step_slug_cbp: step_slug_cbp

In jkennel/hydrorecipes: Hydrogeology steps

step_slug_cbp

Description

Cooper, Bredehoeft and Papadopulos, 1967 Slug test solution

Usage

step_slug_cbp(
  .rec,
  times,
  radius = 1,
  radius_casing = 0.15,
  radius_well = 0.15,
  specific_storage = 1e-06,
  hydraulic_conductivity = 1e-04,
  head_0 = 1,
  thickness = 1,
  n_terms = 16,
  role = "predictor",
  ...
)

Arguments

.rec	the R6 recipe object.
radius	distance from line source or center of well
radius_casing	radius of casing
radius_well	radius of well screen
specific_storage	specific storage of aquifer (L/L)
hydraulic_conductivity	the hydraulic conductivity (L/t)
head_0	initial displacement
thickness	the aquifer thickness (L)
n_terms	number of terms for laplace solution inversion
role	character - the name of the role
...	additional arguments

Value

an updated recipe

References

Cooper, H.H., J.D. Bredehoeft and S.S. Papadopulos, 1967. Response of a finite-diameter well to an instantaneous charge of water, Water Resources Research, vol. 3, no. 1, pp. 263-269.

Examples

# check vs. CBP 1967 table 1
times    <- rep(c(1.0, 2.15, 4.64), 4) * 10^(rep(c(-3, -2, -1, 0), each = 3))
dat <- list(x = times)

formula = formula(x~.)

frec1 = recipe(formula = formula, data = dat) |>
  step_slug_cbp(
    times = x,
    radius = 1.0,
    radius_casing = 1.0,
    radius_well = 1.0,
    specific_storage = 1e-1,
    hydraulic_conductivity = 1.0,
    thickness = 1.0,
    head_0 = 1.0,
    n_terms = 12L
  ) |>
  plate("dt")

jkennel/hydrorecipes documentation built on Dec. 24, 2024, 5:38 p.m.