超声速气流中横向煤油射流的数值模拟

为了对超声速气流中液体横向射流进行数值模拟,采用Eulerian Lagrangian方法描述气液两相流动现象,用随机轨道模型追踪液滴的运动历程,并考虑气流可压缩性、流场不均匀性及液滴变形对液滴运动的影响,开发了相应的计算程序。计算结果发现,液体射流与气流之间存在着强烈的相互作用,液滴在进入气流中不久就破碎成很小的子液滴,受煤油液雾的影响气流速度、温度急剧下降,同时在喷孔上游出现一道弓形激波。与液雾结构的纹影图像和液雾穿透的实验结果对比显示,数值计算结果基本合理。     

颗粒两相流自由剪切层大涡模拟

采用大涡模拟(LES)方法数值仿真了可压缩流自由剪切层内颗粒两相流流场结构.气相控制方程采用二维盒式滤波法滤除小于网格尺度的小涡以得到LES控制方程,而流场中小涡对大涡的影响采用K-ωSST湍流模型模拟.在Lagrange坐标系下用颗粒轨道模型仿真颗粒运动,用双线性插值方法实现Euler坐标系到肌Lagrange坐标系的转换,并考虑了颗粒速度和温度的改变对气相流场的反作用.分析了仿真结果中气相自由剪切层的涡结构特性,不同St数颗粒在流场中的运动特性,颗粒对流场的反作用效果,得到了合理的结论.     

有蒸发稀液雾冷态射流的大涡模拟研究

为研究气相湍流场中液滴弥散和蒸发过程的特点和影响因素,加深对气液耦合作用的认知,用欧拉-拉格朗日相结合的办法数值模拟了稀液雾冷态受限射流和自由射流。欧拉坐标系下的气相湍流场使用大涡模拟技术模拟,离散液相则使用拉格朗日颗粒轨道模型进行描述,考虑相间质量、动量、能量交换,采用相平衡蒸发模型模拟液相的蒸发效应。模拟得到的液相速度场、粒径分布和质量流率等均和实验数据较好吻合。有回流的受限射流中初始惯性对Stokes数大的液滴弥散影响较大,而Stokes数在1附近的小液滴弥散更受气相流场结构影响。液雾单位体积蒸发速率和液滴弥散分布密切相关。自由射流中液滴受初始惯性和受气相流场作用大体一致,大小液滴的弥散分布基本满足自相似性。     

涡流器燃烧室头部两相反应流数值模拟

在任意曲线坐标系中数值模拟带涡流器燃烧室头部的三维两相反应流场,其中,采用偏微分方程法生成块结构网格;气相在Euler坐标系下采用修正的κ-ε双方程模型处理紊流特性,用修正的EDM-Arrhenius紊流燃烧模型模拟化学反应速率,用六通量模型考虑热辐射的影响;液相处理采用在Lagrange坐标系下颗粒群轨道模型.在非交错网格系下采用SIMPLE算法及混分差分格式求解离散方程,采用单元内颗粒源法(PSIC)算法进行两相之间的耦合,计算所得的旋流数与试验数据相当接近.结果表明涡流器几何形状对燃烧室流场影响较大,涡流器的优化设计可以提高燃烧性能.     

航空发动机排气系统红外辐射测量与数值计算比较研究

评估航空发动机红外隐身性能,需获取其排气系统的红外辐射强度。介绍了航空发动机排气系统红外辐射测量方法,通过试验验证得到了3~5mm的红外辐射特征。采用反向蒙特卡洛法对发动机排气系统的红外辐射特征进行数值计算,计算中将固体壁面的反射特性分别设为漫反射体和镜面反射体,并将计算值与试验值进行对比分析。结果表明:两者分布规律一致,但将壁面的反射特性作为镜面反射处理时,计算值与试验值的误差明显减小。     

带壁面射流燃烧室燃烧数值模拟

利用基于非结构化网格有限体积法对三维有壁面射流的燃烧室内两相流动和燃烧进行了数值研究.对气相流动在Euler坐标系下求解,而对液滴相则利用Lagrange方法进行追踪求解.计算区域采用四面体网格进行划分,气相流场用SIMPLEC计算方法,对液滴相采用了欧拉隐式方法.考虑了液滴相与气相的完全双向耦合作用,分别采用了Spalding液滴蒸发模型和涡破碎(EBU)燃烧模型,数值计算结果与文献中实验数据吻合较好.      

超声速气流中液滴破碎的大涡模拟

为了研究超声速气流中的液滴破碎过程,采用基于CLSVOF界面追踪方法的大涡模拟方法。该方法中用可压流和不可压流求解器分别计算超声速气流流场和液体流场,并且在分别计算气相或液相时,另一相作为该相计算时的边界条件参与计算,从而实现可压流求解器和不可压求解的融合。通过去散度化液体速度外推解决界面处大液气密度比所带来的动量误差大的问题。利用该方法对来流Ma=1.358的超声速气流中的水滴破碎进行了数值模拟。模拟结果展示了超声速气流中液滴破碎的几何结构变化细节和基本破碎形态,且与实验吻合较好。数值模拟结果揭示了Rayleigh-Taylor不稳定和气动剪切力对液滴破碎的作用过程。破碎过程中,液滴在垂直来流方向上的最大无量纲宽度(Dmax/D0)约为4.5,无量纲破碎时间(tbU∞/D0)约为29.85。     

超声速混合层液滴两相流动的大涡模拟

为研究超声速气流中液滴与气流的混合及液滴蒸发对混合的影响,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构,气相流场采用亚格子(SGS)模型和切应力输运(k-ωSST)湍流模型,液相模拟采用轨道模型和单液滴蒸发模型。在混合层前缘入口处均匀持续地投放液滴,并在液滴入口处下方添加非周期小扰动,并观察液滴蒸发过程对该小扰动造成的影响。分析了入口小扰动在流场中不同的发展情况,发现液滴的蒸发过程使混合层厚度增加并加速混合层的发展,对气相流场扰动较强,可能导致流动失稳,对混合过程有很大的促进作用。     

带圆筒头部火焰筒两相燃烧流场及性能计算

在三维任意曲线坐标系下采用代数雷诺应力紊流模型模拟紊流粘性,Arrhen inus-EBU紊流模型模拟燃烧室内化学反应速率,采用随机轨道模型模拟气液两相之间的相互作用,数值分析不同进口气流参数对带双圆筒头部环形燃烧室两相紊流燃烧流场的影响,同时采用多维经验分析法预估燃烧室性能。通过两种工况计算表明:在最大工况下所得的燃烧流场和燃烧室性能均优于慢车工况。计算结果与试验数据符合较好,表明燃烧过程的数值模拟以及多维经验分析法可为燃烧室优化设计提供有用的设计依据,适合于工程应用。     

回流环形燃烧室反应流场数值模拟

采用自编计算程序,在任意曲线坐标系下对带扩压器、涡流杯、火焰筒、内外环冷却通道与弯管的回流环形燃烧室三维两相化学反应整体流场进行计算.数值预测火焰筒头部结构参数对燃烧流场的影响.所得计算结果与实验数据符合较好,结果表明:采用的数学模型和计算方法合理,该软件可供回流燃烧室优化设计与研制使用.      

Extension of the low diffusion particle method for near-continuum two-phase flow simulations

The low diffusion (LD) particle method, proposed by Burt and Boyd, is modified for the near-continuum two-phase flow simulations. The LD method has the advantages of easily coupling with the direct simulation Monte Carlo (DSMC) method for multi-scale flow simulations and dramatically reducing the numerical diffusion error and statistical scatter of the equilibrium particle methods. Liquidor solid-phase particles are introduced in the LD method. Their velocity and temperature updating are respectively, calculated from the motion equation and the temperature equation according to the local gas properties. Coupling effects from condensed phase to gas phase are modeled as momentum and energy sources, which are respectively, equal to the negative values of the total momentum and energy increase in liquid or solid phase. The modified method is compared with theoretical results for unsteady flows, and good agreements are obtained to indicate the reliability of the one-way gas-to-particle coupling models. Hybrid LD-DSMC algorithm is implemented and performed for nozzle discharging gas-liquid flow to show the prospect of the LD-DSMC scheme for multi-scale two-phase flow simulations.     

高超声速稀薄流的气粒多相流动DSMC算法建模研究

基于直接模拟Monte Carlo(DSMC)方法,构造适用于DSMC算法的固态和液态颗粒碰撞、聚合和分离模型,发展稀薄条件下双向耦合作用的气粒多相流的DSMC算法,在此基础上初步实现高超声速稀薄流环境中的气粒多相喷流流场数值模拟.算例结果表明该方法能为稀薄过渡区气粒多相流动提供一种新的应用研究手段.     

Investigation on plume interference effect of solid propellant micro-thruster

The three-dimensional numerical simulation of two-phase plume flow of solid propellant micro-thrusters was developed. Then it was used to investigate the plume interference effect by combining the direct simulation Monte Carlo (DSMC) method for multi-component gas flow with the two-way coupling model for two-phase rarefied flow. At different space between the two micro-thrusters and different wall temperature,the plume interference effect was analyzed specifically. The results show that under the plume interference effect the gas is compressed and the flow direction is changed,which resulted in the increasing of gas pressure and temperature; solid phase made no significant effect on the flow parameters of gas phase; with the rising of the space between the two micro-thrusters,the maximum pressure decreased and the maximum temperature increased in the domain under the plume interference effect; the wall temperature could influence the temperature of the gas which is extremely close to the wall,but not the gas pressure.      

三维管道超声速气固两相流动的CFD/DSMC仿真

采用NND格式求解气相N-S方程，采用DSMC方法模拟固相颗粒运动，数值模拟了三维管道超声速气固两相流场，分析了固体颗粒对流场结构产生的影响。计算结果表明，将DSMC方法与CFD方法相结合是解决气固两相流动问题的一种有效途径。     

高超声速平头圆柱绕流化学反应流场的DSMC/EPSM混合算法研究

本文发展了平衡粒子模拟方法(EPSM)，建立了与高温气体化学反应动力学理论相匹配的EPSM耦合模型，并通过混合参数进行流区的自动识别，将EPSM方法与蒙特卡罗直接模拟方法(DSMC)结合，构造了可模拟化学反应流动的DSMC／EPSM混合算法。应用该算法对汲及化学反应的轴对称情况下高超声速平头圆柱绕流流场进行模拟，将结果与DSMC方法的结果进行比较，验证了新算法对求解化学反应流动的可行性。将混合算法的计算效率与DSMC方法的计算效率进行比较，发现混合算法能够大大提高计算效率。     

超声速平头圆柱绕流DSMC/EPSM混合算法研究

本文将EPSM算法与DSMC方法结合，构造了可模拟含近连续流区及过渡流区的DSMC／EPSM混合算法。运用混合算法模拟马赫数等于5．37时超声速平头圆柱绕流，并与DSMC结果进行比较，证明了DSMC／EPSM混合算法的有效性，并对两种算法的计算效率进行了比较。     

模拟真空羽流场的特征线法

介绍了特征线法（MOC法）对真空羽流场进行数值模拟的过程，用MOC法进行了喷管出口流动不均匀情况真空羽流场算例的求解，所得结果与理论分析相符合，与DSMC法计算结果吻合，提高了计算效率，节省了计算时间。结果表明：用MOC法可代替DSMC法对连续介质区进行模拟，并可在用N-S和DSMC耦合数值模拟中确定耦合边界。     

DSMC量子动理学化学反应模型的数值模拟

为评估直接模拟蒙特卡罗（DSMC）方法中新近发展的一种量子动理学（QK）化学反应模型,在自研高超声速稀薄气体流动数值模拟软件RariHV中实现该模型,并进行了单个绝热单元内的热泳平衡/非平衡反应测试。结果表明,数值解与理论解吻合良好,QK模型具有很好的预测化学反应流动的能力。为进一步验证QK模型在真实流动中的表现,还计算了高超声速绕圆柱的化学反应流动,QK模型计算的驻点线上温度和组分分布都与总碰撞能（TCE）模型结果吻合良好。相较于TCE模型,QK模型的优势是不再依赖反应速率系数的实验数据,在深空探测等化学反应数据缺乏的领域具有一定的应用前景。     

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

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

激光光船周围连续/稀薄混合流场的数值模拟

当火箭模式激光推进器在高空稀薄大气中飞行时,需要在激光聚焦区附近喷射稠密气体工质,从而利用激光聚焦击穿气体工质形成等离子体,为光船提供推力。为了分析该光船的推力产生过程,采用连续流体力学计算方法(CFD)和直接数值模拟蒙特卡洛方法(DSMC)耦合的混合算法,数值模拟飞行高度为80 km飞行马赫数为5的激光光船周围高温连续气体和稀薄环境大气耦合的非定常混合流场,并得到推力随时间的演化曲线。结果表明,稀薄效应明显增加了激波厚度,高温效应延迟了推力的产生,并导致冲量耦合系数和比冲的降低。