[論文レビュー] Unconditionally stable exponential time differencing schemes for the mass-conserving Allen-Cahn equation with nonlocal and local effects

Kun Jiang, Lili Ju+2 more|arXiv (Cornell University)|Jul 10, 2021

Solidification and crystal growth phenomena参考文献 39被引用数 9

ひとこと要約

本稿では、局所的および非局所的効果を有する質量保存型アラン＝キャーン方程式に対して、無条件安定な指数時刻差分（ETD）スキームを提案する。線形安定化とETD法を組み合わせることで、最大境界原理（MBP）および質量保存を時間にわたって無条件に保持する1次および2次時間積分スキームを構築した。理論的誤差推定と2次元および3次元における広範な数値検証が行われている。

ABSTRACT

It is well known that the classic Allen-Cahn equation satisfies the maximum bound principle (MBP), that is, the absolute value of its solution is uniformly bounded for all time by certain constant under suitable initial and boundary conditions. In this paper, we consider numerical solutions of the modified Allen-Cahn equation with a Lagrange multiplier of nonlocal and local effects, which not only shares the same MBP as the original Allen-Cahn equation but also conserves the mass exactly. We reformulate the model equation with a linear stabilizing technique, then construct first- and second-order exponential time differencing schemes for its time integration. We prove the unconditional MBP preservation and mass conservation of the proposed schemes in the time discrete sense and derive their error estimates under some regularity assumptions. Various numerical experiments in two and three dimensions are also conducted to verify the theoretical results.

研究の動機と目的

• 質量保存型アラン＝キャーン方程式の無条件安定な時間積分スキームを開発し、最大境界原理（MBP）および質量保存を両方とも時間にわたって保持すること。
• 相場モデルの長時間シミュレーションにおいて、物理的・数学的構造（特に有界性および質量保存）を維持するという挑戦に応えること。
• 正確な質量保存に必要な非局所ラグランジュ乗数項を扱えるように、既存のETDスキームを拡張すること。安定性および精度を保持したまま。
• 正確解および数値スキームの正則性仮定の下で、理論的誤差推定を提示し、2次元および3次元空間において数値的に検証すること。

提案手法

• 数値安定性を向上させるために、質量保存型アラン＝キャーン方程式を線形安定化技術を用いて再定式化する。
• 時間積分のための1次および2次指数時刻差分ルンゲ＝クッタ（ETDRK1およびETDRK2）スキームを構築する。
• 剛性のある非線形項を暗黙的に取り扱うために、指数時刻差分法を適用し、無条件安定性を保証する。
• 積分制約を用いて非局所ラグランジュ乗数項を組み込み、各時間ステップで正確な質量保存を実現する。
• エネルギーに基づくおよびノルムに基づく解析を用いて、時間離散化において最大境界原理（MBP）および質量保存が無条件に保持されることを証明する。
• 正確解および数値スキームの適切な正則性仮定の下で誤差推定を導出する。

実験結果

リサーチクエスチョン

• RQ1 質量保存型アラン＝キャーン方程式に対して、最大境界原理（MBP）を無条件に保持する指数時刻差分スキームを設計できるか？
• RQ2 このようなスキームは、時間離散化設定においても正確な質量保存を保証できるか？
• RQ3 標準的な正則性仮定の下で、提案されたETDスキームの収束性および誤差推定はいかほどか？
• RQ4 2次元および3次元における複雑なダイナミクスを示す相分離の長時間シミュレーションにおいて、これらのスキームはどのように性能を発揮するか？
• RQ5 線形安定化技術は、このクラスの非局所的・非線形PDEに対して、精度と安定性を両立させるために効果的にETD法と組み合わせられるか？

主な発見

• 提案された1次および2次ETDスキームは、最大境界原理（MBP）を無条件に保持し、解がすべての時間において初期境界内で有界であることを保証する。
• 各時間ステップで正確に質量保存が行われ、質量保存型アラン＝キャーン方程式の定式化が要請する条件を満たす。
• 2次元および3次元における数値実験により、エネルギーが単調に減少し、質量が一定に保たれ、解が指定された区間内で有界に保たれることが確認された。
• ETDRK2スキームは時間方向に2次精度を達成し、滑らかな解の下で理論的予測と一致する観察された収束率を示した。
• 膨張する気泡テストでは、気泡の半径が単調に増加し、t ≈ 18で期待される定常状態の半径約0.407に達した。
• 長時間シミュレーションにおいても、スキームは安定的かつ正確に保たれ、安定性のための時間刻み制限が一切不要であることが確認され、無条件安定性が裏付けられた。