[論文レビュー] Unconditional Stability Of A Two-Step Fourth-Order Modified Explicit Euler/Crank-Nicolson Approach For Solving Time-Variable Fractional Mobile-Immobile Advection-Dispersion Equation

Eric Ngondiep|arXiv (Cornell University)|May 6, 2022

Fractional Differential Equations Solutions被引用数 5

ひとこと要約

本稿では、時間変動型の分数マイトリー・イミュービル移動拡散方程式（変数係数Caputo微分を伴う）を解くために、2段階の4次精度を有する修正 explicit Euler/Crank-Nicolson スキームを提案する。この手法は、時間と空間の両方で無条件安定性および $O(k + h^4)$ の収束性を達成しており、既存の1段階スキームと比較して優れた効率性と精度を示している。

ABSTRACT

This paper considers a two-step fourth-order modified explicit Euler/Crank-Nicolson numerical method for solving the time-variable fractional mobile-immobile advection-dispersion model subjects to suitable initial and boundary conditions. Both stability and error estimates of the new approach are deeply analyzed in the $L^{\infty}(0,T;L^{2})$-norm. The theoretical studies show that the proposed technique is unconditionally stable with convergence of order $O(k+h^{4})$, where $h$ and $k$ are space step and time step, respectively. This result indicate that the two-step fourth-order formulation is more efficient than a broad range of numerical schemes widely studied in the literature for the considered problem. Numerical experiments are performed to verify the unconditional stability and convergence rate of the developed algorithm.

研究の動機と目的

• 時間変動型の分数マイトリー・イミュービル移動拡散方程式（変数係数Caputo微分を伴う）に対して、高次精度で無条件安定な数値解法を開発すること。
• 既存の1段階スキームがしばしば低い時間精度や条件付き安定性を示すという限界を克服すること。
• 計算効率を維持したまま、4次精度の空間収束と1次精度の時間収束を達成すること。
• 提案スキームの $L^{ u}(0,T;L^{2})$-ノルムにおける厳密な安定性および誤差解析を提供すること。
• 正確な解を用いた数値実験と収束率解析を通じて、理論的結果の妥当性を検証すること。

提案手法

• 修正 explicit Euler と Crank-Nicolson スキームを組み合わせた2段階4次精度の定式化により、精度と安定性のバランスを図る。
• 時間微分は2次精度の2段階スキームで近似され、空間微分は4次精度のコンact差分法で離散化される。
• 変数係数Caputo分数微分は、特異核を含む定義に基づく数値積分則（四則則）を用いて近似される。
• エネルギー法を用いた $L^{ u}(0,T;L^{2})$-ノルムにおける安定性解析により、時間刻みと空間刻みのサイズに依存しない無条件安定性が証明される。
• 線形方程式系として定式化され、効率的な反復解法で解ける。今後の研究ではプリコンディショニングの可能性も検討される。
• 分数微分項を半陽解的（semi-implicit）に取り扱うことで、安定性が確保され、より大きな時間刻みが可能になる。

実験結果

リサーチクエスチョン

• RQ1 2段階4次精度の数値スキームは、時間変動型の分数マイトリー・イミュービル移動拡散方程式に対して無条件安定性を達成できるか？
• RQ2 提案スキームの時間および空間における収束次数は何か？ また、既存の1段階スキームを上回る性能を示すか？
• RQ3 分数階数が時間的・空間的に変動する場合、本手法の安定性および精度はどのように変化するか？
• RQ4 $L^{ u}(0,T;L^{2})$-ノルムにおける理論的誤差境界が、正確な解を用いた数値実験によって検証可能か？
• RQ5 本スキームは、同種の問題に対して広く用いられている既存スキームと比較して、計算的により効率的か？

主な発見

• 提案スキームは $L^{ u}(0,T;L^{2})$-ノルムにおいて無条件安定であり、時間刻みと空間刻みのサイズにかかわらず安定性が保証される。
• 時間方向に1次収束（$O(k)$）、空間方向に4次収束（$O(h^4)$）を達成し、全体の収束次数は $O(k + h^4)$ となる。
• 数値実験により理論的収束率が確認され、複数のテストケースにおいて空間方向の収束次数が約4.0、時間方向が3.8～4.1の範囲で観測された。
• 例1では、$h = 1/16$ のとき最大誤差は $1.1488 \times 10^{-6}$ であり、収束率は4.1029であり、4次空間精度が確認された。
• 例2（$α = 0.25$）では、$h$ の減少に伴い誤差が $1.0918 \times 10^{-2}$ から $2.5995 \times 10^{-6}$ に減少し、収束率は4.1328であった。これにより高次精度が裏付けられた。
• 理論的解析および数値結果の両方により、広範な既存スキームと比較して、本手法はより効率的かつ高速であることが示された。