A simple predictive three-step chemical model for gaseous stoichiometric hydrogen-oxygen detonation quenching

被引:1
作者
Watanabe, Hiroaki [1 ]
Taileb, Said [2 ]
Veiga-Lopez, Fernando [3 ]
Melguizo-Gavilanes, Josue [1 ,4 ]
Chinnayya, Ashwin [1 ]
机构
[1] Univ Poitiers, Inst PPRIME, UPR 3346, CNRS,ISAE ENSMA, 1 Ave Clement Ader,BP 40109, F-86961 Furoscope Chasseneuil, France
[2] Safran Tech, Energy & Prop Dept, Rue Jeunes Bois, F-78114 Magny Les Hameaux, France
[3] Univ Vigo, Res Inst Phys & Aerosp Sci, Aerolab, Campus Ourense, Orense 32004, Galicia, Spain
[4] Shell Global Solut BV, Major Hazards Management, Energy Transit Campus, NL-1031 HW Amsterdam, Netherlands
关键词
Three-step chemical model; Detonation quenching; Numerical simulation; REDUCED CHEMISTRY; MECHANISM; WAVE; PROPAGATION; INSTABILITY; STABILITY; DYNAMICS; SUBSEQUENT; INITIATION; LIMITS;
D O I
10.1016/j.combustflame.2024.113609
中图分类号
O414.1 [热力学];
学科分类号
摘要
Two types of three-step chain-branching kinetics for gaseous hydrogen-oxygen detonation quenching are proposed. The reaction rates dependency on density and pressure, as well as the change in molecular weight are included and fitted to results of detailed chemistry not only along the Hugoniot curve but also behind transverse waves. The trend of induction and reaction times and reduced activation energy of zero-dimensional constant pressure combustion process as well as the induction, reaction lengths, and peak thermicity of the ZND model are focused on in the calibration procedure. One model (3SM-I) is designed to model the ideal detonation whereas the other (3SM-N) focuses on the D-n - kappa curve during calibration, i.e. the detonation velocity-mean front curvature curves, aiming to capture the non-idealities. The predictive ability using the proposed models is evaluated in a detonation transmission configuration, a layer of 2 H-2 2 - O-2 mixture being confined by N-2, in determining the critical height, reactive layer height at which detonation transmission is not possible. Both approaches get closer to the critical height from detailed chemistry. With the present set of model parameters, 3SM-I is better than 3SM-N for prediction of critical height. Moreover, borrowing analytical ignition criteria and wave hierarchy approach, the comparison of various chemical models that are able to recover the critical height suggests that the slope of the D-n - kappa curve, linked with the chemical response to losses, is the first key parameter for the critical height determination. In addition, the dynamics of the extinction near criticality such as appearance of transverse detonation and quenching distance from 3SM-I are in very good agreement with that of detailed chemistry due to better estimation of the chemical features along the Hugoniot curve and behind transverse waves. The present study provides a guideline to get simplified chemical models from detailed chemistry for predictive simulations. Novelty and significance statement The novelty of research is two-fold. The first is to provide general guideline to get simplified models from detailed chemistry for the simulation of detonation extinction with improved predictive ability. The tested configuration was that of a stoichiometric hydrogen-oxygen mixture semi-confined by nitrogen. The evolution of the induction, reaction times, and the maximum thermicity were used as fitting targets for our 3-step model. Furthermore, they were evaluated not only at the post-shock states along the Hugoniot curve, but also behind the transverse waves, keystones of the cellular structure. Second, borrowing analytical ignition criteria and wave hierarchy approach, the comparison of different chemical models that were able to recover the critical height suggested that the slope of the normal velocity-curvature curve was the first key parameter for the critical height determination. Moreover, the dynamics of extinction near criticality, and quenching distance depended on the reaction rate dependency on density and pressure.
引用
收藏
页数:13
相关论文
共 78 条
  • [1] Bane S.P.M., 2010, GALCIT REPORT FM2010.002
  • [2] BORISOV AA, 1991, PROGR ASTRONAUT AERO, V133, P118
  • [3] The evolution and cellular structure of a detonation subsequent to a head-on interaction with a shock wave
    Botros, Barbara B.
    Ng, Hoi Dick
    Zhu, YuJian
    Ju, Yiguang
    Lee, John H. S.
    [J]. COMBUSTION AND FLAME, 2007, 151 (04) : 573 - 580
  • [4] Continuous spin detonations
    Bykovskii, Fedor A.
    Zhdan, Sergey A.
    Vedernikov, Evgenii F.
    [J]. JOURNAL OF PROPULSION AND POWER, 2006, 22 (06) : 1204 - 1216
  • [5] Development of a virtual optimized chemistry method. Application to hydrocarbon/air combustion
    Cailler, Melody
    Darabiha, Nasser
    Fiorina, Benoit
    [J]. COMBUSTION AND FLAME, 2020, 211 : 281 - 302
  • [6] Effect of the reactor model on steady detonation modeling
    Chatelain, K. P.
    He, Y.
    Mevel, R.
    Lacoste, D. A.
    [J]. SHOCK WAVES, 2021, 31 (04) : 323 - 335
  • [7] Cheevers K, 2022, Arxiv, DOI arXiv:2211.05216
  • [8] Dabora E.K., 1965, 10 S INT COMBUSTION, P817, DOI DOI 10.1016/S0082-0784(65)80225-9
  • [9] Dabora E.K., 1963, Tech. Rep. 3559-E
  • [10] A parallel adaptive method for simulating shock-induced combustion with detailed chemical kinetics in complex domains
    Deiterding, Ralf
    [J]. COMPUTERS & STRUCTURES, 2009, 87 (11-12) : 769 - 783