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

기술 요약: 유체역학적 극한으로서의 리차드(Richards) 방정식

문제 제기
리차드 방정식은 거의 한 세기 동안 불포화 수분 흐름을 위한 표준 거시적 모델로 사용되어 왔다. 그러나 이 방정식은 전통적으로 미시적 이론으로부터 유도되기보다는 거시적 질량 보존과 다르시-벅킹엄(Darcy–Buckingham) 플럭스 법칙으로부터 가정되어 왔다. 결과적으로, 그 유효 영역은 명시되지 않은 상태로 남아 있으며, 우선적 흐름(preferential flow), 이력 현상(hysteresis), 그리고 동적 비평형 효과와 같은 현상을 포착하는 데 실패한다. 이러한 실패는 내부 기공 구조가 국부적 평형 상태에 영구적으로 있다고 가정하여 수분 함량($\theta$)이나 매트릭 포텐셜($\psi$)을 충분한 상태 변수로 간주하는 데서 기인한다. 본 논문은 연속체 운동론적 방정식(CKE)으로부터 리차드 방정식을 엄밀하게 유도함으로써, 미시적 기공 규모의 물리학과 거시적 연속체 기술 사이의 간극을 해결한다.

방법론
유도 과정은 공간적 균질화와 시간적 척도 분리를 분리하는 2단계 극한 과정을 채택하며, 이는 볼츠만 방정식으로부터 나비에-스토크스(Navier–Stokes) 방정식을 유도하는 방식과 유사하다.

1. 공간적 극한 ($\varepsilon \to 0$): 대표 요소 체적(REV)이 점으로 수축된다. 이 순수 운동학적 단계는 중간 척도의 마스터 방정식을 폐쇄된 연속체 운동론적 방정식(CKE)으로 변환한다:
   $$\partial_t g + \nabla \cdot \mathbf{F} = C[g]$$
   여기서 $g(r, x, t)$는 물로 채워진 반지름 $r$인 기공의 비율이다. 플럭스 $\mathbf{F}$는 수송 계수와 구동 구배가 아직 분리되지 않은 "얽힌(entangled)" 사전 폐쇄 플럭스이며, $C[g]$는 REV 내부의 재분배 연산자이다.

2. 시간적 극한 (채프먼-엔스코그(Chapman–Enskog) 축약): 본 논문은 기공 규모의 재분배 시간과 거시적 강제력 시간의 비율인 다름코어 수(Damköhler number), $Da = \tau_{redis}/\tau_{forcing}$에 의해 제어되는 채프먼-엔스코그(CE) 전개를 수행한다. 유도 과정은 $Da \ll 1$임을 가정하며, 이는 거시적 강제력에 비해 REV 내부의 재분산이 빠름을 의미한다.

   • 상태 변수 $g$는 $g = g^{(0)} + Da \, g^{(1)} + O(Da^2)$로 전개된다.
   • 0차 항($g^{(0)}$)은 국부 평형 단계 함수(모든 기공 중 임계 반지름 $r^*$보다 작은 기공은 가득 차 있고, 나머지는 비어 있음)에 해당한다.
   • 1차 항($g^{(1)}$)은 선형화된 재분배 연산자 $J = C'|_{g_{eq}}$를 역전(inverting)함으로써 결정된다.

주요 기여 및 결과

• 리차드 방정식의 유도: 본 논문은 $Da \ll 1$ 조건 하에서 CE 전개를 CKE에 적용했을 때 리차드 방정식이 유일한 거시적 방정식임을 증명한다. 1차 보정 항의 해 존재성 조건(Fredholm alternative)은 거시적 질량 보존 방정식을 도출하여 시스템을 폐쇄한다.
• 수송 계수로서의 수리 전도도: 수리 전도도 $K(\theta)$는 경험적인 구성 관계가 아니라, 기체 운동론의 점성과 구조적으로 유사한 1차 수송 계수로 유도된다. 이는 다음과 같이 표현된다:
  $$K = \phi \langle \kappa | J^{-1} | S_0 C_w \rangle_f$$
  여기서 $\kappa$는 기공 규모의 전도도, $J^{-1}$은 선형화된 재분배 연산자의 역행렬, $S_0 C_w$는 소스 항이다. 이 공식은 $K$가 매질의 완화 스펙트럼(relaxation spectrum)의 특성임을 밝혀낸다.
• 고전적 모델의 회복:
  • 보유 곡선(Retention Curves): 평형 단계 함수는 영-라플라스(Young–Laplace) 법칙을 통해 자연스럽게 보유 곡선 $\psi(\theta)$를 회복한다. 브룩스-코리(Brooks–Corey) 및 반 겐텐(van Genuchten)과 같은 특정 형태들은 기저의 기공 크기 분포(PSD)의 좌표 변환으로 식별된다.
  • Mualem의 불균질성 패널티: 본 논문은 로그 정규 분포를 가진 PSD에 대한 Mualem 인자 $\exp(-4\sigma^2)$를 경험적 적합이 아닌, CE 수송 계수의 직렬 경로 보정(Neumann series)의 정확한 재합(resummation)으로서 유도한다.
• 이중 및 다중 투과율 모델: 본 프레임워크는 완화 스펙트럼이 다중 간극(예: 이봉형 PSD)을 갖는 경우로 확장된다. CKE를 서로 다른 스펙트럼 밴드(빠른 "본체" 기공과 느린 "목" 기공)로 투영하고, $O(1)$ 내에서 교차 밴드 완화를 유지하면서 각 밴드 내에서 CE를 적용함으로써 결합된 리차드 방정식을 유도한다. 이를 통해 Gerke–van Genuchten 이중 투과율 모델, Weiler의 IN3M, 그리고 이동-부동(mobile–immobile) 모델을 제1원리로부터 회복한다. 결정적으로, 밴드 간 교환 계수는 피팅이 아닌 기공 연결성으로부터 계산된다.
• 붕괴 메커니즘: 본 논문은 이 축약의 한계를 명시적으로 규명한다. 전개는 다음과 같은 경우 실패한다:
  1. $Da \sim 1$ (운동론적 강제력)인 경우, 날카로운 전면(front)과 비평형 상태를 초래한다.
  2. 침투 임계값($\theta \to \theta_c$)에서 스펙트럼 간극이 닫히는 경우, 완화 시간이 발산하고 전도도가 소멸한다.
  3. 척도 분리가 상실되는 경우 (연결 길이 가 REV 크기에 근접할 때).

의의
본 논문은 리차드 방정식이 단순히 여러 모델 중 하나가 아니라, 토양 수분의 운동론적 이론으로부터 도출되는 특정한 유체역학적 극한이며, 빠른 국부 재분산($Da \ll 1$) 조건 하에서 엄격하게 유효하다고 주장한다. 공간적 극한과 시간적 극한을 분리함으로써, 이 방정식이 성립하기 위해 필요한 구조적 가정(국부 평형, 통계적 균질성)을 명확히 한다.

이 연구는 흔히 현상론적으로 취급되는 현상들에 대해 제1원리 기반의 토대를 제공한다. 이중 투과율 모델이 기공 네트워크의 다중 밴드 스펙트럼 구조로부터 자연스럽게 발생함을 보여주며, 이력 현상이 운동론에서의 습윤 및 건조 연산자의 비가환성(non-commutativity)에서 기인함을 입증한다. 또한, 수리 전도도를 완화 스펙트럼의 구조적 특성으로 식별함으로써, 침투 임계값 근처에서의 거동과 기공 기하학에 대한 의존성을 엄밀하게 설명한다. 이 유도 과정은 고전적인 토양 물리학 결과들(보유 곡선, 전도도 적분)을 볼츠만 방정식으로부터 나비에-스토크스 방정식이 도출되는 것과 유사하게, 기저의 운동론적 연산자 구조의 결과로서 통합한다.