Chapman-Enskog Closure

Consider a neutral gas consisting of identical hard-sphere molecules of mass and diameter . Admittedly, this is not a particularly physical model of a neutral gas, but we are only considering it for illustrative purposes. The fluid equations for such a gas are similar to Equations (4.47)-(4.49):

Here, is the particle number density, the flow velocity, the scalar pressure, and the acceleration due to gravity. We have dropped the subscript because, in this case, there is only a single species. There is no collisional friction or heating in a single species system. Of course, there are no electrical or magnetic forces in a neutral gas, so we have included gravitational forces instead. The purpose of the closure scheme is to express the viscosity tensor,

The mean-free-path, , for hard-sphere molecules is given by

This formula is fairly easy to understand. The volume swept out by a given molecule in moving a mean-free-path must contain, on average, approximately one other molecule. Observe that is completely independent of the speed or mass of the molecules. The mean-free-path is assumed to be much smaller than the variation lengthscale, , of macroscopic quantities, so that

(4.59) |

In the Chapman-Enskog scheme, the distribution function is expanded, order by order, in the small parameter :

(4.60) |

Here, , , , and so on, are all assumed to be of the same order of magnitude. In fact, only the first two terms in this expansion are ever calculated. To zeroth order in , the kinetic equation requires that be a Maxwellian:

(4.61) |

Recall that . As is well known, there is zero heat flow or viscous stress associated with a Maxwellian distribution function (Reif 1965). Thus, both the heat flux density, , and the viscosity tensor,

It is possible to linearize the kinetic equation, and then rearrange it so as to obtain an integral equation for in terms of . This rearrangement crucially depends on the bilinearity of the collision operator. Incidentally, the equation is integral because the collision operator is an integral operator. The integral equation is solved by expanding in velocity space using Laguerre polynomials (sometimes called Sonine polynomials). It is possible to reduce the integral equation to an infinite set of simultaneous algebraic equations for the coefficients in this expansion. If the expansion is truncated, after terms, say, then these algebraic equations can be solved for the coefficients. It turns out that the Laguerre polynomial expansion converges very rapidly. Thus, it is conventional to keep only the first two terms in this expansion, which is usually sufficient to ensure an accuracy of about percent in the final result. Finally, the appropriate moments of are taken, so as to obtain expression for the heat flux density and the viscosity tensor. Strictly speaking, after evaluating , we should then go on to evaluate , so as to ensure that really is negligible compared to . In reality, this is never done because the mathematical difficulties involved in such a calculation are prohibitive.

The Chapman-Enskog method outlined previously can be applied to any assumed force law between molecules, provided that the force is sufficiently short-range (i.e., provided that it falls off faster with increasing separation than the Coulomb force). For all sensible force laws, the viscosity tensor is given by

whereas the heat flux density takes the form

Here, is the

(4.64) | ||

(4.65) |

where is the

Here,

is the collision frequency, and

(4.69) |

is the thermal velocity. The first two terms in the Laguerre polynomial expansion are shown explicitly (in the round brackets) in Equations (4.66) and (4.67).

Equations (4.66) and (4.67) have a simple physical interpretation. The viscous and thermal diffusivities of a neutral gas can be accounted for in terms of the random-walk diffusion of molecules with excess momentum and energy, respectively. Recall the standard result in stochastic theory that if particles jump an average distance , in a random direction, times a second, then the diffusivity associated with such motion is (Reif 1965). Chapman-Enskog theory basically allows us to calculate the numerical constants and , multiplying in the expressions for and , for a given force law between molecules. Obviously, these coefficients are different for different force laws. The expression for the mean-free-path, , is also different for different force laws.