[論文レビュー] Long time error analysis of the fourth-order compact finite difference methods for the nonlinear Klein-Gordon equation with weak nonlinearity

Yue Feng|arXiv (Cornell University)|Mar 9, 2020

Numerical methods for differential equations参考文献 42被引用数 4

ひとこと要約

本稿では、弱い非線形性 $\varepsilon^p$ を持つ非線形クライン=ゴルドン方程式に対して、第四次精度のコンパクト有限差分法（4cFD）を提案する。長時間 $O(\varepsilon^{-p})$ における厳密な誤差境界が得られ、最適な精度を達成するには $h = O(\varepsilon^{p/4})$ および $\tau = O(\varepsilon^{p/2})$ が必要であることが示され、第二階精度法と比較して優れた空間分解能を示し、数値的に鋭い収束率が確認された。

ABSTRACT

We present the fourth-order compact finite difference (4cFD) discretizations for the long time dynamics of the nonlinear Klein-Gordon equation (NKGE), while the nonlinearity strength is characterized by $\varepsilon^p$ with a constant $p \in \mathbb{N}^+$ and a dimensionless parameter $\varepsilon \in (0, 1]$. Based on analytical results of the life-span of the solution, rigorous error bounds of the 4cFD methods are carried out up to the time at $O(\varepsilon^{-p})$. We pay particular attention to how error bounds depend explicitly on the mesh size $h$ and time step $τ$ as well as the small parameter $\varepsilon \in (0, 1]$, which indicate that, in order to obtain `correct' numerical solutions up to the time at $O(\varepsilon^{-p})$, the $\varepsilon$-scalability (or meshing strategy requirement) of the 4cFD methods should be taken as: $h = O(\varepsilon^{p/4})$ and $τ= O(\varepsilon^{p/2})$. It has better spatial resolution capacity than the classical second order central difference methods. By a rescaling in time, it is equivalent to an oscillatory NKGE whose solution propagates waves with wavelength at $O(1)$ in space and $O(\varepsilon^p)$ in time. It is straightforward to get the error bounds of the oscillatory NKGE in the fixed time. Finally, numerical results are provided to confirm our theoretical analysis.

研究の動機と目的

• 弱い非線形性 $\varepsilon^p$ を特徴とする非線形クライン=ゴルドン方程式（NKGE）に第四次精度コンパクト有限差分法（4cFD）を適用した場合の長時間誤差境界を発展および解析すること。
• 長時間 $O(\varepsilon^{-p})$ にわたり精度を維持するために必要なメッシュ設定（$h$ と $\tau$）を特定すること。
• 4cFD法の $\varepsilon$ スケーラビリティを確立し、古典的な第二階精度中央差分法と比較して優れた空間分解能を示すこと。
• 時間スケーリングを用いてNKGEの振動的形に一般化し、固定時間誤差境界を導出可能にする。
• 理論的誤差推定値の数値的妥当性を検証し、その鋭さを確認すること。

提案手法

• 弱い非線形性 $\varepsilon^p u^{p+1}$ を持つNKGEに4cFD法を適用し、第四次精度の空間ステンシルと第二階精度の時間離散化を用いる。
• 元のNKGEを $O(1)$ の空間波長と $O(\varepsilon^p)$ の時間波長を持つ振動的形に変換するための時間スケーリング変換を適用する。
• 4cFDスキームの長時間 $O(\varepsilon^{-p})$ における厳密な誤差境界を、$h$、$\tau$、$\varepsilon$ に依存する形で導出する。
• エネルギー法と漸近展開を用いて、長時間にわたる誤差の増大を制御する。
• この手法を振動的NKGEに拡張し、固定時間誤差境界を直接導出可能にする。
• 数値的妥当性評価には、高精度の指数波動積分子フーリエ擬スペクトル法を基準解として用い、空間誤差と時間誤差を計算する。

実験結果

リサーチクエスチョン

• RQ1 弱い非線形性を有する非線形クライン=ゴルドン方程式に対して、第四次精度コンパクト有限差分法の長時間誤差挙動はいかなるものか？
• RQ2 長時間 $O(\varepsilon^{-p})$ にわたり精度を維持するためのメッシュサイズ $h$ と時間ステップ $\tau$ は何か？
• RQ3 本問題において、4cFD法は古典的な第二階精度中央差分法と比較して、空間分解能でどのように優れているか？
• RQ4 時間スケーリングを用いて得られる振動的NKGE形に、誤差境界を拡張できるか？
• RQ5 数値実験により、理論的誤差推定値が鋭いかどうかを確認できるか？

主な発見

• 固定された $\varepsilon > 0$ に対して、4cFD法は空間で第四次収束、時間で第二階収束を達成し、設計された精度順序が裏付けられた。
• 非常に高周波数の状態（$\varepsilon \ll 1$）では、$h \lesssim \varepsilon^{p/4}$ および $\tau \lesssim \varepsilon^{p/2}$ の条件下でのみ最適収束が観察され、理論的 $\varepsilon$ スケーリングの妥当性が裏付けられた。
• 必要なメッシュサイズ $h = O(\varepsilon^{p/4})$ は、同程度の精度を達成するための古典的第二階精度中央差分法（$h = O(\varepsilon^{p/2})$ 必要）と比較して、優れた空間分解能を提供する。
• 数値結果により、誤差境界が鋭いことが確認され、異なる $\varepsilon$ 値において理論的予測と一致する収束率が観察された。
• 振動的NKGEの定式化により、固定時間誤差境界を直接導出可能となり、長時間ダイナミクスの解析が簡素化された。
• $p=1$ および $p=2$ の数値実験により、導出されたメッシュ戦略に従うと、4cFD法が $O(\varepsilon^{-p})$ の長時間にわたり高い精度を維持することが確認された。