anaflow.flow.ext_grf¶
- ext_grf(time, rad, S_part, K_part, R_part, dim=2, lat_ext=1.0, rate=- 0.0001, h_bound=0.0, K_well=None, struc_grid=True, lap_kwargs=None)[source]¶
The extended “General radial flow” model for transient flow.
The general radial flow (GRF) model by Barker introduces an arbitrary dimension for radial groundwater flow. We introduced the possibility to define radial dependet conductivity and storage values.
This solution is based on the grf model presented in [Barker88].
- Parameters
time (
numpy.ndarray
) – Array with all time-points where the function should be evaluatedrad (
numpy.ndarray
) – Array with all radii where the function should be evaluatedS_part (
numpy.ndarray
) – Given storativity values for each diskK_part (
numpy.ndarray
) – Given conductivity values for each diskR_part (
numpy.ndarray
) – Given radii separating the disks (including r_well and r_bound).dim (
float
, optional) – Fractional dimension of the aquifer. Default:2.0
lat_ext (
float
, optional) – Lateral extend of the aquifer. Default:1.0
rate (
float
, optional) – Pumpingrate at the well. Default: -1e-4h_bound (
float
, optional) – Reference head at the outer boundary R_part[-1]. Default:0.0
K_well (
float
, optional) – Conductivity at the well. Default:K_part[0]
struc_grid (
bool
, optional) – If this is set toFalse
, the rad and time array will be merged and interpreted as single, r-t points. In this case they need to have the same shapes. Otherwise a structured r-t grid is created. Default:True
lap_kwargs (
dict
orNone
optional) – Dictionary forget_lap_inv
containing method and method_dict. The default is equivalent tolap_kwargs = {"method": "stehfest", "method_dict": None}
. Default:None
- Returns
Array with all heads at the given radii and time-points.
- Return type
References
- Barker88
Barker, J. ‘’A generalized radial flow model for hydraulic tests in fractured rock.’’, Water Resources Research 24.10, 1796-1804, 1988