[论文解读] Algorithm implementation and numerical analysis for the two-dimensional tempered fractional Laplacian

Jing Sun, Daxin Nie+1 more|arXiv (Cornell University)|Feb 7, 2018

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

一句话总结

本文提出了一种结合有限差分格式、双线性插值与加权梯形求积法的二维广义分数阶拉普拉斯算子（∆+λ）β/2的数值方法。该方法对C²解实现了O(h²⁻β)的收敛速率，并通过利用对称块托普利茨系统的FFT加速共轭梯度求解器，高效求解广义分数阶泊松方程，将计算成本从O(N⁶)降低至O(N² log N²)。

ABSTRACT

Tempered fractional Laplacian is the generator of the tempered isotropic L\'evy process [W.H. Deng, B.Y. Li, W.Y. Tian, and P.W. Zhang, Multiscale Model. Simul., 16(1), 125-149, 2018]. This paper provides the finite difference discretization for the two dimensional tempered fractional Laplacian $(\Delta+\lambda)^{\frac{\beta}{2}}$. Then we use it to solve the tempered fractional Poisson equation with Dirichlet boundary conditions and derive the error estimates. Numerical experiments verify the convergence rates and effectiveness of the schemes.

研究动机与目标

• 开发一种针对二维广义分数阶拉普拉斯算子（∆+λ）β/2的鲁棒数值离散化方法。
• 解决分数阶PDE中非局部性、弱奇异性以及高维计算带来的挑战。
• 为带狄利克雷边界条件的广义分数阶泊松方程应用有限差分格式，提供严格的误差估计。
• 设计一种高效迭代求解器，利用刚度矩阵的对称块托普利茨结构。

提出的方法

• 使用加权梯形求积法近似广义分数阶拉普拉斯算子的奇异积分表示。
• 应用双线性插值将求积方法扩展至局部 stencils 之外的全局计算区域。
• 将R²上的积分变换为对称象限，并利用极坐标处理奇异性。
• 通过核函数的数值积分推导离散权重，对中心点(i,j)=(0,0)采用极坐标进行特殊处理。
• 利用快速傅里叶变换（FFT）加速共轭梯度法中的矩阵-向量乘法。
• 通过利用块托普利茨结构，将内存占用从O(N⁴)降低至O(N²)，计算成本从O(N⁶)降低至O(N² log N²)。

实验结果

研究问题

• RQ1 能否为二维广义分数阶拉普拉斯算子构造一种具有可证明收敛性的高阶有限差分格式？
• RQ2 在二维情形下，如何高效处理该算子的非局部性与弱奇异性？
• RQ3 对于光滑解，所提格式的收敛速率是多少？
• RQ4 在实际中，如何高效求解由离散化产生的稠密线性系统？
• RQ5 该方法能否扩展至处理标准分数阶拉普拉斯算子（λ=0）及无界区域？

主要发现

• 所提出的有限差分格式对C²(R²)类解实现了O(h²⁻β)的收敛速率，验证了在给定正则性下的最优收敛阶。
• 误差估计基于精确解在区域闭包上属于C²的假设，而非在整个空间上。
• 离散化产生的刚度矩阵为对称块托普利茨矩阵，其块结构为高效求解策略提供了可能。
• 利用FFT将共轭梯度迭代中矩阵-向量乘法的计算成本从O(N⁶)降低至O(N² log N²)。
• 该方法适用于广义分数阶拉普拉斯算子（λ>0）与标准分数阶拉普拉斯算子（λ=0），相应地调整数值积分技术。
• 数值实验验证了预测的收敛速率，并证实了该格式与求解器在具有与不具有精确解问题上的有效性。