Calculates the well function of the general radial flow model

general_radial_flow_well_function(td, par)

Arguments

td

Numeric vector with the dimensionless time

par

A list with the parameters flow dimension(n), dimensionless radius (rd), and the vector of stehfest coefficients (coeffs)

Value

A numeric vector with the dimensionless drawdown

References

Barker, J. A. A generalized radial flow model for hydraulic tests in fractured rock Water Resources Research, 1988, 24, 1796-1804.

See also

Examples

td <- logseq(-1, 4, 50) par <- list(n = 2, rd = 10, coeffs = stehfest_coeffs(8)) W <- general_radial_flow_well_function(td, par) plot(td, W, type = "l", log = "xy", main = "GRAL RADIAL FLOW: Well function")
#> Warning: 11 y values <= 0 omitted from logarithmic plot