Narrow Results

High-resolution extensions to six Riemann solvers and three flux vector splitting schemes are developed within the framework of a reconstruction-evolution approach. Third-order spatial accuracy is achieved using two different piecewise parabolic reconstructions and a weighted essentially nonoscillatory scheme. A three-stage TVD Runge–Kutta time stepping is employed for temporal integration. The modular development of solvers provides an ease in selecting a reconstruction scheme and/or a Riemann solver/flux vector splitting scheme. The performances of these high-resolution solvers are compared for several one- and two-dimensional test cases. Based on a comprehensive assessment of the solutions obtained with all solvers, it is found that the use of the weighted essentially nonoscillatory reconstruction with the van Leer flux vector splitting scheme provides solutions for a variety of problems with acceptable accuracy.

In this paper, the specific expression for pressure and sound speed in chemical reaction zone of condensed explosives are theoretically deduced, and a new method for deriving the partial derivative of pressure in respect of every conserved quantity is proposed. Combined with the third-order TVD Runge–Kutta method, we develop a parallel solver using the fifth-order high-resolution weighted essentially non-oscillatory (WENO) finite difference scheme to simulate detonation diffraction for two-dimensional condensed explosives. The numerical simulation results revealed the forming reasons of the low-pressure region, the low-density region, the "vortex" region and the "dead zone" in the vicinity of the corner. Furthermore, it demonstrated that the retonation will generate along the inner wall, and it plays an important role in the process of detonation diffraction.