LPP燃烧室燃烧不稳定特征三维模拟

李昊; 刘勇; 张祥; 王旭怀; 杨晨; 刘重阳

doi:10.13224/j.cnki.jasp.20220250

LPP燃烧室燃烧不稳定特征三维模拟

doi: 10.13224/j.cnki.jasp.20220250

李昊^1,,
刘勇¹,
张祥¹,
王旭怀¹,
杨晨¹,
刘重阳^{2, 3}

1.
南京航空航天大学能源与动力学院，南京 210016
2.
西北工业大学动力与能源学院，西安 710072
3.
中国航发四川燃气涡轮研究院，四川绵阳 621000

基金项目: 中国航发自主创新专项资金（ZZCX-2019-009）

详细信息

作者简介:
李昊（1998-），男，硕士生，主要研究方向为燃烧室不稳定燃烧。E-mail：leelihao@126.com

中图分类号: V231.2
计量
- 文章访问数: 115
- HTML浏览量: 63
- PDF量: 54
- 被引次数: 0
出版历程
- 收稿日期: 2022-04-25
- 网络出版日期: 2023-09-27

Three-dimensional simulation of combustion instability characteristics in LPP combustor

1.
College of Energy and Power Engineering，Nanjing University of Aeronautics and Astronautics，Nanjing 210016，China
2.
School of Power and Energy，Northwestern Polytechnical University，Xi’an 710072，China
3.
Sichuan Gas Turbine Research Institute，Aero Engine Corporation of China，Mianyang Sichuan 621000，China

摘要

摘要:
为分析贫油预混预蒸发（LPP）燃烧室的燃烧不稳定（CI）特征，通过对三维亥姆霍兹方程进行了三种不同程度的简化：平均温度场和导入CFD温度场的无源项方程，以及导入燃烧流场特征的有源项方程，分别对单头部LPP燃烧室模型进行了三维频域特征数值仿真。结果表明：燃烧室内的温度分布是燃烧室声学特征频率的重要影响因素，释热率源项对主频无影响。相比于仅设置平均温度场，导入CFD三维温度场可以获得与实验频率更吻合的结果，精度提高了5%。采用解耦方式求解频域方程能够快速建立声学系统与燃烧流场间的联系，释热率和迟滞时间的空间分布特征表现在亥姆霍兹方程的源项，其对预测燃烧室固有频率没有影响，但是能够获得详细声压分布特征。
- 贫油预混预蒸发（LPP）燃烧室 /
- 声学特征频率 /
- 燃烧不稳定（CI） /
- 声学有限元计算 /
- 热声耦合
Abstract:
In order to analyze the combustion instability (CI) characteristics of lean premixed pre-vaporized (LPP) combustor, three-dimensional Helmholtz equations were simplified at three different levels. The average temperature field equation without source term, the imported CFD temperature field equation without source term, and the imported combustion flow field characteristics equation with source term, were respectively simulated for the single-head LPP combustor model in three-dimensional frequency domain. The results showed that the temperature distribution in the combustor was an important factor affecting the acoustic eigenfrequency of the combustor, and the source term of heat release rate had no effect on the main frequency. Compared with only setting of average temperature field, importing the 3D temperature field calculated by CFD can obtain more consistent results with the experimental frequency, and the accuracy was increased by 5%. The relationship between the acoustic system and the combustion flow field can be quickly established by solving the frequency domain equation in a decoupling way. The spatial distribution characteristics of heat release rate and hysteresis time were represented by the source term of the Helmholtz equation, which had no effect on the prediction of the natural frequency of the combustion chamber, but the detailed sound pressure distribution characteristics can be obtained.
- lean premixed pre-vaporized (LPP) combustor /
- acoustic eigenfrequency /
- combustion instability (CI) /
- finite element calculation /
- thermal acoustic coupling

HTML

图 1 燃烧室段横截面示意图

Figure 1. Schematic diagram of combustor cross-section

下载: 全尺寸图片幻灯片

图 2 实验系统示意图

Figure 2. Schematic diagram of experimental system

下载: 全尺寸图片幻灯片

图 3 实验测得CI频率

Figure 3. CI frequency measured from experiment

下载: 全尺寸图片幻灯片

图 4 燃烧室段网格划分

Figure 4. Meshing of combustor

下载: 全尺寸图片幻灯片

图 5 以阶跃函数设置模型温度

Figure 5. Setting model temperature by step function

下载: 全尺寸图片幻灯片

图 6 无源平均温度场模拟与实验结果对比

Figure 6. Comparison between simulation and experimental results of mean temperature field without source term

下载: 全尺寸图片幻灯片

图 7 PIV测量与RANS计算冷态流场轴向速度分布

Figure 7. Axial velocity distribution of cold flow field from PIV measurements and RANS calculation

下载: 全尺寸图片幻灯片

图 8 无量纲温度

Figure 8. Dimensionless temperature

下载: 全尺寸图片幻灯片

图 9 1阶纵向振型声压等值面

Figure 9. Sound pressure iso-surface of the first order longitudinal mode

下载: 全尺寸图片幻灯片

图 10 无源详细温度场模拟与实验结果对比

Figure 10. Comparison between simulation and experimental results of detailed temperature field without source term

下载: 全尺寸图片幻灯片

图 11 无量纲反应速率

Figure 11. Dimensionless rate of reaction

下载: 全尺寸图片幻灯片

图 12 无量纲迟滞时间

Figure 12. Dimensionless delay time

下载: 全尺寸图片幻灯片

图 13 不同频率声压级响应

Figure 13. Sound pressure level response at different frequencies

下载: 全尺寸图片幻灯片

图 14 速度脉动幅值与增益声压级响应

Figure 14. Sound pressure level response of velocity pulsation amplitude and gain

下载: 全尺寸图片幻灯片

表 1 实验工况及结果

Table 1. Experimental conditions and results

工况	进气温度/K	出口温度/K	油气比	振荡频率/Hz
1	345	760	0.036	116
2	345	1120	0.041	141
3	345	1140	0.045	142
4	383	820	0.036	121
5	383	1200	0.040	146
6	383	1230	0.044	148
7	402	850	0.034	123
8	402	1240	0.039	148
9	402	1380	0.045	157

下载: 导出CSV

表 2 网格无关性验证结果

Table 2. Grid independence verification results

网格数量/10⁴	1阶频率/Hz	2阶频率/Hz	3阶频率/Hz
11	96.804	189.35	278.01
16	96.803	189.35	278
21	96.801	189.32	277.95
26	96.801	189.32	277.95
30	96.801	189.32	277.95

下载: 导出CSV

表 3 无源平均温度场模拟与实验结果

Table 3. Simulation and experimental results of mean temperature field without source term

工况编号	特征频率/Hz		误差/%
工况编号	实验	模拟	误差/%
1	116	130.6	12.6
2	141	147.55	4.6
3	142	148.18	4.4
4	121	136.78	13.0
5	146	153.73	5.3
6	148	155.05	4.8
7	123	139.44	13.4
8	148	156.43	5.7
9	157	162.43	3.5

下载: 导出CSV

表 4 无源详细温度场模拟与实验结果

Table 4. Simulation and experimental results of detailed temperature field without source term

工况编号	特征频率/Hz		误差/%
工况编号	实验	模拟	误差/%
1	116	119.07	2.6
2	141	143.19	1.60
3	142	144.8	2.0
4	121	126.67	4.7
5	146	147.9	1.3
6	148	149.79	1.2
7	123	128.71	4.6
8	148	150.01	1.4
9	157	152.77	2.7

下载: 导出CSV

参考文献(27)

[1]	李磊,孙晓峰. 推进动力系统燃烧不稳定性产生的机理、预测及控制方法[J]. 推进技术,2010,31(6): 710-720. LI Lei,SUN Xiaofeng. Mechanism, prediction and control method of combustion instability in propulsion system[J]. Journal of Propulsion Technology,2010,31(6): 710-720. (in Chinese)
[2]	张弛,林宇震,徐华胜,等. 民用航空发动机低排放燃烧室技术发展现状及水平[J]. 航空学报,2014,35(2): 332-350. ZHANG Chi,LIN Yuzhen,XU Huasheng,et al. Development status and level of low emissions combustor technologies for civil aero-engine[J]. Acta Aeronautica et Astronautica Sinica,2014,35(2): 332-350. (in Chinese)
[3]	朱民,DOWLING A P,BRAY K N C. 振荡燃烧过程的计算和诊断[J]. 工程热物理学报,2007,28(1): 165-168. ZHU Min,DOWLING A P,BRAY K N C. Calculations and predications of combustion oscillations[J]. Journal of Engineering Thermophysics,2007,28(1): 165-168. (in Chinese)
[4]	张志浩,刘潇,郑洪涛,等. 燃烧不稳定性实验诊断及数值分析方法研究进展[J]. 燃气轮机技术,2021,34(3): 1-13. ZHANG Zhihao,LIU Xiao,ZHENG Hongtao,et al. Research progress of experimental diagnosis and numerical analysis methods of combustion instability[J]. Gas Turbine Technology,2021,34(3): 1-13. (in Chinese)
[5]	于丹,郭志辉,杨甫江. 贫燃预混燃烧室中的分布式火焰传递函数分析[J]. 推进技术,2016,37(12): 2210-2218. YU Dan,GUO Zhihui,YANG Fujiang. Experimentally study on distribution of flame transfer function in a lean premixed combustor[J]. Journal of Propulsion Technology,2016,37(12): 2210-2218. (in Chinese)
[6]	唐豪杰,高原,朱民. 预混火焰传递函数的测量与分析[J]. 工程热物理学报,2011,32(1): 173-176. TANG Haojie,GAO Yuan,ZHU Min. Measurements and analysis of premixed flame transfer function[J]. Journal of Engineering Thermophysics,2011,32(1): 173-176. (in Chinese)
[7]	PASCHEREIT C, FLOHR P, SCHUERMANS B. Prediction of combustion oscillations in gas turbine combustors[R]. AIAA2001-484, 2001.
[8]	KOPITZ J, HUBER A, SATTELMAYER T, et al. Thermoacoustic stability analysis of an annular combustion chamber with acoustic low order modeling and validation against experiment[C]//Proceedings of ASME Turbo Expo: Power for Land, Sea, and Air. Reno, US: ASME, 2008: 583-593.
[9]	杨甫江,郭志辉,付虓. 基于热声网络法的燃烧不稳定性分析研究[J]. 推进技术,2014,35(6): 822-829. YANG Fujiang,GUO Zhihui,FU Xiao. Analytic study on combustion instability with a thermoacoustic network model[J]. Journal of Propulsion Technology,2014,35(6): 822-829. (in Chinese)
[10]	ZHANG Zhihao,LIU Xiao,GONG Yaozhen,et al. Investigation on flame characteristics of industrial gas turbine combustor with different mixing uniformities[J]. Fuel,2020,259: 116297.1-116297.13.
[11]	YANG Yang. Large eddy simulation calculated flame dynamics of one F class gas turbine combustor[J]. IOP Conference Series:Materials Science and Engineering,2019,643(1): 012141.1-012141.12.
[12]	KREDIET H J,BECK C H,KREBS W,et al. Saturation mechanism of the heat release response of a premixed swirl flame using LES[J]. Proceedings of the Combustion Institute,2013,34(1): 1223-1230. doi: 10.1016/j.proci.2012.06.140
[13]	CAMPA G, COSATTO E, CAMPOREALE S. Thermoacoustic analysis of combustion instability importing RANS data[R]. Milan, Italy: COMSOL Conference, 2012.
[14]	CAMPA G, CAMPOREALE S M, COSATTO E, et al. Thermoacoustic analysis of combustion instability through a distributed flame response function[C]//Proceedings of the ASME Turbo Expo. Copenhagen, Denmark: ASME, 2012: 179-188.
[15]	CAMPOREALE S M,FORTUNATO B,CAMPA G. A finite element method for three-dimensional analysis of thermo-acoustic combustion instability[J]. Journal of Engineering for Gas Turbines and Power,2011,133(1): 011506.1-011506.13.
[16]	CAMPA G, CAMPOREALE S M, GUAUS A, et al. A quantitative comparison between a low order model and a 3D FEM code for the study of thermoacoustic combustion instabilities[C]//Proceedings of ASME 2011 Turbo Expo: Turbine Technical Conference and Exposition. Vancouver, British Columbia, Canada: ASME, 2012: 859-869.
[17]	SHAW V, CLABBERS J, GUTMARK E. Augmentor combustion instability with COMSOL multiphysics[R]. Lausanne, Switzerland: COMSOL Conference, 2018
[18]	高贤智,何沛,冯晓星,等. 单扇区、扇形、全环燃烧室热声不稳定性试验和模拟研究[J]. 航空动力学报,2021,36(8): 1605-1613. GAO Xianzhi,HE Pei,FENG Xiaoxing,et al. Experimental and numerical investigation of thermo-acoustic instability in single-sector, sector and full annular combustors[J]. Journal of Aerospace Power,2021,36(8): 1605-1613. (in Chinese)
[19]	SILVA C F,NICOUD F,SCHULLER T,et al. Combining a Helmholtz solver with the flame describing function to assess combustion instability in a premixed swirled combustor[J]. Combustion and Flame,2013,160(9): 1743-1754. doi: 10.1016/j.combustflame.2013.03.020
[20]	NEZ R, FLESCA M G, MARCHAL D, et al. Experimental and numerical characterizations of acoustic damping rates in a coupled-cavity configuration[R]. Madrid, Spain: the 8th European Conference for Aeronautics and Aerospace Sciences, 2019.
[21]	TAMANAMPUDI G M R. Reduced order modeling and analysis of Combustion Instabilities[D]. West Lafayette, Indiana: Purdue University: 2015.
[22]	DOWLING A P. A kinematic model of a ducted flame[J]. Journal of Fluid Mechanics,1999,394: 51-72. doi: 10.1017/S0022112099005686
[23]	杨旸,余志健. 基于火焰传递函数的旋流预混火焰稳定性分析[J]. 燃气轮机技术,2021,34(1): 14-20. YANG Yang,YU Zhijian. Stability analysis of a swirling premixed combustor based on flame transfer functions[J]. Gas Turbine Technology,2021,34(1): 14-20. (in Chinese)
[24]	LIU Yunpeng,LI Jinghua,HAN Qixiang,et al. Study of combustion oscillation mechanism and flame image processing[J]. AIAA Journal,2019,57(2): 824-835. doi: 10.2514/1.J057614
[25]	吴江海, 吴健. 声固耦合有限元-无限元建模参数选取原则研究[R]. 昆明: 第18届船舶水下噪声学术讨论会论文集, 2021.
[26]	杨宇东,刘爱虢,王鑫慈,等. 燃油喷雾对某小发环形回流燃烧室燃烧特性影响[J]. 航空动力学报,2023,38(4): 830-839. YANG Yudong,LIU Aiguo,WANG Xinci,et al. Effect of fuel spray on combustion characteristics of a small engine annular reverse flow combustor[J]. Journal of Aerospace Power,2023,38(4): 830-839. (in Chinese)
[27]	BAE J,JEONG S,YOON Y. Effect of delay time on the combustion instability in a single-element combustor[J]. Acta Astronautica,2021,178: 783-792. doi: 10.1016/j.actaastro.2020.10.004