Richards' equation as a hydrodynamic limit: Chapman--Enskog reduction of the continuum kinetic equation for unsaturated soil water

Technische Samenvatting: Richards-vergelijking als een Hydrodynamische Limiet

Probleemstelling
De Richards-vergelijking dient al bijna een eeuw als het standaard macroscopische model voor onverzadigde waterstroming. De vergelijking wordt echter traditioneel gepostuleerd vanuit macroscopische massaconservatie en de Darcy-Buckingham fluxwet, in plaats van afgeleid te worden uit een microscopische theorie. Hierdoor blijft het toepasbaarheidsdomein ongespecificeerd, en faalt de vergelijking bij het beschrijven van fenomenen zoals preferentiële stroming, hysterese en dynamische niet-evenwichtseffecten. Deze tekortkomingen komen voort uit de aanname dat de interne porie-schaal configuratie zich permanent in lokaal evenwicht bevindt, waardoor watergehalte ($\theta$) of matriële potentiaal ($\psi$) voldoende toestandsvariabelen zijn. Dit artikel adresseert de kloof tussen de microscopische porie-schaal fysica en de macroscopische continuümbeschrijving door de Richards-vergelijking rigoureus af te leiden als een hydrodynamische limiet van een continuüm kinetische vergelijking (CKE) voor bodemwater.

Methodologie
De afleiding maakt gebruik van een tweestaps limietproces, waarbij ruimtelijke homogenisatie wordt gescheiden van temporele schaalseparatie, analoog aan de afleiding van de Navier-Stokes-vergelijkingen uit de Boltzmann-vergelijking.

1. Ruimtelijke limiet ($\varepsilon \to 0$): Het Representatieve Elementair Volume (REV) wordt gecontracteerd tot een punt. Deze zuiver kinematische stap transformeert de meso-schaal meestervergelijking in een gesloten continuüm kinetische vergelijking (CKE):
   $$\partial_t g + \nabla \cdot \mathbf{F} = C[g]$$
   Hierbij is $g(r, x, t)$ het fractie van poriën met straal $r$ die gevuld zijn met water. De flux $\mathbf{F}$ is een "verstrengelde" pre-closure flux waarbij transportcoëfficiënten en drijvende gradiënten nog niet gescheiden zijn, en $C[g]$ is de intra-REV herverdelingsoperator.

2. Temporele limiet (Chapman-Enskog reductie): Het artikel voert een Chapman-Enskog (CE) expansie uit, gestuurd door het Damköhler-getal, $Da = \tau_{redis}/\tau_{forcing}$, de ratio tussen de herverdelingstijd op porie-schaal en de macroscopische forceringstijd. De afleiding gaat uit van $Da \ll 1$, wat impliceert dat intra-REV herverdeling snel is in vergelijking met macroscopische forcering.

   • De toestandsvariabele $g$ wordt geëxpandeerd als $g = g^{(0)} + Da \, g^{(1)} + O(Da^2)$.
   • De nulde orde ($g^{(0)}$) komt overeen met de lokale evenwichtstapfunctie (alle poriën kleiner dan een kritische straal $r^*$ zijn vol, andere leeg).
   • De eerste orde ($g^{(1)}$) wordt bepaald door het inverteren van een geliniariseerde herverdelingsoperator $J = C'|_{g_{eq}}$.

Kernbijdragen en Resultaten

• Afleiding van de Richards-vergelijking: Het artikel bewijst dat de Richards-vergelijking de unieke macroscopische vergelijking is die wordt verkregen wanneer de CE-expansie wordt toegepast op de CKE onder de conditie $Da \ll 1$. De oplosbaarheidsvoorwaarde voor de eerste-orde correctie (het Fredholm-alternatief) levert de macroscopische massaconservatievergelijking op, die het systeem sluit.
• Hydraulische conductiviteit als transportcoëfficiënt: De hydraulische conductiviteit $K(\theta)$ wordt niet afgeleid als een empirische constitutieve relatie, maar als een eerste-orde transportcoëfficiënt, structureel analoog aan viscositeit in de gas-kinetische theorie. Het wordt uitgedrukt als:
  $$K = \phi \langle \kappa | J^{-1} | S_0 C_w \rangle_f$$
  waarbij $\kappa$ de porie-schaal conductantie is, $J^{-1}$ de inverse van de geliniariseerde herverdelingsoperator, en $S_0 C_w$ de bronterm. Deze formulering onthult $K$ als een eigenschap van het relaxatiespectrum van het medium.
• Herstel van klassieke modellen:
  • Retentiecurves: De evenwichtsstapfunctie herstelt op natuurlijke wijze de retentiecurve $\psi(\theta)$ via de Young-Laplace wet. Specifieke vormen zoals Brooks-Corey en van Genuchten worden geïdentificeerd als coördinaattransformaties van de onderliggende poriegrootteverdeling (PSD).
  • Mualem's heterogeniteitspenaliteit: Het artikel leidt de Mualem-factor $\exp(-4\sigma^2)$ voor log-normale PSD's af als een exacte ressommatie van de seriele-padcorrectie (Neumann-reeks) van de CE-transportcoëfficiënt, in plaats van een empirische fit.
• Dual- en Multi-permeabiliteitsmodellen: Het kader wordt uitgebreid naar gevallen waarin het relaxatiespectrum meerdere openingen bezit (bijv. bimodale PSD's). Door de CKE te projecteren op afzonderlijke spectrale banden (snelle "body" poriën en trage "throat" poriën) en CE toe te passen binnen elke band terwijl cross-band relaxatie behouden blijft bij $O(1)$, leidt het artikel gekoppelde Richards-vergelijkingen af. Dit herstelt het Gerke-van Genuchten dual-permeability model, Weilers IN3M en mobile-immobile modellen vanuit eerste principes. Cruciaal is dat de inter-band uitwisselingscoëfficiënt wordt berekend uit de porie-connectiviteit in plaats van te worden gefit.
• Mechanismen van falen: Het artikel karakteriseert expliciet de grenzen van de reductie. De expansie faalt wanneer:
  1. $Da \sim 1$ (kinetische forcering), wat leidt tot scherpe fronten en niet-evenwechtstoestanden.
  2. De spectrale kloof sluit bij de percolatiedrempel ($\theta \to \theta_c$), waardoor de relaxatietijd divergeert en de conductiviteit verdwijnt.
  3. Schaalseparatie verloren gaat (connectiviteitslengte nadert de grootte van het REV).

Betekenis
Het artikel beweert dat de Richards-vergelijking niet slechts één model onder vele is, maar de specifie beheerste hydrodynamische limiet van de kinetische theorie van bodemwater, strikt geldig onder de voorwaarde van snelle lokale herverdeling ($Da \ll 1$). Door ruimtelijke en temporele limieten te scheiden, verheldert de afleiding de structurele aannames (lokaal evenwicht, statistische homogeniteit) die vereist zijn om de vergelijking geldig te laten zijn.

Het werk biedt een fundament vanuit de eerste principes voor fenomenen die vaak fenomenologisch worden behandeld. Het toont aan dat dual-permeability modellen natuurlijk voortkomen uit multi-band spectrale structuren in het porienetwerk en dat hysterese voortkomt uit de niet-commutativiteit van bevochtigings- en droogoperatoren in de kinetische theorie. Bovendien identificeert het de hydraulische conductiviteit als een structurele eigenschap van het relaxatiespectrum, wat een rigoureuze verklaring biedt voor het gedrag ervan nabij de percolatiedrempel en de afhankelijkheid van de poriegeometrie. De afleiding verenigt klassieke bodemfysische resultaten (retentiecurves, conductiviteitsintegralen) als consequenties van de onderliggende kinetische operatorarchitectuur, vergelijkbaar met hoe de Navier-Stokes-vergelijkingen voortkomen uit de Boltzmann-vergelijking.