Ansys Employee
1) use density-based explicit with explicit time stepping and the Courant number between 0.5 and 1.0. This will automatically adjust time step to the speed of propagating waves,
2) since all properties change abruptly across contact discontinuity, you may want to use 1st order discretization in space (but keep explicit time stepping)