[論文レビュー] Stability analysis for the Implicit-Explicit discretization of the Cahn-Hilliard equation

Dong Li, Tao Tang|arXiv (Cornell University)|Aug 9, 2020

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

ひとこと要約

本稿は、安定化項を追加せず、問題構造を変更せず、Cahn-Hilliard方程式の陰-陽的（IMEX）時間離散化について、きめ細やかな安定性解析を提示する。エネルギーと粘性に自然にスケーリングする最適な時間刻み制約のもとで、条件付きエネルギー安定性を確立し、最大原理の欠如と小さな移動度係数による剛性に対処するため、L∞ 範囲を制御するための新規な「エネルギーとL∞ をトレードする（TEFL）」フレームワークを導入することで、相場分離シミュレーションにおける長年の未解決問題を解決する。

ABSTRACT

Implicit-Explicit methods have been widely used for the efficient numerical simulation of phase field problems such as the Cahn-Hilliard equation or thin film type equations. Due to the lack of maximum principle and stiffness caused by the effect of small dissipation coefficient, most existing theoretical analysis relies on adding additional stabilization terms, mollifying the nonlinearity or introducing auxiliary variables which implicitly either changes the structure of the problem or trades accuracy for stability in a subtle way. In this work we introduce a robust theoretical framework to analyze directly the stability of the standard implicit-explicit approach without stabilization or any other modification. We take the Cahn-Hilliard equation as a model case and prove energy stability under natural time step constraints which are optimal with respect to energy scaling. These settle several questions which have been open since the work of Chen and Shen \cite{CS98}.

研究の動機と目的

• 標準的なIMEXスキームをCahn-Hilliard方程式に適用した場合の安定性を、きめ細やかに証明するという長年の課題を解決すること。
• 理論的解析において、人工的安定化、平滑化、または補助変数に依存しないようにすること。
• 小さな移動度係数の文脈において、エネルギースケーリングに最適な時間刻み制約を導出すること。
• 自然な条件下で直接安定性解析が可能な、強固な理論的フレームワーク（TEFL）を確立すること。
• 長時間・大規模な相場分離シミュレーションにおけるIMEX法の頑健性を確認すること。

提案手法

• 著者らは、方程式を変更せず、安定化を追加せず、解のL∞ 範囲を制御するための新規な「エネルギーとL∞ をトレードする（TEFL）」フレームワークを考案した。
• ChenとShen（2002）が提唱した標準的なIMEXスキームを分析し、非線形項を陽的に、線形拡散項を陰的に取り扱った。
• Cahn-Hilliard方程式の勾配流れ構造を活用した離散的エネルギー散逸の議論により、エネルギー安定性を証明した。
• 主な推論では、L2 および H1 範囲を用いて解のL∞ 範囲を有界化し、スケーリング解析による時間刻みの厳密な制御を実施した。
• エネルギーの減衰とL∞ の増大をバランスさせるために、エネルギースケーリングに最適な時間刻み制約を導出した。
• 2次元および3次元の周期的領域において、フーリエに基づく技術と時間ステップに関する帰納法を用いて解析を実施した。

実験結果

リサーチクエスチョン

• RQ1 安定化を追加せずに、標準的なIMEXスキームがCahn-Hilliard方程式に対してきめ細やかにエネルギー安定性を示せるか？
• RQ2 小さな移動度係数が存在する状況下で、エネルギー安定性のための最適な時間刻み制約は何か？
• RQ3 最大原理の欠如と、小さな粘性係数領域における剛性は、理論的解析でどのように扱えるか？
• RQ4 補助変数を導入したり非線形項を変更したりせずに、直接的な解析が可能か？
• RQ5 TEFLフレームワークは、同様の剛性と非線形性を示す他の相場分離モデルへ一般化可能か？

主な発見

• 本稿は、エネルギースケーリングに最適な時間刻み制約のもとで、IMEXスキームのCahn-Hilliard方程式に対する条件付きエネルギー安定性を確立した。
• 時間刻み制約は τ ≤ 0.0007 min{β₃⁻⁸, β₄⁻⁸ᐟ³, β₅⁻⁸} として導出され、β₃, β₄, β₅ は初期エネルギーとL² 範囲に依存する。
• 安定化、平滑化、または補助変数を回避することで、直接的かつ構造を保った証明が可能となった。
• TEFLフレームワークは、エネルギーおよびH¹ 範囲を用いて解のL∞ 範囲を効果的に制御し、緩い時間刻み制限のもとで安定性を実現した。
• 導出された時間刻み制約は、エネルギースケーリングの観点からほぼ最適であり、相分離ダイナミクスの物理的期待と一致する。
• 本手法は一般化可能であり、高次のIMEXスキームや他の相場分離モデルに対する洗練された安定性解析の道筋を示した。