超声速补燃的数值模拟和实验研究

利用富燃料固体推进剂燃烧产生的高温燃气与侧喷空气混合进行了超声速补燃实验。运用有限体积法求解轴对称N－S方程和组分方程，对实验流场进行了数值模拟，计算中采用了LU时间陷式格式、MUSCL空间离散方法、12组分16反应的化学反应模型，对化学反应源项进行了点隐式处理。结果表明：实验方法可行，计算结果合理、可靠。     

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

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

固体火箭冲压发动机二次燃烧室流场数值计算和试验研究

用二维k-ε湍流模型及简单一步无限快速化学反应模型，对固体火箭冲压发动机二次燃烧室反应流场进行了数值计算，并针对硼推进剂燃烧特性，提出了二次进气的燃烧室设计方案，在此基础上做了试验研究。结果表明，经数值模拟设计的二次燃烧室构型比普通二次燃烧室构型燃烧效率明显提高。     

超声速补燃的三维数值模拟和实验研究

利用贫氧固体推进剂燃烧产生的高温燃气与侧喷空气混合进行了超声速补燃实验。运用有限体积法求解三维N-S方程和组分连续方程组，对实验所涉及的化学反应流场进行了数值模拟。通过考察两种掺混器的穿透深度、流向涡强度等参数，比较了二者的混合及燃烧增强效果，分析了原因。结果表明在所研究的工况下，桶型掺混器性能更优。     

火箭冲压发动机反应流场数值模拟

用偶合算法求解三维定常椭圆型Ｎ－Ｓ方程组,计算了火箭冲压发动机内部的湍流反应流场。用有限体积法离散控制方程,用块隐式法求解离散后的代数方程组。湍流流动用ｋ－ε模型描述,湍流燃烧采用掺混控制的局部瞬时不混合模型。实践表明,尽管密度场变化剧烈,仍能快速收敛。算例中网格点数为３９×１３×１８。计算出的燃烧效率符合实验结果。     

不同飞行状态下固体火箭发动机尾喷焰数值研究

为了研究飞行状态对固体火箭发动机尾喷焰的影响,建立了含化学反应项和组分输运项的N-S(Navier-Stokes)控制方程,并采用MUSCL(monotonic upstream-centered scheme for conversation laws) Roe格式的有限体积法进行求解.在此基础上,结合热力计算结果对某固体火箭发动机在不同飞行状态下的化学反应流和单一组分流进行仿真计算.结果表明:初始倾角和第一马赫波节长度随飞行高度的增加而增大,随来流马赫数的升高而减小.同单一组分流相比,化学反应流对第一马赫波节的影响较初始倾角显著.射流轴线温度因飞行状态的不同有较大涨落差异.初始倾角、第一马赫波节由来流总压决定,射流轴线温度由激波结构、复燃效应、来流条件等因素共同决定.      

用灵敏度法简化化学反应动力模型

介绍了用于简化化学反应动力模型的灵敏度方法的基本思想与具体方法，然后用该方法来对两个化学反应实例：甲醛氧化反应和一氧化碳／氢气在空气中燃烧的化学反应进行简化，均得到了比较好的结果。     

一维喷管化学反应流的吉尔方法研究

建立了一维喷管化学反应流计算模型,并进行了数值求解,给出了喷管内一维化学反应流计算结果。采用吉尔方法解决控制方程组的刚性问题获得满意效果。计算中解决了化学反应流计算中初值难于给定和喉部条件难于确定等问题。并通过对氢氧发动机喷管的分析计算,证明了这种方法的正确性。     

喷嘴结构对非壅塞式固冲发动机结构匹配性的影响

采用雷诺平均N-S方程和k-ε双方程湍流模型，研究了燃气发生器喷嘴构形对补燃效率的影响。通过一次流、二次流流线和燃烧效率分析，提出了一次流切人点的概念。结果表明：按照良好的燃气切人点设计燃气发生器喷嘴，不仅可以获得较高的燃烧效率，而且可以获得对补燃室热结构有利的温度场分布，实现固冲发动机的结构匹配性。     

支板火箭引射冲压发动机引射模态燃烧流动 (Ⅰ)瞬时掺混燃烧流场的数值模拟

针对一种适用于多种工作模态的支板火箭引射冲压组合发动机构型，Mα在0—1来流范围的瞬时掺混(SMC)引射燃烧流场进行了数值模拟。详细分析了低速模态SMC湍流流场的流动结构及特征，并对其Mα0—1范围的性能进行了分析，结果Mα在大于等于0．7范围内获得了推力增强。结论认为该种构型的组合发动机适用于作为机载导弹的动力装置，而更低马赫数范围(包括Mα＝0)内的推力增强取决于多种因素的优化匹配。     

模型加力室大涡模拟的并行计算

采用基于MPI(消息传递库)的并行算法,在贴体网格下对带V形槽稳定器模型加力燃烧室紊流化学反应流场进行数值模拟,湍流模型采用k方程亚网格尺度模型,燃烧模型采用亚网格EBU模型,采用热通量辐射模型估算辐射通量。在程序设计中,采用动态内存分配、分区算法和多点重合交错网格系统,并行计算的结果与单机运行结果的对比表明计算结果是正确的,可以明显的提高运算效率,是解决复杂燃烧流动大规模数值模拟的有效手段。      

预处理AUSM+格式在全速域反应流模拟中的应用

对预处理修正的AUSM+(advection upwind splitting method+)格式进行了详细讨论, 并将其应用于全速域范围内的化学反应流场数值模拟.基于多组分Navier-Stokes方程, 结合控制方程时间导数项的预处理, 采用预处理修正的AUSM+格式发展了一种模拟全速域范围内化学反应流场的有效算法.验证算例包括超声速、亚声速以及极低速情形的化学反应流场, 计算结果表明基于预处理AUSM+格式发展的算法能够有效地用于全速域范围内化学反应流动问题的模拟.      

试验研究旋流数对燃烧室气动性能的影响

为了深入了解双级轴向旋流器的旋流数变化对模型环形燃烧室内气流结构与气动性能的影响,采用粒子图像速度仪(PIV)对燃烧室内冷态和液雾燃烧流场进行测量,试验研究在冷、热态情况下不同旋流数对模型燃烧室内回流区的尺寸、平均速度(u,v)、脉动速度(U<,rms>,V<,rms>)以及雷诺剪切应力u'v'分布的影响.研究结果表明...     

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

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

应用PIV技术测试模型环形燃烧室流场

采用粒子图像测速仪(PIV)测量带双级旋流器模型环形燃烧室内冷态和燃烧流场,试验研究不同进口气流温度对其冷态和燃烧流场内回流区的形成以及气流速度分布的影响,试验结果表明:燃烧室燃烧流场内回流区长度要比冷态流场短,燃烧流场的脉动速度明显大于冷态流场;另外随着进口气流温度增加,燃烧室冷态和燃烧流场的回流区长度都稍有减少.      

喷管化学反应粘性流场的数值模拟

用时间相关法的显格式求解氢氧发动机喷管的粘性化学反应流动,采用６种组分,８个基元反应有限速率的化学反应模型,湍流模型采用Ｂａｌｄｗｉｎ－Ｌｏｍａｘ代数湍流模型,得到流动参数和产物的质量分数在喷管流场中的分布。计算结果表明,采用显格式方法求解喷管化学反应粘性流动是可行的。     

三维贴体坐标系下燃烧室热态流场的大涡模拟

大涡模拟三维贴体坐标系下环形燃烧室火焰简热态紊流瞬态流场。利用椭圆方程方法生成三维贴体网格，计算中采用k方程亚网格尺度模型估算亚网格紊流粘性；亚网格EBU燃烧模型估算化学反应速率；热通量辐射模型估算辐射通量。并在非交错网格系下采用sIMPLE算法和混合差分格式求解离散方程，利用壁面函数处理固壁边界条件。计算结果与实验结果的比较表明，采用大涡模拟方法能更真实反映环形燃烧室火焰简内紊流化学反应流气流结构和燃烧过程。     

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

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

低速发动机短舱压强分布计算

应用第二代面元法对发动机短舱的压强分布进行了计算。计算结果表明其精确度与第一代高阶面元法相当;并可避免第一代面元法计算内流时出现的“泄漏”(leakage)问题。该方法具有节省机时和便于复杂外形应用的特点。     

时空守恒CE/SE数值方法在凹槽燃烧流动中的应用研究

应用时空守恒元和解元CE/SE方法,发展了求解化学反应流场的三维并行计算程序。为验证程序,对超声速冲压喷射燃烧室模型进行了计算,计算结果与实测数据以及相关文献的计算结果吻合较好,验证了本文的处理方法和程序的正确性。采用H2和O2的8组分34基元化学反应模型,对不同长深比的凹槽燃烧流动进行了数值模拟,分析了凹槽增强混合对稳定火焰的作用。计算结果表明,CE/SE方法具有模拟凹槽超声速燃烧流动的能力,数值模拟的结果有助于认识凹槽的火焰稳定机理。