Journal Article

Symplectic integrators in the shearing sheet

Hanno Rein and Scott Tremaine

in Monthly Notices of the Royal Astronomical Society

Published on behalf of The Royal Astronomical Society

Volume 415, issue 4, pages 3168-3176
Published in print August 2011 | ISSN: 0035-8711
Published online August 2011 | e-ISSN: 1365-2966 | DOI: http://dx.doi.org/10.1111/j.1365-2966.2011.18939.x
Symplectic integrators in the shearing sheet

Show Summary Details

Preview

The shearing sheet is a model dynamical system that is used to study the small-scale dynamics of astrophysical discs. Numerical simulations of particle trajectories in the shearing sheet usually employ the leapfrog integrator, but this integrator performs poorly because of velocity-dependent (Coriolis) forces. We describe two new integrators for this purpose; both are symplectic, time-reversible and second-order accurate, and can easily be generalized to higher orders. Moreover, both the integrators are exact when there are no small-scale forces such as mutual gravitational forces between disc particles. In numerical experiments these integrators have errors that are often several orders of magnitude smaller than competing methods. The first of our new integrators (‘SEI’) is well suited for discs in which the typical interparticle separation is large compared to the particles’ Hill radii (e.g., planetary rings), and the second one (‘SEKI’) is designed for discs in which the particles are on bound orbits or the separation is smaller than the Hill radius (e.g. irregular satellites of the giant planets).

Keywords: methods: numerical; celestial mechanics; planets and satellites: dynamical evolution and stability; planets and satellites: formation; planets and satellites: rings

Journal Article.  5153 words.  Illustrated.

Subjects: Astronomy and Astrophysics

Full text: subscription required

How to subscribe Recommend to my Librarian

Users without a subscription are not able to see the full content. Please, subscribe or login to access all content.