涡轮叶栅外换热系数计算

采用求解三维雷诺平均Navier-Stokes(N-S)方程和带转捩模型的二方程Shear stress transport(SST)湍流模型,完成了MARKⅡ,VKI两个高压涡轮导叶叶栅及一个VKI高压涡轮动叶叶栅的外换热系数计算.计算结果与试验数据的对比表明,发展的带转捩模型的数值方法明显地提高了外换热系数的计算精度.      

Halo轨道维持的线性周期控制策略

共线型平动点附近的Halo轨道具有指数不稳定性,轨道维持是必不可少的。推导了基于Halo轨道的误差动力学方程,并证明其一阶近似即为线性周期系统。以一次维持的作用时间为离散步长,并通过定常变换,将所得误差动力学化为线性离散定常系统;则仅需通过极点配置,即可实现Halo轨道镇定。研究结果表明,利用Halo轨道周期性设计的线性周期控制策略,可以满足轨道维持任务的需要。     

地月Halo与DRO支持的往返月球任务轨道

面向未来载人月球与深空探测任务需求,针对地月系统中低能往返环月轨道与平动点轨道的转移设计问题进行研究,在不同三体轨道下系统分析了环月轨道变化对任务飞行时间与燃料消耗等关键参数的影响,并提出月球往返轨道初值猜想搜索策略。为解决设计变量初值敏感性问题,结合空间不变流形与目标平动点轨道构型特性,采用微分修正算法快速构造初始转移轨迹。在同时考虑近月点与燃料最优转移多约束条件下,通过多重配点打靶法与序列二次规划算法对环月轨道与平动点轨道间往返轨迹进一步研究,并推导了约束方程的梯度公式提高计算效率。为分析往返轨道特性与验证设计策略的有效性,针对不同环月轨道倾角、三体轨道幅值等参数变化与转移时间、燃耗的关系进行分析,设计结果对探月飞行器近月空间部署往返轨道设计及参数选取具有重要的参考意义。     

A new stochastic model considering satellite clock interpolation errors in precise point positioning

Precise clock products are typically interpolated based on the sampling interval of the observational data when they are used for in precise point positioning. However, due to the occurrence of white noise in atomic clocks, a residual component of such noise will inevitable reside within the observations when clock errors are interpolated, and such noise will affect the resolution of the positioning results. In this paper, which is based on a twenty-one-week analysis of the atomic clock noise characteristics of numerous satellites, a new stochastic observation model that considers satellite clock interpolation errors is proposed. First, the systematic error of each satellite in the IGR clock product was extracted using a wavelet de-noising method to obtain the empirical characteristics of atomic clock noise within each clock product. Then, based on those empirical characteristics, a stochastic observation model was structured that considered the satellite clock interpolation errors. Subsequently, the IGR and IGS clock products at different time intervals were used for experimental validation. A verification using 179 stations worldwide from the IGS showed that, compared with the conventional model, the convergence times using the stochastic model proposed in this study were respectively shortened by 4.8% and 4.0% when the IGR and IGS 300-s-interval clock products were used and by 19.1% and 19.4% when the 900-s-interval clock products were used. Furthermore, the disturbances during the initial phase of the calculation were also effectively improved.     

高超声速流动中噪声与湍流度的关系

在不可压缩流动中，湍流度对流动失稳与转捩的影响机制有很多研究成果，已发展了一些依据其大小预测流动转捩的模型应用于工程设计。但在高速可压缩流动中，湍流度很难测量，给该速域的流动转捩预测建模带来困难。近年来的研究表明，高速流场（马赫数≥3）中的噪声是影响流动转捩提前与推迟，甚至决定转捩途径的关键因素。但流场噪声与湍流度对流场的激励机制有没有内在联系，或流场噪声与湍流度理论上是否存在相互关联，迄今仍缺少理论分析结果。由于可压缩流场中噪声的测量相对容易，因此探讨流场噪声与湍流度之间的定量关系具有重要的理论意义和实用价值。讨论了Euler系统中各类双曲波特性，重点分析了流场声波中脉动压力与脉动速度的本质联系，从理论上得到了声压级与声致湍流度的关系。同时，为了对影响湍流度的因素有更全面的认识，导出了非声扰动与非声致湍流度的关系。为建立高速流动转捩预测模型和探讨天地一致性问题，奠定了一定的研究基础。     

Satellite availability and positioning performance of uncombined precise point positioning using BeiDou-2 and BeiDou-3 multi-frequency signals

To realize the smooth transition from regional BeiDou Navigation Satellite System (BDS-2) to the global one (BDS-3), the integration of BDS-2 and BDS-3 is important for providing continuous, stable and reliable positioning, navigation and timing (PNT) services for global users. This work used 154 globally distributed multi-GNSS (Global Navigation Satellite System) experiment stations spanning 30 days to analyze the satellite availability and positioning performance of uncombined precise point positioning (UC-PPP) under current BDS-2 and BDS-3 constellations. We focused on three issues: the influence of BDS-3 receiver tracking abilities, the positioning performance among different areas, and the benefit of multi-frequency observations. The results show that the elliptical zone caused by poor BDS-2 satellite visibility is disappeared when the evenly distributed BDS-3 medium earth orbit satellites are introduced. When BDS-3 are integrated with BDS-2, the area with the Position Dilution of Precision (PDOP) less than 2 can be expanded to 75° S-75° N and 30° E-150° W. The positioning performance of BDS-3 and BDS-2/BDS-3 UC-PPP are seriously affected by the receiver tracking abilities of BDS-3 signals. When the maximum pseudo-random noise sequences (PRNs) of BDS-3 satellites tracked by stations are within 30 or 37, the positioning accuracy of static UC-PPP can be improved by 22.94% or 8.27% due to the integration of BDS-2 and BDS-3. Besides, the most improvement of BDS-2 and BDS-3 integration is achieved in Asia-Pacific regions, especially for the kinematic UC-PPP or the poor receiver tracking abilities of BDS-3. Similar to the multi-frequency BDS-2 UC-PPP, the benefit of multi-frequency signals for BDS-3 or BDS-2/BDS-3 UC-PPP is also non-vital. The three-dimensional positioning accuracy of BDS-2/BDS-3 multi-frequency UC-PPP in static mode and kinematic mode are 2.24 cm and 5.39 cm, while the corresponding convergence time are 49.62 min and 73.80 min, respectively. Compared with BDS-2, both the positioning accuracy and the convergence time of BDS-2/BDS-3 joint UC-PPP are improved by approximately over 50%, which indicates that BDS-3 has a great potential to provide high-quality PNT services as other global navigation satellite systems.     

高超声速气动加热关联方法的适应性分析

针对高超声速飞行器热环境预测上有着重要应用的关联方法缺乏构造机制和适应性分析的问题，从驻点热流的Kemp-Riddell公式出发，推导了驻点斯坦顿数和雷诺数之间的经典关联关系，从理论上分析了上述关系的适应性。结果显示斯坦顿数和雷诺数关联关系受马赫数的制约，飞行条件和试验条件马赫数变化范围较大的情况下，更适合采用雷诺数和马赫数的组合变量开展关联分析。当试验来流静焓较低的情况下，过低静焓对热环境的关联特性有不利影响，需要提高试验来流马赫数，远离临界马赫数确保地面数据同飞行条件之间的关联性。在国内4座风洞上开展的风洞试验和DSMC/CFD的数值仿真分析工作，校验了上述结论在连续流区域的正确性。     

Thermal behavior of an isolator with mode transition inducing back-pressure of a dual-mode scramjet

Combustion mode transition is a valuable and challenging research area in dual-mode scramjet engines. The thermal behavior of an isolator with mode transition inducing back-pressure is investigated by direct-connect dual-mode scramjet experiments and theoretical analysis. Combustion experiments are conducted under the incoming airflow conditions of total temperature 1270 K and Mach 2. A small increment of the fuel equivalence ratio is scheduled to trigger mode transition. Correspondingly, the variation of the coolant flow rate is very small. Based on the mea-sured wall pressures, the heat-transfer model can quantify the thermal state variation of the engine with active cooling. Compared with the combustor, mode transition has a greater effect on the iso-lator thermal behavior, and it significantly changes the isolator heat-flux and wall temperature. To further study the isolator thermal behavior from flight Mach 4 to Mach 7, a theoretical analysis is carried out. Around the critical point of combustion mode transition, sudden changes of the isola-tor flowfield and thermal state are discussed.     

Contribution analysis of QZSS to single-frequency PPP of GPS/BDS/GLONASS/Galileo

The Quasi-Zenith Satellite System (QZSS) established by the Japan Aerospace Exploration Agency mainly serves the Asia-Pacific region and its surrounding areas. Currently, four in-orbit satellites provide services. Most users of GNSS in the mass market use single-frequency (SF) receivers owing to the low cost. Therefore, it is meaningful to analyze and evaluate the contribution of the QZSS to SF precise point positioning (PPP) of GPS/BDS/GLONASS/Galileo systems with the emergence of GNSS and QZSS. This study compares the performances of three SF PPP models, namely the GRoup and PHase Ionospheric Correction (GRAPHIC) model, GRAPHIC with code observation model, and an ionosphere-constrained model, and evaluated the contribution of the QZSS to the SF PPP of GPS/BDS/GLONASS/Galileo systems. Moreover, the influence of code bias on the SF PPP of the BDS system is also analyzed. A two-week dataset (DOY 013–026, 2019) from 10 stations of the MGEX network is selected for validation, and the results show that: (1) For cut-off elevation angles of 15, 20, and 25°, the convergence times for the static SF PPP of GLONASS + QZSS are reduced by 4.3, 30.8, and 12.7%, respectively, and the positioning accuracy is similar compared with that of the GLONASS system. Compared with the BDS single system, the convergence times for the static SF PPP of BDS + QZSS under 15 and 25° are reduced by 37.6 and 39.2%, the horizontal positioning accuracies are improved by 18.6 and 14.1%, and the vertical components are improved by 13.9 and 21.4%, respectively. At cut-off elevation angles of 15, 20, and 25°, the positioning accuracy and precision of GPS/BDS/GLONASS/Galileo + QZSS is similar to that of GPS/BDS/GLONASS/Galileo. And the convergence times are reduced by 7.4 and 4.3% at cut-off elevation angles of 20 and 25°, respectively. In imitating dynamic PPP, the QZSS significantly improves the positioning accuracy of BDS and GLONASS. However, QZSS has little effect on the GPS-only, Galileo-only and GPS/BDS/GLONASS/Galileo. (2) The code bias of BDS IGSO and MEO cannot be ignored in SF PPP. In static SF PPP, taking the frequency band of B1I whose multipath combination is the largest among the frequency bands as an example, the vertical component has a systematic bias of approximately 0.4–1.0 m. After correcting the code bias, the positioning error in the vertical component is lower than 0.2 m, and the positioning accuracy in the horizontal component are improved accordingly. (3) The SF PPP model with ionosphere constraints has a better convergence speed, while the positioning accuracy of the three models is nearly equal. Therefore the GRAPHIC model can be used to get good positioning accuracy in the absence of external ionosphere products, but its convergence speed is slower.     

The update of BDS-2 TGD and its impact on positioning

Timing group delay (TGD) is an important parameter that affects the positioning performance of global navigation satellite systems (GNSS). The BeiDou navigation satellite system (BDS) broadcasts TGD corrections from B3I frequency to B1I and B2I frequencies, namely TGD₁ and TGD₂. On July 21, 2017, BDS updated TGD values with a maximum change of more than 4 ns. In this contribution, we explain the motivation for the BDS TGD update, which is due to the systematic bias between narrowly correlated and widely correlated pseudo-ranges in BDS monitoring receivers. To investigate the impact of the updated TGD, BDS signal-in-space range error (SISRE) and user positioning performance regarding single point positioning (SPP) and precise point positioning (PPP) are analyzed. Results show that after the update of TGD, the difference between the new TGD and multi-GNSS experiment (MGEX) differential code bias (DCB) decreases from 1.38 ns to 0.29 ns on TGD₁ and from 0.40 ns to 0.25 ns on TGD₂. With the contribution of more accurate TGD, the systematic bias of BDS radial SISRE no longer exists, and the overall BDS SISRE also reduces from 1.33 m to 0.87 m on B1I/B2I frequency, from 1.05 m to 0.89 m on B1I frequency, from 0.92 m to 0.91 m on B2I frequency, respectively, which proves the similar precision of BDS TGD and MGEX DCB. One week of statistical results from 28 globally distributed MGEX stations shows that the SPP performance improves on non-B3I frequencies after the TGD update, with a maximum improvement of more than 22% for the B1I/B2I or B1I/B3I combination. The new TGD mainly reduces SPP positioning bias in the East component. The updated TGD also slightly improves the PPP convergence performance for the B1I/B3I combination, but mostly contributes to a more accurate estimation of the receiver clock and ambiguities.