[論文レビュー] Efficient schemes with unconditionally energy stability for the anisotropic Cahn-Hilliard Equation using the stabilized-Scalar Augmented Variable (S-SAV) approach

Xiaofeng Yang|arXiv (Cornell University)|Apr 8, 2018

Solidification and crystal growth phenomena参考文献 12被引用数 10

ひとこと要約

本稿では、安定化スカラー拡張変数（S-SAV）アプローチを用いて、異方性Cahn-Hilliard方程式の2つの2次精度で、無条件エネルギー安定な時間積分スキームを提案する。S-AV法と3つの線形安定化項を組み合わせることで、非線形項を半陽的かつ明示的に取り扱い、各時間ステップで定数係数を持つ3つの分離可能な線形方程式を解くのみで、強い異方性に起因する振動を抑えるとともに、大きな時間ステップを許容する。本手法は、2次精度を保ちつつ、定数係数線形方程式のみを解くという高効率性を兼ね備えた、このモデルに対しては初の厳密に無条件エネルギー安定なスキームである。

ABSTRACT

In this paper, we consider numerical approximations for the anisotropic Cahn-Hilliard equation. The main challenge of constructing numerical schemes with unconditional energy stabilities for this model is how to design proper temporal discretizations for the nonlinear terms with the strong anisotropy. We propose two, second order time marching schemes by combining the recently developed SAV approach with the linear stabilization approach, where three linear stabilization terms are added. These terms are shown to be crucial to remove the oscillations caused by the anisotropic coefficients, numerically. The novelty of the proposed schemes is that all nonlinear terms can be treated semi-explicitly, and one only needs to solve three decoupled linear equations with constant coefficients at each time step. We further prove the unconditional energy stabilities rigorously, and present various 2D and 3D numerical simulations to demonstrate the stability and accuracy.

研究の動機と目的

• 強い異方性に起因する数値的振動のため困難である、異方性Cahn-Hilliard方程式に対する効率的かつ無条件エネルギー安定な数値スキームの開発。
• S-AV法と線形安定化を組み合わせることで、従来のスキームに見られるエネルギー安定性の証明不能性や高コストの問題を克服する。
• 各時間ステップで定数係数線形方程式のみを解くことで、時間精度2次を維持しつつ計算効率を確保する。
• 強い異方性下でも不安定化を伴わずに大きな時間ステップを許容できることを示す。
• 提案スキームの無条件エネルギー安定性を、数学的証明により厳密に保証する。これは、先行手法と比較しての重要な前進である。

提案手法

• エネルギー汎関数を新たなスカラー補助変数を用いて二次形式に再定式化するS-AV（Scalar-Augmented-Variable）法を採用し、非線形項の半陽的取り扱いを可能にする。
• 異方性係数γ(∇φ/|∇φ|)に起因する高周波数振動を抑えるために、3つの線形安定化項（S1, S2, S3）を導入する。特に異方性が強い場合に有効である。
• エネルギー散逸則を無条件に保持する2次精度BDF2型時間積分スキーム（Schemes BDF2-S および BDF2-W）を構築する。
• 各時間ステップで、定数係数を持つ3つの独立した線形楕円型方程式系に分離可能であり、標準的な線形ソルバーで解ける。
• 試験関数を試行関数と同じ関数空間にとることで変分式を導出し、エネルギー推定による安定性を保証する。
• 数値的条件数の改善と振動の低減、特に3次元計算において有効なため、Willmore正則化項（βG(φ)）を導入する。

実験結果

リサーチクエスチョン

• RQ1 強い異方性下でも、無条件エネルギー安定性を保証する2次精度時間積分スキームを構築することは可能か？
• RQ2 異方性係数γ(∇φ/|∇φ|)を含む非線形項を、効率的かつ安定的に半陽的かつ明示的に取り扱う方法は何か？
• RQ3 追加の線形安定化項が、異方性係数に起因する数値的振動を抑制する役割を果たすメカニズムは何か？
• RQ4 提案スキームは、強い異方性下でも大きな時間ステップを許容しながら、高い精度と安定性を維持できるか？
• RQ5 各時間ステップで定数係数線形方程式のみを解くという制約のもとで、厳密な無条件エネルギー安定性を達成することは可能か？

主な発見

• 提案されたS-SAVスキームは、変分式に基づくエネルギー推定により、厳密に無条件エネルギー安定であることが証明された。
• 数値実験の結果、S1 = S2 = 2 かつ S3 = 0.001 の組み合わせが高周波数振動を効果的に抑制する一方、S1 = S2 = 0 かつ S3 = 0.001 の組み合わせでは物理的に不自然な振動が生じる。
• δt = 1e−4 の大きな時間ステップでも精度を維持可能であり、δt = 1e−2, 1e−3, 1e−4 の結果は、ほぼ同一の相場分部を示す。
• 異方性パラメータα（0.1から0.3に増加）を大きくすると、平衡状態の形状におけるピラミッドの角が鋭くなる傾向を示し、物理的期待と整合的である。
• Willmore正則化パラメータβを小さくすると、平衡形状の角が鋭くなる傾向が確認され、時間発展プロットからエネルギーの減少がより速くなることも確認された。
• 本手法は一般化可能であり、他の複雑な非線形性を有する勾配流れ問題へも応用可能である。証明フレームワークは、整合的なガラーキン近似に基づく変分式に依存するため。