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