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

技术摘要：作为流体动力学极限的 Richards 方程

问题陈述
近一个世纪以来，Richards 方程一直是描述非饱和水流的标准宏观模型。然而，它传统上是基于宏观质量守恒和 Darcy–Buckingham 通量定律进行假设的，而非从微观理论中推导出来的。因此，其适用范围尚未明确，且无法捕捉诸如优先流（preferential flow）、滞后现象（hysteresis）以及动态非平衡效应等现象。这些失效源于一个假设，即内部孔隙尺度配置始终处于局部平衡状态，从而认为含水率 ($\theta$) 或基质势 ($\psi$) 足以作为状态变量。本文通过将 Richards 方程严格推导为土壤水连续动力学方程（CKE）的流体动力学极限，旨在弥合微观孔隙尺度物理学与宏观连续介质描述之间的鸿沟。

方法论
该推导采用了两步极限过程，将空间同质化与时间尺度分离分开，类似于从玻尔兹曼方程推导纳维-斯托克斯（Navier–Stokes）方程的过程。

1. 空间极限 ($\varepsilon \to 0$)： 代表性单元体积（REV）被收缩为一个点。这一纯运动学步骤将中尺度主方程转化为封闭的连续动力学方程（CKE）：
   $$\partial_t g + \nabla \cdot \mathbf{F} = C[g]$$
   此处，$g(r, x, t)$ 是半径为 $r$ 的孔隙被水填充的分数。通量 $\mathbf{F}$ 是一个“纠缠”的预封闭通量，其中传输系数与驱动梯度尚未分离，$C[g]$ 是单元内（intra-REV）再分配算子。

2. 时间极限（Chapman–Enskog 归约）： 本文通过达姆科勒数（Damköhler number）$Da = \tau_{redis}/\tau_{forcing}$（孔隙尺度再分配时间与宏观强迫时间之比）来控制 Chapman–Enskog (CE) 展开。推导假设 $Da \ll 1$，这意味着相对于宏观强迫过程，单元内的再分配过程是快速的。

   • 变量 $g$ 展开为 $g = g^{(0)} + Da \, g^{(1)} + O(Da^2)$。
   • 零阶项 ($g^{(0)}$) 对应于局部平衡阶跃函数（所有小于临界半径 $r^*$ 的孔隙均为满水状态，其余为空）。
   • 一阶项 ($g^{(1)}$) 通过对线性化再分配算子 $J = C'|_{g_{eq}}$ 求逆来确定。

核心贡献与结果

• Richards 方程的推导： 本文证明了当在 $Da \ll 1$ 的条件下对 CKE 应用 CE 展开时，Richards 方程是唯一的宏观方程。一阶修正项的可解性条件（Fredholm 择优原理）导出了宏观质量守恒方程，从而封闭了系统。
• 水力传导度作为传输系数： 水力传导度 $K(\theta)$ 被推导为一种一阶传输系数，而非经验性的本构关系，其结构类似于气体动力学中的粘度。它表示为：
  $$K = \phi \langle \kappa | J^{-1} | S_0 C_w \rangle_f$$
  其中 $\kappa$ 是孔隙尺度电导率，$J^{-1}$ 是线性化再分配算子的逆，$S_0 C_w$ 是源项。该公式揭示了 $K$ 是介质弛豫谱的一个属性。
• 经典模型的恢复：
  • 持水曲线： 平衡阶跃函数通过 Young–Laplace 定律自然地恢复了持水曲线 $\psi(\theta)$。Brooks–Corey 和 van Genuchten 等特定形式被识别为底层孔径分布（PSD）的坐标变换。
  • Mualem 异质性惩罚： 本文将针对对数正态分布的 Mualem 因子 $\exp(-4\sigma^2)$ 视为 CE 传输系数的串联路径修正（Neumann 级数）的精确重求和，而非经验拟合。
• 双重及多重渗透率模型： 该框架扩展到了弛豫谱具有多个间隙（例如双峰 PSD）的情况。通过将 CKE 投影到不同的谱带（快速“体”孔隙和慢速“喉部”孔隙）上，并在每个谱带内应用 CE 展开，同时保留 $O(1)$ 阶的跨谱带弛豫，本文推导出了耦合的 Richards 方程。这从第一性原理出发恢复了 Gerke–van Genuchten 双重渗透率模型、Weiler 的 IN3M 以及移动-不动（mobile–immobile）模型。至关重要的是，带间交换系数是通过孔隙连通性计算得出的，而非通过拟合获得。
• 失效机制： 本文明确表征了归约过程的极限。当出现以下情况时，展开会失效：
  1. $Da \sim 1$（动力学强迫），导致剧烈的锋面和非平衡态。
  2. 谱间隙在渗透阈值（$\theta \to \theta_c$）处关闭，导致弛豫时间发散且传导度消失。
  3. 尺度分离失效（连通长度接近 REV 大小）。

意义
本文声称，Richards 方程不仅仅是众多模型中的一个，它是土壤水动力学理论在快速局部再分配（$Da \ll 1$）这一严格条件下特有的流体动力学极限。通过分离空间与时间极限，推导过程阐明了使该方程成立所需的结构性假设（局部平衡、统计均匀性）。

这项工作为通常由现象学处理的现象提供了第一性原理的基础。它证明了双重渗透率模型自然产生于孔隙网络中的多谱带结构，并且滞后现象起源于动力学理论中润湿与干燥算子的非交换性。此外，它将水力传导度识别为弛豫谱的结构属性，从而为传导度在渗透阈值附近的性质及其对孔隙几何结构的依赖性提供了严密的解释。该推导将经典的土壤物理学结果（持水曲线、传导度积分）统一为底层动力学算子架构的必然结果，类似于纳维-斯托克斯方程从玻尔兹曼方程中涌现的过程。