Linearly implicit and high-order energy-preserving relaxation schemes for highly oscillatory Hamiltonian systems.

In this paper, a family of novel energy-preserving schemes are presented for numerically solving highly oscillatory Hamiltonian systems. These schemes are constructed by using the relaxation idea in the extrapolated Runge–Kutta (ERK) methods. After obtaining the relaxation parameter, it is shown that the methods can be arbitrarily high order accurate, linearly implicit and keep the original discrete energy conserved, while the previous energy-preserving schemes for highly oscillatory Hamiltonian systems are usually fully implicit. Numerical comparisons with various typical energy-preserving schemes are presented. The numerical results show that the proposed schemes are highly competitive and effective.