[论文解读] A simple finite element method for the boundary value problem with a Riemann-Liouville derivative

Bangti Jin, Raytcho Lazarov+2 more|arXiv (Cornell University)|Feb 27, 2015

Fractional Differential Equations Solutions参考文献 14被引用 3

一句话总结

本文提出了一种针对具有 Riemann-Liouville 分数阶导数（阶数 $α \in (3/2, 2)$）的边值问题的新型有限元方法，通过 Riemann-Liouville 积分变换来平滑解的奇异性，从而在标准有限元方法下实现最优的 $L^2$ 收敛率。该方法实现了最优收敛性，并可自然推广至分数阶 Sturm-Liouville 问题。

ABSTRACT

We consider a boundary value problem involving a Riemann-Liouville fractional derivative of order $α\in (3/2,2)$ on the unit interval $(0,1)$. The standard Galerkin finite element approximation converges slowly due to the presence of singularity term $x^{α-1}$ in the solution representation. In this work, we develop a simple technique, by transforming it into a second-order two-point boundary value problem with nonlocal low order terms, whose solution can reconstruct directly the solution to the original problem. The stability of the variational formulation, and the optimal regularity pickup of the solution are analyzed. A novel Galerkin finite element method with piecewise linear or quadratic finite elements is developed, and $L^2(D)$ error estimates are provided. The approach is then applied to the corresponding fractional Sturm-Liouville problem, and error estimates of the eigenvalue approximations are given. Extensive numerical results fully confirm our theoretical study.

研究动机与目标

• 解决标准 Galerkin 有限元方法在具有 Riemann-Liouville 导数的分数阶边值问题中因解的奇异性导致的收敛速度缓慢的问题。
• 开发一种稳定、高效且精确的有限元方法，即使在存在 $x^{\alpha-1}$ 奇异性的情况下也能保持最优收敛率。
• 将该方法推广至具有 Riemann-Liouville 导数的分数阶 Sturm-Liouville 问题，实现精确的特征值近似。
• 为变换后的变分格式及其有限元离散化提供严格的稳定性与误差分析。

提出的方法

• 通过 Riemann-Liouville 积分算子 ${}_0I_x^{2-\alpha}$ 将原始分数阶边值问题转化为具有非局部低阶项的等价二阶两点边值问题。
• 该变换将奇异性项 $x^{\alpha-1}$ 映射为光滑函数 $x$，从而实现标准有限元逼近。
• 利用加权 Sobolev 空间推导变换后问题的稳定变分格式，并分析其适定性。
• 对变换后的问题应用一致连续分段线性与二次有限元，获得最优的 $L^2(D)$ 误差估计。
• 通过分数阶微分从变换后的解重构原始解，利用 expm1 函数实现稳定数值计算。
• 将该方法扩展至分数阶 Sturm-Liouville 问题，并推导特征值近似值的收敛率。

实验结果

研究问题

• RQ1 基于变换的有限元方法是否能在存在解奇异性的情况下，对具有 Riemann-Liouville 导数的分数阶边值问题实现最优的 $L^2$ 收敛率？
• RQ2 通过 Riemann-Liouville 积分 ${}_0I_x^{2-\alpha}$ 变换是否能有效将奇异性 $x^{\alpha-1}$ 平滑为一个正则函数？
• RQ3 所提出的该方法能否推广至具有可证明收敛性的分数阶 Sturm-Liouville 问题的特征值近似？
• RQ4 变换变量下变分格式的稳定性和正则性行为如何？

主要发现

• 变换后的问题具有稳定的变分格式，具备最优正则性提升，从而为有限元逼近提供了最优收敛率。
• 建立了分段线性元的 $O(h^2)$ 和分段二次元的 $O(h^3)$ 最优 $L^2(D)$ 误差估计。
• 该方法通过变换而非显式减去奇异性项，显著优于标准 Galerkin 有限元方法的收敛性能。
• 该方法可推广至分数阶 Sturm-Liouville 问题，且已建立特征值近似的收敛率。
• 数值结果验证了理论误差估计，并展示了该方法在各种测试用例中的高效性与鲁棒性。
• 通过使用 expm1 函数实现了分数阶微分在解重构过程中的稳定计算，避免了舍入误差。