Astronomy and Astrophysics – Astronomy
Scientific paper
Jan 1998
adsabs.harvard.edu/cgi-bin/nph-data_query?bibcode=1998cemda..70....1p&link_type=abstract
Celestial Mechanics and Dynamical Astronomy, v. 70, Issue 1, p. 1-21 (1998).
Astronomy and Astrophysics
Astronomy
3
Symplectic Integrators, Roundoff Error, Celestial Mechanics
Scientific paper
We investigate the numerical implementation of a symplectic integrator combined with a rotation (as in the case of an elongated rotating primary). We show that a straightforward implementation of the rotation as a matrix multiplication destroys the conservative property of the global integrator, due to roundoff errors. According to Blank et al. (1997), there exists a KAM-like theorem for twist maps, where the angle of rotation is a function of the radius. This theorem proves the existence of invariant tori which confine the orbit and prevent shifts in radius. We replace the rotation by a twist map or a combination of shears that display the same kind of behaviour and show that we are able not only to recover the conservative properties of the rotation, but also make it more efficient in term of computing time. Next we test the shear combination together with symplectic integrator of order 2, 4, and 6 on a Keplerian orbit. The resulting integrator is conservative down to the roundoff errors. No linear drift of the energy remains, only a divergence as the square root of the number of iterations is to be seen, as in a random walk. We finally test the three symplectic integrators on a real case problem of the orbit of a satellite around an elongated irregular fast rotating primary. We compare these integrators to the well-known general purpose, self-adaptative Bulirsch Stoer integrator. The sixth order symplectic integrator is more accurate and faster than the Bulirsch Stoer integrator. The second- and fourth- order integrators are faster, but of interest only when extreme speed is mandatory.
No associations
LandOfFree
Symplectic Integrators: Rotations and Roundoff Errors does not yet have a rating. At this time, there are no reviews or comments for this scientific paper.
If you have personal experience with Symplectic Integrators: Rotations and Roundoff Errors, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Symplectic Integrators: Rotations and Roundoff Errors will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-1654577