[論文レビュー] Characterizing the stabilization size for semi-implicit Fourier-spectral method to phase field equations

Dong Li, Zhonghua Qiao+1 more|arXiv (Cornell University)|Oct 28, 2014

Solidification and crystal growth phenomena参考文献 20被引用数 7

ひとこと要約

本稿は、Cahn-Hilliardや分子線エpitaxyモデルなどの相場面方程式に適用される半陽的フーリエスペクトルスキームに対して、非線形性の制限的なリプシッツ条件を除去することで、無条件エネルギー安定性を確立する。エネルギー安定性が時間刻みに依存しない安定化項の大きさに依存することを示し、時間刻みを大きくとることを可能にするとともに、最小限の仮定のもとで誤差解析を可能にする。

ABSTRACT

Recent results in the literature provide computational evidence that stabilized semi-implicit time-stepping method can efficiently simulate phase field problems involving fourth-order nonlinear dif- fusion, with typical examples like the Cahn-Hilliard equation and the thin film type equation. The up-to-date theoretical explanation of the numerical stability relies on the assumption that the deriva- tive of the nonlinear potential function satisfies a Lipschitz type condition, which in a rigorous sense, implies the boundedness of the numerical solution. In this work we remove the Lipschitz assumption on the nonlinearity and prove unconditional energy stability for the stabilized semi-implicit time-stepping methods. It is shown that the size of stabilization term depends on the initial energy and the perturba- tion parameter but is independent of the time step. The corresponding error analysis is also established under minimal nonlinearity and regularity assumptions.

研究の動機と目的

• 既存の半陽的フーリエスペクトル法のエネルギー安定性解析において、非線形ポテンシャルの制限的なリプシッツ条件を除去すること。
• 非線形性および正則性の仮定を最小限に抑えた条件下で、Cahn-Hilliardや分子線エpitaxyモデルなどの相場面方程式に対して、無条件エネルギー安定性を確立すること。
• 安定化項の大きさが時間刻みに依存しないことを特徴づけ、大きな時間刻みをとることで安定性を維持することを可能にすること。
• 非線形性および正則性の仮定を最小限に抑えた条件下で、厳密な誤差解析を提供し、4階相場面問題の効率的数値シミュレーションの理論的基盤を拡張すること。

提案手法

• 4階相場面方程式に適した、O(Δt)の安定化項を線形化系に追加した安定化半陽的フーリエスペクトル時間積分スキームを提案する。
• 変分的定式化とエネルギー法を用いて、非線形ポテンシャルがリプシッツ条件を満たす必要がない条件でも、無条件エネルギー安定性を証明する。
• 安定化項の大きさを初期エネルギーと摂動パラメータの関数として導出し、時間刻みΔtに依存しないことを示す。
• エネルギー法を用いて事前推定を導出し、最小限の正則性および非線形性仮定のもとで誤差境界を確立する。
• 周期的領域Ω = T²を想定し、スペクトル精度を持つフーリエスペクトル法を用いて空間離散化を処理する。
• 高次スキームにおける一貫性を保証するための射影作用素Π_Nを導入し、解析を高次時間積分スキームへと拡張する。

実験結果

リサーチクエスチョン

• RQ1 非線形ポテンシャルがリプシッツ連続でない場合でも、半陽的フーリエスペクトルスキームに対して無条件エネルギー安定性を証明できるか？
• RQ2 このようなスキームにおいて、安定化項の大きさが初期エネルギー、摂動パラメータ、時間刻みにどのように依存するか？
• RQ3 非線形性および正則性の仮定を最小限に抑えた条件下で、これらの安定化スキームに対して誤差推定をどのように確立できるか？
• RQ4 安定化項を時間刻みに依存させずに、エネルギー安定性および精度を維持できるか？
• RQ5 理論的枠組みを高次時間積分スキームおよびより一般的な相場面モデルへと拡張できるか？

主な発見

• 非線形ポテンシャルがリプシッツ条件を満たさなくても、安定化半陽的フーリエスペクトル法に対して無条件エネルギー安定性が証明された。
• 安定化項の大きさは初期エネルギーと摂動パラメータにのみ依存し、時間刻みΔtとは無関係である。
• エネルギーの減衰は単調的かつ無条件的であり、導出された安定化条件下では任意の時間刻みサイズで安定性が保証される。
• 非線形性および正則性の仮定を最小限に抑えた条件下で誤差解析が確立され、スキームの収束保証が得られた。
• 理論的結果は、計算的証拠と整合的であり、シミュレーションにおいて大きな時間刻みの使用を支持する。
• この枠組みは、高次時間積分スキームおよびより一般的な相場面モデル（6階およびベクトル型モデルを含む）へと拡張可能である。