[論文レビュー] Numerical solutions of some hyperbolic stochastic partial differential equations with mixed derivatives including sine-Gordon equation

Henry C. Tuckwell|arXiv (Cornell University)|Aug 2, 2015

Stochastic processes and financial applications参考文献 28被引用数 3

ひとこと要約

本稿では、空間時間ガウス白色ノイズによって駆動される、混合微分を含む双曲型確率偏微分方程式（SPDE）を解く明示的数値スキームを提示する。この手法は、決定論的および確率的状態の両方で波動的構造を正確に捉え、特に小さなノイズ強度（σ ≤ 0.25）の場合、初期条件がゼロであってもノイズのみによっても波動的構造が出現する強固な収束性と一貫性のあるダイナミクスを示す。

ABSTRACT

We consider linear and nonlinear hyperbolic SPDEs with mixed derivatives with additive space-time Gaussian white noise of the form $Y_{xt}=F(Y) + σW_{xt}.$ Such equations, which transform to linear and nonlinear wave equations, including Klein-Gordon, Liouville's and the sine-Gordon equation, are related to what Zimmerman (1972) called a diffusion equation. An explicit numerical scheme is employed in both deterministic and stochastic examples. The scheme is checked for accuracy against known exact analytical solutions for deterministic equations. In the stochastic case with $F=0$, solutions yield sample paths for the Brownian sheet whose statistics match well exact values. Generally the boundary conditions are chosen to be initial values $Y(x,0)$ and boundary values $Y(0,t)$ on the quarter-plane or subsets thereof, which have been shown to lead to existence and uniqueness of solutions. For the linear case solutions are compared at various grid sizes and wave-like solutions were found, with and without noise, for non-zero initial and boundary conditions. Surprisingly, wave-like structures seemed to emerge with zero initial and boundary conditions and purely noise source terms with no signal. Equations considered with nonlinear $F$ included quadratic and cubic together with the sine-Gordon equation. For the latter, wave-like structures were apparent with $σ\le 0.25$ but they tended to be shattered at larger values of $σ$. Previous work on stochastic sine-Gordon equations is briefly reviewed.

研究の動機と目的

• 混合偏導関数および加法的空間時間白色ノイズを伴う双曲型SPDEを解くための数値スキームの開発と検証。
• 初期および境界条件がゼロであるがノイズが存在する場合に、波動的挙動がどのように出現するかを調査すること。
• 決定論的ケースにおける既知の解析解と比較することで、数値手法の精度および収束性を評価すること。
• ノン線形SPDE、特にサインゴルドン方程式のノイズ強度（σ）の変化に伴う挙動を調査すること。
• 座標変換を用いて、混合微分を含むSPDEの解と古典的波動方程式との関係を検討すること。

提案手法

• 空間時間白色ノイズ $ W_{xt} $ を伴う形 $ Y_{xt} = F(Y) + \sigma W_{xt} $ のSPDEを解くために、明示的有限差分スキームを適用する。
• 線形および非線形決定論的ケース（波動方程式およびリウヴィルの方程式を含む）に対して、正確な解析解と比較してスキームを検証する。
• 確率的ケースでは、ブラウン運動シートのサンプルパスをシミュレートし、統計的性質が理論的期待値と一致するかを検証する。
• 解の一意性を保証するため、定常初期および境界条件 $ Y(x,0) $ および $ Y(0,t) $ を伴う四分円領域 $[0, T] \times [0, X]$ で数値実験を実施する。
• グリッドの細分化による収束性の評価を実施し、異なるグリッドサイズ（例：400×400 と 800×800）および試行回数の違いによる解の比較を行う。
• 二次、三次、およびサインゴルドン $ F(Y) $ を含むノン線形ケースに、ノイズ強度 $ \sigma $ を系統的に変化させながら適用する。

実験結果

リサーチクエスチョン

• RQ1 混合微分および加法的空間時間白色ノイズを伴う双曲型SPDEを、明示的数値スキームが正確にシミュレートできるか？
• RQ2 初期および境界条件がゼロであるがシステムがノイズのみによって駆動される場合、波動的構造が出現するか？
• RQ3 ノイズ強度 $ \sigma $ が、確率的サインゴルドン方程式における波動的解の安定性および一貫性にどのように影響するか？
• RQ4 ノイズ強度 $ \sigma \to 0 $ の極限において、数値解が決定論的解にどの程度収束するか、およびグリッド解像度がその収束に与える影響は？
• RQ5 混合微分を含むSPDE $ Y_{xt} = F(Y) + \sigma W_{xt} $ の解と古典的波動方程式 $ \phi_{tt} - \phi_{xx} = G(\phi) $ の解との関係は？

主な発見

• 明示的数値スキームは、決定論的線形および非線形SPDEの既知の解析解を正確に再現しており、信頼性が確認された。
• 線形SPDE（$ F = 0 $）の場合、シミュレートされたブラウン運動シートのサンプルパスは理論的統計的性質と一致し、ノイズ統合の正当性が裏付けられた。
• 初期および境界条件がゼロであっても、ノイズ駆動により波動的構造が出現し、特に小さな $ \sigma $ の場合に顕著で、内因性のノイズ駆動パターン形成が示された。
• サインゴルドン方程式においては、$ \sigma \leq 0.25 $ の範囲で一貫性のある波動的構造が安定に維持されたが、ノイズ強度が高くなると構造が破壊された。
• グリッドの細分化（例：800×800 対 400×400）により、波動構造の定義が明確になり、標準偏差が低下したため、解像度の向上に伴う収束性が示された。
• $ \sigma = 0.1 $ の場合、特に800×800グリッドと50回の試行を用いた際、確率的シミュレーションの平均と決定論的解との間に良好な一致が観察された。これは、小さなノイズにおいて数値的安定性が保証されていることを示唆している。