Monte Carlo simulation of nonequilibrium shock fronts

Abstract

A direct simulation of a shock-tube experiment carried by AVCO (pressure of 1 Torr and velocity of 6.4 km/sec), repeated by Sharma and Gillespie (1990), and used by Park (1988) to develop his two-temperature model, is presented. Results show that the electronic ground state of N2 is not in rotational nor vibrational equilibrium with that of N2(+). Moreover, a two-temperature model is inadequate to describe nonequilibrium flows behind shocks. The role of impurities is examined. It is shown that the effects of a small fraction of H2O are insignificant. Good agreement with the measurements of Sharma and Gillespie is indicated.