Journal
APPLIED NUMERICAL MATHEMATICS
Volume 68, Issue -, Pages 58-72Publisher
ELSEVIER
DOI: 10.1016/j.apnum.2013.01.003
Keywords
Symplectic integrators; Splitting methods; Near-integrable systems; N-body problems
Categories
Funding
- Ministerio de Ciencia e Innovacion (Spain) [MTM2010-18246-C03]
- FEDER Funds of the European Union
- FP7 GTSnext project
Ask authors/readers for more resources
We present new splitting methods designed for the numerical integration of near-integrable Hamiltonian systems, and in particular for planetary N-body problems, when one is interested in very accurate results over a large time span. We derive in a systematic way an independent set of necessary and sufficient conditions to be satisfied by the coefficients of splitting methods to achieve a prescribed order of accuracy. Splitting methods satisfying such (generalized) order conditions are appropriate in particular for the numerical simulation of the Solar System described in Jacobi coordinates. We show that, when using Poincare Heliocentric coordinates, the same order of accuracy may be obtained by imposing an additional polynomial equation on the coefficients of the splitting method. We construct several splitting methods appropriate for each of the two sets of coordinates by solving the corresponding systems of polynomial equations and finding the optimal solutions. The experiments reported here indicate that the efficiency of our new schemes is clearly superior to previous integrators when high accuracy is required. (C) 2013 IMACS. Published by Elsevier B.V. All rights reserved.
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