首页 | 本学科首页   官方微博 | 高级检索  
     检索      


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 等数据库收录!
设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号