An empirical solar radiation pressure model for satellites moving in the orbit-normal mode

We present a family of empirical solar radiation pressure (SRP) models suited for satellites orbiting the Earth in the orbit normal (ON) mode. The proposed ECOM-TB model describes the SRP accelerations in the so-called terminator coordinate system. The choice of the coordinate system and the SRP parametrization is based on theoretical assumptions and on simulation results with a QZS-1-like box-wing model, where the SRP accelerations acting on the solar panels and on the box are assessed separately. The new SRP model takes into account that in ON-mode the incident angle of the solar radiation on the solar panels is not constant like in the yaw-steering (YS) attitude mode. It depends on the elevation angle of the Sun above the satellite’s orbital plane. The resulting SRP vector acts, therefore, not only in the Sun-satellite direction, but has also a component normal to it. Both components are changing as a function of the incident angle. ECOM-TB has been used for precise orbit determination (POD) for QZS-1 and BeiDou2 (BDS2) satellites in medium (MEO) and inclined geosynchronous Earth orbits (IGSO) based on IGS MGEX data from 2014 and 2015. The resulting orbits have been validated with SLR, long-arc orbit fits, orbit misclosures, and by the satellite clock corrections based on the orbits. The validation results confirm that—compared to ECOM2—ECOM-TB significantly (factor 3–4) improves the POD of QZS-1 in ON-mode for orbits with different arc lengths (one, three, and five days). Moderate orbit improvements are achieved for BDS2 MEO satellites—especially if ECOM-TB is supported by pseudo-stochastic pulses (the model is then called ECOM-TBP). For BDS2 IGSOs, ECOM-TB with its 9 SRP parameters appears to be over-parameterized. For use with BDS2 IGSO spacecraft we therefore developed a minimized model version called ECOM-TBMP, which is based on the same axis decomposition as ECOM-TB, but has only 2 SRP parameters and is supported by pseudo-stochastic parameters, as well. This model shows a similar performance as ECOM-TB with short arcs, but an improved performance with (3-day) long-arcs. The new SRP models have been activated in CODE’s IGS MGEX solution in Summer 2018. Like the other ECOM models the ECOM-TB derivatives might be used together with an a priori model.     

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.     

Rotating machinery fault diagnosis based on convolutional neural network and infrared thermal imaging

Rotating machinery is widely applied in industrial applications. Fault diagnosis of rotating machinery is vital in manufacturing system, which can prevent catastrophic failure and reduce financial losses. Recently, Deep Learning (DL)-based fault diagnosis method becomes a hot topic. Convolutional Neural Network (CNN) is an effective DL method to extract the features of raw data automatically. This paper develops a fault diagnosis method using CNN for InfRared Thermal (IRT) image. First, IRT technique is utilized to capture the IRT images of rotating machinery. Second, the CNN is applied to extract fault features from the IRT images. In the end, the obtained features are fed into the Softmax Regression (SR) classifier for fault pattern identification. The effectiveness of the proposed method is validated using two different experimental data. Results show that the proposed method has a superior performance in identification various faults on rotor and bearings comparing with other deep learning models and traditional vibration-based method.     

27-day variation in solar-terrestrial parameters: Global characteristics and an origin based approach of the signals

The Earth and the near interplanetary medium are affected by the Sun in different ways. Those processes generated in the Sun that induce perturbations into the Magnetosphere-Ionosphere system are called geoeffective processes and show a wide range of temporal variations, like the 11-year solar cycle (long term variations), the variation of ～27?days (recurrent variations), solar storms enduring for some days, particle acceleration events lasting for some hours, etc.In this article, the periodicity of ～27?days associated with the solar synodic rotation period is investigated. The work is mainly focused on studying the resulting 27-day periodic signal in the magnetic activity, by the analysis of the horizontal component of the magnetic field registered on a set of 103 magnetic observatories distributed around the world. For this a new method to isolate the periodicity of interest has been developed consisting of two main steps: the first one consists of removing the linear trend corresponding to every calendar year from the data series, and the second one of removing from the resulting series a smoothed version of it obtained by applying a 30-day moving average. The result at the end of this process is a data series in which all the signal with periods larger than 30?days are canceled.The most important characteristics observed in the resulting signals are two main amplitude modulations: the first and most prominent related to the 11-year solar cycle and the second one with a semiannual pattern. In addition, the amplitude of the signal shows a dependence on the geomagnetic latitude of the observatory with a significant discontinuity at approx. ±60°.The processing scheme was also applied to other parameters that are widely used to characterize the energy transfer from the Sun to the Earth: F10.7 and Mg II indices and the ionospheric vertical total electron content (vTEC) were considered for radiative interactions; and the solar wind velocity for the non-radiative interactions between the solar wind and the magnetosphere. The 27-day signal obtained in the magnetic activity was compared with the signals found in the other parameters resulting in a series of cross-correlations curves with maximum correlation between 3 and 5?days of delays for the radiative and between 0 and 1?days of delay for the non-radiative parameters. This result supports the idea that the physical process responsible for the 27-day signal in the magnetic activity is related to the solar wind and not to the solar electromagnetic radiation.     

高能(450 kV)X射线参考辐射场附加过滤研究

利用MCNP程序对450kV X光机建立模型，计算各参考辐射质下的原级谱、平均光子能量。对ISO 4037-1中已有范围的管电压与附加过滤片厚度关系曲线进行二次拟合，确定模拟厚度范围。利用MCNP计算各附加过滤组合下的平均光子能量值，并与外推法得出的平均光子能量值进行对比，将相对偏差最小时的附加过滤作为最终确定值。从部分能量点的模拟与测量结果的对比可看出，模拟得出的平均光子能量在2%以内与PTB及中国计量院测量结果保持一致，附加过滤最大相对偏差为18%。     

A time-frequency analysis of the thermal radiation background anomalies caused by large earthquakes: A case study of the Wenchuan 8.0 earthquake

Strong earthquakes have an impact on the regional thermal radiation background, which has been both observed and confirmed. This effect produces anomalies in the thermal radiation background (TRBA) and increases the difficulty of extracting a thermal radiation anomaly (TRA) that is associated with an earthquake occurring during the same time period. The extraction and identification of such anomalies has been ignored by previous studies. In this study, we investigate the time-frequency analysis (TFA) method, together with the wavelet filtering of the Daubechies method and the relative power spectrum analysis of the Fourier Welch method to extract and analyse the TRBA caused by the 2008 Wenchuan Ms 8.0 earthquake and the TRA of the 2013 Minxian-Zhangxian Ms 6.6 earthquake using data concerning the brightness temperature of a black body (TBB) from the Fengyun-2 series of geostationary meteorological satellites developed by China. The result showed that this method can effectively extract and analyse the TRBA caused by the Wenchuan earthquake and the TRA of the Minxian-Zhangxian earthquake form a complex background environment. Furthermore, we discussed the impact of the earthquake on the TRBA and segmented the process. The impact is mainly reflected by three aspects; the characteristic period of the TRBA changes, the TRBA occurs at the same time every year, which is identical to the time at which the earthquake anomaly occurred, and the impact process is in stages. We also summarized the correlation between the characteristic parameters of a TRA and the regional thermal radiation background, geography, and climatic factors.     

S弯二元喷管红外辐射特性实验

基于轴对称排气系统设计了一套S弯二元喷管的实验模型,在实验工况下测量了喷管实验件的部分流场参数和远场红外辐射特性,对S弯二元喷管的红外辐射特性进行了分析和研究,并与基准轴对称喷管的红外辐射实验结果进行了对比。结果表明:上方探测平面是S弯二元喷管的主要辐射方向,其最大辐射角度为15°,在该方向角上固体壁面辐射的贡献最大。下方和侧向探测平面各方向角的红外辐射强度明显小于上方探测平面,主要是由于S弯二元喷管遮挡了大部分的内部高温壁面。与基准轴对称喷管对比,S弯二元喷管有着明显的红外抑制作用,在尾部方向S弯二元喷管的积分辐射强度相比轴对称喷管降低了81%,最大辐射强度相比轴对称喷管降低了47%,燃气辐射降低50%。      

高空飞行环境中液体运载火箭底部热环境研究

采用数值模拟和飞行测试验证相结合的方法对液体运载火箭高空对流/辐射耦合换热问题开展系统深入研究。基于燃气多组分输运Navier-Stokes方程、热辐射方程、Realizable k-ε两方程湍流模型，建立了高空含自由流的运载火箭燃气喷流流动模型。辐射模型采用离散坐标法(DOM)，空间离散采用二阶迎风TVD格式，对多个典型飞行高度火箭底部热流进行大型并行计算，将数值结果与试验数据进行广泛对比，验证了计算模型的精度和有效性。数值研究表明，火箭底部辐射热流在刚起飞阶段达到最大值，随着飞行高度上升，辐射热流逐渐降低，火箭底部对流热流表现为先升高后降低的趋势，并在20 km高空达到峰值。本文的预测分析方法对液体运载火箭底部热防护设计具有重要的理论意义和工程应用价值。     

基于SRAM型FPGA单粒子效应综合防护技术

为解决基于SRAM(Static Random Access Memory,静态随机存储器)型FPGA(Field Programmable Gate Array,现场可编程门阵列)的扩频应答机在空间环境中时常发生功能故障的难题,提出了一种低轨卫星基于SRAM型FPGA的应答机的单粒子综合防护方案。该方案在代码级、器件级、单机级、整星级综合采用器件优选、结构屏蔽设计、巡检自愈、自主复位等多种方法,有效杜绝了单粒子事件带来的种种危害,已成功应用于某在轨卫星。在轨测试结果表明:遥测、遥控、跟踪测量等功能、性能指标满足任务要求,且未发生需要地面干预的问题或现象。该技术还可推广应用于其他采用低等级器件及其他类型FPGA、DSP(Digital Signal Processing,数字信号处理)等大规模集成电路的星载单机防护。     

Improved rainfall estimation over the Indian region using satellite infrared technique

The GOES Precipitation Index (GPI) technique (Arkin, 1979) for rainfall estimation has been in operation for the last three decades. However, its applications are limited to the larger temporal and spatial scales. The present study focuses on the augmentation on GPI technique by incorporating a moisture factor for the environmental correction developed by Vicente et al. (1998). It consists of two steps; in the first step the GPI technique is applied to the Kalpana-IR data for rainfall estimation over the Indian land and oceanic region and in the second step an environmental moisture correction factor is applied to the GPI-based rainfall to estimate the final rainfall. Detailed validation with rain gauges and comparison with Tropical Rainfall Measuring Mission (TRMM) merged data product (3B42) are performed and it is found that the present technique is able to estimate the rainfall with better accuracy than the GPI technique over higher temporal and spatial domains for many operational applications in and around the Indian regions using Indian geostationary satellite data. Further comparison with the Doppler Weather Radar shows that the present technique is able to retrieve the rainfall with reasonably good accuracy.