[論文レビュー] Uniform stability for a spatially-discrete, subdiffusive Fokker-Planck equation

William McLean, Kassem Mustapha|arXiv (Cornell University)|Dec 26, 2020

Fractional Differential Equations Solutions参考文献 29被引用数 6

ひとこと要約

本稿は、分数階微分指数 α ∈ (0,1] を持つ時間分数階 Fokker-Planck 方程式の空間的離散化ガラーキン近似に対して、一様安定性を確立する。著者は、α に依存しない定数を有する安定性評価を導出し、非同次項を考慮し、勾配の安定性を提供する。これらの結果は、従来の研究において α に依存する定数の発散を解消することで向上し、有限要素法の誤差解析を簡素化する。

ABSTRACT

We prove stability estimates for the spatially discrete, Galerkin solution of a fractional Fokker-Planck equation, improving on previous results in several respects. Our main goal is to establish that the stability constants are bounded uniformly in the fractional diffusion exponent $\alpha\in(0,1]$. In addition, we account for the presence of an inhomogeneous term and show a stability estimate for the gradient of the Galerkin solution. As a by-product, the proofs of error bounds for a standard finite element approximation are simplified.

研究の動機と目的

• 時間分数階 Fokker-Planck 方程式の半離散ガラーキン解について、すべての α ∈ (0,1] において一様安定性境界を確立すること。
• 安定性解析において非同次項 g を取り入れることで、従来の g ≡ 0 に限った結果を拡張すること。
• ガラーキン解の勾配についての安定性評価を導出し、空間的正則性の解析を強化すること。
• 新しい安定性フレームワークを活用して、標準的な有限要素近似の誤差境界を簡素化すること。
• 不連続ガラーキン時間離散化への解析を拡張し、適切な条件下で α に一様な安定性を証明すること。

提案手法

• 低正則性のデータを扱うために、時間積分された変分形式 (4) を用いて弱解を定式化し、適切な定義を保証する。
• 半離散ガラーキン解 uX ∈ X ⊂ H¹₀(Ω) を定義し、すべての試験関数 χ ∈ X に対して変分方程式 (6) を満たす。
• 分数階積分推定とエネルギー技法を用いて、α → 1 のときの発散を避ける α に一様な境界を導出する。
• Poincaré 不等式と Riemann–Liouville 分数階微分項を制御するための分数階積分不等式 (38) を用い、Ψ(α) を含む。
• 時間方向に区分的多項式近似を用いた不連続ガラーキン法を適用し、α に依存するが有界な定数を用いて完全離散安定性を導出する。
• 安定性フレームワークを活用して、区分的線形有限要素解の最適誤差境界を導出し、α と正則性パラメータ r への明示的依存を示す。

実験結果

リサーチクエスチョン

• RQ1 時間分数階 Fokker-Planck 方程式の半離散ガラーキン解の安定性定数を、分数階指数 α ∈ (0,1] に依存しない形にできるか？
• RQ2 非ゼロの源項 g の存在がガラーキン解の安定性に与える影響は何か？ 一様境界を達成できるか？
• RQ3 ガラーキン解の勾配の振る舞いは、α に一様な安定性を満たすか？
• RQ4 新しい安定性フレームワークは、標準的な有限要素近似の既存の誤差推定を簡素化または改善できるか？
• RQ5 不連続ガラーキン時間離散化は、分数拡散方程式に対して α に一様な安定性を保持するか？

主な発見

• 本稿は、半離散ガラーキン解について一様安定性を確立した：すべての t ∈ (0,T] および α ∈ (0,1] に対して、‖uX(t)‖ + t^{α/2}‖∇uX(t)‖ ≤ C(‖u₀X‖ + ∫₀ᵗ ‖g(s)‖ds + (1/t ∫₀ᵗ ‖sg(g)‖² ds)^{1/2}) が成り立ち、C は α に依存しない。
• α → 1 のとき安定性定数が有界のまま保たれ、以前の研究で観察された Cα ∼ (1−α)^{-1} の発散問題が解消された。
• ガラーキン解の勾配についても同様の α に一様な評価が成立し、正則性解析の向上が可能になる。
• 区分的線形有限要素解 uh については、誤差境界が ∥uh − u∥ = O(t^{-α(2−r)/2} h²) および ∥∇uh − ∇u∥ = O(t^{-α(2−r)/2} h) と与えられ、r ∈ [0,2] は解の正則性を定量化する。
• r = 2 かつ u₀h = Rₕu₀ のとき、誤差境界は時間に一様になる：∥uh(t) − u(t)∥ + h∥∇uh(t) − ∇u(t)∥ ≤ Cα h² (∥u₀∥₂ + M)。
• 不連続ガラーキン時間離散化では、解が ∥Un⁻∥² + ∑ⱼ∥JUKⱼ∥² + ∫₀ᵗ ⟨∂¹⁻ᵃᵗ κ∇U, ∇U⟩ dt ≤ ∥u₀X∥² + CΩ Ψ(α)/κmin tₙ^{1−α} ∫₀ᵗ ‖g(s)‖² ds を満たし、Ψ(α) は α→1 のとき有界であるが、α→0 のとき発散する。