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 VS2DT 2.5 users
beta = parameter
gamma = 1 - 1/beta
Return to: Full migration procedure for data prepared for VS2DT 2.5