Application of Interparticle Contact Algorithm Based on Riemann Solution in SPH

Smooth particle hydrodynamics (SPH) is a new effective numerical method for modeling problems involving large deformation. However when the SPH method is used to simulate shock waves come across contact interface, tensile instability occurs at the contact interface. In this paper we use contact algorithm which describes contact interaction between SPH particles by means of Riemann solution to simulate the propagation of wave in shock tube and projectile impact problems, and compare the numerical results with the analytical solutions. From the compared results, we can conclude that the modified SPH method do well in removing instabilities and improving accuracy of calculation. And due to the introduce of Riemann solution, there is no need to use artificial viscosity and heat flux, which simplifies the program code.