Molecular Simulation of Shocked Materials Using Reaction Ensemble Monte Carlo: Part I. Application to Nitrogen Dissociation

Report No. ARL-TR-3983
Authors: John K. Brennan; Betsy M. Rice
Date/Pages: November 2006; 32 pages
Abstract: We demonstrate the applicability of the Reaction Ensemble Monte Carlo (RxMC) simulation method for calculating the shock Hugoniot properties of a material. The method does not require interaction potentials that simulate bond breaking or bond forming; it requires only the intermolecular potentials and the ideal-gas partition functions for the reactive species that are present. By performing Monte Carlo sampling of forward and reverse reaction steps, RxMC provides information on the chemical equilibria states of the shocked material including the density of the reactive mixture and the mole fractions of the reactive species. We illustrate the methodology for shocked liquid N2, where we find excellent agreement with experimental measurements. This is the first of two reports describing the applicability of the RxMC method to simulating the shock Hugoniot properties of materials.
Distribution: Approved for public release
  Download Report ( 0.234 MBytes )
If you are visually impaired or need a physical copy of this report, please visit and contact DTIC.

Last Update / Reviewed: November 1, 2006