Numerical modelling of Kuiper belt objects’ dynamics – limitations
AbstractThe investigation of KBOs’ dynamics is based on numerical orbital integrations on extremly long time scales due to orbital evolution of particles. The evolution of KBOs to JFCs needs a time-span of order of 109years. Such a long time of integration affects errors. So the question arises what is the boundary of an integration time to distinguish the physical solution from numerical noise and what it depends on. This paper presents numerical integrations of less than 150 massless test particles in the model of the Solar System which consists of 4 giant planets and the central mass. For each test particle computations were repeated at least twice on different computers and using two different methods of integration. The results show that an increase of errors in a solution depends on the eccentricity and the inclination of an orbit. The estimated maximum time-span of integration is of the order of 10 million years for highly elliptic orbits (e 0.6) and up to 125 million years for quasi-circular orbits (for particular model of the Solar System with orbits of massless objects outside Neptune's orbit). After long time-span of integration (120-130 Myrs) the solution can be completely chaotic. It cannot be stated unequivocally that this is one of the possible particle's paths or that this is just a numerical noise. So a different way of studying KBOs’ and SP comets’ dynamical evolution is needed. The integration of equations of motion between particular phases of objects which are considered as comets in different phases of their lives (KBOs − Centaurs − Comets − possibly extinct Comets) could be the new way of studying the dynamical evolution of SP comets.