[论文解读] Long time error analysis of the fourth-order compact finite difference methods for the nonlinear Klein-Gordon equation with weak nonlinearity

Yue Feng|arXiv (Cornell University)|Mar 9, 2020

Numerical methods for differential equations参考文献 42被引用 4

一句话总结

本文提出了一种四阶紧凑有限差分（4cFD）方法，用于求解具有弱非线性的非线性克莱因-戈尔登方程，其中非线性项按 $\varepsilon^p$ 缩放，且 $\varepsilon \in (0,1]$。在长达 $O(\varepsilon^{-p})$ 的长时间范围内，推导出严格的误差界，表明最优精度要求 $h = O(\varepsilon^{p/4})$ 且 $\tau = O(\varepsilon^{p/2})$，展示了该方法在空间分辨率上优于二阶方法，并通过数值实验验证了收敛速率的精确性。

ABSTRACT

We present the fourth-order compact finite difference (4cFD) discretizations for the long time dynamics of the nonlinear Klein-Gordon equation (NKGE), while the nonlinearity strength is characterized by $\varepsilon^p$ with a constant $p \in \mathbb{N}^+$ and a dimensionless parameter $\varepsilon \in (0, 1]$. Based on analytical results of the life-span of the solution, rigorous error bounds of the 4cFD methods are carried out up to the time at $O(\varepsilon^{-p})$. We pay particular attention to how error bounds depend explicitly on the mesh size $h$ and time step $τ$ as well as the small parameter $\varepsilon \in (0, 1]$, which indicate that, in order to obtain `correct' numerical solutions up to the time at $O(\varepsilon^{-p})$, the $\varepsilon$-scalability (or meshing strategy requirement) of the 4cFD methods should be taken as: $h = O(\varepsilon^{p/4})$ and $τ= O(\varepsilon^{p/2})$. It has better spatial resolution capacity than the classical second order central difference methods. By a rescaling in time, it is equivalent to an oscillatory NKGE whose solution propagates waves with wavelength at $O(1)$ in space and $O(\varepsilon^p)$ in time. It is straightforward to get the error bounds of the oscillatory NKGE in the fixed time. Finally, numerical results are provided to confirm our theoretical analysis.

研究动机与目标

• 开发并分析应用于弱非线性非线性克莱因-戈尔登方程（NKGE）的四阶紧凑有限差分（4cFD）方法的长时间误差界，该方程的非线性特征由 $\\varepsilon^p$ 表征。
• 确定为在长达 $O(\\varepsilon^{-p})$ 时间内保持精度所需的网格策略（$h$ 和 $\\tau$）。
• 建立 4cFD 方法的 $\varepsilon$-可扩展性，表明其在空间分辨率上优于经典二阶中心差分方法。
• 通过时间重标度将误差分析扩展至振荡形式的 NKGE，从而实现固定时间的误差界。
• 通过数值实验验证理论误差估计，并确认其精确性。

提出的方法

• 将 4cFD 方法应用于具有弱非线性 $\varepsilon^p u^{p+1}$ 的 NKGE，采用四阶精度的空间模板和二阶时间离散化。
• 应用时间重标度变换，将原始 NKGE 转化为具有 $O(1)$ 空间波长和 $O(\varepsilon^p)$ 时间波长的振荡形式。
• 为 4cFD 格式在长达 $O(\varepsilon^{-p})$ 的长时间内推导出严格的误差界，其显式依赖于 $h$、$\tau$ 和 $\varepsilon$。
• 分析利用能量方法和渐近展开，以控制长时间积分区间内误差的增长。
• 将方法扩展至通过时间重标度获得的振荡 NKGE 形式，从而可直接推导出固定时间的误差界。
• 数值验证采用高精度的指数波积分傅里叶伪谱方法作为参考解，以计算空间和时间误差。

实验结果

研究问题

• RQ1 对于弱非线性的非线性克莱因-戈尔登方程，四阶紧凑有限差分方法的长时间误差行为如何？
• RQ2 为在 $O(\varepsilon^{-p})$ 时间范围内保持精度，所需的网格尺寸 $h$ 和时间步长 $\tau$ 是什么？
• RQ3 与经典二阶中心差分方法相比，4cFD 方法在此问题上的空间分辨率如何？
• RQ4 能否将误差界扩展至通过时间重标度获得的 NKGE 振荡形式？
• RQ5 数值实验是否证实了理论误差估计的精确性？

主要发现

• 对于固定的 $\varepsilon > 0$，4cFD 方法在空间上达到四阶收敛，在时间上达到二阶收敛，验证了其设计精度阶。
• 在高度振荡的参数区域（$\varepsilon \ll 1$）中，仅当 $h \lesssim \varepsilon^{p/4}$ 且 $\tau \lesssim \varepsilon^{p/2}$ 时才能实现最优收敛，验证了理论的 $\varepsilon$-可扩展性。
• 所需网格尺寸 $h = O(\varepsilon^{p/4})$ 提供了优于经典二阶中心差分方法的空间分辨率，后者为达到相似精度需 $h = O(\varepsilon^{p/2})$。
• 数值结果证实误差界是精确的，收敛速率与理论预测在不同 $\varepsilon$ 值下完全匹配。
• 振荡形式的 NKGE 允许直接推导出固定时间的误差界，简化了长时间动力学的分析。
• 针对 $p=1$ 和 $p=2$ 的数值实验表明，在所推导的网格策略下，4cFD 方法在长达 $O(\varepsilon^{-p})$ 的时间范围内仍能保持高精度。