基于单步与多步H<sub>2</sub>/Air机理的爆震波模拟网格尺度研究

李清安; 张永辉; 王可; 范玮; 陈梁霄; 杨海; 韩俊德; 康健

doi:10.13224/j.cnki.jasp.20230481

基于单步与多步H₂/Air机理的爆震波模拟网格尺度研究

doi: 10.13224/j.cnki.jasp.20230481

李清安^1,,
张永辉²,
王可^{2, 3,*, ,},
范玮^{2, 3},
陈梁霄¹,
杨海¹,
韩俊德¹,
康健¹

1.
中国航天科技集团有限公司四川航天系统工程研究所，成都 610100
2.
西北工业大学动力与能源学院，西安 710129
3.
西北工业大学陕西省航空动力系统热科学重点实验室，西安 710129

基金项目: 国家自然科学基金（52176133）；陕西省创新能力支撑计划（2021KJXX-93）

详细信息

作者简介:
李清安（1992-），男，工程师，博士，主要研究方向为燃烧流动与控制。E-mail：mcraecolin@163.com

通讯作者:
王可（1986-），男，教授，博士，主要研究方向为爆震推进理论与工程。E-mail：wangk@nwpu.edu.cn

中图分类号: V231.2
计量
- 文章访问数: 471
- HTML浏览量: 359
- PDF量: 37
- 被引次数: 0
出版历程
- 收稿日期: 2023-07-25
- 网络出版日期: 2024-08-29

Grid resolution study for detonation wave utilizing single-step and multi-step H₂/Air mechanisms

LI Qingan^1
,,
ZHANG Yonghui²,
WANG Ke^{2, 3,*
, ,},
FAN Wei^{2, 3},
CHEN Liangxiao¹,
YANG Hai¹,
HAN Junde¹,
KANG Jian¹

1.
Sichuan Aerospace System Engineering Research Institute，China Aerospace Science and Technology Corporation，Chengdu 610100，China
2.
School of Power and Energy，Northwestern Polytechnical University，Xi’an 710129，China
3.
Shaanxi Key Laboratory of Thermal Sciences in Aero-engine System，Northwestern Polytechnical University，Xi’an 710129，China

摘要

摘要:
为简化爆震推进系统仿真过程机理与网格的遴选工作，提高计算效率，分别采用单步与19步H₂/Air化学反应机理在长宽比为20的计算域中开展了爆震波的二维数值模拟，详细计算了0.5、0.4、0.2、0.1、0.08、0.04、0.02 mm共7种网格尺度下爆震波的传播过程，得到了爆震波胞格结构、爆震波特征参数等的变化规律。研究表明：采用单步或19步机理模拟得到的爆震波胞格结构的临界网格尺度分别为0.1 mm与0.4 mm，且网格尺度越小越容易获得爆震波胞格结构；爆震波胞格结构平均宽度与网格尺度呈线性关系；较大网格尺度会高估爆震波波峰面的静温与静压，但压力平台区的模拟结果几乎不受网格尺度与机理的影响；相较于多步机理，单步机理表现出的较短的点火延迟特性使模拟出的爆震波胞格结构边缘更细腻，使自持传播的爆震波更易获得。
- 爆震波 /
- 单步机理 /
- 胞格结构 /
- 网格无关性 /
- 爆震推进
Abstract:
In order to efficiently simplify the process of mechanism selection and resolution studies for detonation propulsion simulation, a series of detonation wave propagation simulations were conducted in a field with an aspect ratio of 20 using both the single-step and the multi-step H₂/Air mechanisms, all 7 sizes of grid were considered, including 0.5, 0.4, 0.2, 0.1, 0.08, 0.04, 0.02 mm. The cellular structure and the characteristic parameters of a detonation wave were obtained. The study showed that the cellular structure can be obtained when the grid size was smaller than a critical one, which was 0.1 mm and 0.4 mm for the single-step and multi-step mechanism, respectively. And the cellular structure was easier to be obtained when utilizing a smaller grid size. It was found that the averaged cell width was proportional to the grid size. Larger size of grid can enlarge the static pressure and temperature near the detonation peak, while the pressure plateau was weakly affected by both grid size and mechanism. Compared with the multi-step mechanism, the shorter ignition delay characteristics of the single-step mechanism contributed to a more clear edge of cellular structure, making it easier to obtain the self-sustained propagation of detonation wave.
- detonation wave /
- single-step mechanism /
- cellular structure /
- grid independence /
- detonation propulsion

HTML

图 1 计算域示意图

Figure 1. Schematic of the computation zone

下载: 全尺寸图片幻灯片

图 2 单步机理在0.5、0.4、0.2 mm网格尺度下的静压（上）与静温（下）云图

Figure 2. Static pressure （upper） and temperature （lower） contours using the single-step mechanism under0.5，0.4，0.2 mm grid settings

下载: 全尺寸图片幻灯片

图 3 多步机理在0.5、0.4、0.2 mm网格尺度下的静压（上）与静温（下）云图

Figure 3. Static pressure (upper) and temperature (lower) contours using the multi-step mechanism under 0.5，0.4，0.2 mm grid settings

下载: 全尺寸图片幻灯片

图 4 单步机理在0.1、0.08、0.04、0.02 mm网格尺度下的静压（上）与静温（下）云图

Figure 4. Static pressure （upper） and temperature （lower） contours using the single-step mechanism under0.1，0.08，0.04，0.02 mm grid settings

下载: 全尺寸图片幻灯片

图 5 多步机理在0.1、0.08、0.04、0.02 mm网格尺度下的静压（上）与静温（下）云图

Figure 5. Static pressure （upper） and temperature （lower） contours using the multi-step mechanism under0.1，0.08，0.04，0.02 mm grid settings

下载: 全尺寸图片幻灯片

图 6 单步机理在0.1、0.08、0.04、0.02 mm网格尺度下爆震波波锋附近的胞格结构

Figure 6. Cellular structure near the detonation wave peak using the single-step mechanism under 0.1，0.08，0.04，0.02 mm grid settings

下载: 全尺寸图片幻灯片

图 7 多步机理在0.4、0.2、0.1、0.08、0.04、0.02 mm网格尺度下爆震波波锋附近的胞格结构

Figure 7. Cellular structure near the detonation wave peak using the multi-step mechanism under 0.4，0.2，0.1，0.08，0.04，0.02 mm grid settings

下载: 全尺寸图片幻灯片

图 8 单步机理在0.1、0.08、0.04、0.02 mm网格尺度下的胞格结构

Figure 8. Cellular structure using the single-step mechanism under 0.1，0.08，0.04，0.02 mm grid settings

下载: 全尺寸图片幻灯片

图 9 多步机理在0.4、0.2、0.1、0.08、0.04、0.02 mm网格尺度下的胞格结构

Figure 9. Cellular structure using the multi-step mechanism under 0.4，0.2，0.1，0.08，0.04，0.02 mm grid settings

下载: 全尺寸图片幻灯片

图 10 单步机理在0.1、0.08、0.04、0.02 mm网格尺度下的胞格宽度统计结果

Figure 10. Cell width statistics with 0.1，0.08，0.04，0.02 mm grids settings using the single-step mechanism

下载: 全尺寸图片幻灯片

图 11 多步机理在0.4、0.2、0.1、0.08、0.04、0.02 mm网格尺度下的胞格宽度统计结果

Figure 11. Cell width statistics with 0.4，0.2，0.1，0.08，0.04，0.02 mm grids settings using the multi-step mechanism

下载: 全尺寸图片幻灯片

图 12 平均胞格宽度与网格尺度的关系

Figure 12. Averaged cell width versus grid settings

下载: 全尺寸图片幻灯片

图 13 单步机理不同网格尺度的特征静压沿x轴的分布

Figure 13. Characteristic static pressure along the x axis using the single-step mechanism for different grids settings

下载: 全尺寸图片幻灯片

图 14 多步机理不同网格尺度的特征静压沿x轴的分布

Figure 14. Characteristic static pressure along the x axis using the multi-step mechanism for different grids settings

下载: 全尺寸图片幻灯片

图 15 单步机理不同网格尺度的特征静温沿x轴的分布

Figure 15. Characteristic static temperature along the x axis using the single-step mechanism for different grids settings

下载: 全尺寸图片幻灯片

图 16 多步机理不同网格尺度的特征静温沿x轴的分布

Figure 16. Characteristic static temperature along the x axis using the multi-step mechanism

下载: 全尺寸图片幻灯片

图 17 单步机理不同网格尺度的爆震波波速

Figure 17. Propagation speed detonation wave using the single-step mechanism for different grids settings

下载: 全尺寸图片幻灯片

图 18 多步机理不同网格尺度的爆震波波速

Figure 18. Propagation speed detonation wave using the multi-step mechanism for different grids settings

下载: 全尺寸图片幻灯片

图 19 不同机理的点火延迟时间与起始温度的关系

Figure 19. Ignition delay versus initial temperature for different mechanisms under the same initial mixtures

下载: 全尺寸图片幻灯片

图 20 不同机理的点火延迟时间与起始压力的关系

Figure 20. Ignition delay versus initial pressure for different mechanisms under the same initial mixtures

下载: 全尺寸图片幻灯片

参考文献(29)

[1]	FICKETT W,DAVIS W C. Detonation: theory and experiment[M]. New York,US: Dover Publications Incorporation,2000.
[2]	严传俊,范玮. 脉冲爆震发动机原理及关键技术[M]. 西安: 西北工业大学出版社,2005. YAN Chuanjun,FAN Wei. Principles and key technologies of pulse detonation engines [M]. Xi’an: Northwestern Polytechnical University Press,2005. (in Chinese YAN Chuanjun, FAN Wei. Principles and key technologies of pulse detonation engines [M]. Xi’an: Northwestern Polytechnical University Press, 2005. (in Chinese)
[3]	ROY G D,FROLOV S M,BORISOV A A,et al. Pulse detonation propulsion: challenges,current status,and future perspective[J]. Progress in Energy and Combustion Science,2004,30(6): 545-672. doi: 10.1016/j.pecs.2004.05.001
[4]	LU F K,BRAUN E M. Rotating detonation wave propulsion: experimental challenges,modeling,and engine concepts[J]. Journal of Propulsion and Power,2014,30(5): 1125-1142. doi: 10.2514/1.B34802
[5]	WOLAŃSKI P. Detonative propulsion[J]. Proceedings of the Combustion Institute,2013,34(1): 125-158. doi: 10.1016/j.proci.2012.10.005
[6]	KUO K. K. Principles of combustion[M]. 2nd ed. Hoboken,US: John Wiley & Sons Incorporation,2005.
[7]	FAN Weijie,LIU Weidong,PENG Haoyang,et al. Numerical study on ethylene-air continuous rotating detonation in annular combustors with different widths[J]. Journal of Zhejiang University-SCIENCE A,2022,23(5): 388-404. doi: 10.1631/jzus.A2100448
[8]	潘鑫峰,雷雨冰. 爆震波在开缝式弯管内的传播特性研究[J]. 推进技术,2019,40(1): 123-129. PAN Xinfeng,LEI Yubing. Study on propagation characteristics of detonation wave in slotted elbow[J]. Journal of Propulsion Technology,2019,40(1): 123-129. (in Chinese PAN Xinfeng, LEI Yubing. Study on propagation characteristics of detonation wave in slotted elbow[J]. Journal of Propulsion Technology, 2019, 40(1): 123-129. (in Chinese)
[9]	WANG Kuanliang,TENG Honghui,YANG Pengfei,et al. Numerical investigation of flow structures resulting from the interaction between an oblique detonation wave and an upper expansion corner[J]. Journal of Fluid Mechanics,2020,903: A28.1-A28.17.
[10]	孙志鹏,黄玥,栾振业,等. 高温高压进口下甲烷/氧气旋转爆震燃烧室增压特性与流场分析[J]. 推进技术,2022,43(11): 210751. SUN Zhipeng,HUANG Yue,LUAN Zhenye,et al. Pressurization characteristics and flow field analysis of methane/oxygen rotating detonation combustor at high temperature and high pressure inlets[J]. Journal of Propulsion Technology,2022,43(11): 210751. (in Chinese SUN Zhipeng, HUANG Yue, LUAN Zhenye, et al. Pressurization characteristics and flow field analysis of methane/oxygen rotating detonation combustor at high temperature and high pressure inlets[J]. Journal of Propulsion Technology, 2022, 43(11): 210751. (in Chinese)
[11]	覃建秀,朱德华. 受限空间中斜爆震波波系结构及其演化 [J]. 推进技术,2024,45(2):2208039. QIN Jianxiu,ZHU Dehua. Wave structure and structural evolution of oblique detonation waves in confined space [J]. Journal of Propulsion Technology,2024,45(2):2208039. (in Chinese QIN Jianxiu, ZHU Dehua. Wave structure and structural evolution of oblique detonation waves in confined space [J]. Journal of Propulsion Technology, 2024, 45（2）: 2208039. (in Chinese)
[12]	蔡晓东,梁剑寒,林志勇,等. 基于自适应网格加密的超声速可燃气热射流起爆详细反应数值模拟[J]. 航空动力学报,2014,29(10): 2385-2392. CAI Xiaodong,LIANG Jianhan,LIN Zhiyong,et al. Adaptive mesh refinement-based numerical simulation of initiation in supersonic combustible mixtures using hot jet with detailed reaction model[J]. Journal of Aerospace Power,2014,29(10): 2385-2392. (in Chinese CAI Xiaodong, LIANG Jianhan, LIN Zhiyong, et al. Adaptive mesh refinement-based numerical simulation of initiation in supersonic combustible mixtures using hot jet with detailed reaction model[J]. Journal of Aerospace Power, 2014, 29(10): 2385-2392. (in Chinese)
[13]	焦中天,王永佳,李伟,等. 燃料喷孔数对非预混旋转爆震起爆过程的影响[J]. 火箭推进,2021,47(5): 22-34. JIAO Zhongtian,WANG Yongjia,LI Wei,et al. Effects of the number of fuel injection orifices on rotating detonation initiation process under non-premixed conditions[J]. Journal of Rocket Propulsion,2021,47(5): 22-34. (in Chinese JIAO Zhongtian, WANG Yongjia, LI Wei, et al. Effects of the number of fuel injection orifices on rotating detonation initiation process under non-premixed conditions[J]. Journal of Rocket Propulsion, 2021, 47(5): 22-34. (in Chinese)
[14]	李清安,王可,孙田雨,等. 脉冲爆震发动机反压传播规律数值研究[J]. 实验流体力学,2019,33(1): 103-110. LI Qingan,WANG Ke,SUN Tianyu,et al. Numerical study on propagation characteristics of back-pressure in a pulse detonation engine[J]. Journal of Experiments in Fluid Mechanics,2019,33(1): 103-110. (in Chinese LI Qingan, WANG Ke, SUN Tianyu, et al. Numerical study on propagation characteristics of back-pressure in a pulse detonation engine[J]. Journal of Experiments in Fluid Mechanics, 2019, 33(1): 103-110. (in Chinese)
[15]	张颂睿,曾昊,何立明,等. 新型激波聚焦脉冲爆震模型连续起爆探究[J]. 空军工程大学学报(自然科学版),2021,22(6): 42-48. ZHANG Songrui,ZENG Hao,HE Liming,et al. A study of continuous initiation for a new shockwave focus pulse detonation model[J]. Journal of Air Force Engineering University (Natural Science Edition),2021,22(6): 42-48. (in Chinese ZHANG Songrui, ZENG Hao, HE Liming, et al. A study of continuous initiation for a new shockwave focus pulse detonation model[J]. Journal of Air Force Engineering University (Natural Science Edition), 2021, 22(6): 42-48. (in Chinese)
[16]	黄希桥,王丁丁,肖治邑,等. 航空煤油单组分替代燃料DDT过程数值模拟[J]. 西北工业大学学报,2022,40(3): 610-617. HUANG Xiqiao,WANG Dingding,XIAO Zhiyi,et al. Numerical simulation of DDT process of aviation kerosene one-component alternative fuel[J]. Journal of Northwestern Polytechnical University,2022,40(3): 610-617. (in Chinese doi: 10.1051/jnwpu/20224030610 HUANG Xiqiao, WANG Dingding, XIAO Zhiyi, et al. Numerical simulation of DDT process of aviation kerosene one-component alternative fuel[J]. Journal of Northwestern Polytechnical University, 2022, 40(3): 610-617. (in Chinese) doi: 10.1051/jnwpu/20224030610
[17]	沈吴冰懿,黄玥,尤延铖,等. 圆盘旋转爆震燃烧室中爆震波传播模态[J]. 推进技术,2021,42(4): 865-873. SHEN Wubingyi,HUANG Yue,YOU Yancheng,et al. Propagation mode of detonation wave in a plane-radial rotating detonation combustion chamber[J]. Journal of Propulsion Technology,2021,42(4): 865-873. (in Chinese SHEN Wubingyi, HUANG Yue, YOU Yancheng, et al. Propagation mode of detonation wave in a plane-radial rotating detonation combustion chamber[J]. Journal of Propulsion Technology, 2021, 42(4): 865-873. (in Chinese)
[18]	TENG Hong,JIANG Zong,NG H D. Numerical study on unstable surfaces of oblique detonations[J]. Journal of Fluid Mechanics,2014,744: 111-128. doi: 10.1017/jfm.2014.78
[19]	刘朋欣,郭启龙,赵炜,等. 基于旋转爆震三维流场结构分析的计算模型对比研究[J]. 推进技术,2020,41(12): 2757-2765. LIU Pengxin,GUO Qilong,ZHAO Wei,et al. Computational models based on analysis of three dimensional flow field structures in rotating detonation[J]. Journal of Propulsion Technology,2020,41(12): 2757-2765. (in Chinese LIU Pengxin, GUO Qilong, ZHAO Wei, et al. Computational models based on analysis of three dimensional flow field structures in rotating detonation[J]. Journal of Propulsion Technology, 2020, 41(12): 2757-2765. (in Chinese)
[20]	BURKE M P,CHAOS M,JU Y,et al. Comprehensive H₂/O₂ kinetic model for high-pressure combustion[J]. International Journal of Chemical Kinetics,2012,44: 444-474. doi: 10.1002/kin.20603
[21]	MALCOLM W,CHASE J. Journal of physical and chemical reference data: NIST-JANAF thermochemical tables [M]. 4th ed. New York,US: American Chemical Society and the Americal Institute of Physics for the National Institute of Standards and Technology,1998.
[22]	LEE J H S. The detonation phenomenon[M]. Cambridge,UK: Cambridge University Press,2008.
[23]	TAYLOR B D,KESSLER D A,GAMEZO V N,et al. Numerical simulations of hydrogen detonations with detailed chemical kinetics[J]. Proceedings of the Combustion Institute,2013,34(2): 2009-2016. doi: 10.1016/j.proci.2012.05.045
[24]	KANESHIGE M,SHEPHERD J E. Detonation database: GALCIT. FM97-8 [R]. Pasadena,US: Explosion Dynamics Laboratory,1999.
[25]	ENDO T,KASAHARA J,MATSUO A,et al. Pressure history at the thrust wall of a simplified pulse detonation engine[J]. AIAA Journal,2004,42(9): 1921-1930. doi: 10.2514/1.976
[26]	LI Qingan,WANG Ke,WANG Zhicheng,et al. Evaluation of the two-γ model and the energy balance process of a Chapman-Jouguet detonation[J]. Aerospace Science and Technology,2021,117: 106955. doi: 10.1016/j.ast.2021.106955
[27]	SCHULTZ E,SHEPHERD J. Validation of detailed reaction mechanisms for detonation simulation: FM99-5 [R]. Pasadena,US: California Institute of Technology,2000.
[28]	CONAIRE M Ó,CURRAN H J,SIMMIE J M,et al. A comprehensive modeling study of hydrogen oxidation[J]. International Journal of Chemical Kinetics,2004,36(11): 603-622. doi: 10.1002/kin.20036
[29]	MÉVEL R,JAVOY S,LAFOSSE F,et al. Hydrogen-nitrous oxide delay times: shock tube experimental study and kinetic modelling[J]. Proceedings of the Combustion Institute,2009,32(1): 359-366. doi: 10.1016/j.proci.2008.06.171