Computation of the mean hydrodynamic structure of gaseous detonations with losses

Lagrangian dispersion and averaging behind a two-dimensional gaseous detonation front

“…Transverse disturbances then came into play in 2-D configurations. As argued at first by Radulescu et al (2007) and by many others (Maxwell et al 2017;Taileb et al 2018;Reynaud et al 2020;Sow et al 2021), the fluctuations and the induced dispersion explain the differences between the mean quantities from x i at t 0 (mm) numerical simulations and the ZND results. All dispersion quantities (x 2 i , y 2 i ), when non-dimensionalized by (E a /(RT vN ))x ind were found to be self-similar in the time τ/τ ind .…”

“…The grid resolution is approximately 38 and 30 points per CJ induction length for the 2H 2 -O 2 -7Ar mixture and the 2H 2 -O 2 mixture, respectively. This resolution has been shown to be largely sufficient to capture the mean structure (Reynaud et al 2017(Reynaud et al , 2020. In addition, this resolution is enough to reproduce the features of the instantaneous flow fields for weakly unstable mixtures (Mazaheri, Mahmoudi & Radulescu 2012).…”

“…When the leading shock front reached the right-hand boundary during the simulations, the new Lagrangian particles were located to the right of the computational domain and the Lagrangian particles which were located in the discarded left-hand domain were excluded from the simulations. The recycling block technique was successfully utilized to reduce the computational cost by the use of smaller computational domain and to simulate the detonation propagation in the previous studies (Reynaud et al 2017(Reynaud et al , 2020Sow et al 2019;Watanabe et al 2020Watanabe et al , 2021Taileb, Meluguizo-Gavilances & Chinnayya 2021a;Taileb et al 2021b). The length of the propagation for the average procedure is approximately 1000 x ind for 2H 2 -O 2 -7Ar and 1200 x ind for 2H 2 -O 2 .…”

“…Furthermore, Sow, Chinnayya & Hadjadj (2014) proposed the Favre average procedure for the detonation in the non-inertial instantaneous shock frame to take into account the unsteadiness of the shock front. So far, the Favre average procedure to obtain 1-D profiles was applied to planar detonations Radulescu et al 2007;Maxwell et al 2017;Taileb et al 2018;Sow, Lau-Chapdelaine & Radulescu 2021;Taileb, Meluguizo-Gavilances & Chinnayya 2021b), in non-uniform mixtures (Mi, Timofeev & Higgins 2017a;Mi et al 2017b), in mixtures with concentration gradients (Han, Wang & Law 2019), in mixtures with fluctuations in concentrations (Zhou et al 2022), cylindrical detonation (Han et al 2017), also in non-ideal configurations such as detonations bounded by an inert layer (Reynaud, Virot & Chinnayya 2017;Reynaud et al 2020), with wall losses (Chinnayya, Hadjadj & Ngomo 2013;Sow et al 2014;Sow, Chinnayya & Hadjadj 2015, in two-phase detonations with water spray (Watanabe et al 2019(Watanabe et al , 2020(Watanabe et al , 2021 and with fuel spray (Jourdaine, Tsuboi & Hayashi 2022).…”

The influence of the equation of state on the cellular structure of gaseous detonations

“…The governing equations were integrated using our parallel in-house code RESIDENT (REcycling mesh SImulations of DEtoNaTions) [34][35][36] . Details about the numerical methods used for spatial and temporal discretizations as well as the parallelization methodology can 7…”