Geometric two-scale integrators for highly oscillatory system: uniform accuracy and near conservations

Abstract

In this paper, we consider a class of highly oscillatory Hamiltonian systems which involve a scaling parameter $\epsilon \in (0,1]$. The problem arises from many physical models in some limit parameter regime or from some time-compressed perturbation problems. The solution of the model exhibits rapid temporal oscillations with $\mathcal{O}(1)$-amplitude and $\mathcal{O}(1/\epsilon )$-frequency, which makes classical numerical methods inefficient. We apply the two-scale formulation approach to the problem and propose two new time-symmetric numerical integrators. The methods are proved to have the uniform second order accuracy for all $\epsilon$ at finite times and some near-conservation laws in long times. Numerical experiments on a Hénon-Heiles model, a nonlinear Schrödinger equation and a charged-particle system illustrate the performance of the proposed methods over the existing ones.

arXiv comment21 pages

References in corpus (3)

• Exponential integrators preserving first integrals or Lyapunov functions for conservative or dissipative systems
• Comparison of numerical methods for the nonlinear Klein-Gordon equation in the nonrelativistic limit regime
• Uniform error bounds of a time-splitting spectral method for the long-time dynamics of the nonlinear Klein-Gordon equation with weak nonlinearity