[论文解读] Finite difference methods for fractional Laplacians

Yanghong Huang, Adam M. Oberman|arXiv (Cornell University)|Nov 1, 2016

Fractional Differential Equations Solutions参考文献 33被引用 18

一句话总结

本论文提出了一种统一的有限差分框架，用于在一维均匀网格上近似分数阶拉普拉斯算子 $(-\Delta)^{\alpha/2}$，通过基于谱表示和傅里叶分析导出的权重进行离散卷积。主要贡献在于系统比较了多种离散化格式——包括Grünwald-Letnikov、Gorenflo-Mainardi以及基于求积的方法——展示了它们在非光滑解和远场边界条件下下的精度、稳定性与收敛行为。

ABSTRACT

The fractional Laplacian $(-Δ)^{α/2}$ is the prototypical non-local elliptic operator. While analytical theory has been advanced and understood for some time, there remain many open problems in the numerical analysis of the operator. In this article, we study several different finite difference discretisations of the fractional Laplacian on uniform grids in one dimension that takes the same form. Many properties can be compared and summarised in this relatively simple setting, to tackle more important questions like the nonlocality, singularity and flat tails common in practical implementations. The accuracy and the asymptotic behaviours of the methods are also studied, together with treatment of the far field boundary conditions, providing a unified perspective on the further development of the scheme in higher dimensions.

研究动机与目标

• 开发一维分数阶拉普拉斯算子 $(-\Delta)^{\alpha/2}$ 的统一数值框架，用于有限差分离散化。
• 基于分数阶拉普拉斯的不同定义，比较多种现有及新型离散化格式在精度、稳定性与渐近行为方面的表现。
• 解决实际实现中的挑战，如非局部性、奇异权重以及远场边界条件。
• 分析分数阶拉普拉斯算子PDE的收敛性与精度，特别是当解非光滑时的情况。
• 为将该方法推广至高维及具有通用符号的非局部算子奠定基础。

提出的方法

• 将离散分数阶拉普拉斯算子表述为卷积算子：$(-\Delta_h)^{\alpha/2}u_j = \sum_{k=-\infty}^{\infty} (u_j - u_{j-k}) w_k$，其中权重 $w_k$ 由谱表示导出。
• 利用半离散傅里叶分析将离散算子与其符号关联，实现谱分析与稳定性评估。
• 通过多种方法推导权重：傅里叶符号的sinc插值、奇异符号的正则化、Grünwald-Letnikov差分以及有限差分-求积方法。
• 通过缩放符号与权重渐近行为分析精度阶，特别关注其幂律衰减 $|k|^{-\alpha-1}$。
• 实施远场截断策略以处理离散算子的无限支集，确保计算可行性。
• 将该方法应用于非线性PDE，如分数阶热方程、分形Burgers方程和分数阶薄膜方程，以验证其鲁棒性与收敛性。

实验结果

研究问题

• RQ1 当应用于非光滑解时，不同分数阶拉普拉斯有限差分格式在精度、稳定性与收敛性方面如何比较？
• RQ2 分数阶拉普拉斯的谱符号与有限差分格式中所得离散权重之间存在何种关系？
• RQ3 在非局部PDE的数值求解中，如何有效截断远场边界条件而不引入显著误差？
• RQ4 解的非光滑性对分数阶PDE有限差分格式收敛速率有何影响？
• RQ5 能否为具有通用符号的平移不变非局部算子建立统一的有限差分离散化框架？

主要发现

• 离散分数阶拉普拉斯算子 $(-\Delta_h)^{\alpha/2}$ 是谱空间中的乘子，其符号由连续算子的缩放傅里叶符号导出。
• Grünwald-Letnikov权重可实现一致离散化，对 $\alpha \in (0,1)$ 具有 $O(h)$ 精度，通过改进权重可获得更高阶精度。
• 基于正则化傅里叶符号的Gorenflo-Mainardi格式在精度与稳定性方面表现更优，尤其适用于 $\alpha \in (1,2)$，且权重具有更优的渐近衰减特性。
• 数值实验表明，该格式在分数阶热方程中达到预期收敛阶，并在分形Burgers方程与薄膜方程中正确捕捉了长期行为。
• 对于非光滑初值，如 $u_0(x) = -\text{sign}(x)$，当 $\alpha < 1$ 且 $\kappa$ 较小时，解表现出类似冲击的特征，与理论预期一致。
• 该格式对分数阶薄膜方程具有鲁棒性，解在适当初始质量和参数下趋近于预期稳态 $u_\infty(x) = (1 - x^2)^{1 + \alpha/2}$。