[論文レビュー] On Linear and unconditionally energy stable Algorithms for Variable Mobility Cahn-Hilliard Type Equation with Logarithmic Flory-Huggins Potential

Xiaofeng Yang, Jia Zhao|arXiv (Cornell University)|Jan 25, 2017

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

ひとこと要約

本稿では、不変エネルギー二次化（IEQ）手法を用いて、対数的Flory-Hugginsポテンシャルを有する可変移動度Cahn-Hilliard方程式に対して、線形で、無条件エネルギー安定な時間積分スキームを提案する。非線形ポテンシャルを変数変換により二次形式に再定式化することで、非線形項を半陽的処理可能となり、無条件エネルギー安定で、反復解法を用いて効率的に解ける線形対称正定値系が得られる。本稿は、この困難な相場モデルクラスに対して、線形で高次かつ証明可能な無条件エネルギー安定スキームを達成した最初の研究である。

ABSTRACT

In this paper, we consider the numerical approximations for the fourth order Cahn-Hilliard equation with concentration dependent mobility, and the logarithmic Flory-Huggins potential. One challenge in solving such a diffusive system numerically is how to develop proper temporal discretization for nonlinear terms in order to preserve the energy stability at the time-discrete level. We resolve this issue by developing a set of the first and second order time marching schemes based on a novel, called "Invariant Energy Quadratization" approach. Its novelty is that the proposed scheme is linear and symmetric positive definite because all nonlinear terms are treated semi-explicitly. We further prove all proposed schemes are unconditionally energy stable rigorously. Various 2D and 3D numerical simulations are presented to demonstrate the stability, accuracy and efficiency of the proposed schemes thereafter.

研究の動機と目的

• ポリマーおよび相分離系に広く用いられる、対数的Flory-Hugginsポテンシャルを有する可変移動度Cahn-Hilliard方程式に対する、効率的で高精度な時間積分スキームの開発。
• 特に、薄い界面を有する剛性系において、反復的解法を要する非線形スキームや、小さな時間ステップを要する条件付き安定スキームの限界を克服すること。
• 各時間ステップで非線形性を保ちつつ、線形系を維持し、効率的かつ頑健な数値シミュレーションを可能にする無条件エネルギー安定性の維持。
• 熱力学的に整合性のあるモデルに一般的に見られる複雑な有界非線形ポテンシャル、例えば対数的Flory-Hugginsポテンシャルに対して、IEQ手法の適用範囲を拡張すること。
• 複雑な非線形自由エネルギー項を有する勾配流れ問題に対して、高次で線形かつ無条件安定なスキームを構築する一般化可能なフレームワークの提供。

提案手法

• 不変エネルギー二次化（IEQ）手法を用いて、新たな補助変数を導入することで、対数的Flory-Hugginsポテンシャルを二次形式に変換し、非線形項の半陽的処理を可能にする。
• 再定式化された系は元のエネルギー散逸法則を保持するため、離散レベルでもエネルギー安定性が保証される。
• Crank-Nicolson時間積分を用いて、1次および2次時間離散化スキーム（LS1-CNおよびLS2-CN）を構築し、非線形項はIEQ再定式化により半陰的に扱う。
• 得られる線形系は対称正定値となるため、質量巻き込みプリコンディショナを用いたKrylov部分空間法により効率的に解ける。
• 標準的なガレルキン有限要素法を用いて2次元および3次元空間離散化を実施し、変分式による一貫性を保証する。
• 理論的分析により、時間ステップサイズに依存しない、すべての提案スキームの無条件エネルギー安定性が証明された。

実験結果

リサーチクエスチョン

• RQ1 可変移動度Cahn-Hilliard方程式に、対数的Flory-Hugginsポテンシャルを組み合わせた場合に、線形で高次かつ無条件エネルギー安定な時間離散化スキームを構築可能か？
• RQ2 従来の凸分割法や安定化陽解法が失敗するか、計算的に非現実的になる状況において、IEQ手法が安定的で効率的かつ高精度な数値スキームを提供可能か？
• RQ3 IEQに基づくスキームは、非線形解法を回避しつつ、対称正定値線形系を保ち、無条件エネルギー安定性を維持可能か？
• RQ4 提案スキームは、2次元および3次元における相分離ダイナミクスの長時間シミュレーションにおいて、特に異なる初期濃度（φ₀）に対してどのように性能を発揮するか？
• RQ5 IEQフレームワークは、複雑で有界な非線形ポテンシャルを有する他の勾配流れモデルへ一般化可能か？

主な発見

• 提案されたLS1-CNおよびLS2-CNスキームは、変分式に基づくエネルギー安定性解析により、厳密に無条件エネルギー安定であることが証明された。
• スキームは線形で、対称正定値であるため、質量巻き込みプリコンディショナを用いたKrylov部分空間ソルバーにより、効率的に解ける。
• 2次元および3次元における数値シミュレーションは、長時間安定性を示し、自由エネルギーが時間とともに単調に減少する結果が得られ、エネルギー散逸が確認された。
• 3次元における相分離ダイナミクスでは、φ₀ = 0.3およびφ₀ = 0.7の場合は球状ドロップレットの形成が観察され、φ₀ = 0.5の場合は二層構造が形成され、熱力学的期待と整合的であった。
• 2次スキームであるLS2-CNはδt = 0.001で高い精度を達成し、長時間積分においても不安定性や発散が観察されなかった。
• 本研究は、対数的Flory-Hugginsポテンシャルを有する可変移動度Cahn-Hilliard系に対して、線形で高次かつ無条件エネルギー安定なスキームを初めて達成した。これにより、文献における重要な空白が埋められた。