[论文解读] Energy-preserving exponential integrable numerical method for stochastic cubic wave equation with additive noise

Jianbo Cui, Jialin Hong+2 more|arXiv (Cornell University)|Sep 2, 2019

Stochastic processes and financial applications参考文献 34被引用 5

一句话总结

本文提出了一种新颖的能量保持、指数可积的数值方法，用于带加性噪声的随机三次波方程。通过结合谱伽辽金空间半离散化与时间方向的分裂技巧及平均向量场方法，该格式继承了精确解的能量演化规律和指数可积性，实现了强收敛率 $ O(h^2 + \epsilon^8) $，其中 $ h $ 为时间步长，$ \epsilon $ 控制截断误差。

ABSTRACT

In this paper, we present an energy-preserving exponentially integrable numerical method for stochastic wave equation with cubic nonlinearity and additive noise. We first apply the spectral Galerkin method to discretizing the original equation and show that this spatial discretization possesses an energy evolution law and certain exponential integrability property. Then the exponential integrability property of the exact solution is deduced by proving the strong convergence of the semi-discretization. To propose a full discrete numerical method which could inherit both the energy evolution law and the exponential integrability, we use the splitting technique and averaged vector field method in the temporal direction. Combining these structure-preserving properties with regularity estimates of the exact and the numerical solutions, we obtain the strong convergence rate of the numerical method. Numerical experiments coincide with these theoretical results.

研究动机与目标

• 开发一种全离散数值格式，以同时保持带加性噪声的随机三次波方程的能量演化规律与指数可积性。
• 在空间采用谱伽辽金半离散化、时间采用分裂技巧与平均向量场方法的假设下，建立该数值方法的强收敛性。
• 通过利用谱伽辽金逼近的统一指数可积性与正则性估计，证明精确解的指数可积性。
• 弥合带三次非线性的随机波方程在结构保持数值方法方面的空白，特别是保持其内在物理特性。
• 通过与理论预测一致的数值实验验证理论收敛速率。

提出的方法

• 空间半离散化采用谱伽辽金方法，该方法保持能量演化规律，并在半离散系统中表现出指数可积性。
• 时间半离散化采用分裂技巧与平均向量场方法相结合，以确保能量保持并维持时间方向的指数可积性。
• 该方法被设计为继承精确解的能量演化规律与指数可积性，从而确保长期稳定性和精度。
• 推导了精确解与数值解的先验估计与正则性界，以支持收敛性分析。
• 利用高斯尾部估计与索伯列夫嵌入定理，控制数值解中大偏差的概率。
• 应用迭代误差估计推导最终的强收敛速率，误差满足 $ \mathbb{E}[\|u^* - u^{N,k+1}\|_{L^2}^2] \leq C \max(h^2, \epsilon^8) \left( (1/2)^{k+1} + \exp(-\eta / \epsilon^2) \right) $。

实验结果

研究问题

• RQ1 能否构造一种全离散数值格式，使其同时保持带加性噪声的随机三次波方程的能量演化规律与指数可积性？
• RQ2 在适当的正则性与可积性假设下，此类结构保持格式的强收敛速率是多少？
• RQ3 随机三次波方程的精确解是否具有指数可积性，且该性质能否被数值方法继承？
• RQ4 谱伽辽金方法与分裂技巧如何协同作用，以保持方程的内在几何与概率结构？
• RQ5 所提方法能否在保持长期稳定性与能量守恒的同时实现最优收敛速率？

主要发现

• 带加性噪声的随机三次波方程的精确解满足指数可积性：$ \mathbb{E}\left[\exp\left(\int_0^T c\|u(s)\|_{L^6}^2 ds\right)\right] \leq C $，其中 $ C $ 依赖于初始数据、噪声、时间与维度。
• 谱伽辽金半离散化保持了能量演化规律与指数可积性，从而支持全离散格式的强收敛性推导。
• 所提出的全离散方法，结合空间的谱伽辽金与时间的分裂-平均向量场方法，继承了能量演化与指数可积性双重特性。
• 该数值方法的强收敛速率被确立为 $ \mathbb{E}[\|u^* - u^{N,k+1}\|_{L^2}^2] \leq C \max(h^2, \epsilon^8) \left( (1/2)^{k+1} + \exp(-\eta / \epsilon^2) \right) $，表明时间方向为二阶收敛，且误差界呈现高阶指数衰减。
• 数值实验验证了理论收敛速率，证实了该方法的精度与稳定性。
• 该方法确保了能量泛函的矩有界性，并通过高斯尾部估计控制大偏差，支持了该格式的鲁棒性。