超音速气流绕局部透气凹角流动的数值模拟

以具有压力分裂形式的简化Ｎ Ｓ方程为控制方程,数值模拟了超音速来流条件下的激波 边界层干扰被动控制（ｐａｓｓｉｖｅｃｏｎｔｒｏｌｏｆｓｈｏｃｋ ｂｏｕｎｄａｒｙｌａｙｅｒｉｎｔｅｒａｃｔｉｏｎ）。模拟是以预先给定激波前吹气和激波后吸气的流量来实现的。为了定性地确定吹气或吸气对激波 边界层干扰的影响,首先计算了单独吹气和单独吸气两种情况。数值计算时采用了多重扫描法对控制方程差分离散,以反映亚音速区压力对流场的椭圆性影响。     

激波边界层的相互作用对扰动波传播的影响

利用线性稳定性理论和直接数值模拟研究了带有入射斜激波的、来流马赫数Ma=4．5条件下的平板边界层的失稳特性。重点考察了在由于激波边界层相互干扰，平板边界层上形成分离区，又进而产生激波、膨胀波和旋涡等复杂流动现象的流场上游，引入小扰动的TS波后，扰动波传播通过带有这些复杂流动现象的流场时，扰动波的发展变化特点。通过对流场中扰动波（包括基本波和衍生波）演化特征的分析，研究分离和激波等复杂流动现象对平板边界层稳定性的影响特点。数值模拟发现，激波的出现不同于一般的压缩波，在亚声速区与超声速区对扰动波演化的影响是不同的，此外，分离对扰动有稳定作用。     

管内正激波/湍流边界层干扰的数值模拟

本文采用时均Ｎ－Ｓ方程和Ｂａｌｄｗｉｎ／Ｌｏｍａｘ代数湍流模型计算了典型拉伐尔喷管内正激波与湍流边界层干扰流场。计算与实验结果的比较表明，方法可准确地预测激波结构，激波与边界层干扰区流动特征，激波位置、激波前马赫数和壁面压力分布等。     

进气道内激波/边界层干扰及控制研究进展

进气道作为高速航空航天飞行器的重要气动部件,对飞行器的气动力特性、结构重量、隐身性能等有显著影响。激波/边界层干扰现象是高速进气道内普遍存在的一类流动现象,对进气道的性能有突出的影响。发生于进气道内的激波/边界层干扰现象主要可分为正激波/边界层干扰、斜激波/边界层干扰以及三维激波/边界层干扰几类,由于受到侧壁壁面和进气道内背景波系的影响,这些干扰现象偏离了传统基于简化模型的研究结果,具有显著的耦合干扰特征,干扰区间内三维特征明显。概述了发生于进气道内的激波/边界层干扰特性及相关研究进展,并对目前进气道内激波/边界层干扰现象的控制方法进行了总结。     

高马赫数下激波湍流边界层干扰数值研究

应用GAO—YONG可压缩湍流方程组数值模拟了入射斜激波／平板湍流边界层相互干扰现象,计算了来流马赫数为5．0,激波入射角度分别为15．876°、23．287°两种不同激波干扰强度下的流场。计算程序中的对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge—Kutta显式时间推进法求解空间离散后的控制方程。计算较好地模拟了高马赫数下的激波／湍流边界层干扰的流场结构,位移边界层厚度,动量损失厚度等,也比较准确地预测了平板壁面压力、摩阻系数等气动力参数的分布。     

侧板构型对高超侧压进气道起动性能的影响

为了考察侧板构型对高超侧压进气道起动性能的影响,对相同收缩比下侧板分别为前掠和后掠构型的进气道开展了Ma=4来流下的风洞实验及相应的数值模拟研究。实验结果表明,侧板后掠进气道的起动性能优于侧板前掠构型,实验中侧板后掠进气道能够在+2°攻角时实现起动,而侧板前掠进气道仅能在-2°攻角时起动。对流场进行的数值模拟结果表明,侧板后掠进气道不但比侧板前掠进气道具有更高的内收缩段入口马赫数,而且交汇后的侧板激波与底板边界层干扰的强度较弱,使得边界层不易分离,两方面因素共同作用使得侧板后掠进气道的起动性能显著优于侧板前掠构型。     

Sajben跨声速扩压器分离流动中湍流模式数值研究

为研究湍流模式对激波/湍流边界层干扰内流流动的影响,提高数值计算准确度,使用SA,SST k-ω,非线性EASM k-ω,Gao-Yong四个湍流模式对Sajben扩压器内激波/湍流边界层干扰流动进行了数值计算。对流项采用Roe格式离散,扩散项采用二阶中心格式离散,离散后的控制方程用多步Runge-Kutta显示时间推进法求解。文中展示了四个湍流模式计算得到的壁面压力、速度剖面、摩阻系数等分布。计算值与实验值符合很好,四个湍流模式总体上能够较好地模拟扩压器内激波/湍流边界层干扰复杂分离流动。Gao-Yong湍流模式对分离区内的压力、速度型的模拟更加准确,而非线性EASM k-ω模式对分离再附点位置计算最理想。     

进气道中激波-边界层干扰特性研究

针对高超进气道中激波/边界层干扰等现象,选取了24°压缩拐角及二维高超进气道DLR-GK01等算例,比较研究了混合网格并行数值模拟软件中SA、SST、TNT三种湍流模型对激波/边界层干扰等现象的模拟能力,主要考察了不同湍流模型对流场结构、壁面压强以及壁面摩阻的影响。结果表明:三种湍流模型在模拟大分离流场时存在不同的流动分离特性,TNT湍流模型较SST和SA湍流模型模拟的分离起始位置、分离区大小及边界层厚度与实验吻合度最好。     

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

高超声速激波/边界层干扰流动数值模拟研究

在基于Roe格式的全Navier-Stokes方程计算流体力学(CFD)代码中,发展了一种局部熵修正方法,克服了传统熵修正方法在边界层流动模拟中耗散过大的缺点,可用于更加准确的模拟激波/边界层干扰的复杂高超声速流动。对典型高超声速双锥边界层分离与激波干扰的复杂流动进行了模拟,研究了网格密度和熵修正方法耗散性对计算结果的影响。研究表明:高超声速双锥边界层分离与激波干扰流动的数值模拟结果对网格具有很强的敏感性,过稀的网格将产生严重的分离流动预测偏差;低耗散性的局部熵修正方法能更加准确地模拟复杂的高超声速激波与边界层分离流动干扰现象。     

组合发动机可调进气道气动性能

为使涡轮/冲压组合发动机能在宽广的飞行高度、速度内飞行,设计一种组合发动机可调进气道气动结构.编制带黏性修正激波系结构预估程序,实现不同来流马赫数及来流攻角下,进气道内激波系结构布置及强度预估,并指导进气道可调楔板尺寸选取.考虑进气道附面层抽吸后,对来流在马赫数为1.6～4.0范围内,-2°,0°和+2°三种攻角下进气...     

基于LES方法的平板非定常激波/湍流边界层干扰研究

以高超声速发动机进气道湍流分离控制为应用背景,采用大涡模拟(LES)方法进行马赫数为3.0(唇口附近马赫数约为3.0)的激波/湍流边界层干扰(SWTBLI)流场机理研究.利用扰动循环引入的方法,先得到充分发展湍流场,然后根据斜激波关系式引入激波的方法进行激波/湍流干扰模拟.研究结果显示:充分发展湍流场在激波作用下产生逆...     

跨音压气机叶栅的激波结构模型及损失

在分析讨论激波与附面层相互作用和栅后背压引起激波形状及强度变化的基础上, 给出了考虑激波/附面层相互作用及栅后背压的跨音叶栅激波结构的物理数学模型。应用本文所提出的模型分析了跨音叶栅的激波损失, 其结果和实验结果一致。激波损失的精确得出, 使得将激波与附面层相互作用所引起的流动分离损失从流动总损失中分离出来成为可能, 有助于了解激波与附面层相互作用引起流动分离的机理。      

考虑附面层影响的二元混压式进气道设计

采用等激波强度的方法,考虑附面层修正,设计了一种飞行马赫数Ma=3.0的二元混压式进气道.通过数值仿真,模拟了激波-边界层的相互影响,研究了附面层抽吸对内流场的影响,获得了进气道内部复杂的流场分布,以及不同背压下进气道的起动特性.计算表明所设计的进气道性能较好,附面层抽吸对稳定正激波有明显的作用,提高了进气道抗反压能力.给出的方法可用于二元混压式进气道的初步设计和验证.      

超声速压气机叶栅前缘通道激波损失的鼓包控制研究

为了有效减小超声速压气机叶栅变进气马赫数条件下的前缘通道激波损失及由激波诱导的边界层分离,提出了一种带有平直过渡区的新型鼓包结构,并采用数值方法详细分析了新型鼓包结构对激波与激波/边界层相互作用机理以及鼓包几何尺寸与位置对控制效果的影响机制。研究结果表明:新型鼓包在迎风侧凹面产生的压缩波系有效削弱了前缘通道激波的强度,鼓包过渡区产生的膨胀波系使边界层流体加速,明显抑制了局部流动分离,并使分离提前再附。当某一超声速压气机叶栅的前缘通道激波入射在鼓包的过渡区范围内,鼓包高度为0.35倍的边界层厚度且鼓包迎风侧与背风侧长度分别为过渡区长度4倍与5倍时,可以实现较好的控制效果。此外,与无鼓包方案相比,新型鼓包结构可使超声压气机叶栅在设计工况下的总压损失减少4.6%,同时超声速压气机叶栅进气马赫数在1.65~1.8范围内仍能取得较好的气动减损效果。      

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.     

Reflected shock/turbulent boundary layer interaction structure analysis based on large eddy simulation

In order to understand the physical phenomenon of the reflected shock/turbulent boundary layer interaction, the Large Eddy Simulation (LES) is conducted to investigate shock wave and turbulent boundary layer interaction in a $12°$ compression ramp with inlet high Mach number of 2.9. Rescaling/recycling method is used as inflow turbulence generation technique and validated on a supersonic flat plate turbulent boundary layer. The flow field of recycling plane in the plate computation domain is obtained to give the inlet boundary condition for the LES computation. This paper focuses on the reflected shock/turbulent boundary layer interaction region, where the fine flow structure and instantaneous flow field are analyzed in detail. It is found that the unsteady motion of the shock wave leads to the increase of wall pressure fluctuation.     

高超声速内转式进气道型面流场重构

内转式进气道流场参数分布不均，为改善该类进气道的气动性能，本文采用数值仿真方法开展了基于型面的内转式进气道流场重构研究。结果表明：流场重构型面中心线对进气道流场结构及流动特性影响较大，在给定偏距/长度与中心线末端斜率的约束条件下，选取合适的中心线起始角能够大幅提高进气道的气动性能，改善流场参数分布。与进气道原型方案相比，流场重构型面中心线10°起始角的进气道方案总压恢复系数、抗反压能力最大分别提升33.7%、26.4%，自起动马赫数下降1.1。随着流场重构型面中心线起始角增大，唇罩激波/侧壁边界层干扰诱发的流向涡减弱、流向涡传输轨迹向唇罩一侧偏移，低能流向唇罩两侧迁移趋势增强。在研究范围内，随着流场重构型面中心线起始角增大，隔离段出口总压恢复系数先增大后减小，自起动马赫数先下降后不变。     

Numerical study of unsteady starting characteristics of a hypersonic inlet

The impulse and self starting characteristics of a mixed-compression hypersonic inlet designed at Mach number of 6.5 are studied by applying the unsteady computational fluid dynamics (CFD) method. The full Navier-Stokes equations are solved with the assumption of viscous perfect gas model, and the shear-stress transport (SST) k-x two-equation Reynolds averaged Navier- Stokes (RANS) model is used for turbulence modeling. Results indicate that during impulse starting, the flow field is divided into three zones with different aerodynamic parameters by primary shock and upstream-facing shock. The separation bubble on the shoulder of ramp undergoes a generating, growing, swallowing and disappearing process in sequence. But a separation bubble at the entrance of inlet exists until the freestream velocity is accelerated to the starting Mach number during self starting. The mass flux distribution of flow field is non-uniform because of the interaction between shock and boundary layer, so that the mass flow rate at throat is unsteady during impulse starting. The duration of impulse starting process increases almost linearly with the decrease of freestream Mach number but rises abruptly when the freestream Mach number approaches the starting Mach number. The accelerating performance of booster almost has no influence on the self starting ability of hypersonic inlet.     

乘波前体两侧高超声速内收缩进气道一体化设计

为了探索两侧进气系统的流场结构及气动性能,采用吻切锥乘波前体、压升规律可控的一种高超声速内收缩进气道设计了两侧进气布局的高超声速飞行器一体化进气系统,并进行了数值模拟,研究了进气系统的流场结构、速度特性、攻角特性以及侧滑角特性等。结果表明,设计点前体外流场和进气道内流场相互独立,接力点前体前缘激波和进气道前缘激波相互耦合。由于未吞入前体附面层,因而进气道内激波附面层相互作用较弱,没有产生分离;随来流马赫数增大,进气道总压恢复系数减小,增压比增大显著,升阻比几乎不变;随攻角增大,流量系数增大明显,总压恢复系数略有减小,增压比增大明显,升阻比逐渐增大;随侧滑角增大,进气道总体性能逐渐减小,迎风侧进气道性能下降较小,背风侧进气道性能下降明显。