[論文レビュー] Error analysis for a fractional-derivative parabolic problem on quasi-graded meshes using barrier functions

Natalia Kopteva, Xiangyun Meng|arXiv (Cornell University)|May 17, 2019

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

ひとこと要約

本稿は、準グレーデッド時刻メッシュ上でのCaputo時間微分を用いた分数階微分放物型問題に対する、安定性および誤差解析の新規な手法を提示する。バリア関数を用いて、L1型およびアリクハノフ型スキームの両方について、時刻に依存する点ごとの誤差境界を鋭く確立した。また、最適なメッシュグレーディングよりも緩いものでも、正の時刻において最適収束率が達成されることを示しており、実用的意義と理論的新規性を併せ持つ。

ABSTRACT

An initial-boundary value problem with a Caputo time derivative of fractional order $α\in(0,1)$ is considered, solutions of which typically exhibit a singular behaviour at an initial time. For this problem, we give a simple and general numerical-stability analysis using barrier functions, which yields sharp pointwise-in-time error bounds on quasi-graded temporal meshes with arbitrary degree of grading. L1-type and Alikhanov-type discretization in time are considered. In particular, those results imply that milder (compared to the optimal) grading yields optimal convergence rates in positive time. Semi-discretizations in time and full discretizations are addressed. The theoretical findings are illustrated by numerical experiments.

研究の動機と目的

• 時間離散化された分数階微分放物型問題（Caputo微分を用いる）に対して、一般的かつ単純な数値安定性解析を構築すること。
• 任意のグレーディング度数を持つ準グレーデッド時刻メッシュ上での、点ごとの時刻に依存する誤差境界を確立すること。
• 最適でない（より緩い）グレーディングであっても、正の時刻において最適収束率が得られることを示し、実用的意義を有すること。
• バリア関数手法をL1型およびアリクハノフ型時間離散化に拡張し、広範な適用可能性を確保すること。
• 理論的予測を数値実験により検証し、鋭さと収束挙動の妥当性を確認すること。

提案手法

• 離散分数階微分作用素が離散的最大原理を満たすことを保証するバリア関数を用いて、離散作用素の安定性を解析する。
• 解の大きさをソース項の減衰率で評価する重要な安定性推定式（1.2）を導出する。この推定式は、メッシュグレーディングの明示的依存性を有する。
• L1型スキームにバリア関数アプローチを適用する。L1型スキームは、区分的線形補間と畳み込み・四元数法により定義される。
• 分数階サブディフュージョン問題において高次精度を達成するアリクハノフ型スキームに対しても解析を拡張する。
• 離散作用素行列の逆単調性を用いて、バリア関数構築の妥当性を保証する。
• 初期値問題の数値実験を通じて、理論的境界の妥当性を検証する。実験では、分数階数αとメッシュグレーディングパラメータを変化させた。

実験結果

リサーチクエスチョン

• RQ1 準グレーデッド時刻メッシュ上での分数階微分放物型問題に対して、一般的かつ単純な安定性解析を構築可能か？
• RQ2 任意のグレーディングを持つメッシュ上でのL1型およびアリクハノフ型スキームの、点ごとの時刻に依存する誤差挙動はどのように特定できるか？
• RQ3 最適なグレーディングよりも緩い（非最適な）グレーディングでも、正の時刻において最適収束率が達成されるか？
• RQ4 理論的誤差境界は、数値結果と比較して鋭いか、収束率は一致するか？
• RQ5 バリア関数手法は、L1スキームを超えるさまざまな離散分数階微分作用素に対しても体系的に適用可能か？

主な発見

• 提案されたバリア関数手法により、任意のグレーディング度数を持つ準グレーデッドメッシュ上でのL1型およびアリクハノフ型スキームに対して、点ごとの時刻に依存する誤差境界が鋭く得られた。
• 解析により、最適グレーディングよりも緩いグレーディング（L1ではr < (2−α)/α、アリクハノフではr < (3−α)/α）であっても、正の時刻で最適収束率が達成されることを示した。
• L1法では、数値実験により、r=1のとき収束率約1.0、r=(2−α)/0.9のとき約1.4、r=(2−α)/αのとき約1.6の収束率が観測され、理論的予測と整合的であった。
• アリクハノフ法では、r=2のとき収束率約2.0、r=(3−α)/0.95のとき約2.5の収束率が観測され、理論的境界M^{-min{r,3−α}}と一致した。
• アリクハノフ法のノード最大誤差は、M^{-min{αr,3−α}}の速度で減少し、定理18および注意20から得られるグローバル誤差境界を確認した。
• 数値結果により、誤差境界が鋭く、異なるαおよびメッシュパラメータにおいて、点ごとの誤差挙動が理論的推定式（4.5）と密接に一致することが示された。