[論文レビュー] $L^2$ norm error estimates of BDF methods up to fifth-order for the phase field crystal model

Hong-lin Liao, Yuanyuan Kang|arXiv (Cornell University)|Jan 3, 2022

Solidification and crystal growth phenomena被引用数 4

ひとこと要約

本稿は、6次非線形放物型方程式である相場水晶モデルに高次BDF法（3〜5次）を適用した際の最初の$L^2$ノルム誤差推定を確立する。著者らは、新規の離散勾配構造を構築し、直交畳み込みカーネルと離散エネルギー技法を活用することで、非A安定なこれらのスキームにおいても離散レベルでのエネルギー散逸を証明し、最適な誤差境界を導出する。

ABSTRACT

The well-known backward difference formulas (BDF) of the third, the fourth and the fifth orders are investigated for time integration of the phase field crystal model. By building up novel discrete gradient structures of the BDF-$ mk$ ($ mk=3,4,5$) formulas, we establish the energy dissipation laws at the discrete levels and then obtain the priori solution estimates for the associated numerical schemes (however, we can not build any discrete energy dissipation law for the corresponding BDF-6 scheme because the BDF-6 formula itself does not have any discrete gradient structures). With the help of the discrete orthogonal convolution kernels and Young-type convolution inequalities, some concise $L^2$ norm error estimates (with respect to the starting data in the $L^2$ norm) are established via the discrete energy technique. To the best of our knowledge, this is the first time such type $L^2$ norm error estimates of non-A-stable BDF schemes are obtained for nonlinear parabolic equations. Numerical examples are presented to verify and support the theoretical analysis.

研究の動機と目的

• 相場水晶（PFC）モデルの長時間シミュレーションに向けた高次時間積分スキームの開発。
• 非A安定であるが数値的に散逸的であるBDF-3、BDF-4、BDF-5法に対する厳密な$L^2$ノルム誤差推定の確立。
• BDF-6に離散エネルギー散逸則が欠如しているという課題を克服し、高次BDFスキームの限界を明確化すること。
• 勾配流れ問題の長時間粗大化ダイナミクスのシミュレーションにおける高次BDF法の使用の理論的基盤を提供すること。
• 非A安定時間積分法を用いる非線形高次放物型方程式にまで、離散エネルギー解析技法を拡張すること。

提案手法

• BDF-3、BDF-4、BDF-5式のための新規の離散勾配構造を提案し、離散エネルギー散逸則を確立する。
• 時間における誤差伝播を制御するため、直交畳み込みカーネルとヤング型畳み込み不等式を用いる。
• 離散エネルギー技法を適用し、初期データが$L^2$ノルムにある場合の$L^2$ノルム誤差推定を導出する。
• BDFカーネルの再帰的関数的分解を構築し、離散エネルギー汎関数の正定値性を証明する。
• 和分部分積分恒等式および離散テイラー展開を用いて、スキームの整合性と安定性を分析する。
• BDF-6が離散勾配構造を持たないことを示し、その結果、エネルギー散逸則の導出および同様の誤差解析フレームワークの構築が不可能であることを示す。

実験結果

リサーチクエスチョン

• RQ1 非A安定なBDFスキーム（3〜5次）を相場水晶モデルに適用した際、$L^2$ノルム誤差推定を確立できるか。
• RQ2 BDF-3、BDF-4、BDF-5スキームは、新規の離散勾配構造を通じて離散エネルギー散逸則を保持するか。
• RQ3 なぜBDF-6スキームは離散エネルギー散逸則を有しないのか。その影響は誤差解析にどのような意味を持つのか。
• RQ4 直交畳み込みカーネルと離散エネルギー技法を効果的に組み合わせることで、非線形放物型方程式に対する最適な誤差境界を導出できるか。
• RQ5 初期データの$L^2$ノルムが、高次BDFスキームの収束速度に果たす役割は何か。

主な発見

• 著者らは、6次非線形放物型方程式である相場水晶モデルにBDF-3、BDF-4、BDF-5スキームを適用した際の最初の$L^2$ノルム誤差推定を確立する。
• BDF-3からBDF-5に対して、新規の離散勾配構造を用いて離散エネルギー散逸則を効果的に構築し、安定性と収束性を保証する。
• BDF-6スキームは離散勾配構造を持たないため、離散エネルギー散逸則の導出が不可能であり、結果として同様の誤差解析フレームワークは成立しない。
• 直交畳み込みカーネルと離散エネルギー技法を組み合わせることで、最適な$L^2$誤差境界が導出され、誤差推定は初期データの$L^2$ノルムに依存する。
• 定数$\sigma_{L5} = \frac{646631}{1920000}$がBDF-5のエネルギー汎関数に明示的に特定され、離散エネルギー汎関数の正定値性を保証する。
• 数値例は理論的収束速度を確認し、導出された誤差推定の鋭さを検証する。