For gases at low pressure, the ratio of the gas mean free path, λ, to the gap size (known as the Knudsen number:
Kn = λ/
h) grows. For Knudsen numbers greater than 0.1, the gas cannot be treated using the continuum Navier-Stokes equations and the Boltzmann equation must be solved instead.
At steady state, the solutions to the linearized Boltzmann equation for isothermal flow in a narrow gap between parallel plates can be expressed as a combination of Poiseuille and Couette flows. This is analogous to the solutions of the Navier-Stokes equations in the limit of small h0/
l0. Provided that the surfaces of the wall and base are identical (which is normally the case in many practical applications), the Couette contribution to the bulk fluid velocity is unchanged (it remains the mean of the wall and base velocities for identical surfaces). The Poiseuille contribution to the flow is more complicated for a rarefied gas. A practical approach, pioneered by Fukui and Kaneko (
Ref. 3) is to solve the linearized Boltzmann BGK equation over a range of Knudsen numbers and to provide an empirical fit to the flow. This results in the following form for the average velocity of the flow:
where Q(
Kn, αw , αb) is a nondimensional function of the Knudsen number (
Kn) and the tangential momentum accommodation coefficient at the wall (
αw) and base (
αb).
Q(
Kn, αw , αb) is obtained by solving the linearized Boltzmann BGK equation for steady Poiseuille flow with a range of Knudsen numbers and slip coefficients. This approach assumes that stationary solutions of the Boltzmann equation apply inside the gap, that is, that the flow can be treated as quasi-static.
Fukui and Kaneko provided data on Q(
Kn, αw , αb) for the case where
αw=αb (
Ref. 4), which was subsequently fitted to different empirical formulas by Veijola and others (
Ref. 5). Note also that additional, more accurate, data is available in
Ref. 6. Veijola provided two empirical formulas, which apply under different circumstances with various degrees of accuracy:
where μ is the gas viscosity,
p is the gas pressure,
R is the molar gas constant,
T is the temperature, and
Mn is the molar mass of the gas.
Ref. 2 also employs this definition of the mean free path.
In many applications the forces acting on the wall and base are important. The pressure in the gas can be computed correctly by solving Equations 4-10 and
4-11. However, this approach provides only the normal component of the traction acting on the wall and base. To obtain the shear forces, the approach adopted by Torczynski and Gallis (
Ref. 7) is used. They produced an empirical expression for the shear force that has the correct behavior in the free molecular flow and continuum limits as well as in the limits for the accommodation coefficient. Torczynski and Gallis solve the problem of pure Couette flow and derive an empirical function for the slip length that predicts the correct forces for the flow in the gap in several limiting cases. Their empirical expression for the slip length is given by:
where d1≈0.15 and
d2≈0.59. In principle
d1 and
d2 are variables themselves, but they were found to be constant to within the accuracy of the DSMC experiments used to derive their values.
Equation 4-7 gives the following expressions for the shear forces on the wall and base for pure Couette flow:
here Ls is derived from
Equation 4-13 and
p is obtained by solving
Equation 4-10 with
Equation 4-11. Strictly speaking, Torczynski and Gallis’ result applies for Couette flow only and was derived for a more general variable-soft-sphere gas rather than for the linearized BGK equations, using numerical simulations. From a practical perspective, it seems likely that solutions of the linearized BGK equations would also be fitted by these expressions, and in that case it should be possible to combine the forces using superposition. In the absence of a detailed proof
Equation 4-14 is not the default option for the force model, but is available as an additional option.