[论文解读] A simple solver for the fractional Laplacian in multiple dimensions

Victor Minden, Lexing Ying|arXiv (Cornell University)|Feb 11, 2018

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

一句话总结

本论文提出了一种在多维空间中对分数阶拉普拉斯算子进行简单、二阶精度离散化的数值方法，结合奇点减去法与规则网格上的梯形法则，通过快速傅里叶变换（FFT）实现高效计算。该方法基于标准拉普拉斯算子的预条件Krylov求解器，高效求解分数阶扩散问题，展现出最优收敛性与鲁棒性，即使当分数阶参数 α 接近 2 时亦表现良好。

ABSTRACT

We present a simple discretization scheme for the hypersingular integral representation of the fractional Laplace operator and solver for the corresponding fractional Laplacian problem. Through singularity subtraction, we obtain a regularized integrand that is amenable to the trapezoidal rule with equispaced nodes, assuming a high degree of regularity in the underlying function (i.e., $u\in C^6(R^d)$). The resulting quadrature scheme gives a discrete operator on a regular grid that is translation-invariant and thus can be applied quickly with the fast Fourier transform. For discretizations of problems related to space-fractional diffusion on bounded domains, we observe that the underlying linear system can be efficiently solved via preconditioned Krylov methods with a preconditioner based on the finite-difference (non-fractional) Laplacian. We show numerical results illustrating the error of our simple scheme as well the efficiency of our preconditioning approach, both for the elliptic (steady-state) fractional diffusion problem and the time-dependent problem.

研究动机与目标

• 开发一种简单、精确且高效的多维分数阶拉普拉斯算子离散化格式。
• 在具有狄利克雷边界条件的有界区域上，实现空间分数阶扩散问题的快速求解。
• 通过引入有效的预条件策略，缓解当 α → 2 时离散分数阶拉普拉斯系统病态的问题。
• 将简单求积格式的适用范围扩展至解的正则性有限、特别是靠近区域边界的情形。

提出的方法

• 通过泰勒展开实现奇点减去法，对分数阶拉普拉斯算子的超奇异积分核进行正则化处理。
• 在等距结构化笛卡尔网格上，使用梯形法则对正则化后的被积函数进行离散化，得到一个平移不变的离散算子。
• 通过快速傅里叶变换（FFT）高效应用该离散算子，实现 O(N log N) 的计算复杂度。
• 为求解所得线性系统，采用基于标准有限差分拉普拉斯算子的预条件器，并结合Krylov子空间方法。
• 该预条件器利用了离散分数阶拉普拉斯算子可分解为标准拉普拉斯算子加上低秩或稀疏校正矩阵的结构特性。
• 通过移除对应于区域补集的网格点，将该方法扩展至遮蔽的笛卡尔网格，如 L 形区域。

实验结果

研究问题

• RQ1 能否基于奇点减去法，为多维空间中的分数阶拉普拉斯算子开发一种简单且高阶的求积格式？
• RQ2 所得离散算子是否保持平移不变性，从而支持基于FFT的高效计算？
• RQ3 能否通过基于标准拉普拉斯算子的预条件器，有效缓解当 α 接近 2 时离散分数阶拉普拉斯系统的病态性？
• RQ4 对于正则性有限的解，特别是靠近区域边界时，该方法的收敛速率如何？
• RQ5 预条件Krylov求解器在不同维度与分数阶参数 α 下的性能表现如何？

主要发现

• 对于足够光滑的函数（u ∈ C⁶(Rᵈ)），该离散化方法在网格间距 h 上达到二阶收敛率。
• 对于正则性较差的解（通常在有界区域中出现），方法在 ℓ² 范数下表现出一阶收敛性。
• 基于有限差分拉普拉斯算子的预条件器显著减少了Krylov迭代次数，尤其在较大 α 值时效果明显。
• 在二维情况下，带预条件的CG方法在 α = 1.25 至 1.75 范围内仅需 12–29 次迭代，运行时间呈 O(N log N) 量级。
• 在三维情况下，方法依然稳健，迭代次数随 α = 1.75 增加至 30–170，但预条件策略仍具有效性。
• 误差估计的理查德森外推率约为 1.0–1.5，表明在测试网格上实现了最优收敛性。