3차 정확도 기하학적 부피 보존을 위한 Cahn‑Hilliard 모델

초록

본 논문은 기존 Cahn–Hilliard 상변수 모델에서 인터페이스 두께가 유한할 때 발생하는 기하학적 부피 드리프트 문제를 해결한다. Zhou 등(2022)의 개선된 보존 프레임워크를 확장하여, 물리적 시간 스케일에서의 매치드 비대칭 해석을 수행하고, 1차 내부 교정항을 단일 적분 형태로 정리한다. 두 개의 적분 모멘트를 이용한 ‘모멘트 균형’ 조건을 도출해, Q 매핑 커널을 설계하면 인터페이스 두께 ε에 대해 부피 오차가 O(ε³)까지 감소한다. 또한 에너지 소산을 보장하는 무조건적인 수치 스키마와 다양한 실험을 통해 작은 방울 소멸을 방지하고 부피 보존을 거의 완벽히 달성함을 확인한다.

상세 분석

이 연구는 퇴화된 Cahn–Hilliard 방정식이 표면 확산에 의해 구동되는 인터페이스를 묘사할 때, 전통적인 질량 보존(∫Ωϕ) 대신 설계된 단조 함수 Q(ϕ)를 이용해 정확히 보존되는 “Q‑부피” V_Q(t)=½∫Ω(1+Q(ϕ))dx 를 사용함으로써 기하학적 부피 |Ω⁺(t)|와의 차이를 고차 정확도로 줄이는 방법을 제시한다. 먼저 저자는 물리적 시간(스케일링 없는 시간)에서 매치드 비대칭 전개를 수행한다. 이 경우 인터페이스의 정상 속도 V_n은 ε·V₁+O(ε²) 로 전개되며, V₁은 표면 라플라시안에 비례하는 표면 확산 법칙을 만족한다. 내부 변수 z=d/ε를 도입해 전형적인 헤테로클라인 σ(z)=tanh(z/√2) 를 leading profile 로 잡고, ϕ=σ+εHΦ₁+O(ε²) 로 전개한다. 여기서 Φ₁은 LΦ₁=√(2)/3·Q′(σ)−σ′ 를 만족하는 비동차 방정식이며, L은 self‑adjoint 연산자이다. 기존 문헌에서는 Φ₁을 이중 적분 형태로 제시했지만, 저자는 변수 변환 u=σ(z) 를 이용해 Φ₁=σ′·(4/3)∫₀^{σ(z)}