[論文レビュー] Stochastic Solution of Elliptic and Parabolic Boundary Value Problems for the Spectral Fractional Laplacian

Mamikon Gulian, Guofei Pang|arXiv (Cornell University)|Dec 4, 2018

Advanced Mathematical Modeling in Engineering参考文献 23被引用数 5

ひとこと要約

本稿では、非ゼロのディリクレ境界条件を伴うスペクトル的分数ラプラシアン $(-\Delta_{\Omega,g})^{\alpha/2}$ に関する楕円型および放物型境界値問題に対する、新規の確率的解公式（ファインマン＝カック型）を提示する。この手法は、従属停止ブラウン運動を用い、モンテカルロシミュレーションを介した確率的フレームワークを提供することで、有界領域における効率的かつ局所的に並列化可能な計算を可能にし、収束性および正則性解析が厳密に行われる。

ABSTRACT

We prove and implement stochastic solution (or Feynman-Kac) formulas for boundary value problems involving the spectral fractional Laplacian with nonzero Dirichlet boundary condition. The main tools used in the proofs are the abstract Cauchy problem for Feller semigroups together with Balakrishnan's theory of fractional powers. We show the the spectral fractional Laplacian with nonzero Dirichlet boundary conditions is the generator of an appropriate Feller semigroup for subordinate stopped Brownian motion, and obtain a stochastic solution formula for the fractional heat equation in a bounded domain. We then obtain precise regularity and steady-state convergence properties of the parabolic problem using the eigenfunction expansion of the classical solution, which leads to estimates for the survival probability of subordinate stopped Brownian motion. These results allow us to take the long-time limit in the parabolic formula to establish a stochastic solution formula for the Dirichlet boundary value problem. These stochastic solution formulas for the spectral fractional Laplacian with nonzero boundary conditions are novel, and allow for efficient, embarrassingly parallel local solution of the above boundary value problems. We discuss the discretization of these formulas, and verify them in dimensions two and three with benchmark examples. We study the effect of the number of path samples and the time step size for path discretization on the accuracy of the solution.

研究の動機と目的

• 非ゼロのディリクレ境界条件を伴うスペクトル的分数ラプラシアンに対する確率的解公式の開発および検証。
• フェラー半群および従属理論を用いたスペクトル的分数ラプラシアンの確率的解釈の確立。
• 経路に基づくモンテカルロサンプリングを用いた、分数PDEの効率的かつ明白に並列化可能な数値的解法の実現。
• 固有関数展開を用いた、分数熱方程式の正則性および定常収束性の解析。
• 2次元および3次元領域における数値ベンチマークを通じた提案公式の検証。

提案手法

• バラクリシュナンの分数乗の理論を用いて、スペクトル的分数ラプラシアンをフェラー半群の生成作用素として導出。
• プロセス $X^{\Omega,\alpha}_{t} = \sqrt{2}B_{T_{\alpha/2}(t) \wedge \tau_{\Omega}}$ の構築、これは $\alpha/2$-安定な従属過程で時間変更されたブラウン運動であり、領域の境界に達した時点で停止する。
• 非同次分数カウチ問題におけるファインマン＝カック公式の確立：$u(t,x) = \mathbb{E}_x\left[ f(X_t) \chi_{\tau_\Omega > T_{\alpha/2}(t)} + g(X_t) \chi_{\tau_\Omega \leq T_{\alpha/2}(t)} + \int_0^t r(t-s,X_s) \chi_{\tau_\Omega > T_{\alpha/2}(s)} ds \right]$。
• 定常状態のディリクレ問題への解の拡張：$u(x) = \mathbb{E}_x\left[ g(X_{T^{-1}_{\alpha/2}(\tau_\Omega)}) \right] + \mathbb{E}_x\left[ \int_0^{T^{-1}_{\alpha/2}(\tau_\Omega)} r(X_s) ds \right]$。
• 固有関数展開を用いて、従属停止過程の生存確率および正則性の推定を導出。
• 経路サンプリングを用いた数値的実装および検証を、2次元および3次元領域で行い、経路数および時間刻みの大きさに関する収束解析を実施。

実験結果

リサーチクエスチョン

• RQ1 非ゼロのディリクレ境界条件を伴うスペクトル的分数ラプラシアンは、どのように確率的過程として表現可能か？
• RQ2 従属停止ブラウン運動と分数放物型・楕円型PDEの解との間にはどのような関係があるか？
• RQ3 従属停止過程の生存確率は、解の正則性および収束性にどのように影響するか？
• RQ4 確率的解公式は、保証された収束性および並列スケーラビリティを備えて、数値的に実装可能か？
• RQ5 経路数および時間刻みの大きさに関して、モンテカルロ離散化の精度および収束特性はいかなるものか？

主な発見

• スペクトル的分数ラプラシアン $(-\Delta_{\Omega,g})^{\alpha/2}$ が、従属停止ブラウン運動に関連するフェラー半群の生成作用素として示された。
• 抽象的カウチ問題および従属理論を用いて、分数熱方程式の確率的解公式が導出され、証明された。
• 定常状態のディリクレ問題の解は、放物型公式における $t \to \infty$ の極限をとることで得られ、プロセスの領域外への抜ける時間に依存する経路ベースの期待値として表現された。
• 2次元および3次元における数値実験では、経路数の増加および時間刻みの減少に伴い収束が確認され、手法の精度が検証された。
• 本手法は、局所的に並列化可能でグリッドフリーな分数PDEの解法を可能にし、高次元においても良好なスケーリングを示した。
• 固有関数展開を用いて正則性および生存確率の推定が得られ、解フレームワークの適切な定式化が裏付けられた。