超燃冲压发动机是临近空间高超声速飞行器的主要动力形式之一，因而对于超燃冲压发动机的研究一直是国内外关注的热点。超声速燃烧室是超燃冲压发动机的核心部件，用于组织碳氢燃料与空气混合和燃烧，对于超声速燃烧室内的燃烧过程，国内外研究者开展了很多的研究工作。但是受限于当下试验设备和测量手段的限制，很难开展燃烧室流场的高精度测量，尤其是对于这种大尺寸全流场、极短试验时间、强非定常效应的超声速燃烧室，很难获得燃烧流场的三维信息，数值计算方法逐渐成为研究超声速燃烧室湍流燃烧过程的有力工具。超声速燃烧室湍流燃烧的数值模拟需要解决化学反应动力学和湍流燃烧反应相互作用的问题，以往的研究工作主要采用总包反应模型或研究对象为机理比较简单的氢燃料，研究重点也放在燃烧室燃烧稳定后的释热规律等，较少对燃烧室从点火到火焰稳定整个非定常燃烧过程进行研究。因此有必要对超声速燃烧室内湍流燃烧的非定常过程进行数值研究，并考虑反应机理更复杂的碳氢燃料，研究碳氢燃料点火与火焰传播的非定常过程、燃烧与湍流的相互作用规律、燃烧室燃烧不稳定性等规律，为发动机的设计优化提供一定的参考。

本文对超声速燃烧室内碳氢燃料的湍流燃烧过程进行了研究，建立了基于简化反应机理的湍流燃烧非定常数值模拟方法，能够对燃料的点火、火焰传播等非定常现象进行数值研究。主要内容如下：

(1) 针对碳氢燃料复杂的化学反应机理，应用基于误差传递的直接关系图法和敏感性分析方法（DRGEPSA）对乙烯、正十二烷和煤油替代物详细反应机理进行了简化，获得了规模小、精度高的简化反应机理，并通过计算点火延迟时间和层流火焰速率对简化反应机理进行了验证。同时研究表明，机理简化存在下限，当组分数和基元反应数减少到一定数目时，进一步减少组分数和基元反应数会导致简化机理的精度降低。

(2) 采用分离涡模拟（DES）方法对乙烯/空气射流扩散火焰进行了数值计算，湍流/燃烧相互作用采用了涡耗散概念（EDC）模型，对乙烯简化反应机理的适用性进行了数值验证，并与直接数值模拟结果进行了比较，研究了湍流涡旋与火焰面相互作用规律。结果显示分离涡模拟结合乙烯简化反应机理能够有效模拟乙烯扩散燃烧过程，计算得到的温度剖面曲线与火焰推举高度均与直接数值模拟结果吻合得较好；同时，非定常计算结果展示了湍流涡旋对火焰面影响规律，涡旋可以通过拉伸、弯曲和褶皱等方式对火焰面产生影响。

(3) 基于德国宇航中心DLR支板超声速燃烧室对氢燃料扩散火焰进行了数值研究，计算采用了分离涡模拟方法和稳态火焰面模型，并与实验结果进行了对比，对火焰面模型在超声速燃烧数值模拟中的准确性进行检验。计算结果与实验结果吻合得较好，并且结果要明显优于RANS的计算结果；数值结果给出了火焰传播演化的过程，涡旋裹挟着燃料从支板底部脱落，形成离散的反应区“孤岛”结构；随着碳氢化合物组分（乙烯，甲烷）的加入，混合燃料燃烧强度降低，反应区变窄，流场的平均温度也比氢燃料降低约200K；燃烧效率也随之降低。

(4) 对乙烯在带凹腔超声速燃烧室内的燃烧非定常过程进行了数值研究，为了弥补稳态火焰面的局限性，采用了基于火焰面生成流形（FGM）方法的部分预混火焰面模型，利用大涡模拟方法对燃烧室湍流流动进行模拟，并与实验结果进行了对比。为了模拟真实的燃烧室火花塞点火过程，计算加入了火花模型（Spark Model）在冷态流场稳定后开启。冷态流场计算的沿程压力分布曲线、纹影图与实验结果基本一致；火花模型开启后，燃料开始点火，计算得到了与实验拍摄的CH*发光图像定性一致的点火过程图像，燃烧稳定后的沿程压力分布也与实验吻合得较好。计算还展示了中心区和近壁区点火过程的显著差异，中心区火焰结构基本稳定，而壁面火焰结构演化比较复杂，出现了火焰沿边界层逆流传播的现象；随着燃烧趋于稳定，近壁区的这种燃烧现象逐渐消失；释热导致的前传激波会打碎相对有序的湍流涡结构，而燃烧释热会耗散掉小尺度涡旋结构，使流场呈现出大尺度结构。对燃烧室压力和组分浓度脉动进行了提取和分析，发现燃烧室存在较低频率的振荡，分析了这种频率出现的机制，认为是由燃烧室在亚声速燃烧区形成的声学-对流反馈循环导致的热声耦合振荡造成。加入氢、甲烷等燃料后，会使得燃烧的强度变弱，火焰前传距离也变短，同时火焰形态也变得更窄。