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

技術的要約：流体力学的極限としてのリチャーズ方程式

問題提起
リチャーズ方程式は、1世紀近くにわたり、不飽和水流の標準的なマクロモデルとして機能してきた。しかし、これはマクロな質量保存則とダルシー・バックナムのフラックス則から伝統的に仮定されたものであり、微視的な理論から導出されたものではない。その結果、その有効領域は特定されておらず、優先流、ヒステリシス、および動的な非平衡効果といった現象を捉えることができない。これらの失敗は、内部の細孔スケールの構成が局所平衡状態にあると仮定し、含水率（$\theta$）またはマトリックポテンシャル（$\psi$）を十分な状態変数としていることに起因する。本論文は、連続体キネティック方程式（CKE）の流体力学的極限としてリチャーズ方程式を厳密に導出することにより、微視的な細孔スケールの物理学とマクロな連続体記述の間の溝を埋めるものである。

手法
導出には、ボルツマン方程式からナビエ・ストークス方程式を導出するプロセスと同様に、空間的な均質化と時間的なスケール分離を分けた、2段階の極限プロセスを用いる。

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（REV内）の再分配演算子である。

2. 時間極限（チャップマン・エンスコグ展開）: ダムケラー数 $Da = \tau_{\text{redis}}/\tau_{\text{forcing}}$（細孔スケールの再分配時間とマクロな強制時間の比）によって制御されるチャップマン・エンスコグ（CE）展開を行う。導出では $Da \ll 1$、すなわち intra-REV の再分配がマクロな強制に対して高速であることを仮定する。

   • 変数 $g$ は、$g = g^{(0)} + Da \, g^{(1)} + O(Da^2)$ として展開される。
   • 0次項 ($g^{(0)}$) は、局所平衡ステップ関数（臨界半径 $r^*$ より小さいすべての細孔は満たされ、それ以外は空である状態）に対応する。
   • 1次項 ($g^{(1)}$) は、線形化された再分配演算子 $J = C'|_{g_{\text{eq}}}$ を逆転させることで決定される。

主な貢献と結果

• リチャーズ方程式の導出: 本論文は、CE展開を CKE に適用し、$Da \ll 1$ の条件を満たすとき、リチャーズ方程式が唯一のマクロ方程式となることを証明している。1次補正の可解条件（フレドホルムの代替定理）が、マクロな質量保存則を与え、システムを閉じる。
• 輸送係数としての透水係数: 透水係数 $K(\theta)$ は、経験的な構成関係としてではなく、ガス分子運動論における粘性に構造的に類似した、1次の輸送係数として導出される。これは以下のように表される：
  $$K = \phi \langle \kappa | J^{-1} | S_0 C_w \rangle_f$$
  ここで $\kappa$ は細孔スケールのコンダクタンス、$J^{-1}$ は線形化された再分配演算子の逆演算子、$S_0 C_w$ はソース項である。この定式化は、$K$ が媒体の緩和スペクトルの特性であることを明らかにしている。
• 古典的モデルの回収:
  • 保持曲線: 平衡ステップ関数は、ヤング・ラプラスの法則を介して保持曲線 $\psi(\theta)$ を自然に回収する。ブルックス・コリーやファン・ゲンテンのような特定の形式は、基礎となる細孔径分布（PSD）の座標変換として特定される。
  • Mualemの不均一性ペナルティ: 本論文は、対数正規 PSD に対する Mualem の因子 $\exp(-4\sigma^2)$ を、経験的なフィットとしてではなく、CE 輸送係数の直列パス補正（ノイマン級数）の厳密な再総和として導出する。
• 二重および多重浸透率モデル: 緩和スペクトルが複数のギャップを持つ場合（例：バイモーダル PSD）に、この枠組みは拡張される。CKE を異なるスペクトルバンド（高速な「本体」細孔と低速な「スロート」細孔）に投影し、各バンド内で CE を適用しつつ、$O(1)$ でバンド間再分配を保持することで、結合されたリチャーズ方程式を導出する。これにより、Gerke–van Genuchten の二重浸透率モデル、Weiler の IN3M、および移動・不動モデルが第一原理から回収される。決定的なことに、バンド間の交換係数は、フィッティングではなく細孔の連結性から計算される。
• 破綻のメカニズム: 本論文は、この簡約化の限界を明示的に特徴づけている。展開は以下の場合に失敗する：
  1. $Da \sim 1$ （キネティックな強制）の場合、鋭いフロントや非平衡状態が生じる。
  2. パーコレーション閾値 ($\theta \to \theta_c$) においてスペクトルギャップが閉じる場合、緩和時間が発散し、透水係数が消失する。
  3. スケール分離が失われる場合（連結性長が REV サイズに接近する場合）。

意義
本論文は、リチャーズ方程式が単なる一つのモデルではなく、土壌水のキネティック理論の特定の流体力学的極限であり、高速な局所再分配（$Da \ll 1$）の条件下でのみ厳密に有効であることを主張している。空間極限と時間極限を分離することで、方程式が成立するために必要な構造的仮定（局所平衡、統計的均質性）を明確にしている。

本研究は、現象論的に扱われがちな現象に対し、第一原理に基づいた基礎を提供する。二重浸透率モデルが、多バンドのスペクトル構造から自然に生じること、およびヒステリシスがキネティック理論における湿潤および乾燥演算子の非可換性に由来することを実証している。さらに、透水係数を緩和スペクトルの構造的特性として特定しており、パーコレーション閾値付近での挙動や細孔幾何学への依存性を厳密に説明している。この導出は、古典的な土壌物理学の結果（保持曲線、透水係数の積分）を、ボルツマン方程式からナビエ・ストークス方程式が導出されるのと同様に、基礎となるキネティック演算子のアーキテクチャの結果として統一するものである。