Organize the hydrodynamic limit by a small parameter

The Knudsen number compares a molecular mean free path with a macroscopic variation length. The relevant length may be a channel width, a temperature-gradient length, or a wavelength; choosing the largest apparatus dimension can conceal a much steeper local variation. A temporal condition also compares the collision time with the field-evolution time. Small spatial Knudsen number alone does not control a rapidly driven high-frequency disturbance or an initial layer that still contains strong kinetic distortions.

After nondimensionalization, the streaming terms are of order one while collisions carry an inverse small parameter. Write the distribution as a leading term plus successive corrections. The leading collision equation sets the zeroth-order shape to a local Maxwellian, assuming ordinary elastic collision connectivity. Small parameter refers to a scale hierarchy, not automatically to low flow velocity: Mach number and Knudsen number are separate quantities, though some asymptotic regimes couple them.

Linearize the collision operator about the local Maxwellian. At first order, streaming of that Maxwellian drives a correction through the linearized operator. Its invariant moments must vanish for a solution to exist, because collisions cannot create the protected moments. These solvability conditions are the leading conservation equations. In Chapman–Enskog theory, time derivatives are expanded consistently and lower-order balance laws replace them when forming the next correction; one does not invert the collision operator blindly on its null space.

Choose each correction to have zero density, momentum, and energy moments. This matching convention ensures the hydrodynamic fields remain defined by the full distribution rather than being counted again in the correction. The nonzero higher moments of the first correction then supply viscous stress and heat flux. A finite moment truncation and a Chapman–Enskog gradient expansion are related closure tools, but they are not interchangeable derivations and need not have the same stability properties.

If gradients and nonhydrodynamic corrections remain small, inserting the first correction into moment fluxes gives the Navier–Stokes–Fourier constitutive laws. This is a formal asymptotic construction, not a universal error theorem for every initial or boundary-value problem. Near a wall, particles arriving from the boundary have an angular distribution set by accommodation and reflection. A kinetic boundary layer of order a mean free path can coexist with a well-controlled bulk continuum limit.

The heat lab tests one consequence of a diffusive closure. Reset and set mode to 1 and diffusivity to 0.5, then Step once; predict decay exponent 0.25. Reset and set mode to 4 and diffusivity to 0.5, then Step once; predict exponent 4. Faster spatial variation damps faster within that equation. The exact single-mode display contains no mean free path or collision time and cannot establish that the high-mode continuum approximation remains valid.