sncFun {spsh}

Non-capillary Saturation Function to Extend Other Functions

Description

The general purpose method to calculate numerically the effective non-capillary saturation is directly obtained from any arbritrary expression for the rescaled capillary saturation function as described by (Weber et al. 2019). Examples of capillary saturation functions are the well known (van Genuchten 1980), (Kosugi 1996), (Fredlund and Xing 1994) functions.

Usage

sncFun(h, scap)

Arguments

h

A vector of n pressure head values for which scap was calculated

scap

vector of n monotonically decreasing capillary saturation function values calculated by shypFun, rescaled between 0 and 1.

Value

A vector of n elements with calculated saturation content of the non-capillary part.

Author(s)

Tobias KD Weber , [e-mail]

References

van Genuchten MT (1980). “Closed-form equation for predicting the hydraulic conductivity of unsaturated soils.” Soil Science Society of America Journal, 44(5), 892–898. ISSN 03615995 (ISSN).
Kosugi K (1996). “Lognormal distribution model for unsaturated soil hydraulic properties.” Water Resources Research, 32(9), 2697–2703. ISSN 0043-1397.
Fredlund DG, Xing A (1994). “Equations for the soil-water characteristic curve.” Canadian Geotechnical Journal, 31(4), 521–532. ISSN 00083674 (ISSN), doi: 10.1139/t94-061.

Examples

# set variables
p <- c(0.1, 0.4, 0.01, 2, 100, .5)
h <- 10^seq(-2, 6.8, length = 197)

# Calculate the capillary and non-capillary saturation function.
Se <- shypFun(p, h, shpmodel = "01110")$Se     
Snc <- sncFun(Se)