Richards' equation as a hydrodynamic limit: Chapman--Enskog reduction of the continuum kinetic equation for unsaturated soil water
Abstract
Richards' equation for unsaturated water flow is derived from the kinetic theory of the pore-filling distribution g(r,x,t) of a companion paper. It separates two limits that macroscopic theory usually conflates. The spatial limit is purely kinematic: contracting the representative elementary volume (REV) to a point yields the closed continuum kinetic equation (CKE) d_t g + div F = C[g], with F a pore-resolved pre-closure flux and C[g] the occupancy-gated redistribution operator. The dynamics lies in the temporal limit, the subject of this paper: a Chapman-Enskog (CE) reduction of the CKE controlled by the Damkohler number Da (redistribution time over forcing time), the structural analogue of the passage from Boltzmann to Navier-Stokes. The linearized redistribution operator J is self-adjoint and negative semidefinite in the mass inner product, with a one-dimensional kernel fixing a single invariant (water) and hence one macroscopic equation; the CE hierarchy inverts the same J at every order, only the source changing. Four results follow: the equilibrium step defines the retention curve; the linearized water budget plus the inter-REV source set the first-order equation; its solvability is mass conservation, i.e. Richards' equation; and the response function gives the macroscopic flux, identifying the conductivity K as a first-order transport coefficient, the counterpart of viscosity. Mean-field reduction recovers the standard integral formula; the serial-path correction gives Mualem's heterogeneity penalty. A bimodal pore-size distribution opens a gap in the relaxation spectrum; projecting onto the bands and applying CE within each, with cross-band relaxation surviving at O(1), derives the dual- and multiple-permeability models from first principles, the exchange coefficient set by inter-band connectivity. When Da is not small the expansion breaks down and the full CKE is needed.
Related
- Cracking of unsaturated clays using a newly developed small-beam bending apparatus and analysis by digital image correlation
- Fouling, Deactivation, and Self-Healing of Coordinatively Unsaturated Metal Sites in MOF-Derived Catalysts During Real Water Treatment
- A novel soil-water characteristic curve prediction model for unsaturated soil considering pore structure characteristics
- The Statistical physics of unsaturated soil water: kinetic theory and non commutative pore water dynamics
- Evidence for temporal decoupling of olivine weathering and inorganic carbon removal in unsaturated soil columns
BibTeX
@misc{rigon2026richards,
title = {{Richards' equation as a hydrodynamic limit: Chapman--Enskog reduction of the continuum kinetic equation for unsaturated soil water}},
author = {Riccardo Rigon},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2607.17358},
eprint = {2607.17358},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2607.17358},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Richards' equation as a hydrodynamic limit: Chapman--Enskog reduction of the continuum kinetic equation for unsaturated soil water
AU - Riccardo Rigon
PY - 2026
DA - 2026/09/03
PB - arXiv
DO - 10.48550/arxiv.2607.17358
UR - https://arxiv.org/abs/2607.17358
AB - Richards' equation for unsaturated water flow is derived from the kinetic theory of the pore-filling distribution g(r,x,t) of a companion paper. It separates two limits that macroscopic theory usually conflates. The spatial limit is purely kinematic: contracting the representative elementary volume (REV) to a point yields the closed continuum kinetic equation (CKE) d_t g + div F = C[g], with F a pore-resolved pre-closure flux and C[g] the occupancy-gated redistribution operator. The dynamics lies in the temporal limit, the subject of this paper: a Chapman-Enskog (CE) reduction of the CKE controlled by the Damkohler number Da (redistribution time over forcing time), the structural analogue of the passage from Boltzmann to Navier-Stokes. The linearized redistribution operator J is self-adjoint and negative semidefinite in the mass inner product, with a one-dimensional kernel fixing a single invariant (water) and hence one macroscopic equation; the CE hierarchy inverts the same J at every order, only the source changing. Four results follow: the equilibrium step defines the retention curve; the linearized water budget plus the inter-REV source set the first-order equation; its solvability is mass conservation, i.e. Richards' equation; and the response function gives the macroscopic flux, identifying the conductivity K as a first-order transport coefficient, the counterpart of viscosity. Mean-field reduction recovers the standard integral formula; the serial-path correction gives Mualem's heterogeneity penalty. A bimodal pore-size distribution opens a gap in the relaxation spectrum; projecting onto the bands and applying CE within each, with cross-band relaxation surviving at O(1), derives the dual- and multiple-permeability models from first principles, the exchange coefficient set by inter-band connectivity. When Da is not small the expansion breaks down and the full CKE is needed.
ER -