Closed Newton–Cotes trigonometrically-fitted formulae of high order for long-time integration of orbital problems |
| |
Authors: | TE Simos |
| |
Institution: | aLaboratory of Computational Sciences, Department of Computer Science and Technology, Faculty of Science and Technology, University of Peloponnese, GR-221 00 Tripolis, Greece |
| |
Abstract: | The connection between closed Newton–Cotes, trigonometrically-fitted differential methods and symplectic integrators is studied in this paper. Several one-step symplectic integrators have been obtained based on symplectic geometry, as is shown in the literature. However, the study of multi-step symplectic integrators is very limited. The well-known open Newton–Cotes differential methods are presented as multilayer symplectic integrators by Zhu et al. W. Zhu, X. Zhao, Y. Tang, Journal of Chem. Phys. 104 (1996), 2275]. The construction of multi-step symplectic integrators based on the open Newton–Cotes integration methods is investigated by Chiou and Wu J.C. Chiou, S.D. Wu, Journal of Chemical Physics 107 (1997), 6894]. The closed Newton–Cotes formulae are studied in this paper and presented as symplectic multilayer structures. We also develop trigonometrically-fitted symplectic methods which are based on the closed Newton–Cotes formulae. We apply the symplectic schemes in order to solve Hamilton’s equations of motion which are linear in position and momentum. We observe that the Hamiltonian energy of the system remains almost constant as the integration proceeds. Finally we apply the new developed methods to an orbital problem in order to show the efficiency of this new methodology. |
| |
Keywords: | Numerical methods Orbital problems Closed Newton– Cotes differential methods Symplectic integrators Trigonometric fitting Energy preservation |
本文献已被 ScienceDirect 等数据库收录! |
|