A three-dimensional sharp interface Cartesian grid method for solving high speed multi-material impact, penetration and fragmentation problems

This work presents a three-dimensional, Eulerian, sharp interface, Cartesian grid technique for simulating the response of elasto-plastic solid materials to hypervelocity impact, shocks and detonations. The mass, momentum and energy equations are solved along with evolution equations for deviatoric stress and plastic strain using a third-order finite difference scheme. Material deformation occurs with accompanying nonlinear stress wave propagation; in the Eulerian framework the boundaries of the deforming material are tracked in a sharp fashion using level-sets and the conditions on the immersed boundaries are applied by suitable modifications of a ghost fluid approach. The dilatational response of the material is modeled using the Mie-Gruneisen equation of state and the Johnson–Cook model is employed to characterize the material response due to rate-dependent plastic deformation. Details are provided on the treatment of the deviatoric stress ghost state so that physically correct boundary conditions can be applied at the material interfaces. An efficient parallel algorithm is used to handle computationally intensive three-dimensional problems. The results demonstrate the ability of the method to simulate high-speed impact, penetration and fragmentation phenomena in three dimensions.