[論文レビュー] A Variational Lagrangian Scheme for a Phase Field Model: A Discrete Energetic Variational Approach

Chun Liu, Yiwei Wang|arXiv (Cornell University)|Mar 23, 2020

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

ひとこと要約

本稿では、エネルギー散逸構造を離散的エネルギー変分的アプローチによって保存するように、修正されたAllen-Cahn相場モデルに対する変分的ラグランジュスキームを提案する。ラグランジュ写像に正則化項を組み込んだ動的記述により、エネルギー安定性が達成され、最小限のメッシュポイントで薄い拡散界面を正確に捉えることができ、平衡状態の堅牢な計算が可能になる。

ABSTRACT

In this paper, we propose a variational Lagrangian scheme for a modified phase-field model, which can compute the equilibrium states for the original Allen-Cahn type model. Our discretization is based on a prescribed energy-dissipation law in terms of the flow map. By employing a discrete energetic variational approach, this scheme preserves the variational structure of the original energy-dissipation law and is energy stable. Plentiful numerical tests show that, by choosing the initial value properly, our methods can produce the desired equilibrium and capture the thin diffuse interface with a small number of mesh points.

研究の動機と目的

• Allen-Cahn型相場モデルにおける平衡状態を計算する構造保存型数値スキームの開発を目的とする。
• Eulerian法が最小限の計算コストで薄い拡散界面を解像するのを妨げる限界を克服することを目的とする。
• ラグランジュ法を、ラグランジュ写像上に自然な変分構造を欠く$L^2$-勾配流れに拡張することを目的とする。
• エネルギー散逸法則における正則化項を通じてエネルギー安定性とメッシュ品質を保証することを目的とする。
• 本スキームが最小限のメッシュ解像度で鋭い界面を捉え、正しい平衡状態に収束することの有効性を示すこと

提案手法

• 本手法は、流れ写像に基づいて表現された所定のエネルギー散逸法に依拠しており、変分的構造の保存を保証する。
• 離散的エネルギー変分的手法が用いられ、連続的なエネルギー散逸法則が離散レベルでも維持される時間離散スキームが導出される。
• メッシュのねじれを防ぎ、安定性を確保するため、散逸に正則化項$\nu|\nabla\mathbf{u}|^2$が組み込まれている。
• 系はラグランジュ写像$\mathbf{X}(\mathbf{X}_0, t)$を用いて定式化され、相関関数$\varphi$は$\varphi_t + \nabla\varphi \cdot \mathbf{u} = 0$に従って進化する。
• 離散系は、一貫性のある質量行列$\mathbf{M}^e$と剛性行列$\mathbf{K}^e$を用いた有限要素法で解かれる。$\nu > 0$のとき、これらは正定値性を保証する。
• 本手法は自己適応的であり、メッシュの適応的 refinement を必要とせず、メッシュポイントを拡散界面に集中させる。

実験結果

リサーチクエスチョン

• RQ1 変分的ラグランジュスキームは、離散レベルにおいてAllen-Cahn型相場モデルのエネルギー散逸構造を保存できるか？
• RQ2 $L^2$-勾配流れはラグランジュ写像上に自然な変分構造を欠くため、ラグランジュ法はどのようにそれらに適応できるか？
• RQ3 提案されたスキームは、最小限のメッシュ解像度で正確な平衡状態を計算でき、メッシュのねじれを回避できるか？
• RQ4 正則化項$\nu|\nabla\mathbf{u}|^2$は、数値的安定性とメッシュ品質を確保するために果たす役割は何か？
• RQ5 初期条件$\varphi_0$の選択が、ラグランジュシミュレーションにおける正しい平衡状態への収束に与える影響は何か？

主な発見

• 提案されたスキームはエネルギー安定であり、離散的エネルギー変分的定式化を通じて連続的なエネルギー散逸法則の変分的構造を保存する。
• 適切に選ばれた初期条件のもとで、本手法はAllen-Cahnモデルの正しい平衡状態を的確に計算することができる。
• 本スキームは少数のメッシュポイントで薄い拡散界面を捉えることができ、高い解像度効率を示す。
• 数値的実験では、本手法がメッシュのねじれを回避し、時間の経過にかかわらず良好なメッシュ品質を維持する一方、単純なラグランジュ手法とは異なり、劣化しないことが示された。
• 正則化項$\nu|\nabla\mathbf{u}|^2$は、離散的散逸行列$\mathbf{D}$の正定値性を保証するために不可欠であり、これにより数値的崩壊を防げる。
• 不適切な初期条件$\varphi_0$を用いた失敗例では、$\varphi_0$が$[-1,1]$に有界でない場合、ラグランジュ法が正しい平衡状態に到達できないことが示され、初期データ選択の重要性が強調された。