Dynamical evolution of high area-to-mass ratio objects in Molniya orbits

For special demands, some notable orbit types have been developed by human, including the Molniya orbits, which have a relatively high eccentricity up to about 0.7, and a period of 12 h. Considering that space debris with high area-to-mass ratio (A/M) has been discovered, such objects may also exist in Molniya orbits due to spacecraft and upper stages fragmentation events. However, there are not sufficient studies of the complex dynamical phenomena of such orbits. These studies can enrich the knowledge about the long-term evolution of these orbits, be helpful to propose uncatalogued objects observation and identification, and also set the protected region as well as active debris removal. In this paper, the characteristics of 2:1 resonance of Molniya satellite orbits are studied. A large set of numerical simulations, including all the relevant perturbations, is carried out to further investigate the main characteristics, and special attention is payed to the dynamical evolution of objects with high A/M, particularly affected by the direct solar radiation pressure. The long-term dynamical evolution of orbital elements, as well as the dependency of lifetime on the A/M value, is discussed.     

Dynamical evolution of high area-to-mass ratio debris released into GPS orbits

A large set of simulations, including all the relevant perturbations, was carried out to investigate the long-term dynamical evolution of fictitious high area-to-mass ratio (A/M) objects released, with a negligible velocity variation, in each of the six orbital planes used by Global Positioning System (GPS) satellites. As with similar objects discovered in near synchronous trajectories, long lifetime orbits, with mean motions of about 2 revolutions per day, were found possible for debris characterized by extremely high area-to-mass ratios. Often the lifetime exceeds 100 years up to A/M ∼ 45 m²/kg, decreasing rapidly to a few months above such a threshold. However, the details of the evolution, which are conditioned by the complex interplay of solar radiation pressure and geopotential plus luni-solar resonances, depend on the initial conditions. Different behaviors are thus possible. In any case, objects like those discovered in synchronous orbits, with A/M as high as 20–40 m²/kg, could also survive in this orbital regime, with semi-major axes close to the semi-synchronous values, with maximum eccentricities between 0.3 and 0.7, and with significant orbit pole precessions (faster and wider for increasing values of A/M), leading to inclinations between 30° and more than 90°.     

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.     

Semi-analytical theory of mean orbital motion for geosynchronous space debris under gravitational influence

This paper provides a hamiltonian formulation of the equations of motion of an artificial satellite or space debris orbiting the geostationary ring. This theory of order 1 has been formulated using canonical and non-singular elements for eccentricity and inclination. The analysis is based on an expansion in powers of the eccentricity and of the inclination. The theory accounts for the influence of the Earth gravity field expanded in spherical harmonics, paying a particular attention to the resonance occurring for geosynchronous objects. The luni-solar perturbations are also taken into account. We present the resonant motion and its main characteristics: equilibria, stability, fundamental frequencies and width of the resonant area by comparison with a basic analytical model. Finally, we show some results concerning the long term dynamics of a typical space debris under the influence of the gravitational field of the Earth and the luni-solar interactions.     

Single-event upset in geostationary transfer orbit during solar-activity maximum period measured by the Tsubasa satellite

This paper reports single-event upset (SEU) occurrence related to the space radiation environment in geostationary transfer orbit during solar-activity maximum period measured by the Tsubasa satellite. Most SEUs are measured in the inner radiation belt, indicating that they are mainly caused by trapped protons. Thus, the spatial distribution and the temporal variation of the SEU count correlate well with those of trapped protons. The peak SEU rate appears around L = 1.4. The transition point from SEUs caused by trapped protons to those caused by galactic cosmic rays is around L = 2.6. During the experiment period, increased SEU count was sometimes detected due to solar and geomagnetic events outside the inner radiation belt.     

Global dynamics of high area-to-mass ratios GEO space debris by means of the MEGNO indicator

In this paper we provide an extensive analysis of the global dynamics of high-area-to-mass ratios geosynchronous (GEO) space debris, applying a recent technique developed by Cincotta and Simó [Cincotta, P.M., Simó, C.Simple tools to study global dynamics in non-axisymmetric galactic potentials–I. Astron. Astrophys. (147), 205–228, 2000.], Mean Exponential Growth factor of Nearby Orbits (MEGNO), which provides an efficient tool to investigate both regular and chaotic components of the phase space.     

Analytical and semi-analytical investigations of geosynchronous space debris with high area-to-mass ratios

This paper provides a Hamiltonian formulation of the averaged equations of motion with respect to short periods (1 day) of a space debris subjected to direct solar radiation pressure and orbiting near the geostationary ring. This theory is based on a semi-analytical theory of order 1 regarding the averaging process, formulated using canonical and non-singular elements for eccentricity and inclination. The analysis is based on an expansion in powers of the eccentricity and of the inclination, truncated at an arbitrary high order.     

The classical Laplace plane as a stable disposal orbit for geostationary satellites

The classical Laplace plane is a frozen orbit, or equilibrium solution for the averaged dynamics arising from Earth oblateness and lunisolar gravitational perturbations. The pole of the orbital plane of uncontrolled GEO satellites regress around the pole of the Laplace plane at nearly constant inclination and rate. In accordance with Friesen et al. (1993), we show how this stable plane can be used as a robust long-term disposal orbit. The current graveyard regions for end-of-life retirement of GEO payloads, which is several hundred kilometers above GEO depending on the spacecraft characteristics, cannot contain the newly discovered high area-to-mass ratio debris population. Such objects are highly susceptible to the effects of solar radiation pressure exhibiting dramatic variations in eccentricity and inclination over short periods of time. The Laplace plane graveyard, on the contrary, would trap this debris and would not allow these objects to rain down through GEO. Since placing a satellite in this inclined orbit can be expensive, we discuss some alternative disposal schemes that have acceptable cost-to-benefit ratios.     

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.     

大面质比空间碎片在太阳光压和引力作用下的轨道演化

分析了较高轨道(a > 10000km)大面质比空间碎片的轨道动力学演化问题. 重点讨论了位于地球同步轨道的空间碎片轨道演化问题, 并给出轨道偏心率 随时间演化的表达式. 通过进一步分析得出, 倾角大于63°26'的GTO轨 道空间碎片, 仅在J₂和第三体摄动影响下, 会出现轨道偏心率升高; 而对 于大面质比空间碎片, 在J₂项和太阳光压同时作用下, 当近地点指向的角 变率与太阳平黄经变化率接近时, 会出现长期共振现象, 导致轨道偏心率升 高, 近地点降低. 分析还得出, 轨道演化过程中, 偏心率的最大值与初始轨 道近地点的指向有关.     

On the practical exploitation of perturbative effects in low Earth orbit for space debris mitigation

This paper presents the results of a numerical evaluation of the natural lifetime reduction in low Earth orbit, due to dynamical perturbations. The study considers two values for the area-to-mass ratio, a nominal ratio which resembles a typical value of spacecraft in orbit today, and an enhanced ratio which covers the surface augmentation. The results were obtained with two orbit propagators, one of a semi-analytical nature and the second one using non-averaged equations of motion. The simulations for both propagators were set up similarly to allow comparison. They both use the solar radiation pressure and the secular terms of the geopotential ($J_2,J_4$ and $J_6$). The atmospheric drag was turned on and off in both propagators to alternatively study the eccentricity build up and the residual lifetime. The non-averaging case also covers a validation with the full 6?×?6 geopotential. The results confirm the findings in previous publications, that is, the possibility for de-orbiting from altitudes above the residual atmosphere if a solar sail is deployed at the end-of-life, due to the combined effect of solar radiation pressure and the oblateness of the Earth. At near polar inclinations, shadowing effects can be exploited to the same end. The results obtained with the full, non-averaging propagator revealed additional de-orbiting corridors associated with solar radiation pressure which were not found by previous work on space debris mitigation. The results of both tools are compared for specific initial conditions. For nominal values of area-to-mass ratio, instead, it is confirmed that this resonance effect is negligible.The paper then puts the findings in the perspective of the current satellite catalogue. It identifies space missions which are currently close to a resonance corridor and shows the orbit evolution within the resonances with a significantly shorter residual orbital lifetime. The paper finishes with a discussion on the exploitation of these effects with regards to the long-term simulation of the space debris environment and a flux and collision probability comparison.     

Local orbital debris flux study in the geostationary ring

A local orbital debris flux analysis is performed in the geostationary (GEO) ring to investigate how frequently near-miss events occur for each longitude slot in the GEO ring. The current resident space object (RSO) environment at GEO is evaluated, and publicly-available two-line element (TLE) data are utilized in tandem with a geostationary torus configuration to simulate near-miss events incurred by the trackable RSO population at GEO. Methodology for determining near-miss events with this formulation is introduced, and the results of the analysis for a one-year time frame are provided to illustrate the need for active GEO remediation.     

Charge exchange interactions on near-Earth proton radiation for orbit perturbation of high area-to-mass ratio objects

This novel concept expels neutral gas in the presence of geomagnetically-trapped protons in near-Earth orbit. The expelled neutral gas acts to induce charge exchange collisions with the geomagnetically-trapped protons and induce drag on objects which pass through it. The charge exchange collisions between the neutral gas and the geomagnetically-trapped protons create neutrals with similar kinetic energy that are not confined by the geomagnetic field. The charge exchange neutrals are able to collide with orbital objects and perturb their orbits. The delta-v applied by the charge exchange neutral flux is greatest on high area-to-mass objects. Numerical simulation shows charge exchange neutral impacts produce a delta-v on objects on the order of 3.8 x 10⁻¹¹ m/s at a distance of 1 km from the center of the expelled gas in a 1,000 km orbit. The impulse imparted by charge exchange neutral impacts is at least six orders of magnitude smaller than that provided by the induced drag caused by gas expulsion. The localized drag increase can force a majority of small objects into the orbit of the expelled gas cloud, even if that orbit is retrograde to the initial orbit of the objects. This new technique can be applied to the remediation of space debris.     

Practical method to identify orbital anomaly as spacecraft breakup in the geostationary region

Identifying spacecraft breakup events is an essential issue for better understanding of the current orbital debris environment. This paper proposes an observation planning approach to identify an orbital anomaly, which appears as a significant discontinuity in archived orbital history, as a spacecraft breakup. The proposed approach is applicable to orbital anomalies in the geostationary region. The proposed approach selects a spacecraft that experienced an orbital anomaly, and then predicts trajectories of possible fragments of the spacecraft at an observation epoch. This paper theoretically demonstrates that observation planning for the possible fragments can be conducted. To do this, long-term behaviors of the possible fragments are evaluated. It is concluded that intersections of their trajectories will converge into several corresponding regions in the celestial sphere even if the breakup epoch is not specified and it has uncertainty of the order of several weeks.     

A deformation model of flexible,HAMR objects for accurate propagation under perturbations and the self-shadowing effects

A new type of space debris in near geosynchronous orbit (GEO) was recently discovered and later identified as exhibiting unique characteristics associated with high area-to-mass ratio (HAMR) objects, such as high rotation rates and high reflection properties. Observations have shown that this debris type is very sensitive to environmental disturbances, particularly solar radiation pressure, due to the fact that its motion depends on the actual effective area, orientation of that effective area, reflection properties and the area-to-mass ratio of the object is not stable over time. Previous investigations have modelled this type of debris as rigid bodies (constant area-to-mass ratios) or discrete deformed body; however, these simplifications will lead to inaccurate long term orbital predictions. This paper proposes a simple yet reliable model of a thin, deformable membrane based on multibody dynamics. The membrane is modelled as a series of flat plates, connected through joints, representing the flexibility of the membrane itself. The mass of the membrane, albeit low, is taken into account through lump masses at the joints. The attitude and orbital motion of this flexible membrane model is then propagated near GEO to predict its orbital evolution under the perturbations of solar radiation pressure, Earth’s gravity field (J₂), third body gravitational fields (the Sun and Moon) and self-shadowing. These results are then compared to those obtained for two rigid body models (cannonball and flat rigid plate). In addition, Monte Carlo simulations of the flexible model by varying initial attitude and deformation angle (different shape) are investigated and compared with the two rigid models (cannonball and flat rigid plate) over a period of 100?days. The numerical results demonstrate that cannonball and rigid flat plate are not appropriate to capture the true dynamical evolution of these objects, at the cost of increased computational time.     

Investigations of the evolution of optical characteristics and dynamics of proper rotation of uncontrolled geostationary artificial satellites

On the base of the photometrical data an analysis of the periods and light curves of the uncontrolled artificial satellites was performed.     

Optimising filtering of two-line element sets to increase re-entry prediction accuracy for GTO objects

Predicting re-entry epoch of space objects enables managing the risk to ground population. Predictions are particularly difficult for objects in highly-elliptical orbits, and important for objects with components that can survive re-entry, e.g. rocket bodies (R/Bs). This paper presents a methodology to filter two-line element sets (TLEs) to facilitate accurate re-entry prediction of such objects. Difficulties in using TLEs for precise analyses are highlighted and a set of filters that identifies erroneous element sets is developed. The filter settings are optimised using an artificially generated TLE time series. Optimisation results are verified on real TLEs by analysing the automatically found outliers for exemplar R/Bs. Based on a study of 96 historical re-entries, it is shown that TLE filtering is necessary on all orbital elements that are being used in a given analysis in order to avoid considerably inaccurate results.     

Optimal guidance based on receding horizon control for low-thrust transfer to libration point orbits

This paper addresses the design and computation of a guidance law for a transfer mission from an orbit near the Earth to a halo orbit around the libration point L₂ in the Sun–Earth system. The guidance law, which is designed based on receding horizon control and compensates for launch velocity errors that are introduced by inaccuracies of the launch vehicle, is solved using the generating function method. During the design of the closed-loop guidance law, the entire transfer mission, which is considered a nonlinear optimal control problem, is evaluated to obtain a nominal reference trajectory. Using the launch velocity errors and the uncertainty of the model, a spacecraft controlled by the proposed guidance law tracks the reference trajectory. Furthermore, the original Riccati differential equation in the receding horizon control algorithm is replaced by an equivalent convenient form of the Riccati differential equation that is based on the generating function. The high-efficiency solution of the equivalent equation avoids the online direct integration of the original Riccati differential equation, which significantly increases the computational efficiency for the receding horizon control problem. Numerical simulations using a nonlinear bicircular four-body model demonstrate the capabilities of the proposed receding horizon guidance law for the transfer mission. In addition, the generating function method improves the computational efficiency by at least one order of magnitude over the backward sweep method in solving the receding horizon control problem.     

宽高比对微小通道空气流动换热特性影响实验

针对不同宽高比微小矩形通道内流动和换热特性进行了实验研究.微小矩形通道宽0.4mm,宽高比分别为2/5、1/2、2/3和1.实验以空气作为工质.流动雷诺数Re范围是200~1600.实验结果表明:宽高比对微小矩形通道流动换热特性的影响不可忽略.微小矩形通道摩擦因子随着Re增大而减小,泊肃叶数Po随着Re增大而增大,二者随着宽高比增大同时减小;努塞尔数Nu随着Re增大而增大,随着宽高比增大而减小.针对宽高比对微小矩形通道流动换热特性的影响拟合了相应经验关系式.