다상 유동의 다중 접합면을 위한 변분 원리 기반 표면 장력 모델

초록

본 논문은 부피 보존 미분동형사상이라는 무한 차원 리군 위에서 변분 원리를 적용해, 기존 SURFER 방식의 두 상계면 표면 장력식을 일반화한다. 이를 통해 N‑phase( N≥2) 유동에서 다중 접합점(삼중점 등)의 특이성을 자연스럽게 포함한 새로운 오일러 방정식을 도출하고, 접촉각 문제에 대한 수치 실험으로 안정성을 검증한다.

상세 분석

이 연구는 두 단계로 구성된 기존 SURFER(phase‑field) 모델을 수학적 변분 체계에 재배치함으로써, 표면 장력 항을 무한 차원 리군인 부피 보존 미분동형사상(VDiff) 위의 액션에 포함시킨다. 먼저, 두 상계면을 레벨셋 함수 q(x) 로 기술하고, q=0 의 근방을 인공적인 중간 영역(ε) 으로 확장한다. 이 영역에서 표면 장력은 에너지 밀도 σ |∇q| 형태로 정의되며, Noether 정리를 이용해 운동량 보존식과 연계된 응력 텐서를 도출한다. 변분 과정에서 라그랑지안에 운동에너지와 표면 에너지를 더하고, 부피 보존 제약을 라그랑지 승수(압력)로 처리함으로써, 기존 SURFER가 제시한 두 상 오일러 방정식을 정확히 재현한다.

핵심적인 확장은 다중 상을 나타내는 색 함수 φ_a(x) (a=1,…,N) 를 도입하고, 각 상 사이의 인터페이스 에너지를 Σ_{a<b} σ_{ab} |∇φ_a·∇φ_b| 형태로 일반화하는 것이다. 이때 σ_{ab} 는 a‑b 사이의 표면 장력 계수이며, φ_a·φ_b=0 (상호 배타성) 조건을 자연스럽게 만족하도록 설계된다. 변분에 의해 얻어지는 다상 오일러 방정식은 각 상의 질량 보존식, 운동량 방정식, 그리고 인터페이스에 대한 커런시‑조건(연속성) 으로 구성된다. 특히, 다중 접합점에서 발생하는 특이성(예: 삼중점, 사중점)은 중간 영역 ε→0 한계에서도 에너지 밀도가 유한하게 유지되므로, 수치적으로 불안정한 발산을 방지한다.

수치 구현에서는 기존의 VOF/Level‑Set 혼합 기법 대신, φ_a 를 직접 시간 적분하고, 압력-속도 연산을 투영법(Projection)으로 수행한다. 접촉각 문제에서는 고정된 벽면을 추가적인 색 함수 φ_w 로 모델링하고, σ_{aw}, σ_{bw} 와 σ_{ab} 의 비율에 따라 Young‑법칙이 자동으로 만족되는 것을 확인한다. 실험 결과는 주어진 σ_{ij} 값에 대해 접촉각이 이론적 값과 일치하고, 시간 전개가 안정적으로 수렴함을 보여준다.

이 논문의 의의는 (1) 표면 장력 항을 변분 원리와 Noether 정리로 엄밀히 유도함으로써 물리적 일관성을 확보하고, (2) 다중 상 및 다중 접합점을 자연스럽게 포함하는 일반화된 오일러 방정식을 제시했으며, (3) 수치 실험을 통해 제약조건 없이도 안정적인 시뮬레이션이 가능함을 입증했다는 점이다. 향후 복잡한 마이크로플루이드 디바이스, 다상 촉매 반응, 그리고 생물학적 세포막 모델링 등에 바로 적용할 수 있는 강력한 이론적·계산적 토대를 제공한다.