A particle-level set-based sharp interface cartesian grid method for impact, penetration, and void collapse

An Eulerian, sharp interface, Cartesian grid method is developed for the numerical simulation of the response of materials to impact, shocks and detonations. The mass, momentum, and energy equations are solved along with evolution equations for deviatoric stresses and equivalent plastic strain. These equations are cast in Eulerian conservation law form. The Mie–Grüneisen equation of state is used to obtain pressure and the material is modeled as a Johnson–Cook solid. The ENO scheme is employed to capture shocks in combination with a hybrid particle level set technique to evolve sharp immersed boundaries. The numerical technique is able to handle collisions between multiple materials and can accurately compute the dynamics of the immersed boundaries. Results of calculations for axisymmetric Taylor bar impact and penetration of a Tungsten rod into steel plate show good agreement with moving finite element solutions and experimental results. Qualitative agreement with theory is shown for the void collapse phenomenon in an impacted material containing a spherical void.