[论文解读] A modified semi--implict Euler-Maruyama Scheme for finite element discretization of SPDEs with additive noise

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

Structural Health Monitoring Techniques被引用 4

一句话总结

本文提出了一种改进的半隐式 Euler-Maruyama 格式，结合有限元空间半径化，用于求解带有加法型时空噪声的二阶随机偏微分方程。该方法通过引入噪声的线性泛函，相较于标准半隐式 Euler-Maruyama 格式，在均方根 $L^2$ 范数下提升了收敛速率，通过理论分析与二维线性和非线性反应-扩散-对流方程的数值实验得到验证。

ABSTRACT

We consider the numerical approximation of a general second order semi--linear parabolic stochastic partial differential equation (SPDE) driven by additive space-time noise. We introduce a new modified scheme using a linear functional of the noise with a semi--implicit Euler--Maruyama method in time and in space we analyse a finite element method (although extension to finite differences or finite volumes would be possible). We prove convergence in the root mean square $L^{2}$ norm for a diffusion reaction equation and diffusion advection reaction equation. We present numerical results for a linear reaction diffusion equation in two dimensions as well as a nonlinear example of two-dimensional stochastic advection diffusion reaction equation. We see from both the analysis and numerics that the proposed scheme has better convergence properties than the standard semi--implicit Euler--Maruyama method.

研究动机与目标

• 开发一种更精确的数值格式，用于强逼近带有加法噪声的二阶半线性抛物型 SPDE。
• 克服标准半隐式 Euler-Maruyama 格式在存在非线性 Nemytskii 算子时的阶次降低限制。
• 实现通过有限元方法进行空间半径化，而无需显式了解线性算子和噪声协方差算子的特征函数与特征值。
• 在扩散-反应方程和对流-扩散-反应方程中，建立均方根 $L^2$ 范数下的收敛性。
• 通过理论分析与二维数值模拟，展示改进的收敛特性。

提出的方法

• 构建一种基于噪声线性泛函的改进半隐式 Euler-Maruyama 时间步进格式，该泛函源自每个傅里叶模态下随机卷积的精确方差。
• 采用有限元方法进行空间半径化，可在无需谱基知识的前提下灵活处理复杂区域与边界条件。
• 该格式使用随机卷积近似 $O_m = \int_{t_m}^{t_{m+1}} e^{(t_{m+1}-s)A} dW(s)$，其中噪声项通过 $\mathbf{E}[(e_i, O_m)^2] = e^{-2\lambda_i \Delta t} \frac{q_i}{2\lambda_i} (1 - e^{-2\lambda_i \Delta t})$ 计算。
• 通过采用改进的时间步进格式，避免了高阶格式中严格的交换性假设，同时保持了稳定性和精度。
• 利用离散 Gronwall 引理及对有限元投影与时间积分误差项的估计，证明了理论收敛性。
• 数值实现采用具有混合边界条件的矩形区域，噪声采用 Matérn 型协方差，通过无限域积分与精确高斯计算进行近似。

实验结果

研究问题

• RQ1 改进的半隐式 Euler-Maruyama 格式能否在带有加法噪声的 SPDE 中实现优于标准格式的收敛速率？
• RQ2 在无需谱基知识的前提下，使用噪声的线性泛函是否能提升有限元半径化中的稳定性和精度？
• RQ3 在一般边界条件下，该方法能否对非线性对流-扩散-反应 SPDE 保持高阶收敛性？
• RQ4 在均方根 $L^2$ 范数下，改进格式的收敛速率与标准半隐式 Euler-Maruyama 方法相比如何？
• RQ5 当与噪声线性化结合时，空间有限元半径化对整体收敛性有何影响？

主要发现

• 所提格式在均方根 $L^2$ 范数下，相较于标准半隐式 Euler-Maruyama 方法，对线性和非线性 SPDE 均实现了更优的收敛速率。
• 对于扩散-反应方程，该格式表现出收敛阶 $\mathcal{O}(h^\beta + \Delta t^{\beta/2})$（$\beta < 2$），其中 $\beta$ 与解的正则性相关。
• 二维线性反应-扩散方程的数值结果验证了理论收敛速率，并显示误差显著低于标准格式。
• 非线性二维随机对流-扩散-反应方程的测试案例表明，该格式在复杂动力学下仍具有鲁棒性与稳定性。
• 该方法成功处理了具有混合边界条件且噪声协方差不可分的 SPDE，而无需显式进行特征函数分解。
• 理论分析在较弱假设下确认了收敛性，包括 $X(t)$ 为 $H^1(\Omega)$-值过程，且 $F$ 满足 Lipschitz 条件与增长条件。