Stability analysis of interfacial Richtmyer-Meshkov flow of explosion-driven copper interface

Abstract

In this paper, a stability analysis is given to study the unstable mechanism of the Richtmyer-Meshkov flow of explosion-driven copper interface. The Richtmyer-Meshkov flow refers as an interfacial instability growth under shockwave incident loading. Numerical investigations are performed to check the applicability of the two-dimensional hydrocode, which is named AFE2D, and the physical models of detonation waves propagating in the high explosives, equations of state and the constitutive behaviors of solids in the analysis of Richtmyer-Meshkov flow problems. Here we theoretically analyze the two key issues of the unstable mechanism in Richtmyer-Meshkov flow in solids. The unstable mechanism includes temperature related melting mechanism and the plastic evolution related tensile fracture mechanism. In the analysis of the temperature related unstable mechanisms, the calculated temperature increase during the shockwave compression from the shock Hugoniot data in the shockwave physics is not enough to melt the material near the perturbed interface. On the other hand, the temperature increase from the translation of plastic work during perturbation growth which relats to the distribution of the cumulative effective plastic strain is also not enough to supply the thermal energy which is needed to melt the crystal lattice of solid, either. Therefore, the temperature related melting mechanism is not the main factor of the unstable growth of copper interface under explosion driven. In the analysis of the plastic tensile fracture related unstable mechanism, a scaling law between the maximum cumulative effective plastic strain and the scaled maximum amplitude of spikes is proposed to describe the relationship between the plastic deformation of material and the perturbation growth of interface. Combined with a critical plastic strain fracture criterion, the unstable condition of the scaled maximum amplitude of spikes is given. If the spikes grow sufficiently to meet the unstable condition, the interfacial growth will be unstable. Numerical simulations with varying initial configurations of perturbation and yield strength of materials show good agreement with the theoretical stability analysis. Finally, a criterion to judging whether the growth is stable is discussed in the form of competition between the temperature related unstable mechanism and the tensile fracture unstable mechanism.