Abstract



We present an efficient strategy for the numerical propagation of small Solar System objects undergoing close encounters with massive bodies. The trajectory is split in several phases, each of them being the solution of a perturbed twobody problem. Formulations regularized with respect to different primaries are employed in two subsequent phases. In particular, we consider the KustaanheimoStiefel regularization and a novel set of nonsingular orbital elements pertaining to the Dromo family. In order to test the proposed strategy, we perform ensemble propagations in the EarthSun CR3BP using a variable stepsize and order multistep integrator and an improved version of Everhart's RADAU solver of 15th order. By combining the trajectory splitting with regularized equations of motion in shortterm propagations (one year), we gain up to six orders of magnitude in accuracy with respect to the classical Cowell's method for the same computational cost. Moreover, in the propagation of asteroid (99942) Apophis through its 2029 Earth encounter, the position error stays within 100 meters after 100 years. In general, as to improve the performance of regularized formulations, the trajectory must be split between 1.2 and 3 Hill radii from the Earth. We also devise a robust iterative algorithm to stop the integration of regularized equations of motion at a prescribed physical time. The results rigorously hold in the CR3BP, and similar considerations may apply when considering more complex models. The methods and algorithms are implemented in the NAPLES Fortran 2003 code, which is available online as a GitHub repository.  
International

Si 
JCR

Si 
Title

Monthly Notices of the Royal Astronomical Society 
ISBN

00358711 
Impact factor JCR

4,952 
Impact info

Datos JCR del año 2015 
Volume



10.1093/mnras/stx1254 
Journal number


From page

1 
To page

21 
Month

MAYO 
Ranking
