Controlling the Wrapping Effect in the Solution of ODEs for Asteroids

During the last decade, substantial progress has been made in fighting the wrapping effect in self-validated integrations of linear systems. However, it is still the main problem limiting the applicability of such methods to the long-term integration of non-linear systems. Here we show how high-order self-validated methods can successfully overcome this obstacle.We study and compare the validated integration of a Kepler problem with conventional and high-order methods represented by AWA and Taylor models, respectively. We show that this simple model problem exhibits significant wrapping that is particularly difficult to control for conventional first-order methods. It will become clear that utilizing high-order methods with shrink wrapping allows the system to be analyzed in a fully validated context over large integration times. By comparing high-order Taylor model integrations with Taylor model methods subjected to an artificial wrapping effect, we show that utilizing high-order methods to propagate initial conditions is indeed the foremost reason for the successful suppression of the wrapping effect.To further demonstrate that high-order Taylor model methods can be used for the integration of complicated non-linear systems, we summarize results obtained from a fully verified and self-validated orbit integration of the near earth asteroid 1997 XF11. Since this asteroid will have several close encounters with Earth, its analysis is an important application of reliable computations.

[1]  Martin Berz,et al.  5. Remainder Differential Algebras and Their Applications , 1996 .

[2]  Josef Honerkamp,et al.  Klassische Theoretische Physik , 1989 .

[3]  Clifford M. Will,et al.  The theoretical tools of experimental gravitation , 1972 .

[4]  Martin Berz,et al.  Efficient Control of the Dependency Problem Based on Taylor Model Methods , 1999, Reliab. Comput..

[5]  Tullio Levi-Civita.  The n-Body Problem in General Relativity , 1965 .

[6]  Kyoko Makino,et al.  Rigorous analysis of nonlinear motion in particle accelerators , 1998 .

[7]  G. Verschuur.  Impact!: The Threat of Comets and Asteroids , 1996 .

[8]  J. W. Humberston.  Classical mechanics , 1980, Nature.

[9]  Donald K. Yeomans.  Comet and asteroid ephemerides for spacecraft encounters , 1997 .

[10]  Martin Berz,et al.  Computation and Application of Taylor Polynomials with Interval Remainder Bounds , 1998, Reliab. Comput..

[11]  Martin Berz,et al.  Verified Integration of ODEs and Flows Using Differential Algebraic Methods on High-Order Taylor Models , 1998, Reliab. Comput..

[12]  P. Seidelmann.  Explanatory Supplement to the Astronomical Almanac , 2005 .

[13]  William McCrea,et al.  The Theory of Space, Time and Gravitation , 1961 .

[14]  N. Spyrou.  The N-body problem in general relativity , 1975 .

[15]  Wolfgang Kuehn,et al.  Towards an optimal control of the wrapping effect , 1998, SCAN.

[16]  Ramon E. Moore.  Interval arithmetic and automatic error analysis in digital computing , 1963 .

[17]  N. Nedialkov,et al.  Computing rigorous bounds on the solution of an initial value problem for an ordinary differential equation , 1999 .