2021 Vol. 41, No. 6

Display Method:
中国首颗太阳探测科学技术试验卫星羲和号成功发射
2021, 41(6): 845-845.
Abstract:
NEOWISE任务将延期运行至2023年6月
2021, 41(6): 845-845.
Abstract:
NASA选定VIPER月球漫游车着陆点
2021, 41(6): 846-846.
Abstract:
ESA启动Moonlight月球卫星星座建设
2021, 41(6): 846-846.
Abstract:
中国成功发射可持续发展科学卫星
2021, 41(6): 846-846.
Abstract:
20亿年前玄武岩进一步揭示月球演化奥秘
2021, 41(6): 847-847.
Abstract:
国际空间站安全委员会建议延长空间站运行期限
2021, 41(6): 848-848.
Abstract:
NASA资助10项深空探测技术研发青年项目
2021, 41(6): 848-848.
Abstract:
Research Articles
Statistical Analysis on Four Kinds of Space Weather Events during Solar Cycles 21~24
MIAO Juan, LI Zhitao, REN Tingling, WANG Xin
2021, 41(6): 849-857. doi: 10.11728/cjss2021.06.849
Abstract:
Based on the statistics of different intensity levels of solar X-ray flares, solar proton events, geomagnetic storms and relativistic electron flux enhancement of Solar Cycle 21~24, the following conclusions can be drawn. Firstly, the total number of X-ray flares in solar cycle is directly proportional to the number of sunspot maximum, and the correlation coefficients of the total number of flares, X-class flares and sunspot maximum are 0.974 and 0.997, respectively. Secondly, the solar proton events mainly occurred in the first two years before and after the peak year, accounting for 80% of the total number of events. 83% solar proton events with peak flux over 10 pfu (1 pfu=1 cm-2·sr-1·s-1) were mainly accompanied by X-class and M-class flares, but also with a small number of M-class flares. For the solar proton events with the peak flux over 1000 pfu, about 98% of them were accompanied by C-class flares. Thirdly, the most frequent time of geomagnetic storms occurred during Cycle 21, 22, 23 and 24 was in 1982, 1991, 2003 and 2015, respectively, which lagged behind the peak time of sunspot maximum by three years, two years, two years and one year, respectively. Finally, 72% of the relativistic electron flux enhancement occurred in the declining phase of solar cycle, and 24% of them occurred in the rising phase.
Multi-case Study of Current Sheet Flapping Motions Induced by Non-adiabatic Ions
WEI Xinhua, CAI Chunlin
2021, 41(6): 858-868. doi: 10.11728/cjss2021.06.858
Abstract:
In this paper, the y-component of magnetic field line curvature in the plasma sheet was analyzed, and two kinds of shear structures of the flapping current sheet were found, i.e. symmetric and antisymmetric. The alternating bending orientations of the guiding field are exactly corresponding to alternating north-south asymmetries of the bouncing ion population in the sheet center. Those alternating asymmetric plasma sources consequently induce the current sheet flapping motion as a driver. In addition, a substantial particle population with downward motion was observed in the center of a bifurcated current sheet. This population is identified as the quasi-adiabatic particles, and provides a net current opposite to the conventional cross-tail current.
Influence of Continuous Magnetic Activities on the Evolution of the Plasmasphere
TI Shuo, SHEN Chao, CHEN Tao, JI Yong, HE Han, XU Ronglan, HUANG Ya
2021, 41(6): 869-880. doi: 10.11728/cjss2021.06.869
Abstract:
In this paper, the evolution of the plasmasphere's morning side viewing was analyzed under a series of magnetic activities using the extreme ultraviolet logarithmic data observed by the Chang'E-3 extreme ultraviolet camera on 21 February 2014. Based on the phase space distribution of the protons, the evolution of the plasmasphere on the magnetic equatorial plane was simulated when a series of magnetic storm events occurred from 18 to 22 February 2014. Through observation and simulation, it is found that the actual filling speed of the plasmasphere is greater than the filling speed of the plasmasphere during the simulation. It is presumed that the plume structure between the dusk side and the sun side has an affect on the position of the plasmapause on the sun side from side view. In simulation, the overall response of the plasmasphere to the magnetic storm is within 3 hours, but large magnetic storm has a long-term impact on the plasmasphere, which can reach 1~2 days. Continuous magnetic storm events have a significant impact on the plasmasphere, there is a superimposed effect. Refilling of the plasmasphere takes longer than erosion.
Transmission Characteristics of Electromagnetic Waves in Plasma with External Orthogonal Electromagnetic Fields
ZHANG Jie, ZHAO Shanchao, ZHANG Guodong
2021, 41(6): 881-886. doi: 10.11728/cjss2021.06.881
Abstract:
On the basis of Magnetohydrodynamics (MHD) and the theory of electromagnetic wave propagation, a physical model is proposed to solve the blackout in the reentry process of spacecraft. By means of numerical analysis, the variation trend of the electron density around the electromagnetic window is studied with the increase of the normal distance of the aircraft model under the premise of different axial distances of the aircraft model; the variation trend of electron density around the electromagnetic window is analyzed by changing the cross angle of the applied electric field and magnetic field; the variation trend of the ratio of electron density after applying an orthogonal electromagnetic field to electron density without applying electromagnetic field and the variation trend of attenuation of electromagnetic wave. The results show that when the orthogonal electromagnetic field is applied, the electron density around the antenna decreases with the increase of the normal distance of the aircraft, and the ratio of electron density to that without the applied electromagnetic field shows a decreasing trend. Attenuation of the electromagnetic wave decreases with the increase of electromagnetic frequency at different flying heights and magnetic field strength of the applied magnetic field. This provides a new method to alleviate the blackout.
Global Ionospheric TEC and ROTI Variations during a Moderate Geomagnetic Storm
OU Ming, WU Jiayan, CHEN Longjiang, ZHEN Weimin
2021, 41(6): 887-897. doi: 10.11728/cjss2021.06.887
Abstract:
The ionosphere will have strong disturbances on geomagnetic storm days. The physical mechanism causing an ionospheric disturbance is very complex, various dynamics and electrodynamics processes in magnetosphere/thermosphere lead to great differences in ionospheric morphology and response in different locations and different time during a geomagnetic storm. Therefore, the geomagnetic storm has always been a hot and difficult issue in ionospheric research. In August 2014, under the influence of a solar flare explosion and two CME events, a strong magnetic storm was triggered on 27 August, accompanied by an ionospheric storm effect. The maximum Kp index reached 4.7, and the lowest Dst index reached -80 nT, reaching the level of moderate intensity geomagnetic storm. Based on the data of global ionospheric TEC and Rate of TEC Index (ROTI), the global ionospheric effects at different sectors in August 2014 were analyzed. Possible mechanisms responsible for ionospheric disturbances were also discussed. It is found that the ionosphere in the geomagnetic storm has an obvious ionospheric storm effect. During the magnetic storm, the southern hemisphere was dominated by positive storms while the northern hemisphere was characterized by a positive storm followed by a long duration strong negative storm, and the negative phase storm duration was significantly longer than the positive storm; the amplitude variation at low latitude was smaller than that at mid-latitude and the high latitude mainly performed negative storm effects; the peaks of Equatorial Ionospheric Anomaly (EIA) disappeared after the magnetic storm occurs due to the movement of EIA towards the equator. The analysis of ionospheric irregularities during geomagnetic storms shows that ionospheric scintillation occurs before the magnetic storm at European-African sector and disappears after the geomagnetic storm. The results show that there are obvious temporal and spatial differences in the ionospheric changes during the storm.
Ionosphere Peak Height Prediction Performance of International Reference Ionosphere Model
TANG Siyu, ZHOU Chunhong, HUANG Zhi, YUAN Hong
2021, 41(6): 898-904. doi: 10.11728/cjss2021.06.898
Abstract:
Because ionospheric peak height HmF2 is one of the important parameters describing the ionospheric morphology, new models based on digisonde station data and radio occultation data was included in the latest version IRI-2016, and the changes and improvements of the HmF2 prediction would be expected. In this work, Global three dimension variations of HmF2 from the Constellation Observing System for Meteorology, Ionosphere and Climate (COSMIC) are described using the 4- year dataset during the period of 2007-2010. Meanwhile, HmF2 predictions by the model IRI-2016 and IRI-2012 are analyzed in order to evaluate the performance of IRI-2016. Comparisons between COSMIC and IRI show that IRI-derived HmF2 values are high in the middle and high latitude regions and whereas they are relatively low in the equator and low latitude regions. During the period of midnight, HmF2 by IRI-2016 is characteristic of latitudinal gradient, and notably is higher than that by IRI-2012 in the part of the mid-high and the low-latitude regions. However, HmF2 from model IRI-2016 in most equator regions and low-latitudes is lower than IRI-2012 in the daytime. Analysis of ionospheric HmF2 provides valuable references for the future improvement of ionosphere model.
Optimization of Atmospheric Radiative Transfer Model LBLRTM Based on Measured CO2 Data
ZHANG Bingyan, YAN Zhaoai, GUO Wenjie, HU Xiong
2021, 41(6): 905-910. doi: 10.11728/cjss2021.06.905
Abstract:
According to the CO2 observation data of TIMED/SABER from 2002 to 2018, analyze the variation characteristics of CO2 concentration, gives the monthly average fitting formula of CO2 concentration varying with time, height and latitude according to the various characteristics, and uses the nonlinear least square fitting method to fit the CO2 concentration data of different heights and latitudes respectively to generate corresponding fitting parameters, Then, summarize all fitting parameters and generate the fitting parameter file. Combined with the fitting formula, the empirical calculation module of global CO2 concentration is constructed, and the module is applied to the atmospheric radiative transfer model LBLRTM to optimize LBLRTM. Comparing the simulation results of the optimized LBLRTM model with the TIMED/SABER observation data, the root mean square error between the simulation results of the not optimized LBLRTM model and the observation value is 15.4%, while the root mean square error between the simulation results of the optimized LBLRTM model and the observation value is reduced to 8.91%. The results show that this optimization method can further improve the radiation simulation accuracy of the LBLRTM model in the infrared band.
Statistical Characteristics of Stratospheric Mountain Waves over Southern Andes Based on AIRS Observations
HUANG Chouyue, ZHANG Shaodong, HUANG Chunming, HUANG Kaiming, GONG Yun
2021, 41(6): 911-919. doi: 10.11728/cjss2021.06.911
Abstract:
Using the radiation measurement data of the Atmospheric Infrared Sounder (AIRS) on the Aqua satellite from 2013 to 2018, case studies and statistical analysis of mountain waves at 20 km, 27 km, 35 km, and 41 km in the Andes are carried out. Observation results show that the mountain waves over the Andes from 2013 to 2018 mainly occur from May to October, and the monthly average horizontal wavelength, vertical wavelength, and momentum flux have no obvious interannual changes.The horizontal wavelength decreases slightly in May and October compared to June to September. The vertical wavelength gradually increases from May to June and reaches its peak in July. Then it decreases from August to October. The momentum flux shows a similar seasonal changes pattern as the vertical wavelength. In the range of 20 km to 41 km, the horizontal wavelength slowly increases from 43.5~53.9km to 89.3~176.8km. Affected by the background wind field, the vertical wavelength increases with height, from 7.4~14.7 km at 20 km to 7.4~29.7km at 41 km. The momentum flux ranges from 376.0~801.3mPa at 20 km and drops to 10.4~239.3mPa at 41 km. In general, the horizontal wavelength of the mountain waves increases slowly with height as it propagates upward, and the vertical wavelength increases with height under the influence of the background wind field in the case of upstream propagation. The momentum flux decreases significantly with height, indicating that the Andes mountain waves propagate upward while accompanied by strong dissipation. The dissipated energy will be stored in the background atmosphere, which will have an important impact on the upper stratosphere and even the mesosphere.
A Novel Magnetic Configuration for Space Radiation Active Shielding
CAI Minghui, XU Liangliang, YANG Tao, HAN Jianwei
2021, 41(6): 920-927. doi: 10.11728/cjss2021.06.920
Abstract:
Space radiation has been identified as the main health hazard to crews involved in manned Mars missions. Active shielding is more effective than passive shielding to the very energetic particles from cosmic rays. Particle motion in a magnetic field is studied based on the single-particle theory and Monte Carlo method. By comparing the shielding efficiency of different magnetic field configurations, a novel active magnetic shielding configuration with lower mass cost and power consumption is proposed for manned Mars missions. The new magnetic configuration can shield 92.8% of protons and 84.4% of alpha particles with E < 4 GeV·n-1, when considering the passive shielding contribution of 10.0 g·cm-2 Al Shielding, the required magnetic stiffness can be reduced from 27 Tm to 16 Tm. The detailed analysis of mass cost and power consumption shows that active shielding will be a promising means to protect crews from space radiation exposure in manned Mars missions.
Linear Surrogate Uncertainty Analysis Method for Distributed Satellite System
GAO Chen, YANG Zhen, NIU Wenlong
2021, 41(6): 928-935. doi: 10.11728/cjss2021.06.928
Abstract:
Distributed satellite system is widely used in space science mission. The complexity of system brings challenges in conceptual design phase. To evaluate the impact of uncertainty parameters on the system detection result is significant. Traditional analysis method such as differential method and Monte Carlo method have disadvantages in solving this problem. Differential method needs system can be expressed explicit while Monte Carlo method is time consuming which cannot fit fast iteration demand. In this paper, a surrogated linear model method is proposed. The problem is solved by reducing the dimensionality of the results of the uncertainty parameters into a linear combination of coefficients. The coefficients are obtained by weighted regression method using data generated by Monte Carlo process. Simulated with a typical distributed satellite mission, the proposed method is valid compared to Monte Carlo simulation method.
Reliability Analysis of Space Tether Collision Based on Fuzzy Stress-strength Model
HAN Zhichen, LI Aijun, WANG Changqing, LI Xianpeng
2021, 41(6): 936-944. doi: 10.11728/cjss2021.06.936
Abstract:
Space tether system has attracted more and more attention from various countries because of its special structure. The collision reliability research of space tether is an important part of the mission design of space tether. In this paper, the stress-strength model is used to analyze the reliability of space tethers on-orbit collision based on the space debris flux and Poisson distribution on low-Earth orbit. Several tether structure factors those have an impact on the collision reliability of space tethers are proposed, such as tether's diameter and length for single-strand space tether, and tether's distance between ropes, collision angle for double-strand tether. Firstly, the tether residual cross-section at the collision point is modeled after the space debris impacts and cuts the tether. Then, according to the fuzzy stress-strength model analysis, the tether destructive collision probability is calculated. Finally, according to the Poisson method, the reliability of the space tether is calculated over time, and the simulation analysis is carried out to compare the effective on-orbit time of space tether with different structures.
Design and Implementation of a Chirp Transform Spectrometer for Deep Space Exploration
WANG Yehuan, DONG Xiaolong, ZHU Di, LIU Mengwei, GONG Junjie
2021, 41(6): 945-953. doi: 10.11728/cjss2021.06.945
Abstract:
According to the demand of radiation resistance, high resolution and low power consumption for spectrometer in deep space exploration, the design and implementation of a high resolution spectrometer based on Chirp transform are studied. The key problem in the implementation of the high-resolution spectrometer based on a SAW filter with 1 GHz center frequency, 400 MHz bandwidth and 10 µs dispersion time are studied. The system is designed and built by the combination of analog method and digital method, and the basic functions of the spectrometer are realized. The bandwidth of the system is 400 MHz, the frequency resolution is 152 kHz and the power consumption is about 3.6 W. The influence of the nonideal behavior of the compressor device on frequency resolution is analyzed. Amplitude compensation curve and phase compensation curve are carried out based on frequency response characteristics of compressor device. The amplitude mismatch and phase mismatch of the signal generated by the expander are compensated. The simulation results show that the frequency resolution of the system after compensation increases to 108 kHz.
Configuration Keeping Control of Stereo Imaging with Dual-satellite Following Flying Formation
CHEN Gaojie, CHANG Lin, LI Jing, YANG Xiubin, YANG Chunlei, LI Yanbo
2021, 41(6): 954-961. doi: 10.11728/cjss2021.06.954
Abstract:
Aiming at the configuration maintenance problem of stereo imaging with micro double satellites under strong interference and output saturation, a composite anti-disturbance control strategy combining LQR optimal control and PD control is proposed. According to the relative orbit motion mechanism of dual-satellite following flying formation for stereo imaging, the dynamic modeling of relative motion is built with a dynamic analysis of double satellite following formation carried out, and the system motion state equation is obtained. Secondly, an adaptive disturbance observer is designed, which can estimate the system state and disturbance information at the same time. The existence condition of the observer is given by using Lyapunov stability theory and linear matrix inequality. In addition, the pole placement method is used to improve the observer dynamic performance, and the exponential attenuation factor is introduced to improve the convergence speed of the controller. Considering the output saturation characteristics of the actuator, a composite control strategy of weighted PD + LQR feedback and disturbance feedforward compensation is designed to suppress the influence of unknown disturbance and ensure the transient response and steady state performance of the system. Simulation results have verified the effectiveness of the proposed algorithm, which can provide configuration keeping control strategy for dual star stereo imaging.
Combinational Model for Real-time Computing Rocket Falling Point Based on Space Delamination
WANG Qiang, ZHANG Wei, REN Zhenzhen, GONG Jianze, JING Yuhang
2021, 41(6): 962-967. doi: 10.11728/cjss2021.06.962
Abstract:
A mathematical model based on space delamination combining geometric and dynamical models for falling point calculation is proposed to ensure the safety of the space launch mission's flight zone and to achieve high accuracy and fast calculation of rocket falling point in real time. The model stratifies the trajectory into a free segment and a re-entry segment. In the free segment, an elliptical geometric model is used considering Earth oblateness correction and in the re-entry segment, a dynamic model is used adding air resistance. Corrections for factors such as geophysical disturbance and garth rotation are also incorporated. Through validation of actual falling point data, it is found that the model can maximize the integration of accuracy and efficiency, and could be used to calculate rocket falling point for space launch missions with high accuracy and real-time performance.
A Modulation Recognition Algorithm Based on Wavelet Transform Entropy and High-order Cumulant for Satellite Signal Modulation
YAN Wenkang, YAN Yi, FAN Yanan, YAO Xiujuan, GAO Xiang, SUN Wen
2021, 41(6): 968-975. doi: 10.11728/cjss2021.06.968
Abstract:
Modulation recognition is the key link of signal detection and demodulation. Aiming at MAPSK, MQAM, MFSK and MPSK modes used in satellite modulation, a joint modulation recognition algorithm is proposed, which calculates the entropy of wavelet transform and combines highorder cumulant. According to the characteristics of the wavelet transform pair frequency sensitive information, different modulation method to distinguish the results of the calculation of higher-order cumulant and the modulation signal of different complexity entropy results, the calculation results of the above four 4 kinds of modulation signal are analyzed, and entropy is proposed based on wavelet transform and higher-order cumulants joint satellite signal modulation recognition algorithm. Based on the calculated wavelet coefficients of modulation signals, the classification of modulation signals is realized by calculating the entropy value, and the signal classification within the class is realized by using the high-order cumulant. Through simulation analysis, the recognition effect of 0.9 above 8 dB can be achieved. In addition, the method has reference significance for high order (64 order modulation) signals.
Optimization and Implementation of Two-dimensional Fast Acquiring DSSS Burst Datagram Communication for LEO Satellite Constellation
ZHANG Jun, ZHENG Lirong, LIN Baojun, SHEN Yuan, GONG Wenbing, SHUAI Tao
2021, 41(6): 976-982. doi: 10.11728/cjss2021.06.976
Abstract:
Large Doppler shifts and extremely weak signal in the condition of short datagram burst LEO satellite communication are always the technical challenges for fast acquisition receiver. Firstly, the influence of Doppler characteristics on the fast acquisition of LEO satellite spread spectrum signal is analyzed, and the advantages and disadvantages of different acquisition strategies are discussed. Based on double-sampling zero-padding technology of partial matched filter and Fast Fourier Transform (FFT), Optimization and implementation of two-dimensional parallel fast acquiring Direct Sequence Spread Spectrum (DSSS) signal processing method is presented. Different from the single sampling scheme in which the input signal is sampled according to the chip rate, the double sampling scheme has two sampling signals reaching the correlator in each chip time. That is, the correlator works on a clock with twice than the chip rate. The correlator selects a sampling signal at each interval and inputs it to the PN code part correlator for correlation operation, and then outputs the maximum value after passing through the FFT module and energy accumulation module for judgment. Therefore, the maximum sampling signal-to-noise ratio loss is about 0.9 dB, which greatly improves the performance compared with the single sampling method. Finally, combined with the decision oriented carrier forward estimation and PN code feedback synchronization technology, the actual test shows that the average acquisition time of the receiver is 10 ms, which is much lower than the tradition parallel acquisition time of PN sequence code in seconds order. It can better meet the service requirements of short datagram burst communication system with random access for LEO satellite constellation.
Numerical Investigation on Collector of Asteroid's Particles with Various Sizes
DU Yonggang, YAN Chunjie, WANG Chunyong, WANG Jin
2021, 41(6): 983-987. doi: 10.11728/cjss2021.06.983
Abstract:
Collecting and returning samples from the asteroid's surface is one of the significant targets of the current asteroid detection. The key issue of the nitrogen purge sampling technology is the flow field design of the sample container. In order to ensure that the sample container has an efficient collection performance, the flow field design of the sample container must solve two key problems, one problem is the flow field how to form a trap effect on the sample particles, and another problem is the flow field how to effectively classify the sample particles. Therefore, the researchers proposed a scheme of the sample container, which uses the combined effect of centrifugal force and guide blade to collect and classify the sample particles. For verifying the performance and feasibility of this scheme in a microgravity environment, the researchers carried out the numerical simulation with Fluent software. Under the boundary conditions of the container, the simulation program calculated the velocity of nitrogen at the inlet, and the optimal value is 0.6 m·s-1, the efficiency of the container at this velocity reaches 99.7%, the flow field of the sample container also has a good particle classification function, and its capture rate of particles that less than 2mm exceeds 90%. Its pressure loss is also within a reasonable range, the research results of this paper can provide a reference for the upcoming asteroid's sampling task in China.