Assessment of the close approach frequency and collision probability for satellites in different configurations of large constellations

The Earth orbital environment is drastically changing due to an intensification of the space activities. In particular, several projects of large constellations, proposed for the next years for communications purpose like global internet access, Internet of Things, or for Earth observations, will lead to the deployment of several thousands of new satellites at an unprecedented rate. It is a crucial challenge for space traffic management, which will deal with a great number of satellite conjunctions, potentially causing a collision with damaging consequences for the constellation itself and the space environment sustainability.In this paper, we investigate the close approach frequency and the cumulative collision probability for each referenced constellation. For this purpose, we compute the orbital evolution of satellites in different constellations during the lifecycle, from the deployment to the decommissioning, and we apply the CUBE algorithm and the Foster method to assess the collision probability with the background space debris population assuming a constant uncertainty in position. We show the variation of risk defined by the close approach frequency and the cumulative collision probability as a function of the proposed configuration. In particular, satellites of the Iridium and Kuiper constellation, but also satellite of the Telesat constellation on polar orbits are the most exposed at a collision. Moreover, the decommissioning phase contribute for a major part to the final cumulative collision probability.     

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.     

Development and performance evaluation of tropospheric scintillation model on Ku-band satellite link over Akure,Nigeria

Existing amplitude scintillation prediction models often perform less satisfactorily when deployed outside the regions where they were formulated. This necessitates the need to evaluate the performance of scintillation models developed in one region using data data from other regions while documenting their relative errors. Due to its variation with elevation angle, frequency, other link parameters and meteorological factors, we employed three years (January 2016 to December 2018) of concurrently measured satellite radio beacons and tropospheric weather parameters to develop a location-based amplitude scintillation prediction model over the Earth-space path of Akure (7.17^oN, 5.18^oE), South-western Nigeria. The satellite beacon measurement used Tektronix Y400 NetTek Analyzer at 1 s integration time while meteorological parameters, namely; temperature, pressure and relative humidity were measured using Davis Vantage Vue weather station at 1 min integration time. Comparative study of the model’s performance with nine (9) existing scintillation prediction models indicates that the best and worst performing models, in terms of root mean square error (RMSE), are the Statistical Temperature and Refractivity (STN) and direct physical and statistical prediction (DPSP) models with values 11.48 and 51.03 respectively. Also, worst month analysis indicates that April, with respective enhancement and fade values of 0.88 and 0.90 dB for 0.01% exceedance, is the overall worst calendar month for amplitude scintillation.     

A multi-stage regulation strategy of space manipulators with a free-swinging joint failure

To ensure tasks can be completed after a free-swinging joint failure occurs, a multi-stage regulation strategy of space manipulators is proposed. Considering all terms of the dynamics equation, an evaluation model of the regulation ability (EMRA) of active joints over the fault joint is established based on the fuzzy entropy. And then a multi-stage regulation strategy based on the EMRA is designed to regulate the fault joint. The strategy divides the regulation process into several stages, and select a certain active joint to regulative the fault joint in every stage. With this multi-stage regulation strategy, the fault joint can be regulated to the desired angle without huge torque on regulative joints. The simulation is carried out with a 7-DOF space manipulator, verifying the correctness and effectiveness of the multi-stage regulation strategy. The strategy has three advantages: Coriolis and centrifugal terms are both considered for the first time in selecting the regulative joint, making the selection result more in line with the actual regulation process; The influence of the model uncertainty is eliminated in establishing the EMRA, making the evaluation of regulative ability more precise; The fault joint is successfully regulated to the desired angle without huge torque on regulative joints.     

Excavation of artificial caverns inside asteroids by leveraging rotational self-energy

Future space ventures will likely require exploitation of near-Earth asteroid resources. Moreover, it can be envisaged that asteroids may host habitats in their interiors. In fact, a cavern inside an asteroid would be a natural radiation shield against cosmic radiation and may also serve as a confined environment for storage of mined material such as water ice or other processed volatiles such as propellants. To this end, this paper proposes to leverage the asteroid rotational self-energy to remove material from the asteroid interiors and create a spherical cavern, by means of the orbital siphon concept. The siphon is a chain of tether-connected payload masses (the asteroid material), which exploits the rotation of the asteroid for the delivery of mass from the asteroid to escape. Under certain conditions the siphon can be initiated to ensure self-sustained flow of mass from the asteroid to escape. A net orbital siphon effect is generated by connecting new payloads at the bottom of the chain while releasing the upper payloads. Key parameters are discussed, such as the required siphon dimension and the maximum size of the internal cavity that can be excavated, as a function of the asteroid rotational period. Moreover, assuming elastic material behaviour, a closed-form expression for the stress tensor is found and a failure criterion is used to identify regions in the asteroid interiors subjected to the larger stresses. It is shown that the conditions for failure are relaxed as the radius of the internal void increases.     

Estimating the spin axis orientation of the Echostar-2 box-wing geosynchronous satellite

For the first time, the spin axis orientation of an inactive box-wing geosynchronous satellite has been estimated from ground-based optical photometric observations of Echostar-2’s specular reflections. Recent photometric light curves obtained of Echostar-2 over four years suggest that unusually bright and brief specular reflections were occurring twice within an observed spin period. These bright and brief specular reflections suggested two satellite surfaces with surface normals separated by approximately 180°. The geometry between the satellite, the Sun, and the observing location at the time of each of the brightest observed reflections, was used to estimate Echostar-2’s equatorial spin axis orientation coordinates. When considering prograde and retrograde rotation, Echostar-2’s spin axis orientation was estimated to have been located within 30° of either equatorial coordinate pole. Echostar-2’s spin axis was observed to have moved approximately 180° in right ascension, within a time span of six months, suggesting a roughly one year spin axis precession period about the satellite’s angular momentum vector.     

Impact of geomagnetic variation over sub-auroral ionospheric region during high solar activity year 2014

The present work is an attempt to evaluate the impact of changing space weather condition over sub-auroral ionosphere during high solar activity year 2014. In view of this, the GPS based TEC along with Ionosonde data over Indian permanent scientific base “Maitri”, Antarctica (70°46′00″S, 11°43′56″E) has been utilized. The results suggested that the nature of ionospheric responses to the geomagnetic disturbances not only depended upon the status of high latitudinal electro-dynamic processes but also influenced by the seasonal variations. The results revel both negative and positive type of ionospheric response in a single year but during different seasons. The study suggested that the combination of equator-ward plasma transportation along with ionospheric compositional changes causes a negative ionospheric impact especially during summer and equinox seasons. However, the combination of pole-ward contraction of the oval region along with particle precipitation may lead to exhibit positive ionospheric response during the winter season. The plasma transportation direction has been validated with the help of convection boundary (HM boundary) deduced with the help of SuperDARN observations. The ground based ionosonde observations clearly provided the evidence of deep penetration of high energetic particles up to the E-layer heights which results a sudden and strong appearance of E-layer. The strengthening of E-layer is responsible for modification of auroral electrojet and field-aligned current system. Also, the sudden appearance of E-layer along with a decrease in F-layer electron density suggested the dominance of NO⁺ over O⁺ in a considered region under geomagnetic disturbed condition.     

Packaging of thick membranes using a multi-spiral folding approach: Flat and curved surfaces

Elucidating versatile configurations of spiral folding, and investigating the deployment performance is of relevant interest to extend the applicability of deployable membranes towards large-scale and functional configurations.In this paper we propose new schemes to package flat and curved membranes of finite thickness by using multiple spirals, whose governing equations render folding lines by juxtaposing spirals and by accommodating membrane thickness. Our experiments using a set of topologically distinct flat and curved membranes deployed by tensile forces applied in the radial and circumferential directions have shown that (1) the multi-spiral approach with prismatic folding lines offered the improved deployment performance, and (2) the deployment of curved surfaces progresses rapidly within a finite load domain. Furthermore, we confirmed the high efficiency of membranes folded by multi-spiral patterns.From viewpoints of configuration and deployment performance, the multi-spiral approach is potential to extend the versatility and maneuverability of spiral folding mechanisms.     

双模块下颌式内转进气道/圆锥前体一体化布局研究

为实现三维内转进气道的内收缩流场与圆锥前体的外压缩流场的良好匹配，提出了一种双模块下颌式内转进气道/圆锥前体(Double-modules Chin Inward-turning Inlet and Conical Forebody，DCII/CF)一体化设计方法，获得一种新颖的双发并置、侧向安装的DCII/CF一体化布局。针对该布局形式，开展了DCII/CF一体化构型与传统的单模块内转进气道/圆锥前体(Single-module Inward-turning Inlet and Conical Forebody，SII/CF)一体化构型的数值对比研究。结果表明：DCII/CF一体化布局不仅为内转进气道提供了优秀的前体附面层排移效果，还有效避免了传统SII/CF布局中前体附面层与进气道内部流场之间的相互干扰。在Ma∞=6.0设计状态，DCII/CF一体化布局的进气道总压恢复系数相较传统的SII/CF布局有了显著提高，从0.403提高至0.482；但由于前体附面层的排移，该布局的捕获流量略有降低， SII/CF的流量系数为0.956，该布局则为0.917。而在非设计状态，该布局形式同样具备较好的总压恢复性能，在Ma∞=5.0与Ma∞=4.0的总压恢复系数分别达到了0.586和0.682，明显高于SII/CF的总压恢复系数0.507和0.619。     

Substitution of satellite-based land surface temperature defective data using GSP method

Land surface temperature (LST) as an important environmental variable provides valuable information for earth environmental system modelling. Currently, LST is obtained through satellite thermal sensors at various spatial and temporal resolutions. Although spatially continuous satellite-based LST measurements are intended to overcome the shortcomings of sparse ground-based LST measurements, LST images often contain anomalous values due to the existence of clouds or sensor malfunctioning. The problem becomes more serious where the users deal with high spatial resolution characterized by low temporal resolution. This study examines the capability of a newly developed graph signal processing (GSP) method using two-dimensional single-date thermal data. For this purpose, four Landsat/TIRS datasets are analyzed. The data of five elliptical regions on thermal images are eliminated and then reconstructed through the GSP method and using the LST values of the enclosing rectangles containing the ellipsoids. The results indicate that the temperature variation determined by the GSP method generally conforms to the original image LST values. According to a correlation test conducted on the original image LST and those obtained through the GSP method, the values vary from 58% to 95%, which is an above-the-average rate (RMSE from 0.69 to 2.27). The statistical analysis of the original image LST in both the elliptical regions and the enclosing rectangles containing the ellipsoids indicates that an increase in the variance of LST data causes an increased error in the calculation of temperature by the GSP method, and vice versa. The results of the analysis of variance (ANOVA) and Duncan test indicated that an increase in the number of the non-zero spectral bins would result in increased RMSE values for all the dates and the regions. Moreover, the model errors were significant at the 0.05 level across all the image date and five elliptical study regions. Based on the results, the use of this method is recommended for the reconstruction of LST missing values, where dissimilarity of atmospheric conditions limits the use of other methods that depend on the time series data of various dates and a great deal of data calculation.