[论文解读] Numerical solutions of some hyperbolic stochastic partial differential equations with mixed derivatives including sine-Gordon equation

Henry C. Tuckwell|arXiv (Cornell University)|Aug 2, 2015

Stochastic processes and financial applications参考文献 28被引用 3

一句话总结

本文提出了一种显式数值格式，用于求解具有混合偏导数的双曲型随机偏微分方程（SPDEs），包括由时空高斯白噪声驱动的Sine-Gordon方程。该方法在确定性和随机两种情形下均能准确捕捉类似波的结构，即使在零初始条件下，仅由噪声驱动时，解仍表现出鲁棒的收敛性与相干动力学的出现，尤其在小噪声强度（σ ≤ 0.25）时更为显著。

ABSTRACT

We consider linear and nonlinear hyperbolic SPDEs with mixed derivatives with additive space-time Gaussian white noise of the form $Y_{xt}=F(Y) + σW_{xt}.$ Such equations, which transform to linear and nonlinear wave equations, including Klein-Gordon, Liouville's and the sine-Gordon equation, are related to what Zimmerman (1972) called a diffusion equation. An explicit numerical scheme is employed in both deterministic and stochastic examples. The scheme is checked for accuracy against known exact analytical solutions for deterministic equations. In the stochastic case with $F=0$, solutions yield sample paths for the Brownian sheet whose statistics match well exact values. Generally the boundary conditions are chosen to be initial values $Y(x,0)$ and boundary values $Y(0,t)$ on the quarter-plane or subsets thereof, which have been shown to lead to existence and uniqueness of solutions. For the linear case solutions are compared at various grid sizes and wave-like solutions were found, with and without noise, for non-zero initial and boundary conditions. Surprisingly, wave-like structures seemed to emerge with zero initial and boundary conditions and purely noise source terms with no signal. Equations considered with nonlinear $F$ included quadratic and cubic together with the sine-Gordon equation. For the latter, wave-like structures were apparent with $σ\le 0.25$ but they tended to be shattered at larger values of $σ$. Previous work on stochastic sine-Gordon equations is briefly reviewed.

研究动机与目标

• 开发并验证一种用于求解具有混合偏导数和加法型时空白噪声的双曲型SPDE的数值格式。
• 研究在初始条件和边界条件均为零但存在噪声时，随机SPDE中类似波行为的出现机制。
• 通过与已知解析解对比，评估数值方法在确定性情形下的精度与收敛性。
• 研究非线性SPDE（包括Sine-Gordon方程）在不同噪声强度（σ）下的行为特征。
• 通过坐标变换，探讨混合导数SPDE形式的解与经典波动方程之间的关系。

提出的方法

• 采用显式有限差分格式求解形如 $ Y_{xt} = F(Y) + \sigma W_{xt} $ 的SPDE，其中 $ W_{xt} $ 为时空白噪声。
• 通过线性与非线性确定性情形的精确解析解对格式进行验证，包括波动方程与Liouville方程。
• 在随机情形下，模拟布朗曲面的样本路径，并将其统计特性与理论预期进行比较。
• 在四分之一平面区域 $[0, T] \times [0, X]$ 上进行数值实验，初始与边界条件恒定（$ Y(x,0) $ 与 $ Y(0,t) $），以确保解的唯一性。
• 通过网格细化研究评估收敛性，比较不同网格尺寸（如 400×400 与 800×800）下的解，并改变试验次数。
• 将方法应用于非线性情形，包括二次、三次及Sine-Gordon型 $ F(Y) $，系统性地改变噪声强度 $ \sigma $。

实验结果

研究问题

• RQ1 显式数值格式能否准确模拟具有混合导数和加法型时空白噪声的双曲型SPDE？
• RQ2 当初始条件与边界条件均为零但系统仅由噪声驱动时，解中是否会出现类似波的结构？
• RQ3 噪声强度 $ \sigma $ 如何影响随机Sine-Gordon方程中类似波解的稳定性与相干性？
• RQ4 当 $ \sigma \to 0 $ 时，数值解在多大程度上收敛于确定性解？网格分辨率在该过程中起到何种作用？
• RQ5 混合导数SPDE $ Y_{xt} = F(Y) + \sigma W_{xt} $ 的解与经典波动方程 $ \phi_{tt} - \phi_{xx} = G(\phi) $ 之间存在何种关系？

主要发现

• 显式数值格式能准确再现确定性线性与非线性SPDE的已知解析解，验证了其可靠性。
• 对于 $ F = 0 $ 的线性SPDE，模拟的布朗曲面样本路径与理论统计特性一致，验证了噪声积分的正确性。
• 即使初始与边界条件为零，仅由噪声驱动时仍能出现类似波的结构，尤其在小 $ \sigma $ 时更为显著，表明存在内在的噪声驱动模式形成机制。
• 对于Sine-Gordon方程，当 $ \sigma \leq 0.25 $ 时，相干的类似波结构保持稳定，但在更高噪声强度下则出现碎片化。
• 网格细化（如 800×800 与 400×400 对比）使波结构更加清晰，标准差减小，表明随着分辨率提高，解呈现收敛性。
• 当 $ \sigma = 0.1 $ 时，随机模拟的均值与确定性解具有良好一致性，尤其在 800×800 网格与 50 次试验下，表明小噪声下具有良好的数值稳定性。