Astronomy and Astrophysics – Astrophysics
Scientific paper
1999-10-14
Astronomy and Astrophysics
Astrophysics
14 pages, 5 figures. Accepted for publication in the Astronomical Journal
Scientific paper
10.1086/301161
Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly with order in timestep, rendering higher-order methods impractical. However, symplectic integrators are often applied to systems in which perturbations between bodies are a small factor of the force due to a dominant central mass. In this case, it is possible to create optimized symplectic algorithms that require fewer substeps per timestep. This is achieved by only considering error terms of order epsilon, and neglecting those of order epsilon^2, epsilon^3 etc. Here we devise symplectic algorithms with 4 and 6 substeps per step which effectively behave as 4th and 6th-order integrators when epsilon is small. These algorithms are more efficient than the usual 2nd and 4th-order methods when applied to planetary systems.
Chambers John E.
Murison Marc Allen
No associations
LandOfFree
Pseudo-High-Order Symplectic Integrators 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 Pseudo-High-Order Symplectic Integrators, we encourage you to share that experience with our LandOfFree.com community. Your opinion is very important and Pseudo-High-Order Symplectic Integrators will most certainly appreciate the feedback.
Profile ID: LFWR-SCP-O-105597