A high-order simulation method for compressible multiphase flows with condensed-phase explosive detonation in underwater explosions

This study introduces an efficient method designed for the simulation of compressible multiphase flows associated with explosive detonation, primarily in the context of underwater explosions. The proposed approach integrates three core components: the compressible Euler equations, the level-set equation, and the program burn model. Spatial terms of the compressible Euler equations undergo discretization using a fifth-order accuracy weighted essentially non-oscillation reconstruction, while the third-order total variation diminishing Runge-Kutta scheme manages the temporal terms. The level-set method ensures accurate tracking of the multiphase interface. To detail the transition from solid, non-reactive explosives to gaseous detonation products in the condensed charge's detonation reaction zone, the program burn model based on Zeldovich, von Neumann, Doering theory. The efficacy and accuracy of the incorporated program burn model and the multiphase interface capture method are employed through four benchmark tests, exhibiting excellent agreement with previously published data from alternative numerical methods or commercial software. In conclusion, applying the proposed method to four distinct engineering scenarios facilitates a comprehensive understanding of the inherent dynamics associated with detonation and shock wave generation and propagation.