PIC numerical scheme for a two-dimensional Eulerian code

The present work uses the PIC method with operator splitting to hydroelastoplastic multimaterial flow problems in 2D Eulerian coordinates. The operator splitting procedures are carried out in both Lagrangian and remap phases with alternating direction for each time advancement to maintain the accuracy of one-dimensional scheme. A detail description of the numerical approximations is given with a sample calculation of a tungsten rod penetrating an aluminum target.