High-Resolution Numerical Simulation of Detonation and Interface Tracking

The governing equations for the detonation of both gas and explosive are established. Based on the concept of Steger-Warming vector flux splitting, a splitting method is used for the case when the vector flux is not a first order homogenous function. Weighted essentially non-oscillatory (WENO) scheme with high resolution is adopted to identify the discontinuous detonation wave. A two-stage chemical reaction model and an ignition-growth model are employed to describe the chemical reaction process. Numerical simulations of the deflagration-to-detonation transition of premixed gas in the duct, the propagation of detonation waves and the propagation of cylindrical divergent detonation waves for condensed explosive are performed. Moreover, the interaction between shock wave and bubbles by the Level Set interface tracking and ghost fluid method is established, and the results Suggest that the Level Set method is of higher resolution for interface tracking and can be used to Simulate explosion and shock wave. These numerical results are in good agreement with experimental results, supporting the validity of the numerical method and interface tracking.