[논문 리뷰] A Second Order Energy Stable Scheme for the Cahn-Hilliard-Hele-Shaw Equations

Wenbin Chen, Wenqiang Feng외 3명|arXiv (Cornell University)|2016. 11. 09.

Solidification and crystal growth phenomena참고 문헌 12인용 수 9

한 줄 요약

이 논문은 시간 이산화에 수정된 크랭크-니콜슨과 애담스-مول턴 방법을 사용하여 Cahn-Hilliard-Hele-Shaw 방정식에 대해 두 번째 차수 정확도를 가지며, 조건부 에너지 안정성을 보장하는 유한차분 스킴을 제안한다. 비선형 항은 볼록 분할과 반음향적 처리를 통해 다루어지며, 이 스킴은 고유한 해 존재성을 보장하고, $\varepsilon$-독립적인 $\mathbb{Q}^{\infty}(0,T;H_{h}^{1})$ 오차 추정을 통해 최적 수렴률을 달성한다. 수치 시뮬레이션을 통해 검증되었다.

ABSTRACT

We present a second-order-in-time finite difference scheme for the Cahn-Hilliard-Hele-Shaw equations. This numerical method is uniquely solvable and unconditionally energy stable. At each time step, this scheme leads to a system of nonlinear equations that can be efficiently solved by a nonlinear multigrid solver. Owing to the energy stability, we derive an $\ell^2 (0,T; H_h^3)$ stability of the numerical scheme. To overcome the difficulty associated with the convection term $ abla \cdot (ϕ\boldsymbol{u})$, we perform an $\ell^\infty (0,T; H_h^1)$ error estimate instead of the classical $\ell^\infty (0,T; \ell^2)$ one to obtain the optimal rate convergence analysis. In addition, various numerical simulations are carried out, which demonstrate the accuracy and efficiency of the proposed numerical scheme.

연구 동기 및 목표

• Cahn-Hilliard-Hele-Shaw 시스템에 대해 조건부 에너지 안정성을 유지하는 두 번째 차수 시간 정확도를 갖는 수치 스킴을 개발한다.
• 전통적인 $\mathbb{Q}^{\infty}(0,T;\ell^{2})$ 대신 $\mathbb{Q}^{\infty}(0,T;H_{h}^{1})$ 오차 추정을 사용하여, 에너지 안정성과 수렴 분석에서 대류 항 처리의 과제를 해결한다.
• 이 시스템에 대해 두 번째 차수 스킴의 최적 수렴률 분석을 수립한다. 이는 문헌에서 흔치 않은 일이다.
• 각 시간 단계에서 비선형 다중격자 솔버를 통해 고유한 해 존재성과 강건성을 확보한다.
• 2차원 및 3차원 영역에서 광범위한 수치 시뮬레이션을 통해 스킴의 정확성과 효율성을 검증한다.

제안 방법

• 시간에 대한 두 번째 차수 정확도를 확보하기 위해 비선형 화학적 환경 항에 수정된 크랭크-니콜슨 스킴을 적용한다.
• 최고차수 표면 확산 항에 대해 시간 단계 $t^{n+1}$ 및 $t^{n-1}$ 에서 애담스-몰턴 보간 공식을 사용하여 안정성과 수렴 분석을 향상시킨다.
• 에너지의 오목 부분은 애담스-바슈포스 외삽을 통해 명시적으로 처리하여 에너지 안정성을 유지한다.
• 대류 항 $\nabla \cdot (\phi \mathbf{u})$ 는 반음향적으로 다루며, 단계 변수는 외삽되고 속도장은 다르시의 법칙에 따라 암시적으로 갱신된다.
• 에너지의 볼록 부분과 오목 부분으로 분할하는 볼록 분할 전략을 사용하여 조건부 에너지 안정성을 달성한다.
• 각 시간 단계에서 유도된 비선형 연립방정식은 비선형 다중격자 솔버를 통해 효율적으로 해결된다.

실험 결과

연구 질문

• RQ1 Cahn-Hilliard-Hele-Shaw 시스템에 대해 조건부 에너지 안정성과 고유한 해 존재성을 유지하는 두 번째 차수 시간 정확도를 갖는 유한차분 스킴을 구성할 수 있는가?
• RQ2 대류 항 $\nabla \cdot (\phi \mathbf{u})$ 는 어떻게 다뤄야 하며, 이는 전통적인 $\ell^{\infty}(0,T;\ell^{2})$ 보다 더 자연스러운 노름에서 최적 수렴 분석을 가능하게 하는가?
• RQ3 이 시스템에서 두 번째 차수 스킴이 달성할 수 있는 최적 수렴률은 무엇이며, 수정된 시간 이산화 전략을 사용하여 엄밀히 증명할 수 있는가?
• RQ4 스킴이 시간 간격 크기에 관계없이 에너지 안정성을 유지할 수 있는가, 심지어 큰 시간 간격에서도 말이다?
• RQ5 제안된 스킴은 2차원 및 3차원 벤치마크 문제에서 정확성과 효율성 측면에서 어떻게 성능을 발휘하는가?

주요 결과

• 제안된 스킴은 조건부 에너지 안정성을 보장하며, 시간 간격 크기에 관계없이 이산 에너지가 감소한다.
• 각 시간 단계에서 스킴은 고유한 해 존재성을 보장하여 잘 정의된 수치적 해를 보장한다.
• $\ell^{2}(0,T;H_{h}^{3})$ 안정성 추정이 유도되었으며, 이는 수치적 해의 높은 정(regularity)을 나타낸다.
• $\ell^{\infty}(0,T;H_{h}^{1})$ 오차 추정이 확립되어, 전통적인 $\ell^{\infty}(0,T;\ell^{2})$ 대비 완화된 노름에서 최적의 두 번째 차수 수렴을 가능하게 한다.
• 수치 시뮬레이션은 시간에 대한 두 번째 차수 수렴을 확인하였으며, 2차원 및 3차원 영역에서 강건성과 효율성을 입증한다.
• 비선형 다중격자 솔버는 각 시간 단계에서 시스템을 효율적으로 해결하여 스킴의 실용적 타당성을 유지한다.