[논문 리뷰] On the Stability and Accuracy of Partially and Fully Implicit Schemes for Phase Field Modeling

Jinchao Xu, Yukun Li외 2명|arXiv (Cornell University)|2016. 04. 19.

Solidification and crystal growth phenomena인용 수 8

한 줄 요약

이 논문은 상변화 모델에 널리 사용되는 볼록 분할 스킴(CSS)이 시간 간격을 크게 줄인 완전히 암시적 스킴(FIS)과 수학적으로 동치이며, 이로 인해 해의 근사가 지연됨을 보여준다. FIS는 에너지 최소화를 통해 더 높은 정확도와 조건부 에너지 안정성 없이도 무조건적인 에너지 안정성을 확보할 수 있으며, 이는 전통적으로 CSS를 선호하는 데 도전하고, 상변화 시뮬레이션에서 FIS의 재평가를 촉구한다.

ABSTRACT

We study in this paper the accuracy and stability of partially and fully implicit schemes for phase field modeling. Through theoretical and numerical analysis of Allen-Cahn and Cahn-Hillard models, we investigate the potential problems of using partially implicit schemes, demonstrate the importance of using fully implicit schemes and discuss the limitation of energy stability that are often used to evaluate the quality of a numerical scheme for phase-field modeling. In particular, we make the following observations: 1. a convex splitting scheme (CSS in short) can be equivalent to some fully implicit scheme (FIS in short) with a much different time scaling and thus it may lack numerical accuracy; 2. most implicit schemes (in discussions) are energy-stable if the time-step size is sufficiently small; 3. a traditionally known conditionally energy-stable scheme still possess an unconditionally energy-stable physical solution; 4. an unconditionally energy-stable scheme is not necessarily better than a conditionally energy-stable scheme when the time step size is not small enough; 5. a first-order FIS for the Allen-Cahn model can be devised so that the maximum principle will be valid on the discrete level and hence the discrete phase variable satisfies $|u_h(x)|\le 1$ for all $x$ and, furthermore, the linearized discretized system can be effectively preconditioned by discrete Poisson operators.

연구 동기 및 목표

• 앨런-캐언 및 캐언-힐리아드 상변화 모델에 대한 부분적으로 암시적(예: 볼록 분할) 및 완전히 암시적 스킴의 안정성과 정확도를 조사한다.
• 볼록 분할 스킴이 완전히 암시적 스킴보다 본질적으로 더 안정하다는 기존의 통념을 도전한다.
• 에너지 최소화를 통한 완전히 암시적 스킴이 시간 간격 제약 없이 무조건적인 에너지 안정성과 향상된 정확도를 달성할 수 있음을 보여준다.
• 시간 간격 크기가 수치적 해에 미치는 영향, 특히 볼록 분할 스킴에서의 시간 지연 효과를 분석한다.
• 최대 원리 보존과 효과적인 조절을 가능하게 하는 수정된 완전히 암시적 스킴을 제안하고 분석한다.

제안 방법

• 이론적 분석을 통해 앨런-캐언 및 캐언-힐리아드 방정식에 대한 볼록 분할 스킴(CSS)과 완전히 암시적 스킴(FIS)을 비교하여 특정 시간 스케일링 하에서 수학적 동치성을 입증한다.
• 앨런-캐언 모델의 경우, 일阶 CSS가 훨씬 작은 시간 간격을 가진 표준 FIS와 동치임을 증명하며, 이로 인해 해의 근사가 지연됨을 보여준다.
• 캐언-힐리아드 모델의 경우, 표준 CSS가 원래 모델의 비정상적인 변형에 해당함을 보여주며, 이는 안정성 향상의 이유를 설명한다.
• 이산 수준에서 최대 원리를 강제하는 수정된 완전히 암시적 스킴을 제안하여 모든 x에 대해 |u_h(x)| ≤ 1를 보장하고, 이산 푸아송 연산자를 통한 효과적인 조절을 가능하게 한다.
• 에너지 최소화를 각 시간 단계에서 이산 최적화 문제로 재구성하며, 전역 최소화해가 무조건적인 에너지 안정성을 만족함을 보여준다.
• 수치 실험을 통해 다양한 시간 간격에서 두 번째 차수 스킴(CSS2, FIS)을 비교하여 시간 지연 효과와 큰 시간 간격에서의 정확도 저하를 시각화한다.

실험 결과

연구 질문

• RQ1 완전히 암시적 스킴이 원칙적으로 더 정확한 데도 불구하고, 볼록 분할 스킴이 왜 더 안정해 보이는가?
• RQ2 볼록 분할 스킴의 안정성이 본질적인 알고리즘적 이점보다는 시간 간격 크기의 효과적 감소 때문인가?
• RQ3 완전히 암시적 스킴을 시간 간격 제약 없이 무조건적인 에너지 안정성으로 만들 수 있는가, 정확도를 훼손하거나 너무 작은 시간 간격을 요구하지 않도록?
• RQ4 시간 간격 크기가 부분적으로 암시적 및 완전히 암시적 스킴의 정확도와 시간적 일치도에 어떻게 영향을 미치는가?
• RQ5 에너지 최소화를 통해 완전히 암시적 스킴의 시간 간격 제약을 제거하고 물리적 충실도를 유지할 수 있는가?

주요 결과

• 앨런-캐언 모델에 대한 일계 볼록 분할 스킴(CSS)은 시간 간격을 훨씬 줄인 표준 완전히 암시적 스킴(FIS)과 수학적으로 동치이며, 이로 인해 해의 근사가 지연됨을 보여준다.
• 캐언-힐리아드 모델의 경우, 표준 CSS는 원래 방정식의 비정상적인 변형에 해당하며, 이는 정확도를 희생하면서까지 안정성 향상을 설명한다.
• 수치 실험을 통해 시간 간격이 충분히 크지 않은 경우, 무조건적인 에너지 안정성 스킴이 조건부 에너지 안정성 스킴보다 나을 수 있음을 보여주며, 이는 에너지 안정성만으로는 좋은 수치 성능을 보장하지 못함을 시사한다.
• 수정된 일계 완전히 암시적 스킴은 이산 수준에서 최대 원리를 보존하여 모든 x에 대해 |u_h(x)| ≤ 1를 보장하며, 이산 푸아송 연산자를 통한 효과적인 조절을 가능하게 한다.
• 완전히 암시적 스킴의 에너지 최소화 재구성은 시간 간격 제약을 제거하고, 큰 시간 간격에서의 CSS보다 물리적으로 더 정확한 해를 도출한다.
• 수치 결과는 조건부 에너지 안정성 스킴조차도 시간 간격이 너무 클 경우 정확도가 떨어질 수 있음을 보여주며, 에너지 안정성만으로는 좋은 수치 성능을 확보할 수 없음을 강조한다.