[論文レビュー] Mesh-robustness of the variable steps BDF2 method for the Cahn-Hilliard model.

Hong-lin Liao, Bingquan Ji+2 more|arXiv (Cornell University)|Feb 7, 2021

Solidification and crystal growth phenomena参考文献 27被引用数 4

ひとこと要約

本稿では、ユーザーが定義するステップ比制限 $ r_{\text{user}} < 4.864 $ の下でエネルギー散逸とメッシュに頑健な安定性を維持する可変ステップ幅BDF2法をCahn-Hilliardモデルに対して提案する。離散直交畳み込みカーネルと新規の埋め込み不等式を活用することで、時間刻みの変動に依存しない $ L^2 $ ノルムの誤差推定が達成され、メッシュに頑健な収束性が保証され、適応的時間刻みに基づく効率的な長時間シミュレーションが可能になる。

ABSTRACT

The two-step backward differential formula (BDF2) implicit method with unequal time-steps is investigated for the Cahn-Hilliard model by focusing on the numerical influences of time-step variations. The suggested method is proved to preserve a modified energy dissipation law at the discrete levels if the adjoint time-step ratios fulfill a new step-ratio restriction $0<r_k:= au_k/ au_{k-1}\le r_{\mathrm{user}}$ ($r_{\mathrm{user}}$ can be chosen by the user such that $r_{\mathrm{user}}<4.864$), such that it is mesh-robustly stable in an energy norm. We view the BDF2 formula as a convolution approximation of the first time derivative and perform the error analysis by using the recent suggested discrete orthogonal convolution kernels. By developing some novel convolution embedding inequalities with respect to the orthogonal convolution kernels, an $L^2$ norm error estimate is established at the first time under the updated step-ratio restriction $0<r_k\le r_{\mathrm{user}}$.The time-stepping scheme is mesh-robustly convergent in the sense that the convergence constant (prefactor) in the error estimate is independent of the adjoint time-step ratios. On the basis of ample tests on random time meshes, a useful adaptive time-stepping strategy is applied to efficiently capture the multi-scale behaviors and to accelerate the long-time simulation approaching the steady state.

研究の動機と目的

• 可変ステップ幅BDF2スキームにおけるCahn-Hilliardモデルのエネルギー安定性と収束性を維持する課題に対処すること。
• 時間刻みの変動が長時間シミュレーションにおける数値的安定性と精度に与える影響を分析すること。
• 時間刻み比に依存しない誤差境界を持つ、メッシュに頑健な収束理論を確立すること。
• マルチスケールダイナミクスを効果的に捉えることができる実用的な適応的時間刻み戦略の開発。

提案手法

• BDF2法は、可変時間刻みを用いた時間微分の畳み込み近似として定式化される。
• 時間離散化フレームワーク内での正確な誤差解析を実現するために、離散直交畳み込みカーネルが用いられる。
• 時間誤差の増大を制御するため、直交畳み込みカーネルに対して新規の畳み込み埋め込み不等式が導出される。
• エネルギー散逸を離散レベルで維持するため、$ 0 < r_k := \Delta t_k / \Delta t_{k-1} \leq r_{\text{user}} $ というステップ比制限 $ r_{\text{user}} < 4.864 $ が導入される。
• この方法は、修正されたエネルギー散逸法則を満たすことが証明され、エネルギーノルムにおける安定性が保証される。
• ランダムな時間メッシュ上での数値的実験に基づき、安定状態への収束を加速するための適応的時間刻み戦略が実装される。

実験結果

リサーチクエスチョン

• RQ1 可変ステップ幅BDF2法は、任意の時間刻みの変動下でもCahn-Hilliardモデルに対してエネルギー散逸を保持できるか？
• RQ2 どのステップ比制限がエネルギーノルムおよび $ L^2 $ ノルムにおけるメッシュに頑健な安定性と収束性を保証するか？
• RQ3 離散直交畳み込みカーネルをどのように用いることで、非一様時間メッシュにおける鋭い誤差推定を導出できるか？
• RQ4 収束定数は時間刻み比にどの程度依存するか？また、それをそれらに依存しない形にできるか？
• RQ5 適応的時間刻み戦略は、マルチスケールダイナミクスを解像し、長時間シミュレーションを加速するためにどの程度有効か？

主な発見

• ステップ比制限 $ r_k \leq r_{\text{user}} < 4.864 $ の下で、可変時間刻みを用いたBDF2法は離散レベルで修正されたエネルギー散逸法則を保持し、エネルギー安定性が保証される。
• 更新されたステップ比条件の下で、最初の時間ステップにおいて新たな $ L^2 $ ノルム誤差推定が確立され、収束定数が時間刻み比に依存しない。
• 収束はメッシュに頑健であり、誤差境界の係数が隣接する時間刻み比に依存しない。
• 離散直交畳み込みカーネルと導出された埋め込み不等式の活用により、$ L^2 $ ノルムにおける最適収束率が達成される。
• ランダムな時間メッシュ上での数値実験により、適応的時間刻み戦略の頑健性と効率性が確認され、マルチスケール挙動の捉え込みに有効である。
• この手法により、安定性と精度を維持したまま、動的時間刻みの調整によって長時間シミュレーションが加速可能となる。