HPC simulations of shock front evolution for a study of the shock precursor decay in a submicron thick nanocrystalline aluminum

Abstract

The Hugoniot elastic limit (HEL, or the shock precursor) decay phenomenon was investigated under an uniaxial strain condition, in a plate-on-plate impact configuration, using large-scale molecular dynamics (MD) high performance computing (HPC) simulations on a multi-billion 5000 Å thick nanocrystalline aluminum (nc-Al) system with an average grain size of 1000 Å and at five impact velocities ranging from 0.7 to 1.5 km s−1. The averaged stress and strain distributions were obtained in the shock fronts’ travel direction using a material conserving atom slicing method. The loading paths in terms of the Rayleigh lines experienced by the atom system in the evolving shock fronts exhibited a strong dependency on the shock stress levels. This dependency decreased as the impact velocity increased from 0.7 to 1.5 km s−1. By combining the HELs from MD results with plate impact experimental data, the precursor decay for the nc-Al was predicted from nano-to-macro scale thickness range. The evolving shock fronts were characterized in terms of parameters such as the shock front thickness, shock rise time and strain rate. The MD results were further analyzed using a crystal analysis algorithm and a twin dislocation identification method to obtain the densities of the atomistic defects evolving behind the evolving shock fronts. High-fidelity large-scale HPC simulation results showed that certain dislocation partials strongly influenced the elastic–plastic transition response across the HELs. The twinning dislocations increased by more than a factor of 10 during the transition and remained constant under further shock compression.