多孔推进剂装药的动态压缩过程的数值研究

对约束条件下的推进剂多孔床中的动态压缩进行了数值分析。用两组连续混合非平衡流的理论建立了描述动态压缩过程的数学模型。用MacCormack有限差分方法求解了控制方程。以活塞撞击HMX多孔床为例，将数值结果、实验结果以及稳态分析结果进行了比较，证明模型能很好的预测稳态压缩波特性。研究结果有助于弄清高能推进剂的冲击波起爆、燃烧转爆轰的机理。     

高能固体推进剂燃烧转爆轰（DDT）研究综述

以固体火箭发动机装药安全性为背景，从高能固体推进剂的燃烧转爆轰（简称ＤＤＴ）危险性、燃烧转爆轰机理和燃烧转爆轰的模拟三方面对该问题的研究工作进行了综述，重点讨论了ＤＤＴ过程中起关键作用的动态压缩现象、压缩燃烧和热点形成机理，以及ＤＤＴ建模中的本构关系和封闭问题，指出了今后进一步研究的方向。     

发射药在强约束条件下的燃烧转爆轰特性

以Ｄ－４５，Ｓ－１４，３／１三种发射药为研究对象，对钢管约束条件下的高密实火药颗粒床的燃烧转爆轰机理和敏感性进行了实验研究。分别用变片和离子地测量了药床中的压缩波和燃烧波的传播，通过对压缩波阵面和燃烧波阵面的分，发现压缩波形成火焰峰后的已燃区，它将不断加速与加强，最后赶上火焰峰形成强冲击波而诱发爆轰。     

多孔床复合固体推进剂燃烧转爆轰实验研究

详细研究了高限制条件下,多孔床丁羧和丁羟推进剂的燃烧转爆轰过程,考察了推进剂的装药密度及其颗粒尺寸的大小对装药燃烧转爆轰的影响。实验结果表明:多孔床丁羧和丁羟推进剂在一定的条件下是可以产生燃烧转爆轰的,而且构成多孔床推进剂的碎片越细小,则推进剂越易产生燃烧向爆轰的转变。     

脉冲爆轰发动机内三维两相爆轰的数值计算

根据脉冲爆轰发动机内脉冲爆轰的特点,建立三维两相爆轰的理论模型。应用国际上最近发展起来的CE/SE方法,数值模拟爆轰管内三维两相燃烧转爆轰的过程。计算结果表明:当点火起爆后,在燃烧转爆轰过程中,三维效应非常明显。当形成稳定爆轰波后,三维效应不明显。计算压力值和实验压力值符合得较好。     

固体推进剂燃烧转爆轰模拟计算研究

采用一维分离两相流反应模型以及变步长跳步差分计算格式,模拟了固体推进剂的燃烧转爆轰过程。计算结果详细、合理地描述了推进剂在燃烧转爆轰过程中的流场分布情况。文中建立了确定诱导爆轰距离(L_(DDT)的方法,由该法所得L_(DDT)的值与实验值符合得较好,其误差小于15％。     

与爆轰相关的湍流燃烧

爆轰具有较高的燃烧效率,当应用于推进系统时,有较宽广的前景,故引起国内外学者广泛的关注。与爆轰有关的两个重要现象是——爆轰的生成和传播。爆轰的生成分为直接起爆和间接起爆。直接起爆需要极高的能量,因此间接起爆一直是研究的热点,其中爆燃转爆轰(deflagration-to-detonation transition,DDT)则是关注的重点。然而,经过几十年的实验、理论和数值研究,DDT的机理已被理解,但是准确的预测还无法实现,这其中最核心的因素是湍流燃烧。对于爆轰波的传播,特别是高活化能的可燃气体,其波阵面呈强不稳定结构,存在湍流燃烧,波后出现的未燃气团的燃烧也属于湍流燃烧,这使得爆轰波的精确控制也存在较大的困难,但是爆轰波传播的平均速度仍然可以用CJ速度来预测。本文对DDT和爆轰波结构中的湍流燃烧的研究现状进行了分析和评述,并对未来的研究进行了展望。     

曲面激波诱导斜爆轰的数值模拟

为了研究恰当比预混氢气-空气斜爆轰流场的波系结构和流动特征，基于带化学反应的Navier-Stokes方程，对弹头及楔穿越预混气体时诱导的斜爆轰进行了数值模拟。对流项的离散采用Steger-Warming格式，时间项采用二阶Runge-kutta方法。结果表明，对于弹头：（1）在亚爆轰条件下，能够模拟氢气-空气预混诱导爆轰流场的精细结构；（2）在超爆轰条件下，通过精细调整网格，能够很好地分辨强烈耦合的激波和燃烧波，且与Lehr实验吻合良好。对于楔结构：捕获到了清晰的三波点及其复杂精细的斜爆轰流场结构，预测的诱导燃烧距离、激波角和斜爆轰角与实验吻合良好；通过对流场波系结构变化过程的研究，获得了流场三波点随时间的演化过程。     

斜爆轰波系在受限空间内的演变及其临界条件的数值研究

作为一种面向未来高马赫数吸气推进技术的新概念发动机,斜爆轰发动机的运行依赖其稳定的宏观波系结构。本文针对斜爆轰发动机简化模型,采用多组分基元化学反应流动数值模拟技术,数值分析斜爆轰波系在受限空间内的宏观结构特征及其演变,并进一步分析爆轰波系结构转变的临界条件。研究结果发现,随着楔面角度的增加,依次出现四种结构：激波诱导燃烧,斜爆轰波双规则反射,回流区马赫反射,楔面燃烧。对于稳定的波系结构,楔面压缩角同时存在下临界、亚临界以及超临界三种极限条件,在波系演变过程还伴随激波规则反射和马赫反射的转变。     

进气温度和当量比对脉冲爆轰发动机工作过程影响

为了研究进气温度和当量比对脉冲爆轰发动机工作过程的影响,建立了带简单化学反应的气液两相爆轰欧拉-拉格朗日模型,使用二维时空守恒元与求解元(CE/SE)方法和变步长4阶龙格-库塔法分别求解气液两相爆轰方程.计算结果表明:提高进气温度,能加速液滴雾化、蒸发,缩短燃烧转爆轰距离和时间,但是会降低爆轰波的峰值压力;当量比小于1.1时,增加当量比,能缩短燃烧转爆轰距离和时间,提高爆轰波的峰值压力,加快爆轰波传播速度;当量比大于等于1.1时,增加当量比,能小幅提高爆轰波速传播度和缩短燃烧转爆轰距离和时间,但对爆轰波的峰值压力影响很小.      

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

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

球形点火腔对脉冲爆震DDT过程影响的数值研究

为了减小脉冲爆震发动机PDE（pulse detonation engine）单次循环时间，缩短爆燃向爆震转变DDT（deflagration to detonation transition）距离是关键.为此，提出了一种在头部添加球形点火腔的新型爆震管道结构设计方案，并采用丙烷和空气为可爆混合物，通过对新型爆震管DDT过程的二维数值模拟，研究其对DDT距离以及DDT时间的影响.数值模拟结果表明，头部采用球形点火腔后，点火腔中的压缩波经过多次反射后，能够在爆震管中更快地促使爆震波形成，DDT距离大大缩短；当球形点火腔直径为1.5倍等直爆震管直径时，相对于常规等直管爆震管结构，其DDT距离和过程时间分别减少了14%和16.26%.      

环形爆震室中火焰加速数值模拟及验证

对带不同数量孔板的环形爆震室进行数值模拟,并通过试验对数值计算进行验证,来研究火焰加速现象、爆燃向爆震转变过程和不同当量比下起爆距离.数值计算采用二维轴对称非定常Navier-Stokes方程来模拟流体动力学过程.研究发现用较低的点火能量对火焰混合区的可燃气点火产生低速火焰,低速火焰向环形爆震室射流并改变方向向出口传播,火焰在孔板的阻碍作用以及火焰诱导激波和反射波的加速作用下,由层流变为湍流,湍流火焰与其诱导激波相互加强,最终引爆未燃混气;还对爆震波在孔板区的传播过程进行了分析,对不同当量比下的火焰速度和起爆距离进行了模拟研究.      

汽油/空气两相脉冲爆震发动机触发爆震的研究

有效的触发爆震波是脉冲爆震发动机正常工作的关键所在。为了在多循环两相脉冲爆震发动机的缓燃到爆震转变过程中,有效的控制激波反射,成功的触发爆震。通过在爆震室内合理的排放障碍物来促进激波与火焰的相互作用,组织激波反射,适时的触发爆震,从而有效地缩短促发爆震的距离,并最终获得充分发展的爆震波。试验发现爆震管内同样的障碍物在不同的位置起到的作用不同。在直径58 mm,长度1275 mm的汽油/空气PDE爆震管内成功触发频率为50 Hz的爆震。      

激波聚焦诱导气液两相爆震燃烧的数值模拟

对以激波聚焦和增加障碍物方式诱导煤油-空气气液两相爆震燃烧的过程进行了数值模拟.采用欧拉-拉格朗日方法建立了脉冲爆震发动机(PDE)中气液两相流的喷射、雾化、掺混过程.研究发现环形爆震波在爆震管凹腔内经过反射、汇聚后能够引燃可燃混合物.而在障碍物处,激波的反射和再反射聚焦能够形成高温高压点(2700K,25MPa),产生局部爆炸,有助于形成稳定的脉冲爆震燃烧(波面速度为1900m/s,温度为2 950K),有效地缩短由缓燃向爆震转变(DDT)距离至0.45m.     

间接起爆对脉冲爆震发动机热效率影响

建立了基于整个燃烧室内经历不同燃烧方式下所有工质来评估脉冲爆发动机平均循环热效率的方法,通过二维数值方法对两种采用不同数目障碍物的爆震管模型的起爆特性及循环热效率进行了研究,并分析了燃烧波传播的距离与DDT距离之比（定义为ξ）对模型热效率的影响.结果表明:①间接起爆方式及不同DDT强化装置都会影响系统的循环热效率;②就计算模型,当ξ接近1时,其系统平均循环热效率仅为理想爆震循环热效率的60%左右;③随着ξ的增加,间接起爆方式的影响迅速减弱;④当ξ为1.5时,系统平均循环热效率可达爆震燃烧热效率的93%.      

热射流点火对爆震管内火焰加速及爆震波触发影响的数值研究

为研究热射流点火对爆震管内火焰加速及爆震波触发的影响,运用34步26组分丙烷基元反应进行了二维数值模拟,获得了5种不同的热射流发生器几何结构下的爆震管内火焰传播规律及缓燃向爆震转捩(deflagration to detonation transition,DDT)的时间与距离.结果表明:爆震管中湍流在火焰加速阶段起着重要作用,在爆震波触发阶段激波与火焰相互作用占有主导地位.根据DDT过程的定义,得到DDT时间在1.4~2.0ms之间,同时发现热射流发生器长度为150mm,热射流发生器孔径在8mm时DDT时间最短,热射流发生器长度及孔径对DDT距离的影响不大.      

脉冲爆震发动机点火过程离子催化效应数值模拟

利用气体电离理论推导出氢气-空气混合气体电离后组成成分,理论分析活性基团对燃烧速率及剧烈程度的催化效应,以及不同点火能量、活性基团浓度对缓燃转爆震(DDT)过程的影响.结合氢气-空气燃烧23步化学反应动力学机制,采用FLUENT软件对不同工况下的DDT过程进行模拟,与理论分析结果对比.结果表明:点火温度为2000~2500K时,活性基团的加入,可提高燃烧速率,DDT时间可缩短9.91%~21.08%,DDT距离可缩短3.32%~8.08%,DDT时间和DDT距离的改变幅度随点火温度的升高而增大.点火能量较高时应该考虑气体电离效应.