[논문 리뷰] Numerical approximations of the Cahn-Hilliard and Allen-Cahn Equations with general nonlinear potential using the Invariant Energy Quadratization approach

Xiaofeng Yang, Guodong Zhang|arXiv (Cornell University)|2017. 12. 07.

Solidification and crystal growth phenomena참고 문헌 19인용 수 29

한 줄 요약

이 논문은 일반적인 비선형 포텐셜을 갖는 Cahn-Hilliard 및 Allen-Cahn 방정식에 대해 Invariant Energy Quadratization (IEQ) 접근법을 사용하여 일阶 반연산 수치 스킴을 제시한다. 새로운 변수를 도입하여 비선형 포텐셜을 이차형태로 재구성함으로써, 선형이고 조건부 에너지 안정적인 스킴을 도출한다. 주요 기여는 포텐셜의 유계성과 연속성에 대한 온화하고 물리적으로 타당한 조건 하에서 최적의 오차 추정치를 엄밀히 유도한 점이며, 이 조건들은 이중우물 및 정규화된 Flory-Huggins 유형의 일반적인 포텐셜들에 의해 만족된다.

ABSTRACT

In this paper, we carry out stability and error analyses for two first-order, semi-discrete time stepping schemes, which are based on the newly developed Invariant Energy Quadratization approach, for solving the well-known Cahn-Hilliard and Allen-Cahn equations with general nonlinear bulk potentials. Some reasonable sufficient conditions about boundedness and continuity of the nonlinear functional are given in order to obtain optimal error estimates. These conditions are naturally satisfied by two commonly used nonlinear potentials including the double-well potential and regularized logarithmic Flory-Huggins potential. The well-posedness, unconditional energy stabilities and optimal error estimates of the numerical schemes are proved rigorously.

연구 동기 및 목표

• 일반적인 비선형 부피 포텐셜을 갖는 Cahn-Hilliard 및 Allen-Cahn 방정식에 대해 조건부 에너지 안정적인 수치 스킴을 개발하는 것.
• 특히 일반적인 비선형 포텐셜에 대해 IEQ 유형 스킴의 엄밀한 오차 분석 부족 문제를 해결하는 것.
• 포텐셜의 유계성과 연속성에 대한 충분조건을 규명하여 최적의 오차 추정치를 수립하는 것.
• 다양한 비선형 포텐셜을 갖는 다양한 기울기 흐름 모델에 적용 가능한 일반적인 분석 프레임워크를 제공하는 것.
• IEQ 방법의 적용 범위를 에너지 안정성 초과하여 최적 수렴 속도를 갖는 수렴 분석으로 확장하는 것.

제안 방법

• 보조 변수를 사용하여 비선형 포텐셜을 이차형태로 변환하기 위해 Invariant Energy Quadratization (IEQ) 접근법을 적용한다.
• 기존의 PDE를 새로운 변수에 대해 재구성하여 에너지 소산 법칙이 연속 수준에서 그대로 유지되도록 한다.
• 시간 도함수를 일阶 반암시적 스킴으로 이산화하여 선형 연립방정식 시스템을 도출한다.
• 결과로 얻어진 선형 시스템이 구성에 의해 잘 정의되고 조건부 에너지 안정적이도록 보장한다.
• 시험 함수를 시험 함수와 동일한 함수 공간에 놓고 변분 형식을 사용하여 오차 추정치를 유도한다.
• 수학적 귀납법과 이차형 함수의 리프시츠 연속성을 활용하여 새로운 변수에서의 오차를 제한하고 원래 해로의 오차로 환원한다.

실험 결과

연구 질문

• RQ1 일반적인 비선형 포텐셜을 갖는 Cahn-Hilliard 및 Allen-Cahn 방정식에 적용된 IEQ 기반 스킴에 대해 최적의 오차 추정치를 엄밀히 도출할 수 있는가?
• RQ2 최적의 오차 제어를 위해 필요한 포텐셜의 유계성과 연속성에 대해 최소한의 조건은 무엇인가?
• RQ3 IEQ 변환은 비선형성이 존재할 경우 오차 전파 및 수렴 행동에 어떻게 영향을 미치는가?
• RQ4 IEQ 프레임워크는 에너지 안정성 초과하여 특히 일阶 시간 이산화에 대해 오차 추정치를 제공할 수 있는가?
• RQ5 유도된 오차 추정치는 표준 포텐셜인 이중우물 및 정규화된 Flory-Huggins 포텐셜에 대해 어느 정도 유효한가?

주요 결과

• 일阶 IEQ 스킴을 Cahn-Hilliard 및 Allen-Cahn 방정식에 적용한 결과, 오차 추정치의 최적 순서 𝒪(δt)를 확립하였다.
• 오차 추정치의 유계성과 연속성 조건이 충분히 충족되며, 이는 이중우물 및 정규화된 Flory-Huggins 포텐셜에 의해 자연스럽게 만족된다.
• 스킴은 조건부 에너지 안정적이며, 시간 간격이나 메esh 해상도에 관계없이 에너지 감소가 유지됨을 의미한다.
• IEQ 방법에서 유도된 선형 시스템은 잘 정의되어 있으며, Allen-Cahn 방정식의 경우 대칭 양의 정부호 행렬이다.
• 충분히 작은 시간 간격에 대해 ∥φn∥L∞ ≤ κ 라는 균일한 유계성 조건이 확립되었으며, 이는 오차 분석에서 비선형 항을 제어하는 데 필수적이다.
• 변분 형식에 기반한 이 분석 프레임워크는 일致한 갈레르킨 유형의 공간 이산화에 대해 일반적으로 확장 가능하다.