Journal
MONTHLY NOTICES OF THE ROYAL ASTRONOMICAL SOCIETY
Volume 372, Issue 1, Pages 219-223Publisher
BLACKWELL PUBLISHING
DOI: 10.1111/j.1365-2966.2006.10854.x
Keywords
stellar dynamics; methods : N-body simulations; celestial mechanics
Categories
Ask authors/readers for more resources
Algorithmic regularization uses a transformation of the equations of motion such that the leapfrog algorithm produces exact trajectories for two-body motion as well as regular results in numerical integration of the motion of strongly interacting few-body systems. That algorithm alone is not sufficiently accurate and one must use the extrapolation method for improved precision. This requires that the basic leapfrog algorithm be time-symmetric, which is not directly possible in the case of velocity-dependent forces, but is usually obtained with the help of the implicit mid-point method. Here, we suggest an alternative explicit algorithmic regularization algorithm which can handle velocity-dependent forces. This is done with the help of a generalized mid-point method to obtain the required time symmetry, thus eliminating the need for the implicit mid-point method and allowing the use of extrapolation.
Authors
I am an author on this paper
Click your name to claim this paper and add it to your profile.
Reviews
Recommended
No Data Available