Large constellations assessment and optimization in LEO space debris environment

Recent plans for large constellations in Low-Earth Orbit have opened the debate on both their vulnerability and their influence on the already hazardous space debris environment. In fact, given that large constellations normally employ satellites of small size, there might be situations in which cm-size debris could have enough energy to cause fragmentation of a significant part of these spacecraft upon impact, while smaller debris could affect the functionalities of critical subsystems, even compromising the success of disposal operations planned at end-of-life. In this context, this paper investigates: (1) collisions with large objects that could initiate the fragmentation of a significant part of the satellite, and (2) impacts with small debris that might perforate the spacecraft hull thus causing relevant performance/functionality degradation. These two points are merged in a simple statistical tool for risk assessment, which analyses the effects of the main parameters of the constellations on its vulnerability (i.e. operational life, number of satellites, spacecraft cross section, satellites reliability). In more details, the tool relates impact probability (for both small and large debris) to the ballistic response of spacecraft structures and protections, defining the critical configurations that might compromise the expected disposal operations. This method requires a limited knowledge of the spacecraft internal layout, as it is based on a statistical analysis of impact damage instead of a complete evaluation of the vulnerability of each subsystem. In parallel, non-debris related failures are also investigated and statistic models of spacecraft reliability characteristic are proposed. Among the results, it is shown that reducing the lifetime of individual satellites in a constellation might improve the success rate of post-mission disposal, thanks to the reduction of the spacecraft exposure to the space environment with the consequential degradation of its performance. On the other hand, reducing the lifetime would seriously affect the debris environment: the increase in traffic in the most crowded altitudes would be not counterbalanced by the higher post mission disposal success rate, causing an overall increase of the total number of uncontrolled resident objects.     

The situation in the geostationary ring

For more than 25 years, the practice of reorbiting of a geostationary satellite at the end of its mission in order to protect the geostationary orbit (GEO) environment has been recommended and performed by a number of operators. In recent years, an internationally recognised re-orbiting altitude has been defined by the Inter-Agency Space Debris Coordination Committee (IADC). Based on orbital data contained in the DISCOS database, the situation on the geostationary ring is analysed. In January 2004, from 1036 catalogued objects passing through the geostationary region, 340 are controlled within their allocated longitude slots, 395 are drifting above, below or through GEO, and 140 are in a libration orbit. In the periods 1997–2003, 103 spacecraft reached their end of life; 34 were reorbited in compliance with the IADC recommendation, 35 were reorbited below the minimum recommended altitude, and 34 were abandoned without any end-of-life disposal manoeuvre. Apart from these catalogued objects, the ESA 1-m telescope has observed many smaller debris (down to 10–15 cm) in this orbital region representing a collision risk for GEO spacecraft which is difficult to quantify.     

Debris swarms seen by SMEI

The large 3° × 60° fields-of-view of the Solar Mass Ejection Imager (SMEI) instruments are oriented on the stabilized Coriolis satellite to image most of the sky each Sun-synchronous orbit. Besides observing coronal mass ejections, the SMEI mission objective, SMEI also has detected a plethora of Earth-orbiting satellites (resident space objects or RSOs) brighter than ∼8th magnitude at a rate of about 1 per minute. Occasionally, SMEI sees an RSO swarm: a sudden onset of a large number of RSOs, many more than the nominal rate, upto dozens detected in a 4-s frame. These swarms usually last for a few minutes. A sample of six such RSO ensembles is analyzed in this paper in which the distance and the direction of the velocity vector for individual objects are estimated. We present the observational evidence indicating that the swarms must be near-field objects traveling in orbits near that of Coriolis, and that the relatively speeds between the objects and Coriolis are low. Further, analyses indicate that the RSOs are quite close (<20 m) and are generally moving radially away from the satellite. The predicted encounter geometries for Coriolis passing through or near a small debris cloud is, generally, quite inconsistent with the observations. The most likely explanation consistent with the observations is that SMEI is seeing debris being ejected from the Coriolis spacecraft itself. An analysis of distance and brightness for a subset of the RSOs indicates that the median diameter of the debris particles is ∼80 μm.     

Predicting the vulnerability of spacecraft components: Modelling debris impact effects through vulnerable-zones

The space environment around the Earth is populated by more than 130 million objects of 1 mm in size and larger, and future predictions shows that this amount is destined to increase, even if mitigation measures are implemented at a far better rate than today. These objects can hit and damage a spacecraft or its components. It is thus necessary to assess the risk level for a satellite during its mission lifetime. Few software packages perform this analysis, and most of them employ time-consuming ray-tracing methodology, where particles are randomly sampled from relevant distributions. In addition, they tend not to consider the risk associated with the secondary debris clouds. The paper presents the development of a vulnerability assessment model, which relies on a fully statistical procedure: the debris fluxes are directly used combining them with the concept of vulnerable zone, avoiding the random sampling the debris fluxes. A novel methodology is presented to predict damage on internal components. It models the interaction between the components and the secondary debris cloud through basic geometrical operations, considering mutual shielding and shadowing between internal components. The methodologies are tested against state-of-the-art software for relevant test cases, comparing results on external structures and internal components.     

控制力矩陀螺框架控制方法及框架转速测量方法

控制力矩陀螺是一种具有力矩放大特性的惯性执行机构,通常应用于大型航天器姿态控制。近年来,随着相关技术的发展,对基于控制力矩陀螺的中小卫星快速机动平台的需求日益迫切,这不仅需要控制力矩陀螺能够输出大力矩,而且要具有较高的力矩输出精度。本文结合工程实践,提出一种框架转速精度测量方法,以及一种采用正弦永磁同步电机,基于转子磁场定向的矢量控制方案。该框架控制方案中引入摩擦力矩观测器及其补偿算法,在控制回路中通过对摩擦力矩的补偿,可有效提高框架驱动控制系统的稳定性和动态性能。     

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.     

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.     

Simulation of long-term rotational dynamics of large space debris: A TOPEX/Poseidon case study

Accurate knowledge of the rotational dynamics of a large space debris is crucial for space situational awareness (SSA), whether it be for accurate orbital predictions needed for satellite conjunction analyses or for the success of an eventual active debris removal mission charged with stabilization, capture and removal of debris from orbit. In this light, the attitude dynamics of an inoperative satellite of great interest to the space debris community, the joint French and American spacecraft TOPEX/Poseidon, is explored. A comparison of simulation results with observations obtained from high-frequency satellite range measurements is made, showing that the spacecraft is currently spinning about its minor principal axis in a stable manner. Predictions of the evolution of its attitude motion to 2030 are presented, emphasizing the uncertainty on those estimates due to internal energy dissipation, which could cause a change of its spin state in the future. The effect of solar radiation pressure and the eddy-current torque are investigated in detail, and insights into some of the satellite’s missing properties are provided. These results are obtained using a novel, open-source, coupled orbit-attitude propagation software, the Debris SPin/Orbit Simulation Environment (D-SPOSE), whose primary goal is the study of the long-term evolution of the attitude dynamics of large space debris.     

Characterizing space debris longitude-dependent distribution based on RAAN perturbation rate

The space debris environment is one of the major threats against payloads. Space debris orbital distribution is of great importance for space debris environment modeling. Due to perturbation factors, the Right Ascension of Ascending Node (RAAN) of space objects changes consistently, causing regular rotation of the orbit plane around Earth’s axis. Based on the investigation of the RAAN perturbation rate of concerned objects, this paper proposes a RAAN discretization method in order to present the space debris longitude-dependent distribution. Combined with two line element (TLE) data provided by the US Space Surveillance Network, the estimated value from RAAN discretization method is compared with the real case. The results suggest that using only the initial orbital data at the beginning of the time interval of interest, the RAAN discretization method is able to provide reliable longitude distribution of concerned targets in the next following period. Furthermore, spacecraft cumulative flux against space debris is calculated in this paper. The results suggest that the relevance between spacecraft RAAN setup and flux output is much smaller for LEO targets than MEO targets, which corresponds with the theory analysis. Since the nonspherical perturbation is the major factor for RAAN variation, the RAAN perturbation rate has little connection with the size of orbital objects. In other words, the RAAN discretization method introduced in this paper also applies to space debris of different size range, proposing a possible suggestion for the improvement of space debris environment engineering models.     

Impact effects from small size meteoroids and space debris

When the impact risk from meteoroids and orbital debris is assessed the main concern is usually structural damage. With their high impact velocities of typically 10–20 km/s millimeter or centimeter sized objects can puncture pressure vessels and other walls or lead to destruction of complete subsystems or even whole spacecraft. Fortunately chances of collisions with such larger objects are small (at least at present). However, particles in the size range 1–100 μm are far more abundant than larger objects and every orbiting spacecraft will encounter them with certainty. Every solar cell (8 cm² area) of the Hubble Space Telescope encountered on average 12 impacts during its 8.25 years of space exposure. Most were from micron sized particles.     

利用TLE数据分析LEO卫星轨道异常的新方法——————综合判据法

及时准确地发现在轨卫星的轨道异常意义重大. 通过有效的异常算法, 能够找出发生轨道异常的碎片或航天器, 为空间碎片碰撞预警系统分析和验证碰撞事件提供数据支持. 通过对利用TLE (Two Line Elements)数据分析LEO在轨卫星轨道异常的方法研究, 提出了一个利用单个卫星相邻根数时间差控制加综合判据的判别方法. 分析表明, 相对于取单一因素阈值的判别方法, 综合判据法能够最大限度地减少漏判, 并且保持相对较高的判断准确率.      

Dynamics and control for contactless interaction between spacecraft and tumbling debris

Tumbling debris has become a great threat to orbit activities. Contactless interaction is a novel concept for active debris removal, through which the tumbling debris no longer rotates freely but is under control. The contactless interaction method aims to de-tumble the debris and then maintain desired relative states between the spacecraft and debris. The spacecraft is simultaneously stabilized through three-axis attitude control, which makes the de-tumbling and capture operation much safer, more effective and accurate. The dynamics and control for the contactless interaction have been little studied in the past years. This paper considers a generic dynamics and control problem for contactless interaction between a spacecraft and debris. A translational and rotational dynamics model of contactless interaction is proposed and the 6-DOF equations are established. The contactless interaction control law is designed with the backstepping method, and the spacecraft three-axis control law is designed with the PD control. Simulation results show that the angular momentum is transferred from the debris to the spacecraft and the debris is thus de-tumbled. The desired relative states are achieved efficiently. Significantly, the spacecraft and debris no longer rotate in the inertial frame and, hence, the safety and accuracy for capture operation are guaranteed.     

A deorbiter CubeSat for active orbital debris removal

This paper introduces a mission concept for active removal of orbital debris based on the utilization of the CubeSat form factor. The CubeSat is deployed from a carrier spacecraft, known as a mothership, and is equipped with orbital and attitude control actuators to attach to the target debris, stabilize its attitude, and subsequently move the debris to a lower orbit where atmospheric drag is high enough for the bodies to burn up. The mass and orbit altitude of debris objects that are within the realms of the CubeSat’s propulsion capabilities are identified. The attitude control schemes for the detumbling and deorbiting phases of the mission are specified. The objective of the deorbiting maneuver is to decrease the semi-major axis of the debris orbit, at the fastest rate, from its initial value to a final value of about 6471?km (i.e., 100?km above Earth considering a circular orbit) via a continuous low-thrust orbital transfer. Two case studies are investigated to verify the performance of the deorbiter CubeSat during the detumbling and deorbiting phases of the mission. The baseline target debris used in the study are the decommissioned KOMPSAT-1 satellite and the Pegasus rocket body. The results show that the deorbiting times for the target debris are reduced significantly, from several decades to one or two years.     

航天器的低轨道卫星网络接入

随着空间技术的发展,航天器网络化的趋势不断增强。航天器不仅需要将数据实时地传输到地面,而且需要与其他航天器实现信息共享。通过建设LEO卫星网,容易实现航天器的数据传输与共享,但是航天器的网络接入是研究难点之一。如果针对不同系统采用不同的接入方法,将增加接入的复杂度。通过采用区域路由方法,提出一种通用的接入方法,不仅简化了航天器卫星网络的接入,也简化卫星网络系统的设计,实现了接入方法与网络设计的适当分离。计算机模拟表明该方法的可行性。     

Finite element reconstruction approach for on-orbit spacecraft breakup dynamics simulation and fragment analysis

Breakup model is the key area of space debris environment modeling. NASA standard breakup model is currently the most widely used for general-purpose. It is a statistical model found based on space surveillance data and a few ground-based test data. NASA model takes the mass, impact velocity magnitude for input and provides the fragment size, area-to-mass ratio, velocity magnitude distributions for output. A more precise approach for spacecraft disintegration fragment analysis is presented in this paper. This approach is based on hypervelocity impact dynamics and takes the shape, material, internal structure and impact location etc. of spacecraft and impactor, which might greatly affect the fragment distribution, into consideration. The approach is a combination of finite element and particle methods, entitled finite element reconstruction (FER). By reconstructing elements from the particle debris cloud, reliable individual fragments are identified. Fragment distribution is generated with undirected graph conversion and connected component analysis. Ground-based test from literature is introduced for verification. In the simulation satellite targets and impactors are modeled in detail including the shape, material, internal structure and so on. FER output includes the total number of fragments and the mass, size and velocity vector of each fragment. The reported fragment distribution of FER shows good agreement with the test, and has good accuracy for small fragments.     

基于航天器解体事件的空间碎片寿命算法

研究了针对航天器解体事件所生成的空间碎片的寿命计算方法.给出了基于NASA标准航天器解体模型的航天器解体算法.该算法生成的一系列碎片参数,将作为寿命计算的初始条件.总结了现有求解碎片寿命的算法,并提出了一种半分析算法.该算法运用平均根数法的思路,计算了在J2摄动项的影响下,碎片的半长轴和偏心率的变化率;并采用微分积分法预报半长轴和偏心率随时间的变化.为了适应时变大气模型,该算法限制了计算步长.通过与数值法的比较分析了算法的计算速度和精度.选用了3种大气模型:SA76、GOST和MSIS-00,分析了不同大气模型在计算碎片寿命之间的差异.通过与P-78卫星解体事件的实测数据对比验证了整个算法的正确性.      

Analysis of laser ranges and angular measurements data fusion for space debris orbit determination

In the framework of space debris, the orbit determination process is a fundamental step, both, for researchers and for satellite operators. The accurate knowledge of the orbit of space debris objects is needed to allow space debris characterization studies and to avoid unnecessary collision avoidance maneuvers.The accuracy of the results of an orbit determination process depends on several factors as the number, the accuracy, the kind of processed measurements, their distribution along the orbit, and the object-observer relative geometry. When the observation coverage of the target orbit is not homogeneous, the accuracy of the orbit determination can be improved processing different kind of observables. Recent studies showed that the satellite laser ranging technique can be successfully applied to space debris.In this paper, we will investigate the benefits of using laser ranges and angular measurements for the orbit determination process. We will analyze the influence of the number of used observations, of the covered arc of orbit, of each observable, and of the observation geometry on the estimated parameters. Finally, using data acquired on short observation arcs, we analyze the achievable accuracies for the orbital regimes with the highest space debris density, and to the consequences of the data fusion on catalog maintenance operations. The results shown are obtained using only real data (both angular and laser measurements) provided by sensors of the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald owned by the Astronomical Institute of the University of Bern (AIUB) and for some studies also using ranges provided from other stations of the International Laser Ranging Service (ILRS).     

Analysis of the orbit lifetime of CubeSats in low Earth orbits including periodic variation in drag due to attitude motion

It is estimated that more than 22,300 human-made objects are in orbit around the Earth, with a total mass above 8,400,000 kg. Around 89% of these objects are non-operational and without control, which makes them to be considered orbital debris. These numbers consider only objects with dimensions larger than 10 cm. Besides those numbers, there are also about 2000 operational satellites in orbit nowadays. The space debris represents a hazard to operational satellites and to the space operations. A major concern is that this number is growing, due to new launches and particles generated by collisions. Another important point is that the development of CubeSats has increased exponentially in the last years, increasing the number of objects in space, mainly in the Low Earth Orbits (LEO). Due to the short operational time, CubeSats boost the debris population. One of the requirements for space debris mitigation in LEO is the limitation of the orbital lifetime of the satellites, which needs to be lower than 25 years. However, there are space debris with longer estimated decay time. In LEÓs, the influence of the atmospheric drag is the main orbital perturbation, and is used in maneuvers to increment the losses in the satellite orbital energy, to locate satellites in constellations and to accelerate the decay.The goal of the present research is to study the influence of aerodynamic rotational maneuver in the CubeSat?s orbital lifetime. The rotational axis is orthogonal to the orbital plane of the CubeSat, which generates variations in the ballistic coefficient along the trajectory. The maneuver is proposed to accelerate the decay and to mitigate orbital debris generated by non-operational CubeSats. The panel method is selected to determine the drag coefficient as a function of the flow incident angle and the spinning rate. The pressure distribution is integrated from the satellite faces at hypersonic rarefied flow to calculate the drag coefficient. The mathematical model considers the gravitational potential of the Earth and the deceleration due to drag. To analyze the effects of the rotation during the decay, multiple trajectories were propagated, comparing the results obtained assuming a constant drag coefficient with trajectories where the drag coefficient changes periodically. The initial perigees selected were lower than 400 km of altitude with eccentricities ranging from 0.00 to 0.02. Six values for the angular velocity were applied in the maneuver. The technique of rotating the spacecraft is an interesting solution to increase the orbit decay of a CubeSat without implementing additional de-orbit devices. Significant changes in the decay time are presented due to the increase of the mean drag coefficient calculated by the panel method, when the maneuver is applied, reducing the orbital lifetime, however the results are independent of the angular velocity of the satellite.     

柱状弹丸撞击防护屏形成碎片云材料状态特性研究

微流星体及空间碎片的高速撞击威胁着长寿命、大尺寸航天器的安全运行，导致其严重的损伤和灾难性的失效。为精确估计微流星体及空间碎片高速撞击防护屏所产生碎片云对舱壁的损伤，必须确定碎片云中三种状态材料的特性，建立了碎片云特性分析模型，分别计算了柱状弹丸撞击防护屏所产生碎片云以及碎片云中弹丸和防护屏材料三种状态物质的质量分布。通过计算分析可见，弹丸以不同速度撞击防护屏所产生碎片云三种状态物质的质量分布是不同的，速度增大，液化和气化增强，对靶件的损伤小。而在速度小于7km/s时，碎片云以固体碎片的形式存在，对靶件的损伤大。     

A Partially Orthogonal EnKF approach to atmospheric density estimation using orbital debris

A key requirement for accurate trajectory prediction and space situational awareness is knowledge of how non-conservative forces affect space object motion. These forces vary temporally and spatially, and are driven by the underlying behavior of space weather particularly in Low Earth Orbit (LEO). Existing trajectory prediction algorithms adjust space weather models based on calibration satellite observations. However, lack of sufficient data and mismodeling of non-conservative forces cause inaccuracies in space object motion prediction, especially for uncontrolled debris objects. The uncontrolled nature of debris objects makes them particularly sensitive to the variations in space weather. Our research takes advantage of this behavior by utilizing observations of debris objects to infer the space environment parameters influencing their motion.The hypothesis of this research is that it is possible to utilize debris objects as passive, indirect sensors of the space environment. We focus on estimating atmospheric density and its spatial variability to allow for more precise prediction of LEO object motion. The estimated density is parameterized as a grid of values, distributed by latitude and local sidereal time over a spherical shell encompassing Earth at a fixed altitude of 400 km. The position and velocity of each debris object are also estimated. A Partially Orthogonal Ensemble Kalman Filter (POEnKF) is used for assimilation of space object measurements to estimate density.For performance comparison, the scenario characteristics (number of objects, measurement cadence, etc.) are based on a sensor tasking campaign executed for the High Accuracy Satellite Drag Model project. The POEnKF analysis details spatial comparisons between the true and estimated density fields, and quantifies the improved accuracy in debris object motion predictions due to more accurate drag force models from density estimates. It is shown that there is an advantage to utilizing multiple debris objects instead of just one object. Although the work presented here explores the POEnKF performance when using information from only 16 debris objects, the research vision is to utilize information from all routinely observed debris objects. Overall, the filter demonstrates the ability to estimate density to within a threshold of accuracy dependent on measurement/sensor error. In the case of a geomagnetic storm, the filter is able to track the storm and provide more accurate density estimates than would be achieved using a simple exponential atmospheric density model or MSIS Atmospheric Model (when calm conditions are assumed).