난류 내 액적 부피 보존을 위한 곡률‑의존형 보정법

초록

본 논문은 고윌버수(Weber number) 조건에서 기존 위상장 모델이 액적 부피를 유지하지 못하는 문제를 지적하고, 보존형 Allen‑Cahn 방정식에 곡률‑의존형 역확산(term) 을 추가한 새로운 보정법을 제안한다. 직접수치시뮬레이션(DNS) 결과, 제안 모델은 통계적 부피 보존을 달성하면서 수치 불안정, 전역 질량 위배, 인공적 코어싱, 스퓨리어스 속도 증가 등의 부작용을 최소화한다.

상세 분석

이 연구는 위상장(phase‑field) 방법을 이용한 액적‑함유 난류 시뮬레이션에서 가장 오래된 난제인 “액적 부피 손실” 문제를 체계적으로 재조명한다. 기존에는 전역 질량 보정, 국소 확산 억제, 이동도(mobility) 변조 등 다양한 보정 전략이 제안되었지만, 대부분이 정적인 혹은 약한 변형 흐름에서만 검증되었으며, 고윌버수(We ≫ 1) 조건의 완전 발달 난류에서는 여전히 0.5 %~2.5 % 수준의 부피 감소가 관측되었다. 특히, 전역 보정은 작은 액적이 빠르게 소멸하면서 발생하는 질량 재분배를 균일하게 적용함으로써 인공적인 코어싱(coarsening)과 전역 질량 위배 위험을 내포한다. 반면, 보존형 Allen‑Cahn(CAC) 방정식은 인터페이스 정상 방향의 확산을 억제해 일정 정도 부피 보존에 성공했지만, 강한 변형·파열 상황에서는 여전히 손실이 발생한다.

저자들은 이러한 한계를 극복하기 위해 CAC 방정식에 “곡률‑의존형 역확산(counter‑diffusion) 항”을 도입한다. 구체적으로, 인터페이스 곡률 κ에 비례하는 보정 계수를 정의하고, 이를 ∇·(κ · ∇ϕ) 형태의 항으로 추가함으로써, 고곡률(작은 액적) 영역에서의 비물리적 확산을 효과적으로 상쇄한다. 이 보정은 물리적으로는 “강수(precipitation)‑유사 메커니즘”이라 해석될 수 있으며, 전역 질량 보전은 기존 CAC의 발산형(term) 로 유지하면서 국소적으로만 작용한다.

수치 검증은 동일한 격자 해상도와 물성 파라미터를 사용한 3차원 동질등방성 난류(HIT) 환경에서 수행되었다. Weber 수를 10, 30, 50 등으로 변화시켰을 때, 기존 모델들은 1 %~3 % 수준의 부피 손실을 보였으나, 제안 모델은 통계적으로 0 %에 근접한 부피 보존을 달성했다. 또한, 스퓨리어스 속도(spurious velocity) 크기와 인터페이스 두께 변화를 정량화한 결과, 새로운 보정이 기존 모델보다 수치적 안정성을 저해하지 않으며, 계산 비용도 기존 CAC와 동일 수준(추가 연산량 < 5 %)에 머물렀다.

이러한 결과는 곡률‑의존형 보정이 고윌버수 난류에서 작은 액적의 과도한 용해를 억제하고, 전체 액적 집합체의 부피 통계량을 유지하는 데 핵심적인 역할을 함을 시사한다. 또한, 보정 항이 간단한 형태이므로 다양한 위상장 기반 DNS 코드에 손쉽게 적용 가능하다는 실용적 장점도 강조된다.

한계점으로는 현재 연구가 동일한 물성(밀도·점도 차이)과 일정한 인터페이스 두께를 가정했으며, 복합적인 물리(예: 열전달, 화학 반응)와 결합된 경우 보정 항의 파라미터 튜닝이 필요할 수 있다. 향후 연구에서는 다중 물성 차이를 포함한 복합 흐름, 그리고 실험 데이터와의 정량적 비교를 통해 보정 계수의 보편성을 검증할 계획이다.