[論文レビュー] A variational analysis for the moving finite element method for gradient flows

Xianmin Xu|arXiv (Cornell University)|Sep 3, 2020

Advanced Numerical Methods in Computational Mathematics参考文献 29被引用数 4

ひとこと要約

本稿では、エネルギー散逸を保証するOnsager原理を基盤として、勾配流れ方程式に対する移動有限要素法（MFEM）の新規な変分的導出を提示する。MFEMがエネルギー汎関数の局所最小解に収束することを証明し、グローバル最小解が捉えられた場合には最適収束率を示す。線形および非線形拡散問題、特にAllen-Cahn方程式について数値的に検証された。

ABSTRACT

By using the Onsager principle as an approximation tool, we give a novel derivation for the moving finite element method for gradient flow equations. We show that the discretized problem has the same energy dissipation structure as the continuous one. This enables us to do numerical analysis for the stationary solution of a nonlinear reaction diffusion equation using the approximation theory of free-knot piecewise polynomials. We show that under certain conditions the solution obtained by the moving finite element method converges to a local minimizer of the total energy when time goes to infinity. The global minimizer, once it is detected by the discrete scheme, approximates the continuous stationary solution in optimal order. Numerical examples for a linear diffusion equation and a nonlinear Allen-Cahn equation are given to verify the analytical results.

研究の動機と目的

• 勾配流れ系に対する移動有限要素法（MFEM）の新規な変分的導出を、Onsager原理を基盤とするフレームワークとして開発すること。
• 離散MFEM問題が連続勾配流れ系と同一のエネルギー散逸構造を保持することを確立すること。
• 自由ノットを持つ区分的多項式近似を用いて、非線形反応拡散方程式の定常解に対する厳密な誤差解析を実施すること。
• 離散スキームがエネルギー汎関数のグローバル最小解を検出する条件の下で、最適収束率を証明すること。
• 線形および非線形拡散方程式（Allen-Cahn方程式を含む）に対する理論的結果の数値的検証を実施すること。

提案手法

• 散逸とエネルギー変化の項からなるレイリーアン関数を最小化することで、時間発展をOnsager変分原理を用いて導出する。
• 自由ノットを持つ区分的多項式空間における有限要素近似を定式化し、節点値とメッシュ頂点の両方を最適化すべき未知数として扱う。
• 節点値とメッシュ点の時間発展を記述する常微分方程式（ODE）系を導出する。この系は、滑らか化処理や特異関数の取り扱いを要せず、古典的MFEM形式と一致する。
• 連続問題におけるエネルギー散逸構造が離散形式でも保たれることで、安定性および収束解析が可能になる。
• 自由ノットを持つ区分的多項式の非線形近似理論を用いて誤差評価を導出し、特に定常解に対して有効である。
• 数値実験では、1ステップあたりのエネルギー減少量が許容誤差以下に下がるまで離散系を解き、メッシュの適応的 refinement が解の曲率が大きい領域にノードが集中することを観察した。

実験結果

リサーチクエスチョン

• RQ1 Onsager変分原理は、勾配流れ系に対する移動有限要素法を体系的に導出するためのフレームワークとして利用可能か？
• RQ2 離散MFEM形式は、連続勾配流れ系のエネルギー散逸構造を保持するか？
• RQ3 時間無限大におけるMFEMがエネルギー汎関数の局所最小解に収束する条件は何か？
• RQ4 離散空間でグローバル最小解が検出された場合、定常解に対するMFEMの最適収束率は何か？
• RQ5 非線形勾配流れにおけるMFEMにおいて、H¹ノルムの誤差とエネルギー誤差は、メッシュの細分化に伴いどのようにスケーリングされるか？

主な発見

• Onsager原理を用いて導出したMFEMは、連続勾配流れ系のエネルギー散逸構造を正確に保持する。
• 適切な条件下では、MFEM解は時間無限大に近づくにつれてエネルギー汎関数の局所最小解に収束する。
• 離散自由ノット区分的線形空間でグローバル最小解が検出された場合、H¹誤差は要素数Nに関して最適収束率O(N⁻¹)で収束する。
• エネルギー誤差は最適収束率O(N⁻²)で収束し、非線形近似理論から導かれた理論的境界と整合的である。
• ε = 0.05およびε = 0.01のAllen-Cahn方程式に対する数値結果は、最適収束率を確認した：H¹誤差 ≈ O(N⁻¹)、エネルギー誤差 ≈ O(N⁻²)、N = 5の粗いメッシュでも成立。
• 適応的メッシュは、急峻な遷移プロファイルの内部層に自然にノードを集中させ、H¹ノルムにおける誤差が一様に分布していることを示している。