Symplectic integrators: Rotations and roundoff errors

DOI10.1023/A:1008268823420MaRDI QIDQ1286873zbMATH OpenFDO

Authors Jean-Marc Petit

Publication date 2 December 1999

Published in Celestial Mechanics and Dynamical Astronomy (Search for Journal in Brave)

Full work available at URL https://doi.org/10.1023/A:1008268823420

zbMATH Keywords

combination of shearselongated rotating primaryinvariant toriKAM-like theorem for twist mapsKeplerian orbitself-adaptative Bulirsch-Stoer integrator

Mathematics Subject Classification ID

Computational methods for problems pertaining to mechanics of particles and systems (70-08) Celestial mechanics (70F15)

Review text

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. 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 terms of computing time. Next we test the shear combination together with symplectic integrator of order 2, 4, and 6 on a Keplerian orbit.

Recommendations

• Symplectic integrators and their application to dynamical astronomy
• Explicit symplectic integrator for rotating satellites
• Recent progress in the theory and application of symplectic integrators
• scientific article; zbMATH DE number 1453822
• Symplectic integration with floating-point arithmetic and other approximations

Cited in

(6)

• Good rotations
• FOURTH ORDER SYMPLECTIC INTEGRATION WITH REDUCED PHASE ERROR
• The accuracy of symplectic integrators
• Dealing with Parasitic Behaviour in G-Symplectic Integrators
• Explicit symplectic integrator for rotating satellites
• Achieving Brouwer's law with implicit Runge-Kutta methods

This page was built for publication: Symplectic integrators: Rotations and roundoff errors

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q1286873)