Slip Flow Boundary Conditions
Maxwell was the first person to consider the boundary conditions required for the solution of the Boltzmann equation in the case of a slightly rarefied gas (Ref. 1). The boundary conditions he derived can be used as boundary conditions for the Navier Stokes equations in the slip flow regime, although Maxwell considered the particular case of a monatomic hard-sphere gas, so some generalization is required.
Maxwell’s boundary condition takes the form:
(7-1)
where uslip is the slip velocity, n is the boundary normal, τ is the viscous stress tensor, μ is the viscosity of the gas, ρ is its density, and Tg is its temperature. The factor G (which has dimensions of length) is given by:
where λ is the mean free path and av is the tangential momentum accommodation coefficient (for a model in which the surface reflects some molecules diffusely and some specularly, this is equivalent to the fraction of molecules which are reflected diffusely).
A slightly different definition of the mean free path is used here: in terms of Maxwell’s original parameter, l, the mean free path is λ = 2l/3. In Maxwell’s original paper, he uses a coordinate based notation, the vector notation is available in Ref. 2.
In Equation 7-1 the left-hand term represents the phenomena of viscous slip, whilst the right-hand term produces thermal slip or transpiration. Maxwell was aware that the temperature of the gas was not necessarily equal to that of the wall but formulated the boundary condition using only the gas temperatures. When the wall temperature is included, the following equations are obtained (Ref. 2 and Ref. 3):
where Tw is the wall temperature, σs is the viscous slip coefficient, σT is the thermal slip coefficient, and ζT is the temperature jump coefficient. Within a generalized Maxwell’s model the three coefficients σs, σT, and ζT are given by:
(7-2)
where κ is the thermal conductivity of the gas. The mean free path can be computed from the gas properties using the following equation (Ref. 3):
(7-3)
consequently:
In some cases an alternative definition of the mean free path, λ′, is used (as, for example, in Ref. 5):
It is possible to use this definition of the mean free path when entering user defined values for the slip coefficients.
Values of the parameters given in Equation 7-2, are available in Ref. 5, although care should be taken with the slip and temperature jump coefficients to adjust for differences between the author’s definition of the mean free path and that used in COMSOL Multiphysics.