[論文レビュー] Finite difference methods for fractional Laplacians

Yanghong Huang, Adam M. Oberman|arXiv (Cornell University)|Nov 1, 2016

Fractional Differential Equations Solutions参考文献 33被引用数 18

ひとこと要約

本稿では、1次元の均等グリッド上で、スペクトル表現とフーリエ解析を用いて導出された重みを用いた離散畳み込みに基づく、分数ラプラシアン $(-\Delta)^{\alpha/2}$ の統一された有限差分フレームワークを提示する。主な貢献は、グリューンワルド＝レトニコフ、ゴレンフロ＝マイナーディ、および数値積分に基づく手法を含む、さまざまな離散化スキームの体系的比較であり、滑らかでない解や遠方境界条件下での精度、安定性、収束特性を示している。

ABSTRACT

The fractional Laplacian $(-Δ)^{α/2}$ is the prototypical non-local elliptic operator. While analytical theory has been advanced and understood for some time, there remain many open problems in the numerical analysis of the operator. In this article, we study several different finite difference discretisations of the fractional Laplacian on uniform grids in one dimension that takes the same form. Many properties can be compared and summarised in this relatively simple setting, to tackle more important questions like the nonlocality, singularity and flat tails common in practical implementations. The accuracy and the asymptotic behaviours of the methods are also studied, together with treatment of the far field boundary conditions, providing a unified perspective on the further development of the scheme in higher dimensions.

研究の動機と目的

• 1次元における分数ラプラシアン $(-\Delta)^{\alpha/2}$ の有限差分離散化の統一的数値フレームワークの構築を目的とする。
• 分数ラプラシアンの異なる定義に基づく、既存および新規の離散化スキームの精度、安定性、漸近的挙動を比較することを目的とする。
• 非局所性、特異的重み、遠方境界条件といった数値実装上の実用的課題に対処することを目的とする。
• 滑らかでない解を伴うPDEに対して、分数ラプラシアン作用素を有する場合の収束性と精度を分析することを目的とする。
• 一般の記号を持つ非局所作用素へとこのスキームを高次元および一般化する基盤を提供することを目的とする。

提案手法

• 離散分数ラプラシアンを畳み込み作用素として定式化：$(-\Delta_h)^{\alpha/2}u_j = \sum_{k=-\infty}^{\infty} (u_j - u_{j-k}) w_k$ であり、重み $w_k$ はスペクトル表現から導出される。
• 半離散的フーリエ解析を用いて、離散作用素とその記号を関連づけ、スペクトル解析と安定性評価を可能にする。
• 重みを複数のアプローチで導出：フーリエ記号のシンク補間、特異的記号の正則化、グリューンワルド＝レトニコフ差分、有限差分-数値積分法。
• スケーリングされた記号と重みの漸近的挙動を用いて、精度の順序を分析し、特に重みのべき乗則的減衰 $|k|^{-\alpha-1}$ に注目する。
• 離散作用素の無限大支持を扱うための遠方切り捨て戦略を実装し、計算の実行可能性を確保する。
• 非線形PDE、例えば分数熱方程式、フラクタルバーガース方程式、分数薄膜方程式への適用を通じて、スケームの頑健性と収束性を検証する。

実験結果

リサーチクエスチョン

• RQ1 滑らかでない解に適用した場合、分数ラプラシアンの異なる有限差分スキームは、精度、安定性、収束性においてどのように比較されるか？
• RQ2 分数ラプラシアンのスペクトル記号と、有限差分スキームにおける結果得られる離散的重みの間にはどのような関係があるか？
• RQ3 非局所PDEの数値解法において、遠方境界条件をどのように切り捨てれば、顕著な誤差を生じさせずに効果的に処理できるか？
• RQ4 解の滑らかでない性質が、分数PDEの有限差分スキームの収束速度に及ぼす影響は何か？
• RQ5 一般の記号を持つ並進不変非局所作用素の有限差分離散化に対して、統一されたフレームワークを確立できるか？

主な発見

• 離散分数ラプラシアン $(-\Delta_h)^{\alpha/2}$ はスペクトル空間における乗数であり、連続作用素のスケーリングされたフーリエ記号から導出される。
• グリューンワルド＝レトニコフ重みは、$\alpha \in (0,1)$ の場合に $O(h)$ の精度で一貫した離散化を実現し、修正された重みを用いることで高次の精度が達成可能である。
• 正則化されたフーリエ記号に基づくゴレンフロ＝マイナーディスキームは、$\alpha \in (1,2)$ の場合に特に優れた精度と安定性を示し、重みの漸近的減衰が改善されている。
• 数値実験により、分数熱方程式に対して期待通りの収束順序が得られ、フラクタルバーガース方程式および薄膜方程式では正しい長期挙動が正しく捉えられている。
• 初期データが滑らかでない場合、例えば $u_0(x) = -\text{sign}(x)$ の場合、$\alpha < 1$ かつ $\kappa$ が小さいとき、解はショックのような挙動を示し、理論的予測と整合的である。
• 分数薄膜方程式に対してもスケームは頑健であり、適切な初期質量とパrametersのもとで、解は期待される定常状態 $u_\infty(x) = (1 - x^2)^{1 + \alpha/2}$ に収束する。