Author(s):
1. A. Amin:
Department of Mathematics, University of Engineering and Technology,Lahore,Pakistan
2. S. Rehman:
Department of Mathematics, University of Engineering and Technology,Lahore,Pakistan
3. A. Pervaiz:
Department of Mathematics, University of Engineering and Technology,Lahore,Pakistan
Abstract:
The N-body problem is one of the well-known and most central computational problem. The N-body problem of the Solar System is not only a rich source of initial value problems (IVPs) for ordinary differential equations (ODEs), but is also very convenient to understand the orbital evolution of the Solar System; see, for example,[1, 2]. Wide range of numerical integrators have been developed and implemented for performing such N-body simulations. The main objective of this research paper is to analyze and compare the accuracy and efficiency of different ordinary differential equation (ODE) solvers applied to the Kepler's two-body problem for Terrestrial planets. Throughout this paper, the error growth is investigated in terms of global error in position and velocity, and the relative error in terms of angular momentum and total energy of the system. To quantify the quality of different ODE solvers, we performed numerical tests applied to the Kepler's two body problem for Terrestrial planets with local error tolerances ranging from 10-12 to 10-4.
Page(s):
203-211
DOI:
DOI not available
Published:
Journal: Proceedings of the Pakistan Academy of Sciences: A. Physical and Computational Sciences, Volume: 53, Issue: 2, Year: 2016
Keywords:
Keplers two body problem
,
Terrestrial planets
,
ODE solvers
,
Nbody simulations
References:
References are not available for this document.
Citations
Citations are not available for this document.