[論文レビュー] An implicit--explicit second order BDF numerical scheme with variable steps for gradient flows

Dianming Hou, Zhonghua Qiao|arXiv (Cornell University)|Apr 1, 2022

Advanced Numerical Methods in Computational Mathematics被引用数 4

ひとこと要約

本稿では、勾配流れにための可変時間刻みを用いた、スカラーアシスタント変数（SAV）法を用いた、新しい暗黙・陽解法（IMEX）2次BDFスキームを提案する。このスキームは、時間刻み比 ≤ 4.8645 の場合、無条件エネルギー安定性を達成し、非一様メッシュ上での最適誤差推定を伴う厳密な収束解析を提供する。数値実験および粗粒化ダイナミクスのシミュレーションにおける自己適合時間刻み法の効率性によって検証された。

ABSTRACT

In this paper, we propose and analyze an efficient implicit--explicit (IMEX) second order in time backward differentiation formulation (BDF2) scheme with variable time steps for gradient flow problems using the scalar auxiliary variable (SAV) approach. We prove the unconditional energy stability of the scheme for a modified discrete energy with the adjacent time step ratio $γ_{n+1}:=\Dt_{n+1}/\Dt_{n}\leq 4.8645$. The uniform $H^{2}$ bound for the numerical solution is derived under a mild regularity restriction on the initial condition, that is $ϕ(\x,0)\in H^{2}$. Based on this uniform bound, a rigorous error estimate of the numerical solution is carried out on the temporal nonuniform mesh. Finally, serval numerical tests are provided to validate the theoretical claims. With the application of an adaptive time-stepping strategy, the efficiency of our proposed scheme can be clearly observed in the coarsening dynamics simulation.

研究の動機と目的

• 非一様時間メッシュ上での効率的かつ無条件エネルギー安定な数値スキームの開発。
• 特に剛性および長時間シミュレーションにおいて、可変時間刻みを用いた高次・安定なBDF2スキームを構築する課題に対処する。
• 長期間のシミュレーションにおいて、安定性や精度を損なわずに自己適合時間刻み戦略を可能にする。
• やや弱い正則性仮定の下で、提案スキームの厳密な収束解析とエネルギー安定性を確立する。
• 粗粒化ダイナミクスにおける自己適合時間刻みを用いた数値実験を通じて、スキームの効率性と頑健性を示す。

提案手法

• スカラーアシスタント変数（SAV）アプローチを用いて、可変時間刻み（VBDF2）を用いた勾配流れの2次BDFスキームを定式化する。
• IMEX（暗黙・陽解法）時間離散化を導入：線形拡散項を暗黙的に、非線形ポテンシャル項を2次補間作用素を用いて陽的に扱う。
• 隣接する時間刻み比 $\gamma_{n+1} = \Delta t_{n+1}/\Delta t_n \leq 4.8645$ の条件下で、修正された離散エネルギー関数を用いて無条件エネルギー安定性を証明する。
• 時間の中間点 $t_{n+\sigma}$ における時間微分の2次精度近似を構築するために、2次補間作用素 $\Pi_{2,j}\phi(t)$ を用いる。
• SAV構造を保持するため、補助変数 $r^{n+1}$ および $\xi^{n+1}$ を定義し、$r$ と $\xi$ が1次精度であっても、位相変数 $\phi$ における2次精度を保証する。
• 初期データのやや弱い正則性条件 $\phi(\mathbf{x},0) \in H^2$ の下で、非一様メッシュ上での数値解に対する一様な $H^2$ 界を確立し、厳密な誤差解析を可能にする。

実験結果

リサーチクエスチョン

• RQ1 可変時間刻みを用いた非一様時間メッシュ上での、勾配流れに対する無条件エネルギー安定な2次BDFスキームを構築することは可能か？
• RQ2 SAV法を用いた可変ステップBDF2スキームにおいて、無条件エネルギー安定性を保つために許容可能な最大時間刻み比 $\gamma_{n+1}$ は何か？
• RQ3 SAVフレームワークを可変時間刻みを用いたIMEX BDF2スキームに拡張することで、2次精度と安定性を保ちながら、どのように実現できるか？
• RQ4 提案スキームの収束速度は非一様時間メッシュ上ですでにどれほどであり、初期データの正則性にどのように依存するか？
• RQ5 粗粒化ダイナミクスの長時間シミュレーションにおいて、自己適合時間刻み戦略を効果的にサポートできるか？

主な発見

• 提案されたIMEX VBDF2スキームは、時間刻み比 $\gamma_{n+1} \leq 4.8645$ の場合、無条件エネルギー安定である。これは、既存の文献における境界を著しく上回る。
• 初期条件 $\phi(\mathbf{x},0) \in H^2$ のやや弱い仮定の下で、数値解に対する一様な $H^2$ 界が確立され、収束解析に不可欠である。
• 非一様時間メッシュ上での数値解に対して、最適なオーダーの厳密な誤差推定が導出され、時間方向の2次収束が確認された。
• 補助変数 $r^{n+1}$ および $\xi^{n+1}$ が1次精度であっても、SAV定式化の構造のおかげで、位相変数 $\phi$ における2次精度が維持される。
• 数値実験により理論的主張が確認され、予測された収束率と良好に一致し、自己適合時間刻みを用いた場合に顕著な効率向上が示された。
• 粗粒化ダイナミクスのシミュレーションにおける自己適合時間刻みの応用により、長時間シミュレーションにおけるスキームの実用的効率性と頑健性が確認された。