CE/SE方法在脉冲爆轰发动机气液两相爆轰流场计算中应用的研究

摘要

脉冲爆轰发动机是一种利用脉冲式爆轰波产生推力的新概念发动机。与传统的以等压方式燃烧的发动机相比，它具有结构简单、热循环效率高、单位燃料消耗率低、工作范围广等优点，这使得它在军用和民用方面具有广阔的应用前景，将成为未来重要的新型推进系统。本文以液体燃料脉冲爆轰发动机为研究对象，针对两相爆轰流场的特点，建立两相爆轰控制方程，应用近期发展起来的时空守恒元和求解元方法(简称CE/SE方法)数值求解，用Fortran语言编写程序模拟两相爆轰过程。根据时空守恒元与求解元思想，推导了一维、二维和三维CE/SE方法的计算格式，用于对脉冲爆轰发动机管内气液两相爆轰一维、二维、三维流场的数值计算。主要内容包括：
　　 (1)推导了一维CE/SE方法的计算格式，用于脉冲爆轰发动机内汽油/空气两相爆轰流场的一维数值计算，编写了相关计算程序，分析并讨论了一维两相爆轰波传播过程中物理量的变化规律。数值实践表明：CE/SE方法可高效捕获爆轰波等强间断。
　　 (2)推导了二维无粘CE/SE方法的计算格式，编写了相关计算程序，实现了对脉冲爆轰发动机内两相流场的数值模拟，获得了脉冲爆轰发动机内流场的变化规律，分析并讨论了液滴尺寸分布以及不同的起爆温度和起爆压力对两相燃烧转爆轰过程的影响。计算结果表明：起爆温度和起爆压力越大，汽油雾化后液滴尺寸越小，越有利于实现燃烧转爆轰过程。
　　 (3)推导了二维粘性和湍流CE/SE方法的计算格式，分别用于脉冲爆轰发动机内两相粘性和湍流爆轰流场的数值模拟，编写了相关计算程序，得到了脉冲爆轰发动机两相粘性流场和湍流流场的变化规律，分析并讨论了扰流片片数和间距对两相燃烧转爆轰过程的影响。研究发现：相对于无粘模型，应用粘性或湍流模型来描述两相爆轰过程，可更加精细刻画流场特性，更加符合实际过程。计算结果表明：扰流片可有效促进燃烧向爆轰转化的进程。
　　 (4)推导了三维无粘CE/SE方法的计算格式，用于脉冲爆轰发动机内两相流场的数值计算，编写了相关计算程序，分析并讨论了三维无粘两相爆轰流场的变化规律。研究发现，在燃烧转爆轰过程中，存在明显的三维效应。本文详细刻画了两相燃烧转爆轰过程中爆轰波三维结构特征。
　　 (5)推导了三维粘性CE/SE方法计算格式，用于脉冲爆轰发动机内两相流场的数值计算，编写了相关计算程序，分析并讨论了三维两相粘性爆轰流场的变化规律。研究发现：三维粘性数值计算结果比三维无粘数值计算结果更真实地反映爆轰流场特性，尤其在近壁附近，可更加精细刻画爆轰波的三维结构特征。
　　 (6)进行了以汽油/空气为工质的吸气式脉冲爆轰发动机实验研究，实现了33Hz、40Hz、50Hz三种不同工况下的爆轰，得到了三种工况下的爆轰波压力曲线，分析了进气压力变化对脉冲爆轰发动机工作性能的影响，为数值计算提供了实验验证。