[論文レビュー] Weighted and shifted BDF2 methods on variable grids

Minghua Chen, Yu Fan+1 more|arXiv (Cornell University)|Aug 6, 2021

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

ひとこと要約

本稿は、非一様時間グリッド上での放物型PDEを解くために、可変時間ステップを用いた重み付き・シフト付きBDF2（WSBDF2）法を導入する。BDF2とCrank-Nicolsonスキームを重み付きシフトによって統合し、無条件安定性および最適な2次収束性を証明した。前例のBDF2法よりも著しく向上した最適な隣接時間ステップ比を導出し、可変グリッドにおけるPDE解析における理論的空白を埋めた。

ABSTRACT

Variable steps implicit-explicit multistep methods for PDEs have been presented in [17], where the zero-stability is studied for ODEs; however, the stability analysis still remains an open question for PDEs. Based on the idea of linear multistep methods, we present a simple weighted and shifted BDF2 methods with variable steps for the parabolic problems, which serve as a bridge between BDF2 and Crank-Nicolson scheme. The contributions of this paper are as follows: we first prove that the optimal adjacent time-step ratios for the weighted and shifted BDF2, which greatly improve the maximum time-step ratios for BDF2 in [11,15]. Moreover, the unconditional stability and optimal convergence are rigorous proved, which make up for the vacancy of the theory for PDEs in [17]. Finally, numerical experiments are given to illustrate theoretical results.

研究の動機と目的

• PDEにおける可変ステップBDF法の安定性および収束理論の欠如に応えること、特に陰・陽マルチステップスキームに関して。
• 非一様時間グリッド上での重み付き・シフト付きBDF2の定式化を通じて、BDF2とCrank-Nicolsonスキームを統一的フレームワークで結びつけること。
• 安定性と効率性を最大化する最適な隣接時間ステップ比を導出し、従来のBDF2法における境界を改善すること。
• 放物問題にWSBDF2法を適用した場合の無条件安定性および最適な2次収束性を厳密に証明すること。
• スペクトル的空間離散化と可変時間ステップを用いた数値実験を通じて理論的考察を検証すること。

提案手法

• 重み付き・シフト付きBDF2作用素 $ D_2v^n = \theta\tilde{D}_2v^n + (1-\theta)\bar{D}_2v^n $ を定式化し、$ \tilde{D}_2 $ は標準BDF2、$ \bar{D}_2 $ はシフト付きBDF2であり、$ \theta \in [1/2,1] $ とする。
• 時間離散スキームを $ D_2u^n - \theta\Delta u^n - (1-\theta)\Delta u^{n-1} = \theta f^n + (1-\theta)f^{n-1} $ として定義し、可変時間ステップを用いて放物型問題を解く。
• $ L^2 $内積と $ H_0^1 $ テスト関数を用いた弱形式を導出し、安定性および収束推定値を得る。
• エネルギー推定と離散畳み込み恒等式を用いて誤差を評価し、コーシー＝シュワルツの不等式および和分部分積分技術を活用する。
• 誤差伝搬を制御するための離散グローワル不等式を適用し、最適な収束次数が得られる。
• スペクトルコロケーション法を用い、チェビシェフ＝ガウス＝ロバット点を用いて空間離散化を行い、スペクトル的精度と無視できる空間誤差を確保する。

実験結果

リサーチクエスチョン

• RQ1 可変グリッド上でのWSBDF2法の安定性および最適収束性を保証するための最適な隣接時間ステップ比 $ r_k = \tau_k / \tau_{k-1} $ の範囲は何か？
• RQ2 可変時間ステップを用いた放物型PDEに対して、WSBDF2法が無条件安定性および2次収束性を達成できるか？ これは、先行研究における理論的欠落を埋めるものである。
• RQ3 可変時間ステッピングの文脈において、WSBDF2法はBDF2とCrank-Nicolsonスキームの安定性および精度特性をどのように統合するか？
• RQ4 幾何的およびランダムなステップサイズを含む可変時間ステップ戦略が、WSBDF2法の収束挙動に及ぼす影響はどの程度か？
• RQ5 WSBDF2フレームワークを高次法（例：WSBDF3）に非一様グリッド上で拡張可能か？ その理論的課題は何か？

主な発見

• WSBDF2法はすべての $ \theta \in [1/2,1] $ に対して無条件安定性を達成し、$ \theta = 1/2 $ の場合に最適な時間ステップ比 $ r_k \leq 4 $ が得られ、標準BDF2法よりも著しく高い許容可能比を実現した。
• すべての $ \theta \in [1/2,1] $ に対して、離散 $ L^2 $-ノルムにおける最適な2次収束性が厳密に証明され、理論的期待と一致する収束率が得られた。
• 数値実験により、すべてのテストケースで2次収束が確認された：ケースI（$ r_k = 4, 1/4 $ の交互）、ケースII（$ r_k = 2 $）、ケースIII（ランダム $ r_k $）、観測収束率は2.0に近い。
• ケースII（$ r_k = 2 $ の定数）では収束率が2次のままであり、非最適なステップ比でも安定性が保たれていることが示された。
• ケースIII（ランダム時間ステップ）では収束率が約2.0～2.5を維持し、不規則な時間グリッドに対してもロバストであることが示された。
• BDF2（$ \theta = 1 $）とCrank-Nicolson（$ \theta = 1/2 $）の両方を効果的に統合しており、特に $ \theta = 1/2 $ の場合、すべてのテストケースで一貫した収束が得られた。