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.     

Long term dynamics and optimal control of nano-satellite deorbit using a short electrodynamic tether

This paper studies the long term dynamics and optimal control of a nano-satellite deorbit by a short electrodynamic tether. The long term deorbit process is discretized into intervals and within each interval a two-phase optimal control law is proposed to achieve libration stability and fast deorbit simultaneously. The first-phase formulates an open-loop fast-deorbit control trajectory by a simplified model that assumes the slow-varying orbital elements of electrodynamic tethered system as constant and ignores perturbation forces other than the electrodynamic force. The second phase tracks the optimal trajectory derived in the first phase by a finite receding horizon control method while considering a full dynamic model of electrodynamic tether system. Both optimal control problems are solved by direct collocation method base on the Hermite–Simpson discretization schemes with coincident nodes. The resulting piecewise nonlinear programing problems in the sequential intervals reduces the problem size and improve the computational efficiency, which enable an on-orbit control application. Numerical results for deorbit control of a short electrodynamic tethered nano-satellite system in both equatorial and highly inclined orbits demonstrate the efficiency of the proposed control method. An optimal balance between the libration stability and a fast deorbit of satellite with minimum control efforts is achieved.     

Coordinated control of tethered space robot using mobile tether attachment point in approaching phase

Tethered space robots (TSRs) have wide applications in future on-orbit service owing to its flexibility and great workspace. However, the control problem is quite complex and difficult in the phase of approaching target, and the fuel consumption must also be taken into account. Hence, we present a novel scheme of achieving coordinated orbit and attitude control simultaneously for the TSR. Space tether, which can provide greater force compared with the thruster force, is used in the design of the orbit and attitude coordinated controller. A coordinated control mechanism is designed to provide attitude control torques of the pitch and yaw motions by adjusting the position of the mobile tether attachment point, while the roll motion is stabilized by the thruster. In order to guarantee this mechanism to work properly, constant tether tension strategies are utilized to plan an optimal approaching trajectory which is tracked by the coordinated controller of tether force and thruster force. Numerical simulation validates the feasibility of our proposed coordinated control scheme for TSR in the approaching phase. Furthermore, fuel consumption of the orbit and attitude control are both significantly reduced compared with traditional thruster control.     

Spin deployment of Hub-Spoke tethered satellite formation with sliding mode tether tension control

This work develops a tension control strategy for deploying an underactuated spin-stable tethered satellite formation in the hub-spoke configuration. First, the Lagrange equation is used to model the spin-deployment dynamics of the tethered satellite formation. The central spacecraft is modeled as a rigid body, and the tethered subsatellites are simplified as lumped masses. Second, a pure tension controller has been proposed to suppress the tether libration motion in the deployment without thrusting at the subsatellites. A nonlinear sliding mode control is introduced in the tension controller for the underactuated system to suppress the periodic gravitational perturbations caused by the spinning hub-spoke tethered satellite formation. The unknown upper bounds of the perturbations are estimated by adaptive control law. The bounded stability of the closed-loop tension controller has been proved by the Lyapunov theory. Finally, numerical simulations validate the effectiveness and robustness of the proposed controller, i.e., tethers are fully deployed stably to the desired hub-spoke configuration.     

Estimation of flexible space tether state based on end measurement by finite element Kalman filter state estimator

This paper proposes a novel finite element Kalman filter to estimate the unmeasurable state of space tether systems based on the measured state at its ends only. The finite element method calculates the unmeasurable internal state as the virtual measurement based on the dynamic model of the system by imposing the input of measured state at the boundary to the model using the Lagrange multiplier method in the spatial space. Combining the real and virtual measurement into a hybrid measurement model of the system, the full state is reconstructed and propagated in the temporal space by the extended Kalman filter. Two state-space system models, the dynamics-based and kinematics-based state models, in the Kalman filter are explored. The observability and stability of the newly proposed finite element Kalman filter are examined and proved. The advantages of the proposed state estimator are (i) the singularity in the virtual measurement of state caused by the number of internal state greater than the number of state measured at the boundary is eliminated in the statistic meaning by the Kalman filter, and (ii) the effects of noises of the observation data and the uncertainties of model discretization are considered and minimized. The correctness and effectiveness of the proposed state estimator is demonstrated by the numerical analysis of a space tether system orbiting around the Earth. The results show the proposed state estimator with only measured state at the ends of the tether successfully provides an accurate time history estimation of geometric configuration and motion of the entire tether. Moreover, the results also show the difference caused by the dynamics-based and kinematics-based system models in the state estimator is negligible. The kinematics-based system model should be used in the state estimator due to its significantly low computational load. Finally, the proposed method can be easily applied for the state estimation process for other space tethered spacecraft systems.     

Shape and motion estimate of LEO debris using light curves

We have derived a tri-axial ellipsoidal model of an LEO object, a Cosmos 2082 rocket body, including its rotational axis direction, rotation period, precession, and a compositional parameter, using only light curve data from an optical telescope. The brightness of the object was monitored for two days and least-squares fitting was used to determine these values. The derived axial ratios are 100:18:18, the coordinates of the rotational axis direction on the celestial sphere are R.A. = 305.8° and Dec. = 2.6°, and its observed average rotation period is 41 s. When precession is considered, its amplitude and precession period are 30.5° and 29.4 min. These results show that optical light curve data are sufficient to determine the shapes and the motions of LEO objects.     

Estimation of orbital parameters of broken-up objects from in-situ debris measurements

Even sub-millimeter-size debris could cause a fatal damage on a spacecraft. Such tiny debris cannot be followed up or tracked from the ground. Therefore, Kyushu University has initiated IDEA the project for In-situ Debris Environmental Awareness, which conducts in-situ measurements of sub-millimeter-size debris. One of the objectives is to estimate the location of on-orbit satellite fragmentations from in-situ measurements. The previous studies revealed that it is important to find out the right nodal precession rate to estimate the orbital parameters of a broken-up object properly. Therefore, this study derives a constraint equation that applies to the nodal precession rate of the broken-up object. This study also establishes an effective procedure to estimate properly the orbital parameters of a broken-up object with the constraint equation.     

Protecting Earth-orbiting spacecraft against micro-meteoroid/orbital debris impact damage using composite structural systems and materials: An overview

Spacecraft that are launched to operate in Earth orbit are susceptible to impacts by meteoroids and pieces of orbital debris (MMOD). The effect of a MMOD particle impact on a spacecraft depends on where the impact occurs, the size, composition, and speed of the impacting object, the function of the impacted system. In order to perform a risk analysis for a particular spacecraft under a specific mission profile, it is important to know whether or not the impacting particle (or its remnants) will exit the rear of an impacted spacecraft wall. A variety of different ballistic limit equations (BLEs) have been developed for many different types of structural wall configurations. BLEs can be used to optimize the design of spacecraft wall parameters so that the resulting configuration is able to withstand the anticipated variety of on-orbit high-speed impact scenarios. While the level of effort exerted in studying the response of metallic multi-wall systems to high speed particle impact is quite substantial, the extent of the effort to study composite material and composite structural systems under similar impact conditions has been much more limited. This paper presents an overview of the activities performed to assess the resiliency of composite structures and materials under high speed projectile impact. The activities reviewed will be those that have been aimed at increasing the level of protection afforded to spacecraft operating in the MMOD environment, and more specifically, on those activities performed to mitigate the mechanical and structural effects of an MMOD impact.     

Faster algorithm of debris cloud orbital character from spacecraft collision breakup

Space debris is polluting the space environment. Collision fragment is its important source. NASA standard breakup model, including size distributions, area-to-mass distributions, and delta velocity distributions, is a statistic experimental model used widely. The general algorithm based on the model is introduced. But this algorithm is difficult when debris quantity is more than hundreds or thousands. So a new faster algorithm for calculating debris cloud orbital lifetime and character from spacecraft collision breakup is presented first. For validating the faster algorithm, USA 193 satellite breakup event is simulated and compared with general algorithm. Contrast result indicates that calculation speed and efficiency of faster algorithm is very good. When debris size is in 0.01–0.05 m, the faster algorithm is almost a hundred times faster than general algorithm. And at the same time, its calculation precision is held well. The difference between corresponding orbital debris ratios from two algorithms is less than 1% generally.     

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.     

Stabilization of an injected conducting layer for artificially enhancing drag on orbital debris

The evolution of a magnetized conducting medium suspended in magnetic and gravitational fields is examined. In this paper some effects of the influence of velocity fields on the linear stability properties of such layers are investigated. A fully compressible, three-dimensional analysis of the layer is described. The relevant equations are derived and then solved by the MagnetoHydroDynamic SPEctral Compressible Linear Stability (MHDSPECLS) algorithm, a Chebyshev collocation code. The code allows for the computation of magnetic and thermal effects. A complete stabilization of the system is found above a critical velocity of approximately 2500 m/s.     

Assessment of space debris collisions against spacecraft with deorbit devices

Deorbit methods have been employed to remove space debris from orbit. One of these methods is to utilize atmospheric drag. In this method, a membrane loaded into the spacecraft is expanded to increase atmospheric drag. Although this method works without requiring fuel, it has the disadvantage of a high risk of collision with other debris owing to its larger area. Area-time product and energy-to-mass ratio have been used as indices to evaluate the risk of collisions between spacecraft and debris. However, the evaluation criteria were uncertain because these two indices are independent. In this paper, we propose a new evaluation index, single-sheet collision factor (SSCF), that comprehensively evaluates the collision risk based on experiments simulating debris collisions. As a result of the hypervelocity collision experiment, we found that the penetration-area mass of the spacecraft affects the severity of debris collisions. In this paper, the product of the exterior-wall thickness, the exterior-wall density, and the space debris cross-sectional area defines the penetration-area mass of the spacecraft. Furthermore, we compare and evaluate various deorbit methods using SSCF. The comparison showed that the penetration-area mass of the SSCF could be quantitatively determined for the debris-collision severity due to difference in structural materials of spacecraft. SSCF will be used to create rules for space-environment conservation with the expansion of the space-development market.     

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.     

Study of the radiation fields in LEO with the Space Application of Timepix Radiation Monitor (SATRAM)

We present the analysis of data taken by the Space Application of Timepix Radiation Monitor (SATRAM). It is centred on a Timepix detector (300?$\mathrm{\mu }$m thick silicon sensor, pixel pitch 55?$\mathrm{\mu }$m, 256?×?256 pixels). It was flown on Proba-V, an Earth observing satellite of the European Space Agency (ESA) from an altitude of 820?km on a sun-synchronous orbit, launched on May 7, 2013. A Monte Carlo simulation was conducted to determine the detector response to electrons (0.5–7?MeV) and protons (10–400?MeV) in an omnidirectional field taking into account the shielding of the detector housing and the satellite. With the help of the simulation, a strategy was developed to separate electrons, protons and ions in the data. The measured dose rate and stopping power distribution are presented as well as SATRAM’s capability to measure some of the stronger events in Earth’s magnetosphere. The stopping power, the cluster height and the shape of the particle tracks in the sensor were used to separate electrons, protons and ions. The results are presented as well. Finally, the pitch angles for a short period of time were extracted from the data and corrected with the angular response determined by the simulation.     

The bistatic radar capabilities of the Medicina radiotelescopes in space debris detection and tracking

An accurate measurement of the position and trajectory of the space debris fragments is of primary importance for the characterization of the orbital debris environment. The Medicina Radioastronomical Station is a radio observation facility that is here proposed as receiving part of a ground-based space surveillance system for detecting and tracking space debris at different orbital regions (from Low Earth Orbits up to Geostationary Earth Orbits). The proposed system consists of two bistatic radars formed by the existing Medicina receiving antennas coupled with appropriate transmitters. This paper focuses on the current features and future technical development of the receiving part of the observational setup. Outlines of possible transmitting systems will also be given together with the evaluation of the observation strategies achievable with the proposed facilities.     

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.     

Feasibility study of LEO,GEO and Molniya orbit based satellite solar power station for some identified sites in India

The analysis of satellite solar power station (SSPS) is carried out for some specified locations (Delhi, Mumbai, Kolkata and Bengaluru) in India and consequently the performance of the system is evaluated for geostationary earth orbit (GEO) based SSPS, low earth orbit (LEO) based SSPS and Molniya (quasi geostationary) orbit based SSPS for sites located at different latitudes. The analysis of power; received energy over a year and weight of the rectenna array for the same beam intensity showed varied results for Molniya orbit based SSPS, LEO based SSPS and GEO based SSPS. The power delivered by the LEO SSPS was highest which indicated that this SSPS may be efficient for the short term power requirement. However, it is observed from the results of the energy received over a year that the GEO based system is suitable for base load power plant as it is capable of delivering constant energy through out a year. Further, the weight of the rectenna and hence the space required for ground station for same power output is smallest for Molniya orbit based system for a range of rectenna array radius considered. It is thus envisaged that the Molniya orbit based system would be a better choice for commercial use of SSPS. These findings may help for judicious selection of satellite orbit and ground station location for placing the satellite for SSPS for various applications.     

Different space weather effects in anomalies of the high and low orbital satellites

Preliminary results of the EU INTAS Project 00810, which aims to improve the methods of safeguarding satellites in the Earth’s magnetosphere from the negative effects of the space environment, are presented. Anomaly data from the “Kosmos” series satellites in the period 1971–1999 are combined in one database, together with similar information on other spacecraft. This database contains, beyond the anomaly information, various characteristics of the space weather: geomagnetic activity indices (Ap, AE and Dst), fluxes and fluences of electrons and protons at different energies, high energy cosmic ray variations and other solar, interplanetary and solar wind data. A comparative analysis of the distribution of each of these parameters relative to satellite anomalies was carried out for the total number of anomalies (about 6000 events), and separately for high (5000 events) and low (about 800 events) altitude orbit satellites. No relation was found between low and high altitude satellite anomalies. Daily numbers of satellite anomalies, averaged by a superposed epoch method around sudden storm commencements and proton event onsets for high (>1500 km) and low (<1500 km) altitude orbits revealed a big difference in a behavior. Satellites were divided on several groups according to the orbital characteristics (altitude and inclination). The relation of satellite anomalies to the environmental parameters was found to be different for various orbits that should be taken into account under developing of the anomaly frequency models.     

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.     

Analysis and design of Cubesat constellation for the Mediterranean south costal monitoring against illegal immigration

Costal monitoring is focused on fast response to illegal immigration and illegal ship traffic. Especially, the illegal ship traffic has been present in media since April 2015, as the number of reported deaths of immigrants crossing the Mediterranean significantly increased. Satellite images provide a possibility to at least partially control both types of events. This paper defines the principal criteria to select the best satellite constellation architecture for maritime and coastal monitoring, filling the gaps of imagery techniques in term of real-time control. The primary purpose of a constellation is to obtain global measurement improving the temporal resolution. The small size and low-cost are the main factors, which make CubeSats ideal for use in constellations. We propose a constellation of 9 Cubesats distributed evenly in 3 different planes. This reduces the revisit time enhancing the coverage duration. In addition, it also allows observing fire, damage on building and similar disasters. In this analysis, the performance criteria were reported such as the revisit time, the vision duration and the area coverage.