用准稳态方法建立碳氢燃料点火燃烧的简化化学反应动力学模型

基于“准稳态”方法建立了一套碳氢燃料点火燃烧的化学反应动力学模型简化方法和相应的软件SPARCK,并从甲烷点火燃烧的GRI2.11详细基元反应动力学模型出发简化得出了包含16个组分12步总包反应形式的简化化学反应动力学模型,从庚烷点火燃烧的详细基元反应动力学模型出发简化得出了包含25个组分21步总包反应形式的简化化学反应动力学模型。通过其计算结果与CARM软件导出简化模型的计算结果和典型激波管试验结果的对比可以看出,本文简化得到的简化反应动力学模型能较为有效地再现详细基元反应模型的反应机理,简化模型的计算精度与CARM软件导出简化模型的计算精度相当。与详细基元反应动力学模型相比,简化模型有效地减少了反应组分,在工程计算中有比较好的应用前景。     

Ap复合固体推进剂气相点火模型

本文根据异质推进剂的特点,提出了含氧流动热气体点燃Ap复合推进剂的一维气相点火模型。模型中详细考察了点火过程中推进剂表面的分解过程和气相区的化学动力学过程,并利用有限差分法直接求解点火过程的控制方程,获得了点火延迟时间t_(ig)随燃烧室压力P变化的关系,固相区和气相区的温度分布,以及参加反应的各种化学组分在气相区的分布。t_(ig)随P变化的理论计算结果与实验测定曲线比较接近。对于深入研究点火问题有一定参考价值。利用本模型和计算方法还可以从理论上预示其它各种参数对Ap推进剂点火延迟的影响。     

乙烯简化化学动力学模型在HIFiRE燃烧室计算中的应用

为评估化学反应动力学模型对超燃冲压发动机计算结果的影响,首先利用基于"准稳态"假设方法建立了复杂化学反应动力学模型简化软件包(SPARCK),从详细模型出发得到了一个包含20组分16步总包反应的乙烯简化动力学模型.然后采用该简化模型和Princeton大学的简化模型对乙烯点火延迟时间进行了计算,对比显示两个简化模型计算结果基本一致,与实验测量结果均吻合较好,都能准确反映乙烯点火特性.但对HIFiRE燃烧室直连式实验的数值模拟结果显示两个模型计算得到的燃烧流场静温分布差别较大,内流道推力性能相差近12%.相比于Princeton大学的简化模型,SPARCK软件得到的简化模型计算的壁面压力分布与实验结果吻合更好,能够准确描述燃烧室流场现象.     

激波管煤油燃料点火特性研究进展

为了研究煤油的点火特性,国内外采用激波管实验方法,探究了不同因素对煤油的点火延迟时间的影响。本文综述了激波管煤油燃料点火特性研究进展,介绍了激波管的设计与实验要点以及液、气相煤油在激波管中点火特性的研究。结果表明,温度、压力、当量比等因素对煤油的点火特性有着复杂的影响,对目前研究中存在的问题与不足进行了归纳总结,对未来研究的重点进行展望。     

复合推进剂气相点火过程的分析与计算

根据异质推进剂的特点,提出了含氧流动热气体点燃PBAA复合推进剂的气相点火模型。为了详细了解点火过程,进而准确预示点火延迟t-(ig),模型中详细考察了点火期间推进剂表面的分解过程,并利用有限差分法直接求解点火过程的控制方程,获得了点火表面附近的温度和浓度分布与变化曲线。通过与已有实验结果的比较,结果表明,文中提出的理论模型基本合理,计算方法也实用可靠,它们对于深入研究点火问题有一定参考价值。     

爆震波多管点火特性实验

为了研究液体火箭发动机爆震波多管点火的同步性能、多次点火重复性能以及点火火炬性能,组建了氢氧爆震波多管点火实验系统。采用氢气和氧气为工质,常温供气压力0.10.5 M Pa（表压）,混合比1.87.2,进行了多次实验。实验结果表明:爆震波点火技术可以在与液体火箭发动机贮箱压力相适应的较低的供气压力下获得高温（>1300℃）高压（>1 M Pa）爆震产物,并且具备良好的点火重复性能和多管点火同步性能,多次点火重复性时间差和多管点火同步性时间差均小于0.3 m s。爆震波多管点火技术适合用于多燃烧室液体火箭发动机的同步点火。      

高能点火药点燃性能的研究

以 Co₂ 激光为能源对 Al+KClO₄ 高能点火药的点火特性进行了研究。结果表明:Al和KClO₄ 的质量比各为 50 %时,点火延迟时间最小。点火药中掺入少量硫粉能大大减少点火延迟时间。点火延迟时间随点火药的初温的升高而减小。点火延迟时间随点火热流密度的增大而减小。存在一个最小点火热流密度,小于此值点火药将不能被点着。在实验观察和理论分析的基础上,提出一个Al+KClO₄ 高能点火药的点火模型,并对由此模型而导出的控制方程进行了数值求解,计算结果与实验值相符。      

两相爆震波多管同步点火技术的实验研究

为了研究基于液态燃料的爆震波点火技术,进行了一系列液态煤油/氧气爆震波点火实验。实验中氧气和煤油的供给压力分别为1.0 MPa和0.7 MPa,火花塞点火能量为50 mJ。研究了两相爆震波点火技术的基本特性,实验表明:采用低点火能量能够快速产生充分发展的爆震波,煤油/氧气爆震压力可达4.0 MPa,爆震波速度可达1500 m/s到2001 m/s,尾焰温度约为2075 K。开展了两相爆震波由单管向多管传播的实验,验证了两相爆震波多管点火技术的可行性,目前可实现四管同步点火。实验显示两相爆震波点火技术重复性强,多管点火具有较好的同步性,时间差别为几十个微秒量级,适用于多燃烧室液态火箭发动机的同步点火。     

C02激光点燃固体推进剂的研究

选用具有代表性的复合型、双基型推进剂以及金属粉组成的点火药,在新建的CO_2激光点火试验设备上进行了推进剂点火性能试验,获得了重复性较好的试验数据。针对复合推进剂提出一个点火模型及点火判据表示式,计算与试验结果相吻合,证明该模型可以用于预估复合推进剂的点火性能。试验和计算结果都表明点火热流密度对点火延迟时间影响很大,是推进剂点火性能中非常重要的控制量     

弯曲方管内爆燃向爆震转变特性

为缩短爆震管轴向尺寸,采用曲管爆震管替代直管爆震管.基于弯曲方管实验器,进行了二维数值仿真计算,研究了点火位置、点火能量及弯曲段曲率半径对爆震波形成及传播的影响规律.数值模拟结果表明:爆燃向爆震转变(DDT)时间随弯曲段曲率半径增大而减小,DDT距离在弯曲段曲率半径最大时达到最小值;点火能量越大,爆震管内激波强度越大,DDT时间越小;在距离爆震管进口300mm处点火时,爆震管内爆震波压力峰值最大,DDT时间最短.     

强迫对流下硼颗粒燃烧特性影响因素研究

针对冲压发动机燃烧室内强迫对流下的硼颗粒燃烧特性展开了系统研究,考虑气相流动、扩散和表面单步有限化学反应动力作用,建立了强迫对流下硼颗粒燃烧过程的物理和数学模型;采用有限体积法求解含多组分反应流的二维轴对称Navier-Stokes方程,并验证了数值仿真方法的正确性。首先通过数值仿真研究了来流速度、颗粒半径、环境中氧气质量分数和环境压力等因素对硼颗粒燃烧特性的影响,并对其成因展开了详细分析。研究表明,在强迫对流作用下,硼颗粒总的燃烧质量流率和质量流率通量均随来流速度、颗粒半径、环境中氧气质量分数和环境压力的增加而增大。通过深入分析发现,强迫对流下硼颗粒的燃烧质量流率通量随着来流雷诺数的增加而增大。然后基于大量数值仿真结果,对相对静止气氛下的硼颗粒质量流率通量进行了修正,用于描述强迫对流下的硼颗粒燃烧特性。     

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

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

使用GAO—YONG湍流模型数值研究管道凸起流动

用GAO—YONG可压缩湍流方程数值模拟了Delery管道凸起跨音流场中的激波湍流边界层干扰现象。分析了GAO-YONG可压缩湍流方程组对湍流的非平衡、多尺度、各向异性等特性的描述能力。计算中对流项、扩散项分别采用二阶ROE格式和二阶中心差分格式离散,并用多步Runge-Kutta显式时间推进法求解了空间离散后的控制方程。计算很好地模拟到了压力平台区、“入”波结构等典型激波湍流边界层干扰的流动现象,也得到了壁面压力分布、平均速度剖面以及雷诺应力分布等,并与相应的实验数据进行了对比分析,两者符合很好。     

气气喷注器推力室传热特性研究

设计了单喷嘴气气喷注器容热式推力室,进行了0.92～6.1 MPa范围内7个燃烧室压力工况,共17次热试车;采用壁面测温方法获得推力室沿轴向方向的内壁热流分布,得到各室压工况下燃烧室内壁热流分布曲线,不同工况的热流曲线显示出相似的分布;并从中分离出对流传热热流,得到对流传热热流与室压的关系。为拓展范围,采用多组分湍流N-S方程描述推力室内燃烧流动,采用6组分9步反应模型来描述氢氧反应,反应速率由Arrhenius公式计算,进行了5～20 MPa更高室压范围内的燃烧内流场的数值模拟,并耦合计算了各工况燃气与室壁之间的传热,获得了与试验相同规律的结果。     

模型燃气轮机燃烧室中三维气相反应流的数值计算

 本文给出一种模型燃气轮机燃烧室中三维气相反应流流场的数值计算结果。计算中采用了著名的SIMPLE算法,改进的κ-ε双方程湍流模型及湍流扩散燃烧模型。研究了燃烧过程对流动图形的影响。计算结果经过实验验证,两者符合较好。     

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

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

有限元法求解QPNS方程模拟高速旋成体气动力和热

从ＦＮＳ方程中略去流向粘性导数项导出熵变量形式ＱＰＮＳ方程（拟抛物化Ｎ－Ｓ方程）其离散解能自动满足熵增不等式。利用ＧＬＳ（伽辽金／最小二乘法）有限元法构造出ＱＰＮＳ方程的弱解式,所形成的整体方程组用ＧＭＲＥＳ（广义最小余数法）求解。采用空间推进求解ＱＰＮＳ方程,能明显降低计算机内存和运算时间。算例给出典型旋成体的物面压强、剪应力、热流量分布以及流场中密度、温度分布。     

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

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

大速差射流燃烧室内三维湍流气固两相流的实验研究

本文用二维激光多普勒测速仪和探针取样法对大速差射流燃烧室中冷态三维湍流气固两相流动进行了实验研究,测量了气相速度,颗粒速度、颗粒质量流分布及两相各自的湍流特性,探讨了高速射流作用下两相流场的主要物理特征及其对火焰的稳定和强化燃烧的作用。     

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

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