Uniform-in-time Strong Error Estimates of Tamed-FEM to Superlinear SPDEs driven by Multiplicative Noise

技術要約：マルチプリカティブ・ノイズ（乗法的ノイズ）に駆動される超線形SPDEに対する、Tamed-FEMの一様時間強誤差評価

問題設定
本論文は、マルチプリカティブ・ノイズに駆動される、中程度の厚みの界面を持つ確率アレン・カーン方程式などの、ある種の超線形確率偏微分方程式（SPDE）の数値近似を取り扱う。対象となる支配方程式は、超線形成長を示すドリフト項と、リプシッツ連続な拡散項を持つ半線形放物型SPDEである。有限の時間地平におけるこれらの誤差評価については既に確立されているが、無限の時間地平における評価の分析は依然として困難である。標準的な陽解法（オイラー・マルーヤマ法など）は、超線形問題においてモーメントの爆発（blow-up）を起こす可能性がある一方、陰解法は安定ではあるものの、高い計算コストを要する。主な目的は、長時間の無条件安定性を保証し、SPDEのリアプノフ構造を保持する、非線形明示的スキームを開発することである。

手法
著者らは、SPDEの空間離散化に適用される、[19]で導入された**非線形明示的タムド有限要素法（tamed-FEM）**を分析している。この手法は、空間離散化のためのガラーキン有限要素法（FEM）と、時間離散化のためのタムド・オイラー・スキームを組み合わせたものである。「タミング（taming）」メカニズムは、超線形ドリフト $f(\xi)$ を $f_\tau(\xi) = f(\xi)(1 + \tau|\xi|^{2q})^{-1/2}$ （ここで $\tau$ はタイムステップ）へと修正することで、明示的な計算可能性を維持しながら数値的な爆発を防ぐ。

理論的分析は、以下の主要なステップを通じて進行する：

1. 一様時間安定性： 著者らは、$\dot{H}^1$ ノルムにおけるタムド-FEMの厳密な一様時間無条件安定性評価を確立する。これは、タムド関数の新しい評価、具体的には $f'_\tau + \tau|f'_\tau|^2$ の一様上界の導出に依拠しており、これは無限の時間間隔にわたってモーメントの成長を制御するために極めて重要である。
2. 補助プロセス分析： 真の解と数値スキームの間のギャップを埋めるために、補助プロセスを構築する。著者らは、解析的な $C_0$ セミグループの平滑化特性と真の解の正則性を利用して、真の解とこの補助プロセスとの間の強誤差評価を導出する。
3. 誤差分解： 全強誤差は、真の解と補助プロセスとの間の誤差、および補助プロセスとタムド-FEM解との間の誤差に分解される。変分法と伊藤（Itô）解析が、これらの項を一様に抑え込むために用いられる。
4. エルゴード性分析： リアプノフ構造と初期データに対する指数的連続性の評価を活用することで、著者らはタムド-FEMの指数エルゴード性を証明する。これにより、真の不変測度と数値的不変測度の間のワッサースタイン（Wasserstein）-2距離における誤差評価の導出が可能となる。

主な貢献と結果

• 一様時間強誤差評価： 本論文は、マルチプリカティブ・ノイズに駆動される超線形SPDEに適用される非線形明示的スキームに対する、初の一様時間強誤差評価を確立した。
  • マルチプリカティブ・ノイズの場合、強収束次数は $O(h^{1+\gamma} + \tau^{1/2})$ と示される。ここで $h$ は空間メッシュサイズであり、$\gamma \in [0, 1]$ は初期データの正則性に関連する。
  • **アディティブ・ノイズ（加法的ノイズ）**の場合、より高次の収束次数 $O(h^{1+\gamma} + \tau^{(1+\gamma)/2})$ が導かれる。
  • 極めて重要な点として、これらの誤差評価における定数は時間地平に依存せず、$t \to \infty$ においても評価が有効であることを保証している。
• 指数エルゴード性と不変測度の収束： 著者らは、タムド-FEMが指数エルゴード的であることを証明した。その結果、真の不変測度と数値的不変測度の間のエルゴード誤差を、マルチプリカティブ・ノイズの場合は $O(h^{2-\epsilon} + \tau^{1/2})$、アディティブ・ノイズの場合は $O(h^{2-\epsilon} + \tau^{1-\epsilon})$ のワッサースタイン-2距離において導出した。
• 数値検証： 二次元領域における確率アレン・カーン方程式を用いた数値実験により、理論的知見を確認した。テストは以下の事項を実証している：
  • 時間的および空間的離散化の両方に対する、予測された収束率の鋭さ。
  • 誤差定数の時間独立性。すなわち、シミュレーション時間が増加しても、誤差界は劣化しない。
  • スキームの指数エルゴード性。これは、異なる初期状態からのアンサンブル平均の急速な合流によって示されている。

意義
本論文は、マルチプリカティブ・ノイズに駆動される超線形SPDEの一意なエルゴード性を保持する非線形明示的スキームに対して、初の一様時間強誤差評価を提供すると主張している。計算コストの高い陰解法に頼ることなく、明示的スキームにおけるモーメント爆発の課題を克服することで、本研究は相転移モデルやその他の超線形SPDEの長時間シミュレーションのための、堅牢かつ効率的な数値ツールを提供する。本成果は、これらの複雑な確率系における数値近似の長期挙動に関する文献の空白を埋めるものであり、数値的不変測度が真の不変測度へと定量可能な誤差範囲内で収束することを保証するものである。