Long time error analysis of the fourth-order compact finite difference methods for the nonlinear Klein-Gordon equation with weak nonlinearity

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 $\tau$ 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 $\tau = 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.