Molecular Simulations of Hugoniots of detonation products mixtures at chemical equilibrium: Microscopic calculation of the Chapman-Jouguet State
Résumé
In this work, we used simultaneously the Reaction Ensemble Monte Carlo (ReMC) method and the Adaptive Erpenbeck Equation Of State (AE-EOS) method to directly calculate the thermodynamical and chemical equilibrium of mixtures of detonation products on the Hugoniot curve. The ReMC method (W. R. Smith and B. Triska, J. Chem. Phys. 100, pp 3019-3027 (1994)) allows to reach the chemical equilibrium of a reacting mixture, and the AE-EOS method (J. J. Erpenbeck, Phys. Rev. A, 46, p 6406 (1992)) constrains the system to satisfy the Hugoniot relation. Once the Hugoniot curve of the detonation products mixture is established, the CJ state of the explosive can be determined. Performing a NPT simulation at P(CJ) , T(CJ) , we then calculate the direct thermodynamic properties and the following derivative properties of the system using a fluctuation method: calorific capacities, sound velocity and Gruneisen coefficient. As the composition fluctuates, and the number of particles is not necessarily constant in this ensemble, a fluctuation formula has been developed to take into account the fluctuations of mole number and composition. This type of calculation has been applied to several usual energetic materials: nitromethane, tetranitromethane, hexanitroethane, PETN and RDX.
Origine : Fichiers produits par l'(les) auteur(s)
Loading...