T. E. Simos¹*

¹ Laboratory of Computational Sciences, Department of Computer Science and Technology, Faculty of Sciences and Technology, University of Peloponnese, GR-221 00 Tripolis, Greece (tsimos@mail.ariadne-t.gr).

Received 2005 October 5
Accepted 2006 April 18

En este trabajo se investiga la conexión entre formulas Newton-Cotes, métodos diferenciales por ajustes trigonométricos e integradores simplécticos. Se conoce, a través de la literatura, que varios integradores simplécticos de un paso han sido obtenidos basándose en geometría simpléctica. Sin embargo, la investigación de integradores simplécticos multicapa es muy pobre. Zhu et al. (1996) presentaron los conocidos métodos diferenciales Newton-Cotes abiertos como integradores simplécticos multicapa. También Chiou & Wu (1997) investigaron la construcción de integradores simplécticos multicapa basándose en los métodos de integración abierta Newton-Cotes. En este trabajo investigamos las fórmulas cerradas Newton-Cotes y las escribimos como estructuras simplécticas multicapa. Después de esto, construimos métodos simplécticos por ajustes trigonométricos, los cuales se basan en las formulas cerradas Newton-Cotes. Aplicamos los esquemas simplécticos para resolver las ecuaciones de movimiento de Hamilton que son lineales en posición y momento. Observamos que la energía hamiltoniana del sistema permanece casi constante a medida que la integración avanza.

ABSTRACT

The connection between closed Newton-Cotes, trigonometrically-fitted differential methods and symplectic integrators is investigated in this paper. It is known from the literature that several one step symplectic integrators have been obtained based on symplectic geometry. However, the investigation of multistep symplectic integrators is very poor. Zhu et al. (1996) presented the well known open Newton-Cotes differential methods as multilayer symplectic integrators. Also, Chiou & Wu (1997) investigated the construction of multistep symplectic integrators based on the open Newton-Cotes integration methods. In this paper we investigate the closed Newton-Cotes formulae and we write them as symplectic multilayer structures. After this we construct 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 integration procceeds.

Key Words: CELESTIAL MECHANICS - METHODS: NUMERICAL.

DESCARGAR ARTÍCULO EN FORMATO PDF

REFERENCES

Anastassi, Z. A., & Simos, T. E. 2004, NewA, 10(1), 31        [ Links ]

----------. 2005, NewA, 10,(4), 301        [ Links ]

Chiou, J. C., & Wu, S. D. 1997, J. of Chem. Phys., 107, 6894        [ Links ]

Monovasilis, Th., Kalogiratou, Z., & Simos, T. E. 2004, Appl. Num. Anal. Comp. Math., 1(1), 195        [ Links ]

----------. 2005, Appl. Num. Anal. Comp. Math., 2(2), 238        [ Links ]

Psihoyios, G., & Simos, T. E. 2003, NewA, 8(7), 679        [ Links ]

----------. 2004a, Appl. Num. Anal. Comp. Math., 1(1), 205        [ Links ]

----------. 2004b, Appl. Num. Anal. Comp. Math., 1(1), 216        [ Links ]

Sanz-Serna, J. M., & Calvo, M. P. 1994, in Numerical Hamiltonian Problem (London: Chapman and Hall)        [ Links ]

Simos, T. E. 1996, International J. of Modern Phys., C, 7, 825        [ Links ]

----------. 1998a, International J. of Modern Phys., C,9, 1055        [ Links ]

----------. 1998b, International J. of Modern Phys., C, 9, 271        [ Links ]

Simos, T. E. 2000, International J. of Modern Phys., C, 11, 79        [ Links ]

----------. 2002a, in Numerical Methods for 1D, 2D, and 3D Differential Equations Arising in Chemical Problems, Chemical Modelling: Applications and Theory (The Royal Society of Chemistry), 170        [ Links ]

----------. 2002b, NewA, 7(1), 1        [ Links ]

----------. 2003, NewA, 8(5), 391        [ Links ]

----------. 2004a, NewA, 9(6), 409        [ Links ]

----------. 2004b, NewA, 9(1), 59        [ Links ]

----------. 2005, Computing Lett., 1(1), 37        [ Links ]

Stiefel, E., & Bettis, D. G. 1969, Numer. Math., 13, 154        [ Links ]

Van Daele, M., & Vanden Berghe, G. 2004, Appl. Num.Anal. Comp. Math., 1(2), 353        [ Links ]

Vanden Berghe, G., Van Daele, M., & Vande Vyver, H. 2004, Appl. Num. Anal. Comp. Math., 1(1), 49        [ Links ]

Vlachos D. S., & Simos, T. E. 2004, Appl. Num. Anal.Comp. Math., 1(2), 540        [ Links ]

Zhu, W., Zhao, X., & Tang, Y. 1996, J. of Chem. Phys., 104, 2275        [ Links ]

NOTES:

* Active Member of the European Academy of Sciences and Arts, Corresponding Member of the European Academy of Sciences, and Corresponding Member of European Academy of Arts, Sciences, and Humanities.