For each celestial body of multi-planet systems, there are two slowly varying quantities or quasi-integrals, Kepler energy and Laplace integral, which are closely associated with the orbital semimajor axis and eccentricity, respectively. To correct numerical errors of the quantities, we give an extension of Nacozy's approach and develop a new manifold correction method, where corresponding reference values of the quantities at every integration step are obtained from integral invariant relations, and only velocity corrections are used to approximately satisfy the two quasi-integrals. As a result, the scheme does enhance the quality of the integration by significantly raising the accuracy of the two elements. Especially, it is superior to the existing dual scaling method in the improvement of eccentricity in general when the adopted integrator provides a sufficient precision to the eccentricity.

No abstract

When conventional integrators like Runge—Kutta-type algorithms are used, numerical errors can make an orbit deviate from a hypersurface determined by many constraints, which leads to unreliable numerical solutions. Scaling correction methods are a powerful tool to avoid this. We focus on their applications, and also develop a family of new velocity multiple scaling correction methods where scale factors only act on the related components of the integrated momenta. They can preserve exactly some first integrals of motion in discrete or continuous dynamical systems, so that rapid growth of roundoff or truncation errors is suppressed significantly.