预燃级旋流数对双环腔燃烧室流场的影响.pdf

预燃级旋流数对双环腔燃烧室流场的影响李 茂，钟世林，周 雄，马存祥，李逸飞(中国航发四川燃气涡轮研究院，成都 610500)摘 要：采用数值仿真方法，研究了双环腔燃烧室预燃级不同旋流数下的冷态流场和油雾场。结果表明：在旋流器出口会形成 2 个旋涡，在离旋流器出口较远处受主燃孔射流和中心隔离环射流影响也会形成 1 个旋涡。随着旋流数增大，流量系数减小，气流的回流速度及回流量增大，燃油雾化质量改善，从而改善点火性能和小工况燃烧性能；回流区长度基本一致，但回流区径向尺寸增大，旋流器出口燃油质量分数增大，但离旋流器出口 60 mm 处的燃油质量分数减小。受冷态流场 3 涡结构影响，存在 1 个最佳旋流数使得点火电嘴位置处油气比最佳。关键词：航空发动机；双环腔燃烧室；旋流器；旋流数；油雾场；数值仿真；点火中图分类号：V231.2 文献标识码：A 文章编号：1672-2620(2022)06-0000-08Effect of swirl number of pilot swirl device on the flow field in recirculation zone of double-annular combustorLI Mao,ZHONG Shi-lin,ZHOU Xiong,MA Cun-xiang,LI Yi-fei(AECC Sichuan Gas Turbine Establishment,Chengdu 610500,China)Abstract:Numerical simulation was used to study the flow field and fuel mixed field of pilot swirler in the double-annular combustor with different swirl numbers.The results show that two vortices will be formed at the outlet swirler,and a vortex will also be formed at a distance from the outlet wirler due to the influence of the main combustion hole jet and the central isolation ring jet.With the increase of swirl number,the flow coefficient decreases,the axial velocity and return flow of airflow increase which can improve the atomization of fuel and improve the ignition performance and combustion performance under small operating conditions.The length of the recirculation zone under different swirl numbers is basically the same.With the increase of swirl number,the size of the recirculation zone in the radial direction increases and the fuel mass fraction at the outlet swirler increases,but the fuel mass fraction at the outlet swirler 60 mm decreases.Due to the influence of the three vortex structure in the flow field,there is an optimal swirl number that maximizes the fuel-air ratio at the ignition nozzle position.Key words:aero-engine;double-annular combustor;swirler;swirl number;fuel mixed field;numerical simulation;ignition1 引言 双环腔燃烧室在民用航空发动机中应用非常广泛，其主要应用目的是降低污染排放，如GE90、E3 和 CFM56-5B 等发动机中均使用了双环腔燃烧室1-3。双环腔燃烧室由 1 个中心体将火焰筒分隔成内、外 2 个并行的燃烧区，外侧的为预燃级，内侧的为主燃级，从而实现分级分区燃烧。小工况下，仅预燃级喷嘴喷油工作，保持较高的油气比，气流速度低，燃烧时间长，燃烧完全；大工况下，预燃级和主燃级的喷嘴共同工作。对于最大工作马赫数 3 及以上、大体积流量进口、宽工作范围的发动机而言，燃烧室采用径收稿日期：2021-06-08作者简介：李 茂(1996-)，男，四川泸州人，助理工程师，硕士，从事燃烧室性能设计。第 35 卷 第 6 期2022 年 12 月燃气涡轮试验与研究Gas Turbine Experiment and ResearchVol.35,No.6Dec.,2022 12第35卷 燃 气 涡 轮 试 验 与 研 究 13 向分级的双环腔构型，让不同燃烧区域发挥不同的功能，能够满足发动机大体积流量进口、宽工作包线下工作的要求。由于发动机工作包线范围宽，要求风车点火高度在 10 km 以上，且有逆向模态转换过程中能够成功点火的要求，这对燃烧室的点火性能提出了较大挑战。从点火的几个过程中容易发现，燃烧区空间的油气比、燃油分布及速度分布是影响点火性能的关键因素4-5，但速度分布、油气比、回流区结构、燃油雾化质量等又主要受到头部旋流器的影响。为了获得良好的燃油雾化特性，形成良好的流场结构，国内外学者对旋流器做了很多研究，包括旋流器压降、文氏管结构、旋流器几何形状等6-7。PanduRanga Reddy 等8采 用 PIV(粒 子图像测速仪)，研究了单级旋流器矩形燃烧室内的旋流流场，发现随着旋流器气流压降的减少，旋流器后回流区尺寸相应减小；Cai 等9研究了轴向旋流器模型燃烧室油雾燃烧流场，其试验结果表明，油雾速度分布对旋流流场起重要作用，而文氏管结构对油雾速度分布影响较大；蔡文祥等10采用数值模拟方法对旋流器进行了研究，表明旋流器几何形状对燃烧室流场影响较大，优化旋流器设计可提高燃烧性能。但这些研究大多针对单环腔燃烧室，对于工作在大空间、宽工作范围的双环腔燃烧室，对其流场、雾化场、预燃级点火性能及小工况下的燃烧性能研究较少。本文在单级轴向旋流器的基础上，采用数值仿真方法，研究了轴向旋流器旋流数对双环腔燃烧室流场和油雾场的影响，可为双环腔燃烧室的稳定点火及小工况燃烧性能改善提供一定的理论支持。2 数值仿真2.1 物理模型 以某型双环腔燃烧室为研究对象。燃烧室外侧为预燃级，采用单级轴向旋流器；内侧为主燃级，采用双级径向旋流器。本文研究的旋流器均采用直叶片。预燃级头部与主燃级头部数目比为 2:1。双环腔环形燃烧室几何模型如图 1 所示，模型采用 2个预燃级头部和 1 个主燃级头部用于研究。由于本文不考虑温度分布，只计算冷态流场，所以除了保留主燃孔、掺混孔等较大的孔外，去掉了内外火焰筒壁上细小的气膜冷却孔。主燃级头部预燃级头部图 1 双环腔燃烧室及其计算流体域Fig.1 Double-annular combustor and computational fluid domain轴流式旋流器旋流数的计算公式11为：(1)式中：为旋流器旋流数，为旋流器叶片安装角，、分别为旋流器叶片径向位置的内半径和外半径。研究了 4 中旋流器方案：方案 1，采用单级轴向旋流器基准结构(图 2)，即基准旋流器方案；方案 24，以方案 1 为基础，改变旋流器旋流数，并保证文氏管喉道内径、文氏管扩张角和各旋流器与基准旋流器有效流通面积基本不变。旋流数改变主要是通过改变旋流器叶片安装角来实现。各旋流器方案几何参数对比见表 1。其中，ACd 为旋流器有效流通面积，N 为叶片数。与方案 1 相比，方案 2、3 主要是增加旋流器叶片外径，方案 4 是减少叶片数。各方案旋流器的轴向长度、喉道长度及半径不变。各方案叶片结构对比如图 3 所示。图 2 单级轴向旋流器基准结构及其中心剖面图Fig.2 Axial swirler structure and its central section view表 1 各旋流器方案几何参数对比Table 1 Comparison of geometric parameters for various swirler schemes方案SN/()ACd/mm2NSN/%10.60938.0325.89200.0020.79646.0326.182030.7031.00453.0324.332064.8641.24058.6325.4212103.602.2 网格划分 由于燃烧室模型结构较为复杂，因此采用非结构化网格对流体域模型进行网格划分，并在冷却孔、14 李 茂等：预燃级旋流数对双环腔燃烧室流场的影响 第 6 期 由于试验装置有几个速度进口，且在无限制的空间内进行燃烧，因此需对其进行简化。简化后采用二维仿真模型进行仿真，其流体域宽为 130 mm，高为 130 mm，长为 400 mm，以保证其燃烧不受空间限制。采用结构化网格，其网格划分结果如图 6 所示。采用标准和 Realizable 两种湍流模型和增强壁面处理函数进行数值仿真。选取轴向位置X=6.8，10.0，20.0，40.0，70.0，100.0 mm 处沿径向的轴向速度进行对比，结果见图 7。相比于标准湍流模型，除了X=40.0 mm入了与旋转和曲率有关的内容，且方程中不再包含，能够更好地表示光谱的能量转换，即使 k值很小或等于 0 时也不会使得分母为 0，使方程具有奇异性。2.4 计算模型验证 为对比验证两种湍流模型精度，采用 Sydney 大学旋流装置(图 5)12-13实验数据与两种湍流模型计算结果进行对比。选取 N29S054弱旋流模型进行验证。旋流数设定为ss/SNUW=，表 2 给出了 N29S054 弱旋流模型的边界条件。(a)方案 1(c)方案 3(b)方案 2(d)方案 4图 3 各旋流器方案叶片结构对比Fig.3 Comparison of blade structure for each swirler case图 4 燃烧室三维网格划分示意图Fig.4 Schematic of three-dimensional meshing of combustor主燃孔及掺混孔附近进行局部加密。网格数量约为 2 200 万；最小网格质量为 0.14，最大扭曲度为 0.93，满足数值仿真要求。燃烧室三维网格划分见图 4。2.3 计算方法 数值仿真采用商业 FLUENT 软件，选取标准模型和 Realizable 模型进行湍流模型对比。Realizable 模型的湍流动能方程为：(2)(3)式中：表示由平均速度梯度引起的湍流动能 k 的产生项，表示浮力引起的湍流动能 k 的产生项，代表可压湍流中脉动扩张对总的耗散率的影响。与标准模型相比，Realizable 模型引表 2 N29S054 边界条件Table 2 N29S054 boundary conditionsUj/(m/s)Us/(m/s)Ue/(m/s)SNTair/KRejetRes66.0029.7420.000.5427315 70059 000305mm3x Tangential air streamHoneycomb meshesFuel streamAxialair streamWs60mm50mmUeUeUsUjUs15图 5 燃烧装置示意图Fig.5 Schematic of combustion device湍流耗散率方程为：第 35 卷 燃 气 涡 轮 试 验 与 研 究 15 截面处，其他截面处采用 Realizable 湍流模型的轴向速度与实验结果的吻合程度更好。另外，经大量验证，Realizable 湍流模型可有效应用在不同类型流动的数值仿真中，包括旋转均匀剪切流、管道内流动、边界层流动、包含有射流和混合流的自由流动，以及带有分离的流动等14。因此下文的数值计算中，计算模型采用 Realizable 湍流模型，壁面函数采用增强壁面处理函数，用 DPM模型模拟燃油分布，燃油喷点布置在离喷嘴出口0.5 mm 处，避免边界层网格的影响。图 6 二维简化模型结构网格示意图Fig.6 Schematic of two-dimensional simplified model structure grid0.00 0.02 0.04806040200轴向速度/(m/s)径向距离/m(a)X=6.8 mmRealizableEXPStandard0.00 0.02 0.04806040200轴向速度/(m/s)径向距离/m(b)X=10.0 mmRealizableEXPStandard0.00 0.02 0.046040200轴向速度/(m/s)径向距离/m(d)X=40.0 mmRealizableEXPStandard 0.00 0.02 0.046040200轴向速度/(m/s)径向距离/m(e)X=70.0 mmRealizableEXPStandard0.00 0.02 0.046040200轴向速度/(m/s)径向距离/m(f)X=100.0 mmRealizableEXPStandard0.00 0.02 0.04806040200轴向速度/(m/s)径向距离/m(c)X=20.0 mmRealizableEXPStandard图 7 不同轴向截面轴向速度对比图Fig.7 Comparison of axial velocity of different axial sections16 李 茂等：预燃级旋流数对双环腔燃烧室流场的影响 第 6 期3 计算结果与讨论3.1 网格无关性验证 采用基准方案，在参考压力 130 kPa、进口温度400 K 及燃烧室进口流量 0.676 4 kg/s 的点火工况下，对不同网格数量模型进行数值仿真。点火器横截面位置处的轴向速度如图 8(a)所示。可以看出，当网格数在 2 000 万以下时，随着网格数的增加，轴向速度变动较大；当网格数从 2 200 万增加到2 500 万时，2 条曲线基本吻合，轴向速度基本不再变化。不同网格数量下的有效流通面积如图8(b)所示。图中，有效流通面积通过式(4)计算。从图中看，模拟值要大于设计值，网格数量在 2 000 万到2 800 万之间计算出的有效流通面积与设计值较为接近。其原因在于，计算时忽略了下游的动压头，计算出的压降偏小，导致有效流通面积模拟值偏大。综上，对各个方案的数值模拟均采用 2 200 万左右的网格数。式中：表示实际空气流量；为压力降，即上游总压减下游静压；为进口空气密度。(4)3.2 旋流数对冷态流场的影响 针对轴向旋流器，设计时流量系数 Cd采用下式15计算：图 10 为 4 种旋流器结构下燃烧室旋流器中心截面的流线图。从图中可以看到，在各个旋流数下，旋流器出口(区域 1、2)都形成了 2 个旋涡，区域 3 形成 1 个旋涡。但旋流数为 0.609、0.796 时，旋流器出口为对称旋涡；旋流数为 1.004、1.204 时，旋流器出口为非对称旋涡。区域 3 旋涡的形成主要是由于中心隔离环射流与主燃孔射流的阻挡，使得气流速度减小、压力增大，形成回流区。回流区的存在能够形成稳定的点火源，有利于点火。中心隔离环射流的图 9 不同旋流数旋流器的流量系数Fig.9 Flow coefficients of swirler with different swirl numbers旋流数流量系数0.710.750.700.650.600.550.500.450.4 0.6 0.8 1.0 1.2 1.4 1.60.660.590.56-20 0 20 40 600.450.400.350.300.251800 2400 3000轴向速度/(m/s)径向距离/m网格数/104模拟值设计值(a)轴向速度(b)有效流通面积图 8 不同网格数量下的轴向速度和有效流通面积Fig.8 Axial velocity and effective area of swirler under different grid numbers400350300250ACd/mm214761762万万22022520万万现的，而流量系数与叶片安装角表现出反比例关系(式(5)。此外，流量系数还受叶片表面粗糙度、气流分离等因素影响。数值仿真中，流量系数主要通过有效流通面积及设计的物理面积计算得到。图 9 为 4 种不同旋流数旋流器的流量系数。从图中可以看出，随着旋流数的增大，流量系数表现出减小的趋势。这是因为文中旋流数的增大主要是通过增大叶片安装角来实(5)图 10 各旋流器方案中心截面的速度云图及流线图Fig.10 Velocity contour and streamline of the central section of each swirler case(a)SN=0.609(c)SN=1.004(b)SN=0.796(d)SN=1.240-30 -20 -10 0 10 20 30 40 50 60 70 80 90 100 110 120 130 140 150轴向速度/(m/s)对称旋涡中心隔离环射流-20 0 20 0.450.400.350.300.25径向速度/(m/s)SN=0.609SN=0.796SN=1.004SN=1.240Y/m(b)Z=60 mm第 35 卷 燃 气 涡 轮 试 验 与 研 究 17 存在避免了预燃级与主燃级旋流的相互影响。图 11 为旋流器中心截面不同轴向位置(距旋流器出口平面 26 mm，60 mm，91 mm)处的轴向速度。可以看到，随着旋流数增大，旋流器气流产生的中心低压区作用增强，旋流器中心截面轴向位置(径向位置 0.406 m)附近的回流速度增大；随着离旋流器出口距离的增加，回流速度减小。图 12 为旋流器中心截面不同轴向位置处的径向速度。可见，在距旋流器出口平面 26 mm 截面处，径向速度表现出非对称的一正一负 2 个峰值，且随着旋流数增大，旋流器中心轴线(径向位置 0.406 m)下侧的径向速度减小，原因是旋流器产生的旋转角动量会受到中心隔离环轴向动量的影响而被削弱；在距旋流器套筒出口平面 60 mm 及 91 mm 截面处，径向速度基本为负，这也是受中心隔离环射流影响所致。-20 0 20 0.450.400.350.300.25径向速度/(m/s)SN=0.609SN=0.796SN=1.004SN=1.240Y/m(a)Z=26 mm-20 0 20 40 60 80 1000.450.400.350.300.25轴向速度/(m/s)SN=0.609SN=0.796SN=1.004SN=1.240Y/m(a)Z=26 mm-20 0 20 40 60 80 1000.450.400.350.300.25轴向速度/(m/s)SN=0.609SN=0.796SN=1.004SN=1.240Y/m(b)Z=60 mm图 11 旋流器中心截面的轴向速度Fig.11 The axial velocity of the central section of the swirler-20 0 20 40 60 80 1000.450.400.350.300.25轴向速度/(m/s)SN=0.609SN=0.796SN=1.004SN=1.240Y/m(c)Z=91 mm图 12 旋流器中心截面的径向速度Fig.12 The radial velocity of the central section of the swirler-20 0 20 0.450.400.350.300.25径向速度/(m/s)SN=0.609SN=0.796SN=1.004SN=1.240Y/m(c)Z=91 mm点火电嘴位置图 15 燃油质量分数云图Fig.15 Fuel mass fraction contour 0.001 0.002 0.003 0.004 0.005 0.006 0.007 0.008 0.009 0.01 0.011 0.012 0.013(a)SN=0.609(c)SN=1.004(b)SN=0.796(d)SN=1.240燃油质量分数18 李 茂等：预燃级旋流数对双环腔燃烧室流场的影响 第 6 期 旋流器的空气回流量会对回流区的稳定及贫油熄火产生影响，为此为定量分析不同旋流数下的回流量，以旋流器套筒出口平面为参考平面，对 4 个方案的 17 个不同轴向位置的回流量16占旋流器套筒出口平面空气流量的百分比进行了对比，结果如图 13 所示。可见，4 种方案的回流量变化趋势都是先增大后减小、再增大再减小，其主要原因是流场结构表现为非对称的 3 涡结构。回流量在轴向位置455 mm 处达到最大，且随着旋流数的增加而增大，在旋流数 1.004、1.240 时达到 60%左右；后面回流量峰值较小，主要是由于在面积一定的情况下回流速度减小。高温燃气回流量大有利于更多的高温燃气回流，对燃烧室的点火有益。图 14 给出了不同旋流数下回流区长度及大小对比。从图中看，不同旋流数下的回流区长度基本一致，且回流区均止于主燃孔射流位置，这是由于主燃孔射流限制了回流区长度。因此，除了旋流器对回流区有影响外，主燃孔位置对回流区长度也有重要影响。此外，随着旋流数增大，回流区径向尺寸增大。3.3 旋流数对油雾场的影响 图 15 为不同旋流数下主燃区(红色虚线区域)旋流器中心截面的燃油质量分数云图。可以看到，4 种旋流数下，主燃区的燃油分布区别不大，主燃孔之后也有燃油分布；同时燃油沿预燃级旋流器中心轴线呈非对称分布，这是由于中心隔离环射流速度大，压力小，从而使得预燃级主燃区燃油分布存在向下偏移的情况。图 16 为不同轴向截面的燃油质量分数曲线示意图。可以看到，在距旋流器出口 26 mm 截面处，燃油质量分数随旋流数增大而增大；在距旋流器出口图 13 不同旋流器方案沿程回流量对比Fig.13 Comparison of recirculation air along flow channel for different swirler cases回流率/%轴向距离/mm806040200450 500 550SN=0.609SN=0.796SN=1.004SN=1.240图 14 不同旋流器方案的回流区大小Fig.14 Size of recirculation zone for different swirler cases403020100-10-20-300.40 0.45 0.50 0.55 0.60轴向速度/(m/s)(a)回流区长度Z/m0.460.440.420.400.380.360.34 0.45 0.50 0.55 Y/m(b)回流区大小SN=0.609SN=0.796SN=1.004SN=1.240Z/mSN=0.609SN=0.796SN=1.004SN=1.240第 35 卷 燃 气 涡 轮 试 验 与 研 究 1960 mm(点火电嘴中心轴线位置)截面处，随着旋流数增大，燃油质量分数减小。其原因在于，冷态流场为 3 涡结构，在距离旋流器出口稍远处，燃油的雾化蒸发还受到主燃孔射流和中心隔离环射流形成的旋涡影响。在小工况下，燃油流量较小，喷嘴雾化质量不好，而点火电嘴位置处的燃油质量分数受到 3 个旋涡剪切作用的共同影响，因此要合理选择旋流数，使得在点火电嘴位置处油气比更加接近化学恰当比。从上面结果分析，存在 1 个最佳旋流数使得点火电嘴位置处油气比最佳。4 结论 采用数值仿真方法，研究了双环腔燃烧室的流量系数，主燃区的轴向速度、径向速度、回流流量、回流区大小、燃油质量分数等参数，间接分析了其对点火性能、燃烧性能的影响，主要得出如下结论：(1)在本文研究范围内，轴向旋流器的旋流数对冷态流场结构影响较小，整个流场为 3 涡结构，在旋流器出口形成 2 个旋涡，在主燃孔射流及中心隔离环射流之间形成 1 个旋涡。(2)气流的回流速度和回流量随着旋流数的增大而增大。回流速度的增大有利于增强与燃油的剪切作用，在一定程度上改善燃油雾化质量，从而改善点火性能；回流量的增大意味着卷吸了更多的高温燃气回流，有助于形成稳定的回流区火焰，对燃烧室的点火有益。(3)受主燃孔射流限制，不同旋流数下的回流区长度基本一致，但随着旋流数增大，回流区径向尺寸增大。(4)在旋流器出口附近，燃油质量分数随旋流数增大而增大；在离旋流器出口较远处，燃油质量分数随旋流数增大而减小。受冷态流场 3 涡结构影响，存在 1 个最佳旋流数使得点火电嘴位置处油气比最佳。参考文献：1 23456 7891011Walker A D，Denman P A，McGuirk J J.Experimental and computational study of hybrid diffusers for gas turbine combustorsR.ASME GT-2003-38406，2003.赵坚行.民用发动机污染排放及低污染燃烧技术发展趋势 J.航空动力学报，2008，23(6)：986996.尚守堂，程 明，张军峰，等.低排放长寿命燃烧室关键技术分析 C/.大型飞机关键技术高层论坛暨中国航空学会 2007 年学术年会论文集.2007.张 斌.航空发动机燃烧室点火模拟及特性研究 D.南京：南京航空航天大学，2004.江 洋.燃烧室点火过程大涡数值模拟方法研究 D.南京：南京航空航天大学，2017.林宇震，许全宏，刘高恩.燃气轮机燃烧室 M.北京：国防工业出版社，2008.金如山，索建秦.先进燃气轮机燃烧室 M.北京：航空工业出版社，2016.PanduRanga Reddy A，Sujith R I，Chakravarthy S R.Swirlerflow field characteristics in a sudden expansion combustor geometry using PIVR.AIAA 2005-217，2005.Cai J，Jeng S M，Tacina R.The structure of a swirl-stabilized reacting spray issued from an axial swirlerR.AIAA 2005-1424，2005.蔡文祥，胡好生，赵坚行.涡流器燃烧室头部两相反应流数值模拟 J.航空动力学报，2006，21(5)：837842.航空发动机设计手册总编委会.航空发动机设计手册 第 9 册 主燃烧室 M.北京：国防工业出版社，2000.0.0080.0060.0040.00200.26 0.30 0.34 0.38 0.42 0.46 0.50燃油质量分数Y/m(a)Z=26 mm0.26 0.30 0.34 0.38 0.42 0.46 0.50图 16 不同轴向截面的燃油质量分数Fig.16 Fuel mass fraction of different axial sections0.0080.0060.0040.0020(b)Z=60 mm燃油质量分数Y/mSN=0.609SN=0.796SN=1.004SN=1.240SN=0.609SN=0.796SN=1.004SN=1.240(下转第 33 页)第 35 卷 燃 气 涡 轮 试 验 与 研 究 331213141516Tan H P，Shuai Y，Dong SK.Analysis of rocket plume base heating by using backward Monte-Carlo methodJ.Journal of Thermophysics and Heat Transfer，2005，19(1)：125127.Shuai Y，Dong S K，Tan H，P.Simulation of the infrared radiation characteristics of high-temperature exhaust plume including particles using the backward Monte Carlo methodJ.Journal of Quantitative Spectroscopy&Radiative Transfer，2005，95：231240.帅 永，董士奎，刘林华.高温含粒子自由流红外辐射特性的反向蒙特卡罗法模拟J.红外与毫米波学报，2005，24(2)：100104.帅 永，董士奎，刘林华，等.双向蒙特卡罗法模拟参与性介质中的辐射换热 J.工程热物理学报，2006，27(1)：109111.李宏顺.用逆向蒙特卡罗法计算光的辐射强度组分17181920 21 J.华中科技大学学报(自然科学版)，2007，35(2)：113115.黄群星，王 飞，常瑞丽，等.基于逆向蒙特卡罗法的煤粉炉内辐射能传递计算 J.中国电机工程学报，2007，27(8)：8893.吕建伟，王 强.飞机蒙皮红外辐射特征的反向蒙特卡罗计算与分析方法 J.红外与激光工程，2009，38(2)：232237.单 勇，张靖周，邵万仁，等.涡扇发动机排气系统中心锥气膜冷却结构的气动和红外辐射特性实验 J.航空动力学报，2012，27(1)：1015.罗明东.无人机排气系统红外隐身技术研究D.南京：南京航空航天大学，2006.黄 伟，吉洪湖，斯 仁，等.涡扇发动机排气系统红外特征 J.推进技术，2010，31(6)：745750.12 13 14 AI-Abdeli Y M，Masri A R.Stability characteristics and flow fields of turbulent non premixed swirling flamesJ.Combustion Theory and Modeling，2003，7：731766.Kalt P A M，AI-Abdeli Y M，Masri A R，et al.Swirling turbulent non-premixed flames of methane:flowfield and compositional structureC/.Proceedings of the Twenty-Ninth International Symposium on Combustion.2002：19131919.Kim S E，Choudhury D，Patel B.Computations of complex turbulent flows using the commercial code FLUENTM.Dordrecht：Springer Netherlands，1999：259276.齐海帆.燃烧室头部设计对燃烧室性能影响研究 D.南京：南京航空航天大学，2012.代 威，林宇震，张 驰.第 2 级径向旋流器旋流数对燃烧室点火和贫油熄火性能的影响 J.航空动力学报，2015，30(5)：10921098.15 16(上接第 19 页)