液体火箭发动机有无隔板的热态解数值模拟

用数值方法模拟了液体火箭发动机燃烧室内有无隔板两种情况的热态解。气相控制方程用欧拉坐标系下的Navier-Stokes方程组描述,液相控制方程在Lagrangian坐标系下进行描述。      

湍流两相燃烧不稳定性分析

采用数值模拟的方法对液体火箭发动机燃烧室内湍流两相流动与燃烧过程进行了分析,同时研究了燃烧室内高频不稳定燃烧现象。气相控制方程用欧拉坐标系下的Navier-Stokes方程组描述,液相控制方程在Lagrangian坐标系下描述,湍流模型采用高雷诺数的k-ε双方程模型,并采用高压蒸发燃烧模型。通过数值模拟分析了当燃烧室处于稳态燃烧时从喷注面到喉部之间两相混合物的变化规律;当发生高频不稳定燃烧时所研究的控制体内混合气体和喷雾的量发生振荡,并且具有相同的频率和半个周期的相位差。     

液体发动机燃烧室流场模拟的并行SIMPLE算法

用串行和并行ＳＩＭＰＬＥ算法对液体火箭发动机燃烧室内的复杂流动进行了数值模拟。控制方程组用欧拉坐标系下的Ｎ－Ｓ方程组描述。在并行虚拟机环境下研究了交错网格系统中参数的传输和接受模式, 从计算结果可以看出并行计算的效率较高, 是解决此类复杂流动大规模数值模拟问题的有效手段之一。     

旋转部件复杂表面水滴撞击计算

为模拟旋转体结冰问题的过冷水滴运动及撞击过程,基于欧拉方法及单旋转坐标系模型,建立了三维旋转水滴运动模型,并提出了相应的数值求解方法。采用单旋转坐标系对空气及水滴两相流动过程进行处理,通过引入惯性力,将惯性系下的周期转动边界转换为定常流动边界;利用欧拉方法,使用单向耦合形式描述空气—水滴流场;在单旋转坐标系下,向控制方程内引入科里奥利加速度及牵连加速度,进行非惯性系下欧拉方程的修正,从而描述水滴运动过程;采用有限容积求解器对空气及水滴运动的控制方程组进行求解,通过引入源项定义单旋转坐标系下的惯性力,得到空气流场及水滴场的速度、体积分数分布,进而得到表面水滴撞击特性。采用上述方法对旋转帽罩与叶片模型进行算例分析,结果表明所建立的旋转水滴计算方法有效,对比静止状态表面的结果,旋转对帽罩的水滴撞击特性影响甚小,而对桨叶存在显著影响;由于帽罩具有中心对称的特性,因而旋转带来的切向速度变化对其水滴撞击特性影响不明显;桨叶表面水滴收集系数随旋转角速度增大而增大,同时收集系数在表面的分布会向迎风方向偏移,较大的角速度对应了更为显著的收集系数增幅与偏移现象。     

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

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

有隔板和无隔板的液体火箭发动机冷态全流场数值模拟

针对液体火箭发动机燃烧室内有隔板和无隔板的情况，用数值模拟的方法对发动机冷态全流场进行了计算，探讨了燃烧室内有隔板时边界条件的处理、计算区域的离散以及隔板对流场参数的影响等问题，用欧拉坐标系下的Navier－Stokes方程组描述气相控制方程，通过数值模拟成功地对燃烧室内有隔板和无隔板的流场进行了数值仿真，表明了对隔板体的数值处理方法的可行性。     

液体箭发动机燃烧不稳定性分析模型

通过一个简化的物理模型。运用计算流体动力学分析方法，研究了液体火简发动机的不稳定燃烧现象，验证了该模型及分析方法的简捷和有效性，描述发动机燃烧室内气相流动过气液两控制方程是非稳准一维欧拉方程组，液相则用拉格朗日方法描述，方程组中的源项反映了气液两相耦合及液雾蒸发燃烧的影响。针对不稳定燃烧的机理，考察了压力脉动随着时间的推进对声学扰动的敏感性，计算格式采用两步积分的预测校正法进行迭代求解。     

液体火箭发动机燃烧不稳定性分析模型

通过一个简化的物理模型，运用计算流体动力学分析方法，研究了液体火箭发动机的不稳定燃烧现象，验证了该模型及分析方法的简捷和有效性。描述发动机燃烧室内气相流动过程的控制方程是非稳态准一维欧拉方程组，液相则用拉格朗日方法描述，方程组中的源项反映了气液两相耦合及液雾蒸发燃烧的影响。针对不稳定燃烧的机理，考察了压力脉动随着时间的推进对声学扰动的敏感性。计算格式采用两步积分的预测校正法进行选代求解。     

化学非平衡解耦算法精度分析及其改进

模拟化学非平衡流问题涉及到流动方程和化学反应方程两部分的求解,流动方程组中包含所有组元的偏微分方程,导致求解变量成数量级增加,点隐算法和全隐算法在求解源项刚性问题时,通常会涉及到矩阵求逆运算,这些因素带来的巨大计算量限制了隐式算法在复杂工程问题中的应用。1993年刘君提出了化学非平衡流解耦算法,将控制方程组分解为流动和化学反应两部分,流动方程的求解采用冻结流假设来描述流体微团沿流线的运动过程,化学反应方程的求解描述流体微团在随体坐标系下发生绝热、定容的爆炸过程。结合解耦算法和有限体积法的特点对这种解耦算法进行改进,提出的优化算法不需要求解组元变量所对应的偏微分方程组,只求解由5个基本变量构成的,形式上与量热完全气体近似的偏微分方程组,通过对比计算结果发现优化算法可以显著地提高计算效率。同时,将精细积分方法应用于化学非平衡流问题的求解中,通过与传统的VODE方法和α-QSS方法对比发现,精细积分方法的鲁棒性更优、精度对时间步长不敏感,适当的选取时间步长,可以充分发挥精细积分方法的优势。     

一种基于执行坐标系的导弹运动建模方法

为解决导弹制导控制系统建模中弹体模型与系统的匹配问题,提出了一种在弹体执行坐标系上构建导弹六自由度运动模型的方法。首先在弹体执行坐标系对导弹进行受力分析,建立导弹动力学方程组,然后利用弹体执行坐标系与弹体坐标系的关系,直接列写运动学方程组。该方法简化了弹体建模过程中气动数据的处理,统一了动力学和运动学方程组以及制导控制方程组的参考坐标系,从而使制导控制系统建模过程更为简单,同时有效提高了仿真效率。最后,通过仿真验证了该方法的有效性。     

回流燃烧室燃烧过程的三维数值模拟

在三维任意曲线坐标系下数值模拟回流燃烧室火焰筒内两相燃烧过程,采用RNG k-ε模型模拟紊流粘性,EBU-Arrhenius模型模拟燃烧反应速率、离散坐标模型以及六通量模型考虑辐射传热,液相采用颗粒轨道模型,气相采用SIMPLE算法求解,并用PSIC算法考虑气液两相之间的相互作用的影响,计算得到燃烧室内速度、温度等各气流参数分布.通过将计算与实验结果对比表明,计算方法可靠,离散坐标模型优于六通量辐射模型,更适用于模拟火焰筒内两相燃烧流场.      

液氢/液氧火箭发动机喷雾燃烧过程三维数值模拟

用三维湍流Ｎ－Ｓ方程及拟流体模型描述液体火箭发动机内部喷雾两相燃烧流动过程,两相之间的质量交换采用蒸发模型计算,气相化学反应速率采用Ａｒｒｈｅｎｉｕｓ公式计算。采用ＳＩＭＰＬＥＲ与ＩＰＳＡ算法求解两相控制方程并对氢氧发动机燃烧流动过程进行了三维数值分析,得到了发动机推力室内燃气参数的详细分布情况。结果表明对氢氧发动机而言,蒸发过程非常快,推进剂混合过程是燃烧过程的决定因素。     

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

采用欧拉-拉格朗日方法对燃气-蒸汽发射动力装置内高温高压超音速燃气中的横向喷雾的气液两相流进行了数值模拟研究，气相采用RK-AUSM+格式求解N-S方程，液相应用颗粒轨道模型，两相之间的耦合通过在气相各守恒方程中添加源项来实现。针对不同的喷水孔径、喷水压差以及水汽质量比进行了数值模拟计算，分析了各种条件下的液滴蒸发情况以及对气相流场和能量输出的影响，为燃气-蒸汽发射动力装置的设计优化和能量调节技术提供了理论依据     

雾化燃料在超声速气流中横向喷射混合的数值模拟

本文用双流体模型对液态雾化燃料在直简形燃烧室中横向喷射问题进行了数值研究,气相用二阶迎负ＴＶＤ格式求解三维全Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程,液相用二阶预估、校正ＮＮＤ格式求解三维Ｅｕｌｅｒ方程。描述相间相互作用的常微分方程用预估、校正的二阶Ｒｕｎｇｅ－Ｋｕｔｔａ法求解。分别考察了喷口处气相压力、相间速度滑移、液滴直径和燃烧室入口处预置的超声速流向涡对雾化燃料混合效果的影响。结果表明：气相较液相的扩     

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

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