Transient motion of a confined rarefied gas due to wall heating or cooling

Abstract

The transient motion that arises in a confined rarefied gas as a container wall is rapidly heated or cooled is simulated numerically. The Knudsen number based on nominal gas density and characteristic container dimension is varied from near-continuum to highly rarefied conditions. Solutions are generated with the direct simulation Monte Carlo method. Comparisons are made with finite-difference solutions of the Navier–Stokes equations, the limiting free-molecular values, and (continuum) results based on a small perturbation analysis. The wall heating and cooling scenarios considered induce relatively large acoustic disturbances in the gas, with characteristic flow speeds on the order of 20% of the local sound speed. Steady-state conditions are reached after on the order of 5 to 10 acoustic time units, here based on the initial speed of sound in the gas and the container dimension. As rarefaction increases, the initial gas response time is decreased. For the case of a rapid increase in wall temperature, transient rarefaction effects near the wall greatly alter gas response compared to the continuum predictions, even at relatively small nominal Knudsen number. For wall cooling, the continuum solution agrees well with direct simulation at that same Knudsen number. A local Knudsen number, based on the mean free path and the scale length of the temperature gradient, is found to be a more suitable indicator of transient rarefaction effects.

Gas-flow animation by unsteady heating in a microchannel

“…To consider shorter time scales, several researchers have investigated the problem in the limit of sudden temperature variations by examining the kinetic response of a dilute gas to a step function change in wall temperatures. [8][9][10][11][12] Close agreement between analytical predictions and molecular simulations has been demonstrated in Ref. 12.…”

“…In nondimensional formulation, the heating signal is given by F͑t͒ = F j ͑t͒ =H͑t͒, where H is the Heaviside step function. The problem has been investigated in several works, including DSMC simulation studies 11 and theoretical analyses concentrating on the free-molecular limit. [8][9][10]12 The latter supplies the early time response ͑for times not exceeding the mean collision time͒ of arbitrarily large systems.…”

Gas-Surface Interaction Model Evaluation for DSMC Applications

“…The accuracy of the Nocilla form depends on the constraints used to fit the parameters. Comparisons of Navier-Stokes and DSMC predictions of the response of a similar system are <nven in [19]. °…”

“…Intermolecular collisions are not considered (i.e., Kn-> =°), so that the evolution of the system is determined solely by the GSI model. The characteristic time for transport across the system z = H/v mp , where v mp is the most probable molecular speed at T, steady state conditions are achieved on the order of I0z [19].…”

Unsteady Heat Transfer In a Gas Mixture

“…The transient behaviour in collisionless regime was studied in Ref. [2], while the oscillatory boundary conditions were stated in Ref. [3].…”