Analytical orbit predictions with air drag using KS uniformly regular canonical elements

A new nonsingular analytical theory for the motion of near Earth satellite orbits with the air drag effect is developed for long term motion in terms of the KS uniformly regular canonical elements by a series expansion method, by assuming the atmosphere to be symmetrically spherical with constant density scale height. The series expansions include up to third order terms in eccentricity. Only two of the nine equations are solved analytically to compute the state vector and change in energy at the end of each revolution, due to symmetry in the equations of motion. Numerical comparisons of the important orbital parameters semi major axis and eccentricity up to 1000 revolutions, obtained with the present solution, with KS elements analytical solution and Cook, King-Hele and Walker's theory with respect to the numerically integrated values, show the superiority of the present solution over the other two theories over a wide range of eccentricity, perigee height and inclination.     

Orbit computations with the Eulerian topological regular elements with oblateness

In this paper, the Eulerian set of topological regular elements was utilized to develop an orbit computation package for the initial value problem for any conic motion of artificial satellites in the Earth's gravitational field with axial symmetry. Applications for the two types of short- and long-term predictions are considered. The numerical results proved the high efficiency and flexibility of the package.     

Contraction of satellite orbits using KS uniformly regular canonical elements in an oblate diurnally varying atmosphere

A new non-singular analytical theory for the motion of near Earth satellite orbits with the air drag effect is developed in terms of the Kustaanheimo and Stiefel (KS) uniformly regular canonical elements, by assuming the atmosphere to be oblate diurnally varying with constant density scale height. The series expansions include up to third-order terms in eccentricity and c (a small parameter dependent on the flattening of the atmosphere). Only two of the nine equations are solved analytically to compute the state vector and change in energy at the end of each revolution, due to symmetry in the equations of motion. Numerical comparisons of the important orbital parameters semimajor axis and eccentricity up to 1000 revolutions, obtained with the present solution, with the third-order analytical theories of Swinerd and Boulton and in terms of the KS elements, with respect to the numerically integrated values, show the superiority of the present solution over the other two theories over a wide range of eccentricity, perigee height and inclination.     

Contraction of near-Earth satellite orbits using uniformly regular KS canonical elements in an oblate atmosphere with density scale height variation with altitude

A new non-singular analytical theory for the contraction of near-Earth satellite orbits under the influence of air drag is developed in terms of uniformly regular Kustaanheimo and Stiefel (KS) canonical elements using an oblate atmosphere with variation of density scale height with altitude. The series expansions include up to fourth power in terms of eccentricity and c (a small parameter dependent on the flattening of the atmosphere). Only two of the nine equations are solved analytically to compute the state vector and change in energy at the end of each revolution, due to symmetry in the equations of motion. It is observed that the analytically computed values of the semi-major axis and eccentricity are consistent with the numerically integrated values up to 500 revolutions over a wide range of the drag-perturbed orbital parameters. The theory can be effectively used for re-entry of near-Earth objects.     

Analytical approach using KS elements to short-term orbit predictions including J 2

Analytical solutions using KS elements are derived. The perturbation considered is the Earth's zonal harmonic J ₂. The series expansions include terms of fourth power in the eccentricity. Only two of the nine KS element equations are integrated analytically due to the reasons of symmetry. The analytical solution is suitable for short-term orbit computations. Numerical studies show that reasonably good estimates of the orbital elements can be obtained in one step of 10 to 30 degrees of eccentric anomaly for near-Earth orbits of moderate eccentricity. For application purposes, the analytical solution can be effectively used for onboard computation in the navigation and guidance packages, where the modelling of J ₂ effect becomes necessary.     

Analytical short-term orbit predictions withJ 3 andJ 4 in terms of KS elements

Analytical theory for short-term orbit motion of satellite orbits with Earth's zonal harmonicsJ ₃ andJ ₄ is developed in terms of KS elements. Due to symmetry in KS element equations, only two of the nine equations are integrated analytically. The series expansions include terms of third power in the eccentricity. Numerical studies with two test cases reveal that orbital elements obtained from the analytical expressions match quite well with numerically integrated values during a revolution. Typically for an orbit with perigee height, eccentricity and inclination of 421.9 km, 0.17524 and 30 degrees, respectively, maximum differences of 27 and 25 cm in semimajor axis computation are noted withJ ₃ andJ ₄ term during a revolution. For application purposes, the analytical solutions can be used for accurate onboard computation of state vector in navigation and guidance packages.     

Prediction of trajectories in the Earth's gravitational field with axial symmetry using KS uniform regular elements

In this paper, a motion prediction algorithm based on the KS regular elements is developed for the motion in the Earth's gravitational field with axial symmetry. The algorithm is of recursive nature and general in the sense that it could be applied for any conic motion whatever the number of zonal harmonic coefficientsN 2 may be. Applications of the algorithm for the problem of the final state prediction are illustrated by numerical examples of eight typical ballistic missiles for geopotential model with zonal harmonic terms up toJ ₃₆. A final state of any desired accuracy is obtained for each case study, a result which shows the efficiency and the flexibility of the algorithm.     

Halo orbit of regularized circular restricted three-body problem with radiation pressure and oblateness

In this paper, computation of the halo orbit for the KS-regularized photogravitational circular restricted three-body problem is carried out. This work extends the idea of Srivastava et al. (Astrophys. Space Sci. 362: 49, 2017) which only concentrated on the (i) regularization of the 3D-governing equations of motion, and (ii) validation of the modeling for small out-of-plane amplitude (\(A_z =110000\) km) assuming the third-order analytical approximation as an initial guess with and without differential correction. This motivated us to compute the halo orbits for the large out-of-plane amplitudes and to study their stability analysis for the regularized motion. The stability indices are described as a function of out-of-plane amplitude, mass reduction factor and oblateness coefficient. Three different Sun–planet systems: the Sun–Earth, Sun–Mars and the Sun–Jupiter are chosen in this study. Stable halo orbits do not exist around the \(L_{1}\) point, however, around the \(L_{2}\) point stable halo orbits are found for the considered systems.     

Prediction of satellite orbits contraction due to diurnally varying oblate atmosphere and altitude-dependent scale height using KS canonical elements

A new non-singular analytical theory for the motion of near-Earth satellite orbits with the air drag effect is developed in terms of uniformly regular KS canonical elements. Diurnally varying oblate atmosphere is considered with variation in density scale height dependent on altitude. The series expansion method is utilized to generate the analytical solutions and terms up to fourth-order terms in eccentricity and c (a small parameter dependent on the flattening of the atmosphere) are retained. Only two of the nine equations are solved analytically to compute the state vector and change in energy at the end of each revolution, due to symmetry in the equations of motion. The important drag perturbed orbital parameters: semi-major axis and eccentricity are obtained up to 500 revolutions, with the present analytical theory and by numerical integration over a wide range of perigee height, eccentricity and inclination. The differences between the two are found to be very less. A comparison between the theories generated with terms up to third- and fourth-order terms in c and e shows an improvement in the computation of the orbital parameters semi-major axis and eccentricity, up to 9%. The theory can be effectively used for the re-entry of the near-Earth objects, which mainly decay due to atmospheric drag.