Deorbiter CubeSat mission design

This paper presents the mission design for a CubeSat-based active debris removal approach intended for transferring sizable debris objects from low-Earth orbit to a deorbit altitude of 100 km. The mission consists of a mothership spacecraft that carries and deploys several debris-removing nanosatellites, called Deorbiter CubeSats. Each Deorbiter is designed based on the utilization of an eight-unit CubeSat form factor and commercially-available components with significant flight heritage. The mothership spacecraft delivers Deorbiter CubeSats to the vicinity of a predetermined target debris, through performing a long-range rendezvous maneuver. Through a formation flying maneuver, the mothership then performs in-situ measurements of debris shape and orbital state. Upon release from the mothership, each Deorbiter CubeSat proceeds to performing a rendezvous and attachment maneuver with a debris object. Once attached to the debris, the CubeSat performs a detumbling maneuver, by which the residual angular momentum of the CubeSat-debris system is dumped using Deorbiter’s onboard reaction wheels. After stabilizing the attitude motion of the combined Deorbiter-debris system, the CubeSat proceeds to performing a deorbiting maneuver, i.e., reducing system’s altitude so much so that the bodies disintegrate and burn up due to atmospheric drag, typically at around 100 km above the Earth surface. The attitude and orbital maneuvers that are planned for the mission are described, both for the mothership and Deorbiter CubeSat. The performance of each spacecraft during their operations is investigated, using the actual performance specifications of the onboard components. The viability of the proposed debris removal approach is discussed in light of the results.     

Reentry survival analysis of tumbling metallic hollow cylinder

The survival of orbital debris reentering the Earth’s atmosphere is considered. The numerical approach of NASA’s Object Reentry Survival Analysis Tool (ORSAT) is reviewed, and a new equation accounting for reradiation heat loss of hollow cylindrical objects is presented. Based on these, a code called Survivability Analysis Program for Atmospheric Reentry (SAPAR) has been developed, and the new equation for reradiation heat loss is validated. Using this equation in conjunction with the formulation used in ORSAT, a comparative case study on the Delta-II second stage cylindrical tank is given, demonstrating that the analysis using the proposed equation is in good agreement with the actual recovered object when a practical value for thermal emissivity is used. A detailed explanation of the revised formulation is given, and additional simulation results are presented. Finally, discussions are made to address the applicability of the proposed equation to be incorporated in future survival analyses of orbital debris.     

基于雷诺数的大气阻力模型在飞行器再入预报中的应用

随着对空间技术服务需求的增加和空间碎片主动移除技术的实现,未来空间碎片将以数量多、质量大、难分解等特点频繁再入大气层,给地面人员和财产安全造成更多威胁.因此,亟需对火箭体等大型航天器的大气再入进行预警,然而因缺乏合适的大气阻力系数模型难以实现高精度的大气再入预报.为此,在简化航天器模型的基础上引入基于雷诺数的大气动力模...     

An adaptive threshold method for improving astrometry of space debris CCD images

Optical survey is a main technique for observing space debris, and precisely measuring the positions of space debris is of great importance. Due to several factors, e.g. the angle object normal to the observer, the shape as well as the attitude of the object, the variations of observed characteristics for low earth orbital space debris are distinct. When we look at optical CCD images of observed objects, the size and brightness are varying, hence it’s difficult to decide the threshold during centroid measurement and precise astrometry. Traditionally the threshold is given empirically and constantly in data reduction, and obviously it’s not suitable for data reduction of space debris. Here we offer a solution to provide the threshold. Our method assumes that the PSF (point spread function) is Gaussian and estimates the signal flux by a directly two-dimensional Gaussian fit, then a cubic spline interpolation is performed to divide each initial pixel into several sub-pixels, at last the threshold is determined by the estimation of signal flux and the sub-pixels above threshold are separated to estimate the centroid. A trail observation of the fast spinning satellite Ajisai is made and the CCD frames are obtained to test our algorithm. The calibration precision of various threshold is obtained through the comparison between the observed equatorial position and the reference one, the latter are obtained from the precise ephemeris of the satellite. The results indicate that our method reduces the total errors of measurements, it works effectively in improving the centering precision of space debris images.     

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).     

Apparent rotation properties of space debris extracted from photometric measurements

Knowledge about the rotation properties of space debris objects is essential for the active debris removal missions, accurate re-entry predictions and to investigate the long-term effects of the space environment on the attitude motion change. Different orbital regions and object’s physical properties lead to different attitude states and their change over time.Since 2007 the Astronomical Institute of the University of Bern (AIUB) performs photometric measurements of space debris objects. To June 2016 almost 2000 light curves of more than 400 individual objects have been acquired and processed. These objects are situated in all orbital regions, from low Earth orbit (LEO), via global navigation systems orbits and high eccentricity orbit (HEO), to geosynchronous Earth orbit (GEO). All types of objects were observed including the non-functional spacecraft, rocket bodies, fragmentation debris and uncorrelated objects discovered during dedicated surveys. For data acquisition, we used the 1-meter Zimmerwald Laser and Astrometry Telescope (ZIMLAT) at the Swiss Optical Ground Station and Geodynamics Observatory Zimmerwald, Switzerland. We applied our own method of phase-diagram reconstruction to extract the apparent rotation period from the light curve. Presented is the AIUB’s light curve database and the obtained rotation properties of space debris as a function of object type and orbit.     

一种基于TLE数据的轨道异常分析方法

空间在轨物体的轨道异常是航天工程及预警领域普遍关注的问题,及时发现轨道异常意义重大,通过分析空间物体的轨道异常,可以及时发现和识别规避事件或碰撞事件,还可以了解监测网的能力.本文提出一种基于TLE数据的简单的轨道异常分析方法——长半轴变化法.该方法快速有效,应用到低轨在用卫星和美俄解体碎片的异常分析中,异常物体正确识别率可达到100%;对美俄解体碎片进行轨道异常分析后得出,美国空间监视网可以稳定探测90%以上的解体碎片.      

基于相控阵雷达波束篱笆的空间碎片数量与分布估计方法

随着载人航天与空间站等航天活动的增多,不能有效防护、也无法定期跟踪和编目的小尺寸(尤其是1～10 cm)碎片的危害越来越受到关注,这些碎片信息的获取依赖于统计采样技术.针对简化的相控阵雷达波束篱笆空间碎片探测模式,提出了一种采用统计技术估计空间碎片总数量以及高度和倾角分布的方法.将碎片穿越波束篱笆的过程用Poisson分布来建模,根据观测时段内穿越波束篱笆目标的平均到达率及测量的轨道高度和倾角数据来估计给定轨道高度范围或倾角范围内碎片的数量,进而得到碎片的总数量以及碎片数量随轨道高度或倾角的分布.在获取雷达散射截面信息时,该方法还可用于估计碎片数量随尺寸的分布.通过仿真实验验证了该方法的有效性.      

The fast debris evolution model

The ‘particles-in-a-box’ (PIB) model introduced by Talent [Talent, D.L. Analytic model for orbital debris environmental management. J. Spacecraft Rocket, 29 (4), 508–513, 1992.] removed the need for computer-intensive Monte Carlo simulation to predict the gross characteristics of an evolving debris environment. The PIB model was described using a differential equation that allows the stability of the low Earth orbit (LEO) environment to be tested by a straightforward analysis of the equation’s coefficients. As part of an ongoing research effort to investigate more efficient approaches to evolutionary modelling and to develop a suite of educational tools, a new PIB model has been developed. The model, entitled Fast Debris Evolution (FADE), employs a first-order differential equation to describe the rate at which new objects ?10 cm are added and removed from the environment. Whilst Talent [Talent, D.L. Analytic model for orbital debris environmental management. J. Spacecraft Rocket, 29 (4), 508–513, 1992.] based the collision theory for the PIB approach on collisions between gas particles and adopted specific values for the parameters of the model from a number of references, the form and coefficients of the FADE model equations can be inferred from the outputs of future projections produced by high-fidelity models, such as the DAMAGE model.     

A debris image tracking using optical flow algorithm

This study proposes a motion detection and object tracking technique for GEO debris in a sequence of images. A couple of techniques (called the “stacking method” and “line-identifying technique”) were recently proposed to address the same problem. Although these techniques are effective at detecting the debris position and motion in the image sequences, there are some issues concerned with computational load and assumed debris motion. This study derives a method to estimate motion vectors of objects in image sequence and finally detect the debris locations by using a computer vision technique called an optical flow algorithm. The new method detects these parameters in low computational time in a serial manner, which implies that it has an advantage to track not only linear but also nonlinear motion of GEO debris more easily than the previous methods. The feasibility of the proposed methods is validated using real and synthesized image sequences which contain some typical debris motions.     

Space debris cumulative flux considering the Interval Distance-based method

The availability of engineering models to estimate the risk from space debris is essential for space missions. According to current research, cumulative flux calculation is mostly carried out based on the equal-width interval discretization. The method discretizes the volume around the Earth into cells defined in earth centered inertial coordinates. The resulting debris flux onto a target object is shown to depend on the chosen size of the cells. To avoid a discretization error, this must be accounted for. In order to present reliable flux predictions for space mission, the algorithm improvement is an ongoing topic for the related research field. The aim of this study was to examine the discretization error during the cumulative flux determination process. Both the effect of interval step length and the orbital boundary are under investigation. Several typical orbits are selected as examples here and the 2018/01/03 TLE (Two Line Element) data published by the US Space Surveillance Network is used as the debris background in this paper. Furthermore, the Interval Distance-Based method for Discretization (IDD) is adopted in this paper. A position-centered flux determination method is introduced based on the IDD method. According to the example analysis, the IDD used in the flux calculation process provides results which are less affected by the interval step-size setup; and the orbital boundary has no effect on the calculation process. In other words, the discretization error is significantly reduced. The position-centered method provided a possible suggestion for the improvement of space debris environment models.     

满足多终端约束的二次曲线迭代制导方法研究

某些空间载荷会对入轨精度和入轨姿态同时提出很高的要求,应用于运载火箭控制系统的摄动制导方法的入轨精度无法满足要求,而传统的迭代制导方法无法约束终端的入轨姿态。为此提出了一种满足多终端约束的二次曲线迭代制导方法,该方法通过二次曲线形式规划整个真空段的飞行制导程序角,实时满足位置、速度与终端姿态约束,从而可以使得火箭以期望姿态角实现高精度入轨。算例结果表明,该方法能同时保证高精度的入轨指标和入轨姿态,并能适应偏差状态、约束姿态角变化与轨道根数小幅变化,具有较高的工程应用价值。     

星间洛仑兹力编队飞行的平衡点及零速度曲面

提出了一种新的利用星间洛仑兹力控制卫星相对运动的方法:使主星产生自旋磁场,副星带电,通过控制副星所受的星间洛仑兹力进行编队。假设主星产生的人造磁场表现为偶极子并且运动在一条开普勒圆形轨道上,副星恒定带电并在主星附近运动,同时假设星间洛仑兹力只影响副星的运动而不影响主星的运动。推导了副星在主星HCW(当地垂直当地地平坐标系)坐标系下的相对运动动力学方程。针对偶极子与HCW坐标轴X轴重合的情况下,推导了动力系统下的平衡点(并采用稳定性分析方法分析其线性化意义下的稳定性)、积分常数和零速度曲面,证明了有界相对运动的存在。最后用数值仿真验证了上述结论。     

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.     

Radioactive satellites: Intact reentry and breakup by debris impact

There is a substantial mass of radioactive material in nuclear reactors or radioisotope thermal generators (RTGs) in orbit about the earth. This paper examines the reentry of intact nuclear fuel cores and RTGs and the fragmentation and subsequent radioactive debris cloud deposition and evolution resulting from the impact of orbital debris upon an orbiting reactor, fuel core, or RTG. To assess the intact reentry, decay rates and a predicted decay date using historical and projected orbital decay data, are estimated. The current NASA debris environment model is utilized to estimate impact rates and debris cloud evolution of a fragmentation event. Results of these analyses are compared and concepts are tendered which would tend to minimize the radiological debris hazard to personnel and structures both on the earth's surface and in low earth orbit.     

提高卫星自主导航精度的滤波算法和仿真

如何实时准确地对系统误差进行估计和校准提高基于“地-日-月”测量的卫星自主导航精度的一个重要途径。针对这一问题,将系统误差作为新增的状态向量,设计了轨道参数状态与系统误差状态解耦估计的两步滤波算法,该算法在估计轨道参数的同时校准系统误差,数学仿真结果表明,采用文章所提出的滤波算法可以将自主导航的位置精度从10km(3σ)量级提高到300m(3σ)左右。     

Practical guidelines for electro-dynamic tethers to survive from orbital debris impacts

Japan Aerospace Exploration Agency (JAXA) has proposed an active debris removal using electro-dynamic tether to reduce large space debris in the low-Earth orbit. However, a tether strand is thin but long enough to have a large area so that it is vulnerable to small particles. This vulnerability might be the weakest point of a tether system against orbital debris. In order to overcome this weakest point, a double tether system, in which two tether strands are tied together at even intervals to form equally spaced loops, has been suggested as one of the promising candidates. This paper provides a mathematical approach to estimate the survival probability of a double tether system and then apply the approach to evaluate the mission success rate of the active debris removal using electro-dynamic tether that JAXA has proposed. It can be concluded the countermeasure to get enough success rate can be obtained. The result is simulated for Advanced Earth Observing Satellite II (ADEOS-II) re-entry from 800 km sun synchronized orbit to atmosphere. The simulation shows that mission success rate over 90% can be obtained with number of loops over 1000 and 10 mm clearance between two strands.     

空间碎片云演变过程的阶段划分

根据碎片云从破碎点开始向空间扩散过程中碎片密度和形状的变化规律,以几何形状和起主要作用的因素为特征,定义了球形、椭球形、绳形、螺旋线形、全方位弥漫直至球壳形六个演变阶段.论述了在各个阶段的主要特征和对演变过程起主要作用的因素.总结了与演变过程相关的轨道运动理论和研究方法,分析了各个阶段演变的动力学原理.在球形阶段起主要作用的是分离速度;椭球形阶段可以利用线性化相对运动方程进行分析;绳形与螺旋线形在几何上有质变,但都有结点和结线,并可以利用速度增量理论分析和解释其存在的原因.轨道摄动力消除了结点和结线,导致碎片云的全方位弥漫,并最终使碎片云趋于球壳形.推导和罗列了各阶段转换标志点时刻的计算公式,利用计算机仿真的方法,给出了近地轨道各个阶段碎片云分布示意图,验证了演变过程阶段划分的合理性.     

基于解体事件的空间碎片轨道演化算法研究

针对航天器解体事件所生成的空间碎片的演化过程,进行了数学分析,确定了新生成的空间碎片的速度增量,在该增量作用下碎片轨道会发生变更,本文根据该增量得出了空间碎片在轨道变更后的轨道根数,分析了在大气阻力摄动作用下,空间碎片的数目和轨道分布的演化情况,给出了相关结果,结果表明此算法可行。     

空间碎片的空间密度算法研究

由于空间碎片的影响, 空间环境日益恶劣,有必要建立空间碎片环境工程模型对空间碎片撞击航天器进行风险评估.本文研究了空间碎片环境模型中的重要环节, 即碎片的空间密度问题,在统计理论与椭圆轨道理论基础上,分析及推导了空间碎片在空间中的分布状况及空间密度,并得出了碎片空间密度的空间坐标函数解析表达式.同时对引入假设条件的合理性进行了讨论,并利用双行元数据对结论进行了验证, 利用本文方法得出的结果与双行元数据吻合.