^{1}

^{1}

^{*}

This paper deals with the re-entry predictions of the space objects from the low eccentric orbit. Any re-entering object re-enters the Earth’s atmosphere with a high orbital velocity. Due to the aerodynamic heating the object tends to break into multiple fragments which later pose a great risk hazard to the population. Here a satellite is considered as the space object for which the re-entry prediction is made. This prediction is made with a package where the trajectory path, the time of re-entry and the survival rate of the fragments is done. The prediction is done using DRAMA 2.0—ESA’s Debris Risk Assessment and Mitigation Analysis Tool suite, MATLAB and Numerical Prediction of Orbital Events software. The predicted re-entry time of OSIRIS 3U was found to be on 7th March 2019, 7:25 (UTC), whereas the actual re-entry time was on 7th March 2019, 7:03 (UTC). The trajectory path found was 51.5699 deg. (Lat), −86.5738 deg. (Long.) with an altitude of 168.643 km. But the actual trajectory was 51.76 deg. (Lat), −89.01deg. (Long.) with an altitude of 143.5 km.

Space objects refer to astronomical objects as well as the artificial space objects, i.e. naturally occurring or man-made objects in space. Both the astronomical and artificial objects tend to enter the Earth’s atmosphere, especially the man-made objects. Man-made objects like satellites and rockets are being sent to space for communication, navigation and other space missions. They tend to decay after a certain orbital lifetime. Spent rocket stages, old satellites and fragments from disintegration, erosion and collisions are considered as debris. In particular, an accurate estimation of orbital decay of objects during the final stages of re-entry is of considerable importance. This helps to predict the re-entry time and location and thus plan proper hazard assessment and mitigation strategies. The database available for the prediction of orbital lifetime and re-entry of debris objects is the set of two line elements (TLEs). In general, the physical parameters of the objects like mass, area of cross section, shape and dimensions are not available accurately. Further, the atmosphere in which the objects decay varies significantly. A low Earth orbit lies between the altitudes 150 and 2000 kilometer, with a period of about 88 minutes to 127 minutes. Some of the important studies carried out in the area reported in this study are in References [

In this paper, a method to carry out the re-entry predictions of space debris entering from low eccentricity Earth orbit, employing the orbital data in the form of TLEs, is presented. The ballistic coefficient and eccentricity of the re-entering objects are considered as uncertain parameters. The Earth’s zonal harmonic terms J_{2} to J_{6} are included along with the drag perturbation [

The re-entry time estimation in each case is computed using NPOE software and MATLAB. The influence of luni-solar perturbations, Earth’s oblateness and atmospheric drag are considered to predict the re-entry time. The re-entry trajectory path along with the debris risk analysis is determined using Debris Risk Assessment and Mitigation Analysis (DRAMA) software. The orbital data and other details of the space objects which re-enter are taken from space-track.org maintained by US Air Force.

NPOE is an interactive computer program for computers which can model important orbital events and predict the long-term evolution of satellites in Earth orbits. Program NPOE implements a special perturbation solution of orbital motion using a variable step size Runge-Kutta-Fehlberg (RKF78) integration method to numerically integrate Cowell’s form of the system of differential equations. Orbital events are predicted using Brent’s method for finding the root of a nonlinear equation.

MATLAB has also been used to predict the long-term behaviour of the Earth’s satellites subjected to various perturbations. The program actualizes a special perturbation solution of orbital motion as same as NPOE using the variable step size Runge-Kutta-Fehlberg (RKF78) integration method to numerically solve Cowell’s form of the system of differential equation subjected to the central body gravity and other external forces which is otherwise called as orbital initial value problem (IVP).MATLAB is used to plot the graphs of the orbital events.

DRAMA is a comprehensive tool [

The following tools that are available within DRAMA are:

1) ARES: Assessment of Risk Event Statistics

The ARES provides an assessment of collision-related events between an operational spacecraft and trackable objects orbiting the Earth, the statistical collision probability, the mean number of conjunction avoidance manoeuvres and the fuel consumption associated to those manoeuvres.

2) MIDAS: MASTER (Based) Impact Flux and Damage Assessment

In MIDAS user-defined BLEs can be provided and flux computations are now performed using the MASTER-2009 model. Debris and meteoroid collision flux and damage analysis are done using MIDAS.

3) OSCAR: Orbital Spacecraft Active Removal

OSCAR offers the possibility to select between different standardized methods (ISO, ECSS) to generate forecasts of future solar and geomagnetic activity. OSCAR allows for the simulation of drag augmentation devices as a new disposal system.

4) CROC: Cross Section of Complex Bodies

Computes the cross-section of complex bodies for different aspect angle conditions.

5) SARA: (Re-Entry) Survival and Risk Analysis

It is used for the simulation of the re-entry into the Earth’s atmosphere.

This prediction package is used to study the re-entry of OSIRIS 3U satellite. OSIRIS-3U is a CubeSat mission launched on 14 August 2017. Orbital Satellite for Investigating the Response of the Ionosphere to Stimulation and Space Weather is the acronym of OSIRIS which is a three-unit CubeSat developed by the students of the Penn State University. The mission of OSIRIS-3U is to investigate the radio wave interaction in the ionosphere, particularly the interaction of high-power radio waves.

Newton’s law [

F = − G M m r 2

where,

G = 6.673 × 10 − 11 m 3 ⋅ kg − 1 ⋅ s − 2 , the gravitational constant;

M = the mass of centre object;

m = the mass of orbiting object;

r = the distance between the two masses.

Transforming into the vector form,

r ¨ = − G M r 3 r

The mass of the satellite (m) is neglected due to its small mass.

In this study eccentricity and ballistic coefficient are considered as uncertain parameters [^{2} [

As mentioned in the introduction, NPOE and MATLAB script was used to implement a special perturbation solution of orbital motion using a variable step size Runge-Kutta-Fehlberg (RKF78) integration method to numerically solve Cowell’s form of the system of differential equation subject to the central body gravity and other external forces.

a ( r , v , t ) = r ¨ ( r , r ˙ , t ) = a g ( r ) + a d ( r , v , t ) + a s m ( r , t ) + ( r , t )

The satellite’s acceleration due to Earth’s gravity field is calculated with a vector equation derived from the gradient of the potential function expressed as

a g ( r , t ) = ∇ ϕ ( r , t )

This acceleration vector is a combination of pure two-body or point mass gravity acceleration and the gravitational acceleration due to higher-order non-spherical terms in the Earths geo-potential.

The acceleration experienced by the satellite due to atmospheric drag [

a d ( r , v , t ) = − 1 2 ρ ( r , t ) | v r | v r C d A m

The acceleration contribution of the Sun and Moon represented by point masses is given by

a s m ( r , t ) = − μ m ( r m − b | r m − b | 3 + r e − m | r e − m | 3 ) − μ s ( r s − b | r s − b | 3 + r e − s | r e − s | 3 )

Generally, in the low earth orbit, the space objects are more prone to the gravity and atmospheric drag. But when the altitude exceeds approximately 1000 km, the solar effects from the Sun increase. Here, OSIRIS 3U is at an altitude between 180 km - 210 km and so the solar radiation pressure along with the luni-solar perturbation is being neglected.

In this case, the variation of re-entry time is determined based on varying the ballistic coefficient (BC) [

It is found based on the data plotted (Figures 1-4) that the ballistic coefficient does not show much variation in the re-entry time as expected. The graphs show that the orbital lifetime with BC = 60 kg/m^{2} and BC = 80 kg/m^{2} show very less variation.

Hence, it was neglected as being an uncertain parameter and BC = 70 kg/m^{2} was taken for the prediction of re-entry time of OSIRIS 3U. The observed semi-major axis and the osculating eccentricity from the TLEs are given in

TLE Epoch (UTC) | Semi Major Axis (km) | Eccentricity |
---|---|---|

04.02.19 20:31 | 6669.964476 | 0.0005899031 |

09.02.19 13:25 | 6662.310370 | 0.0002481707 |

10.02.19 21:01 | 6664.590117 | 0.0004735118 |

11.02.19 22:34 | 6661.128219 | 0.0007581420 |

13.02.19 21:10 | 6656.298443 | 0.0011430093 |

15.02.19 19:43 | 6653.712256, | 0.0011598911 |

18.02.19 19:44 | 6641.503586 | 0.0012821558 |

22.02.19 12:08 | 6649.046280 | 0.0011109638 |

The osculating eccentricity and ballistic coefficient taken as uncertain parameters were found to be 0.0080354104 and 70 kg/m^{2}, respectively. Using the estimated values, the re-entry time was predicted to be on 7th March 7:25 (UTC) as shown in the graphs (

The Annual Collision Probability (ACP) [

c = v D A c Δ t

When simulation was performed on OSIRIS 3U with the ARES tool for the ACP functionality, the functionality returned with the value as in

Functionalities | Values |
---|---|

ACP_d | 0.7141E−05 |

ACP_w | 0.6424E−03 |

Flux_d | 0.6910E−01 |

Flux_w | 0.1011E+02 |

where, ACP_w—probability of collision with any object of the whole population; ACP_d—probability of collision with detectable objects; Flux_d—flux due to the detected population [1/km^{2}/year]; Flux_w—flux due to the whole population [1/km^{2}/year].

Mean number of Manoeuvres per Year

The probability of collision for an object-to object encounter is computed by

P = 1 2 π d e t ( C ) ∫ − R R ∫ − R 2 − x 2 R 2 − x 2 e − 1 2 δ r T C − 1 δ r d y d x

where,

R—sum of the two object radii;

δ r —vector between a point in the integration area and the point where the near-miss is predicted.

Risk Reduction and False Alarm Rate

Finally,

The required Δ V j (for each population group) is obtained by integration of ∆V F j over the manoeuvring area. The total ∆V ( Δ V T ) is obtained by adding the contribution of each population group:

Δ V T = ∑ j Δ V j

It is observed that when sooner the manoeuvres are performed the more the cheaper they are. However, manoeuvres that performed much before the close approach have large uncertainties associated. The less the time to an event, the better the uncertainties are known, and therefore, the better the prediction.

The propellant mass fraction burned during the expected avoidance manoeuvres, during the satellite lifetime, is linked to the required ∆V to perform those manoeuvres and the propulsion system characteristics.

The computation of the propellant mass fraction requires the specific impulse ( I s p ). The ratio between the propellant mass ( m p ) to be burned by a known ∆V and the initial satellite mass ( m o ) is given by:

P M F a m = m p m o = [ 1 − e − Δ V I s p g ]

For the given simulation time ∆t and the impact flux F, the number of impacts ( n i m p ) and the probability of collision ( P c o l l ) was computed.

Number of Impacts (NOI)— n i m p = F ⋅ A ⋅ Δ t

Probability of collision (POC)— P c o l l = 1 − e − n i m p

The orbit of OSIRIS 3U is the target orbit and is defined by the following parameters:

Semi-major Axis (a) = 6646.4 km.

Eccentricity (e) = 8.39E-3.

Inclination (i) = 51.6305 deg.

Right ascension of the ascending node (Ω) = 162.247 deg.

Argument of Perigee (ω) = 350.489 deg.

And the objects considered are of range:

0.1000E+00 m—Lower Threshold

0.1000E+03 m—Upper Threshold

The flux was considered for a diameter range 1 mm < d < 20 cm for a time span between May 1st, 2001 and May 1st, 2050. All sources of debris and the meteoroid were also considered.

m p > m

where,

mp = particle mass.

d p > d

where,

dp = particle diameter.

On comparing the overall number of impacts on the reference sphere with those on the Sun oriented surface there are more impacts on the oriented surface than on the sphere.

The residual lifetime of OSIRIS 3U was estimated by OSCAR for the minimum cross-section in flight direction. From

The cross-sectional [

The cross section of a satellite with respect to the aspect angle θ = 0 , φ = 0 was found to be 19,900.0 mm^{2} and the perturbation of atmospheric drag acts on this surface.

For an angle of θ = 45 , φ = 15 the surface area is 45411 mm^{2}.

Due to the large orbital velocity, the spacecraft experiences high mechanical loads and heating rates while entering the atmosphere and during the aero-braking process the potential and kinetic energy of the spacecraft is converted into thermal energy that is consumed by the spacecraft, resulting in thermal and mechanical loads destroying the spacecraft either completely or partially [

Object-oriented method was used to analyse and the breakup altitude was set for a range between 75 km and 85 km. The solar panel breakup altitude was assumed to be at 95 km.

Along with the re-entry trajectory dynamics, the aerodynamic, aero thermodynamic and thermal analyses were also performed for OSISRIS 3U re-entry. The initial conditions of the trajectory [

The trajectory path was found to be 51.5699 deg. (Lat), −86.5738 deg. (long.) with an altitude of 168.643 km at a velocity of 7.51313 km/s with an approximate temperature of 300 K.

It was found that no objects survived upon re-entry and so there was no ground impact.

The predicted re-entry time was found to be on 7th March 2019, 7:25 (UTC), whereas the actual re-entry time was on 7th March 2019, 7:03 (UTC). The trajectory path found was 51.5699 deg. (Lat), −86.5738 deg. (Long.) with an altitude of 168.643 km. But the actual trajectory was 51.76 deg. (Lat), −89.01 deg. (Long.) with an altitude of 143.5 km.

The authors declare no conflicts of interest regarding the publication of this paper.

Lawrence, C.S. and Sharma, R.K. (2019) Re-Entry of Space Objects from Low Eccentricity Orbits. International Journal of Astronomy and Astrophysics, 9, 200-216. https://doi.org/10.4236/ijaa.2019.93015