Overview

In the previous chapters, the SPH method is introduced for solving CFD problems including problems with shock waves. This chapter extends and applies the SPH method to simulate the complicated explosion process of high explosives (HE). The combination of adaptive, meshfree and Lagrangian nature of the SPH method plus an explicit algorithm make the SPH method very attractive in treating highly dynamic phenomena with large deformations and large inhomogeneities that occur in the extremely transient HE explosion process.

A variable smoothing length model is employed in this chapter to ensure a minimal and sufficient number of neighboring particles to contribute to the discrete particle approximations at the current time step, so as to ensure a proper adaptability of the SPH to the drastic movement of particles.

The SPH formulation is performed based on the Euler equations, as the explosion is an extremely fast phenomenon. The JWL equation of state for high explosives is incorporated into the SPH equations.

Two numerical examples, a one-dimensional TNT slab detonation and a two-dimensional TNT explosive gas expansion, are investigated with comparisons to the results from other sources. The simulation results show that the SPH method can give a very good prediction for both the magnitude and the form of the detonation waves as well as the pressure distributions in the explosion process. The major physics of the HE explosion can be well captured in the simulation.

The SPH method is also applied to investigate the shaped charge explosion. The effects of different detonation...