AECSC-IBM航空发动机燃烧室数值模拟软件研发与检验

王煜栋; 王方; 周佳伟; 金捷

doi:10.13224/j.cnki.jasp.20220216

AECSC-IBM航空发动机燃烧室数值模拟软件研发与检验

doi: 10.13224/j.cnki.jasp.20220216

王煜栋^1,,
王方^{1, 2, 3,*, ,},
周佳伟¹,
金捷^{1, 2, 3}

1.
北京航空航天大学能源与动力工程学院航空发动机数值仿真研究中心，北京 100191
2.
北京航空航天大学成都航空动力创新研究院，成都 611930
3.
北京航空航天大学江西研究院，南昌 330096

基金项目: 国家自然科学基金（91741125）; 国家科技重大专项（2017-Ⅰ-0004-0005)

详细信息

作者简介:
王煜栋（1998-），男，硕士生，主要从事两相湍流燃烧理论及其应用研究。E-mail: wangyudong@buaa.edu.cn

通讯作者:
王方（1972-），女，副教授、硕士生导师，博士，主要从事两相湍流燃烧理论及其应用研究。E-mail: fwang@buaa.edu.cn

中图分类号: V231.2
计量
- 文章访问数: 1128
- HTML浏览量: 325
- PDF量: 121
- 被引次数: 0
出版历程
- 收稿日期: 2022-04-15
- 网络出版日期: 2022-09-14

Development and inspection of aero-engine combustor simulation software AECSC-IBM

WANG Yudong^1
,,
WANG Fang^{1, 2, 3,*
, ,},
ZHOU Jiawei¹,
JIN Jie^{1, 2, 3}

1.
Aeroengine Numerical Simulation Research Center，School of Energy and Power Engineering，Beihang University，Beijing 100191，China
2.
Chengdu Innovation Research Institute on Aircraft Power，Beihang University，Chengdu 611930，China
3.
Jiangxi Research Institute, Beihang University，Nanchang 330096，China

摘要

摘要:
针对航空发动机燃烧室的高保真数值模拟需求，基于浸没边界方法（IBM）及大涡模拟-输运概率密度函数湍流燃烧模型（LES-TPDF）开发软件AECSC-IBM，用网格标记映射燃烧室真实几何结构。通过模拟双旋流燃烧室算例和Sandia射流火焰算例检验湍流流动和燃烧的模拟精度。在双旋流燃烧室模拟中，旋流器出口时均轴向、径向、切向速度平均误差分别为15.7%、23.8%和15.0%。在射流火焰的模拟中，Flame-E和Flame-F的温度、燃料质量分数平均相对误差分别为14.69%、5.22%和14.18%、5.54%。进一步将AECSC-IBM软件应用于某真实结构单头部燃烧室算例，模拟得到出口温度与实验数据相比方均根误差为11.66%。算例检验表明AECSC-IBM软件能快速精确映射几何模型，大幅减少复杂几何高质量网格生成工作量，高效准确地模拟航空发动机燃烧室内的两相湍流燃烧现象，模拟结果可为燃烧室精细化研发提供燃烧场数据参考，具有工程实用价值。
- AECSC-IBM软件 /
- 大涡模拟(LES) /
- 输运概率密度函数方程(TPDF) /
- 浸没边界方法(IBM) /
- 燃烧室数值模拟
Abstract:
To meet the needs of high-fidelity numerical simulation of aero-engine combustors, the AECSC-IBM software was developed based on the immersion boundary method (IBM) and the large-eddy simulation-transported probability density function combustion model (LES-TPDF). The original geometric structure of the combustor was mapped with grid markers. The simulation accuracy of turbulent flow and combustion was verified by a simulation example of a twin-swirl combustor and Sandia jet flame. In the simulation of twin cyclone combustors, the average errors of axial,radial and tangential velocities at the outlet of the swirler were 15.7%, 23.8%, and 15.0%, respectively. The average relative errors of temperature and fuel mass fraction for flame-E, flame-F were 14.69% and 5.22%, 14.18% and 5.54%, respectively. Furthermore, AECSC-IBM software was applied to a single head combustor of a real structure, and the root mean square error of the simulated exit temperature was 11.66% compared with the experimental data. Example tests showed that AECSC-IBM software can map geometric models quickly and accurately, reduce the workload of complex geometry high-quality mesh generation greatly, and simulate the two-phase turbulent combustion phenomenon in the aero-engine combustion chamber efficiently and accurately. The simulation results can provide a reference for combustion field data in combustion chamber refinement research and development,presenting practical engineering value.
- AECSC-IBM software /
- large-eddy simulation (LES) /
- transported probability density function (TPDF) /
- immersed boundary method (IBM) /
- combustor numerical simulation

HTML

图 1 GTMC几何结构示意图^[19]

Figure 1. Schematic diagram of the geometric structure of the GTMC^[19]

下载: 全尺寸图片幻灯片

图 2 GTMC算例网格与IBM标记

Figure 2. Grid and IBM mark for the GTMC

下载: 全尺寸图片幻灯片

图 3 GTMC算例中央剖面瞬态速度矢量图

Figure 3. Central section transient velocity vector diagram of the GTMC

下载: 全尺寸图片幻灯片

图 4 GTMC算例三维瞬态流线图

Figure 4. 3D transient streamline diagram of GTMC

下载: 全尺寸图片幻灯片

图 5 AECSC-IBM模拟结果与商业软件FLUENT模拟结果中央界面时均轴向速度图

Figure 5. Time-averaged axial velocity diagram of the center interface by AECSC-IBM simulation and commercial software FLUENT simulation

下载: 全尺寸图片幻灯片

图 6 旋流器出口时均轴向速度分布

Figure 6. Time averaged axial velocity distribution at the outlet of the swirler

下载: 全尺寸图片幻灯片

图 7 旋流器出口时均切向速度分布

Figure 7. Time averaged tangential velocity distribution at the outlet of the swirler

下载: 全尺寸图片幻灯片

图 8 旋流器出口时均径向速度分布

Figure 8. Time averaged radial velocity distribution at the outlet of the swirler

下载: 全尺寸图片幻灯片

图 9 射流火焰网格图

Figure 9. Grid diagram of jet flame

下载: 全尺寸图片幻灯片

图 10 射流火焰瞬态速度云图

Figure 10. Transient velocity diagram of jet flame

下载: 全尺寸图片幻灯片

图 11 射流火焰湍流动能云图

Figure 11. Turbulent kinetic energy diagram of jet flame

下载: 全尺寸图片幻灯片

图 12 射流火焰瞬态温度云图

Figure 12. Transient temperature diagram of jet flame

下载: 全尺寸图片幻灯片

图 13 射流火焰瞬态甲烷质量分数云图

Figure 13. Transient methane mass fraction diagram of jet flame

下载: 全尺寸图片幻灯片

图 14 射流火焰时均温度沿径向的分布

Figure 14. Time-averaged temperature distribution of jet flame along the radial direction

下载: 全尺寸图片幻灯片

图 15 射流火焰时均甲烷质量分数沿径向的分布

Figure 15. Time-averaged methane mass fraction distribution of jet flame along the radial direction

下载: 全尺寸图片幻灯片

图 16 单头部燃烧室几何模型

Figure 16. Geometry of the single-head combustor

下载: 全尺寸图片幻灯片

图 17 单头部燃烧室三角面网格

Figure 17. Triangular mesh of the single-head combustor

下载: 全尺寸图片幻灯片

图 18 计算域背景网格与IBM网格标记

Figure 18. Computational domain background grid and IBM grid mark

下载: 全尺寸图片幻灯片

图 19 三维IBM网格标记

Figure 19. 3D IBM grid mark

下载: 全尺寸图片幻灯片

图 20 中央截面瞬态速度矢量图

Figure 20. Central section transient velocity vector diagram

下载: 全尺寸图片幻灯片

图 21 三维流线图

Figure 21. 3D streamline diagram

下载: 全尺寸图片幻灯片

图 22 旋流器轴向时均速度与回流区

Figure 22. Time-averaged axial velocity and recirculation zone of the swirler

下载: 全尺寸图片幻灯片

图 23 燃烧室中央剖面三维瞬态温度云图

Figure 23. 3D transient temperature diagram of the central section of the combustor

下载: 全尺寸图片幻灯片

图 24 燃烧室内时均气相煤油质量分数分布

Figure 24. Time-averaged kerosene mass fraction distribution in the combustor

下载: 全尺寸图片幻灯片

图 25 燃烧室中央剖面三维时均温度分布

Figure 25. 3D time-averaged temperature distribution on the central section of the combustor

下载: 全尺寸图片幻灯片

图 26 燃烧室三维时均高温区分布

Figure 26. 3D time-averaged high temperature zone distribution in the combustor

下载: 全尺寸图片幻灯片

图 27 燃烧室出口测点位置径向温度分布

Figure 27. Radial temperature distribution of the measuring points at the exit of the combustor

下载: 全尺寸图片幻灯片

表 1 GTMC时均速度平均相对误差

Table 1. Average relative error of GTMC time-averaged velocity

高度/mm	速度方向	平均相对误差/%
高度/mm	速度方向	AECSC-IBM	FLUENT
2	轴向	13.8	15.1
5	轴向	12.0	11.4
10	轴向	21.2	15.8
2	径向	36.6	34.5
5	径向	16.9	11.9
10	径向	18.0	10.2
2	切向	15.4	11.2
5	切向	11.6	11.2
10	切向	18.0	14.3

下载: 导出CSV

表 2 各横截面时均径向温度分布平均相对误差

Table 2. Time-averaged radial temperature distribution relative error for each cross section

算例	$ z/d $	平均相对误差/%	算例	$z/d $	平均相对误差/%
Flame-E	1	13.30	Flame-F	1	19.57
Flame-E	2	11.51	Flame-F	2	21.23
Flame-E	3	12.77	Flame-F	3	18.60
Flame-E	45	22.38	Flame-F	45	7.16
Flame-E	60	13.79	Flame-F	60	8.61
Flame-E	75	14.37	Flame-F	75	9.91

下载: 导出CSV

表 3 各横截面时均径向甲烷质量分数分布平均相对误差

Table 3. Time-averaged radial methane mass fraction distribution relative error each cross section

算例	$ z/d $	平均相对误差/%	算例	$z/d $	平均相对误差/%
Flame-E	1.0	4.98	Flame-F	1.0	4.44
Flame-E	2.0	4.43	Flame-F	2.0	5.14
Flame-E	3.0	3.84	Flame-F	3.0	3.93
Flame-E	7.5	6.40	Flame-F	7.5	4.96
Flame-E	15.0	6.43	Flame-F	15.0	9.25

下载: 导出CSV

参考文献(27)

[1]	尚守堂,林宏军,程明,等. 航空发动机燃烧室数值仿真技术工程应用分析[J]. 航空动力,2021,19(2): 66-69. SHANG Shoutang,LIN Hongjun,CHENG Ming,et al. Engineering applications of numerical simulation technology for aero engine combustor[J]. Aerospace Power,2021,19(2): 66-69. (in Chinese)
[2]	王方,王煜栋,姜胜利,等. AECSC-JASMIN湍流燃烧仿真软件研发和检验[J]. 航空学报,2021,42(12): 128-140. WANG Fang,WANG Yudong,JIANG Shengli,et al. Development and inspection of AECSC-JASMIN turbulent combustion simulation software[J]. Acta Aeronautica et Astronautica Sinica,2021,42(12): 128-140. (in Chinese)
[3]	JONES W P,MARQUIS A J,VOGIATZAKI K. Large-eddy simulation of spray combustion in a gas turbine combustor[J]. Combustion and Flame,2014,161(1): 222-239.
[4]	李威,张哲巅. 某航空发动机燃烧室天然气湿燃烧数值模拟研究[J]. 热能动力工程,2021,36(5): 126-133. doi: 10.16146/j.cnki.rndlgc.2021.05.019 LI Wei,ZHANG Zhedian. Numerical investigation on a jet-engine combustor fueled by humidified natural gas[J]. Journal of Engineering for Thermal Energy and Power,2021,36(5): 126-133. (in Chinese) doi: 10.16146/j.cnki.rndlgc.2021.05.019
[5]	莫妲,程明,万斌,等. 三旋流燃烧室的数值模拟与实验[J]. 航空动力学报,2017,32(11): 2568-2575. MO Da,CHENG Ming,WAN Bin,et al. Numerical simulation and experiment of tri-swirl combustor[J]. Journal of Aerospace Power,2017,32(11): 2568-2575. (in Chinese)
[6]	丁勇能,田勇,王波,等. 某重型燃气轮机DLN燃烧室数值模拟-热态分析[J]. 燃气轮机技术,2019,32(2): 40-43. doi: 10.16120/j.cnki.issn1009-2889.2019.02.007 DING Yongneng,TIAN Yong,WANG Bo,et al. Numerical simulation of DLN combustor for a heavy-duty gas turbine[J]. Gas Turbine Technology,2019,32(2): 40-43. (in Chinese) doi: 10.16120/j.cnki.issn1009-2889.2019.02.007
[7]	ZHOU Y,LE J L,HUANG Y. LES of combustion flow field in a practical aeroengine combustor with two-stage counter-rotating swirler[J]. Journal of Propulsion Technology,2018,39(7): 141-154.
[8]	阎超,于剑,徐晶磊,等. CFD模拟方法的发展成就与展望[J]. 力学进展,2011,41(5): 562-589. doi: 10.6052/1000-0992-2011-5-lxjzJ2010-082 YAN Chao,YU Jian,XU Jinglei,et al. On the achievements and prospects for the method of computational fluid dynamics[J]. Advances in Mechanics,2011,41(5): 562-589. (in Chinese) doi: 10.6052/1000-0992-2011-5-lxjzJ2010-082
[9]	PESKIN C S. Flow patterns around heart values:a numerical method[J]. Journal of Computational Physics,1972,10(2): 252-271. doi: 10.1016/0021-9991(72)90065-4
[10]	ROMA A M.A multilevel self-adaptive version of the immersed boundary method[D].New York:New York University,1996.
[11]	LEVEQUE R J,LI Z. Immersed interface method for stokes flow with elastic boundaries or surface tension[J]. SIAM Journal on Scientific Computing,1997,18(3): 709-735. doi: 10.1137/S1064827595282532
[12]	BEYER R P,LEVEQUE R J. Analysis of a one-dimensional model for the immersed boundaries method[J]. SIAM Journal on Numerical Analysis,1992,29(2): 332-364. doi: 10.1137/0729022
[13]	张漫,王铮钧,王晶,等. 航空发动机内流全场流动的大涡模拟[J]. 航空动力,2021,19(2): 59-60. ZHANG Man,WANG Zhengjun,WANG Jing,et al. Large eddy simulation on internal flow of aero engine[J]. Aerospace Power,2021,19(2): 59-60. (in Chinese)
[14]	金捷,刘邓欢. 航空发动机燃烧室湍流两相燃烧模型发展现状[J]. 南京航空航天大学学报,2016,48(3): 303-309. doi: 10.16356/j.1005-2615.2016.03.002 JIN Jie,LIU Denghuan. Recent advances in turbulent two-phase combustion models[J]. Journal of Nanjing University of Aeronautics and Astronautics,2016,48(3): 303-309. (in Chinese) doi: 10.16356/j.1005-2615.2016.03.002
[15]	JONES W P,PRASAD V N. Large eddy simulation of the sandia flame series (D-F) using the Eulerian stochastic field method[J]. Combustion and Flame,2010,157(9): 1621-1636. doi: 10.1016/j.combustflame.2010.05.010
[16]	JONES W P,MARQUIS A J,WANG F. Large eddy simulation of a premixed propane turbulent bluff body flame using the Eulerian stochastic field method[J]. Fuel,2015,140: 514-525. doi: 10.1016/j.fuel.2014.06.050
[17]	WANG F,LIU R,DOU L,et al. A dual timescale model for micromixing and its application in LES/TPDF simulations of turbulent nonpremixed flames[J]. Chinese Journal of Aeronautics,2019,32(4): 52-64.
[18]	LUTZ A E,KEE R J,GRCAR J F,et al.OPPDIF:a FORTRAN program for computing opposed flow diffusion flames[R].Livermore,US:Sandia National Lab,1997.
[19]	MEIER U,HEINZE J,FREITAG S,et al. Spray and flame structure of a generic injector at aeroengine conditions[J]. Journal of Engineering for Gas Turbines and Power,2012,134(3): 031503.1-031503.9.
[20]	JABERI F A,COLUCCI P J,JAMES S,et al. Filtered mass density function for large-eddy simulation of turbulent reacting flows[J]. Journal of Fluid Mechanics,1999,401(2): 85-121.
[21]	LAGAE A,DUTRE P. An efficient ray-quadrilateral intersection test[J]. Journal of Graphics Tools,2005,10(4): 23-32. doi: 10.1080/2151237X.2005.10129208
[22]	STOKIE J M,WETTON B T R. Stability analysis for the immersed fiber problem[J]. SIAM Journal on Applied Mathematics,1995,55(6): 1577-1591. doi: 10.1137/S0036139994267018
[23]	STOKIE J M,WETTON B R. Analysis of stiffness in the immersed boundary method and implications for time-stepping schemes[J]. Journal of Computational Physics,1999,154(1): 41-64. doi: 10.1006/jcph.1999.6297
[24]	曾家,金捷,张晟,等. 基于LES-PDF方法的双旋流模型燃烧室数值模拟[J]. 气体物理,2019,4(5): 52-64. doi: 10.19527/j.cnki.2096-1642.0789 ZENG Jia,JIN Jie,ZHANG Sheng,et al. Numerical simulation of double-swirled model combustor based on LES-PDF[J]. Physics of Gases,2019,4(5): 52-64. (in Chinese) doi: 10.19527/j.cnki.2096-1642.0789
[25]	BARLOW R S,FRANK J H.Effects of turbulence on species mass fractions in methane/air jet flames[J].Symposium (International) on Combustion,1998,27(1): 1087-1095.
[26]	KUNDU K P,DEUR J M.A simplified reaction mechanism for calculation of emissions in hydrocarbon (Jet-A) combustion[R].AIAA 93-2341,1993.
[27]	WANG J,MAO R,TAO W,et al. Numerical simulation of combustor effusion cooling flow based on source term method[R].ASME Paper GT2019-90885,2019.