[论文解读] Full discretisation of semi-linear stochastic wave equations driven by multiplicative noise

Rikard Anton, David Cohen+2 more|arXiv (Cornell University)|Feb 28, 2015

Stochastic processes and financial applications参考文献 17被引用 10

一句话总结

本文提出了一种用于带乘性噪声的半线性随机波动方程的全离散数值格式，结合线性有限元空间离散化与随机三角方法（STM）进行时间积分。主要贡献在于均方收敛性，其误差界与空间网格尺寸无关，同时保持了迹公式（期望能量的线性漂移），优于传统格式如欧拉-马鲁亚马方法在长时间能量行为上的表现。

ABSTRACT

A fully discrete approximation of the semi-linear stochastic wave equation driven by multiplicative noise is presented. A standard linear finite element approximation is used in space and a stochastic trigonometric method for the temporal approximation. This explicit time integrator allows for mean-square error bounds independent of the space discretisation and thus do not suffer from a step size restriction as in the often used Störmer-Verlet-leap-frog scheme. Furthermore, it satisfies an almost trace formula (i.e., a linear drift of the expected value of the energy of the problem). Numerical experiments are presented and confirm the theoretical results.

研究动机与目标

• 为由乘性噪声驱动的半线性随机波动方程开发一种全离散数值格式。
• 在均方意义下实现强收敛，且误差界与空间离散化无关。
• 在数值解中保持几何性质，如迹公式（期望能量的线性漂移）。
• 克服传统格式如斯特默-维尔特-跃迁法中固有的时间步长限制。

提出的方法

• 在具有多边形边界的有界凸域上，使用标准线性有限元方法进行空间离散化。
• 时间积分采用随机三角方法（STM），一种专为随机波动方程设计的显式时间积分格式。
• 将STM应用于空间半离散化后得到的半离散系统，确保稳定性与精度。
• 通过希尔伯特空间中的能量估计与插值技术进行分析，特别关注加权索伯列夫范数。
• 为精确解、有限元半离散解与全离散解分别推导迹公式，表明期望能量呈现线性漂移。
• 数值实验验证了理论收敛速率以及在不同时间步长与网格尺寸下迹公式的保持性。

实验结果

研究问题

• RQ1 能否为带乘性噪声的半线性随机波动方程构造一种全离散格式，使其在均方收敛性上与空间网格尺寸无关？
• RQ2 所提出的随机三角方法是否能保持精确解的迹公式（期望能量的线性漂移）？
• RQ3 在长时间能量演化模拟中，STM与经典格式如欧拉-马鲁亚马法和克兰克-尼科尔森-马鲁亚马法相比表现如何？
• RQ4 全离散格式在均方意义下的收敛速率是多少？其对空间与时间离散化参数的依赖关系如何？
• RQ5 该方法在大时间步长下是否仍能保持如能量漂移行为等几何性质？

主要发现

• 所提出的全离散格式实现了均方收敛，且误差界与空间网格尺寸无关，避免了斯特默-维尔特-跃迁法中常见的时间步长限制。
• 随机三角方法保持了迹公式，确保数值解的期望能量呈现与精确解一致的线性漂移。
• 数值实验确认了理论收敛速率，观测到的斜率与时间与空间方向上预测的收敛阶数一致。
• 在长时间模拟中，STM优于经典欧拉-马鲁亚马类格式，能正确捕捉期望能量的线性漂移，而欧拉格式则表现出错误的能量增长。
• 即使在大时间步长（如k=0.1）与长时域（如[0,250]）下，STM仍能保持期望能量演化的优异行为，证实其鲁棒性。
• 数值实验中的蒙特卡洛误差可忽略不计（如<0.01），验证了所观测收敛性与迹公式行为的可靠性。