Numerical simulations of flame propagation and DDT in obstructed channels filled with hydrogen–air mixture

We study flame acceleration and deflagration-to-detonation transition (DDT) in channels with obstacles using 2D and 3D reactive Navier–Stokes numerical simulations. The energy release rate for the stoichiometric H2–air mixture is modeled by a one-step Arrhenius kinetics. Computations show that at initial stages, the flame and flow acceleration is caused by thermal expansion of hot combustion products. At later stages, shock–flame interactions, Rayleigh–Taylor, Richtmyer–Meshkov, and Kelvin–Helmholtz instabilities, and flame–vortex interactions in obstacle wakes become responsible for the increase of the flame surface area, the energy-release rate, and, eventually, the shock strength. Computations performed for different channel widths d with the distance between obstacles d and the constant blockage ratio 0.5 reproduce the main regimes observed in experiments: choking flames, quasi-detonations, and detonations. For quasi-detonations, both the initial DDT and succeeding detonation reignitions occur when the Mach stem, created by the reflection of the leading shock from the bottom wall, collides with an obstacle. As the size of the system increases, the time to DDT and the distance to DDT increase linearly with d2. We also observe an intermediate regime of fast flame propagation in which local detonations periodically appear behind the leading shock, but do not reach it.