Influence of Earth’s shadowing effects on space debris stability

Solar radiation pressure affects the evolution of high area-to-mass geostationary space debris. In this work, we extend the stability study of Valk et al. (2009) by considering the influence of Earth’s shadows on the short- and long-term time evolutions of space debris. To assess the orbits stability, we use the Mean Exponential Growth factor of Nearby Orbits (MEGNO), which is an efficient numerical tool to distinguish between regular and chaotic behaviors. To reliably compute long-term space debris motion, we resort to the Global Symplectic Integrator (GSI) of Libert et al. (2011) which consists in the symplectic integration of both Hamiltonian equations of motion and variational equations. We show how to efficiently compute the MEGNO indicator in a complete symplectic framework, and we also discuss the choice of a symplectic integrator, since propagators adapted to the structure of the Hamiltonian equations of motion are not necessarily suited for the associated variational equations. The performances of our method are illustrated and validated through the study of the Arnold diffusion problem. We then analyze the effects of Earth’s shadows, using the adapted conical and cylindrical Earth’s shadowing models introduced by Hubaux et al. (2012) as the smooth shadow function deriving from these models can be easily included into the variational equations. Our stability study shows that Earth’s shadows greatly affect the global behaviour of space debris orbits by increasing the size of chaotic regions around the geostationary altitude. We also emphasize the differences in the results given by conical or cylindrical Earth’s shadowing models. Finally, such results are compared with a non-symplectic integration scheme.     

Semi-analytical investigations of high area-to-mass ratio geosynchronous space debris including Earth’s shadowing effects

This paper investigates the long-term perturbations of the orbits of geosynchronous space debris influenced by direct radiation pressure including the Earth’s shadowing effects. For this purpose, we propose an extension of our homemade semi-analytical theory [Valk, S., Lemaître, A., Deleflie, F. Semi-analytical theory of mean orbital motion for geosynchronous space debris under gravitational influence. Adv. Space Res., submitted for publication], based on the method developed by Aksnes [Aksnes, K. Short-period and long-period perturbations of a spherical satellite due to direct solar radiation. Celest. Mech. Dyn. Astron. 13, 89–104, 1976] and generalized into a more convenient non-singular formalism. The perturbations accounting for the direct radiation pressure with the Earth’s shadow are computed on a revolution-by-revolution basis, retaining the original osculating Hamiltonian disturbing function. In this framework, we compute the non-singular mean longitude at shadow entry and shadow exit at every orbital revolution in opposition to classical approaches where the singular eccentric anomalies at shadow entry and shadow exit are computed. This new algorithm is developed using non-singular variables. Consequently, it is particularly suitable for both near-circular and near-equatorial orbits as well as orbits which transit periodically around null eccentricities and null inclinations.The algorithm is tested by means of numerical integrations of the equations, averaged over the short periods, including radiation pressure, J₂, the combined Moon and Sun third body attraction as well as the long-term effects of the 1:1 resonance occurring for geosynchronous objects. As an extension of [Valk, S., Lemaître, A., Anselmo, L. Analytical and semi-analytical investigations of geosynchronous space debris with high area-to-mass ratios influenced by solar radiation pressure. Adv. Space Res., doi:10.1016/j.asr.2007.10.025, 2007b], we especially apply our analysis to space debris with area-to-mass as high as 20 m²/kg. This paper provides numerical and semi-analytical investigations leading to a deep understanding of the long-term evolution of the semi-major axis. Finally, these semi-analytical investigations are compared with accurate numerical integrations of the osculating equations of motion over time scales as high as 25 years.     

Long-term predictability of orbits around the geosynchronous altitude

For a family of simulated geostationary, geosynchronous and super-geostationary orbits we compute the indicator called MEGNO in order to find out which orbits are chaotic and the timescale of their exponential divergence. A symplectic integrator of the Wisdom–Holman type was used for this purpose, with an integration span of 40 years. The results indicate that chaotic orbits exist only at the separatrix between geostationary libration and circulation and that they are relatively rare. The super-geostationary region seems to be entirely regular and quasi-periodic on the timescale of few decades.     

Influence of the 2:1 resonance in the orbits of GPS satellites

Due to the characteristics of their orbits the GPS satellites are submitted to the following main perturbations: terrestrial gravitational field, luni-solar gravitational attraction and solar radiation pressure (including the effects of the Earth's shadow). An additional perturbation arises due to the 2:1 commensurability of the orbital period of the satellite with the period of the Earth's rotation. An analytical theory is briefly presented to solve the equations of motion including the previously mentioned effects. The analytical solution, based on the Lie-Hori method, is compared with a numerical integration of the equations.     

Resonant motion, ballistic escape, and their applications in astrodynamics

A special set of solutions governing the motion of a particle, subject to the gravitational attractions of the Earth, the Moon, and, eventually, the Sun, is discussed in this paper. These solutions, called resonant orbits, correspond to a special motion where the particle is in resonance with the Moon. For a restricted set of initial conditions the particle performs a resonance transition in the vicinity of the Moon. In this paper, the nature of the resonance transition is investigated under the perspective of the dynamical system theory and the energy approach. In particular, using a new definition of weak stability boundary, we show that the resonance transition mechanism is strictly related to the concept of weak capture. This is shown through a carefully computed set of Poincaré surfaces, at different energy levels, on which both the weak stability boundary and the resonant orbits are represented. It is numerically demonstrated that resonance transitioning orbits pass through the weak stability boundaries. In the second part of the paper the solar perturbation is taken into account, and the motion of the resonant orbits is studied within a four-body dynamics. We show that, for a wide class of initial conditions, the particle escapes from the Earth–Moon system and targets an heliocentric orbit. This is a free ejection called a ballistic escape. Astrodynamical applications are discussed.     

基于GLONASS星历的预报轨道的误差分析

推导了协议地球坐标系下的卫星运动方程.通过分析由GLONASS(Global Navigation Satellite System)广播星历参数确定的卫星预报轨道的拟合精度,指出了摄动力模型的简化、积分器的选择,以及忽略了极移影响等因素是引起拟合误差的主要因素,其中摄动力模型的简化起最主要的作用.通过对卫星轨道运动方程积分30?min,可知由摄动力模型的简化、积分器和忽略极移影响等因素引起的拟合误差分别为0.827?m,0.224?m和0.025?m.要提高预报轨道拟合的精度,关键是要对摄动力简化特别是地球引力摄动高阶项的截断以及日月引力场简化造成的轨道预报精度损失加以控制.     

Two-dimensional phase plane structure and the stability of the orbital motion for space debris in the geosynchronous ring

Based on the orbital resonance model, we study the two-dimensional phase plane structure of the motion of space debris orbiting the geosynchronous ring under the combined effects of the tesseral harmonics J₂₂, J₃₁ and J₃₃ of the Earth’s gravitational field. We present the main characteristic parameters of the two-dimensional phase plane structure. We also analyze the stability of the two-dimensional phase plane structure with numerical method. Our main findings indicate that the combined effects of the tesseral harmonics J₂₂, J₃₁ and J₃₃ fully determine the two-dimensional phase plane structure of the space debris, and it remains robust under the effect of the Earth’s actual gravitational field, the luni-solar perturbations and the solar radiation pressure with the normal area-to-mass ratios.     

Impact of the time span selected to calibrate the ballistic parameter on spacecraft re-entry predictions

The trajectory modeling of satellites that are re-entering the Earth’s atmosphere, as a result of natural orbital decay, has always been a challenging task. Residual lifetime estimations and re-entry predictions are affected by substantial uncertainties, associated with atmospheric density models, with the forecasts of the relevant solar and geomagnetic activity indices and with tracking data, which for uncontrolled re-entries are usually sparse and not particularly accurate. Furthermore, modeling the aerodynamic forces that act on low altitude satellites is a formidable task, especially for objects of a complex shape and unknown attitude evolution.     

Dynamics and control of spacecraft hovering using the geomagnetic Lorentz force

To achieve hovering, a spacecraft thrusts continuously to induce an equilibrium state at a desired position. Due to the constraints on the quantity of propellant onboard, long-time hovering around low-Earth orbits (LEO) is hardly achievable using traditional chemical propulsion. The Lorentz force, acting on an electrostatically charged spacecraft as it moves through a planetary magnetic field, provides a new propellantless method for orbital maneuvers. This paper investigates the feasibility of using the induced Lorentz force as an auxiliary means of propulsion for spacecraft hovering. Assuming that the Earth’s magnetic field is a dipole that rotates with the Earth, a dynamical model that characterizes the relative motion of Lorentz spacecraft is derived to analyze the required open-loop control acceleration for hovering. Based on this dynamical model, we first present the hovering configurations that could achieve propellantless hovering and the corresponding required specific charge of a Lorentz spacecraft. For other configurations, optimal open-loop control laws that minimize the control energy consumption are designed. Likewise, the optimal trajectories of required specific charge and control acceleration are both presented. The effect of orbital inclination on the expenditure of control energy is also analyzed. Further, we also develop a closed-loop control approach for propellantless hovering. Numerical results prove the validity of proposed control methods for hovering and show that hovering around low-Earth orbits would be achievable if the required specific charge of a Lorentz spacecraft becomes feasible in the future. Typically, hovering radially several kilometers above a target in LEO requires specific charges on the order of 0.1 C/kg.     

双小行星系统引力场建模方法研究

针对弱引力双小行星系统的引力场建模问题,本文采用复杂度和精度依次递增的球体–球体模型、椭球体–球体模型和改进的限制性椭球体–椭球体模型来进行引力场建模,并分别采用椭圆积分以及无积分环节、计算效率高的二阶二次球谐函数来表征引力势,从而比较精确地刻画双小行星系统和探测器构成的限制性全三体问题的动力学模型;针对双小行星系统1999KW4,对其不同的引力场模型进行了仿真研究,分别给出了不同模型下的等效势能函数曲面及零速度曲线,比较了不同模型下的平动点位置坐标偏差。结果表明,二阶二次球谐函数计算引力势的椭球体-椭球体模型计算精度高,复杂程度低,计算量更少,计算速度更快,能够较精确的对双小行星系统进行引力场建模。     

Low-thrust transfer to Backflip orbits

The aim of the work is to design a low-thrust transfer from a Low Earth Orbit to a “useful” periodic orbit in the Earth–Moon Circular Restricted Three Body Model (CR3BP). A useful periodic orbit is here intended as one that moves both in the Earth–Moon plane and out of this plane without any requirements of propellant mass. This is achieved by exploiting a particular class of periodic orbits named Backflip orbits, enabled by the CR3BP. The unique characteristics of this class of periodic solutions allow the design of an almost planar transfer from a geocentric orbit and the use of the Backflip intrinsic characteristics to explore the geospace out of the Earth–Moon plane. The main advantage of this approach is that periodic plane changes can be obtained by performing an almost planar transfer. In order to save propellant mass, so as to increase the scientific payload of the mission, a low-powered transfer is considered. This foresees a thrusting phase to gain energy from a departing circular geocentric orbit and a second thrusting phase to match the state of the target Backflip orbit, separated by an intermediate ballistic phase. This results in a combined application of a low-thrust manoeuvre and of a periodical solution in the CR3BP to realize a new class of missions to explore the Earth–Moon neighbourhoods in a quite inexpensive way. In addition, a low-thrust transit between two different Backflip orbits is analyzed and considered as a possible extension of the proposed mission. Thus, also a Backflip-to-Backflip transfer is addressed where a low-powered probe is able to experience periodic excursions above and below the Earth–Moon plane only performing almost planar and very short transfers.     

A high order method for orbital conjunctions analysis: Sensitivity to initial uncertainties

A high order method to quickly assess the effect that uncertainties produce on orbital conjunctions through a numerical high-fidelity propagator is presented. In particular, the dependency of time and distance of closest approach to initial uncertainties on position and velocity of both objects involved in a conjunction is studied. The approach relies on a numerical integration based on differential algebraic techniques and a high-order algorithm that expands the time and distance of closest approach in Taylor series with respect to relevant uncertainties. The modeled perturbations are atmospheric drag, using NRLMSISE-00 air density model, solar radiation pressure with shadow, third body perturbation using JPL’s DE405 ephemeris, and EGM2008 gravity model. The polynomial approximation of the final position is used as an input to compute analytically the expansion of time and distance of closest approach. As a result, the analysis of a close encounter can be performed through fast, multiple evaluations of Taylor polynomials. Test cases with objects ranging from LEO to GEO regimes are considered to assess the performances and the accuracy of the proposed method.     

Properties of the high area-to-mass ratio space debris population at high altitudes

In the framework of its space debris research activities ESA established an optical survey program to study the space debris environment at high altitudes, in particular in the geostationary ring and in the geostationary transfer orbit region. The Astronomical Institute of the University of Bern (AIUB) performs these surveys on behalf of ESA using ESA’s 1-m telescope in Tenerife. Regular observations were started in 1999 and are continued during about 120–140 nights per year. Results from these surveys revealed a substantial amount of space debris at high altitudes in the size range from 0.1 to 1 m. Several space debris populations with different dynamical properties were identified in the geostationary ring. During the searches for debris in the geostationary transfer orbit region a new population of objects in unexpected orbits, where no potential progenitors exist, was found. The orbital periods of these objects are clustered around one revolution per day; the eccentricities, however, are scattered between 0 and 0.6. By following-up some of these objects using the ESA telescope and AIUB’s 1-m telescope in Zimmerwald, Switzerland, it was possible to study the properties of this new population. One spectacular finding from monitoring the orbits over time spans of days to months is the fact that these objects must have extreme area-to-mass ratios, which are by several orders of magnitudes higher than for ‘normal-type’ debris. This in turn supports the hypothesis that the new population actually is debris generated in or near the geostationary ring and which is in orbits with periodically varying eccentricity and inclination due to perturbations by solar radiation pressure. In order to further study the nature of these debris, multi-color and temporal photometry (light curves) were acquired with the Zimmerwald telescope. The light curves show strong variations over short time intervals, including signals typical for specular reflections. Some objects exhibit distinct periodic variations with periods ranging from 10 to several 100 s. All this is indicative for objects with complicated shapes and some highly reflective surfaces.     

Study of waves in the magnetotail region with cluster and DSP

The study of the neutral sheet is of fundamental importance in understanding the dynamics of the Earth’s magnetosphere. From the earliest observation of the magnetotail, it has been found that the neutral sheet frequently appears to be in motion due to changing solar wind conditions and geomagnetic activity. Multiple crossings of the neutral sheet by spacecraft have been attributed to a flapping motion of the neutral sheet in the north–south direction, a wavy profile either along the magnetotail or the dawn–dusk direction. Cluster observations have revealed that the flapping motions of the Earth’s magnetotail are of internal origin and that kink-like waves are emitted from the central part of the tail and propagate toward the tail flanks. This flapping motion is shown here to propagate at an angle of ∼45° with x_GSM. A possible assumption that the flapping could be created by a wake travelling away from a fast flow in the current sheet is rejected. Other waves in the magnetotail are found in the ULF range. One conjunction event between Cluster and DoubleStar TC1 is presented where all spacecraft show ULF wave activity at a period of approximately 5 min during fast Earthward flow. These waves are shown to be Kelvin–Helmholtz waves on the boundaries of the flow channel. Calculations show that the conversion of flow energy into magnetic energy through the Kelvin–Helmholtz instability can contribute to a significant part of flow breaking between Cluster and DoubleStar TC1.     

Evolution of the orbital elements for objects with high area-to-mass ratios in geostationary transfer orbits

A new population of uncatalogued objects in geosynchronous Earth orbits (GEO), with a mean motion of about 1 rev/day and eccentricities up to 0.6, has been identified recently. The first observations of this new type of objects were acquired in the framework of the European Space Agency’s (ESA) search for space debris in GEO and the geostationary transfer orbit (GTO) using the ESA 1-m telescope on Tenerife. Earlier studies have postulated that the perturbations due to the solar radiation pressure can lead to such large eccentricities for GEO objects with a high area-to-mass ratio (A/M). The simulations showed that the eccentricities of GEO objects with large A/M exhibit periodic variations with periods of about one year and amplitudes depending on the value of A/M. The findings of these studies could be confirmed by observations from the ESA 1-m telescope on Tenerife.     

Heliotropic dust rings for Earth climate engineering

This paper examines the concept of a Sun-pointing elliptical Earth ring comprised of dust grains to offset global warming. A new family of non-Keplerian periodic orbits, under the effects of solar radiation pressure and the Earth’s J₂ oblateness perturbation, is used to increase the lifetime of the passive cloud of particles and, thus, increase the efficiency of this geoengineering strategy. An analytical model is used to predict the orbit evolution of the dust ring due to solar-radiation pressure and the J₂ effect. The attenuation of the solar radiation can then be calculated from the ring model. In comparison to circular orbits, eccentric orbits yield a more stable environment for small grain sizes and therefore achieve higher efficiencies when the orbit decay of the material is considered. Moreover, the novel orbital dynamics experienced by high area-to-mass ratio objects, influenced by solar radiation pressure and the J₂ effect, ensure the ring will maintain a permanent heliotropic shape, with dust spending the largest portion of time on the Sun facing side of the orbit. It is envisaged that small dust grains can be released from a circular generator orbit with an initial impulse to enter an eccentric orbit with Sun-facing apogee. Finally, a lowest estimate of 1 × 10¹² kg of material is computed as the total mass required to offset the effects of global warming.     

On the stabilizing effect of Solar Radiation Pressure in the Earth-Moon system

Solar sails change the natural dynamics of systems: Trajectories that are driven by gravitational forces can be displaced and changed because of the effect of Solar Radiation Pressure (SRP). Moreover, if the lightness number of the sail is large enough, the instability of certain orbits can be diminished and even removed. In this paper we modify two models for the motion of a probe in the Earth-Moon system that include the effect of Sun’s gravity to take also into account the effect of SRP. These models, the Bicircular Problem (BCP) and the Quasi-Bicircular Problem (QBCP), are periodic perturbations of the Earth-Moon Restricted Three Body Problem (RTBP). The models are modified to consider the effect of the SRP upon a solar sail. We provide examples of periodic orbits that are stabilized (or made less unstable) due to the effect of SRP.     

Multiple Sun–Earth Saddle Point flybys for LISA Pathfinder

LISA Pathfinder is an ESA mission due to be launched in the next two years. The gravity gradiometer onboard has the sensitivity required to test predictions by gravitational theories proposed as alternatives to Dark Matter such as TeVeS. Within the Solar System measurable effects are predicted only in the vicinity of gravitational saddle points (SP). For this reason it has been proposed to fly LPF by the Earth–Sun SP, at some 259,000 km from Earth. This could be done in an extension to the nominal mission which uses a Lissajous orbit about the Earth–Sun L1 point. The responsibility for LPF mission design lies with ESA/ESOC, who have designed the transfer trajectories, orbits about L1, and station keeping strategies. This article describes an analysis performed by Astrium to support a suggestion for a possible mission extension to a saddle point crossing. With only very limited fuel availability, reaching the saddle point is a significant challenge. In this article, we present recent advances in the work on trajectory design. It is demonstrated that reaching the SP is feasible once the LPF mission is completed. Furthermore, in a significant enhancement, it is demonstrated that trajectories including more than one SP flyby are possible, thus improving the science return for this proposed mission extension.     

小行星远距离抵近轨道

讨论了小行星引力一阶项可被忽略情况下的小行星远距离轨道设计及动力学。此时,航天器的运动受太阳引力和太阳光压的影响。航天器和小行星的加速度之差在这两者之间形成的独特的相对动力学,为航天器在小行星附近停驻与观测提供特定轨道。完整解决了小行星处于圆形日心轨道这一较简单情况,也考虑和阐述了椭圆轨道情况,并取得了一些初步结果。