前缘钝化和化学非平衡效应对斜爆震发动机进气道特性的影响

为了探究前缘钝化、化学非平衡效应对斜爆震发动机进气道的工作性能及出口温度边界层分布的影响,采用热完全气体、化学非平衡气体两种模型对不同顶板、唇口前缘钝化半径下斜爆震发动机进气道进行数值模拟,结果表明：相比基准进气道,钝化后进气道上壁面温度边界层厚度较大,下壁面温度边界层厚度较小;在化学非平衡气体模型下,顶板前缘钝化半径(R1)>2mm时进气道的顶板附近和分离区内离解反应较为明显,唇口前缘钝化半径(R2)>2mm时进气道的唇罩、唇口板附近离解反应较为明显;当钝化半径≥4mm时,两种气体模型下进气道出口总压恢复系数和静温的相对变化量绝对值大于0.5%,有必要考虑化学非平衡气体效应对进气道出口性能的影响。     

基于边界层转捩的高超声速进气道特性研究

为了探索边界层非强迫转捩对进气道性能的影响,采用数值计算的方法开展了边界层转捩对轴对称混压式高超声速进气道流场特性的研究。研究表明:随着进气道中心锥锥尖钝化半径增大,边界层转捩先推迟。当锥尖钝度大到一定程度时,边界层转捩位置前移。随着钝化半径进一步增大,边界层转捩再次推迟,转捩位置逐渐后移。来流湍流度越大,边界层越不稳定,边界层转捩越易发生。与湍流边界层相比,考虑边界层转捩时进气道的总压恢复系数及流量系数较高、热载荷及阻力系数较小,Ma=6.5时喉道处总压恢复系数最高上升17.3%,进气道阻力最大下降17.4%。边界层转捩对壁面热流密度分布影响较大,但对壁面压力分布影响较小。钝化影响进气道的自起动性能,随着钝化半径增大,自起动马赫数升高,而边界层转捩对进气道自起动性能影响较小。     

壁面温度控制对平板边界层影响的数值研究

通过对零压力梯度的平板边界层流动施加温度控制,展开壁面温度控制对平板层流边界层和湍流边界层影响的研究,探索温度控制对平板转捩雷诺数和壁面摩擦阻力的影响规律。采用带有转捩模式的三方程湍流模型对平板边界层流动进行数值模拟,重点考察了壁面摩阻系数、平板转捩雷诺数、湍流边界层流动随壁面温度变化的规律。计算结果表明在壁面温度从288 K 增大到432 K 时,边界层转捩雷诺数增大约36%,表面摩擦阻力减少约9.6%。研究分析表明：加热控制使层流区域温度边界层内粘性作用增强,雷诺切应力和湍动能减小,流动更加稳定;而湍流区域边界层内粘性底层中速度梯度和粘性切应力减小,导致壁面处摩擦切应力减小。因此壁面加热控制可以延迟边界层转捩,减小湍流区摩阻系数,并减小平板摩擦阻力。     

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

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

高超声速飞行器-进气道一体化热流数值计算

采用CFD(计算流体动力学)技术, 开展了飞行器前体/发动机一体化气动热环境分析.对层流区、转捩区和湍流区分别采用计算模型, 在湍流区利用压缩性修正的SSGZ-Jk-ε湍流模型, 在转捩区引入代数型转捩因子模型描述边界层由层流逐渐过渡为完全湍流的流动过程.计算了前体和内通道的表面热流, 并与实验结果进行了对比.结果表明所采用的计算方法可以较好地预测前体及发动机内通道热流率, 流动状态、几何结构及激波入射对热流值影响较大.      

粗糙壁面诱导的流动转捩数值模拟方法

为预测粗糙壁面诱导的流动转捩,对Langtry提出的γ-Re_θ转捩模型增加粗糙壁面的比耗散率和涡黏性的边界条件,并对模型中的经验关联函数——转捩动量厚度雷诺数进行修正,引入等效沙粒表面粗糙度作为变量,使模型具有一定的预测粗糙壁面诱导的流动转捩的能力.对粗糙平板自然转捩算例和变压力梯度平板绕流算例进行数值模拟,计算结果与风洞实验数据符合的较好.主要结论如下:表面粗糙度一般会增加边界层内的传热系数和阻力系数,同时使层流到湍流的转捩位置提前;自然转捩的转捩位置受表面粗糙度影响较大,与光滑壁面相比,平板算例中0.15mm的表面粗糙度使转捩位置提前40%;分离诱导转捩的转捩位置受表面粗糙度影响较弱,随着表面粗糙度逐渐增加,转捩位置和分离泡位置略有后移,且分离泡强度逐渐减弱,分离泡之后的阻力系数增加.      

超燃进气道激波/湍流边界层干扰

针对超燃进气道湍流边界层/激波干扰引起的分离问题,采用基于5阶WENO数值格式的大涡模拟(LES)方法开展流场湍流非定常预测,旨在分析进气道湍流化技术实现进气道起动的可行性。研究表明,平板激波/湍流边界层干扰(STBLI)问题,LES方法能够清晰、可靠预测反射、分离激波形成过程及激波与充分发展湍流边界层的相互干扰,定量结果与试验一致;进气道研究方面,层流状态下,激波干扰产生强分离,导致进气道堵塞,而采用湍流化控制后试验和计算均表明流场分离明显减小,流场稳定且无明显堵塞现象,进气道可以起动,总压恢复系数达到要求,该结果表明,利用强湍流化减弱分离,实现进气道起动思想是可行的。      

TSTO并联分离激波/边界层干扰流动特性分析

针对两级入轨飞行器的缩比模型，通过试验与数值模拟相结合的方式，在马赫数6条件下开展典型级间距状态的激波/边界层干扰流场研究，详细分析干扰区壁面及空间的流动结构与特性。结果表明：试验中模型壁面边界层在激波入射之前为层流状态，在强激波干扰后迅速转捩为湍流状态，因此试验结果在第一道激波作用结束之前与层流计算结果吻合，而在第一道激波作用结束之后与湍流计算结果一致；激波/边界层干扰呈现复杂的三维流动特征和明显的开放结构，强激波在壁面形成的高压区呈弧状向下游展开，轨道级头部产生的入射激波在级间来回反射，强度依次递减；同时，干扰区内存在展向弯曲的主分离线与再附线、沿流向排列的二次分离线与再附线、流动剪切形成的旋涡结构以及包括鞍点、结点、焦点在内的临界点；层流边界层受到激波作用形成的分离区明显大于湍流边界层，同时开放特征更为显著。     

真实气体效应对双楔绕流流动特性的影响

为研究高马赫数下真实气体效应对双楔绕流流动特性的影响，采用量热完全气体模型和化学非平衡气体模型两种气体模型对层流流态、转捩流态、全湍流流态的双楔绕流流场进行了数值模拟，其中转捩流态计算采用基于间歇因子的γ-Re_θt转捩模型，全湍流计算中采用k-ωSST湍流模型。结果表明：双楔拐角附近发生流动分离，与量热完全气体模型相比，相同流态下，化学非平衡气体模型计算得到的分离区尺寸更小。在边界层转捩位置的分析中发现，双楔拐角附近的分离区对边界层转捩有重要影响，分离激波产生的位置与边界层开始转捩的位置高度吻合，与量热完全气体模型相比，化学非平衡气体模型时的边界层转捩起始位置更靠后。针对真实气体效应下单位雷诺数、拐角角度对边界层转捩的影响规律研究发现，真实气体效应在一定程度上抑制了分离区内的转捩流动，随着单位雷诺数减小、拐角角度增大，真实气体效应对边界层转捩的影响随之增大。     

转捩模式与转捩准则预测高超声速边界层流动

对原始的k-ω-γ转捩模式和"层流+转捩准则"模型进行了改进,在2种方法中分别增加了横流模态时间尺度和横流转捩准则用于预测横流失稳诱导转捩。通过对网格预处理可并行计算获得边界层外缘信息以及边界层内横流速度。采用不同雷诺数条件下的0°攻角尖锥以及HIFiRE-5外形对2种方法预测高超声速边界层转捩的性能进行了比较分析。研究结果表明,2种方法均能正确反映高超声速边界层转捩起始位置和转捩区长度随雷诺数的变化趋势,但不能捕捉转捩区热流峰值;"层流+转捩准则"模型计算得到的传热系数在全湍流区较k-ω-γ转捩模式偏高。对于同时存在流向不稳定和横流不稳定的HIFiRE-5外形,改进的k-ω-γ转捩模式和改进的"层流+转捩准则"模型相比于原始的模型均能更加准确地预测中心线两侧横流失稳诱导形成的转捩;对于中心线附近因速度剖面拐点引起的边界层转捩,"层流+转捩准则"模型由于与边界层厚度相关,预测得到的转捩位置较试验结果靠前,k-ω-γ转捩模式与试验结果吻合很好。     

压气机叶栅叶片表面附面层流态变化影响因素探讨

以平面叶栅中的二元叶栅模型为试验对象,测量了在不同来流条件下叶片表面流场分布情况及栅后气流参数,分析了不同来流条件下叶片表面附面层流动状态的变化。并借助数值模拟手段重点研究了在不同来流马赫数和冲角下,叶片表面压力梯度对层流附面层向紊流附面层转捩过程的影响,通过利用实验数据分析研究来流条件对转捩过程的影响,为从机理上更深刻地认识叶片表面粘性附面层转捩机制提供了科学参考依据。     

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

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

一种CARET进气道喉部附面层抽吸数值研究

以一种CARET(后掠双斜切双压缩面)进气道为研究对象,设计喉道附面层抽吸槽以控制流动分离。采用CFD数值计算软件对进气道在设计点工况下（马赫数为2.0)下内、外流场进行计算,以总压恢复系数和进气道出口总压畸变为评价指标分析不同抽吸方案的效果。结果表明:喉道附面层抽吸能够稳定结尾正激波,削弱激波/附面层干扰,抑制流动分离,显著改善流场,提高总压恢复系数,减小出口畸变;喉道段抽吸槽位置靠前能够明显降低出口畸变;随着抽吸量的增大,附面层抽吸对进气道内特性性能提升的贡献越来越小。      

来流条件对叶片表面附面层转捩影响的数值模拟研究

通过对二维压气机叶栅叶片表面附面层流态的数值模拟分析,研究了在不同来流马赫数和冲角下,叶片表面层流附面层向紊流附面层的转捩过程,转捩区域的大小等。归纳分析了部分来流条件对转捩过程的影响机制。结果表明:来流马赫数和冲角对转捩起始点和结束点的位置,即转捩所占区域的大小有重要影响。研究工作取得的进展不仅为进一步从机理上深刻认识来流条件对叶片表面附面层流态的影响机理提供了一定的科学依据,而且对今后指导叶片表面附面层流态模似实验工作的开展具有重要参考价值。      

喷注壁面楔状流层流边界层速度与温度相似解

建立了壁面有喷注的楔状流层流边界层冷却数学模型,求解经相似变换得到的描述无量纲流函数和无量纲温度的常微分方程,获得了楔状流层流边界层无量纲速度和温度的相似解,给出了考虑壁面喷注边界条件的形式简洁、物理意义明确的无量纲流函数拟合式;用Runge-Kutta法求解无量纲速度与温度的常微分方程,获得了壁面有喷注的楔状流层流边界层无量纲速度和温度随楔形角、吹风参数、冷却介质温度的变化规律.通过计算0°,18°和36°楔形角的楔状流层流边界冷却速度与温度分布得到以下结论:楔状流速度边界层和温度边界层厚度都随楔形角的增大而变小,随着吹风参数的增大而增大;当冷却介质温度越低、吹风参数越小、楔形角越大时,靠近壁面处楔状流层流边界层内温度梯度越大.      

温度梯度对平板边界层转捩影响的数值模拟

为了研究温度梯度对边界层转捩的影响,在对有/无压力梯度的T3系列平板转捩实验进行数值模拟获得满意结果的基础上,对原型实验进行重新设计,以在平板边界层内形成不同温度梯度并进行数值模拟。计算结果表明随着温度梯度的增大,转捩位置向下游推延发生,平板相同位置边界层的形状参数和壁面切应力系数也相应大幅增高,边界层湍流脉动得到抑制,这是因为温度梯度形成的密度分层,湍流能量为了克服因为密度分层形成的法向浮力而被逐渐耗散,但是温度梯度对转捩影响与其它流动特性对转捩的影响相比较弱。     

Experimental investigation on aero-heating of rudder shaft within laminar/turbulent hypersonic boundary layers

The aero-heating of the rudder shaft region of a hypersonic vehicle is very harsh, as the peak heat flux in this region can be even higher than that at the stagnation point. Therefore, studying the aero-heating of the rudder shaft is of great significance for designing the thermal protection system of the hypersonic vehicle. In the wind tunnel test of the aero-heating effect, we find that with the increase of the angle of attack of the lifting body model, the increasement of the heat flux of the rudder shaft is larger under laminar flow conditions than that under turbulent flow conditions. To understand this, we design a wind tunnel experiment to study the effect of laminar/turbulent hypersonic boundary layers on the heat flux of the rudder shaft under the same wind tunnel freestream conditions. The experiment is carried out in the ?2 m shock tunnel(FD-14 A) affiliated to the China Aerodynamics Research and Development Center(CARDC). The laminar boundary layer on the model is triggered to a turbulent one by using vortex generators, which are 2 mm-high diamonds. The aero-heating of the rudder shaft(with the rudder) and the protuberance(without the rudder) are studied in both hypersonic laminar and turbulent boundary layers under the same freestream condition. The nominal Mach numbers are 10 and 12, and the unit Reynolds numbers are2.4 × 10~6 m~(-1) and 2.1 × 10~6 m-1. The angle of attack of the model is 20°, and the deflection angle of the rudder and the protuberance is 10°. The heat flux on the model surface is measured by thin film heat flux sensors, and the heat flux distribution along the center line of the lifting body model suggests that forced transition is achieved in the upstream of the rudder. The test results of the rudder shaft and the protuberance show that the heat flux of the rudder shaft is lower in the turbulent flow than that in the laminar flow, but the heat flux of the protuberance is the other way around,i.e., lower in the laminar flow than in the turbulent flow. The wind tunnel test results is also validated by numerical simulations. Our analysis suggests that this phenomenon is due to the difference of boundary layer velocities caused by different thickness of boundary layer between laminar and turbulent flows, as well as the restricted flow within the rudder gap. When the turbulent boundary layer is more than three times thicker than that of the laminar boundary layer, the heat flux of the rudder shaft under the laminar flow condition is higher than that under the turbulent flow condition. Discovery of this phenomenon has great importance for guiding the design of the thermal protection system for the rudder shaft of hypersonic vehicles.     

二维进气道不启动流场非定常特性的混合LES/RANS模拟

采用一种混合大涡模拟/雷诺平均Navier-Stokes(LES/RANS)方程模拟方法结合5阶WENO(weighted essentially non-oscillatory)格式对马赫数为3的来流中、内收缩比为1.5的不启动状态下的二维进气道进行了计算,再现了不启动进气道中的非定常流场.计算结果表明:所采用的模拟方法对入口处的平均绝热壁温、摩擦速度和雷诺应力的计算精度较好,进气道不启动流场中激波波系和分离区存在大时空尺度的低频运动,其占主导的特征频率和典型的激波/湍流边界层干扰问题中激波和分离区的低频频率接近,且进气道出现了间歇性的启动状态.