Simulation of stationary shock waves in porous copper with smoothed particle hydrodynamics

The parameters of the shock adiabat of a porous material are calculated in the mesoscopic setting with the developed method of a moving window. The porous material is considered a framework of solid material that fills the space between pores, the mechanical properties and shock adiabat of which are often well known. The essence of the method is as follows: uncompressed material flows into the computational box with a constant velocity, while the outflow velocity from the box is chosen based on iterations so as to make the wave front immobile relative to the window, because the stationary mode of shock-wave propagation must be achieved to calculate the shock adiabats. The simulation of shock waves is performed both in the standard setting with an immobile piston (inverse-motion method) and in the system of a moving window. It is demonstrated that the wave profiles obtained with both methods are identical after the stationary mode is achieved. As an example, the shock adiabats of porous copper, which adequately reproduce experimental data for different porosities, are calculated. The proposed mesoscopic method of calculating the response of porous materials to shock compression in a moving window enables direct calculation of desired shock adiabats for porous materials that have not been studied experimentally.