Collision detection for N-body Kepler systems

Visser, P.M.

doi:10.1051/0004-6361/202243754

Collision detection for N-body Kepler systems

Title

Collision detection for N-body Kepler systems

Author

Visser, P.M. (TU Delft Mathematical Physics)

Date

2023

Abstract

Context. In a Keplerian system, a large number of bodies orbit a central mass. Accretion disks, protoplanetary disks, asteroid belts, and planetary rings are examples. Simulations of these systems require algorithms that are computationally efficient. The inclusion of collisions in the simulations is challenging but important. Aims. We intend to calculate the time of collision of two astronomical bodies in intersecting Kepler orbits as a function of the orbital elements. The aim is to use the solution in an analytic propagator (N-body simulation) that jumps from one collision event to the next. Methods. We outline an algorithm that maintains a list of possible collision pairs ordered chronologically. At each step (the soonest event on the list), only the particles created in the collision can cause new collision possibilities. We estimate the collision rate, the length of the list, and the average change in this length at an event, and study the efficiency of the method used. Results. We find that the collision-time problem is equivalent to finding the grid point between two parallel lines that is closest to the origin. The solution is based on the continued fraction of the ratio of orbital periods. Conclusions. Due to the large jumps in time, the algorithm can beat tree codes (octree and k-d tree codes can efficiently detect collisions) for specific systems such as the Solar System with N < 108. However, the gravitational interactions between particles can only be treated as gravitational scattering or as a secular perturbation, at the cost of reducing the time-step or at the cost of accuracy. While simulations of this size with high-fidelity propagators can already span vast timescales, the high efficiency of the collision detection allows many runs from one initial state or a large sample set, so that one can study statistics.

Subject

gravitation
methods: analytical
methods: statistica
celestial mechanics
planets and satellites: formation
protoplanetary disks

To reference this document use:

http://resolver.tudelft.nl/uuid:0c612b70-93ac-4e80-be59-0cc99d10e30e

DOI

https://doi.org/10.1051/0004-6361/202243754

ISSN

0004-6361

Source

Astronomy & Astrophysics, 669, 1-15

Part of collection

Institutional Repository