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

Resumo Técnico: A Equação de Richards como um Limite Hidrodinâmico

Definição do Problema
A equação de Richards tem servido como o modelo macroscópico padrão para o fluxo de água não saturada há quase um século. No entanto, ela é tradicionalmente postulada a partir da conservação de massa macroscópica e da lei de fluxo de Darcy–Buckingham, em vez de ser derivada de uma teoria microscópica. Consequentemente, seu domínio de validade permanece não especificado, e ela falha em capturar fenômenos como o fluxo preferencial, a histerese e efeitos de não-equilíbrio dinâmico. Essas falhas decorrem da suposição de que a configuração interna em escala de poro está permanentemente em equilíbrio local, tornando o conteúdo de água ($\theta$) ou o potencial matricial ($\psi$) variáveis de estado suficientes. Este artigo aborda a lacuna entre a física de escala de poro microscópica e a descrição contínua macroscópica, derivando rigorosamente a equação de Richards como um limite hidrodinâmico de uma equação cinética contínua (CKE) para a água no solo.

Metodologia
A derivação emprega um processo de limite em duas etapas, separando a homogeneização espacial da separação de escala temporal, analogamente à derivação das equações de Navier–Stokes a partir da equação de Boltzmann.

1. Limite Espacial ($\varepsilon \to 0$): O Volume Elementar Representativo (REV) é contraído a um ponto. Este passo puramente cinemático transforma a equação mestra de mesoescala em uma equação cinética contínua (CKE) fechada:
   $$\partial_t g + \nabla \cdot \mathbf{F} = C[g]$$
   Aqui, $g(r, x, t)$ é a fração de poros de raio $r$ preenchidos com água. O fluxo $\mathbf{F}$ é um fluxo de pré-fechamento "emaranhado", onde os coeficientes de transporte e os gradientes de força ainda não foram separados, e $C[g]$ é o operador de redistribuição intra-REV.

2. Limite Temporal (Redução de Chapman–Enskog): O artigo realiza uma expansão de Chapman–Enskog (CE) controlada pelo número de Damköhler, $Da = \tau_{redis}/\tau_{forcing}$, a razão entre o tempo de redistribuição na escala de poro e o tempo de força macroscópica. A derivação assume que $Da \ll 1$, implicando que a redistribuição intra-REV é rápida em comparação com a força macroscópica.

   • A variável de estado $g$ é expandida como $g = g^{(0)} + Da \, g^{(1)} + O(Da^2)$.
   • A ordem zero ($g^{(0)}$) corresponde à função degrau de equilíbrio local (todos os poros menores que um raio crítico $r^*$ estão cheios, outros vazios).
   • A primeira ordem ($g^{(1)}$) é determinada pela inversão de um operador de redistribuição linearizado $J = C'|_{g_{eq}}$.

Principais Contribuições e Resultados

• Derivação da Equação de Richards: O artigo prova que a equação de Richards é a única equação macroscópica obtida quando a expansão CE é aplicada à CKE sob a condição $Da \ll 1$. A condição de solvabilidade para a correção de primeira ordem (alternativa de Fredholm) resulta na equação de conservação de massa macroscópica, fechando o sistema.
• Condutividade Hidráulica como um Coeficiente de Transporte: A condutividade hidráulica $K(\theta)$ é derivada não como uma relação constitutiva empírica, mas como um coeficiente de transporte de primeira ordem, estruturalmente análogo à viscosidade na teoria cinética dos gases. Ela é expressa como:
  $$K = \phi \langle \kappa | J^{-1} | S_0 C_w \rangle_f$$
  onde $\kappa$ é a condutância na escala de poro, $J^{-1}$ é o inverso do operador de redistribuição linearizado, e $S_0 C_w$ é o termo de fonte. Esta formulação revela $K$ como uma propriedade do espectro de relaxação do meio.
• Recuperação de Modelos Clássicos:
  • Curvas de Retenção: A função degrau de equilíbrio recupera naturalmente a curva de retenção $\psi(\theta)$ via lei de Young–Laplace. Formas específicas como Brooks–Corey e van Genuchten são identificadas como transformações de coordenadas da distribuição de tamanho de poros (PSD) subjacente.
  • Penalidade de Heterogeneidade de Mualem: O artigo deriva o fator de Mualem $\exp(-4\sigma^2)$ para PSDs log-normais como uma ressumação exata da correção de caminho serial (série de Neumann) do coeficiente de transporte CE, em vez de um ajuste empírico.
• Modelos de Permeabilidade Dupla e Múltipla: O arcabouço se estende a casos onde o espectro de relaxação possui múltiplos hiatos (ex: PSDs bimodais). Ao projetar a CKE em bandas espectrais distintas (poros de "corpo" rápidos e poros de "gargalo" lentos) e aplicar CE dentro de cada banda enquanto mantém a relaxação entre bandas em $O(1)$, o artigo deriva equações de Richards acopladas. Isso recupera o modelo de permeabilidade dupla de Gerke–van Genuchten, o IN3M de Weiler e modelos de móvel–imóvel a partir de primeiros princípios. Crucialmente, o coeficiente de troca entre bandas é computado a partir da conectividade dos poros, em vez de ser ajustado.
• Mecanismos de Ruptura: O artigo caracteriza explicitamente os limites da redução. A expansão falha quando:
  1. $Da \sim 1$ (força cinética), levando a frentes nítidas e estados de não-equilíbrio.
  2. O hiato espectral se fecha no limiar de percolação ($\theta \to \theta_c$), fazendo com que o tempo de relaxação divirja e a condutividade desapareça.
  3. A separação de escala é perdida (comprimento de conectividade aproxima-se do tamanho do REV).

Significância
O artigo afirma que a equação de Richards não é meramente um modelo entre muitos, mas o limite hidrodinâmico específico da teoria cinética da água no solo, válido estritamente sob a condição de redistribuição local rápida ($Da \ll 1$). Ao separar os limites espacial e temporal, a derivação esclarece as suposições estruturais (equilíbrio local, homogeneidade estatística) necessárias para que a equação seja válida.

O trabalho fornece uma base de primeiros princípios para fenômenos frequentemente tratados de forma fenomenológica. Demonstra que modelos de permeabilidade dupla surgem naturalmente de estruturas espectrais de múltiplas bandas na rede de poros e que a histerese se origina da não-comutatividade dos operadores de molhagem e secagem na teoria cinética. Além disso, identifica a condutividade hidráulica como uma propriedade estrutural do espectro de relaxação, oferecendo uma explicação rigorosa para seu comportamento próximo ao limiar de percolação e sua dependência da geometria dos poros. A derivação unifica resultados clássicos da física do solo (curvas de retenção, integrais de condutividade) como consequências da arquitetura do operador cinético subjacente, de forma semelhante a como as equações de Navier–Stokes emergem da equação de Boltzmann.