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

ひとこと要約

本稿では、2次元および3次元単位球体において、重み付き$L^2$空間におけるジャコビ多項式から導かれる直交基底関数を用いて、分数ラプラシアンおよび分数ポisson問題を解くスペクトル法を提示する。この手法は、すべての$\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.

研究の動機と目的

• 境界領域、特に2次元および3次元単位球体における分数ラプラシアンの高次精度数値解法の開発を目的とする。
• 分数PDEの解における正則性の低下という課題に、従来の差分法および有限要素法では対応しきれないことに対処することを目的とする。
• 単位球体における分数ラプラシアンの固有関数および固有値構造を活用し、スペクトル的に正確な離散化を構築することを目的とする。
• 重み付き$L^2$空間における直交基底関数の使用を通じて、計算の効率性を保証することを目的とする。
• 2次元および3次元の例、特に時間依存問題を用いた数値的検証を通じて、高い精度と収束性を示す。

提案手法

• 本手法は、単位球体における分数ラプラシアンの固有関数を用い、$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 2次元および3次元単位球体における分数ラプラシアンに対して、すべての$\alpha \in (0,2)$においてスペクトル精度を達成するスペクトル法を構築可能か？
• RQ2 分数PDEの解における正則性の低下は、数値スキームにおいてどのように効果的に扱えるか？
• RQ3 重み付き$L^2$空間における基底関数の直交性を活用することで、計算の効率性と高精度を保証できるか？
• RQ4 スペクトル法を用いた分数拡散方程式の時間離散化スキームの時間収束率はどの程度か？
• RQ5 空間および時間方向の数値結果は、理論的誤差推定とどのように一致するか？

主な発見

• 数値収束解析により、すべての$\alpha \in (0,2)$に対して空間方向にスペクトル精度が達成されていることが確認された。
• 基底関数の直交性のおかげで、スパースかつ良好に条件付けられた線形系が得られ、計算が効率的である。
• 陰的エイラー時間積分を用いた時間離散化スキームは、理論的期待と整合的に$\mathcal{O}(\Delta t)$の収束を示した。
• 初期条件が径対称である分数拡散方程式の解は時間経過に伴い径対称を保ち、$t=1$における解のプロファイルを正確に捉えている。
• 解の誤差は時間離散化に起因しており、スペクトル精度のおかげで空間誤差は無視できるほど小さい。
• $\alpha = 0.5, 1.0, 1.5$のさまざまな値に対しても、一貫した収束率が観察され、本手法の頑健性が裏付けられた。