[論文レビュー] Sharp $H^1$-norm error estimates of two time-stepping schemes for reaction-subdiffusion problems

Jincheng Ren, Hong-lin Liao+2 more|arXiv (Cornell University)|Nov 20, 2018

Fractional Differential Equations Solutions参考文献 9被引用数 21

ひとこと要約

本稿は、キャプートー微分を伴う反応-サブディファージョン問題に対する非一様時間刻みの2つのスキーム—L1および分数階クランク・ニコルソン—について、鋭いH¹ノルム誤差推定を確立する。初期特異性とキャプートー微分の畳み込み構造に起因する古典的解析で失われる最適収束率を回復するため、改良された離散グローバル不等式を導入し、時間・空間誤差分解技術とグローバル一致性解析を組み合わせることで、最適なグレーディングパラメータを有する非一様メッシュ上でも、両スケームに対してH¹ノルム誤差境界𝒪(τ^{2−α})を達成する。

ABSTRACT

Due to the intrinsically initial singularity of solution and the discrete convolution form in numerical Caputo derivatives, the traditional $H^1$-norm analysis (corresponding to the case for a classical diffusion equation) to the time approximations of a fractional subdiffusion problem always leads to suboptimal error estimates (a loss of time accuracy). To recover the theoretical accuracy in time, we propose an improved discrete Grönwall inequality and apply it to the well-known L1 formula and a fractional Crank-Nicolson scheme. With the help of a time-space error-splitting technique and the global consistency analysis, sharp $H^1$-norm error estimates of the two nonuniform approaches are established for a reaction-subdiffusion problems. Numerical experiments are included to confirm the sharpness of our analysis.

研究の動機と目的

• 初期特異性とキャプートー微分の畳み込み構造に起因する、分数階サブディファージョン問題における時間方向H¹ノルム誤差推定の非最適性を解消すること。
• 通常の古典的解析で失われる非一様時間離散化におけるH¹ノルムの最適時間収束率を回復すること。
• リャオらの理論枠組みを拡張し、一般の非一様メッシュ上でのL1および分数階クランク・ニコルソンスキームに対して、鋭いH¹ノルム誤差推定を達成すること。
• 変化するメッシュグレーディングパラメータおよび分数階数を用いた数値実験を通じて、理論的推定の鋭さを検証すること。

提案手法

• 初期特異性とキャプートー微分の畳み込み構造を扱えるように適合させた改良された離散グローバル不等式を提案する。
• 時間・空間誤差分解技術を適用し、時間誤差と空間誤差を分離することで、各成分を独立に解析可能にする。
• グローバル一貫性解析を用いて、特にt=0近傍で局所的打切り誤差を全体の時間区間で制御する。
• 初期特異性を解消するため、パラメータγを用いたグレーディング時間メッシュを用い、τk ∝ tk−1+γ−1とする。
• 改良された離散グローバル不等式と安定性・一貫性の議論を組み合わせることで、L1および分数階クランク・ニコルソンスキームの両方に対して鋭いH¹ノルム誤差推定を確立する。
• 有限差分法による空間離散化と変化するメッシュグレーディングパラメータを用いた数値実験を通じて、理論的結果を検証する。

実験結果

リサーチクエスチョン

• RQ1 初期特異性と畳み込み構造に起因するキャプートー微分を伴う反応-サブディファージョン問題の非一様時間離散化において、最適なH¹ノルム誤差推定を回復できるか。
• RQ2 t=0近傍での解の滑らかでない性質が、従来のH¹ノルム解析によって最適収束率を達成できない原因となり、その是正策は何か。
• RQ3 改良された離散グローバル不等式は、L1および分数階クランク・ニコルソンスキームの両方において、H¹ノルムでの最適時間収束率を回復するために効果的に利用可能か。
• RQ4 これらのスキームのH¹ノルムで理論的収束率𝒪(τ^{2−α})を達成するための最適グレーディングパラメータγは何か。
• RQ5 変化するメッシュパラメータおよび分数階数を用いた数値実験により、理論的H¹ノルム誤差境界が鋭いかどうかを確認できるか。

主な発見

• 本稿は、初期特異性に起因する非一様時間メッシュ上でも、L1および分数階クランク・ニコルソンスキームの両方に対して、H¹ノルム誤差推定𝒪(τ^{2−α})を達成する鋭い推定値を確立した。
• 改良された離散グローバル不等式は、H¹ノルム解析における精度の損失を効果的に緩和し、古典的手法が失敗する場面でも最適誤差境界を達成可能にした。
• 数値実験により、メッシュグレーディングパラメータγ ≥ γ_optのとき、収束率が𝒪(τ^{2−α})に達することが確認され、理論的推定の鋭さが裏付けられた。
• L1スキームではγσ ≥ 2−αのとき、分数階クランク・ニコルソンスキームではγσ ≥ 2のとき最適収束が達成され、γ_optはそれぞれに応じて決定される。
• 時間・空間誤差分解技術は、時間誤差と空間誤差を効果的に分離し、時間離散化によるH¹ノルム誤差寄与を精密に制御・解析可能にした。
• 数値結果から、非一様メッシュは均等メッシュに比べて収束率を著しく向上させ、特にγ ≥ γ_optのとき理論的予測と一致することが確認された。