Integrators for highly oscillatory Hamiltonian systems: an homogenization approach

Abstract : We introduce a systematic way to construct symplectic schemes for the numerical integration of a large class of highly oscillatory Hamiltonian systems. The bottom line of our construction is to consider the Hamilton-Jacobi form of the Newton equations of motion, and to perform a two-scale expansion of the solution, for small times and high frequencies. The approximation obtained for the solution is then used as a generating function, from which the numerical scheme is derived. Several options for the derivation are presented, some of them also giving rise to non symplectic variants. The various integrators obtained are tested and compared to several existing algorithms. The numerical results demonstrate their efficiency.