Loading...
The adaptive Verlet method
Huang, Weizhang Leimkuhler, Benedict J.
Huang, Weizhang
Leimkuhler, Benedict J.
Citations
Altmetric:
Abstract
We discuss the integration of autonomous Hamiltonian systems via dynamical rescaling of the vector field (reparameterization of time). Appropriate rescalings (e.g., based on normalization of the vector field or on minimum particle separation in an N-body problem) do not alter the time-reversal symmetry of the flow, and it is desirable to maintain this symmetry under discretization. For standard form mechanical systems without rescaling, this can be achieved by using the explicit leapfrog--Verlet method; we show that explicit time-reversible integration of the reparameterized equations is also possible if the parameterization depends on positions or velocities only. For general rescalings, a scalar nonlinear equation must be solved at each step, but only one force evaluation is needed. The new method also conserves the angular momentum for an N-body problem. The use of reversible schemes, together with a step control based on normalization of the vector field (arclength reparameterization), is demonstrated in several numerical experiments, including a double pendulum, the Kepler problem, and a three-body problem.
Description
This is the published version, also available here: http://dx.doi.org/10.1137/S1064827595284658.
Date
1997-01-01
Journal Title
Journal ISSN
Volume Title
Publisher
Society for Industrial and Applied Mathematics
Collections
Research Projects
Organizational Units
Journal Issue
Keywords
time-reversible methods, symplectic methods, Hamitonian systems, variable stepsize methods, Verlet, leapfrong, N-Body problems
Citation
Huang, Weizhang & Leimkuhler, Benedict. "The adaptive Verlet method." SIAM J. Sci. Comput., 18(1), 239–256. (18 pages). http://dx.doi.org/10.1137/S1064827595284658.