Spacecraft formation flying control around L2 sun-earth libration point using on–off SDRE approach

In this paper, on–off SDRE control approach is presented for spacecraft formation flying control around sun-earth L2 libration point. Orbits around libration points are significant targets for many space missions mainly because of efficient fuel consumption. Furthermore, less propellant usage can be achieved by considering optimal control approaches in spacecraft formation flying control design. Among various nonlinear and optimal control methods, SDRE has shown to be a popular controller in various missions due to the privileges including efficiency, accuracy and robustness. The spacecraft are assumed to have on–off thrusters as actuators. It requires them to be fed with a sequence of on–off pulses which is regarded as a challenge for spacecraft designers. Hence, the main contribution of this paper is designing an on–off SDRE approach for the formation flight around sun-earth L2 point with uncertainty with energy and accuracy considerations. Including on–off input as a constraint is not feasible for SDRE implementation because it makes the system non-affine. An alternative is utilizing an integral action technique and an auxiliary control to make the system affine which leads to on–off SDRE approach. It has also been shown that the proposed method is robust against parametric uncertainties of the states. Present study aims to design an energy-beneficial, simple and attractive controller for a complex nonlinear system with on–off inputs and uncertainty in CRTBP. Simulation results show that the on–off SDRE control could provide the formation flight around L2 point with high accuracy using less energy consumption.     

Ionospheric total electron content response to annular solar eclipse on June 21, 2020

The effects of physical events on the ionosphere structure is an important field of study, especially for navigation and radio communication. The paper presents the spatio-temporal ionospheric TEC response to the recent annular solar eclipse on June 21, 2020, which spans across two continents, Africa and Asia, and 14 countries. This eclipse took place on the same day as the June Solstice. The Global Navigation Satellite System (GNSS) based TEC data of the Global Ionosphere Maps (GIMs), 9 International GNSS Service (IGS) stations and FORMOSAT-7/COSMIC-2 (F7/C2) were utilized to analyze TEC response during the eclipse. The phases of the TEC time series were determined by taking the difference of the observed TEC values on eclipse day from the previous 5-day median TEC values. The results showed clear depletions in the TEC time series on June 21. These decreases were between 1 and 9 TECU (15–60%) depending on the location of IGS stations. The depletions are relatively higher at the stations close to the path of annular eclipse than those farther away. Furthermore, a reduction of about ?10 TECU in the form of an equatorial plasma bubble (EPB) was observed in GIMs at ~20° away from the equator towards northpole, between 08:00–11:00 UT where its maximum phase is located in southeast Japan. Additionally, an overall depletion of ~10% was observed in F7/C2 derived TEC at an altitude of 240 km (hmF2) in all regions affected by the solar eclipse, whereas, significant TEC fluctuations between the altitudes of 100 km ? 140 km were analyzed using the Savitzky-Golay smoothing filter. To prove TEC depletions are not caused by space weather, the variation of the sunspot number (SSN), solar wind (V_SW), disturbance storm-time (Dst), and Kp indices were investigated from 16th to 22nd June. The quiet space weather before and during the solar eclipse proved that the observed depletions in the TEC time series and profiles were caused by the annular solar eclipse.     

Equatorial anomaly according to the Interkosmos-19 data and IRI model: A comparison

The comparison of the IRI model with the foF2 distribution in the equatorial anomaly region obtained by topside sounding onboard the Interkosmos-19 satellite has been carried out. The global distribution of foF2 in terms of LT-maps was constructed by averaging Intercosmos-19 data for summer, winter, and equinox. These maps, in fact, represent an empirical model of the equatorial anomaly for high solar activity F10.7 ~ 200. The comparison is carried out for the latitudinal foF2 profiles in the characteristic longitudinal sectors of 30, 90, 210, 270, and 330°, as well as for the longitudinal variations in foF2 over the equator. The largest difference between the models (up to 60%) for any season was found in the Pacific longitudinal sector of 210°, where there are a few ground-based sounding stations. Considerable discrepancies, however, are sometimes observed in the longitudinal sectors, where there are many ground-based stations, for example, in the European or Indian sector. The discrepancies reach their maximum at 00 LT, since a decay of the equatorial anomaly begins before midnight in the IRI model and after midnight according to the Interkosmos-19 data. The discrepancies are also large in the morning at 06 LT, since in the IRI model, the foF2 growth begins long before sunrise. In the longitudinal variations in foF2 over the equator at noon, according to the satellite data, four harmonics are distinguished in the June solstice and at the equinox, and three harmonics in the December solstice, while in the IRI model only two and one harmonics respectively are revealed. In diurnal variations in foF2 and, accordingly, in the equatorial anomaly intensity, the IRI model does not adequately reproduce even the main, evening extremum.     

微等离子体转化二氧化碳发射光谱诊断

火星大气层的主要成分为二氧化碳，如果能够利用低温等离子体方法高效分解二氧化碳，使其转化为氧气和一氧化碳加以利用，可以大幅降低航天员生命保障相关载荷长途运输的成本，进一步提高生命保障能力。低温等离子体放电过程中会产生大量活性组分，可以在数百度温度下实现二氧化碳的高效解离，是具有很大潜力的二氧化碳解离与转化方式。设计了一种尺度在亚毫米级、功率输入为数瓦的直流微槽等离子体反应器，可以在较低气体温度下实现二氧化碳分解。测量了反应器电流、功率等放电参数，采用发射光谱确定了体系中激发态组分，分析了激发态粒子谱线强度随输入电压、稀释气体比例等反应器工作参数变化，利用氮气分子振转谱带测量了等离子体放电区振动温度和气体温度。研究表明，添加氩、氦、氮气均可以增强二氧化碳的分解，添加氦气可以促进二氧化碳的电离过程。稀释气体激发态因具有高能量，可以通过潘宁解离通道增强二氧化碳分解。氦组分激发态的能量高于二氧化碳电离能，可以促进二氧化碳的电离反应。微等离子体内存在强烈的振动 平动非平衡现象：振动温度约为4400~4800K，而气体温度仅为450 ~600K，表明可以通过合理的放电和结构设计，定向将能量注入到振动态，从而进一步促进二氧化碳的振动解离。     

Evaluation of the IRI-2016 and NeQuick electron content specification by COSMIC GPS radio occultation,ground-based GPS and Jason-2 joint altimeter/GPS observations

We examined performance of two empirical profile-based ionospheric models, namely IRI-2016 and NeQuick-2, in electron content (EC) and total electron content (TEC) representation for different seasons and levels of solar activity. We derived and analyzed EC estimates in several representative altitudinal intervals for the ionosphere and the plasmasphere from the COSMIC GPS radio occultation, ground-based GPS and Jason-2 joint altimeter/GPS observations. It allows us to estimate a quantitative impact of the ionospheric electron density profiles formulation in several altitudinal intervals and to examine the source of the model-data discrepancies of the EC specification from the bottom-side ionosphere towards the GPS orbit altitudes. The most pronounced model-data differences were found at the low latitude region as related to the equatorial ionization anomaly appearance. Both the IRI-2016 and NeQuick-2 models tend to overestimate the daytime ionospheric EC and TEC at low latitudes during all seasons of low solar activity. On the contrary, during high solar activity the model results underestimated the EC/TEC observations at low latitudes. We found that both models underestimated the EC for the topside ionosphere and plasmasphere regions for all levels of solar activity. For low solar activity, the underestimated EC from the topside ionosphere and plasmasphere can compensate the overestimation of the ionospheric EC and, consequently, can slightly decrease the resulted model overestimation of the ground-based TEC. For high solar activity, the underestimated EC from the topside ionosphere and plasmasphere leads to a strengthening of the model underestimation of the ground-based TEC values. We demonstrated that the major source of the model-data discrepancies in the EC/TEC domain comes from the topside ionosphere/plasmasphere system.     

Plasma transport process in the equatorial/low-latitude ionosphere

The behaviour of the equatorial/low-latitude ionosphere and the transport processes during magnetic disturbed and quiet periods of a high solar activity year, 2014, in the American sector are investigated. Parameters used include vertical drift (Vz), transport term (W), NmF2, hmF2 and scale-height (H). The F2 plasma variations followed the diurnal local solar pattern, being higher at daytime. The sunset maximum and sunrise minimum peaks of hmF2 were directly opposite to the scale height (H) pattern. The plasma distribution was basically controlled by combined actions of the electrodynamic convection/thermospheric composition, which is geomagnetic activity dependent. The annual, semi-annual and winter-anomalies of the F2 parameters were higher at the dip equator in comparison with the low-latitude. The Vz pre-reversal peak magnitude coincided with hmF2 peak and the effects are more pronounced during geomagnetic disturbed conditions. The transport term pattern was similar to that of the scale height and it is suggested as a proxy parameter for quantifying low-latitude plasma irregularities and distribution of thermospheric composition.     

The adjusted optical properties for Galileo/BeiDou-2/QZS-1 satellites and initial results on BeiDou-3e and QZS-2 satellites

Solar Radiation Pressure (SRP) is the dominant non-gravitational perturbation for GNSS (Global Navigation Satellite System) satellites. In the absence of precise surface models, the Empirical CODE Orbit Models (ECOM, ECOM2) are widely used in GNSS satellite orbit determination. Based on previous studies, the use of an a priori box-wing model enhances the ECOM model, especially if the spacecraft is a stretched body satellite. However, so far not all the GNSS system providers have published their metadata. To ensure a precise use of the a priori box-wing model, we estimate the optical parameters of all the Galileo, BeiDou-2, and QZS-1 (Quasi Zenith Satellite System) satellites based on the physical processes from SRP to acceleration. Validation using orbit prediction proves that the adjusted parameters of Galileo and QZS-1 satellites exhibit almost the same performance as the corresponding published and “best guess” values. Whereas, the estimated parameters of BeiDou-2 satellites demonstrate an improvement of more than 60% over the initial “guess” values. The resulting optical parameters of all the satellites are introduced into an a priori box-wing model, which is jointly used with ECOM and ECOM2 model in the orbit determination. Results show that the pure ECOM2 model exhibits better performance than the pure ECOM model for Galileo, BeiDou-2 GEO and QZS-1 orbits. Combined with the a priori box-wing model the ECOM model (ECOM+BW) results in the best Galileo, BeiDou-2 GEO and QZS-1 orbits. The standard deviation (STD) of satellite laser ranging residuals reduce by about 20% and 5% with respect to the pure ECOM2 model for Galileo and BeiDou-2 GEO orbits, while the reductions are about 40% and 60% for QZS-1 orbits in yaw-steering and orbit-normal mode respectively. BeiDou-2 IGSO and MEO satellite orbits do not benefit much from the a priori box-wing model. In summary, we suggest setting up a unified SRP model of ECOM+BW for Galileo, QZS-1, and BeiDou-2 orbits based on the adjusted metadata. In addition, we estimate the optical parameters of BeiDou-3e and QZS-2 satellites using a limited number of tracking stations. Results regarding the unified SRP model indicate the same advantages, the STD of satellite laser ranging residuals reduces by about 30% and 20% for QZS-2 and BeiDou-3e orbits respectively over orbit products without a priori model. The estimation procedure is effective and easy to apply to the new emerging satellites in the future.     

螺旋桨滑流对发动机进气道气动性能的影响

针对螺旋桨/短舱/进气道一体化构型,采用多重参考系(MRF)法和滑移网格非定常法进行了数值模拟研究,实现了螺旋桨/短舱/进气道一体化内外流耦合流场仿真,研究了不同工况下螺旋桨滑流对进气道气动性能的影响,并与无桨构型的进气道性能进行了对比,最后通过数值分析解释了螺旋桨桨盘后总压上升的原因。结果表明:多重参考系法和滑移网格非定常法对进气道出口的时均参数的计算结果吻合较好,多重参考系法可用于螺旋桨及进气道的气动性能参数计算;螺旋桨做功使进气道总压恢复系数增高,螺旋桨转速越大,总压恢复系数的增量就越大,且螺旋桨滑流使进气道出口畸变指数明显增大;通过数值分析推导的桨盘前后总压关系公式与仿真结果的误差均在3%以内。      

逆向进气复合角发散冷却孔流量系数试验

针对回流燃烧室外环壁面冷气与燃气流动方向相反的情况,试验研究了发散孔进气方向与环腔通道气流流动方向相反的条件下,复合角发散冷却孔流量系数随发散孔倾角(25°~55°)、偏角(90°~180°)及压降系数等参数的变化规律。研究结果表明:复合角发散孔的流量系数随倾角的增大而增大,随偏角的增大而减小；随着压降系数的增大,复合角发散孔的流量系数先呈近似线性增长,当压降系数大于0.95后,呈非线性增长且增长速度减缓。      

过氧化氢/紫外线/臭氧氧化技术处理肼类污水的应用

采用放大试验规模（处理量1 m³/h）的H₂O₂（过氧化氢）/UV（紫外线）/O₃（臭氧）氧化技术处理肼类推进剂污水，在30.0±1.6 ℃，pH=9.0±0.2的条件下，研究肼类的过氧化氢增强光解臭氧化降解。对比不同氧化技术的协同效应以及对污水降解的影响，重点考察过氧化氢、紫外线、臭氧和初始浓度等工艺参数对降解效果的影响，对氧化技术的应用进行了优化。研究结果表明：该技术可使COD去除率提高27.66%，肼类的降解速率随着过氧化氢投加量、紫外线辐射强度、O₃投加速率和水质地提高而升高，随着初始质量浓度的提高而下降。在最佳工艺下，处理5 000 mg/L质量浓度的废水，处理60 min时，COD去除率分别为98.62%（偏二甲肼）、99.17%（甲基肼）、99.94%（肼）和93.25%（单推-3）。