High-Order Numerical Methods and Reaction Rate Law for Non-Ideal Detonation Modeling

This study demonstrates the ability of high-order, shock-capturing methods to simulate detonation wave propagation in right circular cylinders of non-ideal explosives. 2-D axisymmetric, reactive Euler equations were solved for the detonation of the highly non-ideal explosive ammonium nitrate-fuel oil (ANFO) using a simple pressure-dependent rate law and perfect gas equation of state. The fractional step method was applied to the convection and reaction operators using a 2nd-order Strang splitting. The flow solver consisted of a finite difference scheme where the flux is reconstructed by a modified fifth-order weighted essentially non-oscillatory (WENO5-Z) with local characteristic decomposition and local Lax-Friedrichs (LLF) flux splitting. Time integration was performed by a third order total variation diminishing (TVD) Runge-Kutta time discretization method, also known as the strong stability preserving (SSP) method