Despite its fundamental importance for a broad range of applications, little is understood about the behaviour of metals during the initial phase of shock compression. Here, we present molecular dynamics (MD) simulations of shock-wave propagation through a metal allowing a detailed analysis of the dynamics of high strain-rate plasticity. Previous MD simulations have not seen the evolution of the strain from one- to three-dimensional compression that is observed in diffraction experiments. Our large-scale MD simulations of up to 352 million atoms resolve this important discrepancy through a detailed understanding of dislocation flow at high strain rates. The stress relaxes to an approximately hydrostatic state and the dislocation velocity drops to nearly zero. The dislocation velocity drop leads to a steady state with no further relaxation of the lattice, as revealed by simulated X-ray diffraction.

The key to perfect radiation endurance is perfect recovery. Since surfaces are perfect sinks for defects, a porous material with a high surface to volume ratio has the potential to be extremely radiation tolerant, provided it is morphologically stable in a radiation environment. Experiments and computer simulations on nanoscale gold foams reported here show the existence of a window in the parameter space where foams are radiation tolerant. We analyze these results in terms of a model for the irradiation response that quantitatively locates such window that appears to be the consequence of the combined effect of two length scales dependent on the irradiation conditions: (i) foams with ligament diameters below a minimum value display ligament melting and breaking, together with compaction increasing with dose (this value is typically ∼5 nm for primary knock on atoms (PKA) of ∼15 keV in Au), while (ii) foams with ligament diameters above a maximum value show bulk behavior, that is, damage accumulation (few hundred nanometers for the PKA's energy and dose rate used in this study). In between these dimensions, (i.e., ∼100 nm in Au), defect migration to the ligament surface happens faster than the time between cascades, ensuring radiation resistance for a given dose-rate. We conclude that foams can be tailored to become radiation tolerant.

We describe the first of two large-scale atomic simulation projects on materials failure performed on the 12-teraflop ASCI (Accelerated Strategic Computing Initiative) White computer at Lawrence Livermore National Laboratory. This is a multimillion-atom simulation study of crack propagation in rapid brittle fracture where the cracks travel faster than the speed of sound. Our finding centers on a bilayer solid that behaves under large strain like an interface crack between a soft (linear) material and a stiff (nonlinear) material. We verify that the crack behavior is dominated by the local (nonlinear) wave speeds, which can be in excess of the conventional sound speeds of a solid.

Abstract-Analysis of the results obtained from material simulations is important in the physical sciences. Our research was motivated by the need to investigate the properties of a simulated porous solid as it is hit by a projectile. This paper describes two techniques for the generation of distance fields containing a minimal number of topological features, and we use them to identify features of the material. We focus on distance fields defined on a volumetric domain considering the distance to a given surface embedded within the domain. Topological features of the field are characterized by its critical points. Our first method begins with a distance field that is computed using a standard approach, and simplifies this field using ideas from Morse theory. We present a procedure for identifying and extracting a feature set through analysis of the MS complex, and apply it to find the invariants in the clean distance field. Our second method proceeds by advancing a front, beginning at the surface, and locally controlling the creation of new critical points. We demonstrate the value of topologically clean distance fields for the analysis of filament structures in porous solids. Our methods produce a curved skeleton representation of the filaments that helps material scientists to perform a detailed qualitative and quantitative analysis of pores, and hence infer important material properties. Furthermore, we provide a set of criteria for finding the "difference" between two skeletal structures, and use this to examine how the structure of the porous solid changes over several timesteps in the simulation of the particle impact.

We describe the second of two large-scale atomic simulation projects on materials failure performed on the 12-teraflop ASCI (Accelerated Strategic Computing Initiative) White computer at the Lawrence Livermore National Laboratory. This investigation simulates ductile failure by using more than one billion atoms where the true complexity of the creation and interaction of hundreds of dislocations are revealed.