Design method with controllable velocity direction at throat for inward-turning inlets

In the design of a hypersonic inward-turning inlet by applying the traditional basic flowfield, a reflected shock-wave is formed in the isolator due to the continuous reflection of the cowlreflected shock wave in the basic flow-field, which interacts with the boundary layer to produce a considerable influence on the performance of the inlet. Here, a basic flow-field design method that can control the velocity direction at the throat section is developed, and numerical simulations are conducted to demonstrate the effectiveness of this method. The method presented in this paper can achieve the absorption of the reflected waves at the shoulder of the basic flow-field by adjusting the variation law of the center radius in the basic flow-field, and a smooth transition between the compression surface and the isolator can also be produced. The Mach number and total pressure recovery coefficient of the inlet designed according to this method are 3.00 and 0.657, respectively, at design point(the incoming flow Mach number Ma1= 6.0). The results show that with this method, the inlet can efficiently weaken both the reflection of the shock wave and the interaction between the boundary layer and the reflected shock waves, which improves the aerodynamic performance of the inlet.     

基于描述函数的燃烧振荡建模与仿真

针对贫预混燃烧在燃气轮机燃烧室中的燃烧振荡问题进行主动控制技术研究，一般需依赖于较为准确的火焰动力学模型，但是通过仿真或者试验测量得到完整、准确的非线性火焰动力学响应特性都存在较大的困难。针对该问题，本文基于声学分析，结合非线性建模方法的应用，提出了一种基于描述函数的热声耦合系统分析方法，既能够体现火焰非线性动力学过程的影响，又便于实验测量与系统辨识。本文针对旋流预混燃烧室的热声耦合特性进行了实验测量，通过系统辨识得到热声耦合系统的描述函数模型，在此基础上开展热声耦合系统时域仿真分析。实验测量结果表明，包含火焰动力学过程的声波传播过程具有较强的非线性特征，并且在一定频率范围内描述函数的幅值和相位随扰动量幅值变化规律体现出较一致的规律；仿真分析结果证明了辨识得到的描述函数可以对燃烧室中的热声耦合形成过程和极限环振荡特征进行描述。     

Shock/shock interactions between bodies and wings

This paper examines the Shock/Shock Interactions (SSI) between the body and wing of aircraft in supersonic flows. The body is simplified to a flat wedge and the wing is assumed to be a sharp wing. The theoretical spatial dimension reduction method, which transforms the 3D problem into a 2D one, is used to analyze the SSI between the body and wing. The temperature and pressure behind the Mach stem induced by the wing and body are obtained, and the wave configurations in the corner are determined. Numerical validations are conducted by solving the inviscid Euler equations in 3D with a Non-oscillatory and Non-free-parameters Dissipative (NND) finite difference scheme. Good agreements between the theoretical and numerical results are obtained. Additionally, the effects of the wedge angle and sweep angle on wave configurations and flow field are considered numerically and theoretically. The influences of wedge angle are significant, whereas the effects of sweep angle on wave configurations are negligible. This paper provides useful information for the design and thermal protection of aircraft in supersonic and hypersonic flows.     

Flow characteristics around airfoils near transonic buffet onset conditions

In transonic flow, buffet is a phenomenon of flow instability caused by shock wave/boundary layer interaction and flow separation. The phenomenon is common in transonic flow, and it has serious impact on the structural strength and fatigue life of aircraft. In this paper, three typical airfoils: the supercritical OAT15A, the high-speed symmetrical NACA64A010, and the thin, transonic/supersonic NACA64A204 are selected as the research objects. The flow fields of these airfoils under pre-buffet and buffet onset conditions are simulated by Unsteady Reynolds Averaged Navier-Stokes (URANS) method, and the mode analysis of numerical results is carried out by Dynamic Mode Decomposition (DMD). Qualitative and quantitative analysis of the shock wave motion, shock wave intensity, shock foot bubble and trailing edge separation, and pressure coefficient fluctuation were performed to attain deep insight of transonic buffet flow features of different airfoils near buffet onset conditions. The results of DMD analysis show that the energy proportion of the steady mode of these airfoils decreases dramatically when approaching the buffet onset angle of attack, while the growth rate of the primary mode increases inversely. It was found that at the onset of buffet, there exist different degrees of merging behavior between shock foot bubble and trailing edge separation during one buffet cycle, and the instability of shock wave and separation induced shear layer are closely related to the merging behavior.     

Ⅱ型PIO反馈与前馈抑制系统对比

 当舵机速率限制饱和时会产生附加时延,恶化飞行品质,触发Ⅱ型驾驶员诱发振荡（PIO）。运用描述函数法研究速率限制反馈（RLF）和死区增稳（DASA）Ⅱ型PIO抑制系统的开环相位补偿能力,提出RLF抑制系统的参数优化选择方法,分析了DASA抑制系统中死区和滤波器参数对相位补偿能力的影响。研究了两型抑制系统对正弦与偏差输入信号的开环响应特点。建立了含速率限制器的人机系统数学模型,基于描述函数法探究了Ⅱ型PIO的产生机理,推导了抑制Ⅱ型PIO发生的公式,应用连续信号相位补偿法抑制Ⅱ型PIO的发生。对阶跃、离散和正弦3种易于诱发PIO现象的跟踪任务进行数值仿真,对比了两种抑制系统的Ⅱ型PIO抑制能力。结果表明,反馈型RLF抑制系统对Ⅱ型PIO的抑制能力强于前馈型DASA抑制系统,有效降低了Ⅱ型PIO发生的可能性。     

对称交叉激波和湍流边界层相互作用的数值研究

针对两种双尖鳍外形的对称交叉激波与湍流边界层相互作用,采用N-S方程和两种湍流模型进行了计算。研究了网格收敛性、鳍的角度和湍流模型对壁面压强、Stanton数和壁面摩擦力线的影响。弱相互作用的计算结果较好,强相互作用的壁面压强和摩擦力线的计算结果与试验吻合较好,而Stanton数的计算结果较差,峰值高达试验的2.5倍左右。随着鳍的角度的增加,壁面压强和Stanton数的分布从单调分布发展为M型分布,两者的峰值不在相同的位置。湍流模型对壁面压强和壁面摩擦力线影响很小,对Stanton数计算的影响很大,SST模型比BSL模型表现好一些。     

Heating heavy ions in the polar corona by collisionless shocks: A one-dimensional simulation

Recently a new model for explaining the observations of preferential heating of heavy ions in the polar solar corona was proposed ( and ). In that model the ion energization mechanism is the ion reflection off supercritical quasi-perpendicular collisionless shocks in the corona and the subsequent acceleration by the motional electric field E = −V × B/c. The mechanism of heavy ion reflection is based on ion gyration in the magnetic overshoot of the shock. The acceleration due to the motional electric field is perpendicular to the magnetic field, giving rise to large temperature anisotropy with T_⊥ ? T_∥, in agreement with SoHO observations. Such a model is tested here by means of a one dimensional test particle simulation where ions are launched toward electric and magnetic profiles representing the shock transition. We study the dynamics of O⁵⁺, as representative of coronal heavy ions for Alfvénic Mach numbers of 2–4, as appropriate to solar corona. It is found that O⁵⁺ ions are easily reflected and gain more than mass proportional energy with respect to protons.     

轴对称后向台阶不稳定流动及压强振荡数值研究

使用翼柱形装药的战术发动机,在翼槽烧完时会出现轴对称的后向台阶结构,易产生转角涡脱落引起的声涡耦合现象。根据战术发动机中后向台阶流动特点设计简化模型,从线性稳定性理论出发,应用声涡耦合模态预测方法进行分析。采用湍流流动的大涡模拟(LES)方法,进行了较宽来流Re数条件下后向台阶不稳定流动的数值模拟,获得了转角涡脱落产生的一般过程,分析了不同来流条件下的旋涡多尺度运动规律。通过分析流场中关键位置的压强振荡,研究了气动声学现象。结果表明,声涡耦合模态预测方法可对不同速度下可能诱发的声场固有频率进行合理预测。流场头部压强波动的振幅最大,频率与声场轴向声模态相耦合。再发展区出现低于轴向基频的宽峰振荡主要受大尺度旋涡运动的影响,该频率与基频接近时,会造成基频的多峰现象。当近壁面小尺寸旋涡形成频率接近声场某阶固有频率时,该频率及更低阶频率就有可能被激发。     

三维流场乘波体快速设计方法及多目标优化

为扩大乘波体的设计空间,精确地计算设计流场,引入高超声速激波装配法,结合基于多重网格的搜索技术加速寻点效率,扩展流线追踪方法到三流流场,建立普适于一般三维流场的乘波体快速设计方法。对设计曲线使用B样条方法建模并提取设计变量,结合改进的Pareto遗传算法对升阻比和容积率进行多目标优化设计,获得了性能较高的乘波外形。研究表明此方法效率高,设计空间大,是一种乘波体设计的新方法和新思路。     

前向空腔的弓形激波特性研究

在FD-14A激波风洞中Ma=10流场对前向空腔构型开展试验研究,应用高速阴影技术捕捉弓形激波的平均位置及振荡幅值,利用压力传感器测量空腔底部的脉动压力。在现有无空腔钝头体激波脱体距离预测方法的基础上发展了前向空腔构型的激波脱体距离预测方法,结合国外的试验测量结果与Organ-pipe理论,验证了这种方法的有效性和适用性,且该方法对激波脱体距离的预测结果与FD-14A风洞试验结果一致。此外,基于这种方法讨论了空腔振荡频率预测方法存在的争议。最后,研究了Ma=10流场下球锥体-前向空腔构型的脱体激波振荡幅值与平均速度的规律。