Conservation properties of numerical integrators for highly oscillatory Hamiltonian systems

Modulated Fourier expansion is used to show long-time near-conservation of the total and oscillatory energies of numerical methods for Hamiltonian systems with highly oscillatory solutions. The numerical methods considered are an extension of the trigonometric methods. A brief discussion of conservation properties in the continuous problem and in the multi-frequency case is also given.

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

[2]  Uri M. Ascher,et al.  On Some Difficulties in Integrating Highly Oscillatory Hamiltonian Systems , 1999, Computational Molecular Dynamics.

[3]  N. H. McClamroch,et al.  Dynamics and control of an elastic dumbbell spacecraft in a central gravitational field , 2003, 42nd IEEE International Conference on Decision and Control (IEEE Cat. No.03CH37475).

[4]  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 .

[5]  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 .

[6]  E. Hairer,et al.  Solving Ordinary Differential Equations I , 1987 .

[7]  Uri M. Ascher,et al.  The Midpoint Scheme and Variants for Hamiltonian Systems: Advantages and Pitfalls , 1999, SIAM J. Sci. Comput..

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

[9]  J. Izaguirre.  Longer Time Steps for Molecular Dynamics , 1999 .

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

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