N-body problem solution with composition numerical integration methods

The subject of this paper is the experimental study of the composition numerical integration methods while solving the N-body problem. Various composition methods based on symmetric second order Verlet method are considered. Step size control algorithm for these methods is described. The principles...

Full description

Saved in:
Bibliographic Details
Published in2016 XIX IEEE International Conference on Soft Computing and Measurements (SCM) pp. 196 - 198
Main Authors Andreev, V. S., Goryainov, S. V., Krasilnikov, A. V.
Format Conference Proceeding
LanguageEnglish
Published IEEE 01.05.2016
Subjects
Adaptation models
adaptive stepsize
Algorithm design and analysis
composition method
dynamical systems simulation
Earth
gravitational N-body problem
Heuristic algorithms
Mathematical model
Numerical models
reversibitity
Stability analysis
Verlet method
Online AccessGet full text

Cover

More Information
Summary:The subject of this paper is the experimental study of the composition numerical integration methods while solving the N-body problem. Various composition methods based on symmetric second order Verlet method are considered. Step size control algorithm for these methods is described. The principles for experimental estimation of the time reversibility for composition methods with fixed and adaptive integration steps is described. The results six-body problem simulation by composition method of accuracy order 6 and Runge-Kutta 45 method are compared. Conclusions about the effectiveness of proposed algorithms for long-term simulation of the chosen problem are given. Study shows that semi-implicit composition methods have greater numerical stability for N-body problem simulation than the explicit Runge-Kutta methods. It is shown that the semi-implicit composition methods with fixed integration step are time-reversible, but they lose this property in a case of adaptive step size.
DOI:10.1109/SCM.2016.7519726