New families of symplectic splitting methods for numerical integration in dynamical astronomy
[EN] 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...
| Autores: | , , , , , |
|---|---|
| Tipo de recurso: | artículo |
| Fecha de publicación: | 2013 |
| País: | España |
| Institución: | Universitat Politècnica de València (UPV) |
| Repositorio: | RiuNet. Repositorio Institucional de la Universitat Politécnica de Valéncia |
| Idioma: | inglés |
| OAI Identifier: | oai:riunet.upv.es:10251/38064 |
| Acceso en línea: | https://riunet.upv.es/handle/10251/38064 |
| Access Level: | acceso abierto |
| Palabra clave: | Symplectic integrators Splitting methods Near-integrable systems N-body problems MATEMATICA APLICADA |
| Sumario: | [EN] 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 Poincaré 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. |
|---|