高马赫数来流条件下斜激波与流向涡对相互作用

针对斜激波入射流向涡对形成的复杂流动干扰现象,在来流马赫数6条件下,利用后掠斜坡式涡流发生器生成流向涡对,采用数值模拟结合激波风洞实验,研究了斜激波与流向涡对相互作用引起的激波变形和旋涡演变。结果表明,斜激波与流向涡对相互作用的剧烈程度大致分为强、中、弱三种。在强相互作用下,激波面向上游凸起的程度最剧烈,呈现类“T”字型特征,其对称面处接近于一道正激波,波后出现回流区和局部流动滞止,使得流向涡对发生破碎;在中等相互作用下,对称面处存在局部正激波,但没有出现回流区;在弱相互作用下,不再出现局部正激波。进一步地分析发现,斜激波与流向涡对相互作用产生的流向涡对破碎现象,不符合经典的激波作用下单个流向涡破碎理论的预测。采用流向涡对对称面处的气流参数,对经典的流向涡破碎理论进行修正后,能够预测本文流向涡对在斜激波作用下是否发生破碎。本文为掌握激波作用下流向涡对破碎现象的规律,提供了重要参考。     

带舵旋成体侧向喷流流场特性分析

针对带舵旋成体外形,来流马赫数为2.0,迎角为0°,10°和20°,喷口位于×字舵中间,首先开展了喷主流干扰流场的数值模拟与风洞实验校验,得到相符的喷口附近子午线压力分布,验证了数值模拟方法的可行性;然后根据数值模拟得到的喷主流干扰流场结构和流动参数分布,分析了喷主流干扰流场结构随迎角的变化特点;最后从喷口附近旋成体和尾舵表面压力分布的变化分析了喷主流干扰对气动力/力矩放大因子的影响.研究结果表明:随迎角增大,喷流前激波向喷口靠近,分离区减小,喷主流压比增大,喷流穿透高度增大;气动力/力矩放大因子均大于1.0,随迎角增大,放大因子减小.喷口位于×字舵中间时,喷主流干扰对喷流直接力有增益作用;×字舵的存在减弱了喷流的包裹作用,对流场起到了有利干扰.      

超声速飞行器侧向喷流干扰流场传统数值模拟方法的误差分析

使用CFD方法,分别就真实喷管边界和简化喷口边界,计算超声速飞行器侧向喷流干扰流场,研究边界条件对干扰流场及气动力的影响.使用k-ε湍流模型封闭雷诺平均N-S方程,利用非结构网格对流场进行空间离散.通过对比,计算结果与实验值吻合良好,证明该方法具有一定可靠性.进一步研究表明喷流边界条件对喷流干扰流场具有一定影响:相对于简化喷口边界,真实喷管边界喷流出口的非均匀性导致喷口上游分离涡和激波位置较为靠前,从而引起附加气动力和力矩的变化;由于摩擦阻力的作用,真实喷管静推力存在损失;喷流压比为500时,总法向力和总俯仰力矩在两种边界条件之间的误差分别为8.21%和22.4%,误差较大.在进行侧向喷流干扰流场的精确计算时,需要考虑边界条件的影响.      

超声速侧向多喷流干扰特性数值模拟

为了研究多个喷流喷管对导弹控制力的干扰影响,本文通过数值求解N-S方程来模拟超声速外流场中横向喷流的干扰流场,采用分块对接网格和"O"型网格技术,精确模拟喷口截面及弹翼形状,生成高质量的贴体计算网格。通过对多种喷管控制组合的超声速横向喷流干扰流场的数值模拟,研究和分析了喷口附近流场的涡系结构和波系结构,并将喷管几种排列组合对导弹喷流干扰力放大因子的影响进行了分析研究,得出一些多喷流干扰的结论。     

超声速主流中横向喷流干扰流场数值模拟

本文用二阶迎风 TVD格式求解层流 N- S方程 ,数值模拟了超声速主流与横向喷流相互作用复杂干扰流场 ,得到了合理的涡系结构 ,计算波系结构及沿物面压力分布也与实验结果一致。计算还表明 ,TVD格式不但具有较强的捕捉激波的能力 ,而且用来模拟复杂涡系结构也是有效的     

尖拱弹身横向喷流数值模拟

本文通过数值求解层流Navier—Stakes方程来模拟超声速流场中横向超声速喷流的干扰流场，采用分块对接网格和“O”型网格技术，精确模拟喷口截面形状并生成高质量的贴体计算网格，对尖拱弹身的横向超声速喷流的干扰流场进行了数值模拟，研究了喷口附近的流场结构、涡系结构、波系结构及流动分离等流场特性，得到喷流流场的干扰效应，计算结果与实验值吻合较好。     

导弹侧向喷流干扰及多喷口耦合效应数值模拟

侧向喷流与外部来流的干扰流场相当复杂,流场内会出现弓形激波、再附激波和分离旋涡等复杂的物理现象。通过数值求解NS方程,对导弹的侧向喷流干扰流场进行了数值模拟研究,重点讨论了采用空气冷喷流进行喷流干扰模拟的相似模拟准则,通过与燃气喷流的对比,验证了喷流干扰模拟准则在导弹侧向喷流干扰数值模拟中的可靠性。利用建立的侧向喷流模拟方法,对某导弹外形的多喷口耦合效应进行了数值模拟研究,分析了侧向多喷口耦合干扰下的放大因子及流场结构,相关结论可为导弹喷流控制系统设计提供参考依据。     

高压涡轮封严流与主流相互作用的机理研究

为探究高压涡轮轮毂封严流与主流相互作用的机理，借助经过实验校核的数值模拟方法，详细分析了封严流对主流端区二次流结构的影响以及封严流与主流的相互作用过程。研究发现：一方面，主流从叶片前缘位置侵入封严结构内部，在封严出口处形成封严回流涡，并在封严结构内部诱导出一个反向涡，这两个涡直接影响封严结构的封严效率；另一方面，封严出口处封严回流涡与叶片通道内的马蹄涡压力面分支在流向上旋转方向一致，互相融合并增强通道涡强度。封严结构决定了封严回流涡流出的位置和速度方向，直接影响封严回流涡与马蹄涡压力面分支的相互作用过程，从而决定了损失的大小。研究还发现，当封严流和主流在封严出口交界面上流量相当且存在一定的周向速度差时，封严出口会发生Kelvin-Helmholtz不稳定现象。此时伴随大量边界层低能流体进入封严结构内，封严流周向速度减小，马蹄涡的压力面分支和封严回流涡随之减弱，继而使端区二次流损失减小。     

叶片式涡流发生器对压缩拐角流动分离的控制

采用基于k-ω湍流模型的雷诺平均Navier-Stokes（RANS）方程方法,研究了叶片式涡流发生器（VG）对于马赫数Ma_∞=2.9时24°压缩拐角边界层分离的控制作用。计算结果表明：叶片式涡流发生器诱发的流向涡,是控制拐角处边界层分离的主要因素,流向涡强度越大控制效果越好。流向涡增大了主流与边界层内的动量输运,沿壁面法向速度型更加饱满,并使得压缩拐角处的二维分离转变为三维分离,改变了激波边界层干扰的结构,分离区长度减小了39.68%。相比于相向旋转,同向旋转叶片式涡流发生器改善了分离区内的压力分布,分离区总长度减小量相当,但分离点距转折点处的长度更短,且系统阻力增量更小。对于相向旋转叶片式涡流发生器,后缘高度增大,分离区总长度减小,系统阻力增量先减小后增大;相向旋转叶片间距越大,分离区总长度越小,系统阻力增量越大;同向旋转叶片间距越大,分离区总长度越大,系统阻力增量越小。高度对叶片式涡流发生器诱发的流向涡强度起主要作用,异向与同向叶片间距的影响较小。     

稀薄过渡流区横向喷流干扰效应数值模拟研究

为了准确预测稀薄过渡流区横向喷流与稀薄大气的干扰流动特征,本文在建立直角与表面非结构网格混合结构的DSMC数值算法以及碰撞网格自适应算法基础上,提出基于MPI的静态随机负载平衡技术,构建了DSMC并行计算代码。计算的不同压力比条件下三维平板模型横向喷流与稀薄大气干扰的分离长度与低密度风洞试验有较好的一致性,验证了本文算法的可靠性。开展了细长钝双锥外形高超声速稀薄来流与超声速/高超声速横向喷流干扰效应的研究,计算分析了不同飞行高度、不同飞行速度、不同飞行攻角、不同喷流推力下复杂流场结构和对气动力特性的影响规律。考察了RCS喷管出口参数不同(均匀/非均匀)对喷口附近分离涡和分离长度的影响。     

交叉射流特性的试验研究

本文研究了稳燃腔中交叉射流两相流动。使用三维激光粒子动态分析仪 (PDA)测量了正三角形截面副射流与主射流速度比为 2 .5时的交叉射流流动 ,获得了粒子湍流脉动强度及粒子在交叉射流尾迹中的浓度分布     

进气横流对前缘冲击流量系数的影响

对靶面为圆弧面的单孔前缘冲击流量系数进行了实验测量,通过改变进气横流和单孔射流间的速比大小来模拟一排冲击孔中各孔实际所感受横流的不同。研究了孔雷诺数和进气横流速比对流量系数的影响,孔雷诺数范围:1.5×104~4.0×104,横流速比范围:0~0.3。结果显示冲击孔流量系数随孔雷诺数的增大而增大,最大增幅7%左右;随横流速比的增大而线性减小,最大减幅15%左右。由实验数据得到了冲击孔流量系数分别与孔雷诺数和横流速比的经验关系式。     

非均匀横流作用下冲击射流冷却的数值研究

运用数值计算的方法对非均匀横流作用下冲击射流冷却的流动和传热过程进行了三维数值研究,揭示出射流冲击孔和横流孔相对排列方式、射流雷诺数和横流雷诺数相对大小、射流孔中心至横流孔出口截面距离等因素对冲击射流冷却流动传热特性的影响。结果表明:若保持总的冷却流量不变,当横流雷诺数与射流雷诺数之比等于或大于0.5后,有横流时的冲击射流冷却局部努塞尔数比无横流时有较为显著的下降;在非均匀的横流作用下,对于叉排和顺排方式,由于射流和横流的相互作用,壁面温度分布呈现出复杂的特征;射流孔和横流孔呈顺排布置时,峰值努塞尔数所对应的射流驻点区向下游偏移。      

高精度ENO格式在射流数值模拟中的应用

本文分析了ENO格式的特点，并应用于全NS方程的迁移项和压力顶，模拟了轴对称和三维射流。首先对轴对称自由射流为算例比较了高阶TVD格式和ENO格式，证实了ENO格式捕获复杂激波波系的能力；其次计算了四喷管湍射流干扰流场，给出了相应的计算结果；最后计算了四股射流流入圆管内的流场，给出了物理量沿壁面和截面的分布。     

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

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

凹腔前缘角对超声速燃烧室性能的影响

针对带有不同前缘角的凹腔内流动和燃烧过程,分别在冷态和燃烧条件下探讨了前缘角对凹腔内流动损失及阻力特性的影响.研究表明:在壁面垂直喷射的喷口上游和凹腔内部均会形成低速、高温回流区,有利于点火及火焰稳定,燃烧反压通过边界层的亚声速区域上传,形成激波/边界层干扰结构.减小前缘角,可使剪切层分离位置提前,更偏向凹腔内部,导致凹腔后壁面再附激波增强,进而增大了总压损失,降低了总压恢复系数;亦可导致凹腔前、后壁面压差阻力增大,阻力系数上升.进一步认识了凹腔内部流场及稳焰增混机理,进而为优化凹腔结构设计提供依据.      

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

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

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

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

支板喷射超声速湍流燃烧的大涡模拟

为了发展可行的超声速湍流燃烧大涡模拟方法,将设定型PDF(Probability Density Function)模型与LES(Large Eddy Simulation)相结合以封闭亚格子湍流-燃烧相互作用,并将模型用于支板喷射超声速湍流燃烧流场的数值模拟。分别对冷流及燃烧流场进行了模拟,计算结果与实验测量符合较好,表明了所采用方法及模型的可行性。冷流条件下,大尺度湍流涡通过卷吸、拉伸运动主导支板尾迹区的近场混合,并通过破碎过程影响远场混合。燃烧条件下,回流区尺度扩大,剪切层中形成的高温燃烧产物通过大涡卷吸以及回流区末端对流作用进入回流区并与其中的燃料喷流相互作用,使部分燃料预热升温并进入回流区两侧剪切层与主流空气混合、燃烧,从而实现火焰稳定。在薄反应剪切层及大尺度反应涡的边界区域,LES网格不足以直接求解湍流与燃烧的相互作用,PDF模型给出了较强的亚格子脉动。     

中心分级燃烧室头部冷态流场特性实验研究

为了深入研究中心分级燃烧室的流场特性,采用PIV方法对其头部冷态流场开展了实验研究,重点分析了中心分级燃烧室头部流场的结构特征以及主燃级旋流数对头部流场的影响。实验结果表明：中心分级燃烧室头部典型流场结构形成一个较大的中心主回流区、较小的端部回流区及角回流区,主、预燃级气流的相互耦合过程分为独立射流、气流掺混和气流合并三个阶段。主燃级旋流数对头部流场结构影响较大,随主燃级旋流数增大,主、预燃级气流耦合作用减弱,回流区分布形态发生较大变化,新的流场结构特征为中心主回流区范围明显收缩,端部回流区向下游延伸,在主、预燃级之间剪切层形成较长回流区。针对该中心分级燃烧室头部,主、预燃级气流由耦合流动变为解耦流动的临界主燃级旋流数为0.8~1.0。