固体火箭发动机喷管气固两相流动的数值模拟

对颗粒相采用颗粒轨道模型,气相求解可压缩N-S方程组,计算方法采用显式Runge-Kutta时间推进法与有总变差衰减(TVD)性质的高精度MUSCL-Roe格式;自主开发了曲线坐标系下二维轴对称可压缩N-S方程组的解算器Solve2D,研究了固体火箭发动机喷管中颗粒相对流场的影响以及不同尺寸颗粒运动规律.结果表明:颗粒相对流场的影响主要表现在喷管喉部以及扩张段,和单相流场相比,沿轴线马赫数减小,且颗粒尺寸越小减少得越多;沿轴线气相温度升高,且颗粒尺寸越小温度升高越多;颗粒尺寸越小,无粒子区越小;颗粒越大与收缩段壁面碰撞越剧烈,无粒子区越大.      

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

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

高温高压超声速气液两相流场数值模拟

采用欧拉-拉格朗日方法对燃气-蒸汽发射动力装置内高温高压超声速燃气中横向喷雾的气液两相流进行了数值模拟研究,气相采用RK(Runge-Kutta)-AUSM+(advection upstream splitting method+)格式求解Navier-Stokes方程,液相应用颗粒轨道模型,两相之间的耦合通过在气相...     

变截面管道中激波诱导的两相流动

本文研究扩张型与收缩型变截面管道中激波诱导的非定常两相流动。对于稀相气固悬浮体采用了双连续介质模型,对于两相流动按照准一维近似处理。控制方程则是利用算子分裂技术和二阶GRP方法数值求解。文中给出了波后气固两相的流动结构,并讨论了粒子对两相流场的影响。     

大功率氧煤燃烧器内湍流气固两相流动数值模拟

本文对新型同轴射流缩口式氧煤枪所组成的大功率氧煤燃烧器内冷态轴对称二维湍流气固两相流动用ｋ－ｅ－ＡＰ两相湍流模型进行了数值预报，揭示了燃烧器内两相流场特征，预报结果和实验结果相符。进一步对相应高炉实际工况高速条件下的冷态两相流场进行了预报，其结果为优选结构尺寸和流动参数提供了重要依据。     

水中弹药引信外置涡轮高速旋转两相流动数值仿真

为预测空化条件下的涡轮转速,获得涡轮叶片附近空化情况、压力系数、阻力系数等数据,采用SchnerrSauer空化模型,进行了水中弹药引信涡轮高速旋转两相流动的三维流场数值仿真研究。计算得到不同水深(10m、30m及50m)和弹速(5~100m/s)条件下的流场数据。结果表明:引信涡轮表面在弹速30m/s左右时空化初生,且空化尺度随着弹速的增加而增大;弹速100m/s时,头部粘性阻力占总阻力的18.7%;涡轮转速与水深、静压、弹速相关,高速时涡轮转速随水深减小明显,与弹速基本成正比线性关系,非线性度约为0.998。     

基于气固两相双流体模型研究火箭发动机斜切喷管流场特性

为了研究火箭发动机(SRM)斜切喷管的两相流动特性,采用气体-颗粒相双流体模型,并结合多区域混合网格技术,对发动机斜切喷管内气相与颗粒相的相互作用规律进行研究,探索颗粒直径与颗粒质量分数变化对发动机喷管气固两相流动特性的影响。结果表明:固体颗粒相的存在,对发动机斜切喷管的流场结构产生重要影响,导致喷管轴线附近存在一个燃气流动速度较低,温度较高的区域。同时,喷管壁面附近存在无粒子区,随着颗粒直径的增加,无粒子区域的范围逐渐扩大。并且,颗粒直径越大,其运动速度越小,在喷管内的滞留时间越长。颗粒直径与质量分数的变化同样会影响发动机喷管的流场结构,随着颗粒直径的增加,发动机喷管轴线处气相马赫数先减小后增大,而燃气温度则先增大后减小;发动机推力的变化趋势与马赫数变化趋势相同,但两者并不同时达到极值点。颗粒相的质量分数越大,沿喷管轴线方向的气相马赫数和发动机推力越小,喷管两相流损失越大。     

可压稀相两相流场的数值模拟

针对可压两相流动，在阻力系数中考虑了可压缩性及颗粒聚积所造成的影响，用TVD格式对沙层，平板及JPL喷管内的气固稀相两相流进行了数值模拟，所获结果和已有的实验数据及有关文献吻合很好，在JPL喷管中，随着质量分数的增加，由于气体和颗粒间的相互作用逐渐增强，气相的温度增大，马赫数减小，但压力变化很小，而随着颗粒直径的增加，气相轴线上的温度减小。     

前后翼型装药结构SRM气固两相流动特性

基于可压缩流动强守恒型Navier-Stokes方程,运用颗粒轨道模型(PTM)和shear stress transfer(SST)k-ω湍流模型,采用计算单元内颗粒源法(PSIC)进行气固两相耦合计算,建立了某型前后翼型装药结构固体火箭发动机(SRM)工作初期的三维两相内流场数值计算模型.对比分析了纯气相和气固两相...     

气／固两相流测量系统

面粉机提粉管中有流动着的气体和面粉，称之两仃流，采用由二个极板组成的电容在和电容－电压变换电路组成的电容式传感器测量管道中流动着的面粉密度。用压力传感器测量管道中的风压并计算出风速，由风速和测得的面粉密度可算出管道的面粉流量，根据流量大小，计算机可控制每个磨口的工作状态。     

H2引燃雾化煤油超燃混合的数值研究

为研究H2引燃的煤油超燃混合问题，采用流体模型描述两相流动，气相反应系统的全N－S方程风迎风TVD格式求解，液相扩散Euler方程用NND格式求解，化学反应源项采用点隐处理，对单喷嘴喷H2超声速燃烧、H2引燃的煤油超燃混合问题进行了数值研究。获得了压力、密度、温度和组元浓度场的分布。结果表明：与单喷嘴喷H2相,比下游喷嘴H2的射流穿透浓度和扩散范围增大，燃烧区也变大；虽然煤油穿透深度大，但其扩散较H2差，且下游出现无煤油区。     

计算两相流的一种改进方法

针对两相流动中密度及颗粒浓度的变化，建立了一种考虑颗粒相体积浓度值变化的颗粒相压力修正方程及考虑可压缩性影响的流体相压力修正方程，在不可压两相流的计算方法上，发展了一种计算可压两相流的方法，以速度协变分量为求解变量，对不可压两相密相流及可压稀相两相流进行了计算，并和已有的实验数据及有关文献进行了比较，吻合性较好，表明计算方法可行。     

H_2引燃雾化煤油超燃混合的数值研究

为研究H2 引燃的煤油超燃混合问题 ,采用流体模型描述两相流动 ,气相反应系统的全N S方程风迎风TVD格式求解 ,液相扩散的Euler方程用NND格式求解 ,化学反应源项采用点隐处理 ,对单喷嘴喷H2 超声速燃烧、H2 引燃的煤油超燃混合问题进行了数值研究。获得了压力、密度、温度和组元浓度场的分布。结果表明 :与单喷嘴喷H2 相比 ,下游喷嘴H2 的射流穿透深度和扩散范围增大 ,燃烧区也变大 ;虽然煤油穿透深度大 ,但其扩散较H2 差 ,且下游出现无煤油区。     

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

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

气液喷流与超声速凹腔流场的相互作用

为了解凹腔火焰稳定器的工作过程,通过试验和数值仿真手段,对超声速条件下凹腔的流动特性进行了深入细致的研究,探讨了气/液喷流与超声速凹腔流动的相互作用机理。研究结果表明,无喷流时超声速凹腔流场具有五个典型的特征;引入喷流会引入新的流场特征,同时流场结构会发生机理性的变化;气/液喷流对凹腔整体流动特征的影响是一致的;有气体喷流时,不同压降下凹腔流场结构是类似的;而对液体喷流而言,提高喷注压降会增加雾化距离、射流穿透度、喷流厚度,提高来流速度则会使其减小。     

星孔型装药发动机三维两相流场的数值模拟

为了研究颗粒在星孔型装药固体火箭发动机燃烧室和喷管中的运动轨迹以及颗粒与发动机壁面的碰撞情况,针对可压两相流动,采用了高雷诺数下的k ε湍流模型和欧拉 拉格朗日两相流模型,用全速度SIMPLE方法对方程组进行求解,并用PSIC方法进行气固耦合计算。计算得出了流场内两相的速度、温度等参数的分布及多种情况下固体颗粒的运动轨迹。在燃气生成量确定的情况下,从距离喷管较近的某些位置进入流场的颗粒比较容易撞击壁面;颗粒的尺寸和局部产生的旋涡对颗粒的轨迹和碰撞也会产生较大的影响。     

二维凹腔超声速流动的混合RANS/LES模拟

基于Menter的k-ω SST湍流模型构建了一种混合RANS/LES模拟方法,通过采用一个与到壁面距离相关的衔接函数将处理近壁区的SST模型过渡到处理主流区的Yoshizawa一方程亚格子模型.利用此方法对用于超燃冲压发动机的凹腔的二维超声速流动进行模拟,模拟结果再现了二维凹腔剪切层的拟序结构,计算得到的凹腔自激振荡频率、时均统计的压力分布以及压力剖面与实验结果吻合较好.     

煤油超燃冲压发动机两相流场数值模拟(I)数值校验及总体流场特征

针对一种带三维导流型凹槽的煤油超燃冲压发动机,考虑两相(颗粒随机轨道方法)、湍流(Menter的SSTk ω模型)和多步燃烧机理(代用燃料C12H23的17组分30步化学反应机理),对飞行马赫数M0为4和6的典型工况开展了初步数值研究。结果流道壁面压强分布与试验值取得了较好的吻合,不过M0=6工况的第一凹槽组附近的存在约23%的数值误差。给出了误差原因。进一步分析表明(1)M0=6工况的分段燃料喷注对整体燃烧产生了有利影响,导致了较高的燃烧效率。(2)M0=4工况采用了不合理喷注参数,结果煤油穿透深度不足,获得了较低的燃烧效率。(3)三维导流型凹槽附近的湍流效应明显,对掺混产生了有利影响。不过超声速主流的湍流强度远大于亚声速主流情况。     

Fluid structure interaction simulation of supersonic parachute inflation by an interface tracking method

An Arbitrary Lagrangian–Eulerian(ALE) approach with interface tracking is developed in this paper to simulate the supersonic parachute inflation. A two-way interaction between a nonlinear finite element method and a finite volume method is accomplished. In order to apply this interface tracking method to problems with instantaneous large deformation and self-contact, a new virtual structure contact method is proposed to leave room for the body-fitted mesh between the contact structural surfaces....