[論文レビュー] Numerical approximation of the stochastic heat equation with a distributional reaction term

Ludovic Goudenège, El Mehdi Haress+1 more|arXiv (Cornell University)|May 13, 2024

Advanced Thermodynamics and Statistical Mechanics被引用数 1

ひとこと要約

本稿は、1次元空間における分布的反応項を伴う確率的熱方程式に、tamed Euler有限差分スキームを適用した際の収束速度を確立する。確率的ぬいぐるみ技術と駆動項のBesov正則性を活用することで、L^m収束を証明し、特異な駆動項の正則性に応じた収束速度を示す。駆動項が有界可測関数の場合、空間でほぼ最適な(1/2 − ε)、時間で(1/4 − ε)の収束速度を達成する。また、駆動項が有限測度である場合でさえも正の収束速度が得られる。

ABSTRACT

We study the numerical approximation of the stochastic heat equation with a distributional reaction term. Under a condition on the Besov regularity of the reaction term, it was proven recently that a strong solution exists and is unique in the pathwise sense, in a class of H\"older continuous processes. For a suitable choice of sequence $(b^k)_{k\in \mathbb{N}}$ approximating $b$, we prove that the error between the solution $u$ of the SPDE with reaction term $b$ and its tamed Euler finite-difference scheme with mollified drift $b^k$, converges to $0$ in $L^m(\Omega)$ with a rate that depends on the Besov regularity of $b$. In particular, one can consider two interesting cases: first, even when $b$ is only a (finite) measure, a rate of convergence is obtained. On the other hand, when $b$ is a bounded measurable function, the (almost) optimal rate of convergence $(\frac{1}{2}-\varepsilon)$-in space and $(\frac{1}{4}-\varepsilon)$-in time is achieved. Stochastic sewing techniques are used in the proofs, in particular to deduce new regularising properties of the discrete Ornstein-Uhlenbeck process.

研究の動機と目的

• 有界可測駆動項を超える、Besov空間に属する分布的駆動項を伴う確率的熱方程式の数値収束結果を拡張すること。
• 経路ごとの意味で、滑らか化された駆動項を用いたtamed Euler有限差分スキームの強収束速度を確立すること。
• 反応項bのBesov正則性に基づいて、L^m(Ω)における収束速度を定量化すること。
• bが関数でない有限測度である場合でも収束が達成可能であることを示すこと。

提案手法

• 空間方向に中央差分、時間方向に明示的刻み幅を用いたtamed Euler有限差分スキームを構築し、滑らか化された駆動項列(bk)k∈Nを導入する。
• 特に離散オーナイズ・ウーレンプソン過程とその正則化性を扱うために適応された確率的ぬいぐるみ技術を用いて解析を行う。
• 真の解uと数値解un,kとの間の誤差を、モーメント推定と駆動項差分‖b − bk‖B^{γ−1}_pのBesovノルム制御を用いて評価する。
• 弱い正則性仮定のもとで数値誤差の増大を制御するため、対数因子を含むグローバル型の補題を用いる。
• 特異な駆動項の存在下でも誤差制御を可能にするために、臨界指数を備えた修正版の確率的ぬいぐるみ補題に依拠する。
• モーメントバウンド、離散半群推定、および解の経路ごとの正則性を組み合わせることで収束を確立する。

実験結果

リサーチクエスチョン

• RQ1 反応項bがBesov空間B^γ_{p,∞}に属する分布である場合、確率的熱方程式の数値解の収束速度はどの程度達成可能か？
• RQ2 bが関数ではなく有限測度である場合、tamed Euler有限差分スキームはL^m(Ω)で収束するか？
• RQ3 bが有界可測関数である場合、時間と空間における最適な収束速度は何か？ また、既知の鋭い収束速度（時間で1/4、空間で1/2）にどの程度近づけるか？
• RQ4 確率的ぬいぐるみ技術と離散オーナイズ・ウーレンプソン過程の正則化性は、特異な駆動項の下で誤差推定にどのように寄与するか？
• RQ5 対数項を含む臨界グローバル補題を用いることで、低正則性駆動項を伴う数値スキームにおける誤差伝搬を制御できるか？

主な発見

• 数値スキームは、駆動項bのBesov正則性、特に‖b − bk‖B^{γ−1}_pに依存したL^m(Ω)収束を示す。
• bが有限測度（すなわちγ − 1/p ≤ −1）である場合でさえも、正の収束速度が得られる。
• 有界可測駆動項（γ = 0, p = ∞）の場合、空間で(1/2 − ε)、時間で(1/4 − ε)の（ほぼ）最適収束速度が達成される。
• 弱い正則性下でも離散解誤差を制御可能な、臨界指数を備えた確率的ぬいぐるみ補題の新規応用に依拠する。
• 離散オーナイズ・ウーレンプソン過程に新たな正則化性が導出され、数値スキームの誤差バウンドにおいて中心的な役割を果たす。
• 対数項を含む臨界グローバル補題の使用により、駆動項の非リプシッツ性を扱い、スキーム内の誤差増大を制御できる。