用双流体—轨道模型模拟四角喷燃模型炉内三维湍流两相流动和煤粉燃烧

目前炉内两相流动和煤粉燃烧数值模拟中多半用颗粒随机轨道模型和单流体无滑移模型，这些模型都难以完整地给出三维空间内颗粒速度，浓度，湍流度分布的信息。主采用双流体－轨道模型（颗粒相连续介质－轨道模型）对一个四角喷燃模型炉内三维湍流两相流动及煤粉燃烧进行了模拟。此模型基于欧拉气相方程组、欧拉颗粒连续方程组和动量方程组以及拉氏颗粒能量方程和质量变化的方程，并使用k-ε-kp两相湍流模型，EBU－Arrhenius湍流燃烧模型，离散坐标辐射传热模型，煤粉颗粒的水分蒸发，热解挥发模型和焦炭燃烧的扩散－动力模型等。热态模拟中，为了减小为信散造成的影响，采用了扭转坐标法（将坐标扭转一定的角度使之与煤粉射流方程一致）。为了检验数值模拟，采用三维相位移普勒测速仪（PDPA）对于冷态模型炉内湍流两相流场进行了测量，得到了两相速度，湍流脉动及颗粒浓度的分布。分别对冷态模炉内两相流动和热态模型炉内三维两相流动和煤粉进行了模拟，冷态两相流动的计算与实验结果的对比表明预报的两相流场是合理的，热态模拟的结果给两相速度，气相温度、组分浓度及壁面热，显示出靠近出口处气相速度和温度分布不对称，造成一个局部高温区。     

固体火箭发动机内气粒两相流动的SPH-FVM耦合方法数值模拟

为解决现有数值模型在求解固体火箭发动机内颗粒的燃烧流动等问题时所遇到的难题,基于拟流体模型,采用光滑粒子流体动力学方法 (SPH)求解离散颗粒相,追踪颗粒运动轨迹,采用有限体积方法 (FVM)对气体连续相进行描述,捕捉流场特性,两相间通过曳力、压力梯度、热传导、体积分数等参量进行耦合,建立了两种不同坐标系下方法间的耦合框架。耦合方法模拟了喷气推进实验室(JPL)喷管中气粒两相流动过程,得到了气相和颗粒相的参数分布规律,与相关实验结果及离散颗粒模型(DPM)数值模拟结果对比,吻合较好。针对Φ315mm实验发动机工作状态和结构特点,进行了发动机燃烧室内过载条件下的气粒两相流动数值模拟,分析了纵、横加速度载荷对发动机燃烧室内粒子场和聚集带的影响,与已有的数值模拟结果基本吻合。结果表明,新方法求解发动机气粒两相流问题准确可靠,适用于发动机内更为复杂的两相流问题数值模拟。     

超声速流场中侧向射流的数值研究

通过数值方法研究拱柱体上侧向射流与超声速来流的相互干扰,为侧向射流喷管出口横截面几何形状的设计提供依据,实现高速飞行器机动飞行的姿态控制。采用多块嵌套网格并行算法,数值求解Navier-Stokes方程以及k-ω湍流模型,分别对四种不同横截面形状(相同面积的圆形、椭圆形、楔形以及倒置楔形)的射流喷管出口形成的流场进行了研究。圆形喷管出口的流场数值计算结果与实验结果相一致,表明本文采用的数值方法可靠以及网格分布合理。对比物面压力分布和射流干扰放大因子发现,对于长轴与来流垂直的椭圆形喷管出口,在射流上游一侧的物面附近形成较大的分离区以及较强的激波,由此产生的高压区有利于侧向射流对飞行器姿态的控制。     

喷管二维跨声速两相湍流流场的数值模拟

采用ＳＩＭＰＬＥ算法,结合分散颗粒轨道（ＰＳＩＣ）模型,在同位网格上实现了喷管中跨声速两相湍流流场的数值模拟。根据Ｒｈｉｅ,Ｃｈｏｗ提出的动量插值概念,以协变物理速度分量作为求解变量,完成跨声速气相流场的数值计算。湍流模式采用ｋ－ε双方程模型,计算结果与提供的实验数据相符合。     

湍流模型在单/双缝引射矢量喷管中的验证分析研究

利用数值模拟方法,分别求解以SA、SSTk-ω、EASMk-ω或k-ζ湍流模型封闭的RANS方程,针对单/双缝引射的激波控制射流推力矢量喷管展开研究,在多个主喷流压比NPR(4.6,7.0,8.78,10.0)和次主流压比SPR(0.7,1.0)下,系统考察了四种不同湍流模型对射流推力矢量喷管内外流场的流场结构,性能参数及主喷管内壁面压力分布的预测能力,探讨了由激波控制产生的推力矢量随不同参数的变化规律。数值计算结果表明,四种湍流模型都能比较准确地预测出射流推力矢量喷管的性能参数和喷管流场结构,双缝引射的推力矢量效率并不一定优于单缝引射。     

双钟型喷管高度补偿特性的数值分析

针对双钟型喷管流场中存在的激波、激波与附面层相互作用、大尺度分离流动等复杂的物理现象，应用RNGκ-ε湍流模型封闭二维粘性可压缩N-S方程组，采用二阶迎风格式及隐式Gauss-Seidel迭代方法进行耦合求解，对双钟型喷管流场进行数值模拟研究。结果表明，双钟型喷管在低空出现壁面可控流动分离，避免了过膨胀现象，性能优于大面积比钟型喷管；在高空下为完全附着流动，性能优于小面积比钟型喷管，实现了高度补偿特性。     

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

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

气动塞式喷管的流场数值模拟和设计参数探讨

针对气动塞式喷管流场中存在的激波、激波与附面层相互作用、大尺度分离流动等复杂的物理现象 ,本文采用扩展的压强校正法求解用紊流模型封闭的可压缩湍流平均 N-S方程组 ,利用交错网格系统抹平迭代过程中的数值振荡 ,并通过高分辨率的 TVD格式进一步改善压强校正法的激波捕获效果 ,对气动塞式喷管的流场进行了详细的数值模拟。并进一步探讨了喷管倾角、塞锥长度和二次流流量比对喷管流场结构和性能的影响。     

突扩回流与大速差射流回流湍流气－固两相流动的数值模拟

本文建立了气－固两相流的多连续介质模型，并分别对突扩回流与大速差射流回流这两种与实际工程应用密切相关的复杂流动条件下的受限湍流气一固两相流动进行了数值模拟。单相湍流回流流场的模拟结果与实验符合良好。两相流动的模拟结果表明，与突扩回流相比，大速差射流回流更有利于实现较高的颗粒浓度与气相回流区及低速尾流区的匹配，增加两相滑移速度，延长颗粒停留时间及强化两相间的混合。     

固体火箭发动机喷管燃气舵铰链力矩参变分析

以固体火箭发动机喷管燃气舵绕流流场为研究对象,采用数值模拟技术,通过离散求解流动方程及相关湍流模型方程,对包含舵基、舵片和发动机喷管在内的复杂流场进行了变参数分析,获得并揭示了舵基、舵片附近的流动现象和流动特点.验证了所引入的数值模拟方法的有效性,分析结果对舵基、舵片的结构安装以及控制机构的设计具有借鉴意义.     

固体火箭发动机内轴对称两相流场的一体化数值模拟

将固体火箭发动机燃烧室和喷管作为一个整体，应用时间相关法进行了轴对称两相流场的一体化数值模拟。对气相采用ＭａｃＣｏｒｍａｃｋ显式差分格式，对颗粒相则采用特征线法进行计算。计算结果表明该方法可以对带潜入式喷管的发动机内复杂的两相流场进行有效的数值模拟。     

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

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

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

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

全耦合的一维两相喷管流的数值解

本文对控制方程组作了无因次处理,用龙格-库塔-基尔法进行了数值解。计算了五种粒度及两种微粒-气体质量流率比下微粒速度滞后及温度滞后;气相和凝相参数沿喷管长度的变化。分析讨论了所获得的结果,并简要地评论了几种一维两相喷管流模型。      

N2O/HTPB固液火箭发动机喷管两相流计算

利用二维轴对称N-S方程对选用氧化亚氮/丁羟基燃料推进剂的固液混合火箭发动机的喷管两相流进行了计算.计算采用MacCormack时间推进预报校正二步格式,采用了Baldwin-Lomax代数湍流模型和两相平衡流模型.计算了三种氧燃比下4个不同喷管的喷管流场参数,并计算了喷管性能,通过比较两相流和气相流的计算结果,分析了不同氧燃比和喷管形状对喷管性能的影响,认为固液火箭发动机的性能主要受氧燃比的影响,为固液混合火箭发动机的设计提供了依据.      

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

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

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

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

固体发动机喷管内两相流的有限元数值模拟

本文给出了固体发动机喷管内两相流的数值模似。从非守恒型两相流甚本方程出发,采用加权积分技术,导出了离散点的流动物理量的有限元数值计算公式。数值模拟分别与JPL喷管的两相流壁面静压分布测量结果以及与差分法计算结果作了比较,比较结果表明:二者是吻合的。本文还以中型喷管为例,计算了粒子的轨迹,粒子的速度滞后和燃气的温度滞后,计算结果表明;入口粒子冲击角对粒子轨迹有明显的影响。     

轴对称激波诱导矢量喷管三维两相流数值模拟

运用颗粒轨道模型对轴对称激波诱导矢量喷管内流场进行了数值模拟,并与纯气相条件下的计算结果进行了对比分析,研究了1～30μm不同直径颗粒的运动轨迹和内流场参数的分布特征,并对矢量喷管的推力性能进行了研究。结果表明,颗粒相的加入对射流推力矢量喷管内流场和推力性能都带来很大影响,影响规律与颗粒直径密切相关。为两相流条件下射流矢量喷管的试验研究提供了一定依据。     

大涡模拟超声速两相流混合层的非线性分析

为研究超声速两相流混合机理,采用大涡模拟(LES)方法数值仿真了超声速混合层内液滴两相流流场结构.用二维盒式滤波法处理可压缩气相N-S方程,得到亚格子形式的LES控制方程;用轨道模型仿真液滴运动,并用单液滴蒸发模型模拟液滴蒸发过程.利用伪邻近法和互信息量法确定了相空间重构的嵌入维数和时间延迟,对混合层压力时间序列进行了相空间重构.分析了仿真结果中混合层的涡结构特性和液滴对流场的反作用效果,并讨论了混合层失稳过程的非线性特征.发现混合层发展的不同阶段与其系统的非线性特征值密切相关.