Abstract
The Newton-Störmer/Verlet-leapfrog method (S/V) is a symplectic and symmetric one of order two, which, when applied to separable Hamiltonian dynamical systems, becomes explicit and conserves quadratic first integrals, e.g., the angular momentum in the N-body problem.
As for high accuracy long-time computations required by dynamical systems coming from astronomy, the order two is too low, we consider composed S/V methods up to order 10. Beyond all these remarkable qualities, our numerical experiments on the outer as well as inner Solar system (N-body problem for N larger than 5) show that, in practice, the conservation of the Hamiltonian of the flow, up to a prescribed order, still remains an open task. Specifically, the order at which the Hamiltonian is conserved seems to be quasi independent of the order of the composed method in spite of the fact that the compensated summation, as a specific technique to reduce the round-off errors, was used. In order to clear up this aspect we perform a round-off error analysis of S/V method as a one-step method applied to an arbitrary vector field.
Despite of the fact that we succeeded in improving the classical round-off error estimates, numerical experiments lead to the opinion that the accumulation of round-off errors bears the main responsibility for the lack of accuracy in long time energy conservation.
Authors
Calin-Ioan Gheorghiu,
Adrian-Cristian Muresan
Keywords
References
See the expanding block below.
Cite this paper as
C.I. Gheorghiu, A.C. Mureşan, On the cccuracy of Stoermer/Verlet method as the numerical integrator of the n-body problem – Application to solar system, Rom. Astron. J., 16 (2006) 93-102.
About this paper
Journal
Romanian Astronomical Journal
Publisher Name
Romanian Academy
Paper on the journal website
Print ISSN
1220-5168
Online ISSN
2285-3758
Google Scholar Profile
[1] Applegate, J. H. Douglas, M. R. Guersel, Y. Sussman, G. J. Wisdom, J: 1986, Astron, J., 92, 176
[2] Benetin, G., Giorgilli, A. : 1994, J. Statist. Phys.,74, 1117
[3] Butcher, J.: 2003, Numerical Methods for Ordinary Differential Equations, J. Wiley, Chichester
[4] Dekker, T.J.: 1971, Numer. Math. 18, 224.
[5] Fukushima, T.: 2001, Astron. J., 121, 1768
[6] Gladman, B., Duncan, M., Candy, J: 1991, Celest. Mech., 52, 221
[7] Heirer, E., Lubicj, C. Wanner, G.: 2002, Geometric Numerical Integration – Structure Preserving Algorithms for Ordinary Differential Equations, Springer Verlag, Berlin, Heidelberg.
[8] Hairer, E., Mairer, M.: 2002, http://www.unige.ch/math/folks
[9] Hairer, E., Lubich, C., Wanner, G.: Acta Numerica, 1
[10] Henrici, P: 1962, Discrete Variable Methods for Ordinary Differential Equations, J. Wiley, New York, London
[11] Higham, N.J.: 1993, SIAM J. Sci. Comput. 14, 787
[12] Huang, W., Leimkuhler, B.: 1997, SIAM J. Sci. Comput., 18, 239
[13] Quinn, T., Tremaine, S: 1990, Astron. J., 99, 1016
[14] Quinn, T., Tremaine, S., Duncan, M: 1991, Asstron. J., 101, 2287
[15] Sanz-Serna, J.M., Calvo, P.M.: 1994, Numerical Hamiltonian Problems, Chapman and Hall, London.
[16] Stuart, A.M., Humphries, A.R.: 1996, Dynamical systems and Numerical Analysis, Cambridge University Press. Cambridge, UK.
[17] Sussman, G.J., Wisdom, J: 1992, Science, 257, 56
[18] Wisdom, J., Holman, M.: Astron. J., 102, 1528
[19] Yosida, H.: 1990, Phys. Lett., A, 150, 262.
[20] Yosida, H: 1993, Celest Mech. Dyn. Astron., 56, 27
[21] Zhang, M., Skeel, R. D.: 1997, Appl. Num. Meth., 25, 297