[論文レビュー] Energy-preserving exponential integrable numerical method for stochastic cubic wave equation with additive noise

Jianbo Cui, Jialin Hong+2 more|arXiv (Cornell University)|Sep 2, 2019

Stochastic processes and financial applications参考文献 34被引用数 5

ひとこと要約

本稿では、加法的ノイズを伴う確率的三次波動方程式に対する、エネルギー保存性と指数的可積分性を併せ持つ新しい数値解法を提案する。空間的分離にはスペクトルガレルキン法を、時間的分離には分割技法と平均化ベクトル場法を組み合わせることで、解のエネルギー保存則と指数的可積分性を保持する。強収束率は $ O(h^2 + \epsilon^8) $ であり、$ h $ は時間刻み、$ \epsilon $ は切断誤差を制御するパラメータである。

ABSTRACT

In this paper, we present an energy-preserving exponentially integrable numerical method for stochastic wave equation with cubic nonlinearity and additive noise. We first apply the spectral Galerkin method to discretizing the original equation and show that this spatial discretization possesses an energy evolution law and certain exponential integrability property. Then the exponential integrability property of the exact solution is deduced by proving the strong convergence of the semi-discretization. To propose a full discrete numerical method which could inherit both the energy evolution law and the exponential integrability, we use the splitting technique and averaged vector field method in the temporal direction. Combining these structure-preserving properties with regularity estimates of the exact and the numerical solutions, we obtain the strong convergence rate of the numerical method. Numerical experiments coincide with these theoretical results.

研究の動機と目的

• 加法的ノイズを伴う確率的三次波動方程式のエネルギー保存則と指数的可積分性を両方保持する完全離散数値スキームの構築を目的とする。
• 空間的スペクトルガレルキン離散化と時間的分割法に平均化ベクトル場法を組み合わせた仮定の下で、数値解法の強収束性を確立することを目的とする。
• スペクトルガレルキン近似の一様な指数的可積分性と正則性推定を活用して、正確な解の指数的可積分性を証明することを目的とする。
• 特に、内在的な物理的性質を保存する構造保存型数値解法のギャップを埋め、三次非線形性を有する確率的波動方程式に対する手法の発展を目的とする。
• 理論的予測と整合する数値実験を通じて、理論的収束率の妥当性を検証することを目的とする。

提案手法

• 空間的離散化には、エネルギー保存則を保持し、半離散系において指数的可積分性を示すスペクトルガレルキン法を用いる。
• 時間的離散化には、エネルギー保存性と時間的指数的可積分性を維持するため、分割技法と平均化ベクトル場法を組み合わせる。
• 本手法は、正確な解の両方の性質（エネルギー保存則と指数的可積分性）を継承するように構築されており、長期的な安定性と精度を保証する。
• 収束解析を支援するため、正確な解および数値解のための事前推定と正則性バウンドを導出する。
• 確率的偏差の大きさを制御するために、ガウス尾部推定とソボレフ埋め込み定理を用いる。
• 反復的誤差推定を適用して最終的な強収束率を導出し、誤差は $ \mathbb{E}[\|u^* - u^{N,k+1}\|_{L^2}^2] \leq C \max(h^2, \epsilon^8) \left( (1/2)^{k+1} + \exp(-\eta / \epsilon^2) \right) $ で抑えられる。

実験結果

リサーチクエスチョン

• RQ1 加法的ノイズを伴う確率的三次波動方程式に対して、エネルギー保存則と指数的可積分性を両方保持する完全離散数値スキームを構築可能か？
• RQ2 適切な正則性および可積分性仮定の下で、このような構造保存型スキームの強収束率は何か？
• RQ3 確率的三次波動方程式の正確な解は指数的可積分性を有するか？ また、その性質は数値解法に継承可能か？
• RQ4 スペクトルガレルキン法と分割技法は、方程式の内在的な幾何学的・確率的構造をどのように保存するか？
• RQ5 提案手法は、最適な収束率を達成しつつ、長期的安定性とエネルギー保存性を維持できるか？

主な発見

• 加法的ノイズを伴う確率的三次波動方程式の正確な解は、指数的可積分性を満たす：$ \mathbb{E}\left[\exp\left(\int_0^T c\|u(s)\|_{L^6}^2 ds\right)\right] \leq C $ であり、$ C $ は初期データ、ノイズ、時間、次元に依存する。
• スペクトルガレルキン半離散化は、エネルギー保存則と指数的可積分性の両方を保持し、完全離散化の強収束の導出を可能にする。
• 空間にスペクトルガレルキン法、時間に分割-平均化ベクトル場法を組み合わせた本稿で提案する完全離散法は、エネルギー保存則と指数的可積分性の両方を継承する。
• 数値解法の強収束率は $ \mathbb{E}[\|u^* - u^{N,k+1}\|_{L^2}^2] \leq C \max(h^2, \epsilon^8) \left( (1/2)^{k+1} + \exp(-\eta / \epsilon^2) \right) $ として確立され、時間方向で2次収束であり、誤差バウンドに高次指数的減衰が見られる。
• 数値実験により理論的収束率が確認され、提案手法の精度と安定性が裏付けられる。
• 本手法はエネルギー関数のモーメント有界性を保証し、ガウス尾部推定により大きな偏差の確率を制御し、スキームの頑健性を支持する。