[論文レビュー] On a space discretization scheme for the Fractional Stochastic Heat Equations

Latifa Debbi, Marco Dozzi|arXiv (Cornell University)|Feb 23, 2011

Fractional Differential Equations Solutions参考文献 25被引用数 5

ひとこと要約

本稿は、乗法的円筒形ウィナー過程によって駆動される分数確率的熱方程式のための、非局所的分数ラプラシアン作用素を離散的ラプラシアン作用素の分数乗によって近似する新しい空間離散化スキームを提案する。著者らは収束速度を確立し、収束を保証する分数階数 $\alpha$ の臨界値を特定し、非局所的作用素を伴う方程式のための厳密な数値フレームワークを提供する。

ABSTRACT

In this work, we introduce a new discretization to the fractional Laplacian and use it to elaborate an approximation scheme for fractional heat equations perturbed by a multiplicative cylindrical white noise. In particular, we estimate the rate of convergence.

研究の動機と目的

• 乗法的ノイズを伴う分数確率的熱方程式のための新しい空間離散化スキームを開発すること。
• 確率的偏微分方程式における非局所的分数ラプラシアン作用素の数値的近似の課題に取り組むこと。
• 提案されたスキームの収束速度を推定し、収束を保証する分数階数 $\alpha$ の臨界値を同定すること。
• ラプラシアンの離散的分数乗を用いた分数作用素の数値的近似のための厳密な数学的基盤を提供すること。

提案手法

• 均等グリッド上での有限差分法を用いてラプラシアン作用素を離散化し、その後その離散作用素を分数乗 $\alpha/2$ に上げることで分数ラプラシアンを近似する。
• 積分公式 $A_{\alpha} = \frac{\sin(\alpha\pi/2)}{\pi}\int_0^\infty z^{\alpha/2 - 1} A(Iz + A)^{-1} dz$ を用いて分数ラプラシアンのスペクトル表現を導出する。
• 無限次元の円筒形ウィナー過程を有限次元部分空間に射影することにより、確率的熱方程式の有限次元近似を構築する。
• 離散初期データおよびノイズ項に行列指数関数 $e^{-t\mathbb{A}_{{\mathcal{V}}_{n-1}}^{\alpha/2}}$ を適用することで、離散解を定義する。
• 有限次元における伊藤積分を用いて、確率的畳み込み項を $\sum_{j=1}^{n-1} \int_0^t E_n e^{-\mathbb{A}_{{\mathcal{V}}_{n-1}}^{\alpha/2}(t-s)} P_n g(u^n(s)) \, dB_j(s)$ と表現する。
• 補間理論および作用素ノルムの見積もりを用いて、ソボレフ空間における近似誤差を制御し、収束速度を導出する。

実験結果

リサーチクエスチョン

• RQ1 提案された有限差分スキームが分数確率的熱方程式に適用された際の最適収束速度は何か？
• RQ2 収束はラプラシアンの分数階数 $\alpha$ にどのように依存するか？
• RQ3 $L^2$-ノルムにおける数値スキームの収束を保証するための $\alpha$ の臨界値は何か？
• RQ4 離散ラプラシアンを分数乗に上げることで分数ラプラシアンを効果的に近似できるか？理論的根拠は何か？
• RQ5 円筒形ウィナー過程によって駆動される確率的項は、スキームの収束挙動にどのように影響を与えるか？

主な発見

• 提案されたスキームは、初期データの正則性および分数階数 $\alpha$ に依存するが、$O(n^{-\delta})$ の収束速度（$\delta > 0$）を達成する。$\alpha$ に明示的な依存関係が存在する。
• 収束は $\alpha > 1$ の場合に保証され、臨界閾値 $\alpha = 1$ は収束が失敗するか遅延する可能性がある境界ケースとして特定される。
• 離散的分数ラプラシアンは、標準ラプラシアンの有限差分離散化の分数乗として定義され、一貫性と安定性を有する近似を提供する。
• 有限次元設定における解は行列指数関数と伊藤積分を用いて構築され、適切に定義されており、パスワイズ収束を保証する。
• 誤差見積もりは補間理論および作用素ノルムの境界を用いて導出され、$\mathbb{E}[\sup_{t\in[0,T]} \|u^n(t) - u(t)\|_{L^2}^p] \leq C_{T,n,\|g\|}(1 + \mathbb{E}[\|u_0\|_{L^2}^p])$ が成り立つ。ここで $C$ は $n$ およびノイズ係数 $g$ に依存する。
• 離散作用素 $\mathbb{A}_{{\mathcal{V}}_{n-1}}^{\alpha/2}$ の固有関数が連続的分数ラプラシアンの固有関数と一致することが示され、スペクトル的一致性が保証される。