The direct simulation Monte Carlo method has evolved over 50 years into a powerful numerical technique for the computation of complex, nonequilibrium gas flows. In this context, "nonequilibrium" means that the velocity distribution function is not in an equilibrium form due to a low number of intermolecular collisions within a fluid element. In hypersonic flow, nonequilibrium conditions occur at high altitude and in regions of flowfields with small length scales. In this paper, the theoretical basis of the direct simulation Monte Carlo technique is discussed. In addition, the methods used in direct simulation Monte Carlo are described for simulation of high-temperature, real gas effects and gas-surface interactions. Several examples of the application of direct simulation Monte Carlo to flows around blunt hypersonic vehicles are presented to illustrate current capabilities. Areas are highlighted where further research on the direct simulation Monte Carlo technique is required.
Nomenclature
A= surface element of area, m 2 C = particle velocity vector, m∕s e r = specific rotational energy, J∕kg f = probability density function g = relative velocity, m∕s Kn = Knudsen number k = Boltzmann's constant, J∕kg L = characteristic length scale, m m = mass of a particle, kg N = average number of particles in a cell N p = number of particles in a cell N t = number of iterations n = number density, m −3 P c = collision probability r = particle position vector, m T = temperature, K t = time, s V = cell volume, m 3 Δf = rate of change due to collision processes ε act = activation energy, J ε r = rotational energy, J ε tot = total collision energy, J λ = mean free path, m ζ = number of rotational degrees of freedom σ = collision cross section, m 2 τ r = rotational relaxation time, s ω = viscosity temperature exponent