某重型燃机燃烧室热态流场数值模拟
王嘉奇;冯兆兴;王国峰;刘洋;史俊瑞
【摘 要】为了更好地了解燃机燃烧室中燃烧过程的细节以及流场特性,以重型燃机为研究对象,建立燃烧室的三维模型,采用单步化学反应机理对该燃机燃烧室热态流场进行数值模拟.结果表明:燃烧室轴向中心截面的最高温度和最大速度约为2500 K和120 m/s;燃烧室的速度场、温度场、反应物、生成物分布能较好地反映出该燃机燃烧室的实际情况,计算得到的出口温度分布与理论值近似,为同类燃机燃烧室的数值模拟研究提供了参考价值.
【期刊名称】《沈阳工程学院学报（自然科学版）》 【年(卷),期】2018(014)004 【总页数】5页(P305-309)
【关键词】燃烧室;热态流场;数值研究 【作 者】王嘉奇;冯兆兴;王国峰;刘洋;史俊瑞
【作者单位】沈阳工程学院研究生部,辽宁 沈阳 110136;沈阳工程学院能源与动力学院,辽宁 沈阳 110136;沈阳工程学院能源与动力学院,辽宁 沈阳 110136;沈阳工程学院研究生部,辽宁 沈阳 110136;山东理工大学交通与车辆工程学院,山东 淄博 255000
【正文语种】中 文 【中图分类】TK16
1 物理模型
燃烧室为沿周向分布20个火焰筒的逆流式环管型燃烧室，单个火焰筒的整体结构如图1所示。扩压器与压气机后气缸相连，内含出口静叶、整流叶片及长短不一的支板。燃烧室外壳直径最大值为3 100 mm，内腔直径为2 968 mm，全长为900 mm，质量约为1 320 kg。 图1 某型燃机单个火焰筒的整体结构 2 数值模型 2.1 控制方程
求解重型燃机内部的三维流动过程主要是求解Reynolds时均化Navier-Stokes方程及关联量的输运方程的湍流模型，即标准k-ε湍流模型。计算过程中必须遵循4个基本控制方程。 1）连续性方程
该燃气轮机燃烧室的流动问题属于稳态燃烧问题 2）动量方程
3）能量方程
式中，keff为有效热传导系数；Jj为第j种组分的扩散通量；E为总能，是势能和内能之和；Sh为化学反应热和其他用户定义的体积热源项。 4）组分方程
式中，Yi代表第i组分的质量分数；Ri代表化学反应中组分i的净生成速率；Si代
表组分i由离散项和其他用户定义的源项引起额外的生成速率。 2.2 模型建立和网格划分
针对该燃烧室，运用Pro-E软件创建三维仿真分析模型。对于火焰筒上的细小气膜孔，采用气膜带简化处理，这样既能减少网格的数量，又能保证空气的冷却效果。对火焰筒头部的主旋流器和8个微型旋流器采用增加空气入口旋流的方法进行简化，利用空气入口切向、轴向、径向3个方向的速度分量来实现旋流，保证气体经过旋流器后，过量空气系数达到设计要求。
对火焰筒环形区的微型旋流器以及的主旋流器等部位进行加密处理，在温度和速度骤变的地方，采用了细密的网格；对于流体流动影响不大的地方，采用较疏的网格，减少网格数量，从而减少计算量。燃烧室头部网格如图2所示。 图2 燃烧室头部网格
该模型选择的网格质量指标是EquiSize Skew，是通过单元大小计算歪斜度以检查网格的质量，EquiSize Skew在0～0.4之间（lower-upper）的网格数占比为96.70%，说明网格划分质量良好。 2.3 边界条件 1）空气进口
该燃烧室设置了带有旋流且供燃料燃烧的空气入口、15排气膜带入口和4个冷却孔入口。所有的空气入口都设置为质量流量进口，给定每个进口的质量流量数值、方向、湍流度和水力直径。给定空气的温度为666 K，空气的总进气量为6.59 kg/s。 2）燃料进口
燃料进口包括8个微型旋流器、16个值班燃料小孔入口以及8个燃料喷杆入口。所有燃料入口边界也设置为质量流量进口，给定质量流量的数值、湍流度及水力直径。设定燃料入口温度为288 K，燃气的总进气量为102 g/s。
3）出口
出口边界设置为压力出口，给定湍流度、回流度和水力直径。 4）壁面
热边界类型设置为Radiation类型，假设壁面处无速度滑移，热生成率设为0.3 W/m3，壁厚为0 mm，内部发射率为0.6，外部发射率为0.8。 3 计算结果与分析 3.1 速度场分布特性
图3为数值分析得到的速度分布云图及矢量图。燃烧室在喷嘴的出口处形成了回流区，这是由于高速旋转的气流对喷嘴出口中心部位的气流产生卷吸效应，从而在背后形成一个低压区域，这对于旋流燃烧室的稳定燃烧至关重要。同时，回流也延长了燃料在火焰筒内的停留时间，可以提高燃料的燃烧效率。回流区下游截面的平均速度在60 m/s左右，并且轴向速度由前往后逐渐增加，这是由于随着燃料与空气燃烧反应的进行，燃气平均温度逐渐升高，使得火焰筒内气体的体积流量增加，而且火焰筒出口通流面积渐缩，导致轴向速度增加，出口平均速度达到100 m/s左右。
图3 燃烧室内轴向中心截面流场分布特性 3.2 组分浓度场分布特性
图4 燃烧室内中心截面的组分分布特性
图4为数值分析得到的各组分的浓度分布图。从图4中可以看出，CH4燃料喷入火焰筒后，进行快速燃烧，浓度急剧下降。O2的浓度场则与燃烧速度分布有关，在燃烧的中心区域，由于燃烧作用，O2消耗多，浓度低；外围区域，燃烧反应不剧烈，O2的浓度高。火焰筒头部的主燃区O2浓度最低，环形燃烧区后半部也出现O2浓度低的区域，此区域达到燃料燃烧的温度，O2因燃烧而消耗掉。 两种燃烧产物CO2和H2O浓度场分布规律基本相同，在主燃区浓度大，最大的
质量分数相差不多。CO2和H2O分布大致以中轴线对称，从第一级束腰环开始，CO2和H2O浓度开始逐渐降低，由掺混孔进入的大量空气使得CO2和H2O浓度迅速下降。出口截面，CO2的质量分数降到0.04左右，H2O的质量分数降到0.034，这是因为此区域是燃烧剧烈区域。反应物和生成物的分布均能较好地反映出该型燃气轮机燃烧室的实际燃烧情况。 3.3 温度场分布特性
图5为利用不同燃烧反应机理模型得到的温度分布。CH4与空气在环形燃烧室的前半部进行预混，在后半部进行燃烧。因此，环形燃烧室的前半部温度低，而后半部的温度高。燃烧主要稳定在主燃区的头部，热点温度达到2 500 K左右。出口截面的温度分布梯度较小，说明温度分布较均匀。出口截面的平均温度为1 295.6 K，热点温度达到1 625.8 K。单步化学反应机理、简化机理、详细化学反应机理和利用GR13.0计算的温度场的温度分布基本一致，火焰形状基本相同，燃烧室内的温度分布大致以轴线为对称轴分布。该数据满足出口温度涡轮导向叶片的温度要求，模拟得到的燃烧室出口的高温区的分布位置与给定实验得到的高温区分布位置相同，温度略高于实验的数据。 图5 燃烧室轴向截面的温度分布 4 结 论
1）采用单步化学反应机理模拟得到的燃烧室结果表明，旋流和回流位置分布合理，燃气与空气的混合情况良好，各组分分布均能较好地反映出燃烧室的实际燃烧情况。 2）与不同燃烧反应机理的模拟结果对比后，验证了单步化学反应机理模拟的温度场分布特性合理，出口温度场与给定实验出口温度场吻合。
3）该数值模拟分析方法计算速度快，对模拟网格要求低，是一种能够适合燃烧室进行快速分析的有效方法。 参考文献
【相关文献】
［1］姜 伟，赵士杭.燃气轮机原理、结构与应用［M］.北京：科学出版社，2002. ［2］张宝诚.航空发动机实验和测试技术［M］.北京：北京航空航天大学出版社，2005. ［3］VAUCHELLES D，TAUPIN B，BOUKHALFAA.Experimental study of lean premixed turbulent combustion in a scale gas turbine chamber［J］.Experimental Thermal and Fluid Science，2004，28（7）：683-690.
［4］TYLISZCZAK A，BOGUSLAWSKI A，NOWAK D.Numerical simulations of combustion process in a gas turbine with a single and multi-point fuel injection system［J］.Applied Energy，2016，174：153-165.
［5］ BULAT G，JONES W P，MARQUIS A J.NO and CO formation in an industrial gas-turbine combustion chamber using LES with the Eulerian sub-grid PDF method［J］.Combustion and Flame，2014，161（7）：1804-1825.
［6］史 丹.我国能源经济的总体特征、问题及展望［J］.中国能源，2007，29（1）：5-12。 ［7］张晨曦.三旋流燃烧室头部方案设计与燃烧流场数值模拟［D］.哈尔滨：哈尔滨工程大学，2013.
［8］王成军，艾延廷，张宝诚.航改燃气轮机环形燃烧室三维全流场数值模拟［J］.航空发动机，2009，35（6）：22-24.
［9］董 刚，黄 鹰，陈义良.不同化学反应机理对甲烷射流湍流扩散火焰计算结果影响的研究［J］.燃料化学学报，2000，28（1）：49-54.
［10］范珍涔，程 波，曹宗华，等.反应机理对燃烧室出口温度分布仿真结果的影响［J］.燃气涡轮试验与研究，2016，29（4）：17-21.
［11］梁 双.甲烷燃料的燃烧机理及其在燃气轮机中的应用［D］.沈阳：沈阳航空航天大学，2013. ［12］潘宏刚，郑飞飞，夏永放，等.汽轮机转子临界转速实验研究［J］.沈阳工程学院学报：自然科学版，2016，12（4）：314-318.
