On prediction of re-entry time of an upper stage from GTO

The evolution of objects in geostationary transfer orbit (GTO) is determined by a complex interplay of atmospheric drag and luni-solar gravity. These orbits are highly eccentric (eccentricity >0.7) and have large variations in velocity and perturbations during a revolution. The periodic changes in the perigee altitudes of these orbits are mainly due to the gravitational perturbations of the Sun and the Moon. The re-entry time of the objects in such orbits is sensitive to the initial conditions. The aim of this paper is to study the re-entry time of the cryogenic stage of the Indian geo-synchronous launch vehicle, GSLV-F04/CS, which has been decaying since 2 September 2007 from initial orbit with eccentricity equal to 0.706. Two parameters, initial eccentricity and ballistic coefficient, are chosen for optimal estimation. It is known that the errors are more in eccentricity for the observations based on two line elements (TLEs). These two parameters are computed with response surface method using a genetic algorithm for the selected eight different zones, based on rough linear variation of the mean apogee altitude during 200 days orbit evolution. The study shows that the GSLV-F04/CS will re-enter between 5 December 2010 and 7 January 2011. The methodology is also applied to study the re-entry of six decayed objects (cryogenic stages of GSLV and Molniya satellites). Good agreement is noticed between the actual and the predicted re-entry times. The absolute percentage error in re-entry prediction time for all the six objects is found to be less than 7%. The present methodology is being adopted at Vikram Sarabhai Space Centre (VSSC) to predict the re-entry time of GSLV-F04/CS.