[논문 리뷰] Convergence Analysis of a Linear, Unconditionally Energy-Stable SAV Finite Element Method for the Cahn-Hilliard Equation

Na Li, Yongchao Zhao|arXiv (Cornell University)|2026. 02. 25.

Solidification and crystal growth phenomena인용 수 0

한 줄 요약

논문은 Cahn-Hilliard 방정식에 대해 선형이고 에너지 안정적인 SAV 기반 유한요소 스킴을 개발하고, 무조건 에너지 안정성과 시간 및 공간에서의 최적 차수 수렴을 수치 검증으로 입증한다.

ABSTRACT

This paper proposes a finite element scheme, based on the Scalar Auxiliary Variable (SAV) approach, for the Cahn-Hilliard equation--a model that possesses significant physical relevance and a rich mathematical structure. A convergence analysis of the fully discrete scheme is conducted under suitable regularity assumptions, confirming optimal-order convergence in both time and space for the phase variable, chemical potential, and auxiliary variable in the H1-norm. Furthermore, the scheme is proven to be unconditionally energy stable. Finally, a numerical example is presented to demonstrate the effectiveness of the method and to confirm the theoretical convergence rates.

연구 동기 및 목표

• 에너지 소실을 이산 수준에서도 보존하는 numerically efficient scheme for the Cahn-Hilliard equation를Motivate 하고 분석한다.
• 비선형성을 선형화하면서 안정성을 유지하는 SAV 기반 재정의를 도입한다.
• fully discrete scheme에 대해 무조건 에너지 안정성과 고유 해를 보장한다.
• 1st-order 시간 수렴 및 2nd-order 공간 수렴을 가지는 a priori 오차 추정치를 도출한다.
• 수치 실험으로 이론적 결과와 수렴 속도를 확인한다.

제안 방법

• SAV 기반 시스템을 얻기 위해 r(t)라는 스칼라 보조 변수로 CH 방정식을 재정의한다.
• Backward Euler 시간 이산화와 공간에서의 선형 유한요소를 적용하여 매 단계에서 선형 시스템을 얻는다.
• 이산 에너지 법칙을 통해 무조건 에너지 안정성을 보이고 fully discrete scheme의 고유 해를 보장한다.
• φ, μ, r에 대한 H1-노름에서의 엄밀한 오차 추정을 도출하기 위해 타원 사영(elliptic projection)을 도입한다.
• phase 변수와 화학 퍼텐셜에 대해 시간에서 1차, 공간에서 2차의 수렴을 보이는 수렴 결과를 확립한다.

(a) Variation of the mass with time

실험 결과

연구 질문

• RQ1 SAV 기반 fully discrete scheme이 모든 시간 스텝에 대해 무조건적 에너지 안정성을 유지하는가?
• RQ2 SAV 유한요소 프레임워크에서 φ, μ, 및 보조 변수 r에 대한 a priori 오차 추정은 무엇인가?
• RQ3 적절한 규칙성하에서 SAV-FE 스킴의 시간 및 공간의 수렴 차수는 얼마인가?
• RQ4 타원 사영 기법이 SAV 기반 스킴의 수렴 분석에 어떻게 기여하는가?

주요 결과

• SAV 유한요소 스킴은 이산 에너지가 시간에 따라 단조롭게 감소하는 무조건 에너지 안정성을 가지며 안정적이다.
• 이산화의 선형 구조로 인해 각 시간 스텝에서 고유해를 가진다.
• 최적 차수 수렴이 확립되며: φ와 μ의 H1 노름에서 시간은 1st-order이고 공간은 2nd-order이며 r에 대한 수렴도 적절하다.
• 타원 사영은 엄밀한 오차 추정치를 얻고 투영 오차와 이산화 오차를 분리하는 데 사용된다.
• 수치 실험은 제안된 방법의 이론적 수렴 속도와 안정성을 입증한다.

(b) Variation of the total free energy with time