[논문 리뷰] Stabilized exponential-SAV schemes preserving energy dissipation law and maximum bound principle for the Allen-Cahn type equations

Lili Ju, Xiao Li외 1명|arXiv (Cornell University)|2022. 03. 14.

Solidification and crystal growth phenomena인용 수 4

한 줄 요약

이 논문은 앨런-클로우 타입 방정식에 대해 이산 에너지 소산 법칙과 최대 경계 원칙(MBP)을 동시에 유지하는 최초의 1차 및 2차 안정화된 지수 스칼라 보조변수(sESAV) 스킴을 제안한다. 지수 시간 적분과 새로운 안정화 기법을 조합함으로써, 스킴은 조건부 에너지 안정성과 стр한 해 경계를 보장하며, 최적의 오차 추정치와 광범위한 수치적 검증을 제공한다.

ABSTRACT

It is well-known that the Allen-Cahn equation not only satisfies the energy dissipation law but also possesses the maximum bound principle (MBP) in the sense that the absolute value of its solution is pointwise bounded for all time by some specific constant under appropriate initial/boundary conditions. In recent years, the scalar auxiliary variable (SAV) method and many of its variants have attracted much attention in numerical solution for gradient flow problems due to their inherent advantage of preserving certain discrete analogues of the energy dissipation law. However, existing SAV schemes usually fail to preserve the MBP when applied to the Allen-Cahn equation. In this paper, we develop and analyze new first- and second-order stabilized exponential-SAV schemes for a class of Allen-Cahn type equations, which are shown to simultaneously preserve the energy dissipation law and MBP in discrete settings. In addition, optimal error estimates for the numerical solutions are rigorously obtained for both schemes. Extensive numerical tests and comparisons are also conducted to demonstrate the performance of the proposed schemes.

연구 동기 및 목표

• 앨런-클로우 타입 방정식의 시간 이산화 설정에서 에너지 소산 법칙과 최대 경계 원칙(MBP)을 동시에 유지하는 수치적 스킴을 개발한다.
• 일반적으로 MBP를 유지하지 못하는 표준 SAV 스킴의 한계를 극복하기 위해 안정화된 지수-SAV 프레임워크를 도입한다.
• 제안된 스킴에 대한 엄밀한 안정성 및 수렴 분석을 수립하여 최적의 오차 한계를 보장한다.
• 장기간의 응집 동역학과 상분리 과정에 대한 광범위한 수치 실험을 통해 스킴의 효능을 입증한다.
• 새로운 안정화 메커니즘을 통해 SAV 기반 방법의 적용 가능성을 기울기 흐름으로까지 확장하여 에너지 안정성과 해의 유한성 양쪽을 보장한다.

제안 방법

• 에너지 기능을 재구성하기 위해 지수 스칼라 보조변수(ESAV) 접근법을 채택하여 이산 스킴에서 비선형 항에 대해 양수 계수를 보장한다.
• 스킴에 선형 항을 더하고 빼는 방식으로 안정화 항을 도입하며, 이는 MBP 유지 증명에 핵심적이며 원래 에너지 구조를 변화시키지 않는다.
• 일단 단계 방법 프레임워크를 사용해 1차 및 2차 시간 이산화를 구성함으로써 적응형 시간 스텝 전략의 사용을 가능하게 한다.
• 에너지 소산과 MBP 유지가 서로 독립적이며, 고전적 안정화 스킴과 달리 동시에 증명될 수 있음을 증명한다.
• 안정성과 정확도를 향상시키기 위해 비선형 항을 지수 시간 적분을 사용한 반음성적 처리 방식을 적용한다.
• 엄밀한 분석을 통해 $L^2$ 및 $L^\frown$ 노름에서 최적의 오차 추정치를 유도하며, 이를 수치적으로 검증한다.

실험 결과

연구 질문

• RQ1 이산 시간에서 에너지 소산 법칙과 최대 경계 원칙(MBP)을 동시에 유지할 수 있는 SAV 기반 스킴을 설계할 수 있는가?
• RQ2 해 경계 유지에 있어 표준 SAV 스킴의 본질적 불안정성을 유지하면서도 에너지 안정성을 유지하기 위해 어떻게 수정할 수 있는가?
• RQ3 지수 시간 적분과 안정화 항이 에너지 소산과 MBP 유지 증명을 독립적으로 가능하게 하는 데서 수행하는 역할은 무엇인가?
• RQ4 이러한 안정화된 지수-SAV 스킴에 대해 최적의 오차 추정치를 엄밀히 도출할 수 있는가?
• RQ5 제안된 스킴은 상분리 및 응집 동역학의 장기간 시뮬레이션에서 어떻게 성능을 발휘하는가?

주요 결과

• 제안된 sESAV 스킴은 1차 및 2차 스킴 모두에 대해 이산 에너지 소산 법칙을 조건 없이 유지한다.
• 스킴은 최대 경계 원칙(MBP)을 엄격히 준수하여, 초기 자료가 $|u_{\text{init}}| \leq \beta$ 를 만족할 경우 모든 시간에 걸쳐 해가 절대값으로 $\beta$ 이내에 유지됨을 보장한다.
• 1차 스킴에 대해 $\mathcal{O}(\tau + h^2)$, 2차 스킴에 대해 $\mathcal{O}(\tau^2 + h^2)$ 의 최적 오차 추정치가 엄밀히 도출되었다.
• 수치 결과는 $\tau = 0.01$ 인 sESAV2 스킴이 장기간의 응집 동역학을 정확히 포착함을 보여주며, 이중우물 잠재력의 경우 $t \approx 604$, 플로리-후긴스 잠재력의 경우 $t \approx 602$ 에 안정 상태에 도달한다.
• 모든 시뮬레이션 기간 동안 해가 $L^\infty$ 노름에서 유한하게 유지됨을 확인하여 MBP 유지가 확인되었고, 에너지가 단조 감소함을 통해 에너지 소산이 검증되었다.
• 스킴은 강건하며 비물리적 진동이나 해 경계 위반 없이 IFRK4 스킴과 유사한 결과를 생성한다.