[論文レビュー] Stabilized exponential-SAV schemes preserving energy dissipation law and maximum bound principle for the Allen-Cahn type equations

Lili Ju, Xiao Li+1 more|arXiv (Cornell University)|Mar 14, 2022

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）の両方を保存する数値スキームの開発を目的とする。
• 通常のSAVスキームがMBPを保存できないという限界を克服するため、安定化指数関数的SAVフレームワークを導入する。
• 提案されたスキームに対して、厳密な安定性および収束解析を確立し、最適な誤差境界を保証する。
• 長時間にわたる粗大化ダイナミクスおよび相分離過程の広範な数値実験を通じて、スキームの有効性を示す。
• 新規の安定化機構を用いて、勾配流れへのSAVベースの手法の適用範囲を拡大し、エネルギー安定性と解の有界性の両方を保証する。

提案手法

• エネルギー関数を再定式化するため、指数関数的スカラーアシスタント変数（ESAV）手法を採用し、離散スキームにおける非線形項の係数が正であることを保証する。
• スキームに線形項を加算・減算することで安定化項を導入し、MBPの保存を証明するために不可欠であるが、元のエネルギー構造を変更しない。
• 1段階法フレームワークを用いて1次および2次時間離散化を構築し、適応的時間刻み幅戦略の使用を可能にする。
• エネルギー散逸とMBP保存が独立して同時に確立できることを証明し、従来の安定化スキームとは異なり、並列に証明可能であることを示す。
• 指数時間積分を用いた非線形項の半 implicitly 記述により、安定性と精度を向上させる。
• 厳密な解析を通じて、$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スキームと同等の結果を生成した。