ISSN:
1572-9478
Source:
Springer Online Journal Archives 1860-2000
Topics:
Physics
Notes:
Abstract The main limit to the time span of a numerical integration of the planetary orbits is no longer set by the availability of computer resources, but rather by the accumulation of the integration error. By the latter we mean the difference between the computed orbit and the dynamical behaviour of the real physical system, whatever the causes. The analysis of these causes requires an interdisciplinary effort: there are physical model and parameters errors, algorithm and discretisation errors, rounding off errors and reliability problems in the computer hardware and system software, as well as instabilities in the dynamical system. We list all the sources of integration error we are aware of and discuss their relevance in determining the present limit to the time span of a meaningful integration of the orbit of the planets. At present this limit is of the order of 108 years for the outer planets. We discuss in more detail the truncation error of multistep algorithms (when applied to eccentric orbits), the coefficient error, the method of Encke and the associated coordinate change error, the procedures used to test the numerical integration software and their limitations. Many problems remain open, including the one of a realistic statistical model of the rounding off error; at present, the latter can only be described by a semiempirical model based upon the simpleN 2 ∈ formula (N=number of steps, ∈=machine accuracy), with an unknown numerical coefficient which is determined only a posteriori.
Type of Medium:
Electronic Resource
URL:
http://dx.doi.org/10.1007/BF01234550
Permalink