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.     

Validation of non-linear split window algorithm for land surface temperature estimation using Sentinel-3 satellite imagery: Case study; Tehran Province,Iran

In recent years, land surface temperature (LST) has become critical in environmental studies and earth science. Remote sensing technology enables spatiotemporal monitoring of this parameter on large scales. This parameter can be estimated by satellite images with at least one thermal band. Sentinel-3 SLSTR data provide LST products with a spatial resolution of 1 km. In this research, direct and indirect validation procedures were employed to evaluate the Sentinel-3 SLSTR LST products over the study area in different seasons from 2018 to 2019. The validation method was based on the absolute (direct) evaluation of this product with field data and comparison (indirect) evaluation with the MODIS LST product and the estimated LST using the non-linear split-window (NSW) algorithm. Also, two emissivity estimation methods, (1) NDVI thresholding method (NDVI-THM) and (2) classification-based emissivity method (CBEM), were used to estimate the LST using the NSW method according to the two thermal bands of Sentinel-3 images. Then, the accuracy of these methods in estimating LST was evaluated using field data and temporal changes of vegetation, which the NDVI-THM method generated better results. For indirect evaluation between the Sentinel-3 LST product, MODIS LST product, and LST estimated using NSW, four filters based on spatial and temporal separates between pairs of pixels and pixel quality were used to ensure the accuracy and consistency of the compared pairs of a pixel. In general, the accuracy results of the LST products of MODIS and Sentinel-3, and LST estimated using NSW showed a similar trend for LST changes during the seasons. With respect to the two absolute and comparative validations for the Sentinel-3 LST products, summer with the highest values of bias (?1.24 K), standard deviation (StDv = 2.66 K), and RMSE (2.43 K), and winter with the lowest ones (bias of 0.14 K, StDv of 1.13 K, and RMSE of 1.12 K) provided the worst and best results for the seasons in the period of 2018–2019, respectively. According to both absolute and comparative evaluation results, the Sentinel-3 SLSTR LST products provided reliable results for all seasons on a large temporal and spatial scale over our studied area.     

GPS phase difference variation statistics: A comparison between phase scintillation index and proxy indices

Scintillated GPS phase observations are traditionally characterized by the phase scintillation index, derived from specialized GPS receivers usually tracking at 50 Hz. Geodetic quality GPS receivers, on the other hand, are normally tracking at frequencies up to 1 Hz. However, availability of continuously operating geodetic receivers both in time and geographical location are superior to scintillation receiver’s coverage in many parts of the world. This motivates scintillation studies using regional and global geodetic GPS networks. Previous studies have shown the usefulness of GPS estimated total electron content variations for detecting ionospheric irregularities. In this paper, collocated geodetic and scintillation receivers are employed to compare proxy indices derived from geodetic receivers with the phase scintillation index during quiet and moderately disturbed ionospheric conditions. Sensitivity of the phase scintillation indices at high latitude stations to geomagnetic activity is discussed. Global mapping of ionospheric disturbances using proxy indices from real-time 1 Hz GPS stations are also presented.     

Aircraft robust multidisciplinary design optimization methodology based on fuzzy preference function

This paper presents a Fuzzy Preference Function-based Robust Multidisciplinary Design Optimization (FPF-RMDO) methodology. This method is an effective approach to multidisciplinary systems, which can be used to designer experiences during the design optimization process by fuzzy preference functions. In this study, two optimizations are done for Predator MQ-1 Unmanned Aerial Vehicle (UAV): (A) deterministic optimization and (B) robust optimization. In both problems, minimization of takeoff weight and drag is considered as objective functions, which have been optimized using Non-dominated Sorting Genetic Algorithm (NSGA). In the robust design optimization, cruise altitude and velocity are considered as uncertainties that are modeled by the Monte Carlo Simulation (MCS) method. Aerodynamics, stability and control, mass properties, performance, and center of gravity are used for multidisciplinary analysis. Robust design optimization results show 46% and 42% robustness improvement for takeoff weight and cruise drag relative to optimal design respectively.     

Temporal extrapolation of TEC using WNN during 2007–2018 and comparison against GIM,IRI2016 and Kriging

In this paper, a new method of temporal extrapolation of the ionosphere total electron content (TEC) is proposed. Using 3-layer wavelet neural networks (WNNs) and particle swarm optimization (PSO) training algorithm, TEC time series are modeled. The TEC temporal variations for next times are extrapolated with the help of training model. To evaluate the proposed model, observations of Tehran GNSS station (35.69°N, 51.33°E) from 2007 to 2018 are used. The efficiency of the proposed model has been evaluated in both low and high solar activity periods. All observations of the 2015 and 2018 have been removed from the training step to test the proposed model. On the other hand, observations of these 2 years are not used in network training. According to the F10.7, the 2015 has high solar activity and the 2018 has quiet conditions. The results of the proposed model are compared with the global ionosphere maps (GIMs) as a traditional ionosphere model, international reference ionosphere 2016 (IRI2016), Kriging and artificial neural network (ANN) models. The root mean square error (RMSE), bias, dVTEC = |VTEC_GPS ? VTEC_Model| and correlation coefficient are used to assess the accuracy of the proposed method. Also, for more accurate evaluation, a single-frequency precise point positioning (PPP) approach is used. According to the results of 2015, the maximum values of the RMSE for the WNN, ANN, Kriging, GIM and IRI2016 models are 5.49, 6.02, 6.34, 6.19 and 13.60 TECU, respectively. Also, the maximum values of the RMSE at 2018 for the WNN, ANN, Kriging, GIM and IRI2016 models are 2.47, 2.49, 2.50, 4.36 and 6.01 TECU, respectively. Comparing the results of the bias and correlation coefficient shows the higher accuracy of the proposed model in quiet and severe solar activity periods. The PPP analysis with the WNN model also shows an improvement of 1 to 12 mm in coordinate components. The results of the analyzes of this paper show that the WNN is a reliable, accurate and fast model for predicting the behavior of the ionosphere in different solar conditions.     

Flow visualization study of an airship model using a water towing tank

The work presented investigates the aerodynamic drag and flow physics of an airship configuration. The investigation is carried out experimentally in a water towing tank. The purpose of this project is to improve the understanding of the aerodynamic characteristics of an airship through quantitative and qualitative approaches. Qualitative measurements of integral drag force and moment are measured in the presence of a crosswind. Flow visualization using colored and fluorescent dyes has been applied to improve knowledge of the aerodynamics of airship configurations with an emphasis on the boundary layer separation, transition and hull–appendage interference effects.     

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.     

Lineament mapping and fractal analysis using SPOT-ASTER satellite imagery for evaluating the severity of slope weathering process

The present study describes a remote sensing approach for preparing lineament map that subsequently indicates the influence of lineament density in the severity of weathering development. In this study, SPOT-5 data, the integration of SPOT-ASTER and Digital Elevation Model (DEM) data were used and processed. The existence of an active fault system in the south of Mashhad city, NE Iran and presence of schistose rocks in this area result in the development of numerous lineament features. This region was selected for this research. Lineament features including fractures, bedding plane, cleavage, shear zones and schistosity were mapped in the study area. The results indicate that the highest concentration of lineaments occurred in the central-western and south-eastern parts of the study area, which coincide with metamorphic outcrops and NW-SE trending fault system. A comparison of lineament statistical analysis and field survey demonstrated that the structural discontinuities have a significant effect on forming and distribution of weathering profiles. It was observed that increasing the number, length and density of structural discontinuities led to strong severity in weathering, which can produce deep residual soils susceptible to landslide occurrence. The remote sensing approach developed in this study can be applicable for preparing lineament maps and evaluating the severity of weathering development in other active fault zones around the world.     

Experimental and numerical investigation of threedimensional vortex structures of a pitching airfoil at a transitional Reynolds number

This research examines the vortex behaviors and aerodynamic forces in dynamic stall phenomena at a transitional Reynolds number(Re = 90000) using experimental and numerical approaches.Periodic sinusoidal pitching motion at two different reduced frequencies is used to achieve the dynamic stall of a NACA 0012 airfoil.Several leading edge vortices form and detach in the dynamic stall stage.The flow then quickly transitions to a full separation zone in the stall stage when the angle of attack starts to decrease.There is discrepancy between the phaseaveraged and instantaneous flow field in that the small flow structures increased with angle of attack, which is a characteristic of the flow field at the transitional Reynolds number.The interaction between the streamwise vortices in the three-dimensional numerical results and the leading edge vortex are the main contribution to the turbulent flow.In addition, the leading edge vortex that supplies vortex lift is more stable at higher reduced frequency, which decreases the lift fluctuation in the dynamic stall stage.The leading edge vortex at higher reduced frequency is strong enough to stabilize the flow, even when the airfoil is in the down-stroke phase.