Numerical Study of Detonation Wave Di ↵ raction in Hydrogen-based Mixtures
Stany Gallier,Rémy Mével,Dmitry Davidenko,Florian Pintgen,Joseph Shepherd
2013-01-01
Abstract:As a detonation wave propagates from a confined tube to an unconfined space, it undergoes di↵raction, which is characterized by a change of geometry from a quasi-planar wave to a cylindrical or spherical wave [1–4]. Depending on the mixture composition, thermodynamic conditions, detonation velocity at the tube exit, the geometry of the area expansion and tube cross section, the detonation wave can be either extinguished (sub-critical regime) or re-initiated (super-critical regime) [2]. Detonation di↵raction has been extensively studied both experimentally and through numerical simulations. Since Zeldovich et al.’s pioneering work [5], numerous studies have been performed on detonation di↵raction using a variety of experimental geometry and techniques, as summarized by Schultz [4] in 2000. Since then, a number of investigations have been performed, for example by Khasainov et al. [6], Meredith et al. [7], and Nagura et al. [8]. The most advanced experimental investigation is the study of Pintgen [1, 2] who performed simultaneous imaging of the shock front and reaction zone, using schlieren and Planar Laser Induced Fluorescence (PLIF), as well as multi-exposure and stereoscopic chemiluminescence imaging. Jones et al. [10–13] and Oran et al. [14, 15] used numerical simulations to study the di↵raction of detonation and reignition process at an abrupt tube diameter increase (180 angle). Layered and uniform H2-O2(-Ar) mixtures were modeled. The chemical schemes were either 1-step or 2-step models. Arienti et al. [9] and Nagura et al. [8] performed simulations of detonation di↵raction from a tube to an unconfined volume (90 angle). Arienti used a 1-step reaction model characterized by di↵erent activation energies whereas Nagura employed a 2-step scheme to model a stoichiometric H2-air mixture. Khasainov et al. [6] studied the di↵raction of detonation from a tube to a cone. A stoichiometric C2H2-O2 mixture was considered using a 1-step pressure dependent reaction model. Deiterding [16] studied the failure and reignition of di↵racting detonations in a stoichiometric H2-O2 mixture diluted with 70% argon. High-resolution numerical simulations with mesh refinement were performed using Euler equations. The reactivity of the mixture was described using a detailed reaction model composed of 34 irreversible reactions and 9 species. From this literature review, we see that most numerical studies have been performed using an inviscid gas model (Euler equations) and global (1or 2-step) kinetic schemes. The present study aims at performing numerical simulation of detonation di↵raction using a viscous gas model (Navier-Stokes equations) and accurate chemical and thermodynamic models to allow for a more realistic comparison with the experimental data of Pintgen in terms of di↵racting detonation behaviour and structure.