[论文解读] Jacobi-Galerkin spectral method for eigenvalue problems of Riesz fractional differential equations

Lizhen Chen, Zhiping Mao+1 more|arXiv (Cornell University)|Mar 9, 2018

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

一句话总结

本文提出了一种雅可比-伽辽金谱方法，用于求解具有齐次狄利克雷边界条件的里兹分数阶微分特征值问题。通过在特殊构造的索伯列夫空间中以对称弱形式表述问题，该方法保持了矩阵的对称性与正定性，从而实现了严格的误差估计并证明了指数收敛性。关键贡献在于数值验证了 Weyl 型渐近律 $\lambda_n = \mathcal{O}(n^{2\alpha})$ 和条件数缩放 $\chi_N = \mathcal{O}(N^{4\alpha})$。通过大量数值实验，该方法实现了谱精度，并确认了理论收敛速率。

ABSTRACT

An efficient Jacobi-Galerkin spectral method for calculating eigenvalues of Riesz fractional partial differential equations with homogeneous Dirichlet boundary values is proposed in this paper. In order to retain the symmetry and positive definiteness of the discrete linear system, we introduce some properly defined Sobolev spaces and approximate the eigenvalue problem in a standard Galerkin weak formulation instead of the Petrov-Galerkin one as in literature. Poincaré and inverse inequalities are proved for the proposed Galerkin formulation which finally help us establishing a sharp estimate on the algebraic system's condition number. Rigorous error estimates of the eigenvalues and eigenvectors are then readily obtained by using Babuška and Osborn's approximation theory on self-adjoint and positive-definite eigenvalue problems. Numerical results are presented to demonstrate the accuracy and efficiency, and to validate the asymptotically exponential oder of convergence. Moreover, the Weyl-type asymptotic law $ λ_n=\mathcal{O}(n^{2α})$ for the $n$-th eigenvalue $λ_n$ of the Riesz fractional differential operator of order $2α$, and the condition number $N^{4α}$ of its algebraic system with respect to the polynomial degree $N$ are observed.

研究动机与目标

• 开发一种高阶、稳定的数值方法，用于求解具有狄利克雷边界条件的里兹分数阶微分方程特征值问题。
• 保持离散系统的对称性与正定性，确保良好的条件数和广义特征值问题的高效求解。
• 基于巴布什卡与奥斯伯恩的自伴正定特征值问题逼近理论，严格建立特征值与特征向量的误差估计。
• 通过数值方法验证 Weyl 型渐近律 $\lambda_n = \mathcal{O}(n^{2\alpha})$ 和条件数缩放 $\chi_N = \mathcal{O}(N^{4\alpha})$。
• 展示该方法在不同分数阶 $2\alpha \in (0,2)$ 及高达 $5.6$ 的高阶情况下的指数收敛速率与鲁棒性。

提出的方法

• 该方法在专门构造的索伯列夫空间中采用伽辽金弱形式，以容纳特征函数在区域边界附近的奇异性。
• 使用雅可比基函数 $\{(1-x^2)^\alpha J^{\alpha,\alpha}_n(x)\}$ 近似特征函数，以捕捉在 $x = \pm 1$ 处的 $\sim (1-x^2)^\alpha$ 奇异性。
• 该公式确保了所得代数系统的对称性与正定性，从而能够高效稳定地求解广义特征值问题。
• 为所提出的伽辽金空间导出了庞加莱不等式与逆不等式，这对控制离散系统的条件数至关重要。
• 基于巴布什卡与奥斯伯恩的理论，推导了特征值与特征向量的误差估计，适用于自伴正定特征值问题。
• 通过改变多项式阶数 $N$ 和分数阶 $2\alpha$，利用高精度求解器进行大量数值实验，以验证收敛性与渐近行为。

实验结果

研究问题

• RQ1 所提出的雅可比-伽辽金谱方法是否对里兹分数阶微分方程的特征值实现了指数收敛？
• RQ2 随着 $N$ 增大，离散系统的条件数是否准确地由 $\mathcal{O}(N^{4\alpha})$ 预测？
• RQ3 计算得到的特征值是否在 $n$ 较大时遵循 Weyl 型渐近律 $\lambda_n = \mathcal{O}(n^{2\alpha})$？
• RQ4 该方法在更高分数阶 $2\alpha > 2$（如 $3.6$ 和 $5.6$）下是否仍能保持高精度与稳定性？
• RQ5 数值特征函数在端点附近是否能有效捕捉预期的奇异性行为 $\sim (1-x^2)^\alpha$？

主要发现

• 该方法对特征值实现了渐近指数收敛，数值结果表明随着多项式阶数 $N$ 增加，误差迅速衰减。
• 特征值按 $\lambda_n = \mathcal{O}(n^{2\alpha})$ 缩放，数值结果表明在 $2\alpha \in (1.2, 2.0)$ 范围内，$\lambda_n$ 与 $n^{2\alpha}$ 呈强线性关系。
• 离散系统的条件数按 $\chi_N = \mathcal{O}(N^{4\alpha})$ 缩放，通过 $\chi_N$ 与 $N$ 的对数-对数图验证了预测的幂律关系。
• 当 $2\alpha = 2$ 时，第一特征值 $\lambda_1$ 收敛至精确值 $\left(\frac{\pi}{2}\right)^2$，且随着 $n$ 增大，$\rho_n = \lambda_n / \left(\frac{\pi n}{2}\right)^{2\alpha}$ 收敛于 1。
• 该方法在更高分数阶下（包括 $2\alpha = 3.6$ 和 $5.6$）仍保持高精度与高效性，且观察到一致的收敛行为。
• 数值结果验证了理论预测：第一特征函数 $\psi_1(x)$ 在边界附近的行为为 $\sim (1-x^2)^\alpha$，而雅可比基函数能有效捕捉该行为。