An Averaging Technique for Highly Oscillatory Hamiltonian Problems

In this paper, we are concerned with the numerical solution of highly oscillatory Hamiltonian systems with a stiff linear part. We construct an averaged system whose solution remains close to the exact one over bounded time intervals, possesses the same adiabatic and Hamiltonian invariants as the original system, and is nonstiff. We then investigate its numerical approximation through a method which combines a symplectic integration scheme and an acceleration technique for the evaluation of time averages developed in [E. Cances et al., Numer. Math., 100 (2005), pp. 211-232]. Eventually, we demonstrate the efficiency of our approach on two test problems with one or several frequencies.

[1]  V. Arnold SMALL DENOMINATORS AND PROBLEMS OF STABILITY OF MOTION IN CLASSICAL AND CELESTIAL MECHANICS , 1963 .

[2]  Antonio Giorgilli,et al.  Realization of holonomic constraints and freezing of high frequency degrees of freedom in the light of classical perturbation theory. Part I , 1987 .

[3]  V. Arnold Mathematical Methods of Classical Mechanics , 1974 .

[4]  J. Moser On invariant curves of area-preserving mappings of an anulus , 1962 .

[5]  P. Deuflhard A study of extrapolation methods based on multistep schemes without parasitic solutions , 1979 .

[6]  F. Legoll,et al.  Integrators for highly oscillatory Hamiltonian systems: an homogenization approach , 2009 .

[7]  Antonio Giorgilli,et al.  Realization of holonomic constraints and freezing of high frequency degrees of freedom in the light of classical perturbation theory. Part II , 1987 .

[8]  Ernst Hairer,et al.  Modulated Fourier Expansions of Highly Oscillatory Differential Equations , 2003, Found. Comput. Math..

[9]  E. Hairer,et al.  Geometric Numerical Integration , 2022, Oberwolfach Reports.

[10]  W. Gautschi Numerical integration of ordinary differential equations based on trigonometric polynomials , 1961 .

[11]  Ernst Hairer,et al.  Asymptotic expansions and backward analysis for numerical integrators , 2000 .

[12]  Marlis Hochbruck,et al.  Error analysis of exponential integrators for oscillatory second-order differential equations , 2006 .

[13]  E. Hairer,et al.  Geometric Numerical Integration: Structure Preserving Algorithms for Ordinary Differential Equations , 2004 .

[14]  Ernst Hairer,et al.  Long-Time Energy Conservation of Numerical Methods for Oscillatory Differential Equations , 2000, SIAM J. Numer. Anal..

[15]  Robert D. Skeel,et al.  Long-Time-Step Methods for Oscillatory Differential Equations , 1998, SIAM J. Sci. Comput..

[16]  Marlis Hochbruck,et al.  A Gautschi-type method for oscillatory second-order differential equations , 1999, Numerische Mathematik.

[17]  S. Reich Backward Error Analysis for Numerical Integrators , 1999 .

[18]  E. Hairer,et al.  Numerical Energy Conservation for Multi-Frequency Oscillatory Differential Equations , 2005 .