Thiem (1906) radial flow
- class thiem_1906.ThiemWell
Bases:
objectThiem (1906) radial flow solution.
The default ThiemWell object adopts the Aq2dConf aquifer class and SteadyWell class. A single method (.hd) is defined representing the head difference between two radii under a pumping stress; for example, representing the head difference between two observation boress at different radial distances from the pumping bore, or representing the head difference between the pumping bore at its well radius and an observation bore. Default radii are r1=0.1 and r2=1.
- Attributes:
aq (obj) : Aq2dConf aquifer object. well (obj) : SteadyWell object.
- class Aq2dConf(K=1, Ss=0.0001, B=10, bot=0, name='Aq2dConf class')
Bases:
Aquifer2D confined aquifer class.
Subclass of the Aquifer class defining a horizontal, 2D confined aquifer with infinite lateral extent and confined storage.
The default Aq2dConf object has hydraulic conductivity K=1, specific storativity Ss=0.0001, aquifer saturated thickness B=10 and aquifer bottom (datum) elevation bot=0. Exceptions occur if invalid values are provided for K, Ss or B.
The .info and .draw methods display the aquifer information and diagram.
- Attributes:
Ss (float) : Aquifer specific storativity (units 1/L, default 1.0e-4).
- property B
float : Aquifer thickness (units L, default 10.0).
Setter method checks for valid values and triggers an exception if invalid values are specified.
- property D
float : Aquifer diffusivity (units L2/T, default 1.0e+4).
- property K
float : Aquifer hydraulic conductivity (units L/T, default 1.0).
Setter method checks for valid values and triggers an exception if invalid values are specified.
- property S
float : Aquifer storage coefficient (units 1, default 1.0e-3).
- property Ss
float : aquifer specific storativity (units 1/L, default 1.0e-4).
Setter method checks for valid values and triggers an exception if invalid values are specified.
- draw(dw=6)
Display a drawing of the aquifer.
- Args:
dw (float) : Width of figure (default 6.0).
- info()
Print the aquifer information.
- is_1d = False
- is_2d = True
- is_confined = True
- is_finite = False
- is_heterogeneous = False
- is_homogeneous = True
- is_infinite = True
- is_leaky = False
- is_semifinite = False
- is_unconfined = False
- class SteadyWell(x=0.0, y=0.0, r=0.05, q=0.0, pf=1, name='Steady state flowing well')
Bases:
objectSteady state flowing well class.
The default SteadyWell object has coordinates x=0.0 and y=0.0, well radius r=0.05, penetration fraction pf=1.0 and well rate q=0.0. Exceptions occur if invalid values are provided for r or pf. The well rate can be negative, positive or zero and is used to set the .state attribute as extract, inject or off.
- Attributes:
x (float) : Well x coordinate (units L, default 0.0). y (float) : Well y coordinate (units L, default 0.0). r (float) : Well radius (units L, default 0.05). q (float) : Well rate (units L3/T, default 0.0). pf (float) : Well penetration depth (fraction of aquifer depth,
default 1.0).
name (str) : Well name (default ‘’).
- info()
Print the well information.
- is_steady = True
- is_transient = False
- property pf
float : Well penetration depth.
Setter method checks for valid values and triggers an exception if invalid values are specified.
- property r
float : Well radius.
Setter method checks for valid values and triggers an exception if invalid values are specified.
- property state
str : Well state.
- draw(dw=8)
Display the definition diagram.
- Args:
dw (float) : Width of figure (default 8.0).
- hd(r1=0.1, r2=1)
Head difference between two radii.
- info()
Print the solution information.