Influence of interatomic bonding potentials on detonation properties.

The dependences of the macroscopic detonation properties of a two-dimensional (2D) diatomic (AB) molecular system on the fundamental molecular properties were investigated. This includes examining the detonation velocity, reaction zone thickness, and critical width as functions of the exothermicity (Q) of the gas-phase reaction [AB --> (1/2)(A(2) + B(2))] and the gas-phase dissociation energy (D(e)(AB)) for AB --> A + B . Following previous work, molecular dynamics (MD) simulations with a reactive empirical bond-order potential were used to characterize the shock-induced response of a diatomic AB molecular solid, which exothermically reacts to produce A2 and B2 gaseous products.

Shock waves in polycrystalline iron.

The propagation of shock waves through polycrystalline iron is explored by large-scale atomistic simulations. For large enough shock strengths the passage of the wave causes the body-centered-cubic phase to transform into a close-packed phase with most structure being isotropic hexagonal-close-packed (hcp) and, depending on shock strength and grain orientation, some fraction of face-centered-cubic (fcc) structure. The simulated shock Hugoniot is compared to experiments.

Direct observation of the alpha-epsilon transition in shock-compressed iron via nanosecond x-ray diffraction.

In situ x-ray diffraction studies of iron under shock conditions confirm unambiguously a phase change from the bcc (alpha) to hcp (epsilon) structure. Previous identification of this transition in shock-loaded iron has been inferred from the correlation between shock-wave-profile analyses and static high-pressure x-ray measurements. This correlation is intrinsically limited because dynamic loading can markedly affect the structural modifications of solids.

Nanohydrodynamics simulations: an atomistic view of the Rayleigh-Taylor instability.

Nanohydrodynamics simulations, hydrodynamics on the nanometer and nanosecond scale by molecular dynamics simulations for up to 100 million particles, are performed on the latest generation of supercomputers. Such simulations exhibit Rayleigh-Taylor instability, the mixing of a heavy fluid on top of a light in the presence of a gravitational field, initiated by thermal fluctuations at the interface, leading to the chaotic regime in the long-time evolution of the mixing process. The early-time behavior is in general agreement with linear analysis of continuum theory (Navier-Stokes), and the late-time behavior agrees quantitatively with experimental observations.

Microscopic view of structural phase transitions induced by shock waves.

Multimillion-atom molecular-dynamics simulations are used to investigate the shock-induced phase transformation of solid iron. Above a critical shock strength, many small close-packed grains nucleate in the shock-compressed body-centered cubic crystal growing on a picosecond time scale to form larger, energetically favored grains. A split two-wave shock structure is observed immediately above this threshold, with an elastic precursor ahead of the lagging transformation wave.