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

Resumen Técnico: La Ecuación de Richards como Límite Hidrodinámico

Planteamiento del Problema
La ecuación de Richards ha servido como el modelo macroscópico estándar para el flujo de agua no saturada durante casi un siglo. Sin embargo, tradicionalmente se postula a partir de la conservación de la masa macroscópica y la ley de flujo de Darcy–Buckingham, en lugar de derivarse de una teoría microscópica. En consecuencia, su dominio de validez permanece sin especificar, y falla al capturar fenómenos tales como el flujo preferencial, la histéresis y los efectos de no equilibrio dinámico. Estos fallos derivan de la suposición de que la configuración interna a escala de poro está permanentemente en equilibrio local, lo que hace que el contenido de agua ($\theta$) o el potencial mátrico ($\psi$) sean variables de estado suficientes. Este artículo aborda la brecha entre la física de la escala de poro microscópica y la descripción de continuo macroscópica mediante la derivación rigurosa de la ecuación de Richards como un límite hidrodinámico de una ecuación cinética de continuo (CKE) para el agua en el suelo.

Metodología
La derivación emplea un proceso de límite de dos pasos, separando la homogeneización espacial de la separación de escalas temporales, de forma análoga a la derivación de las ecuaciones de Navier–Stokes a partir de la ecuación de Boltzmann.

1. Límite Espacial ($\varepsilon \to 0$): El Volumen Elemental Representativo (REV) se contrae hasta un punto. Este paso puramente cinemático transforma la ecuación maestra mesoscópica en una ecuación cinética de continuo (CKE) cerrada:
   $$\partial_t g + \nabla \cdot \mathbf{F} = C[g]$$
   Aquí, $g(r, x, t)$ es la fracción de poros de radio $r$ llenos de agua. El flujo $\mathbf{F}$ es un flujo de pre-cierre "entrelazado" donde los coeficientes de transporte y los gradientes de fuerza aún no están separados, y $C[g]$ es el operador de redistribución intra-REV.

2. Límite Temporal (Reducción de Chapman–Enskog): El artículo realiza una expansión de Chapman–Enskog (CE) controlada por el número de Damköhler, $Da = \tau_{redist}/\tau_{forcing}$, la relación entre el tiempo de redistribución a escala de poro y el tiempo de forzamiento macroscópico. La derivación asume que $Da \ll 1$, lo que implica que la redistribución intra-REV es rápida en comparación con el forzamiento macroscópico.

   • El estado variable $g$ se expande como $g = g^{(0)} + Da \, g^{(1)} + O(Da^2)$.
   • El orden cero ($g^{(0)}$) corresponde a la función de paso de equilibrio local (todos los poros más pequeños que un radio crítico $r^*$ están llenos, los otros vacíos).
   • El primer orden ($g^{(1)}$) se determina invirtiendo un operador de redistribución linealizado $J = C'|_{g_{eq}}$.

Contribuciones Clave y Resultados

• Derivación de la Ecuación de Richards: El artículo demuestra que la ecuación de Richards es la única ecuación macroscópica obtenida cuando se aplica la expansión CE a la CKE bajo la condición $Da \ll 1$. La condición de solubilidad para la corrección de primer orden (el alternativo de Fredholm) produce la ecuación de conservación de masa macroscópica, cerrando el sistema.
• Conductividad Hidráulica como Coeficiente de Transporte: La conductividad hidráulica $K(\theta)$ se deriva no como una relación constitutiva empírica, sino como un coeficiente de transporte de primer orden, estructuralmente análogo a la viscosidad en la teoría cinética de gases. Se expresa como:
  $$K = \phi \langle \kappa | J^{-1} | S_0 C_w \rangle_f$$
  donde $\kappa$ es la conductancia a escala de poro, $J^{-1}$ es el inverso del operador de redistribución linealizado, y $S_0 C_w$ es el término fuente. Esta formulación revela que $K$ es una propiedad del espectro de relajación del medio.
• Recuperación de Modelos Clásicos:
  • Curvas de Retención: La función de paso de equilibrio recupera naturalmente la curva de retención $\psi(\theta)$ mediante la ley de Young–Laplace. Se identifican formas específicas como Brooks–Corey y van Genuchten como transformaciones de coordenadas de la distribución de tamaño de poros (PSD) subyacente.
  • Penalización de Heterogeneidad de Mualem: El artículo deriva el factor de Mualem $\exp(-4\sigma^2)$ para PSD log-normal como una resumación exacta de la corrección de trayectoria en serie (serie de Neumann) del coeficiente de transporte CE, en lugar de un ajuste empírico.
• Modelos de Doble y Múltiple Permeabilidad: El marco se extiende a casos donde el espectro de relajación posee múltiples brechas (por ejemplo, PSD bimodal). Al proyectar la CKE sobre bandas espectrales distintas (poros de "cuerpo" rápidos y poros de "cuello" lentos) y aplicar CE dentro de cada banda manteniendo la relajación entre bandas a $O(1)$, el artículo deriva ecuaciones de Richards acopladas. Esto recupera el modelo de doble permeabilidad de Gerke–van Genuchten, el IN3M de Weiler y los modelos de móvil–inmóvil a partir de primeros principios. Crucialmente, el coeficiente de intercambio entre bandas se calcula a partir de la conectividad de los poros en lugar de ser ajustado.
• Mecanismos de Ruptura: El artículo caracteriza explícitamente los límites de la reducción. La expansión falla cuando:
  1. $Da \sim 1$ (forzamiento cinético), lo que conduce a frentes agudos y estados de no equilibrio.
  2. La brecha espectral se cierra en el umbral de percolación ($\theta \to \theta_c$), causando que el tiempo de relajación diverja y la conductividad desaparezca.
  3. Se pierde la separación de escalas (la longitud de conectividad se aproxima al tamaño del REV).

Significancia
El artículo sostiene que la ecuación de Richards no es simplemente un modelo entre muchos, sino el límite hidrodinámico específico de la teoría cinética del agua en el suelo, válido estrictamente bajo la condición de una redistribución local rápida ($Da \ll 1$). Al separar los límites espaciales y temporales, la derivación clarifica los supuestos estructurales (equilibrio local, homogeneidad estadística) necesarios para que la ecuación sea válida.

El trabajo proporciona una base de primeros principios para fenómenos a menudo tratados de forma fenomenológica. Demuestra que los modelos de doble permeabilidad surgen naturalmente de estructuras de banda espectral múltiple y que la histéresis se origina de la no conmutatividad de los operadores de humectación y secado en la teoría cinética. Además, identifica la conductividad hidráulica como una propiedad estructural del espectro de relajación, ofreciendo una explicación rigurosa de su comportamiento cerca del umbral de percolación y su dependencia de la geometría de los poros. La derivación unifica los resultados clásicos de la física del suelo (curvas de retención, integrales de conductividad) como consecuencias de la arquitectura del operador cinético subyacente, de manera similar a cómo las ecuaciones de Navier–Stokes emergen de la ecuación de Boltzmann.