二维喷管的初始流动

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度混合WENO/TCD格式,对Ma=1.4的超声速平面射流初始流场进行了数值研究.数值结果清晰地描述了超声速平面射流初始流场的结构特征,包括主涡环与激波结构以及它们演变过程.因主涡环内存在涡导激波对,激波与涡相互作用加速射流剪切层失稳,使剪切层首次卷起形成小涡的位置出现在涡导激波对后,此与亚声速射流情况不同.小涡串卷起成后,相继与涡导激波对相互作用,使激波出现明显的变形并加速主涡环失稳.      

超声速来流与侧向射流作用下的三维流场结构

研究超声速流动与侧向射流作用下的三维流场结构对飞行器表面流体分离控制以及超燃发动机燃料混合等具有重要作用。基于大涡模拟(LES)方法和高阶TCD/WENO混合格式,对来流马赫数Ma=4,平板上侧向射流的流场结构特性进行了数值模拟。模拟结果表明,流场内的主要波系由弓形激波、分离诱导的激波以及桶形激波组成,这与相关实验结果相符。同时,计算结果还清晰显示了射流与主流相互作用下,射流上游的分离区流动呈现高度瞬态特性,分离区内旋涡结构不断变化,数目由2个到4个再到6个不断发展,平板表面流线复杂。射流尾涡摆动,涡管上扬,呈现非定常特性。     

超声速高度欠膨胀冲击射流的大涡模拟

冲击射流广泛应用于短距起飞垂直降落飞行器(SVTOL)等航空航天领域.本文采用大涡模拟方法对高度欠膨胀的超声速冲击射流的流场进行了数值模拟.本文数值模拟得到了高度欠膨胀冲击射流流场中的激波结构和内外剪切层中不同尺度的涡结构.数值结果观察到了马赫盘的振荡,以及在斜激波、马赫盘及大尺度涡结构的共同作用下,射流内外剪切层之间的环形激波的生成与消失的周期过程.并对流场内剪切层的涡结构的演化进行了研究,数值结果显示内剪切层的大尺度涡结构的形成与马赫盘的振荡相关,在内外剪切层的作用下形成了壁射流区内外交错的涡结构.     

超声速弱欠膨胀冲击射流流场结构

为了对超声速弱欠膨胀冲击射流的流场结构细节进行研究,使用大涡模拟方法对其进行了数值模拟。利用三阶迎风和四阶对称紧致格式对无量纲化轴对称可压缩滤波N-S方程进行空间离散,时间上推进采用的是三阶精度的TVD型Rugge-kutta法。亚格子尺度模型采用的是修正Sm agorinsky涡粘性模型。通过与经典的冲击射流实验比较,证明了程序的可靠性。数值模拟得到了剪切层以及壁面射流中的涡结构和主射流中的激波结构,并且在此基础上对涡合并和板前激波和涡干扰现象进行了深入研究。发现涡合并现象主要出现在流场的上游,越往下游出现的几率越小;涡和板前激波的相互作用会引起激波位置和强度以及冲击平板上冲击区的压强的显著变化,同时也会导致涡的变形。     

射流涡发生器对激波边界层作用诱导的流体分离控制大涡模拟研究

为了研究射流涡发生器对激波边界层作用所诱导的流动分离控制机理及其流场特性,基于大涡模拟(Large eddy simulation)方法和高阶TCD/WENO混合格式,对来流马赫数Ma=2.5情况下,平板上射流涡发生器对激波与边界层相互作用所诱导流场进行了数值模拟。结果表明,射流涡发生器对激波边界层的流体分离有一定的抑制作用,与无控制情况相比,射流作用下进出口总压恢复系数由85.9%提高到94.6%。射流尾涡主要集中于一环状区域内,在该区域内,入射激波与马蹄涡、桶形激波上方的涡管以及剪切涡相互作用,导致整体尾流被激波往下压缩。同时在激波的压缩下,各涡之间相互缠绕、挤压合并,形成多个流向小涡结构,将边界层外的高速流体卷入边界层内,从而增加边界层底层能量,达到抑制流动分离的目的。     

强欠膨胀圆管射流中二次涡环的形成机理

基于大涡模拟方法与高精度混合格式,对出口压力比(n)为1.4与4.0的弱与强欠膨胀圆管射流初始流场进行了数值模拟。结果清晰地描述了两种压力比条件下欠膨胀射流初始流场的结构特征,包括涡与激波结构及它们演变过程,发现强欠膨胀射流中二次涡环是由三波点向下游形成的滑移层失稳而直接卷起,此与弱欠膨胀情况不同。另外,还讨论了二次涡环在绕主涡环外缘向后翻转过程中的演变过程及其对主涡核的影响。当n=4.0时,得到射流域内桶形激波长度与马赫盘直径分别为24mm与0.84mm,与经验公式的计算结果相符。     

亚声速等膨胀平面射流的初始流场结构

基于大涡模拟方法与5阶精度混合TCD/WENO格式,对等膨胀射流在马赫数为Ma=0.6条件下的初始流场结构进行了数值模拟。数值结果描述了射流等膨胀过程中,射流初始流场的形成、变化及发展过程。并发现了射流初始流场中涡核为射流失稳的起始点,研究了射流混合层与涡核剪切层上涡的产生、发展与破碎以及它们卷吸环境气体并与射流气体混合的整个过程,揭示了射流初始过程中涡核与涡环直径的变化趋势,得到了射流轴向速度与压力的分布。数值结果可为相关发动机的设计提供重要参考。     

超声速涡轮叶栅超声速气膜冷却数值研究

为了研究超声速涡轮叶栅通道内的超声速气膜冷却,采用数值计算的方法,对主流压比2.33~4、冷气入射角度15°~45°条件下的涡轮叶栅超声速气膜流动和传热进行了研究。计算结果表明:超声速气膜射流与主流作用后产生的斜激波与尾缘激波交汇,形成两道反射激波,其中一道反射激波作用在气膜孔下游的叶片表面又形成了反射;在不同的主流压力下,超声速气膜射流在叶片法向和展向上展现出不同的发展特征,对转涡对(CVP)在展向上相互挤压,扼制了高温主流卷入叶片壁面;主流压比增加到4,气膜射流区在法向拉长,在展向相对较弱,导致主流在对转涡对(CVP)的作用下被卷入气膜射流的底层,壁面冷却效率降低;气膜入射角从15°增大到45°,冷却效率整体上呈先上升后下降趋势,在入射角30°时冷却效率相对最大,这与射流的穿透能力、冷却气流再覆壁面特征有关。     

基于自适应网格的超声速液体射流破碎过程研究

为研究超声速气流中液体横向射流的破碎过程,采用脉冲背景光方法和VOF方法开展了实验和数值研究。为提高液体横向射流中气液界面和气流场特征捕捉的精确性,采用自适应网格技术对于气液界面、激波出现位置进行网格细化,计算得到了较为精细的气液界面、激波特征及涡系结构。研究结果表明:在低成本仿真模拟条件下,利用自适应网格计算得到的射流轨迹和轮廓与实验吻合较好,射流轨迹的最大误差为10%;射流初始段在超声速气流条件下,仍然存在一段高度约为1.9倍喷孔直径且圆柱形态保持较好的连续光滑液柱。随着喷注压降的升高,液柱的长度逐渐增大;主流气体流经液柱发生三维绕流,在射流附近和近壁面区域形成不断演化的反转涡对,反转涡对的形成加速了液体射流一次破碎过程。     

超声速燃烧气动斜坡喷注器研究

为增强超声速燃烧过程的燃料混合,设计了通过喷嘴阵列形成流向涡强化混合的气动斜坡实验件.在直连式实验台上对气动斜坡实验件进行了纹影及油流谱等冷试实验研究.并采用大涡模拟(LES)方法数值模拟了超声速流场中气动斜坡的流动及燃烧特性.结果表明:气动斜坡喷注器的喷嘴阵列设计能有效加速喷流后流向涡的卷起,提高掺混特性;其燃烧过程也能使喷流形成的剪切层上抬,增进混合;由于喷流方向与主流方向夹角较小,总压损失也较物理斜坡、横向喷流等混合增强方式小.      

亚声速等膨胀方管射流轴置换现象的数值研究

基于可压缩Navier-Stokes方程,采用大涡模拟方法与高精度WENO/TCD混合格式,对Ma=0.6的亚声速等膨胀方管射流的初始流场进行了数值研究。数值结果清晰地描述了初始主涡环的形成与三维演变过程。方管射流具有平面管与圆管射流的一般流场特征,但方管周向曲率的不一致导致轴向流动速度不均匀,使方形主涡环出现复杂的Biot-Savart自诱导变形,并在方形涡环截面上诱导形成4对反向流向涡对。在主涡环的Biot-Savart自诱导变形与流向涡对的诱导速度共同作用下致使主涡环截面形状相对于初始时旋转45°,完成轴置换。另外,在亚声速等膨胀射流中,流向涡对的诱导速度在轴置换中占主要作用。     

超声速冲击射流的PIV实验研究

冲击射流广泛应用于短距、垂直起降飞行器（SVTOL）等航空航天领域，然而却伴随着流场与噪声等诸多方面的问题，笔者采用PIV（粒子图像测速）技术对超声速冲击射流的流场结构和涡结构进行了深入研究。实验发现在不同的冲击工况下，冲击射流流场结构呈轴对称和螺旋结构，其瞬时流场的主涡结构也有类似特征，脉动流场在瞬时流场主涡结构的基础上会附加与主涡旋转方向相反的次涡结构，在对称和螺旋两种模态下，由于涡结构的影响，冲击射流的近壁速度存在较强的脉动。     

高雷诺数湍流横侧射流的大涡模拟

采用大涡模拟(LES)方法研究横向主流与壁面射流均为高雷诺数的壁面横侧射流(JICF),采用动态Smagorinsky涡粘模型对亚格子尺度进行封闭。数值模拟结果中,近场和远场的平均速度、均方根值均与实验结果符合很好。瞬时压力等值面在射流迎风侧与射流迹线成正交关系,表明JICF近场迎风涡是由于射流出口上游剪切层Kelvin Helmholtz不稳定性引起的。基于长度尺度RD,Broadwell和Breidenthal提出了射流迹线公式,通过拟合数值模拟结果得到公式中常数A=1.35,B=0.3。采用平均流线、平均场涡量与Q准则分析三维涡旋结构。横向主流遇到射流后发生夹带现象,在射流出口下游受涡旋结构的影响发生卷起;平均涡量等值面表明,反向旋转涡对(CVP)的涡旋方向与射流出口下游马蹄涡两个分支的涡旋方向相反;由于横向主流和射流雷诺数较高,Q准则表征的涡旋结构在射流出口下游一倍射流直径位置开始发生破碎。     

低雷诺数下横流-射流中剪切涡的试验研究

为深入分析横流-射流(JICF)的流动特性及其中的复杂涡系结构,从流动机理上研究燃机叶片气膜冷却,揭示高温燃气流与冷却流的掺混机理,本文对横向流中单孔射流所形成的剪切涡开展了试验研究。主要研究了速度比、雷诺数及射流角对JICF所形成剪切涡的影响。结果表明:速度比、雷诺数以及射流角会改变主流与射流之间的掺混程度,从而改变射流轨迹的曲率、高度及垂向渗透能力,最终改变剪切涡的特性;迎风涡与背风涡分别是由射流边界层涡与主流边界层涡形成的,当主流边界层涡强度大于射流边界层涡时,背风涡是流场的主导结构,反之,迎风涡将成为流场的主要涡系结构。     

火花放电合成射流与超声速来流相互干扰特性数值模拟研究

为了指导火花放电式合成射流激励器在超声速流动控制中的应用,数值模拟研究了火花放电合成射流与超声速来流的相互干扰特性。研究表明火花放电式合成射流在超声速流场中产生强烈扰动,产生较强的激波结构;随着射流的喷出,激励器上游分离区和流场中激波呈先增强后减弱的趋势,激波由弓形激波逐渐弱化为斜激波,并且随着放电能量的增加射流与主流的动量通量比不断增大,射流的干扰和控制能力显著增强。由于超声速流的较大惯性及其对腔内气体的引射作用,激励器的腔体回填速率大幅下降、回填时间明显增长,使得激励器的工作频率受到很大限制。     

脉冲横向射流的大涡模拟

为了研究脉冲气膜冷却的流动特性,采用基于高阶对称加权本质无振荡(WENO)-Roe格式和隐式牛顿迭代时间推进、动态亚格子模型和预处理技术下的大涡模拟(LES)计算方法,对吹气比为1.0下的斜管横向射流进行了稳态和脉冲两种不同工况下的数值模拟.计算结果显示了射流流场中存在发卡涡结构,研究表明射流孔进口位置的脉冲改变了射流出口拟序结构,即改变了射流与主流的掺混过程.      

钝头双锥喷流致冷流场结构及密度脉动特性

采用大涡模拟方法对钝头双锥喷流致冷流场开展了数值模拟,研究了超声速喷流混合流场结构特征及密度脉动特性。大涡模拟方法基于隐式亚格子模型,空间离散采用高精度通量限制型紧致格式,时间推进采用显式Runger-Kutta方法。数值模拟清晰地捕捉到了流场波系结构,精细地预测了流动发生失稳、转捩以及发展为充分发展湍流的物理过程,直接获得了流场密度脉动特性。通过有、无喷流状态对称面流场的对比,发现超声速喷流能够有效冷却光学窗口;喷流与主流形成的混合层不稳定,很快发生失稳和转捩,形成大尺度湍流结构,进而引起强烈的密度脉动。此外,获得了钝头双锥整体模型喷流致冷流场的空间发展形态特征。      

三维超声速不等密度平面混合层的小激波形成与演变

基于大涡模拟方法与高精度混合格式,对两种密度比(1.0与0.138)条件下的高对流马赫数(Mac=1.1)三维空间发展超声速平面混合层进行了数值模拟。数值结果清晰描述了扰动激励下超声速平面混合层的演变过程,讨论了三维混合层中小激波的结构及其形成过程。分析了混合层两侧密度变化对混合层失稳及小激波形成的影响,发现当上、下层密度不等时,混合层失稳加速,小激波强度减弱,而密度较低侧的小激波会更弱。     

小孔辅助射流提高冷效的机理分析

本文采用数值计算手段,模拟普通的圆柱形单孔结构和小孔辅助射流结构在不同吹风比下的流动和换热,侧重通过直观演示气膜孔下游反向对旋涡对（肾形涡）的生成、发展以及相互作用过程,揭示小孔辅助射流改善气膜冷却效果的机理。结果表明：与单孔结构相比,小孔辅助射流结构,由于小孔射流的干涉作用,主孔射流形成的肾形涡的尺度和强度均有较大程度的减小,冷气射流与主流的掺混减弱,对冷气的向上抬升作用减小,避免了冷气穿透主流脱离壁面,大大提高冷效。随着吹风比的增加,与圆柱形单孔相比,气膜冷却效果改善更加明显。     

脉冲能量沉积对超声速射流/激波相互作用掺混的控制研究

为探究超声速射流掺混增强的有效控制方法,基于二维非稳态雷诺平均Navier-Stokes方程,研究了脉冲能量沉积对超声速射流与斜激波相互作用后增强掺混的控制机理和规律。宽度为5 mm的低密度平面射流同向完全膨胀射入马赫数为2.5的主流中,并与20°压缩斜面所产生的斜激波相互作用。在流场中引入脉冲能量沉积对射流掺混进行控制,并考虑射流马赫数、射流长度、能量沉积位置及激励频率对掺混效果的影响。对比激励流场与基础流场,结果表明:脉冲能量沉积的加入诱导射流形成大尺度涡结构,显著提高射流的掺混效果;能量沉积位置对于不同马赫数射流的掺混增强具有不同的控制效果;激励前后射流宽度的对比表明,脉冲能量沉积的无量纲激励频率在0.22附近时,可使得射流的掺混效果最佳。