| A hybrid weighted essentially non-oscillatory (WENO)/centered-difference (CD) numerical method, with low numerical dissipation, high-order shock-capturing, and structured adaptive mesh refinement (SAMR), has been developed for the direct numerical simulation (DNS) of the multicomponent, compressive, reactive Navier-Stokes equations. The method enables accurate resolution of diffusive processes within reaction zones. This numerical method is verified with a series of one- and two-dimensional test problems, including a convergence test of a two-dimensional unsteady reactive double Mach reflection problem. Validation of the method is conducted with experimental comparisons of three applications all of which model multi-dimensional, unsteady reactive flow: an irregular propane detonation, shock and detonation bifurcations, and spark ignition deflagrations.;The numerical approach combines time-split reactive source terms with a high-order, shock-capturing scheme specifically designed for diffusive flows. A description of the order-optimized, symmetric, finite difference, flux-based, hybrid WENO / centered-difference scheme is given, along with its implementation in a high-order SAMR framework. The implementation of new techniques for discontinuity flagging, scheme-switching, and high-order prolongation and restriction is described. In particular, the refined methodology does not require upwinded WENO at grid refinement interfaces for stability, allowing high-order prolongation and thereby eliminating a significant source of numerical diffusion within the overall code performance.;A minimally reduced irregular detonation mixture mechanism (22 species and 53 reversible reactions) is developed and combined with the WENO-CD numerical method to accurately model two-dimensional hydrocarbon (propane) detonations with detailed chemistry and transport. First of its kind, resolved double Mach reflection (DMR) detonation simulations with a large hyrdocarbon mixture are presented. Detailed discussions and comparisons of the influence of grid resolution, lower-order numerical methods, and inviscid approximations are made in addition to the detailed presentation of fluid dynamics found in an unsteady, highly unstable, reactive DMR simulation. Also conducted are direct experimental comparisons to soot foils and schlieren images with an unresolved large-scale propane detonation channel simulation.;The numerical method is also applied to the DNS of two other problems, detonation/shock bifurcations and spark ignited deflagrations. Through the resolution of viscous/diffusive scales, new insights into how a bifurcated foot develops after a detonation end wall reflection, and how geometry can influence the development of a flame kernel after spark ignition are found. |
