[論文レビュー] $H^1$-norm stability and convergence of an L2-type method on nonuniform meshes for subdiffusion equation

Chaoyu Quan, Xu Wu|arXiv (Cornell University)|May 12, 2022

Differential Equations and Numerical Methods被引用数 4

ひとこと要約

本稿では、一般非一様メッシュ上でサブディフュージョン方程式のL2型スキームについて、H¹ノルムの安定性と収束性を確立する。時間刻み比にやや緩い制約（0.457 ≤ ρₖ ≤ 3.562）を課した場合、重要な双一次形式の半正定値性を証明し、これにより、r > 5/α − 1を満たす修正型グレーデッドメッシュ上でのH¹ノルムにおける収束順序(5−α)/2の長期間安定性と収束性を確立する。

ABSTRACT

This work establishes $H^1$-norm stability and convergence for an L2 method on general nonuniform meshes when applied to the subdiffusion equation. Under mild constraints on the time step ratio $ρ_k$, such as $0.4573328\leq ρ_k\leq 3.5615528$ for $k\geq 2$, the positive semidefiniteness of a crucial bilinear form associated with the L2 fractional-derivative operator is proved. This result enables us to derive long time $H^1$-stability of L2 schemes. These positive semidefiniteness and $H^1$-stability properties hold for standard graded meshes with grading parameter $15/α-1$. To the best of our knowledge, this study is the first work on $H^1$-norm stability and convergence of L2 methods on general nonuniform meshes for the subdiffusion equation.

研究の動機と目的

• 一般非一様時刻メッシュ上でのサブディフュージョン方程式に対するL2型スキームのH¹ノルムの安定性と収束性を確立すること。
• 特に低正則性解の下でL2法のH¹安定性解析が不足していることに対処すること。
• 時間刻み比にやや緩い制約の下で、L2分数階微分作用素に関連する双一次形式の半正定値性を証明すること。
• 特定のグレーディングパラメータを有するグレーデッドメッシュおよび修正型グレーデッドメッシュへ安定性と収束結果を拡張すること。
• 厳密なH¹ノルム誤差解析を提供し、数値的に収束速度を確認すること。

提案手法

• k ≥ 2 の場合に時間刻み比制約 ρₖ ∈ [0.4573328, 3.5615528] の下で、双一次形式 $\mathcal{B}_n(v,w) = \sum_{k=1}^n \langle L_k^\alpha v, \delta_k w \rangle$ の半正定値性を証明する。
• 長期間のH¹安定性を導出するため、離散分数階グロワールド不等式とエネルギー技法に基づく新しい解析フレームワークを用いる。
• 標準的グレーデッドメッシュ（グレーディングパラメータ r ≤ 3.2016538）に対しても、H¹安定性を確立する。
• 一般非一様メッシュ上でのH¹ノルム誤差推定を導出し、r > 5/α − 1 を満たす修正型グレーデッドメッシュ上では収束順序 (5 − α)/2 を証明する。
• 空間方向には、チェビシェフ＝ガウス＝ロバッチョフ点を用いたスペクトルコロケーション法を適用し、高い空間精度を達成する。
• 最大H¹誤差を測定することで、変化する α、r、N の条件下で、修正型グレーデッドメッシュ上での数値実験により理論的結果を検証する。

実験結果

リサーチクエスチョン

• RQ1 時間刻み比にやや緩い制約の下で、一般非一様メッシュ上でのL2分数階微分作用素に関連する双一次形式は半正定値であるか？
• RQ2 サブディフュージョン方程式に対して、非一様メッシュ上での陰的L2スキームの長期間H¹安定性は確立可能か？
• RQ3 修正型グレーデッドメッシュ上でのL2法におけるH¹ノルムの最適収束速度は何か？また、これはグレーディングパラメータ r にどのように依存するか？
• RQ4 数値的検証において、双一次形式の半正定値性結果は鋭いか？
• RQ5 理論的H¹ノルム収束速度 (5 − α)/2 は実際の数値実験でも観察可能であり、観測された収束順序と整合するか？

主な発見

• 時間刻み比制約 0.4573328 ≤ ρₖ ≤ 3.5615528（k ≥ 2）の下で、双一次形式 $\mathcal{B}_n$ が半正定値であることが証明され、これによりH¹安定性解析が可能となる。
• 一般非一様メッシュ上でのサブディフュージョン方程式に対して、同次ディリクレ境界条件のもとでL2スキームの長期間H¹安定性が確立される。
• 標準的グレーデッドメッシュでは、グレーディングパラメータが 1 < r ≤ 3.2016538 を満たす場合にH¹安定性が成立する。
• r > 5/α − 1 を満たす修正型グレーデッドメッシュ上では、H¹ノルムにおける収束順序 (5 − α)/2 が証明される。
• 数値結果により、H¹ノルムにおける観測収束順序が、理論的期待と整合的に $\min\{r\alpha, 3 - \alpha\}$ に近いことが確認された。
• 最小固有値の下界 $\lambda_{\min} \geq C N^{\alpha}$ が数値的に最適であることが確認され、最大固有値の下界 $\lambda_{\max} \geq C N^{r\alpha}$ も同様に成立する。