超声速边界层三维扰动引起小激波的数值研究

通过直接数值模拟的方法,研究了M=4.5的超声速边界层中三维扰动的演化。以某一剖面为入口,加入一个及一对三维T-S波,发现随展向波数的增加,扰动幅值的增长率逐渐减小,证实了M=4.5的超声速边界层中,当三维扰动达到一定幅值时会有小激波出现,为建立可压缩流动稳定性理论提供了依据。     

非轴对称静子对压气机气动性能的影响

为降低进气畸变对压气机气动性能的影响,设计了 1种压气机非轴对称静子,并对设计方案开展数值模拟研究。仿真结果表明:在最高效率工况下,非轴对称静子能减小畸变区静叶的流动分离,缩小叶尖低密流区域,提升通道的流通能力,压气机的最高效率约增加 0.46%,此外,畸变区叶片进口气流角得到改善,在 90%叶高处的峰值气流角降低 2.5°;在近失速工况下,非轴对称静子能降低畸变区静叶上半叶高的扩压因子,缩小分离范围,虽略微恶化了叶根区域流场,但压气机整体气动性能与流通能力有所提升,能够在更低的流量下工作,稳定裕度增加 31.5%。     

阳极层霍尔推力器的阳极分段形式对内磁极刻蚀速率的影响

由于推力器内径向电势和等离子体的分布会影响轰击推力器器壁离子的运动,本文通过阳极的不同分段形式来改变推力器内电势分布,进而研究推力器阳极的不同分段形式对内磁极的刻蚀速率分布。首先利用PIC方法研究入射离子的数量和能量分布,接着对内磁极内表面和上表面的刻蚀速率进行仿真计算。由结果可知：入射到内磁极两个面上的离子数量和能量,都是在三分段阳极时达到最大值,其他三种情况偏低。对于内磁极内表面和上表面的刻蚀速率也是在三分段阳极时最大,其次是单阳极,而二、四分段阳极时的刻蚀速率最小。最后通过单阳极和二分段阳极时刻蚀速率的试验结果对比,验证了此研究结果的正确性,为下一步提高推力器的寿命研究提供了参考数据和研究方法。     

Strong interactions of incident shock wave with boundary layer along compression corner

The coherent structure and instability of the interaction of incident shock wave with boundary layer developing on a compression corner are experimentally studied. The experiments are carried out in a supersonic wind tunnel of Mach number 2. Particular attention is paid to shock patterns and unsteady shock motions induced by the separation bubble. The high-speed schlieren is used to visualize the flowfield evolution and to characterize the instability. The snapshot proper orthogonal decomposition of schlieren sequences is applied to investigate the primary coherent structure in the flowfield. Fast Fourier transform and continuous wavelet transformation are applied to characterize the instability. The results show that there are large-scale low-frequency oscillations of the shock waves and small-scale high-frequency pulsations in the separation region. The peak frequency of shock oscillation is mainly concentrated in the range of 100–1000 Hz. The pulsation of the small flow structure in the separation bubble is mainly concentrated above 12.5 kHz. Based on the results of experimental analysis, the preliminary mechanism of the large-scale instability of such interaction is obtained.     

A regional model for the prediction of M(3000)F2 over East Asia

Research on empirical or physical models of ionospheric parameters is one of the important topics in the field of space weather and communication support services. To improve the accuracy of predicting the monthly median ionospheric propagating factor at 3000 km of the F2 layer (identified as M(3000)F2) for high frequency radio wave propagation, a model based on modified orthogonal temporal–spatial functions is proposed. The proposed model has three new characteristics: (1) The solar activity parameters of sunspot number and the 10.7-cm solar radio flux are together introduced into temporal reconstruction. (2) Both the geomagnetic dip and its modified value are chosen as features of the geographical spatial variation for spatial reconstruction. (3) A series of harmonic functions are used to represent the M(3000)F2, which reflects seasonal and solar cycle variations. The proposed model is established by combining nonlinear regression for three characteristics with harmonic analysis by using vertical sounding data over East Asia. Statistical results reveal that M(3000)F2 calculated by the proposed model is consistent with the trend of the monthly median observations. The proposed model is better than the International Reference Ionosphere (IRI) model by comparison between predictions and observations of six station, which illustrates that the proposed model outperforms the IRI model over East Asia. The proposed method can be further expanded for potentially providing more accurate predictions for other ionospheric parameters on the global scale.     

Integrated supersonic wind tunnel nozzle

In supersonic wind tunnels, the airflow at the exit of a convergent-divergent nozzle is affected by the connection between the nozzle and test section, because the connection is a source of disturbance for supersonic flow and the source of disturbance generated by this disturbance propagates downstream. In order to avoid the disturbance, the test can only be carried out in the rhombus area. However, for the supersonic nozzle, the rhombus region is small, limiting the size and attitude angle of the test model. An integrated supersonic nozzle is a nozzle and a test section as a whole, which is designed to weaken or eliminate the disturbance. The inviscid contour of the supersonic nozzle is based on the method of characteristics. A new curve is formed by the smooth connection between the inviscid contour and test section, and the boundary layer is corrected for the overall curve. Integrated supersonic nozzles with Mach number 1.5 and 2 are designed, which are based on this method. The flow field is validated by numerical and experimental results. The results of the study highlight the importance of the connection about the nozzle outlet and test section. They clearly show that the wave system does not exist at the exit of the supersonic nozzle, and the flow field is uniform throughout the test section.     

Influence of non-equilibrium reactions on the optimization of aerothrust aeroassisted maneuver with orbital change

The spaceplane is perspective vehicle due to wide maneuverability in comparison with a space capsule. Its maneuverability is expressed by the larger flight range and also by a possibility to rotate orbital inclination in the atmosphere by the aerodynamic and thrust forces. Orbital plane atmospheric rotation maneuvers can significantly reduce fuel costs compared to rocket-dynamic non-coplanar maneuver. However, this maneuver occurs at Mach numbers about 25, and such velocities lead to non-equilibrium chemical reactions in the shock wave. Such reactions change a physicochemical air property, and it affects aerodynamic coefficients. This paper investigates the influence of non-equilibrium reactions on the aerothrust aeroassisted maneuver with orbital change. The approach is to solve an optimization problem using the differential evolution algorithm with a temperature limitation. The spaceplane aerodynamic coefficients are determined by the numerical solution of the Reynolds-averaged Navier-Stokes equations. The aerodynamic calculations are conducted for the cases of perfect and non-equilibrium gases. A comparison of optimal trajectories, control laws, and fuel costs is made between models of perfect and non-equilibrium gases. The effect of a chemically reacting gas on the finite parameters is also evaluated using control laws obtained for a perfect gas.     

Highly efficient computation method for hazard quantification of uncontained rotor failure

Uncontained Engine Rotor Failure (UERF) can cause a catastrophic failure of an aircraft, and the quantitative assessment of the hazards related to UERF is a very important part of safety analysis. However, the procedure for hazard quantification of UERF recommended by the Federal Aviation Administration in advisory circular AC20-128A is cumbersome, as it involves building auxiliary lines and curve projections. To improve the efficiency and general applicability of the risk angle calculation, a boundary discretization method is developed that involves discretizing the geometry of the target part/structure into node points and calculating the risk angles numerically by iterating a particular algorithm over each node point. The improved efficiency and excellent accuracy for the developed algorithm was validated through a comparison with manual solutions for the hazard quantification of the engine nacelle structures of a passenger aircraft using the guidance in AC20-128A. To further demonstrate the applicability of the boundary discretization method, the proposed algorithm was used to examine the influence of the target size and the distance between the target and rotor on the hazard probability.     

高超声速有攻角锥飞行工况下迎风面波包的演化

飞行工况下由于壁温与来流温度之比较低，Mack模态的不稳定性会得到显著增强，因此Mack模态占主导的有攻角锥迎风面相对侧面可能会提前转捩。本文采用高分辨率直接数值模拟研究了高超声速有攻角锥在飞行工况下迎风面Mack模态的演化规律，Mack模态由迎风中心线附近的一个短时局部壁面吹吸激发。波包的空间分布和不同模态幅值沿流向的演化过程表明在有攻角条件下，Mack模态的演化过程与零度攻角锥边界层中的类似，基频共振是Mack模态最可能的转捩形式。