Function to evaluate the Theis well function: $$ W(u) = \int_{u}^{\infty} \frac{\exp{-u}}{u} du$$

theis_well_function(td, par = NULL)

Arguments

td

Dimensionless time \(t_n\) defined as $$ t_d=\frac{4Tt}{r^{2}S} $$ where \(r\) is the distance between the pumping and observation well in m, \(S\) is the storage coefficient (dimensionless), \(T\) is the Transmisivity in m2/s and \(t\) is the time in s

par

additional parameters

Value

This function returns the value of the Theis well function

References

Theis, C. The relation between the lowering of the piezometric surface and the rate and duration of discharge of a well using groundwater storage. Transactions of the American Geophysical Union, 1935, 2, 519-524

See also

Examples

u <- seq(0.05,30,0.05) td <- 1.0/u W <- theis_well_function(td, par) plot(td, W, type = "l", log = "xy")
#> Warning: 154 y values <= 0 omitted from logarithmic plot