[논문 리뷰] Implicit-explicit BDF $k$ SAV schemes for general dissipative systems and their error analysis

Fukeng Huang, Jie Shen|arXiv (Cornell University)|2021. 03. 10.

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

한 줄 요약

이 논문은 일반적인 소산 시스템, 특히 Allen-Cahn 및 Cahn-Hilliard 유형 방정식을 포함하여 조건부 안정성 없이 고차수(k=1에서 5까지)의 음성-양성 BDF$k$ 스칼라 보조변수(SAV) 스킴을 제안한다. 특별한 제어 인자 $\eta_k^{n+1}$ 를 도입함으로써 에너지 기반 노름에 대한 균일한 유계를 확보하여, 비선형 항에 대한 전역 리프시츠 조건을 요구하지 않고도 $L^\infty(0,T;H^2)$ 노름에서 엄밀한 최적 오차 추정을 가능하게 한다.

ABSTRACT

We construct efficient implicit-explicit BDF$k$ scalar auxiliary variable (SAV) schemes for general dissipative systems. We show that these schemes are unconditionally stable, and lead to a uniform bound of the numerical solution in the norm based on the principal linear operator in the energy. Based on this uniform bound, we carry out a rigorous error analysis for the $k$th-order $(k=1,2,3,4,5)$ SAV schemes in a unified form for a class of typical Allen-Cahn type and Cahn-Hilliard type equations. We also present numerical results confirming our theoretical convergence rates.

연구 동기 및 목표

• 일般적인 소산 시스템에 대해 고차수, 조건부 안정성 없는 시간 적분 스킴을 개발하는 것.
• 기존의 SAV 방법이 고차수 스킴에 대해 엄밀한 오차 분석을 부족하게 하는 한계를 극복하는 것.
• Allen-Cahn 및 Cahn-Hilliard 유형 방정식에 적용 가능한 $k$차수 BDF$^k$ SAV 스킴에 대한 통합된 오차 분석 프레임워크를 수립하는 것.
• 각 시간 단계에서 상수 계수를 가진 선형 연립방정식 한 번만 풀면 되도록 계산 효율성을 확보하는 것.
• $L^\infty(0,T;H^2)$ 노름에서 비선형 항에 대한 전역 리프시츠 가정 없이 최적 오차 추정을 제공하는 것.

제안 방법

• 에너지 기반 노름에 대한 균일한 유계를 확보하기 위해 SAV 설정에 새로운 제어 인자 $\eta_k^{n+1}$ 를 도입한다.
• BDF$k$ 시간 이산화와 비선형 항에 대한 암시-명시(IMEX) 처리를 조합한다.
• SAV 방법을 일반적인 소산 시스템에 적용하여 보조변수를 통해 원래의 에너지 함수를 이차형으로 변환한다.
• BDF$k$ 스킴에 대한 안정성 결과 [13]와 균일한 $H^1$ 유계를 조합하여 오차 분석을 가능하게 한다.
• 엄밀한 귀납적 추론을 사용하여 $L^\infty(0,T;H^2)$ 노름에서 최적 오차 추정을 유도한다.
• 수치적 검증을 통해 $k=1$에서 $5$까지의 이론적 수렴률이 확인된다.

실험 결과

연구 질문

• RQ1 일반적인 소산 시스템에 대해 조건부 안정성 없이 고차수(k=1에서 5까지)의 BDF$^k$ SAV 스킴을 구성할 수 있는가?
• RQ2 고차수 SAV 스킴에 대해 에너지 기반 노름에 대한 균일한 유계를 확보하기 위해 어떤 방법을 사용할 수 있는가?
• RQ3 각 시간 단계에서 상수 계수를 가진 선형 연립방정식 한 번만 풀면 되도록 하여 계산 효율성을 유지할 수 있는가?
• RQ4 $L^\infty(0,T;H^2)$ 노름에서 비선형 항의 전역 리프시츠 연속성 가정 없이 최적 오차 추정을 유도할 수 있는가?
• RQ5 $k=1$에서 $5$까지의 수치 결과가 이론적 수렴률을 확인하는가?

주요 결과

• 제안된 $k$차수 BDF$^k$ SAV 스킴은 조건부 안정성을 가지며, 수치 해의 $H^1$-노름에 대해 균일한 유계를 확보한다.
• $L^\infty(0,T;H^2)$ 노름에서 Allen-Cahn 및 Cahn-Hilliard 유형 방정식에 대해 최적 오차 추정이 확립된다.
• 스킴은 각 시간 단계에서 상수 계수를 가진 선형 연립방정식 한 번만 풀면 되므로 계산 효율성이 보장되며, 표준 IMEX 스킴과 유사한 성능을 보인다.
• 비선형 항의 전역 리프시츠 연속성 가정 없이 오차 추정이 도출되어 SAV 방법의 적용 범위가 확장된다.
• $k=1$에서 $5$까지의 수치 결과는 이론적 수렴률을 확인하여 스킴의 정확성과 안정성을 검증한다.
• 이 방법은 고차수 정확도와 조건부 안정성을 확보하면서도 적응형 시간 스텝에 적합하다.