[論文レビュー] Freezing Stochastic Travelling Waves

Gabriel J. Lord, Vera Thümmler|arXiv (Cornell University)|Jun 2, 2010

Stochastic processes and financial applications参考文献 37被引用数 4

ひとこと要約

本稿では、乗法的および加法的ノイズを伴うナグモ方程式における確率的移動波解および波速の計算のための新しい数値的手法を提案する。波を固定化する目的で、レベル集合追跡法とL²ノルム最小化法を比較する。波形の安定化は、基準関数からのずれを最小化することで実現され、波速が時間に依存する確率的変数として効率的に計算可能となる。主な発見として、ストラトノヴィッチノイズは波速と波の支持集合を増加させ、イトノイズは波速を減少させ、相関長にあまり依存しないことが示された。

ABSTRACT

We consider in this paper travelling wave solutions to stochastic partial differential equations and corresponding wave speed. As a particular example we consider the Nagumo equation with multiplicative noise which we mainly consider in the Stratonovich sense. A standard approach to determine the position and hence speed of a wave is to compute the evolution of a level set. We compare this approach against an alternative where the wave position is found by minimizing the $L^2$ norm against a fixed profile. This approach can also be used to stop (or freeze) the wave and obtain a stochastic partial differential algebraic equation that we then discretize and solve. Although attractive as it leads to a smaller domain size it can be numerically unstable due to large convection terms. We compare numerically the different approaches for estimating the wave speed. Minimization against a fixed profile works well provided the support of the reference function is not too narrow. We then use these techniques to investigate the effect of both \Ito and Stratonovich noise on the Nagumo equation as correlation length and noise intensity increases.

研究の動機と目的

• 確率的偏微分方程式における波速の計算のための数値的手法を開発・比較すること、特にノイズが存在する際の移動波に焦点を当てる。
• イトおよびストラトノヴィッチノイズがナグモ方程式における波速、幅、支持集合に与える影響を調査すること。
• 固定されたプロファイルに対するL²ノルム最小化に基づく波固定化手法を導入・検証し、計算領域を小さくし、シミュレーションを効率化可能とする。
• ノイズの強度および空間的相関長が、確率的系における波動の伝播および安定性に与える影響を分析すること。

提案手法

• 標準的手法では、波の位置をレベル集合により追跡し、固定されたuレベル（例：u=0.5）の位置を用いて波速を計算する。
• 代替的手法では、時間発展する解と固定された基準プロファイルとの間のL²ノルムを最小化することで、波を計算領域内で「固定」する。
• 固定化技術は、波動の運動を相殺する対流項を導入する。この項はL²最小化条件から得られ、確率的偏微分方程式と代数方程式の組み合わせ（SPDAE）を形成する。
• SPDAEは離散化され、数値的に解かれる。波速は最小化条件から時間に依存する変数として計算される。
• ノイズはイトおよびストラトノヴィッチ形式で実装され、乗法的ノイズはuに依存し、加法的ノイズは定数である。
• 数値的シミュレーションにより、両手法からの波速推定値を比較し、平均波速、分布、波面ダイナミクスなどの統計的性質を分析する。

実験結果

リサーチクエスチョン

• RQ1 レベル集合法とL²最小化法は、確率的移動波の波速推定において、どのように比較されるか？
• RQ2 ノイズ強度および空間的相関長は、ナグモ方程式における波速および波形幅にどのような影響を与えるか？
• RQ3 イトとストラトノヴィッチノイズの選択は、波動の伝播および安定性にどのように影響するか？
• RQ4 L²最小化に基づく波固定化手法は、強い対流項が存在する場合でも安定的かつ正確な結果をもたらすか？
• RQ5 加法的ノイズ下での波速および波面ダイナミクスの統計的挙動は何か、特に核生成および消滅の観点からどうか？

主な発見

• ストラトノヴィッチノイズは波速と波の空間的サポートを両方とも増加させ、ノイズ強度が高くなるほど波面が広がる。
• ストラトノヴィッチノイズの空間的相関長を増加させると波速が低下し、波のサポートも減少する。これは波形の形状に対する安定化効果を示唆する。
• イトノイズは波速を低下させ、相関長にほとんど感受性がない。これはストラトノヴィッチノイズとは異なるメカニズムによる波の抑制を示唆する。
• 加法的ノイズはノイズ強度が高くなるにつれて波速を増加させ、乗法的ノイズと同様の傾向を示す。これは波面前方の微小な摂動が伝播を促進するためと推察される。
• L²最小化手法は、基準関数のサポートが波の実際のサポートと著しく狭くない限り信頼性があるが、それ以外の場合は数値的失敗が生じる。
• 波固定化SPDAEアプローチにより、計算領域のサイズとノイズパス生成コストが削減されるが、強い対流項による数値的不安定性は依然として課題のままである。