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

Sintesi Tecnica: L'equazione di Richards come Limite Idrodinamico

Enunciato del Problema
L'equazione di Richards ha servito come modello macroscopico standard per il flusso d'acqua insatura per quasi un secolo. Tuttavia, essa viene tradizionalmente postulata dalla conservazione della massa macroscopica e dalla legge del flusso di Darcy-Buckingham piuttosto che essere derivata da una teoria microscopica. Di conseguenza, il suo dominio di validità rimane non specificato, e essa non riesce a catturare fenomeni quali il flusso preferenziale, l'isteresi e gli effetti di non-equilibrio dinamico. Questi fallimenti derivano dall'assunzione che la configurazione a scala di poro interna sia permanentemente in equilibrio locale, rendendo il contenuto d'acqua ($\theta$) o il potenziale matriciale ($\psi$) variabili di stato sufficienti. Questo articolo affronta il divario tra la fisica microscopica a scala di poro e la descrizione del continuo macroscopico, derivando rigorosamente l'equazione di Richards come limite idrodinamico di un'equazione cinetica del continuo (CKE) per l'acqua nel suolo.

Metodologia
La derivazione impiega un processo di limite a due fasi, separando l'omogeneizzazione spaziale dalla separazione delle scale temporali, analogamente alla derivazione delle equazioni di Navier-Stokes dall'equazione di Boltzmann.

1. Limite Spaziale ($\varepsilon \to 0$): Il Volume Elementare Rappresentativo (REV) viene contratto in un punto. Questo passaggio puramente cinematico trasforma l'equazione master mesoscopica in un'equazione cinetica del continuo (CKE) chiusa:
   $$\partial_t g + \nabla \cdot \mathbf{F} = C[g]$$
   Qui, $g(r, x, t)$ è la frazione di pori di raggio $r$ riempiti d'acqua. Il flusso $\mathbf{F}$ è un flusso di pre-chiusura "intrecciato" dove i coefficienti di trasporto e i gradienti di forza motrice non sono ancora separati, e $C[g]$ è l'operatore di ridistribuzione intra-REV.

2. Limite Temporale (Riduzione di Chapman-Enskog): L'articolo esegue un'espansione di Chapman-Enskog (CE) controllata dal numero di Damköhler, $Da = \tau_{redis}/\tau_{forcing}$, ovvero il rapporto tra il tempo di ridistribuzione a scala di poro e il tempo di forzante macroscopica. La derivazione assume che $Da \ll 1$, implicando che la ridistribuzione intra-REV sia veloce rispetto alla forzante macroscopica.

   • La variabile di stato $g$ è espansa come $g = g^{(0)} + Da \, g^{(1)} + O(Da^2)$.
   • L'ordine zero ($g^{(0)}$) corrisponde alla funzione a gradino di equilibrio locale (tutti i pori più piccoli di un raggio critico $r^*$ sono pieni, gli altri vuoti).
   • L'ordine uno ($g^{(1)}$) è determinato invertendo un operatore di ridistribuzione linearizzato $J = C'|_{g_{eq}}$.

Contributi Chiave e Risultati

• Derivazione dell'Equazione di Richards: L'articolo dimostra che l'equazione di Richards è l'unica equazione macroscopica ottenuta quando l'espansione CE viene applicata alla CKE sotto la condizione $Da \ll 1$. La condizione di solvibilità per la correzione del primo ordine (l'alternativa di Fredholm) produce l'equazione di conservazione della massa macroscopica, chiudendo il sistema.
• Conducibilità Idraulica come Coefficiente di Trasporto: La conducibilità idraulica $K(\theta)$ viene derivata non come una relazione costitutiva empirica, ma come un coefficiente di trasporto del primo ordine, strutturalmente analogo alla viscosità nella teoria cinetica dei gas. È espressa come:
  $$K = \phi \langle \kappa | J^{-1} | S_0 C_w \rangle_f$$
  dove $\kappa$ è la conduttanza a scala di poro, $J^{-1}$ è l'inverso dell'operatore di ridistribuzione linearizzato, e $S_0 C_w$ è il termine sorgente. Questa formulazione rivela $K$ come una proprietà dello spettro di rilassamento del mezzo.
• Recupero dei Modelli Classici:
  • Curve di Ritenzione: La funzione a gradino di equilibrio recupera naturalmente la curva di ritenzione $\psi(\theta)$ tramite la legge di Young-Laplace. Forme specifiche come Brooks-Corey e van Genuchten sono identificate come trasformazioni di coordinate della distribuzione delle dimensioni dei pori (PSD) sottostante.
  • Penalità di Eterogeneità di Mualem: L'articolo deriva il fattore di Mualem $\exp(-4\sigma^2)$ per le PSD log-normali come una risommazione esatta della correzione del percorso in serie (serie di Neumann) del coefficiente di trasporto CE, piuttosto che come un adattamento empirico.
• Modelli a Doppia e Multipla Permeabilità: Il framework si estende ai casi in cui lo spettro di rilassamento possiede molteplici gap (ad esempio, PSD bimodali). Proiettando la CKE su distinti bande spettrali (pori "corpo" veloci e pori "collo" lenti) e applicando la CE all'interno di ogni banda mantenendo il rilassamento tra le bande a $O(1)$, l'articolo deriva equazioni di Richards accoppiate. Questo recupera il modello a doppia permeabilità di Gerke-van Genuchten, l'IN3M di Weiler e i modelli mobile-immobile da primi principi. Crucialmente, il coefficiente di scambio tra le bande è calcolato dalla connettività dei pori anziché essere adattato.
• Meccanismi di Rottura: L'articolo caratterizza esplicitamente i limiti della riduzione. L'espansione fallisce quando:
  1. $Da \sim 1$ (forzante cinetica), portando a fronti netti e stati di non-equilibrio.
  2. Il gap spettrale si chiude al livello di percolazione ($\theta \to \theta_c$), causando la divergenza del tempo di rilassamento e la scomparsa della conducibilità.
  3. La separazione delle scale viene meno (la lunghezza di connettività si avvicina alla dimensione del REV).

Significatività
L'articolo sostiene che l'equazione di Richards non è semplicemente un modello tra i tanti, ma è il limite idrodinamico specifico della teoria cinetica dell'acqua nel suolo, valido strettamente sotto la condizione di rapida ridistribuzione locale ($Da \ll 1$). Separando il limite spaziale da quello temporale, la derivazione chiarisce le assunzioni strutturali (equilibrio locale, omogeneità statistica) necessarie affinché l'equazione sia valida.

Il lavoro fornisce una base dai primi principi per fenomeni spesso trattati fenomenologicamente. Dimostra che i modelli a doppia permeabilità derivano naturalmente da strutture spettrali a più bande (distribuzioni delle dimensioni dei pori) e che l'isteresi ha origine dalla non-commutatività degli operatori di bagnatura e asciugatura nella teoria cinetica. Inoltre, identifica la conducibilità idraulica come una proprietà strutturale dello spettro di rilassamento, offrendo una spiegazione rigorosa del suo comportamento vicino alla soglia di percolazione e della sua dipendenza dalla geometria dei pori. La derivazione unifica i risultati classici della fisica del suolo (curve di ritenzione, integrali di conducibilità) come conseguenze dell'architettura dell'operatore cinetico sottostante, in modo simile a come le equazioni di Navier-Stokes emergono dall'equazione di Boltzmann.