Analysis of high explosive burn methodologies in cylinder test simulations

High explosive (HE) detonation is typically modeled at the continuum scale with either a programmed burn or reactive burn method. Reactive burn methods use a shock induced reaction rate which controls the transition from explosive reactants to products. This requires a fine computational mesh which ideally resolves the reaction zone. In programmed burn methods, the HE lighting time and detonation velocity are determined from an evolution equation which should then be consistent with energy release from the equation of state (EOS). These methods can utilize mesh resolutions much larger than the reaction zone, but at the cost of decoupling the shock front from the reaction. This work investigates simulations of the cylinder test, which measures the expansion of an explosive-filled copper tube, using different programmed and reactive burn models. These include programmed burn with a Huygen's construction wavefront, Detonation Shock Dynamics (DSD) with a velocity-adjusted EOS or the dynamic small resolved heat release model (DASHER) model, and the Arrhenius Wescott-Stewart-Davis (AWSD) reactive burn model. DASHER is a relatively new programmed burn method which assumes that most of the chemical energy is instantaneously released at the shock and the remainder is released via the AWSD burn model. Key advantages and disadvantages of each model are discussed, with particular attention paid to the DASHER model.