Uniform-in-time Strong Error Estimates of Tamed-FEM to Superlinear SPDEs driven by Multiplicative Noise

技术摘要：针对受乘性噪声驱动的超线性 SPDE 的一致时间强误差估计（Tamed-FEM）

问题陈述
本文研究了一类由乘性噪声驱动的超线性随机偏微分方程（SPDE）的数值近似问题，具体为具有中等厚度界面的随机 Allen–Cahn 方程。该控制方程是一个半线性抛物型 SPDE，其漂移项表现出超线性增长，而扩散项是 Lipschitz 连续的。虽然这类方程在有限时间范围内的强误差估计已有成熟结论，但在无限时间范围内的此类估计分析仍具挑战性。标准显式格式（如 Euler–Maruyama 方法）在处理超线性问题时往往会出现矩爆炸（moment blow-up）现象，而隐式格式虽然稳定，却会带来高昂的计算成本。主要目标是开发一种非线性显式格式，以确保长时期的无条件稳定性，保持 SPDE 的 Lyapunov 结构，并提供尖锐的一致时间强误差估计。

方法论
作者分析了一种此前在 [19] 中引入的非线性显式被驯化有限元法（tamed-FEM），该方法应用于 SPDE 的空间离散。该方法结合了用于空间离散的 Galerkin 有限元法（FEM）和用于时间离散的被驯化 Euler 格式。其“被驯化”（taming）机制通过修改超线性漂移项 $f(\xi)$ 为 $f_\tau(\xi) = f(\xi)(1 + \tau|\xi|^{2q})^{-1/2}$（其中 $\tau$ 为时间步长），以防止数值爆炸，同时保持显式计算的可行性。

理论分析通过以下几个关键步骤进行：

1. 一致时间稳定性： 作者建立了被驯化-FEM 在 $\dot{H}^1$ 范数下的严格一致时间无条件稳定性估计。这依赖于对被驯化函数推导新的估计，特别是对于 $f'_\tau + \tau|f'_\tau|^2$ 的一致上界，这对于控制无限时间间隔内矩的增长至关重要。
2. 辅助过程分析： 构建了一个辅助过程来弥合精确解与数值格式之间的差距。作者利用解析 $C_0$-半群的平滑性质和精确解的正规性，推导了精确解与该辅助过程之间的强误差估计。
3. 误差分解： 总强误差被分解为精确解与辅助过程之间的误差，以及辅助过程与被驯化-FEM 解之间的误差。通过变分法和 Itô 微积分来界定这些项在时间上的一致性。
4. 遍历性分析： 利用 Lyapunov 结构和关于初始数据的指数连续性估计，作者证明了被驯化-FEM 的指数遍历性。这使得能够推导出精确不变测度与数值不变测度在 Wasserstein-2 距离下的误差估计。

主要贡献与结果

• 一致时间强误差估计： 本文为应用于受乘性噪声驱动的超线性 SPDE 的非线性显式格式建立了首个一致时间强误差估计。
  • 对于乘性噪声情况，强收敛阶被证明为 $O(h^{1+\gamma} + \tau^{1/2})$，其中 $h$ 是空间网格尺寸，$\gamma \in [0, 1]$ 与初始数据的正则性相关。
  • 对于加性噪声情况，推导出了更高阶的收敛阶 $O(h^{1+\gamma} + \tau^{(1+\gamma)/2})$。
  • 至关重要的是，这些误差界中的常数与时间范围无关，确保了估计在 $t \to \infty$ 时依然有效。
• 指数遍历性与不变测度收敛： 作者证明了被驯化-FEM 具有指数遍历性。因此，他们推导了在 Wasserstein-2 距离下精确不变测度与数值不变测度之间的遍历误差估计，显示出在乘性噪声下收敛阶为 $O(h^{2-\epsilon} + \tau^{1/2})$，在加性噪声下为 $O(h^{2-\epsilon} + \tau^{1-\epsilon})$。
• 数值验证： 在二维区域上的随机 Allen–Cahn 方程的数值实验证实了上述理论发现。测试表明：
  • 预测的关于时间与空间离散的收敛阶是尖锐的。
  • 误差常数具有时间无关性，即随着模拟时间从 $T=1$ 增加到 $T=50$，误差界并未恶化。
  • 该格式具有指数遍历性，这通过不同初始状态下的系综均值（ensemble means）快速合并得到了证实。

意义
本文声称提供了首个针对受乘性噪声驱动的超线性 SPDE 的非线性显式格式的一致时间强误差估计，且该格式保留了唯一的遍历性。通过在不诉诸高计算成本的隐式方法的情况下克服显式格式中矩爆炸的挑战，这项工作为相变模型及其他超线性 SPDE 的长时模拟提供了一种稳健且高效的数值工具。研究结果填补了关于此类复杂随机系统数值近似长期行为的文献空白，确保了数值不变测度能以可量化的误差界向真实不变测度收敛。