van Genuchten Model

For h < 0:

se = 1 / (1 + |alpha*h|beta)gamma

Kr = ( 1 - C * D (-gamma) )2 / D (gamma/2)

C = |alpha*h| (beta-1)

D = 1 + |alpha*h| beta

For h >= 0:

se = 1

Kr = 1

where

se = effective saturation = (VMC - RMC)/(porosity - RMC)

Kr = relative hydraulic conductivity

h = pressure head

alpha = parameter, see Note to VS2DH 1.0 users

beta = parameter

gamma = 1 - 1/beta

Return to: Contents