三维化学非平衡无粘流场的数值模拟

本文用ＮＮＤ格式求解了化学非平衡流场，利用点隐式处理组分连续方程的刚性，采用全流场激波捕捉办法和特征理论，成功地进行了物面边界条件的处理，模拟了高马赫数来流条件下捻激波流场。化学反应采用了五组分五反应模型，建立起来的三维流场计算程序可用于复杂外形化学非平衡无粘流场的数值模拟。     

液氢液氧发动机喷管化学反应流数值模拟

用时间相关的半隐格式求解液氢液氧发动机喷管中的化学反应流动，考虑了６种燃烧产物，８个基元反应，得到了流动参数和产物的质量分数在流场中的分布，为发动机化学动力学损失计算和喷管热结构设计提供了必要数据。计算结果表明，采用半隐格式解决组分连续方程刚性引起的数值不稳定是有效的。     

喷管化学反应流的对角化点隐式法求解

采用弱耦合对角化点隐式方法的ＭａｃＣｏｒｍａｃｋ格式,求解喷管的粘性化学反应流动,解决了控制方程组的刚性问题,又大大减少了常规点隐式方法中求逆矩阵的大量工作。对液氢液氧火箭发动机,采用６种组分、８个反应有限速率的化学反应模型和Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流场参数在喷管中的分布。计算结果表明,采用对角化点隐式方法求解化学反应喷管粘性流动能提高计算效率。     

非平衡流解耦方法及其计算激波诱导燃烧的应用验证

介绍并发展了一种改进的化学非平衡流动解耦方法.将流动控制方程与化学反应生成源项解耦处理,组分对流项的矢通量分裂基于流动方程密度对流项的分裂形式给出.对流项采用五阶WENO格式求解,化学反应源项的刚性采用简化的隐式方法进行处理.将该方法应用于预混H₂/O₂和H₂/Air激波诱导燃烧计算,得到的定常及非定常燃烧过程与实验观测相符,振荡频率与实验测量较吻合,结果表明该方法应用于多组分多步反应流计算是可行的.      

预爆管设置方式对主爆轰点火的影响

应用含组分守恒方程的Navier-Stokes方程,对H2+O2预混气体的二维爆轰进行了数值计算.在建立方法、及验证的基础上,对比研究了化学动力学模型、预爆轰管几何位置对爆震燃烧室流场及引爆性能的影响.在计算过程中,对方程中的对流项采用了空间为二阶精度的AUSMPW+格式,扩散项则采用了二阶中心差分离散.采用流体力学时间步内的当地积分法克服了非平衡源项的刚性问题,保持了LU-SSOR隐式算法的高效率.非定常计算采用了双时间方法.预爆轰管与主爆震室平行,采用两种设置方案:一种正向安装推力壁上,另一种在推力壁附近与其相对设置.研究结果表明,反射激波及过爆轰驱动在爆轰点火和传播中起重要的作用,与推力壁相对设置的预爆轰管更有利于爆震的引发.     

固液混合火箭发动机喷管流动计算

利用二维轴对称 N-S方程和组分方程对选用液氧 /端羟基聚丁二烯推进剂的固液混合火箭发动机喷管流动进行了计算。计算采用 LU时间隐式格式、MUSCL空间离散和 Van Leer矢通量分裂方法 ,8组分 1 0反应的化学反应模型 ,对化学源相进行了点隐式处理。喷管入口条件通过燃速计算和热力计算得到 ,计算了多个氧化剂流率和装药长度情况下的喷管流场 ,分析了真空比冲和推力与氧化剂流率的关系 ,为固液混合火箭发动机喷管设计提供了依据      

跨音速喷管化学反应流数值模拟

用时间相关的半隐格式有限差分数值方法求解了化学非平衡反应跨音速喷管流场，在喷管收敛段，流动接近化学平衡状态，控制方程的刚性问题严重，数值积分困难。通过对时间差分项隐式离散、对空间差分项显式离散，流场边界采用参考平面上的特征线法计算等，成功地解决了由于化学反应有限速率带来的数值解不稳定问题。该格式简单、需要计算机存贮空间少。本文完成了一维和轴对称非平衡化学反应喷管流动计算，并与化学平衡流和冻结流的计算结果做了比较。     

冲压燃烧室内的燃料扩散性能研究

应用含组分守恒方程的质量平均Navier-Stokes方程和B-L代数湍流模型,本文数值模拟了带支板引射器的模型燃烧室内部的流场结构,对比研究了支板引射器数量以及燃料注入角度对燃料混合性能的影响.在计算过程中,对方程中的对流项采用了空间为二阶精度的TVD和AUSMPW格式,扩散项则采用了二阶中心差分离散.采用流体力学时间步内的当地积分法克服了非平衡源项的刚性问题,保持了LU-SSOR隐式算法的高效率.研究结果表明,燃料有角度注入及采用多支板构型可提高燃料的混合性能.     

高超声速化学非平衡粘性激波层方程的数值解

本文讨论了高超声速粘性激波层方程数值计算时差分格式引起的物理失真问题。具体分析了全隐格式格式粘性的影响,并作了数值试验。为了验证隐式结果的可靠性。在超声速激波风洞中测量了钝锥的表面压力分布,并与计算结果作了比较,两者基本一致。 本文采用隐式有限差分法数值计算了高超声速化学非平衡粘性激波层绕细长球锥的流动。计算时采用连续方程和法向动量方程耦合求解的方法以解决细长体远后身区计算中的问题。应用网格技术和加强系数矩阵主对角元素优势的方法提高了化学非平衡流计算的雷诺数范围。文中给出了高超声速化学非平衡流的计算结果,并与其它文献的结果作了比较。     

碳氢燃料超声速燃烧的化学动力学研究：（Ⅰ）甲烷超声速燃烧的简化化？…

采用雷诺平均Ｎ－Ｓ方程对碳氢燃料甲烷在超声速流动中混合燃烧过程进行了数值研究邮简化的甲烷燃烧的反应机理及反应动力学模型,建立了Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数涡粘性湍流模型,采用了高精度的ＴＶＤ隐式格式和化学源项隐式解法,计算结果和氢气／空气超燃的计算结果进行了比较表明所采用的反应机理及计算模型合理,算法可行,就算例结果而言碳氢化合物作为一种超前其的优势,值得进一步深入研究。     

化学氧碘激光(COIL)三维混合反应流场数值模拟研究

通过求解层流Navier-Stokes方程和组分输运方程,数值模拟了化学氧碘激光(COIL)三维扩压器和光腔中10组分21反应的内流流场,计算中没有考虑水蒸气在低温下的冷凝以及激光能量的输出。通过与国外文献计算结果比较可知,物性系数、化学反应速率常数、边界条件处理等对组分浓度影响较大,缺乏这些数据库将会影响计算结果。计算表明,三维副流射流基本穿透主流到达中心线,主副流混合较强,激发态粒子高浓度区主要在中心线附近,因此,也只有在中心线附近有较强的小信号增益系数,进一步研究混合增强技术可以提高原材料的使用效率以及激光的出光功率。本研究可为强激光武器系统之新型压力恢复系统的设计提供技术储备。     

化学非平衡流动解耦算法的回顾与新进展

超声速化学非平衡流动的数值模拟一直是计算流体力学领域的难点,主要体现在如下几个方面：物理过程非常复杂,存在着激波、燃烧波等各种复杂波系的相互作用;超声速化学非平衡流动属于典型的时空多尺度物理问题,其控制方程存在严重的刚性,给数值求解带来了很大困难。对国内外的解耦算法研究现状简单综述后,主要介绍1993年刘君提出的解耦算法的理论基础,流动方程采用冻结流模型,源项方程模拟流体微团在当地绝热、定容的热力学系统内发生的化学反应过程。通过引入两个中间变量,即等效内能和等效比热比,将与温度无关的生成焓从流动方程组能量项中分离出去,源项方程组中包含等效内能,使用不同算子对流动方程和源项方程解耦求解。与传统解耦算法相比,源项方程的求解过程中包含状态参数和组元同时变化。结合刘君解耦算法机理和有限体积法空间平均特性,介绍近期在提高算法计算效率方面的研究进展,包括流动方程优化算法和耦合过程优化。采用优化算法对经典的激波诱导燃烧算例进行数值模拟,与不同文献结果进行对比,验证了优化算法的时空精度。通过总结经验发现化学非平衡反应仅发生在流场局部区域,提出质量生成率判据,结合相应模拟结果验证了方法的可行性,可以进一步提高化学反应算子的计算效率。     

电弧喷射推力器化学非平衡数值模拟

为了准确把握电弧喷射推力器工作过程物理机理与特征,采用化学非平衡模型对不同压缩室直径下低功率氮电弧喷射推力器工作情况进行了数值模拟。模型考虑了工作过程中的分解反应、电离反应和复合反应,化学动力学模型为4组分,4反应的有限速率化学反应模型。采用二阶精度NND格式数值求解耦合电磁源项和化学反应源项的N S方程组,采用有限控制容积积分方法求解由麦克斯韦方程组推导出的电磁场方程。数值模拟的结果揭示了推力器内部电弧能量输入作用和高温电离气体的离解电离状况,给出了不同压缩室直径下推力器的推力、比冲和推进效率。结果分析表明,压缩室直径对推力器性能具有较大影响。     

共用尾喷管多管脉冲爆震发动机数值模拟研究

采用基元反应模型和显式迎风TVD差分格式数值模拟方法,研究共用尾喷管多管脉冲爆震发动机(Pulse Detonation Engine,简称PDE)各爆震管之间相互影响.采用点隐算法解决化学反应引起的刚性问题,模拟了共用尾喷管内的流场,以及爆震波退化为激波遇到尾喷管收敛斜面反射及反射波向爆震管上游传播的过程.研究结果为多管PDE的共用尾喷管设计提供了理论基础.     

时空守恒元和解元方法的爆震波一维数值模拟

运用基元化学反应模型和流体力学的基本原理，建立了脉冲爆震发动机一维非定常反应流数学模型。应用一种改进的时空守恒元和解元(CESE)方法建立控制方程的差分格式，针对化学反应源项的刚性问题，在对非线性源项进行线性化处理的基础上，应用隐式梯形方法实现源项的积分。对氢氧爆震波的形成和传播过程进行数值计算，结果表明，应用CESE方法和基元化学反应模型能够成功地模拟氢氧爆震波的精细结构。     

跨声速操纵面嗡鸣的数值研究

我国的一些高速歼击机，在飞行中曾多次出现操纵面嗡鸣问题，它严重影响了飞机的性能和飞行安全，本文采用数值计算方法，对此问题做较深入的分析和研究。气动力计算采用Ｎａｖｉｅｒ－Ｓｔｏｋｅｓ方程，并采用ＬＵ－ＴＶＤ混合差分格式求解由此气动力解操纵面的弹性振动方程，求出现嗡鸣的Ｍ数及振动的频率、振幅等参数。计算中采用Ｂ－Ｌ代数湍流模型，网格生成采用代数方法及保角变换方法，文中讨论了隐式算法在非定常气动力计算     

激波绕射双方块加速诱导爆轰的数值研究

基于二维单步化学反应Euler方程,应用Roe/HLL混合格式以及自适应网格加密(AMR)技术,对方管中弱激波绕射双方块障碍物加速诱导爆轰现象进行了数值研究。计算结果表明,弱激波与双方块相互作用并发生绕射、碰撞可有效加速诱导爆轰,且在管道阻塞比为45%左右时可达到最佳效果,与相关实验结果相符。另外,研究了方块个数及排列方式对激波绕射加速诱导爆轰的影响,得到了最佳方块数量与方块排列方式。     

用CSCM格式求解激波反射问题

本文对Euler方程采用CSCM格式,数值求解了激波反射问题,与目前流行的交替方向隐式(ADI)格式以及总变差递减(TVD)格式的计算结果进行了比较,并作了简要分析讨论。