slug_cbp: slug_cbp

In jkennel/aquifer: Well solutions for aquifer tests

slug_cbp

Description

This function calculates H/H0 for a well of finite diameter using equation 7 from Cooper, Bredehoeft and Papdopulos, 1967 Response of a Finite-Diameter Well to an Instantaneous Charge Water.

This function uses Stehfest method for the Inverse Laplace Transform.

Usage

slug_cbp(time, radius_casing = NULL, radius_screen = NULL,
  radius = NULL, storativity = 1e-05, transmissivity = 10,
  head_0 = 1, n = 12L)

Arguments

time	elapsed time since slug introduction
radius_casing	radius of the casing where the water level is
radius_screen	radius of the screened interval
radius	radius from borehole center to calculate the response
storativity	aquifer storativity
transmissivity	aquifer transmissivity
head_0	height of introduced head change
n	coefficient for Stehfest algorithm (should be < 20)

Value

vector of H/H0

Examples

t <- 5.776 * 10^(seq(-4, 3, 0.01))
s <- slug_cbp(t, radius_screen = 7.6)
plot(s~t, log = 'x', type = 'l')

jkennel/aquifer documentation built on July 31, 2022, 11:33 p.m.