In this study, an Eulerian, finite-volume method is developed for the numerical simulation of elastic-plastic response of compressible solid materials under impact loading. The governing equations of mass, momentum, and energy along with evolution equations for deviatoric stresses are solved in Eulerian conservation law form. A fourth-order CWENO shock-capturing method that was developed for gas dynamics has been extended to high strain rate solids problems. In this method fluxes are determined on a staggered grid at places where solution is smooth. As a result, the method does not rely on the solution of Riemann problems and enjoys the flexibility of using any type of equation of state. Boundary conditions at material interfaces are also treated by a special ghost cell approach. The position of material interfaces is advanced to the new time using a particle level set method. A fifth-order Godunov-type central scheme is used to solve the Hamilton–Jacobi (HJ) equation of level sets in two space dimensions. The capabilitie of the proposed method is evaluated at the end by comparing numerical results with the experimental results and the reported benchmark solutions for the Taylor rod impact, spherical groove jetting and void collapse problems. Key Words Taylor rod impact, particle level set, central weighted essentially non oscillatory.
