On the Dominant Lunisolar Perturbations for Long-Term Eccentricity Variation: The Case of Molniya Satellite Orbits

The aim of this work is to investigate the main dominant terms of lunisolar perturbations that affect the orbital eccentricity of a Molniya satellite in the long term. From a practical point of view, these variations are important in the context of space situational awareness—for instance, to model the long-term evolution of artificial debris in a highly elliptical orbit or to design a reentry end-of-life strategy for a satellite in a highly elliptical orbit. The study assumes a doubly averaged model including the Earth’s oblateness effect and the lunisolar perturbations up to the third-order expansion. The work presents three important novelties with respect to the literature. First, the perturbing terms are ranked according to their amplitudes and periods. Second, the perturbing bodies are not assumed to move on circular orbits. Third, the lunisolar effect on the precession of the argument of pericenter is analyzed and discussed. As an example of theoretical a application, we depict the phase space description associated with each dominant term, taken as isolated, and we show which terms can apply to the relevant dynamics in the same region.

Semi-Analytical Estimates for the Orbital Stability of Earth’s Satellites

Normal form stability estimates are a basic tool of Celestial Mechanics for characterizing the long-term stability of the orbits of natural and artificial bodies. Using high-order normal form constructions, we provide three different estimates for the orbital stability of point-mass satellites orbiting around the Earth. (i) We demonstrate the long-term stability of the semimajor axis within the framework of the $$J_2$$ J 2 problem, by a normal form construction eliminating the fast angle in the corresponding Hamiltonian and obtaining $${\mathcal {H}}_{J_2}$$ H J 2 . (ii) We demonstrate the stability of the eccentricity and inclination in a secular Hamiltonian model including lunisolar perturbations (the ‘geolunisolar’ Hamiltonian $${\mathcal {H}}_\mathrm{gls}$$ H gls ), after a suitable reduction of the Hamiltonian to the Laplace plane. (iii) We numerically examine the convexity and steepness properties of the integrable part of the secular Hamiltonian in both the $${\mathcal {H}}_{J_2}$$ H J 2 and $${\mathcal {H}}_\mathrm{gls}$$ H gls models, which reflect necessary conditions for the holding of Nekhoroshev’s theorem on the exponential stability of the orbits. We find that the $${\mathcal {H}}_{J_2}$$ H J 2 model is non-convex, but satisfies a ‘three-jet’ condition, while the $${\mathcal {H}}_\mathrm{gls}$$ H gls model restores quasi-convexity by adding lunisolar terms in the Hamiltonian’s integrable part.

Impact of orbital orientations and radii on TianQin constellation stability

TianQin is a proposed space-based gravitational-wave observatory mission to be deployed in high circular Earth orbits. The equilateral-triangle constellation, with a nearly fixed orientation, can be distorted primarily under the lunisolar perturbations. To accommodate science payload requirements, one must optimize the orbits to stabilize the configuration in terms of arm-length, relative velocity, and breathing angle variations. In this paper, we present an efficient optimization method and investigate how changing the two main design factors, i.e. the orbital orientation and radius, impacts the constellation stability through single-variable studies. Thereby, one can arrive at the ranges of the orbital parameters that are comparatively more stable, which may assist future refined orbit design.

Searching for Orbits with Minimum Fuel Consumption for Station-Keeping Maneuvers: An Application to Lunisolar Perturbations

The present paper has the goal of developing a new criterion to search for orbits that minimize the fuel consumption for station-keeping maneuvers. This approach is based on the integral over the time of the perturbing forces. This integral measures the total variation of velocity caused by the perturbations in the spacecraft, which corresponds to the equivalent variation of velocity that an engine should deliver to the spacecraft to compensate the perturbations and to keep its orbit Keplerian all the time. This integral is a characteristic of the orbit and the set of perturbations considered and does not depend on the type of engine used. In this sense, this integral can be seen as a criterion to select the orbit of the spacecraft. When this value becomes larger, more consumption of fuel is required for the station keeping, and, in this sense, less interesting is the orbit. This concept can be applied to any perturbation. In the present research, as an example, the perturbation caused by a third body is considered. Then, numerical simulations considering the effects of the Sun and the Moon in a satellite around the Earth are shown to exemplify the method.

Semianalytic Integration of High-Altitude Orbits under Lunisolar Effects

The long-term effect of lunisolar perturbations on high-altitude orbits is studied after a double averaging procedure that removes both the mean anomaly of the satellite and that of the moon. Lunisolar effects acting on high-altitude orbits are comparable in magnitude to the Earth’s oblateness perturbation. Hence, their accurate modeling does not allow for the usual truncation of the expansion of the third-body disturbing function up to the second degree. Using canonical perturbation theory, the averaging is carried out up to the order where second-order terms in the Earth oblateness coefficient are apparent. This truncation order forces to take into account up to the fifth degree in the expansion of the lunar disturbing function. The small values of the moon’s orbital eccentricity and inclination with respect to the ecliptic allow for some simplification. Nevertheless, as far as the averaging is carried out in closed form of the satellite’s orbit eccentricity, it is not restricted to low-eccentricity orbits.