Open main menu

Chapman–Enskog theory provides a framework in which equations of hydrodynamics for a gas can be derived from the Boltzmann equation. The technique justifies the otherwise phenomenological constitutive relations appearing in hydrodynamical descriptions such as the Navier–Stokes equations. In doing so, expressions for various transport coefficients such as thermal conductivity and viscosity are obtained in terms of molecular parameters. Thus, Chapman–Enskog theory constitutes an important step in the passage from a microscopic, particle-based description to a continuum hydrodynamical one.

The theory is named for Sydney Chapman and David Enskog, who introduced it independently in 1916 and 1917.[1]

DescriptionEdit

The starting point of Chapman–Enskog theory is the Boltzmann equation for the 1-particle distribution function  :

 

where   is a nonlinear integral operator which models the evolution of   under interparticle collisions. This nonlinearity makes solving the full Boltzmann equation difficult, and motivates the development of approximate techniques such as the one provided by Chapman–Enskog theory.

Given this starting point, the various assumptions underlying the Boltzmann equation carry over to Chapman–Enskog theory as well. The most basic of these requires a separation of scale between the collision duration   and the mean free time between collisions  :  . This condition ensures that collisions are well-defined events in space and time, and holds if the dimensionless parameter   is small, where   is the range of interparticle interactions and   is the number density.[2] In addition to this assumption, Chapman–Enskog theory also requires that   is much smaller than any extrinsic timescales  . These are the timescales associated with the terms on the left hand side of the Boltzmann equation, which describe variations of the gas state over macroscopic lengths. Typically, their values are determined by initial/boundary conditions and/or external fields. This separation of scales implies that the collisional term on the right hand side of the Boltzmann equation is much larger than the streaming terms on the left hand side. Thus, an approximate solution can be found from

 

It can be shown that the solution to this equation is a Gaussian:

 

where   is the molecule mass and   is the Boltzmann constant.[3] A gas is said to be in local equilibrium if it satisfies this equation.[4] The assumption of local equilibrium leads directly to the Euler equations, which describe fluids without dissipation, i.e. with thermal conductivity and viscosity equal to  . The primary goal of Chapman–Enskog theory is to systematically obtain generalizations of the Euler equations which do incorporate dissipation. This is achieved by expressing deviations from local equilibrium as a perturbative series in Knudsen number  , which is small if  . Conceptually, the resulting hydrodynamic equations describe the dynamical interplay between free streaming and interparticle collisions. The latter tend to drive the gas towards local equilibrium, while the former acts across spatial inhomogeneities to drive the gas away from local equilibrium.[5] When the Knudsen number is of the order of 1 or greater, the gas in the system being considered cannot be described as a fluid.

To first order in  , one obtains the Navier–Stokes equations. Second and third orders give rise to the Burnett equations and super-Burnett equations.

Mathematical formulationEdit

Since the Knudsen number does not appear explicitly in the Boltzmann equation, but rather implicitly in terms of the distribution function and boundary conditions, a dummy parameter   is introduced to keep track of the appropriate orders in the Chapman–Enskog expansion:

 

It can be seen that small   implies the collisional term   dominates the streaming term  , which is the same as saying the Knudsen number is small. Thus, the appropriate form for the Chapman–Enskog expansion is

 

Solutions which can be formally expanded in this way are known as normal solutions to the Boltzmann equation.[6] Clearly, this class of solutions excludes non-perturbative contributions (such as  ), which appear in boundary layers or near internal shock layers. Thus, Chapman–Enskog theory is restricted to situations in which such solutions are negligible.

Substituting this expansion and equating orders of   leads to the hierarchy

 

where   is an integral operator, linear in both its arguments, which satisfies   and  . The solution to the first equation is a Gaussian:

 

for some functions  ,  , and  . It is tempting to equate these functions with the physical hydrodynamic fields defined as moments of  :

 

From a purely mathematical point of view, however, the two sets of functions are not necessarily the same for   (for   they are equal by definition). Indeed, proceeding systematically in the hierarchy, one finds that similarly to  , each   also contains arbitrary functions of   and   whose relation to the physical hydrodynamic fields is a priori unknown. One of the key simplifying assumptions of Chapman–Enskog theory is to assume that these otherwise arbitrary functions can be written in terms of the exact hydrodynamic fields and their spatial gradients. In other words, the space and time dependence of   enters only implicitly through the hydrodynamic fields. This statement is physically plausible, since for small Knudsen numbers, one expects to enter the hydrodynamic regime in which the state of the gas is determined solely by the hydrodynamic fields. In the case of  , the functions  ,  , and   are assumed exactly equal to the physical hydrodynamic fields.

While these assumptions are physically plausible, there is the question of whether solutions which satisfy these properties actually exist. More precisely, one must show that solutions exist satisfying

 

Moreover, even if such solutions exist, there remains the additional question of whether they span the complete set of normal solutions to the Boltzmann equation, i.e. do not represent an artificial restriction of the original expansion in  . One of the key technical achievements of Chapman–Enskog theory is to answer both of these questions in the positive.[6] Thus, at least at the formal level, there is no loss of generality in the Chapman–Enskog approach.

With these formal considerations established, one can proceed to calculate  . The result is[1]

 

where   is a vector and   a tensor, each a solution of a linear inhomogeneous integral equation that can be solved explicitly by a polynomial expansion. Note that the colon denotes the double dot product,   for tensors  ,  .

PredictionsEdit

To first order in the Knudsen number, the heat flux   is found to obey Fourier's law of heat conduction,[7]

 

and the momentum-flux tensor   is that of a Newtonian fluid,[7]

 

with   the identity tensor. Here   and   are constants which we now identify with the thermal conductivity and viscosity. They can be calculated explicitly in terms of molecular parameters by solving a linear integral equation; the table below summarizes the results for a few important molecular models (  is the molecule mass and   is the Boltzmann constant).[8]

Table 1: Predicted expressions for thermal conductivity and viscosity.
Model     Notes
Rigid elastic spheres of diameter       Correct to 3 decimal places.
Molecules with repulsive force         denotes the Gamma function, and   is a numerical factor. Chapman and Cowling list several values of the latter, e.g.   and  .[9]
Lennard-Jones potential:         is a function of   which can be calculated numerically. It varies from   for   to   for  .[10]


With these results, it is straightforward to obtain the Navier–Stokes equations. Taking velocity moments of the Boltzmann equation leads to the exact balance equations for the hydrodynamic fields  ,  , and  :

 

As in the previous section the colon denotes the double dot product,  . Substituting the Chapman–Enskog expressions for   and  , one arrives at the Navier–Stokes equations.

Comparison with experimentEdit

An important prediction of Chapman–Enskog theory is that viscosity is independent of density (this can be seen for each molecular model in table 1, but is actually model-independent). This surprising result traces back to James Clerk Maxwell, who inferred it in 1860 on the basis of more elementary kinetic arguments.[11] It is well-verified experimentally for gases at ordinary densities.

Table 2: Experimentally measured values of   for the first five noble gases.[12]
Helium 2.45
Neon 2.52
Argon 2.48
Krypton 2.535
Xenon 2.58

On the other hand, the theory predicts that   does depend on temperature. For rigid elastic spheres, the predicted scaling is  , while other models typically show greater variation with temperature. For instance, for molecules repelling each other with force   the predicted scaling is  , where  . Taking  , corresponding to  , shows reasonable agreement with the experimentally observed scaling for helium. For more complex gases the agreement is not as good, most likely due to the neglect of attractive forces.[13] Indeed, the Lennard-Jones model, which does incorporate attractions, can be brought into closer agreement with experiment (albeit at the cost of a more opaque   dependence; see the Lennard-Jones entry in table 1).[14]

Chapman–Enskog theory also predicts a simple relation between   and   in the form  , where   is the specific heat at constant volume and   is a purely numerical factor. For spherically symmetric molecules, its value is predicted to be very close to   in a slightly model-dependent way. For instance, rigid elastic spheres have  , and molecules with repulsive force   have   (the latter deviation is ignored in table 1). The special case of Maxwell molecules (repulsive force  ) has   exactly.[15] Since  ,  , and   can be measured directly in experiments, an easy experimental test of Chapman–Enskog theory is to measure   for the spherically symmetric noble gases. Table 2 shows that there is reasonable agreement between theory and experiment.[12]

ExtensionsEdit

The basic principles of Chapman–Enskog theory can be extended to more diverse physical models, including gas mixtures and molecules with internal degrees of freedom. In the high-density regime, the theory can be adapted to account for collisional transport of momentum and energy, i.e. transport over a molecular diameter during a collision, rather than over a mean free path (in between collisions). Including this mechanism predicts a density dependence of the viscosity at high enough density, which is also observed experimentally.

One can also carry out the theory to higher order in the Knudsen number. In particular, the third-order contribution   has been calculated by Burnett.[16] In general circumstances, however, these high-order corrections should be approached with caution, given that the Chapman–Enskog expansion may not always converge.[17] (On the other hand, the expansion is thought to be at least asymptotic to solutions of the Boltzmann equation, in which case truncating at low order still gives accurate results.)[18] Even if the higher order corrections do afford improvement in a given system, the interpretation of the corresponding hydrodynamical equations is still debated.[19]

See alsoEdit

NotesEdit

  1. ^ a b Chapman, Sydney; Cowling, T.G. (1970), The Mathematical Theory of Non-Uniform Gases (3rd ed.), Cambridge University Press
  2. ^ Balescu, Radu (1975), Equilibrium and Nonequilibrium Statistical Mechanics, John Wiley & Sons, ISBN 978-0-471-04600-4
  3. ^ Cercignani, Carlo (1975), Theory and Application of the Boltzmann Equation, Elsevier, pp. 78–79, ISBN 978-0-444-19450-3
  4. ^ Balescu, p. 450
  5. ^ Balescu, p. 451
  6. ^ a b Grad, Harold (1958), "Principles of the Kinetic Theory of Gases", in Flügge, S. (ed.), Encyclopedia of Physics, XII, Springer-Verlag, pp. 205–294
  7. ^ a b Bird, R. Bryon; Armstrong, Robert C.; Hassager, Ole (1987), Dynamics of Polymeric Liquids, Volume 1: Fluid Mechanics (2nd ed.), John Wiley & Sons, pp. 10–11
  8. ^ Chapman & Cowling, chapter 10
  9. ^ Chapman & Cowling, p. 172
  10. ^ Chapman & Cowling, p. 185
  11. ^ Maxwell, James (1860), "V. Illustrations of the dynamical theory of gases.—Part I. On the motions and collisions of perfectly elastic spheres", Philosophical Magazine, 19 (124): 19–32, doi:10.1080/14786446008642818
  12. ^ a b Chapman & Cowling p. 249
  13. ^ Chapman & Cowling, pp. 230-232
  14. ^ Chapman & Cowling, pp. 235-237
  15. ^ Chapman & Cowling, pp. 247
  16. ^ Burnett, D. (1936), "The Distribution of Molecular Velocities and the Mean Motion in a Non-Uniform Gas", Proceedings of the London Mathematical Society, 40: 382, doi:10.1112/plms/s2-40.1.382
  17. ^ Santos, Andres; Brey, J. Javier; Dufty, James W. (1986), "Divergence of the Chapman–Enskog Expansion", Physical Review Letters, 56 (15): 1571–1574, doi:10.1103/PhysRevLett.56.1571, PMID 10032711
  18. ^ Grad, Harold (1963), "Asymptotic Theory of the Boltzmann Equation", The Physics of Fluids, 6 (2): 147, doi:10.1063/1.1706716
  19. ^ García-Cólin, L.S.; Velasco, R.M.; Uribe, F.J. (2008), "Beyond the Navier–Stokes equations: Burnett hydrodynamics", Physics Reports, 465 (4), doi:10.1016/j.physrep.2008.04.010

ReferencesEdit

The classic monograph on the topic:

  • Chapman, Sydney; Cowling, T.G. (1970), The Mathematical Theory of Non-Uniform Gases (3rd ed.), Cambridge University Press

Contains a technical introduction to normal solutions of the Boltzmann equation:

  • Grad, Harold (1958), "Principles of the Kinetic Theory of Gases", in Flügge, S. (ed.), Encyclopedia of Physics, XII, Springer-Verlag, pp. 205–294