Function to calculate the drawdown using a general_radial_flow solution

general_radial_flow_solution(ptest, a, t0, n, t)

Arguments

ptest

A pumping_test object

a

Slope of the straight line fitted to the drawdown data using the Cooper-Jacob approach

t0

Intercept of the straight line fitted to the drawdown data using the Cooper-Jacob approach

n

Flow dimension

t

Numeric vector with time

Value

A numeric vector with the calculated 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

data("general_radial_flow1") ptest.grf1 <- pumping_test("Well1", Q= 0.02322, r = 26.2, t = general_radial_flow1$t, s = general_radial_flow1$s) par <- list(rw = 0.1, rc = 0.1, rd = 1) ptest.grf1$additional_parameters <- par sol0 <- general_radial_flow_solution_initial(ptest.grf1) sol <- general_radial_flow_solution(ptest.grf1, sol0$a, sol0$t0, sol0$n, ptest.grf1$t) print(sol)
#> [1] 4.196671e-05 4.528718e-05 4.170410e-03 1.910431e-02 3.715632e-02 #> [6] 6.428579e-02 9.221252e-02 1.060167e-01 1.328903e-01 1.554152e-01 #> [11] 1.799381e-01 2.031719e-01 2.304936e-01 2.535721e-01 2.800244e-01 #> [16] 3.003670e-01 3.154535e-01 3.279113e-01 3.458372e-01 3.666092e-01 #> [21] 3.845094e-01 4.015599e-01 4.210028e-01 4.364319e-01 4.512190e-01 #> [26] 4.654139e-01 4.777203e-01 4.896124e-01 4.998563e-01 5.110362e-01 #> [31] 5.242348e-01 5.358193e-01 5.481350e-01 5.589683e-01 5.715418e-01 #> [36] 5.796799e-01 5.925145e-01 6.020767e-01 6.132171e-01 6.213456e-01 #> [41] 6.318978e-01 6.429807e-01 6.512730e-01 6.577678e-01 6.641446e-01 #> [46] 6.719548e-01 6.803487e-01 6.885460e-01 6.965546e-01 7.050870e-01 #> [51] 7.120412e-01 7.181836e-01 7.228872e-01 7.288432e-01 7.353436e-01 #> [56] 7.429876e-01 7.467471e-01