anaflow.flow.thiem

thiem(rad, r_ref, transmissivity, rate=-0.0001, h_ref=0.0)

The Thiem solution.

The Thiem solution for steady-state flow under a pumping condition in a confined and homogeneous aquifer. This solution was presented in [Thiem1906].

Parameters:

• rad (numpy.ndarray) – Array with all radii where the function should be evaluated.

• r_ref (float) – Reference radius with known head (see h_ref).

• transmissivity (float) – Transmissivity of the aquifer.

• rate (float, optional) – Pumpingrate at the well. Default: -1e-4

• h_ref (float, optional) – Reference head at the reference-radius r_ref. Default: 0.0

Returns:

head – Array with all heads at the given radii.

Return type:

numpy.ndarray

References

[Thiem1906]

Thiem, G., ‘’Hydrologische Methoden, J.M. Gebhardt’’, Leipzig, 1906.

Notes

The parameters rad, r_ref and transmissivity will be checked for positivity. If you want to use cartesian coordiantes, just use the formula r = sqrt(x**2 + y**2)

Examples

>>> thiem([1,2,3], 10, 0.001, -0.001)
array([-0.3664678 , -0.25615   , -0.19161822])