Three-dimensional Elastic-plastic Hydrodynamics Adaptive Euler Method

The "MEPH3D" is one of those three dimensional multi-component elastic and plastic hydrodynamics programs. It can be used for computing high-speed impact and explode problems, and this program has already been applied to many practical problems. Based on three dimensional elastic-plastic program MEPH3D, we do some modifications to make it much stronger and more all-purpose, so that we can combine this program with the parallel adaptive mesh refinement package "PARAMESH" which is is a package of Fortran 90 subroutines designed to provide an application developer with an easy route to extend an existing serial code which uses a logically cartesian structured mesh into a parallel code with adaptive mesh refinement(AMR). The routine to advance the solution is done in three orientations by turns. In each orientation we use Lagrangian-Remapping method to compute: in Lagrangian routine multi-component fluid in cell method (MFIC) is applied; in transport routine we apply high-resolution method "MUSCL" to single-component grid cells, and at multi-component grid cells "Youngs" interface method is used. We define many rules to test refinement or derefinement such as the interface of the matters, the surface of detonation wave, the gradient of the flow field, and so on. To keep conservation laws a new interpolation function is introduced. And we choose to use no-permanent guard cell method, so as to save the memory. The "morton curve" method manages the distribution of grid blocks among the processors in an effort to achieve load balance, and by improving data locality to lower communication costs. And then we give out the results of some high-speed impact and explode problems to validate the effectiveness.