Detonation propagation in weakly confined gases

TL;DR

This paper combines CFD simulations and analytical models to study how detonations propagate in weakly confined layered gases, revealing different regimes and transition mechanisms relevant to rotating detonation engines.

Contribution

It introduces a unified numerical-theoretical framework for layered detonation dynamics, including criteria for overdriven and underdriven regimes based on impedance and layer thickness.

Findings

Identification of flowfield regimes depending on impedance and layer thickness

Development of an analytical phase map for detonation behavior

Validation of models with CFD results

Abstract

This study investigates the propagation of detonations along a layered configuration where a reactive gas is weakly confined by a hotter inert layer. CFD simulations are performed using a single-step, non-Arrhenius reaction model designed to suppress cellular instabilities, enabling formulation of a theoretical framework directly compared with simulation results. The simulations reach a quasi-steady state, revealing distinct flowfield regimes that depend on the acoustic-impedance ratio and relative layer thicknesses, with some detonations exhibiting velocity deficits while others propagate above the ideal Chapman-Jouguet (CJ) speed. Analytical models are developed to interpret these regimes. When a precursor shock is observed in the inert layer, the detonation is overdriven; this is modeled using shock-polar analysis and velocity estimates based on the approach of Mitrofanov (Acta Astronaut. 3:995-1004, 1976). An analytical criterion for precursor shock onset is proposed. In underdriven scenarios, the detonation front exhibits positive curvature, analyzed using a geometric construction wherein the relationship between wave speed and front curvature is evaluated a priori. A simplified characteristic-based model captures the decay of the shock wave in the inert layer, after which shock-polar analysis determines the resulting wave interaction. Predictions from these models are assembled into a phase map delineating regions of overdriven and underdriven behavior, along with corresponding shock interactions, in the space of acoustic impedance and area ratios. This map is compared directly with CFD results. The combined numerical-theoretical framework clarifies transition mechanisms governing layered detonations and provides insights into detonation dynamics relevant to rotating detonation engines in which the detonation is bounded by hotter combustion products from a previous cycle.