Document Type


Publication Date


Digital Object Identifier (DOI)


Large-scale molecular dynamics (MD) simulations are performed to investigate shock propagation in single crystal copper. It is shown that the P-V plastic Hugoniot is unique regardless of the sample's orientation, its microstructure, or its length. However, the P-V pathway to the final state is not, and depends on many factors. Specifically, it is shown that the pressure in the elastic precursor (the Hugoniot elastic limit (HEL)) decreases as the shock wave propagates in a micron-sized sample. The attenuation of the HEL in sufficiently-long samples is the main source of disagreement between previous MD simulations and experiment: while single crystal experiments showed that the plastic shock speed is orientation-independent, the simulated plastic shock speed was observed to be orientation-dependent in relatively short single-crystal samples. Such orientation dependence gradually disappears for relatively long, micrometer-sized, samples for all three low-index crystallographic directions 100, 110, and 111, and the plastic shock velocities for all three directions approach the one measured in experiment. The MD simulations also demonstrate the existence of subsonic plastic shock waves generated by relatively weak supporting pressures.

Rights Information

Creative Commons License
This work is licensed under a Creative Commons Attribution 3.0 License.

Was this content written or created while at USF?


Citation / Publisher Attribution

Journal of Physics: Conference Series, v. 500, art. 172008