Shock simulation by the particle method SPH (Q1062556)

The particle method SPH (smoothed particle hydrodynamics), which can be derived from the exact equations of motion by multiplying each term by an appropriately chosen kernel and then integrating over the given domain, is applied to one-dimensional shock tube problems. The basic idea is to use a new kind of artificial viscosity in order to get a satisfactory shock simulation, since standard artificial viscosities like bulk viscosity or von Neumann-Richtmyer viscosity lead to under- and overshooting effects or smearing of the shock front. It is further shown that a super Gaussian interpolating kernel with fourth-order errors gives significantly better shock resolutions than the standard Gaussian kernel with second-order errors. Numerical results are given and compared with those obtained by \textit{G. A. Sod} [ibid. 27, 1-31 (1978; Zbl 0387.76063)] and \textit{B. van Leer} [ibid. 32, 101-136 (1979)].