[论文解读] 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

一句话总结

本文为时间分数阶 Fokker-Planck 方程（分数阶扩散指数 α ∈ (0,1]）的空间离散 Galerkin 逼近建立了统一稳定性。作者推导出稳定性估计，其常数在 α 上一致有界，考虑了非齐次源项，并提供了梯度稳定性。这些结果改进了先前的工作，消除了稳定性常数中与 α 相关的发散问题，简化了有限元方法的误差分析。

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 方程的半离散 Galerkin 解在所有 α ∈ (0,1] 上建立统一稳定性界，包括经典极限 α = 1。
• 在稳定性分析中考虑非齐次源项 g，扩展了以往仅限于 g ≡ 0 的结果。
• 推导 Galerkin 解梯度的稳定性估计，增强空间正则性分析。
• 通过新稳定性框架简化标准有限元逼近的误差界。
• 将分析扩展至间断 Galerkin 时间格式，在适当条件下证明 α-一致稳定性。

提出的方法

• 通过时间积分变分形式（4）表述弱解，以处理低正则性数据，确保适定性。
• 定义半离散 Galerkin 解 uX ∈ X ⊂ H¹₀(Ω)，满足对所有测试函数 χ ∈ X 的变分方程（6）。
• 利用分数阶积分估计与能量技巧推导 α-一致有界界，避免在 α → 1 时发生发散。
• 使用 Poincaré 不等式与涉及 Ψ(α) 的分数阶积分不等式（38），控制 Riemann–Liouville 分数阶导数项。
• 在时间方向采用间断 Galerkin 方法，结合分段多项式逼近，推导出在 α-依赖但有界的常数下的全离散稳定性。
• 利用稳定性框架，推导出分段线性有限元解的最优误差界，显式依赖于 α 和正则性参数 r。

实验结果

研究问题

• RQ1 时间分数阶 Fokker-Planck 方程的半离散 Galerkin 解的稳定性常数能否与分数阶指数 α ∈ (0,1] 无关？
• RQ2 非零源项 g 的存在如何影响 Galerkin 解的稳定性？能否仍实现统一有界性？
• RQ3 Galerkin 解的梯度在 α-一致稳定性方面表现如何？
• RQ4 新稳定性框架能否简化或改进标准有限元逼近的现有误差估计？
• RQ5 间断 Galerkin 时间格式是否能保持分数阶扩散方程的 α-一致性稳定性？

主要发现

• 本文建立了半离散 Galerkin 解的统一稳定性：对所有 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} 随 α→1 而发散的问题。
• Galerkin 解的梯度满足类似的 α-一致估计，支持更优的正则性分析。
• 对于分段线性有限元解 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)。
• 对于间断 Galerkin 时间格式，解满足 ∥Un⁻∥² + ∑ⱼ∥JUKⱯ∥² + ∫₀ᵗ ⟨∂¹⁻ᵃᵗ κ∇U, ∇U⟩ dt ≤ ∥u₀X∥² + CΩ Ψ(α)/κmin tₙ^{1−α} ∫₀ᵗ ‖g(s)‖² ds，其中 Ψ(α) 在 α→1 时有界，但在 α→0 时发散。