航空发动机两级反向旋流燃烧室燃烧流场大涡模拟研究

推进技术 ›› 2018, Vol. 39 ›› Issue (7): 1576-1589.

航空发动机两级反向旋流燃烧室燃烧流场大涡模拟研究

周瑜¹,乐嘉陵^1,2,黄渊²

西北工业大学动力与能源学院，陕西西安 710072,西北工业大学动力与能源学院，陕西西安 710072;中国空气动力研究与发展中心超高速空气动力研究所高超声速冲压发动机技术重点实验室，四川绵阳 621000,中国空气动力研究与发展中心超高速空气动力研究所高超声速冲压发动机技术重点实验室，四川绵阳 621000

发布日期:2021-08-15
基金资助:
国家自然科学基金( 9164120019)。通讯作者:周瑜，男，博士生，研究领域为航空发动机燃烧室数值模拟。 E-mail: zhouyu1984@126.com

LESof Combustion Flow Fieldina Practical Aeroengine Combustorwith Two-Stage Counter-Rotating Swirler

School of Power and Energy，Northwestern Polytechnical University，Xi’an 710072，China,School of Power and Energy，Northwestern Polytechnical University，Xi’an 710072，China;Science and Technology on Scramjet Laboratory，Hypervelocity Aerodynamics Institute of China Aerodynamics Researchand Development Center，Mianyang 621000，China and Science and Technology on Scramjet Laboratory，Hypervelocity Aerodynamics Institute of China Aerodynamics Researchand Development Center，Mianyang 621000，China

Published:2021-08-15

摘要/Abstract

摘要： 为深入了解真实航空发动机内燃烧流场，采用动态亚网格模型结合单步快速化学反应(FC)、火焰面( FLM)及反应进度变量( FPV)等三种燃烧模型对径向两级反向旋流燃烧室单头部构型进行了大涡模拟。为避免模型简化误差，对燃烧室包括全部气膜冷却孔在内的精细结构进行了完全仿真。在已达到统计定常状态的冷态流场基础上首先模拟了燃料喷注和掺混过程，约2.6ms后燃料到达真实的点火位置，随后采用FPV模型在半径3mm的球形区域数值模拟了点火，展示了在主燃孔横向射流作用下初始火焰沿化学恰当比混合分数等值线传播并充满整个火焰筒的发展过程，结果显示火焰到达燃烧室出口的耗时约为6～7ms。不同模型算法预测的平均温度场与CARS测量结果作了对比，LES-FPV，RANS-FPV，LES-FLM 以及参考文献中RANS-FLM 计算平均误差分别为3.47%，4.17%，7.76% 和11.22%，表明LES改进了模拟精度，且FPV模型显著优于FLM模型。RANS-FPV预测的出口存在严重热斑，导致其给出的出口温度分布因子(OTDF) 及最大径向温度分布因子(RTDF) 值分别达到0.593和0.313；LES-FPV结果均匀性最好，其预测值分别为0.284和0.193。

关键词: 航空发动机燃烧室；反向双旋流；大涡模拟；动态亚网格模型；快速化学反应模型；火焰面模型；反应进度变量模型

Abstract: To improve knowledge of the reacting flow in real aeroengines，large eddy simulation(LES) with dynamic Smagorinsky subgrid model is used to explore the complex unsteady flowfield in a single burner of atypical aeroengine combustor with radial two-stage counter-rotating swirler. Three combustion models namelyfast chemistry model(FC)，flamelet model(FLM)and progress-variable model(FPV)are applied. The com.plex geometric configuration including all film cooling holes is fully simulated without any conventional simplifica.tion in order to essentially reduce the modeling errors. Based on the non-reacting flowfield which has been ad.vanced to the statistically stationary state in early work，unsteady process of fuel injection is first simulated with.out reaction. It cost about 2.6ms for adequate fuel reaching the real ignition location，then pseudo ignition is nu.merically conducted using FPV model in the given sphere region with a radius of 3mm. The process that initialflame propagates along the isoline of stoichiometric mixture and fills the whole chamber under the effect of trans.verse flow is simulated. It costs about 6～7ms for the flame arriving at the outlet. Compared to coherent anti-Stokes Raman scattering(CARS)measurement of mean temperature on central profile，mean error of LES-FPV， RANS-FPV，LES-FLM and RANS-FLM in the reference are 3.47%，4.17%，7.76% and 11.22% respectively， indicating that LES improves the accuracy and FPV performs better than FLM. Due to the severe hot spot at theoutlet，RANS-FPV gives the poorest outlet temperature distribution factor(OTDF)and maximum radial tempera. ture distribution factor(RTDF)as 0.593 and 0.313 respectively. Corresponding values predicted by LES-FPV are 0.284 and 0.193，presenting the best uniformity among all cases.

Key words: Aeroengine combustor；Counter-rotating swirler；Large eddy simulation；Dynamic Smago. rinsky model；Fast chemistry model；Flamelet model；Progress variable model

周瑜,乐嘉陵,黄渊. 航空发动机两级反向旋流燃烧室燃烧流场大涡模拟研究[J]. 推进技术, 2018, 39(7): 1576-1589.

ZHOU Yu¹,Le Jia-ling^1,2,HUANG Yuan². LESof Combustion Flow Fieldina Practical Aeroengine Combustorwith Two-Stage Counter-Rotating Swirler[J]. Journal of Propulsion Technology, 2018, 39(7): 1576-1589.

参考文献

[1] Yang Z.Large-Eddy Simulation: Past,Present and the Future[J].Chinese Journal of Aeronautics,2015,28(1): 11-24.
[2] Kim S,Menon S.Large-Eddy Simulation of a High-Pressure Single-Element Lean Direct-Injected Gas-Tur.bine Combustor[R].AIAA 2014-0131.
[3] Bulat G,Jones W P,Navarro-Martinez S.Large EddySimulations of Isothermal Confined Swirling Flow in anIndustrial Gas-Turbine[J].International Journal of Heat and Fluid Flow,2015,51:50-64.
[4] Filosa A,NollB E,Di Domenico M,et al.Numerical Investigations of a Low Emission Gas Turbine CombustorUsing Detailed Chemistry[R].AIAA 2014-3916.
[5] Ilie M.Numerical Computations of Turbulent SwirlingReacting Flows Using URANS and LES[R].AIAA 2009-5407.
[6] Widenhorn A,Noll B,Aigner M.Numerical Study of a Non -Reacting Turbulent Flow in a Gas Turbine ModelCombustor[R].AIAA 2009-647.
[7] Javadi A,Nilsson H.Time-Accurate Numerical Simula.tions of Swirling Flow with Rotor-Stator Interaction[J].Flow Turbulence Combust,2015,95: 755-774.
[8] See Y C,Ihme M.LES Investigation of Flow Field Sen.sitivity in a Gas Turbine Model Combustor[R].AIAA 2014-0621.
[9] Kim J C,Sung H G.Large Eddy Simulation of the Tur.bulent Flow Field in a Swirl Stabilized Annular Combus.tor[R].AIAA 2009-645.
[10] Fureby C.LES of a Multi-Burner Annular Gas Turbine Combustor[J].Flow Turbulence Combust,2010,84: 543-564.
[11] Tachibana S,Saito K,Yamamoto T,et al.Experimen.tal and Numerical Investigation of Thermo-Acoustic In.stability in a Liquid-Fuel Aero-Engine Combustor at El.evated Pressure: Validity of Large-Eddy Simulation ofSpray Combustion[J].Combustion and Flame,2015, 162: 2621-2637.
[12] Hermeth S,Staffelbach G,Gicquel L Y M,et al.Bi.stable Swirled Flames and Influence on Flame Transfer Functions[J].Combustion and Flame,2014,161: 184-196.
[13] Zeng W,Li H,Chen B,et al.Kinetic Simulation of Combustion Process in The Combustor of the Aero-En.gine[J].Combustion Science and Technology,2014, 186(8): 1097-1114.
[14] Ruggles A,Kelman J.A Gas Turbine Combustor for In.stability Research and LES Validation: Methods andMean Results[J].Combustion Science and Technology, 2014,186(3): 313-331.
[15] Worth N A,Dawson J R.Modal Dynamics of Self-Excit.ed Azimuthal Instabilities in an Annular Combustion Chamber[J].Combustion and Flame, 2013, 160: 2476-2489.
[16] Motheau E,Nicoud F,Poinsot T.Mixed Acoustic-En.tropy Combustion Instabilities in Gas Turbines[J].Jour.nal of Fluid Mechanics,2014,749: 542-576.
[17] Brankovic A,Ryder R C Jr,Petersen E L.Time-De.pendent Simulation of Fuel-Flexible Gas-Turbine Burn.ers[R].AIAA 2006-806.
[18] Neophytou A,Cuenot B.Large-Eddy Simulation of Ig.nition and Flame Propagation in a Trisector Combustor[J].Journal of Propulsion and Power,2016,32(2): 345-359.
[19] Masquelet M M,Cao S,Pai M.Ignition and FlamePropagation in a Liquid-Fueled Swirling Combustor[R].AIAA 2016-4687.
[20] Barré D,Esclapez L,Cordier M,et al.Flame Propaga.tion in Aeronautical Swirled Multi-Burners: Experimen.tal and Numerical Investigation[J].Combustion and Flame,2014,161: 2387-2405.
[21] Koo H,Raman V.Large-Eddy Simulation of a Turbu.lent Sooting Flame in a Swirling Combustor[R].AIAA 2015-0167.
[22] Fedina E,Fureby C,Borzov S A,et al.Combustion LES of CESAR Multi-Burner Annular Combustor[R].AIAA 2011-785.
[23] Zettervall N,Fedina E,Nordin-Bates K,et al.Com.bustion LES of a Multi-Burner Annular Aero-EngineCombustor Using a Skeletal Reaction Mechanism for Jet-A Air Mixtures[R].AIAA 2015-4020.
[24] Lietz C,Heye C,Raman V.Flame Stability Analysis inan Ultra Compact Combustor Using Large-Eddy Simula.tion[R].AIAA 2014-1022.
[25] Yang S,Wang X,Yang V,et al.Comparison of Flame.let-Progress-Variable and Finite-Rate Chemistry LESModels in a Preconditioning Scheme[R].AIAA 2017-0605.
[26] See Y C,Ihme M.Large Eddy Simulation of a Gas Tur.bine Model Combustor[R].AIAA 2013-0172.
[27] Pierce C D.Progress-Variable for Large-Eddy Simula.tion of Turbulent Combustion[D].California: Stanford University,2001.
[28] YAN Ying-wen,ZHAO Jian-xing,ZHANG Jing-zhou, et al.Large-Eddy Simulation of Turbulent ReactingFlow in the Combustor by Three-Dimensional Body-Fit.ted Grid[J].Journal of Propulsion Technology,2005, 26(3): 219-222.(颜应文,赵坚行,张靖周,等 .三维贴体坐标系下燃烧室热态流场的大涡模拟[ J].推进技术, 2005,26(3): 219-222.)
[29] YAN Ying-wen,LIU Yong,ZHAO Jian-xing,etal.A Parallel Algorithm for Large Eddy Simulation of Two-Phase Reacting Flows in Model Combustor[J].Journal of Propulsion Technology,2008,29(4): 431-437.(颜应文,刘勇,赵坚行,等 .模型燃烧室两相燃烧大涡模拟的并行计算[ J].推进技术, 2008,29(4): 431-437.)
[30] 何跃龙,邓远灏,颜应文,等 .大涡模拟模型燃烧室燃烧性能计算[ J].航空动力学报, 2012,27(9): 1939-1947.
[31] ZHOU Yu,LE Jia-ling,CHEN Liu-jun,et al.LargeEddy Simulation of Aeroengine Combustor with Counter-Rotating Swirler[J].Journal of Propulsion Technology, 2017,38(4): 909-917.(周瑜,乐嘉陵,陈柳君,等 .径向双旋流燃烧室流场结构大涡模拟研究[ J].推进技术, 2017,38(4): 909-917.)
[32] 陈敏,宋文艳,肖隐利,等 .斜切径向旋流燃烧室主燃区光学测量与特性分析[ J].航空动力学报, 2013,28(8): 1727-1735.* 收稿日期: 2017-11-02；修订日期: 2017-12-23。基金项目:国家自然科学基金( 9164120019)。通讯作者:周瑜,男,博士生,研究领域为航空发动机燃烧室数值模拟。 E-mail: zhouyu1984@126.com (编辑:刘萝威)