The Statistical physics of unsaturated soil water: kinetic theory and non commutative pore water dynamics

Resumen Técnico: Física Estadística del Agua en Suelos No Saturados

Planteamiento del Problema
La dinámica de los medios porosos no saturados es inherentemente irreversible, caracterizada por configuraciones fuera del equilibrio que dependen de la historia del forzamiento (humedecimiento frente a secado). Los modelos clásicos, principalmente la ecuación de Richards, asumen un equilibrio termodinámico local y una única variable de estado escalar (contenido de agua $\theta$ o potencial mátrico $\psi$). Este marco falla al explicar:

1. Histéresis: Estados macroscópicos idénticos que corresponden a configuraciones internas diferentes dependiendo de la trayectoria de humedecimiento/secado.
2. Flujo Preferencial: La coexistencia de la redistribución lenta, controlada por capilaridad, y el flujo rápido de derivación (bypass) a través de poros grandes durante forzamientos intensos, lo cual viola los supuestos de dominio único de la ecuación de Richards.
3. Dependencia de la Tasa: La incapacidad de los modelos independientes de la tasa para capturar la dependencia de la histéresis y los regímenes de flujo respecto a la velocidad del forzamiento.

Los modelos existentes de dominio dual o permeabilidad dual abordan estos problemas de manera fenomenológica, introduciendo variables de estado adicionales y términos de intercambio sin derivarlos de una descripción física común.

Metodología
El artículo desarrolla una teoría de mecánica estadística pasando a través de tres escalas, tratando la fracción de ocupación de poros $g(r, x, t)$ (la fracción de poros de radio $r$ llenos con agua en la posición $x$ y tiempo $t$) como la variable de estado fundamental.

1. Microescala (Tasas Cinéticas):

   • Las transferencias de agua entre poros están gobernadas por tasas de Hagen–Poiseuille $\kappa(r)$ y un potencial de conducción $\Phi(r, r')$, definido como la diferencia en los potenciales químicos entre clases de poros (incorporando términos capilares y gravitatorios).
   • Se deriva una constante de tasa dependiente del estado $\kappa(r, g)$, que transiciona de la dinámica de Lucas–Washburn (invasión inicial rápida) al flujo de Hagen–Poiseuille (estado estacionario) a medida que aumenta la fracción de llenado $g$.

2. Mesoescala (Ecuación Maestra):

   • El promedio sobre un Volumen Elemental Representativo (VER) produce una ecuación maestra para $g$.
   • La ecuación incluye términos de ganancia ($G$) y pérdida ($L$) que representan el intercambio de agua entre clases de poros. Estos términos utilizan un "Stosszahlansatz" (hipótesis de caos molecular), asumiendo la independencia estadística de las ocupaciones del donante y el receptor, lo que conduce a una forma bilineal $g(r')[1-g(r)]$.
   • Se demuestra que el sistema es un flujo de gradiente de Onsager, que desciende un funcional de energía libre de Gibbs $F[g]$ con una movilidad dependiente del estado. Esto asegura el teorema H ($dF/dt \leq 0$) para sistemas isotérmicos y sin forzamiento.

3. Escala del Continuo (Ecuación de Campo):

   • La contracción del VER a un punto promueve a $g$ a un campo continuo.
   • La ecuación fundamental resultante es una ley de balance:
     $$\partial_t g + \nabla \cdot \mathbf{F} = C[g] - E - T$$
     donde $C[g]$ es el operador de redistribución intra-VER, $\nabla \cdot \mathbf{F}$ es la divergencia de transporte (flujo inter-VER), y $E, T$ son sumideros (evaporación, absorción radicular).
   • El flujo $\mathbf{F}$ es no local en el radio del poro, acoplando diferentes clases de poros mediante la matriz de conectividad $C(r, r')$.

Contribuciones Clave y Resultados

• Unificación de Regímenes de Flujo: La teoría introduce un único parámetro de control adimensional, el número de Damköhler resuelto por poro $Da(r, x) = \tau_{redis}(r) / \tau_{forcing}(x)$, que compara el tiempo de relajación interno con la escala de tiempo del forzamiento.

  • $Da \ll 1$ (Cuasi-estático): El sistema sigue el equilibrio local. La reducción de Chapman–Enskog de la ecuación cinética recupera la ecuación de Richards, donde el potencial mátrico $\psi$ y la conductividad hidráulica $K(\theta)$ emergen como límites.
  • $Da \gtrsim 1$ (Cinético): La redistribución no puede seguir el ritmo del forzamiento. El agua ocupa vías geométricamente accesibles (poros grandes) en lugar de aquellas termodinámicamente óptimas, reproduciendo naturalmente el flujo preferencial sin supuestos de dominio dual ad hoc.

• Origen Geométrico de la Histéresis: La histéresis no se atribuye a la biestabilidad por poro (ej. efectos de botella de tinta), sino a la holonomía de un "haz de forzamiento" (forcing bundle).

  • Los operadores de humedecimiento ($L_W$) y secado ($L_D$) no conmutan ($[L_W, L_D] \neq 0$) debido a la asimetría entre el llenado no selectivo (precipitación) y el vaciado selectivo (drenaje/evaporación restringido por la topología).
  • Esta no conmutatividad crea una fase geométrica (estado residual) tras un ciclo de forzamiento cerrado.
  • Predicción: El área del lazo de histéresis $H$ escala cuadráticamente con la intensidad del forzamiento ($H \propto I^2$) para ciclos lentos, una firma ausente en los modelos independientes de la tasa.

• Leyes Constitutivas desde la Microestructura:

  • La conductividad hidráulica $K(\theta)$ se deriva como un coeficiente de transporte (límite espectral de la densidad de la movilidad) en lugar de un parámetro ajustado.
  • La teoría recupera los límites clásicos: modelos de haces capilares (límite diagonal), modelos de camino crítico (límite espectral cerca de la percolación) y modelos de presión capilar dinámica (aproximación de un momento).
  • La capacidad de campo se identifica como el umbral de percolación donde la red de agua se fragmenta ($\theta_c$), causando que la conductividad macroscópica desaparezca.

• Variable de Estado Irreducible: Fuera del límite cuasi-estático, el contenido de agua escalar $\theta$ es insuficiente para caracterizar el sistema. La distribución completa $g(r, x, t)$ es la variable de estado irreducible.

Significancia y Reivindicaciones
El artículo afirma proporcionar un marco constructivo y fundamentado termodinámicamente que unifica fenomenologías previamente distintas (flujo de Richards, flujo preferencial, histéresis) en una única ecuación cinética.

• Rigor Teórico: La teoría se presenta como un análogo de los medios porosos de la ecuación de Boltzmann, donde $g(r)$ corresponde a la distribución de velocidades $f(v)$, y la reducción de Chapman–Enskog refleja la derivación de las ecuaciones de Navier–Stokes.
• Poder Predictivo: Todos los insumos (distribución de tamaño de poros $f(r)$, conectividad $C(r, r')$, tortuosidad $\bar{\tau}$) son propiedades geométricas medibles mediante imágenes de micro-CT, lo que no requiere calibración contra datos macroscópicos de retención o conductividad.
• Falsabilidad: La teoría ofrece predicciones específicas y contrastables, notablemente el escalamiento $H \propto I^2$ del área del lazo de histéresis con la intensidad del forzamiento y la existencia de un umbral de intensidad derivado para el inicio del flujo preferencial basado en la estructura del suelo y las propiedades del fluido.
• Alcance: El marco es generalizable a otros sistemas con transiciones restringidas y dependientes de la historia (ej. deformación plástica, histéresis magnética) y puede extenderse para incluir términos de adsorción, osmóticos o térmicos mediante el refinamiento del potencial de conducción $\Phi$ sin alterar el álgebra de operadores central.

El artículo concluye que la dicotomía entre el flujo de equilibrio y el de no equilibrio no es una cuestión de elegir diferentes modelos, sino de reconocer el cruce continuo controlado por el número de Damköhler dentro de una única teoría cinética unificada.