[論文レビュー] On the Stability and Accuracy of Partially and Fully Implicit Schemes for Phase Field Modeling

Jinchao Xu, Yukun Li+2 more|arXiv (Cornell University)|Apr 19, 2016

Solidification and crystal growth phenomena被引用数 8

ひとこと要約

本稿は、相場モデルに広く用いられる凸分割スキーム（CSS）が、著しく小さい時間刻みを用いた完全陰解法（FIS）と数学的に同等であることを示しており、その結果、解の近似が遅延する。FISにエネルギー最小化を組み合わせることで、時間刻みの制限なしに優れた精度と無条件エネルギー安定性を達成できることを示し、従来のCSSの優位性を疑問視し、相場モデルのシミュレーションにおけるFISの再評価を提唱する。

ABSTRACT

We study in this paper the accuracy and stability of partially and fully implicit schemes for phase field modeling. Through theoretical and numerical analysis of Allen-Cahn and Cahn-Hillard models, we investigate the potential problems of using partially implicit schemes, demonstrate the importance of using fully implicit schemes and discuss the limitation of energy stability that are often used to evaluate the quality of a numerical scheme for phase-field modeling. In particular, we make the following observations: 1. a convex splitting scheme (CSS in short) can be equivalent to some fully implicit scheme (FIS in short) with a much different time scaling and thus it may lack numerical accuracy; 2. most implicit schemes (in discussions) are energy-stable if the time-step size is sufficiently small; 3. a traditionally known conditionally energy-stable scheme still possess an unconditionally energy-stable physical solution; 4. an unconditionally energy-stable scheme is not necessarily better than a conditionally energy-stable scheme when the time step size is not small enough; 5. a first-order FIS for the Allen-Cahn model can be devised so that the maximum principle will be valid on the discrete level and hence the discrete phase variable satisfies $|u_h(x)|\le 1$ for all $x$ and, furthermore, the linearized discretized system can be effectively preconditioned by discrete Poisson operators.

研究の動機と目的

• アレン＝キャーンおよびカーン＝ヒリャードの相場モデルに対する部分的陰解法（例：凸分割）と完全陰解法の安定性および精度を調査すること。
• 凸分割スキームが完全陰解法よりも本質的に安定であるという一般的な認識に疑問を呈すること。
• エネルギー最小化を組み合わせた完全陰解法が、時間刻みの制限なしに無条件エネルギー安定性を達成し、精度を向上させられることを実証すること。
• 時間刻みの大きさが数値解に与える影響、特に凸分割スキームにおける時間遅延効果を分析すること。
• 最大原理を保持する修正された完全陰解法を提案・分析すること。

提案手法

• 理論的分析により、アレン＝キャーンおよびカーン＝ヒリャード方程式に対する凸分割スキーム（CSS）と完全陰解法（FIS）の比較を行い、特定の時間スケーリング下で数学的に同等であることを示す。
• アレン＝キャーンモデルに関して、1次精度のCSSが、はるかに小さな時間刻みを用いたFISと同等であることを証明し、解の近似が遅延することを示す。
• カーン＝ヒリャードモデルに関して、標準的なCSSが元のモデルの非自明な摂動版に対応することを示し、安定性の向上を説明する。
• 離散レベルで最大原理を満たすように制御する修正された完全陰解法を提案し、|u_h(x)| ≤ 1 をすべての x に対して満たすようにし、離散ラプラシアン作用素による有効なプリコンディショニングを可能にする。
• エネルギー最小化を各時間ステップで離散最適化問題として再定式化し、グローバル最小解が無条件エネルギー安定性を満たすことを示す。
• 時間刻みを変化させた第二階層のスキーム（CSS2, FIS）を比較する数値実験を行い、CSSにおける時間遅延効果と精度の低下を示す。

実験結果

リサーチクエスチョン

• RQ1 完全陰解法が原理的により正確であるにもかかわらず、なぜ凸分割スキームは完全陰解法よりも安定性に優れているのか？
• RQ2 凸分割スキームの安定性は、アルゴリズム的利点よりも、実効的な時間刻みの縮小に起因している程度はどの程度か？
• RQ3 完全陰解法を、精度を損なわず、極めて小さな時間刻みを必要とせずに、無条件エネルギー安定性にできるか？
• RQ4 時間刻みの大きさが、部分的および完全陰解法の両方における精度と時間的整合性にどのように影響するか？
• RQ5 エネルギー最小化を用いることで、完全陰解法における時間刻みの制限を排除し、物理的整合性を維持できるか？

主な発見

• アレン＝キャーンモデルにおける1次精度の凸分割スキーム（CSS）は、はるかに小さな時間刻みを用いた標準的完全陰解法（FIS）と数学的に同等であり、解の近似が遅延することを示す。
• カーン＝ヒリャードモデルに関して、標準的なCSSは元の式の非自明な摂動版に対応しており、安定性の向上は精度の低下を伴う。
• 時間刻みが十分に小さくない場合、無条件エネルギー安定性のスキームが条件付きエネルギー安定性のスキームよりも優れているとは限らないことが、数値実験で示された。
• 修正された1次精度完全陰解法は、離散レベルで最大原理を保持し、すべての x に対して |u_h(x)| ≤ 1 を満たすようにし、離散ラプラシアン作用素による有効なプリコンディショニングを可能にする。
• 完全陰解法のエネルギー最小化再定式化により、時間刻みの制限が解消され、大時間刻みのCSSよりも物理的により正確な解が得られる。
• 数値結果から、無条件エネルギー安定性のスキームですら、時間刻みが大きすぎる場合には不正確であることが判明し、エネルギー安定性のみでは良好な数値性能を保証できないことが強調された。