DSMC and R13 modeling of the adiabatic surface

Abstract

Adiabatic wall boundary conditions for rarefied gas flows are described with the isotropic scattering model. An appropriate sampling technique for the direct simulation Monte Carlo (DSMC) method is presented, and the corresponding macroscopic boundary equations for the regularized 13-moment system (R13) are obtained. DSMC simulation of a lid driven cavity shows slip at the wall, which, as a viscous effect, creates heat that enters the gas while there is no heat flux in the wall. Analysis with the macroscopic equations and their boundary conditions reveals that this heat flux is due to viscous slip heating, and is the product of slip velocity and shear stress at the adiabatic surface. DSMC simulations of the driven cavity with adiabatic walls are compared to R13 simulations, which both show this non-linear effect in good agreement for Kn < 0.3.

DSMC Simulation of Rarefied Gas Flow over a Backward- Facing Step: Effect of Expansion ratio

“…From the pressure contours in figures (10)(11)(12) variations in pressure along the flow direction is similar for all ER larger than 2. The pressure changes in the transverse direction are more pronounced in the vicinity of the step, with regions of higher pressure near the walls and pressure reducing sharply away from the wall.…”

“…As known, the density of air decreases gradually as altitude increases, and the rarefaction effects become more evident near outer space as the flow changes from near continuum to free molecular regime. As the degree of rarefaction rises, the continuum solvers which use Navier-Stokes equations lose their credibility [12], whereas the Boltzmann equation can aptly describe the behavior of gas flow at every degree of rarefaction [12]. The Boltzmann equation is given as [13], [14]:…”

Thermoacoustic waves in pulsed Pirani sensor MEMS. Numerical Investigation

“…This gas-surface interaction model can be described by keeping the particle velocity magnitude invariant and the velocity direction of the reflected particle set based on isotropic scattering boundary conditions in the half unite sphere. We discuss this problem first in the spherical coordinate system   1 2 , , [13,18]. The three velocity components on the inner cylinder for adiabatic case reflected by: …”

Numerical simulation of the pulsed Pirani gauges

“…This gas-surface interaction model can be described by keeping the particle velocity magnitude invariant and the velocity direction of the reflected particle set based on isotropic scattering boundary conditions in the half unite sphere. We discuss this problem first in the spherical coordinate system , , c 0, , 0, 2 2 [13,18]. The three velocity components on the inner cylinder for adiabatic case reflected by:…”