超声速欠膨胀冲击射流流场振荡的数值模拟

为了研究超声速欠膨胀冲击射流自激振荡与喷嘴进出口压比的关系,对轴对称冲击射流进行了数值模拟,计算中采用了雷诺应力模型。结果发现:当喷嘴压比较小时,流场振荡随喷嘴压比的增大而加剧。当压比处于某一范围内时,流场振幅最大,振荡最剧烈。随着压比的继续增大,流场振幅又减小。流场振幅随喷嘴压比的变化与冲击单音强度随喷嘴进出口压比的变化相似,说明流场振荡与冲击单音有内在的关联。     

三维亚声速冲击射流流场的数值模拟

利用有限体积法对三维可压缩的N－S方程进行离散，对喷嘴的亚声速垂直冲击射流和斜冲击射流进行了数值模拟。网络的划分采用非结构性网络，湍流模型为RNGκ-ε湍流模型，在近壁处采用壁面函数进行修正，流场计算结果与实验值吻合一致。计算可得到垂直冲击流场结构的三个区域：自由射流、滞止区和壁面射流区域，并显示出垂直冲击点附近区域的卷吸回流、速度分布。冲击射流近壁处存在着回流现象，由于边界层流动和近壁回流引起距离固壁2D内的速度梯度变化较大。斜冲击在近壁处平行壁面速度的最大值，以及壁面上的总压最大值都偏离冲击点，并且相对应。     

超声速欠膨胀冲击射流的数值模拟

利用有限体积法对轴对称可压缩的N－S方程进行离散，对喷嘴的超声速垂直冲击射流进行了数值模拟，网格的划分采用非结构性网络，得到详细的波系结构，激波形状和位置与实验相吻俣，大压比下，壁面上冲击中心区域压力相对低，分布较平坦，压力最大点所在位置偏离中心冲击点，并且在近壁处有回流。     

计算流体动力学在航天飞行器射流逆流中的应用

超声速射流逆流通常用于导弹、航天飞机、卫星和飞船等飞行器运动状态的控制。欠膨胀超声速射流逆流的流场包含有多激波（如弓形激波和马赫盘）、接触间断和剪切层，其结构非常复杂。本文采用激波高分辨率有限差分（ＴＶＤ）格式，对恒定自由流条件，各种不同射流出出口压比的超声速轴对称逆射流进行了数值模拟，且对各种条件下的物理现象给予了分析。计算的马赫盘和弓形激波位置与实验值相吻合，为此类流动问题提供了一种有效的预测     

气液交叉射流动力学行为的基础研究

为了模拟燃油在高速喷嘴内的快速蒸发混合过程,基于Eulerian—Eulerian方法,用RNGk一8湍流模型建立了高速喷嘴内两相流流动的数学模型,数值模拟了高速喷嘴内气液两相交叉射流的流动规律。改变气液动量比和韦伯数等条件下的两相流场计算结果表明：动量比增大3倍时横向穿透深度略有增大;而较大的韦伯数能起到强化横向射流的破碎效果。计算结果与同条件下的试验结果相符合。结果为研究HiTAC快速速喷嘴内可燃混合气的快速形成奠定了基础。     

通过矩形喷嘴的非稳定液体射流

介绍了关于用矩形喷嘴进行的瞬态液体射流的实验研究结果。流动显示演示了在喷流过程中液体射流的种种不稳定现象。笔者推广了Ryhming的理论模型，计算了喷嘴出口的射流速度，并与实验结果进行了比较。     

内管射流占优时轴对称共轴射流结构分析

以轴对称Reynolds-averaged Navier-Stokes(RANS)方程为控制方程模拟了内流占优条件下轴对称共轴喷嘴射流流场,并基于湍动能、速度和压强等流场参数的变化规律,将这类流场划分为交汇段、初始段和主流段等三个区段,给出三段流场的结构特征和流动机理,发现初始段内湍流动能大大低于单管射流,将导致共轴射流噪声低于单管射流;主流段内湍流动能则高于单管射流,将加强射流的掺混作用.将流量计算结果与试验结果进行了对比,证明使用的算法是有效的.      

冲击射流流动换热超大涡模拟研究

为了准确预测发动机热端部件中广泛采用的冲击射流冷却复杂的流动和换热特性,发展了基于BSL k-ω模型的超大涡模拟(VLES)高精度模拟方法,并对高雷诺数Re=4×104,两种不同射流距离2和6的单孔冲击射流及三孔冲击射流这一经典的流动传热问题进行三维非稳态高精度数值计算。同时,将分离涡方法 (DDES)和k-ωSST,RNG,Transition SST三种RANS方法的数值模拟和开发的超大涡模拟(VLES)方法进行对比。研究表明,VLES方法均能够准确捕捉冲击射流流场的复杂非稳态流动及传热特征,包括自由射流区、壁面射流区小尺度涡系和大尺度湍流结构的演化和破碎,同时冲击壁面的换热系数计算结果与实验值吻合较好。DDES方法未能准确捕捉流场复杂的小尺度湍流结构,壁面换热计算结果与实验值差异较大。RANS方法计算的换热结果与实验数据差异最大,基本未能预测到壁面换热特性。在相同的计算网格和计算方法下,VLES方法计算结果优于DDES方法,DDES方法一般好于RANS方法。这表明新开发的VLES方法能够准确地计算冲击射流相关的流动及换热问题。     

三维湍流冲击射流流动与传热特性的数值研究

采用SIMPLE算法和RNGk-ε湍流模型,通过求解三维Navier-Stokes方程和能量方程,对矩形管和方管湍流冲击射流的流动及传热特性进行了数值研究。模拟结果显示出在矩形管和方管射流的冲击面附近分别存在主流方向速度的两个和四个偏心峰值,偏心速度峰值的形成是当射流冲击壁面时,沿壁面射流方向的涡量扩散和旋涡拉伸作用而引起的;偏心速度峰值的形成与射流Re数和喷嘴-冲击板间距的大小密切相关,而偏心速度峰值的个数则取决于射流喷嘴的几何形状;冲击射流的传热特性受流场结构的控制,偏心速度峰值的存在导致了冲击壁面上局部Nu数偏心峰值的形成。     

不同拓扑结构喷嘴内交叉射流的数值研究

利用CFD软件对五种不同拓扑三维燃料—空气快速混合喷嘴进行了数值计算;应用两相流模型和RNGk-ε湍流模型,描述了快速混合喷嘴内两相流动过程。计算得出5种拓扑燃料—空气交叉射流喷嘴中湍流强度分布、总压损失系数、在流向方向各个截面温度分布和出口径向温度分布、出口速度分布和出口径向质量分布系数。通过对这些数据的分析,结果表明,射流下游添加肋结构可以有效增强混合气的快速形成,肋添加越多,对燃料的破碎和燃料—空气的混合越有利。     

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

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

复杂几何形状喷管内外三维流场的数值模拟

研究了球面收敛矩形扩张的二元喷管（ＳＣＦＮ）流场的数值模拟,为了有效地捕获激波与滑流面及较好地模拟分离,采用了Ｒｏｅ格式、ＭＵＳＣＬ方法及Ｖｅｎｄａｔａｋｒｉｓｈｎａｎ限制器、Ｓｐａｌａｒｔ－Ａｌｌｍａｒａｓ方程紊流模型。研究中给出了两个喷管模型四个流动状态的数值模拟结果。结果表明,除了有严重分离的流动情况以外,该算法能得到相当满意的计算结果。     

喷流-外流干扰流场数值模拟

 研究了一种适用于推力矢量的内外流干扰复杂流场的数值模拟方法。主要讨论了能灵活处理复杂外形的重叠-搭接网格技术;适于低速大迎角流动及重叠 搭接网格技术的湍流模型以及是否考虑内流的喷流边界条件的定义。在数值模拟方法研究基础上对矢量推力飞机喷流-外流干扰流场进行了分析研究。包括计及喷管内流的内外流干扰流场数值模拟;给定喷管出口边界条件的喷流-外流干扰流场数值模拟;复杂外形引入喷流边界的数值模拟。以上计算反映了喷流对不同外流情况时的流扬的影响和气动力的影响。数值模拟说明,采用重叠-搭接混合网格处理复杂外形,使用恰当的湍流模型可以较好地模拟喷流-外流的干扰流场。     

适于黏性计算的自适应笛卡儿网格生成及其应用

复杂外形的黏性网格生成已经成为计算流体力学(CFD)在工程应用中的主要瓶颈,高效而鲁棒的网格生成技术研究显得尤为迫切。使用基于物面外形自适应加密的笛卡儿网格方法生成计算网格,投影方法拟合壁面边界,将投影得到的柱形单元沿法向分层获得适于黏性计算的网格。应用特征恢复技术恢复凹面特征并改善凸面处的网格质量。流场解算中应用最小二乘法重构获得2阶精度,Roe格式计算无黏通量,Holmes-Connell方法计算黏性通量。加入Spalart-Allmaras一方程湍流模型以获得湍流解。上下对称高斯 赛德尔迭代(LU-SGS)格式隐式时间推进,对网格重排序保证格式的稳定性和平衡性。算例显示网格生成方法是快速的和鲁棒的,流场计算精度满足工程需求,因此该方法非常适合应用于实际工程计算。     

合成射流与横向主流作用的流场数值分析

建立了将合成射流激励器腔体、出口喉道、外部流场作为单连域计算处理的吹/吸型边界模型。在此基础上,对横向主流偏转角及上下壁面压力变化进行了讨论。计算结果表明,激励器布置位置是影响主流偏转角的重要因素,在主流通道内侧,激励器与主流出口之间的上下壁面压力差较大。文中还对两种不同模式下合成射流对主流的控制效果做了比较。     

相邻激励器合成射流流场数值模拟及机理研究

建立了将合成射流激励器腔体、出口喉道、外部流场作为单连域计算处理的吹／吸型边界模型。在此基础上,对不同相位差、不同振幅、不同频率的相邻激励器相互作用形成的合成射流流场进行了数值分析。计算结果表明：相邻激励器工作时的相位差、振幅不同、驱动频率不同对其形成的合成射流流场有很大影响,合成射流不再对称分布,流动将发生偏转。其机理是由于两激励器吸入和排出流体流动不同(不同相、不同幅值、不同频率),使得两列旋涡对不对称,因此在两列旋涡对之间存在涡量强度不同和压强梯度,从而引起旋涡对向低压侧和强涡量区偏转。     

斜切喷管欠膨胀超声速射流的数值模拟

采用三阶MUSCL TVD格式解三维可压缩平均雷诺纳维尔-斯托克斯方程组，湍流模型为Spalart-Allmaras代数模型，数值模拟了两种斜切角喷管在喷口压比为2.0，2．8，3．4和4．0条件下的欠膨胀超声速射流场，获得了射流场流谱和参数分布，计算结果和实验数据符合良好，并给出了喷管斜切角和喷口压强比与射流偏转角和扩张角的分布曲线。     

用双流体模型模拟超声速气固两相流动

采用了双流体模型对JPL(Jet Propulsion Laboratory)喷管中气固两相流动以及超声速射流两相流动进行了数值模拟,并研究了可压缩两相流动中气相与颗粒的相互作用规律.自主开发的一般曲线坐标系下二维轴对称可压缩双流体程序Solve2D,对气相求解Navier-Stokes(N-S)方程组,采用k-ε湍流模型,颗粒相求解Euler方程组.对JPL喷管内的两相流场和湍流两相射流流场进行了数值模拟,研究了不同颗粒质量百分数以及不同颗粒直径时的气固两相流场的流动规律.      

改进S-A湍流模型对横向射流的CFD模拟

通过在Spalart-Allmaras(S-A)湍流模型中加入应变率张量，来改善模型对横向射流(JICF)的计算精度。模拟了速度比分别为0.5和1.5工况下的横向射流。经实验值验证，结果表明:改进后的S-A模型预测精度有所提升；并对湍流模型改进部分对横向射流的作用机制做出了分析。将改进后S-A模型的计算结果与k-ε模型，雷诺应力模型(RSM)和大涡模拟(LES)的结果进行比较，发现S-A模型在速度比为0.5时精度最高，而在速度比为1.5工况下没有明显优势；并对这一现象做出了合理分析。