[論文レビュー] A PDE approach to fractional diffusion in general domains: a priori error analysis

Ricardo H. Nochetto, Enrique Otárola+1 more|arXiv (Cornell University)|Feb 4, 2013

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

ひとこと要約

本稿では、一般の有界領域における楕円型作用素の分数乗を含む分数拡散問題を解くためのPDEベースの手法を提案する。半無限円筒に問題を拡張し、重み付きソボレフ空間の枠組みを用いることで、有限要素離散化の事前誤差推定を導出する。結果として、準均等メッシュでは最適な正則性は得られるが収束順序は最適でないことが示され、非等方的メッシュでは準最適な収束順序が達成される。

ABSTRACT

The purpose of this work is the study of solution techniques for problems involving fractional powers of symmetric coercive elliptic operators in a bounded domain with Dirichlet boundary conditions. These operators can be realized as the Dirichlet to Neumann map for a degenerate/singular elliptic problem posed on a semi-infinite cylinder, which we analyze in the framework of weighted Sobolev spaces. Motivated by the rapid decay of the solution of this problem, we propose a truncation that is suitable for numerical approximation. We discretize this truncation using first degree tensor product finite elements. We derive a priori error estimates in weighted Sobolev spaces. The estimates exhibit optimal regularity but suboptimal order for quasi-uniform meshes. For anisotropic meshes, instead, they are quasi-optimal in both order and regularity. We present numerical experiments to illustrate the method's performance.

研究の動機と目的

• 非局所的作用素（例えば分数ラプラシアン）に支配される分数拡散問題を解くための堅牢な数値解法の開発。
• 分数作用素の非局所性に対処するため、拡張領域における局所的で退化／特異な楕円型問題に変換する。
• 拡張問題の有限要素離散化の収束特性を分析し、特に重み付きソボレフ空間における誤差推定に焦点を当てる。
• 準均等メッシュと非等方的メッシュの両方を用いた場合の、拡張問題の解法における収束率の性能を比較する。
• Caffarelli-Silvestre拡張フレームワークを用いて、ラプラシアンに限らず一般の2階、対称的、一様楕円型作用素へへの一般化。

提案手法

• Caffarelli-Silvestre拡張を用いて、半無限円筒 $\mathcal{C} = \Omega \times (0, \infty)$ 上のディリクレ・トゥ・ノイマン写像問題として分数拡散問題を定式化する。
• 重み $y^\alpha$, $\alpha = 1 - 2s$ を含む退化／特異な楕円型PDEとして拡張問題を定式化し、重み付きソボレフ空間で扱う。
• 解の $y$ 方向における指数的減衰を利用するため、無限円筒を有限高さ $\mathcal{C}_\mathcal{Y}$ で切断して数値近似を可能にする。
• 拡張領域の形状正則な三角形分割上での1次テンソル積有限要素を用いて、切断問題を離散化する。
• 拡張問題の重み付き $H^1$-ノルムにおける事前誤差推定を導出し、誤差がメッシュサイズおよび要素数に依存することを示す。
• 準均等および非等方的メッシュの両方における収束挙動を分析し、非等方的メッシュでは収束率が向上することを示す。

実験結果

リサーチクエスチョン

• RQ1 一般の有界領域における楕円型作用素の分数乗を、局所的PDE定式化を用いてどのように数値的に近似できるか？
• RQ2 準均等メッシュと非等方的メッシュを用いた場合、拡張問題の有限要素近似の収束挙動はどのように異なるか？
• RQ3 半無限円筒の切断が数値解の精度に及ぼす影響はどの程度か？
• RQ4 この手法によって、元の分数拡散問題の解において、最適な正則性と準最適な収束順序を達成できるか？
• RQ5 誤差推定は、基礎となる楕円型作用素の固有値特性および重み $y^\alpha$ にどのように依存するか？

主な発見

• 拡張問題において、誤差推定は最適な正則性を達成しており、解は重み付き $H^1$-空間 $\widetilde{H}^1_L(\mathcal{C}, y^\alpha)$ に属する。
• 準均等メッシュでは収束順序が最適でないことが示され、$|\log(\#\mathscr{T}_\mathcal{Y})|^s (\#\mathscr{T}_\mathcal{Y})^{-1/(n+1)}$ のスケーリングに従う。
• 非等方的メッシュで $\gamma > 3/(1 - \alpha)$ を満たす場合、誤差推定は順序および正則性の両面で準最適となり、最良の可能な収束レートに一致する。
• 切断誤差は、切断高さ $\mathcal{Y}$ に従って指数的に減少し、$\|\nabla(\mathfrak{u} - v)\|_{L^2(\mathcal{C}, y^\alpha)} \lesssim e^{-\sqrt{\lambda_1}\mathcal{Y}/4}\|f\|_{\mathbb{H}^s(\Omega)'}$ と評価できる。
• 有限要素解 $V_{\mathscr{T}_\mathcal{Y}}$ は、重み付きエネルギーノルムにおいて真の解 $\mathfrak{u}$ に収束し、誤差はメッシュ依存項によって上界が与えられる。
• この手法は、ラプラシアンに限らず、任意の2階、対称的、一様楕円型作用素 $\mathcal{L}$ に対しても、一般化されたCaffarelli-Silvestre拡張を用いて一般化可能である。