Decreasing Computing Time with Symplectic Correctors in Adaptive Timestepping Routines
arXiv:1011.3830 · doi:10.1088/0004-6256/141/1/3
Abstract
It has previously been shown that varying the numerical timestep during a symplectic orbital integration leads to a random walk in energy and angular momentum, destroying the phase space-conserving property of symplectic integrators. Here we show that when altering the timestep symplectic correctors can be used to reduce this error to a negligible level. Furthermore, these correctors can also be employed to avoid a large error introduction when changing the Hamiltonian's partitioning. We have constructed a numerical integrator using this technique that is nearly as accurate as widely used fixed-step routines. In addition, our algorithm is drastically faster for integrations of highly eccentricitic, large semimajor axis orbits, such as those found in the Oort Cloud.
Accepted to AJ, 29 pages, 8 figures
References in corpus (4)
Cited by in corpus (12)
- Oort Cloud and Scattered Disc formation during a late dynamical instability in the Solar System
- TTVFast: An efficient and accurate code for transit timing inversion problems
- Sedna and the Oort Cloud Around a Migrating Sun
- An Oort cloud origin for the high-inclination, high-perihelion Centaurs
- A Possible Divot in the Size Distribution of the Kuiper Belt's Scattering Objects
- The Canada-France Ecliptic Plane Survey (CFEPS) - High Latitude Component
- Reassessing the formation of the Inner Oort cloud in an embedded star cluster II: Probing the inner edge
- An updated estimate of the number of Jupiter-family comets using a simple fading law
- Observational Signatures of a Massive Distant Planet on the Scattering Disk
- Exploring the role of the Sun's motion in terrestrial comet impacts
- An Oort cloud origin of the Halley-type comets
- The Influence of Outer Solar System Architecture on the Structure and Evolution of the Oort Cloud