Implementation of a dislocation-density based single-crystal model into a continuum shock hydrodynamics code

The dynamic thermomechanical responses of polycrystalline materials under shock loading are often dominated by the interaction of defects and interfaces. For example, within metals, a prescribed deformation associated with a shock wave may be accommodated by crystallographic slip, provided a sufficient population of mobile dislocations is available. However, if the deformation rate is large enough, there may be an insufficient number of freely mobile dislocations. In this case, additional dislocations may be nucleated, or alternate mechanisms (e.g. twinning, damage) activated in order to accommodate the deformation. Direct numerical simulation at the mesoscale offers insight into these physical processes that can be invaluable to the development of macroscale constitutive theories, if the mesoscale models adequately represent the anisotropic nonlinear thermomechanical response of individual crystals and their interfaces. The paper briefly outlines a continuum mesoscale modeling framework founded upon a nonlocal dislocation-density based crystal plasticity theory. The nonlocal theory couples continuum dislocation transport with an otherwise local single-crystal model employing nonlinear thermoelasticity and crystallographic plasticity. Dislocation transport is modeled by enforcing dislocation conservation at a slip-system level through the solution of advection-diffusion equations. The configuration of geometrically necessary dislocation density gives rise to a back-stress that inhibits or accentuates the flow of dislocations. In particular, this paper emphasizes recent implementation of the coupled nonlocal model into a 3D shock hydrocode and simulation results for the dynamic response of polycrystalline copper in two and three dimensions.