[논문 리뷰] Efficient schemes with unconditionally energy stability for the anisotropic Cahn-Hilliard Equation using the stabilized-Scalar Augmented Variable (S-SAV) approach

Xiaofeng Yang|arXiv (Cornell University)|2018. 04. 08.

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

한 줄 요약

이 논문은 안정화된 스칼라 보조 변수(S-SAV) 방법을 사용하여 비등방성 Cahn-Hilliard 방정식에 대해 두 가지 이阶, 조건 없이 에너지 안정적인 시간 이산화 방법을 제안한다. SAV 방법과 세 가지 선형 안정화 항을 조합함으로써, 모든 비선형성을 반명시적으로 처리하고 각 시간 스텝에서 상수 계수를 가진 세 개의 분리된 선형 방정식만을 풀게 되어, 강한 비등방성으로 인한 진동을 제거하면서도 큰 시간 스텝을 허용할 수 있다. 이 방법은 이 모델에 대해 조건 없이 에너지 안정성이 보장되고 이阶 정확도를 가지며 높은 효율성을 갖는 첫 번째 증명 가능한 방법이다.

ABSTRACT

In this paper, we consider numerical approximations for the anisotropic Cahn-Hilliard equation. The main challenge of constructing numerical schemes with unconditional energy stabilities for this model is how to design proper temporal discretizations for the nonlinear terms with the strong anisotropy. We propose two, second order time marching schemes by combining the recently developed SAV approach with the linear stabilization approach, where three linear stabilization terms are added. These terms are shown to be crucial to remove the oscillations caused by the anisotropic coefficients, numerically. The novelty of the proposed schemes is that all nonlinear terms can be treated semi-explicitly, and one only needs to solve three decoupled linear equations with constant coefficients at each time step. We further prove the unconditional energy stabilities rigorously, and present various 2D and 3D numerical simulations to demonstrate the stability and accuracy.

연구 동기 및 목표

• 강한 비등방성으로 인한 수치적 진동으로 인해 도전적인 비등방성 Cahn-Hilliard 방정식에 대해 효율적이고 조건 없이 에너지 안정적인 수치적 방법을 개발한다.
• 기존 방법의 한계—예를 들어 증명 가능한 에너지 안정성 부족이나 높은 계산 비용—을 극복하기 위해 SAV 방법과 선형 안정화를 조합한다.
• 각 시간 스텝에서 상수 계수를 가진 선형 방정식만을 풀음으로써 시간에 대한 이계 정확도를 유지하면서도 계산 효율성을 확보한다.
• 강한 비등방성 조건에서도 안정성 없이 큰 시간 스텝을 허용할 수 있음을 입증한다.
• 제안된 방법에 대해 조건 없이 에너지 안정성의 엄밀한 수학적 증명을 제공한다. 이는 이전 방법들과의 핵심적 발전이다.

제안 방법

• 에너지 함수를 새로운 스칼라 보조 변수를 도입하여 이차형으로 재구성함으로써 SAV(스칼라 보조 변수) 방법을 적용하여 비선형 항을 반명시적으로 처리할 수 있도록 한다.
• 특히 비등방성 계수 γ(n)에 의해 유도되는 고주파수 진동을 억제하기 위해 세 가지 선형 안정화 항(S1, S2, S3)을 도입한다.
• 에너지 소산 법칙을 조건 없이 유지하는 두 가지 이계 BDF2 유형의 시간 이산화 방법(Schemes BDF2-S 및 BDF2-W)을 구성한다.
• 각 시간 스텝에서 상수 계수를 가진 세 개의 독립된 선형 타원형 방정식으로 시스템을 분리하여 표준 선형 해법기로 풀 수 있도록 한다.
• 시험 함수를 시험 함수 공간과 동일한 공간에 놓고 변분 형식을 사용하여 약한 형태를 유도하고, 에너지 추정을 통해 안정성을 확보한다.
• 특히 3차원 시뮬레이션에서 수치 조건을 향상시키고 진동을 감소시키기 위해 Willmore 정규화 항(βG(φ))을 적용한다.

실험 결과

연구 질문

• RQ1 강한 비등방성으로 인해 조건 없이 에너지 안정성을 보장하면서도 이계 시간 이산화 방법을 비등방성 Cahn-Hilliard 방정식에 적용할 수 있는가?
• RQ2 비등방성 계수 γ(∇φ/|∇φ|)를 포함하는 비선형 항은 어떻게 효율적이고 안정적으로 반명시적으로 처리할 수 있는가?
• RQ3 추가적인 선형 안정화 항은 비등방성 계수로 인한 수치적 진동 억제에 어떤 역할을 하는가?
• RQ4 제안된 방법은 강한 비등방성 조건에서도 큰 시간 스텝을 사용하면서도 높은 정확도와 안정성을 유지할 수 있는가?
• RQ5 각 시간 스텝에서 상수 계수를 가진 선형 방정식만을 풀면서도 조건 없이 에너지 안정성이 보장되는 것이 가능한가?

주요 결과

• 제안된 S-SAV 방법은 변분 형식에 기반한 에너지 추정을 통해 조건 없이 에너지 안정성이 보장됨을 엄밀히 증명하였다.
• 수치 실험 결과, S1 = S2 = 2이고 S3 = 0.001인 조합은 고주파수 진동을 효과적으로 억제하는 반면, S1 = S2 = 0이고 S3 = 0.001인 조합은 비물리적 진동을 유도한다.
• 시간 스텝 δt = 1e−4를 유지함으로써 정확도를 확보할 수 있으며, δt = 1e−2, 1e−3, 1e−4의 결과는 거의 동일한 계면 영역 프로파일을 보였다.
• 비등방성 매개변수 α를 0.1에서 0.3으로 증가시킬수록 평형 형태에서 피라미드 모양의 모서리가 더 날카롭게 나타나며, 물리적 기대와 일치한다.
• Willmore 정규화 매개변수 β를 감소시킬수록 평형 형태의 모서리가 더 날카롭게 나타나고, 시간 진동도 플롯을 통해 β가 작을수록 에너지 감소 속도가 더 빠르게 나타남을 확인하였다.
• 이 방법은 복잡한 비선형성을 가진 다른 그래디언트 플로우 문제로도 일반화되고 확장 가능하며, 증명 프레임워크가 변분 형식과 일관된 갈레르킨 근사에 기반하기 때문이다.