Efficient Orbit Integration by Linear Transformation for Consistency of Kepler Energy, Full Laplace Integral, and Angular Momentum Vector
Efficient Orbit Integration by Linear Transformation for Consistency of Kepler Energy, Full Laplace Integral, and Angular Momentum Vector
复制标题
DOI:
10.1086/420808
复制
发表时间:
2004-06
期刊:
影响因子:
--
通讯作者:
T. Fukushima
中科院分区:
文献类型:
--
作者:
T. Fukushima
By adopting a general linear transformation as the method of manifold correction, we modify our dual scaling method to integrate quasi-Keplerian orbits numerically. The new method adjusts the integrated position and velocity at each integration step in order to exactly satisfy the relations for the Kepler energy, angular momentum vector, and the full Laplace vector. In the case of no perturbation, the integration errors in all the orbital elements except the mean longitude at the epoch, which grows linearly with time, are reduced to the level of the machine epsilon throughout the integration. For perturbed orbits, the integration errors in position are smaller than with the previous methods of manifold correction. Since its wide applicability is unchanged and the cost of additional computation is similarly negligible, we recommend the new method as the best of our methods of manifold correction.