Hantush-Jacob leaky aquifer

Hantush-Jacob leaky aquifer (Hantush-Jacob, 1955)

time <- 1:2000
flow_rate <- c(runif(1000, 0, 0.003),
               rep(0.0, 1000))

# high
dat <- data.frame(time = 1:2000, flow_rate = flow_rate)
hj_100 <- recipe(flow_rate~time, dat) |>
  step_aquifer_leaky(time, 
                     flow_rate, 
                     leakage = 100, 
                     radius = 100, 
                     storativity = 1e-6,
                     transmissivity = 1e-4) |> 
  plate()

# medium
hj_200 <- recipe(flow_rate~time, dat) |>
  step_aquifer_leaky(time, 
                     flow_rate, 
                     leakage = 200, 
                     radius = 100, 
                     storativity = 1e-6,
                     transmissivity = 1e-4) |> 
  plate()

# low
hj_1000 <- recipe(flow_rate~time, dat) |>
  step_aquifer_leaky(time, 
                     flow_rate, 
                     leakage = 1000, 
                     radius = 100, 
                     storativity = 1e-6,
                     transmissivity = 1e-4) |> 
  plate()

References

Hantush, M.S. and C.E. Jacob, 1955. Non-steady radial flow in an infinite leaky aquifer, Am. Geophys. Union Trans., vol. 36, no. 1, pp. 95-100.

Prodanoff, J.H.A., Mansur, W.J. and Mascarenhas, F.C.B., 2006. Numerical evaluation of Theis and Hantush-Jacob well functions. Journal of hydrology, 318(1-4), pp.173-183.