火星表面热环境对航天器热控影响分析

火星大气对太阳辐射产生吸收和散射作用，同时还将与火星表面航天器发生对流换热。热设计时难以直接评估对流、辐射和导热三种换热对航天器的影响，从而确定主要的控温途径。在调研火星表面辐射、大气等热环境的基础上，从线性化传热系数和对流辐射比的角度对比分析了辐射、对流和导热对航天器的影响。器表辐射传热系数随光学属性和温度的变化范围为0.3～1.4W/(m2·℃)，对流传热系数随风速变化为0.2～1.5W/(m2·℃)，器内导热传热系数可控制在0.25W/(m2·℃)以下。结果表明，太阳辐射较火星表面和天空辐射而言是主要外热源，航天器表面的辐射和对流换热为两条并联换热途径，两者均可成为主要换热途径，器内导热传热是控制航天器内外隔热的主要可控因素。     

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.     

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

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

Research on navigation of satellite constellation based on an asynchronous observation model using X-ray pulsar

Pulsar navigation is a promising navigation method for high-altitude orbit space tasks or deep space exploration. At present, an important reason for restricting the development of pulsar navigation is that navigation accuracy is not high due to the slow update of the measurements. In order to improve the accuracy of pulsar navigation, an asynchronous observation model which can improve the update rate of the measurements is proposed on the basis of satellite constellation which has a broad space for development because of its visibility and reliability. The simulation results show that the asynchronous observation model improves the positioning accuracy by 31.48% and velocity accuracy by 24.75% than that of the synchronous observation model. With the new Doppler effects compensation method in the asynchronous observation model proposed in this paper, the positioning accuracy is improved by 32.27%, and the velocity accuracy is improved by 34.07% than that of the traditional method. The simulation results show that without considering the clock error will result in a filtering divergence.     

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.     

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

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

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.     

飞机封闭腔室减重降噪优化方法

噪声水平是现代飞机封闭腔室设计的重要指标,设计中单纯地对壁板进行加肋处理虽能显著降低噪声但会大大增加结构质量,为此对由封闭腔室构成的结构-声耦合系统进行了减重降噪优化研究。基于结构-声耦合有限元模型,利用有限元软件ACTRAN计算了频谱加载时舱内的声压响应。通过试验对简化处理及数值计算进行了验证,并修正了相关模型参数。为了降低结构质量,以加强肋为边界对舱门壁板进行了分区,通过对各个区域壁板厚度及肋条截面积的优化设计,使系统动刚度分配更趋合理,降低了声辐射能量以及结构-声腔的耦合性,从而实现了在满足噪声约束条件下减轻结构质量的目标。本文的工作对实际工程中由加肋壁板所构成的类似结构的减重降噪设计有着较好的工程指导价值。     

Numerical investigation on flow nonuniformity-induced hysteresis in scramjet isolator

Numerical simulation and theoretical analysis were conducted to study the hysteresis inside scramjet isolator during the reciprocating process of back pressure variation. It is revealed that only a regular reflection is theoretically possible for two leading shocks when the inflow Mach number is greater than 2.0, and no hysteresis can occur in the transition between shock reflection types. Nevertheless, wall suction, gas injection, and background waves cause non-uniformity of the incoming flow and would make hysteresis possible. Besides the classical hysteresis in the transition between shock reflection, new kinds of hysteresis were found in both the deflection angle of separated boundary layer and the location of the shock train. Moreover, the occurrence of hysteresis in the deflection angle of the separated boundary layer is accompanied with the shock reflection hysteresis. In the case with background waves or gas injection, hysteresis in the starting position of leading shock was observed too. As back pressure decreases, the leading shock does not follow the same path as that as the back pressure increases, and it is anchored at the location where the background shock or the injection interacts with the leading shock. It is inferred that, if two strong adverse pressure gradient regions move towards and interact with each other, hysteresis will take place when they start to separate.     

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.