[論文レビュー] Numerical approximations for fractional elliptic equations via the method of semigroups

Nicole Cusimano, Félix del Teso+1 more|arXiv (Cornell University)|Dec 4, 2018

Differential Equations and Boundary Problems被引用数 7

ひとこと要約

本稿では、有界領域におけるスペクトル分数ラプラシアン $(-\Delta_{\mathcal{B}})^s$ を含む分数級楕円型方程式を解くための新しい数値的手法を提示する。半群アプローチを用い、熱半群の有限要素離散化と特異積分の単調な求積則を組み合わせることで、滑らかなデータでは $O(h^{4})$、$L^2$ データでは $O(h^{2s})$ の収束速度を達成し、ディリクレ、ノイマン、ロビン境界条件に対して高い精度と柔軟性を提供する。

ABSTRACT

We provide a novel approach to the numerical solution of the family of nonlocal elliptic equations $(-\\Delta)^su=f$ in $\\Omega$, subject to some homogeneous boundary conditions $\\mathcal{B}(u)=0$ on $\\partial \\Omega$, where $s\\in(0,1)$, $\\Omega\\subset \\mathbb{R}^n$ is a bounded domain, and $(-\\Delta)^s$ is the spectral fractional Laplacian associated to $\\mathcal{B}$ on $\\partial \\Omega$. We use the solution representation $(-\\Delta)^{-s}f$ together with its singular integral expression given by the method of semigroups. By combining finite element discretizations for the heat semigroup with monotone quadratures for the singular integral we obtain accurate numerical solutions. Roughly speaking, given a datum $f$ in a suitable fractional Sobolev space of order $r\\geq 0$ and the discretization parameter $h>0$, our numerical scheme converges as $O(h^{r+2s})$, providing super quadratic convergence rates up to $O(h^4)$ for sufficiently regular data, or simply $O(h^{2s})$ for merely $f\\in L^2(\\Omega)$. We also extend the proposed framework to the case of nonhomogeneous boundary conditions and support our results with some illustrative numerical tests.

研究の動機と目的

• 有界領域におけるスペクトル分数ラプラシアンを用いた非局所楕円型問題を高精度かつ高効率に解く数値スキームの開発。
• $(-\Delta_{\mathcal{B}})^s u = f$ を同次および非同次境界条件で解く課題に対処すること。
• 逆分数ラプラシアンの半群表現を活用して、安定かつ高次精度の離散化を可能とすること。
• データ $f$ の正則性仮定が異なる場合の、提案手法の厳密な収束速度の確立。

提案手法

• 解は半群公式 $(-\Delta_{\mathcal{B}})^{-s}f(x) = \frac{1}{\Gamma(s)}\int_0^\infty e^{t\Delta_{\mathcal{B}}}f(x) \frac{dt}{t^{1-s}}$ で表現され、分数級逆ラプラシアンと熱半群を結びつける。
• 有限要素法を用いて熱半群 $e^{t\Delta_{\mathcal{B}}}f$ を離散化し、各時刻 $t$ における解の空間的近似を実現する。
• 時間方向の特異積分を単調な求積則で近似し、安定性と高次収束性を確保する。
• 解を同次問題と補正項に分解することで、非同次境界条件への拡張を実現する。
• 有限要素誤差見積もりと求積誤差バウンドを組み合わせた収束解析により、グローバルな収束速度を導出する。
• ディリクレ、ノイマン、ロビン境界条件に適用可能であり、より一般的な作用素への拡張も可能である。

実験結果

リサーチクエスチョン

• RQ1 有界領域におけるスペクトル分数ラプラシアンに対して、半群法を用いて高次精度の数値スキームを導出できるか？
• RQ2 半群表現における時間積分に高次求積則を適用し、空間的有限要素離散化と組み合わせた場合、得られる収束速度は何か？
• RQ3 $L^2(\Omega)$ に属する低正則性のデータに対して、この手法はどのように動作し、収束性は保たれるか？
• RQ4 非同次境界条件に拡張しても、最適な収束速度を維持できるか？
• RQ5 非凸多角形領域でもこのアプローチは頑健か？特異性付近でのメッシュの細分化により、最適収束が回復可能か？

主な発見

• データ $f \in \mathbb{H}_{\mathcal{B}}^r(\Omega)$ に対して、$O(h^{r+2s})$ の収束速度を達成し、十分に滑らかなデータでは最大 $O(h^4)$ の収束が得られる。
• $f \in L^2(\Omega)$ の場合、$O(h^{2s})$ の収束速度を示し、与えられた正則性において最適である。
• 熱半群の有限要素近似と単調な求積則の組み合わせにより、安定性と高次精度が保証される。
• 解を調和補正項と同次問題に分解することで、非同次境界条件への拡張が実現される。
• 理論的収束結果は、さまざまな境界条件とデータ正則性において予測された収束速度を示す数値実験によって裏付けられる。
• このフレームワークは、他の2階楕円型作用素や境界条件へも一般化可能であり、求積則および有限要素部の最小限の修正で対応可能である。