[论文解读] A priori estimates of a finite element method for fractional diffusion problems by energy arguments

Samir Karaa, Kassem Mustapha+1 more|arXiv (Cornell University)|May 30, 2016

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

一句话总结

该论文为在凸域上具有可变扩散系数的时间分数阶扩散方程，建立了伽辽金分段线性有限元方法的先验误差估计。通过采用带有 $t^m$-加权范数的精细能量分析方法，并结合分数阶导数的广义莱布尼茨公式，该研究建立了最优的 $L^2$ 和 $H^1$ 误差界，以及对光滑和非光滑初始数据均适用的拟最优 $L^\infty$ 误差界。

ABSTRACT

In this article, the Galerkin piecewise-linear finite element (FE) method is applied to approximate the solution of time-fractional diffusion equations with variable diffusivity on bounded convex domains. Standard energy arguments do not provide satisfactory results for such a problem due to the low regularity of its exact solution. Using a delicate energy analysis, {\it a priori} optimal error bounds in $L^2(\Omega)$-, $H^1(\Omega)$-, and quasi-optimal in $L^{\infty}(\Omega)$-norms are derived for the semidiscrete FE scheme for cases with smooth and nonsmooth initial data. The main tool of our analysis is based on a repeated use of an integral operator and use of a $t^m$ type of weights to take care of the singular behavior at $t=0.$ The generalized Leibniz formula for fractional derivatives is found to play a key role in our analysis. Numerical experiments are presented to illustrate some of the theoretical results.

研究动机与目标

• 推导具有可变扩散系数的时间分数阶扩散问题中半离散有限元方法的先验误差界。
• 解决因 $t=0$ 处解的正则性较低而削弱标准能量技术的问题。
• 为光滑和非光滑初始数据建立最优的 $L^2$ 和 $H^1$ 收敛速率，以及拟最优的 $L^\infty$ 误差界。
• 提出一种新颖的能量分析框架，结合 $t^m$-加权范数和分数阶莱布尼茨公式，以捕捉初始时刻的奇异性行为。

提出的方法

• 将伽辽金分段线性有限元方法应用于有界凸域上的半离散时间分数阶扩散问题。
• 通过重复应用积分算子来处理分数阶导数的结构。
• 在能量分析中引入 $t^m$ 类型的权重，以控制 $t=0$ 处解的奇异性行为。
• 将分数阶导数的广义莱布尼茨公式作为关键分析工具。
• 进行精细的能量估计，以考虑精确解的低正则性。
• 通过系统的分部积分和范数控制，推导出 $L^2(\Omega)$、$H^1(\Omega)$ 和 $L^\infty(\Omega)$ 范数下的误差界。

实验结果

研究问题

• RQ1 能否对标准能量论证进行调整，以在具有低正则性解的分数阶扩散问题中实现最优误差估计？
• RQ2 如何利用 $t^m$-加权范数改善具有初始层奇异性的时间分数阶PDE的分析？
• RQ3 分数阶导数的广义莱布尼茨公式在推导精确误差界中起到什么作用？
• RQ4 在半离散有限元设置下，能否为光滑和非光滑初始数据均实现最优的 $L^2$ 和 $H^1$ 误差估计？
• RQ5 时间分数阶扩散方程的有限元逼近能否实现拟最优的 $L^\infty$ 收敛？

主要发现

• 为半离散有限元格式推导出对光滑和非光滑初始数据均有效的最优 $L^2(\Omega)$-范数误差界。
• 通过结合 $t^m$-权重的精细能量分析，建立了最优的 $H^1(\Omega)$-范数误差界。
• 获得了拟最优的 $L^\infty(\Omega)$-范数误差界，表明在最大范数下实现了近乎最优的收敛性。
• 分数阶导数的广义莱布尼茨公式在能量估计中处理分数阶导数乘积法则时至关重要。
• 重复使用积分算子与 $t^m$-加权范数的结合，有效捕捉了 $t=0$ 处的初始层奇异性。
• 数值实验验证了理论收敛速率在 $L^2$、$H^1$ 和 $L^\infty$ 范数下的表现。