[论文解读] A MODIFIED SEMI{IMPLICT EULER-MARUYAMA SCHEME FOR FINITE ELEMENT DISCRETIZATION OF SPDES

Gabriel J. Lord, Antoine Tambue|arXiv (Cornell University)|Apr 12, 2010

Stochastic processes and financial applications参考文献 26被引用 22

一句话总结

本文提出了一种改进的半隐式 Euler-Maruyama 格式，结合有限元空间半离散化，用于求解带有空间-时间噪声的二阶半线性抛物型随机偏微分方程（SPDEs）。与标准半隐式 Euler-Maruyama 格式相比，该方法在均方根 L² 范数和 H¹ 范数下的收敛率均有提升，尤其在具有空间粗糙噪声的扩散方程和对流-扩散-反应方程中表现更优。

ABSTRACT

We consider the numerical approximation of a general second order semi{linear parabolic stochastic partial dierential equation (SPDE) driven by additive space-time noise. We introduce a new scheme using in time a linear functional of the noise with a semi{implicit Euler{ Maruyama method and in space we analyse a nite element method although extension to nite dierences or nite volumes would be possible. We consider noise that is white in time and either in H 1 or H 2 in space. We give the convergence proofs in the root mean square L 2 norm for a diusion reaction equation and in root mean square H1 norm in the presence of advection. We examine the regularity of the initial data, the regularity of the noise and errors from projecting the noise. We present numerical results for a linear reaction diusion equation in two dimensions as well as a nonlinear example of two-dimensional stochastic advection diusion reaction equation. We see from both the analysis and numerics that we have better convergence properties over the standard semi{implicit Euler{Maruyama method.

研究动机与目标

• 开发一种更精确的数值格式，用于求解由加法型空间-时间噪声驱动的半线性抛物型 SPDE。
• 在存在空间粗糙噪声的条件下，改进现有半隐式 Euler-Maruyama 方法的收敛性质。
• 分析噪声正则性（空间上为 H¹ 或 H²）以及初值正则性对数值收敛性的影响。
• 研究数值格式中由噪声离散化引起的投影误差。
• 通过二维线性和非线性 SPDE 的数值实验验证理论结果。

提出的方法

• 该格式采用半隐式 Euler-Maruyama 时间离散化，并引入空间-时间噪声的线性泛函，以增强稳定性和精度。
• 空间离散化采用有限元方法，尽管理论上可扩展至有限差分或有限体积方法。
• 该方法在反应-扩散方程上分析了均方根 L² 范数下的收敛性，在对流-扩散-反应方程上分析了均方根 H¹ 范数下的收敛性。
• 噪声在时间上建模为白噪声，在空间上为 H¹ 或 H² 空间，对噪声投影误差进行了仔细处理。
• 理论收敛性证明综合考虑了初值和噪声的正则性，确保了稳健的误差界。
• 在二维线性反应-扩散方程和二维非线性随机对流-扩散-反应方程上进行了数值实验，以验证该方法。

实验结果

研究问题

• RQ1 所提出的改进半隐式 Euler-Maruyama 格式相较于标准方法，在带有加法型空间-时间噪声的 SPDE 中，收敛率有何提升？
• RQ2 空间噪声正则性（H¹ 或 H²）对数值格式收敛行为有何影响？
• RQ3 初值正则性和噪声投影对格式整体误差有何影响？
• RQ4 在存在对流项的情况下，该方法是否能在 H¹ 范数下保持更优的收敛性？
• RQ5 数值结果在多大程度上验证了线性和非线性 SPDE 测试案例中的理论收敛率？

主要发现

• 所提出的格式在均方根 L² 范数和 H¹ 范数下均优于标准半隐式 Euler-Maruyama 方法的收敛率。
• 对于二维线性反应-扩散方程，在相同噪声和初值条件下，该方法表现出更优的收敛行为。
• 在对流-扩散-反应情形下，该方法在 H¹ 范数下仍保持更优的收敛性，与理论预期一致。
• 分析表明，噪声在 H¹ 或 H² 空间中的正则性显著影响收敛率，正则性越高，误差控制越佳。
• 数值结果证实，噪声离散化引起的投影误差在该格式中得到有效控制，有助于实现稳定且精确的模拟。
• 该方法在线性和非线性 SPDE 上均表现出鲁棒性，在所有测试场景中均一致优于标准方法。