[论文解读] Uniform error bounds of a time-splitting spectral method for the long-time dynamics of the nonlinear Klein-Gordon equation with weak nonlinearity

Weizhu Bao, Yue Feng+1 more|arXiv (Cornell University)|Jan 29, 2020

Advanced Mathematical Physics Problems参考文献 57被引用 7

一句话总结

该论文为应用于弱非线性非线性克莱因-戈尔登方程的时间分裂谱方法建立了统一的误差界，证明了即使在大时间步长下，当非线性强度为 $\varepsilon^p$ 且模拟时间达到 $O(\varepsilon^{-\beta})$ 时，时间方向仍保持二阶精度。该方法在空间方向实现谱精度，并在所有 $\varepsilon \in (0,1]$ 范围内保持统一收敛性，数值验证了实数与复数振荡情形下的有效性。

ABSTRACT

We establish uniform error bounds of time-splitting Fourier pseudospectral (TSFP) methods for the nonlinear Klein--Gordon equation (NKGE) with weak power-type nonlinearity and $O(1)$ initial data, while the nonlinearity strength is characterized by $\varepsilon^{p}$ with a constant $p \in \mathbb{N}^+$ and a dimensionless parameter $\varepsilon \in (0, 1]$, for the long-time dynamics up to the time at $O(\varepsilon^{-β})$ with $0 \leq β\leq p$. In fact, when $0 < \varepsilon \ll 1$, the problem is equivalent to the long-time dynamics of NKGE with small initial data and $O(1)$ nonlinearity strength, while the amplitude of the initial data (and the solution) is at $O(\varepsilon)$. By reformulating the NKGE into a relativistic nonlinear Schrödinger equation, we adapt the TSFP method to discretize it numerically. By using the method of mathematical induction to bound the numerical solution, we prove uniform error bounds at $O(h^{m}+\varepsilon^{p-β}τ^2)$ of the TSFP method with $h$ mesh size, $τ$ time step and $m\ge2$ depending on the regularity of the solution. The error bounds are uniformly accurate for the long-time simulation up to the time at $O(\varepsilon^{-β})$ and uniformly valid for $\varepsilon\in(0,1]$. Especially, the error bounds are uniformly at the second order rate for the large time step $τ= O(\varepsilon^{-(p-β)/2})$ in the parameter regime $0\leβ

研究动机与目标

• 分析弱非线性特征为 $\varepsilon^p$ 且初始数据为 $O(1)$ 的非线性克莱因-戈尔登方程（NKGE）的长时间动力学。
• 为时间分裂傅里叶伪谱（TSFP）方法建立统一的误差界，确保对所有 $\varepsilon \in (0,1]$ 和长达 $O(\varepsilon^{-\beta})$ 的长时间区间均有效，其中 $0 \leq \beta \leq p$。
• 将分析扩展至具有空间波长 $O(1)$ 和时间波长 $O(\varepsilon^\beta)$、波速为 $O(\varepsilon^{-\beta})$ 的高度振荡复数 NKGE。
• 证明在 $0 \leq \beta < p$ 的情形下，TSFP 方法即使在时间步长 $\tau = O(\varepsilon^{-(p-\beta)/2})$ 下仍保持时间方向的二阶收敛性。

提出的方法

• 将 NKGE 重新表述为相对论性非线性薛定谔方程，以简化对弱非线性动力学的分析。
• 应用时间分裂傅里叶伪谱（TSFP）方法对重述后的方程进行空间和时间的离散化。
• 使用数学归纳法控制数值解，并在 $H^1$-范数下推导统一的误差估计。
• 建立误差界为 $O(h^m + \varepsilon^{p-\beta}\tau^2)$，其中 $h$ 为空间网格尺寸，$\tau$ 为时间步长，$m \geq 2$ 取决于解的正则性。
• 证明误差界对所有 $\varepsilon \in (0,1]$ 和长时间区间 $T = O(\varepsilon^{-\beta})$ 均保持一致，确保在弱非线性区域的鲁棒性。
• 将分析扩展至波速为 $O(\varepsilon^{-\beta})$、空间波长为 $O(1)$、时间波长为 $O(\varepsilon^\beta)$ 的复数高度振荡 NKGE。

实验结果

研究问题

• RQ1 能否为弱非线性非线性克莱因-戈尔登方程的长时间模拟中使用的时间分裂谱方法建立统一的误差界？
• RQ2 当时间步长按 $\tau = O(\varepsilon^{-(p-\beta)/2})$ 缩放时，TSFP 方法是否仍保持时间方向的二阶精度，且适用于 $0 \leq \beta < p$ 的情形？
• RQ3 在 $\beta$ 和 $p$ 的不同参数区间下，误差行为如何变化，特别是在时间步长大小和总步数方面？
• RQ4 TSFP 方法及其误差界能否扩展至波速为 $O(\varepsilon^{-\beta})$、时间振荡为 $O(\varepsilon^\beta)$ 的复数高度振荡 NKGE？
• RQ5 误差界是否对所有 $\varepsilon \in (0,1]$ 均保持一致，包括 $\varepsilon \to 0^+$ 的极限情形，适用于长时间模拟？

主要发现

• TSFP 方法在长达 $O(\varepsilon^{-\beta})$ 时间范围内对 NKGE 实现了统一的误差界，大小为 $O(h^m + \varepsilon^{p-\beta}\tau^2)$，且对所有 $\varepsilon \in (0,1]$ 均有效。
• 当 $0 \leq \beta < p$ 时，即使采用大时间步长 $\tau = O(\varepsilon^{-(p-\beta)/2})$，方法仍保持时间方向的二阶收敛性，证实了在弱非线性区域的鲁棒性。
• 数值结果验证了在所有测试的 $\varepsilon$ 值（包括 $\varepsilon \ll 1$）下，时间方向保持二阶收敛性，空间方向保持谱精度，适用于实数与复数 NKGE。
• 当 $\beta = 1$ 时，方法在所有 $\varepsilon$ 下均显示统一的二阶时间收敛性；而当 $\beta = 2$ 和 $\beta = 3$ 时，仅在 $\tau \lesssim \varepsilon^{3/2}$ 和 $\tau \lesssim \varepsilon^3$ 时观察到二阶收敛性。
• 方法的总计算成本在 $\beta \leq p/2$ 时为 $O(N\ln N)$，但在 $\beta > p/2$ 时增加至 $O(N\ln N / \varepsilon^{(3\beta-p)/2})$，反映出高度振荡区域的挑战。
• 误差界对所有 $\varepsilon \in (0,1]$ 均保持一致，确保了在弱非线性和小数据极限下的可靠性。