[论文解读] Spectral Method for the Fractional Laplacian in 2D and 3D

Kailai Xu, Eric Darve|arXiv (Cornell University)|Dec 18, 2018

Numerical methods in engineering参考文献 20被引用 9

一句话总结

本文提出了一种谱方法，用于在二维和三维单位球中求解分数阶拉普拉斯算子和分数阶泊松问题，采用基于雅可比多项式构造的正交基函数，并在加权 $L^2$ 空间中进行。该方法对所有 $\alpha \in (0,2)$ 实现谱精度，且由于基函数的正交性，计算效率高，数值实验验证了其高精度和 $\mathcal{O}(\Delta t)$ 的时间收敛率。

ABSTRACT

A spectral method is considered for approximating the fractional Laplacian and solving the fractional Poisson problem in 2D and 3D unit balls. The method is based on the explicit formulation of the eigenfunctions and eigenvalues of the fractional Laplacian in the unit balls under the weighted $L^2$ space. The resulting method enjoys spectral accuracy for all fractional index $α\in (0,2)$ and is computationally efficient due to the orthogonality of the basis functions. We also proposed a numerical integration strategy for computing the coefficients. Numerical examples in 2D and 3D are shown to demonstrate the effectiveness of the proposed methods.

研究动机与目标

• 开发一种在有界区域（特别是二维和三维单位球）中求解分数阶拉普拉斯算子的高阶数值方法。
• 解决分数阶PDE解的正则性降低问题，该问题限制了传统有限差分法和有限元法的应用。
• 利用单位球中分数阶拉普拉斯算子的显式特征函数与特征值结构，构建谱精度的离散化方法。
• 通过在加权 $L^2$ 空间中使用正交基函数，确保计算效率。
• 通过二维和三维算例（包括时间依赖问题）对方法进行数值验证，展示高精度与收敛性。

提出的方法

• 该方法使用单位球中分数阶拉普拉斯算子的特征函数，表达式为 $d_{m,0}^{\alpha,d}(1-r^2)^{\alpha/2}P_m^{(\alpha/2,1/2)}(2r^2-1)$，其中 $P_m^{(\alpha/2,1/2)}$ 为雅可比多项式。
• 在具有权 $r^{d-1}$ 的加权 $L^2$ 空间中，将解展开为这些正交函数的谱基，确保谱精度。
• 采用伽辽金弱形式，通过与测试函数积分导出线性系统，其矩阵为对称正定矩阵。
• 提出一种数值积分策略，用于计算谱展开中的系数，确保离散设置下的精度。
• 对于时间依赖问题，采用隐式欧拉时间推进格式，得到形式为 $(I + \Delta t B^{-1}AD)c^{k+1} = c^k$ 的线性系统。
• 矩阵 $A$ 为基函数的格拉姆矩阵，$B$ 为对角矩阵，用于处理权重，$D$ 用于缩放基函数。

实验结果

研究问题

• RQ1 能否在二维和三维单位球中为分数阶拉普拉斯算子构造一种谱方法，使其对所有 $\alpha \in (0,2)$ 均实现谱精度？
• RQ2 如何在数值格式中有效处理分数阶PDE解的正则性降低问题？
• RQ3 能否利用加权 $L^2$ 空间中基函数的正交性，以确保计算效率与高精度？
• RQ4 在使用谱方法求解分数阶扩散方程时，时间离散化格式的时间收敛率是多少？
• RQ5 数值结果在空间和时间上与理论误差估计相比如何？

主要发现

• 数值收敛性研究证实，该谱方法对所有 $\alpha \in (0,2)$ 均在空间上实现谱精度。
• 由于基函数的正交性，该方法计算效率高，导致线性系统稀疏且条件数良好。
• 采用隐式欧拉时间推进的时间离散化格式表现出 $\mathcal{O}(\Delta t)$ 的收敛率，与理论预期一致。
• 具有径向初始条件的分数阶扩散方程的解保持径向对称，该方法能准确捕捉 $t=1$ 时刻的解分布。
• 解的误差主要由时间离散化引起，由于谱精度的存在，空间误差可忽略不计。
• 该方法对不同 $\alpha$ 值（包括 $\alpha = 0.5, 1.0, 1.5$）均表现稳健，且观察到一致的收敛速率。