Department of Computer Science, University of Illinois at Urbana-Champaign, 1304 W. Springfield Ave., Urbana, Illinois 61801-2987, USA
Abstract:
A Hamiltonian system possesses dynamics (e.g. preservation of volume in phase space and symplectic structure) that call for special numerical integrators, namely canonical methods. Recent research on this aspect have shown that canonical numerical integrators may be needed for Hamiltonian systems. In this paper, we focus on numerical experiments that compare canonical and non-canonical numerical integrators. Test problems are taken from different areas in physical sciences. These experiments help to buttress the claims that canonical numerical integrators give results that mimic the qualitative behavior of the original system and that canonical numerical integrators are suitable for long time integrations. Our experiments indicate that higher-order canonical methods allow for larger timestep than lower-order canonical methods.