Nonequilibrium-Direction-Dependent Rotational Energy Model for Use in Continuum and Stochastic Molecular Simulation

A new rotational energy exchange model for direct simulation Monte Carlo and multi-temperature Navier–Stokes methods is presented. The direct simulation Monte Carlo model is based only on collision quantities and reduces to a rotational collision number in the continuum limit, applicable for use with the Jeans relaxation equation. The model is formulated based on recent molecular dynamics simulations of rotational relaxation in nitrogen (Valentini et al., Physics of Fluids, Vol. 24, No. 10, 2012, p. 106101) and accounts for the dependence of the relaxation rate on the direction to the equilibrium state. This enables a single parameterization of the model to accurately simulate rotational relaxation in both compressing and expanding flows, unlike the widely used Parker model. The direct simulation Monte Carlo model is simple to implement, numerically efficient, and accurately reproduces a range of pure molecular dynamics solutions, including isothermal relaxations, normal shock waves, and expansions. A gen...