[論文レビュー] 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

ひとこと要約

本稿では、同次Dirichlet境界条件を満たすRiesz分数階微分固有値問題を解くためのジャコビ・ガレルキンスペクトル法を提案する。特に定義されたソボレフ空間内での対称的弱形式として問題を定式化することにより、行列の対称性と正定値性を保ち、厳密な誤差推定と指数的収束の証明が可能になる。主な貢献は、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.

研究の動機と目的

• Riesz分数階微分方程式の固有値問題に対して高次の安定な数値解法を開発すること。
• 離散系の対称性と正定値性を保つことにより、良好な条件数と代数的固有値問題の効率的かつ安定な解法を確保すること。
• 自己随伴かつ正定値固有値問題に対するBabuškaおよびOsbornの近似理論を用いて、固有値および固有関数の誤差推定を厳密に確立すること。
• 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$ の特異性を捉える。
• この定式化により、得られる代数的系の対称性と正定値性が保証され、一般化固有値問題の効率的かつ安定な解法が可能になる。
• 提案されたガレルキン空間に対して、Poincaré不等式および逆不等式を導出する。これらは離散系の条件数の上界を求める上で不可欠である。
• 自己随伴かつ正定値固有値問題に対するBabuškaおよびOsbornの理論を用いて、固有値および固有関数の誤差推定を導出する。
• 多項式次数 $N$ と分数階数 $2\alpha$ を変化させた数値実験を実施し、高精度ソルバーを用いて収束性および漸近的挙動の妥当性を検証する。

実験結果

リサーチクエスチョン

• RQ1 提案されたジャコビ・ガレルキンスペクトル法は、Riesz分数階微分方程式の固有値に対して指数的収束を達成するか？
• RQ2 離散系の条件数は、$N$ の増加に伴い $\mathcal{O}(N^{4\alpha})$ として正確に予測されるか？
• RQ3 計算された固有値は、$n$ が大きい場合に $\lambda_n = \mathcal{O}(n^{2\alpha})$ のWeyl型漸近法則に従うか？
• 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$ のように振る舞うという理論的予測を裏付けており、ジャコビ基底がこれを効果的に捉えていることが確認された。