[论文解读] 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 分数阶导数的边值问题建立了变分格式，建立了强制性、稳定性及转移定理，从而为有限元方法实现最优收敛速率提供了保障。主要贡献在于对有限元误差进行了严格分析，表明收敛速率仅依赖于右端项的正则性，并通过数值实验得到验证。

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$ 光滑度的伽辽金有限元逼近的收敛速率。
• 通过计算实验对理论误差估计进行数值验证。

提出的方法

• 使用分数阶 Sobolev 空间（Riemann-Liouville 情况下为 $H_0^{\alpha/2}(0,1)$，Caputo 情况下为不同的试函数/测试函数空间）建立弱问题。
• 在 $H_0^{\alpha/2}(0,1)$ 上证明 Riemann-Liouville 变分格式的强制性，确保问题适定。
• 建立转移定理，表明解的正则性与源项 $f$ 的光滑度成正比，从而实现最优有限元收敛速率。
• 在均匀网格上应用 Galerkin 有限元方法，采用分段多项式逼近。
• 推导 $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}$-范数收敛速率，观测到的收敛率约为 1.0–1.5（$\alpha \in (1,2)$ 时）。
• 在 $L^2$-范数下，实际观测的收敛速率比理论估计高出约半个阶次，表明仍存在需进一步研究的差距。
• 对于如 $f(x) = x^{-1/4}$ 的奇异源项，当 $\alpha \to 1^+$ 时收敛速度减慢，尤其在 Caputo 情况下更为明显，但收敛速率仍与理论预期一致。