[論文レビュー] Stability and Error estimates of the SAV Fourier-spectral method for the Phase Field Crystal Equation

Xiaoli Li, Jie Shen|arXiv (Cornell University)|Jul 17, 2019

Solidification and crystal growth phenomena参考文献 16被引用数 4

ひとこと要約

本稿では、相場結晶（PFC）方程式に対する安定化スカラーサービス変数（SAV）フーリエスペクトル法を提案する。時間方向にSAVアプローチを適用し、空間方向にスペクトル離散化を組み合わせることで、無条件エネルギー安定性で、かつ線形なスキームを達成する。本研究は、完全離散線形スキームに対する最初の厳密な最適誤差推定値を提供し、時間方向に2次収束、空間方向にスペクトル精度を示す。安定化は、時間ステップが大きい場合の精度を維持するために不可欠であるが、安定性のためには必要ではない。

ABSTRACT

We consider fully discrete schemes based on the scalar auxiliary variable (SAV) approach and stabilized SAV approach in time and the Fourier-spectral method in space for the phase field crystal (PFC) equation. Unconditionally energy stability is established for both first- and second-order fully discrete schemes. In addition to the stability, we also provide a rigorous error estimate which shows that our second-order in time with Fourier-spectral method in space converges with order $O(Δt^2+N^{-m})$, where $Δt$, $N$ and $m$ are time step size, number of Fourier modes in each direction, and regularity index in space, respectively. We also present numerical experiments to verify our theoretical results and demonstrate the robustness and accuracy of the schemes.

研究の動機と目的

• 相場結晶（PFC）方程式に対する無条件エネルギー安定性で、かつ線形な数値スキームの開発を通じて、非線形スキームの複雑さを克服すること。
• 特に線形スキームに関して、既存のPFC数値法における収束解析と誤差推定の欠如に対処すること。
• SAVフレームワークとスペクトル空間離散化を用いて、完全離散スキームの厳密な誤差推定を確立すること。
• 安定化の役割が、安定性を損なわずに大きな時間ステップでも精度を維持する上で果たす役割を調査すること。
• 数値実験を通じて、提案スキームの頑健性と精度を実証すること。

提案手法

• PFC方程式を線形系に再定式化するため、スカラーサービス変数（SAV）アプローチを適用し、エネルギー安定性を達成するとともに、時間離散化を簡素化する。
• 周期的境界条件を活用して、フーリエスペクトル法を空間離散化に用い、空間方向にスペクトル精度を達成する。
• 安定化項を追加することで、S-SAV（安定化SAV）定式化を導入し、大きな時間ステップでも精度を向上させる。
• Crank-Nicolson法または後退差分公式（BDF）を時間方向に用いて完全離散スキームを構築し、無条件エネルギー安定性を保証する。
• エネルギー安定性に基づく離散解の一様有界性を導出し、厳密な誤差解析を可能にする。
• 離散エネルギー安定性および正則性仮定に基づく誤差解析を実施し、最適収束率を導出する。

実験結果

リサーチクエスチョン

• RQ1 SAVアプローチをPFC方程式に拡張することで、最適収束率を達成する線形で、無条件エネルギー安定なスキームが得られるか？
• RQ2 S-SAV定式化における安定化項が、大きな時間ステップでも正確な結果を得るために果たす役割は何か？
• RQ3 PFC方程式に完全離散線形スキームを適用した場合、高次の非線形性を考慮しても、厳密な誤差推定値を導出可能か？
• RQ4 時間ステップサイズの選択が、特に安定化を適用した・しない場合に、SAV-フーリエスペクトルスキームの精度に与える影響は何か？
• RQ5 スペクトル空間離散化が、スキーム全体の収束速度および数値的効率に果たす貢献度はどの程度か？

主な発見

• 提案された2次S-SAV-フーリエスペクトルスキームは、それぞれ時間方向に $ O( au^2 + N^{-m}) $、空間方向に $ O( au^2 + N^{-m}) $ の最適収束率を達成する。ここで $ \tau $ は時間ステップ、$ N $ はフーリエモードの数である。
• 1次および2次完全離散スキームの両方について、無条件エネルギー安定性が証明され、時間ステップサイズに依存しない頑健性が保証される。
• S-SAV定式化における安定化項は、安定性や収束性のためには必要ではないが、大きな時間ステップでも妥当な精度を達成する上で不可欠である。
• 数値実験により、時間方向に2次収束、空間方向にスペクトル収束が確認され、理論的誤差推定値の妥当性が裏付けられた。
• 安定化を適用した場合、大時間ステップでも複雑な結晶成長ダイナミクス（樹状成長、不配位、界面移動など）を正確に捉えることができた。
• 過冷却液体の相転移およびエpitaxial成長のシミュレーションにおいて、数値的拡散を最小限に抑えつつ、頑健で効率的な性能を示した。