The Use of Kepler Solver in Numerical Integrations of Quasi-Keplerian Orbits

Chen Deng,Xin Wu,Enwei Liang
DOI: https://doi.org/10.1093/mnras/staa1753
IF: 4.8
2020-01-01
Monthly Notices of the Royal Astronomical Society
Abstract:A Kepler solver is an analytical method used to solve a two-body problem. Inthis paper, we propose a new correction method by slightly modifying the Keplersolver. The only change to the analytical solutions is that the obtainment ofthe eccentric anomaly relies on the true anomaly that is associated to a unitradial vector calculated by an integrator. This scheme rigorously conserves allintegrals and orbital elements except the mean longitude. However, the Keplerenergy, angular momentum vector and Laplace-Runge-Lenz vector for perturbedKepler problems are slowly-varying quantities. However, their integralinvariant relations give the quantities high-precision values that directlygovern five slowly-varying orbital elements. These elements combined with theeccentric anomaly determine the desired numerical solutions. The newly proposedmethod can considerably reduce various errors for a post-Newtonian two-bodyproblem compared with an uncorrected integrator, making it suitable for adissipative two-body problem. Spurious secular changes of some elements orquasi-integrals in the outer solar system may be caused by short integrationtimes of the fourth-order Runge-Kutta algorithm. However, they can beeliminated in a long integration time of 10^8 years by the proposed method,similar to Wisdom-Holman second-order symplectic integrator. The proposedmethod has an advantage over the symplectic algorithm in the accuracy but givesa larger slope to the phase error growth.
What problem does this paper attempt to address?