Lie-poisson integrators for a rigid satellite on a circular orbit

Abstract

In the last two decades, many structure preserving numerical methods like Poisson integrators have been investigated in numerical studies. Since the structure matrices are different in many Poisson systems, no integrator is known yet to preserve the Poisson structure of any Poisson system. In the present paper, we propose Lie-Poisson integrators for Lie{Poisson systems whose structure matrix is different from the ones studied before. In particular, explicit Lie-Poisson integrators for the equations of rotational motion of a rigid body (the satellite) on a circular orbit around affixed gravitational center have been constructed based on the splitting. The splitted parts have been composed by a first, a second and a third order compositions. It has been shown that the proposed schemes preserve the quadratic invariants of the equation. Numerical results reveal the preservation of the energy and agree with the theoretical treatment that the invariants lie on the sphere in long-term with different orders of accuracy.