[論文レビュー] Variational formulation of problems involving fractional order differential operators

Bangti Jin, Raytcho Lazarov+1 more|arXiv (Cornell University)|Jul 17, 2013

Fractional Differential Equations Solutions参考文献 10被引用数 7

ひとこと要約

本稿では、$\alpha \in (1,2)$ の順序をもつ Caputo および Riemann-Liouville 分数マージン微分を含む境界値問題の変分形式を構築し、コercivity、安定性、シフト定理を確立することで、有限要素法における最適収束率の達成を可能にする。主な貢献は、右辺の正則性にのみ依存する $H^{\alpha/2}$-および $L^2$-ノルムにおける厳密な有限要素誤差解析であり、数値実験によって裏付けられている。

ABSTRACT

In this work, we consider boundary value problems involving Caputo and Riemann-Liouville fractional derivatives of order $α\in(1,2)$ on the unit interval $(0,1)$. These fractional derivatives lead to non-symmetric boundary value problems, which are investigated from a variational point of view. The variational problem for the Riemann-Liouville case is coercive on the space $H_0^{α/2}(0,1)$ but the solutions are less regular, whereas that for the Caputo case involves different test and trial spaces. The numerical analysis of these problems requires the so-called shift theorems which show that the solutions of the variational problem are more regular. The regularity pickup enables one to establish convergence rates of the finite element approximations. Finally, numerical results are presented to illustrate the error estimates.

研究の動機と目的

• 順序 $\alpha \in (1,2)$ の Caputo および Riemann-Liouville 分数マージン微分を含む分数マージン境界値問題の変分形式を構築すること。
• 解の正則性が向上することを保証する変分的安定性およびシフト定理を確立すること。
• 有限要素近似の収束率を、右辺 $f$ の滑らかさにのみ依存する形で導出すること。
• 理論的誤差推定値の数値的妥当性を、計算実験を通じて提供すること。

提案手法

• 分数マージンソボレフ空間、特に Riemann-Liouville では $H_0^{\alpha/2}(0,1)$ を用い、Caputo では別個の試行／試験関数空間を用いる弱問題の定式化。
• Riemann-Liouville 変分形式が $H_0^{\alpha/2}(0,1)$ 上でコーエルシビティであることを証明し、適切な定式化を保証。
• 解がソース項 $f$ の滑らかさに比例して正則性を獲得することを示すシフト定理を確立し、最適な有限要素収束を可能にする。
• 一様メッシュ上での区分的多項式近似を用いたガラーキン有限要素法の適用。
• $L^2$ および $H^{\alpha/2}$-ノルムにおける誤差推定値を導出し、$f$ の正則性に依存する収束率を明示。
• 特異ソース項を用いた数値実験により、$\alpha$ や $f$ の正則性が変化する条件下での理論的収束率の妥当性を検証。

実験結果

リサーチクエスチョン

• RQ1 順序 $\alpha \in (1,2)$ の Caputo および Riemann-Liouville 分数マージン微分を含む変分形式を、一貫してどのように導出できるか。
• RQ2 これらの変分問題の解に成り立つ正則性特性（シフト定理）は何か、そしてソース項 $f$ にどのように依存するか。
• RQ3 これらの分数マージン問題に対するガラーキン有限要素近似の収束率を、どのように厳密に確立できるか。
• RQ4 $L^2$-ノルム誤差推定値が実際の収束率を下回る理由は何か、その乖離はどのように説明できるか。
• RQ5 ソース $f$ が特異的で $\alpha \to 1^+$ のとき、解が原点付近でどのように振る舞い、どのような是正策が提案されるか。

主な発見

• Riemann-Liouville 変分形式は $H_0^{\alpha/2}(0,1)$ 上でコーエルシビティを示し、適切な定式化を保証するが、解は $x^{\alpha-1}$ の形の特異性を示し、一様メッシュ上では振動を引き起こす。
• Caputo の形式は別個の試行・試験関数空間を必要とし、適切な関数空間設定のもとで変分形式は安定である。
• シフト定理により、解の正則性がソース項 $f$ の滑らかさに比例して向上することが示され、最適な有限要素収束率が達成可能となる。
• 数値結果により、理論的予測と一致する $H^{\alpha/2}$-ノルムにおける収束率が確認され、$\alpha \in (1,2)$ の範囲で観察される収束率は約 1.0～1.5 であった。
• $L^2$-ノルムでは、理論的推定値よりも約半次の高い収束率が観測されたが、これは未解決のギャップを示しており、さらなる調査を要する。
• 例えば $f(x) = x^{-1/4}$ のような特異ソースでは、$\alpha \to 1^+$ に近づくと収束が遅くなるが、特に Caputo の場合に顕著であり、それでも理論的期待値と一致する収束率を示した。