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

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

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

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

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

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

超声速横向射流三维流场结构特征

为了研究超声速燃烧室内燃料与空气快速掺混过程的流场特性，基于可压缩Navier-Stokes方程，采用大涡模拟方法和高精度WENO-TCD混合格式对来流马赫数为2.68，喷压比为36的超声速横向射流流场结构进行数值研究。数值结果清晰描述了超声速主流与横向射流相互作用过程的流场结构特征，得到了三维激波形态的演变规律以及它们在强化混合过程中的作用。另外，因桶形激波背风面低压区处的斜压效应，射流气体在桶形激波背风面形成一对螺旋向上的反向涡对，反向涡对的卷吸作用诱导进入壁面边界层的主流向上运动，形成冲击射流。冲击射流以v=557m/s的法向速度向上冲击桶形激波背风面，因而在桶形激波背风面留下类三角锥面凹痕。     

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

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

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

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

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

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

波瓣混合器涡系结构及射流掺混机理的数值研究

借助流体力学软件ANSYS CFX,对波瓣混合器射流掺混流场进行了全三维定常数值模拟,研究了流场中各涡系结构的形成机理及发展过程,并详细探讨了其加速射流掺混过程的作用机制.结果表明:基于SST(shear stress transport)模型的封闭N-S方程能较好地模拟波瓣混合器射流掺混过程,波瓣特殊几何外形诱导产生的流向涡主要通过扭曲内外涵交界面的间接方式加速射流掺混过程,波瓣下游剪切层中K-H(Kelvin-Helmholtz)不稳定性发展而成的正交涡是直接加速射流掺混的关键因素,波谷附近二次流之间的相互作用所产生的通道涡对该区域内的射流掺混有明显的加速作用,受波瓣前缘切割的边界层在径向压力梯度作用下沿波瓣表面卷起而形成的马蹄涡对射流掺混的影响不是特别明显.      

马赫数1.95超声速欠膨胀喷流声辐射数值分析

针对马赫数为1.95的欠膨胀超声速喷流声辐射特性,采用高精度计算格式的大涡模拟（LES）方法进行数值研究。通过对喷流平均流及湍流脉动统计结果的对比,确定数值方法的精确性。细致分析了超声速喷流流动特征,特别关注了其中激波胞格和湍流相互作用现象。基于流动与声场辐射间的关联性,分析了超声速喷流马赫波、宽频激波噪声辐射特性及形成机理。数值结果表明超声速喷流的剪切层演化及激波胞格与湍流间的相互作用构成欠膨胀超声速喷流的主要噪声源。     

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

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

横向喷流与超音速主流干扰流场的数值研究

 采用欧拉方程与可压缩涡对模型相结合的数值方法,研究了超音速主流与横向喷流干扰流场的气动特性。喷流横截面用一对反向旋转的可压缩涡来模拟。采用有限元Osher格式及非结构网格数值求解二维欧拉方程以得到非等熵外流场的气动特性及弓波位置。将欧拉方程数值解结果与可压缩涡对流场解结果相结合以研究喷流的轨迹及喷流横截面的变化。最后给出了数值模拟结果并与工程估算结果进行了比较,表明本文所述方法是可行的。     

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

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

等离子体合成射流对钝头激波的控制与减阻

为了验证等离子体合成射流对超声速流动的流动控制和减阻效果,在考虑热完全气体效应的情况下,工程拟合等离子体热力学属性和输运系数,利用能量源项模型对超声速平板和球头等典型流场结构进行了数值模拟,并取得了与实验较为一致的结果。研究结果表明：对于马赫数为2的超声速平板流动,等离子体合成射流能有效干扰边界层的发展,并诱导产生一系列大尺度结构;对于马赫数为3的球头流场,等离子体合成射流能显著改变激波脱体距离与球头的阻力特性;在放电后第1个周期内,合成射流能使球头平均阻力降低6.3%,而在射流峰值情况下使阻力降低32.0%,激波脱体距离增加2倍,实现了激波控制和流动减阻的预期效果。     

超声速射流离散频率噪声的屏蔽抑制及对流动的影响

根据对超声速射流离散频率噪声产生的反馈机理的认识,提出了一种能够有效地破坏反馈环的形成,从而抑制射流离散频率噪声的喷嘴屏蔽方法.这种方法是通过阻隔反馈声波使其不能到达喷嘴唇口、同时屏蔽物不与射流接触来实现降噪的目的.实验结果表明,对于合适的屏蔽罩参数,正常降噪效果可达5分贝以上.应用激光多普勒测速技术测量了超声速自由射流在采取屏蔽降噪措施前后轴线速度的分布,发现屏蔽后射流势核长度增加近50%,表明这种降噪方法有助于射流推力的提高,但不利于射流的掺混.     

超声速来流稳焰凹腔上游气体燃料横向喷注的流动混合特征

对Ma=1.7的空气超声速来流稳焰凹腔上游气体燃料横向喷注的流动混合特征进行了实验和数值模拟研究。在喷注的模拟燃料氦气中添加丙酮,采用平面激光诱导荧光技术(PLIF)观测了设置开式凹腔时燃料在展向和流向的空间分布。针对下游布置不同凹腔时燃料分布相似的实验结果,选取L/D=7的凹腔对其燃料混合流场进行了大涡模拟。实验和仿真结果表明:实验来流条件下喷注燃料大部分直接由喷流引起的上升反转旋涡对带走,只有小部分燃料经由喷流与凹腔流动涡结构的相互作用进入凹腔剪切层,并随剪切层运动对流进入凹腔内部。     

Simulation of underexpanded supersonic jet flows with chemical reactions

To achieve a detailed understanding of underexpanded supersonic jet structures influenced by afterburning and other flow conditions, the underexpanded turbulent supersonic jet with and without combustions are investigated by computational fluid dynamics(CFD) method.A program based on a total variation diminishing(TVD) methodology capable of predicting complex shocks is created to solve the axisymmetric expanded Navier–Stokes equations containing transport equations of species. The finite-rate ratio model is employed to handle species sources in chemical reactions. CFD solutions indicate that the structure of underexpanded jet is typically influenced by the pressure ratio and afterburning. The shock reflection distance and maximum value of Mach number in the first shock cell increase with pressure ratio. Chemical reactions for the rocket exhaust mostly exist in the mixing layer of supersonic jet flows. This tends to reduce the intensity of shocks existing in the jet, responding to the variation of thermal parameters.     

Subsonic impulsively starting flow at a high angle of attack with shock wave and vortex interaction

Impulsively starting flow, by a sudden attainment of a large angle of attack, has been well studied for incompressible and supersonic flows, but less studied for subsonic flow. Recently, a preliminary numerical study for subsonic starting flow at a high angle of attack displays an advance of stall around a Mach number of 0.5, when compared to other Mach numbers. To see what happens in this special case, we conduct here in this paper a further study for this case, to display and analyze the full flow structures. We find that for a Mach number around 0.5, a local supersonic flow region repeatedly splits and merges, and a pair of left-going and right-going unsteady shock waves are embedded inside the leading edge vortex once it is sufficiently grown up and detached from the leading edge. The flow evolution during the formation of shock waves is displayed in detail. The reason for the formation of these shock waves is explained here using the Laval nozzle flow theory. The existence of this shock pair inside the vortex, for a Mach number only close to 0.5, may help the growing of the trailing edge vortex responsible for the advance of stall observed previously.