numMualem {spsh}

Function to Numerically Compute the Mualem Integral

Description

This function will calculate numerically Mualems Integral (Mualem 1976) and return Hydraulic Conductivity Values.

Usage

numMualem(h, scap, pcon = NA)

Arguments

h

vector of length l pressure head values.

scap

Capillary saturation [cm3 cm-3].

pcon

vector of soil hydraulic conductivity model parameters, the first argument is q and the second r.

Details

The numerical solution of Mualems integral relies on the trapezoidal rul of integration.

Value

returns a vector of length l of calulcated conductivity values at h.

Author(s)

Tobias KD Weber , [e-mail]

References

Mualem Y (1976). “New Model for predicting hydraulic conductivity of unsaturated porous media.” Water Resources Research, 12(3), 513–522. ISSN 0043-1397, http://www.scopus.com/inward/record.url?eid=2-s2.0-0016961814&partnerID=40&md5=e3773c228d43852dad22d656738389c9.

Examples

h <- 10^seq(-3, 6.8, length = 501)
p = c(.05, .5, .01, 1.8, 100, .5)
shyp.L <- shypFun.01110(p, h)

Ks <- p[5]
tau <- p[6]
Se <- shyp.L[['Se']]
Khrnum <- numMualem(h, pcon = tau, scap = Se) 

Khnum <- Ks * Se^tau * Khrnum

plot(log10(h), log10(shyp.L[['Kh']]), ylim = c(-10, 2.3), 
     xlim = c(-1,6), ylab = "log10 Kunsat [ cm/d ]", xlab = "pF [ - ]", type = "l", lwd = 8)
lines(log10(h), log10(Khnum), col = "red", lwd = 2)