Define aquifer — define_aquifer • anem

Define an aquifer as confined or unconfined and with saturated hydraulic conductivity. Other optional parameters can be defined including boundaries and recharge.

define_aquifer(aquifer_type, Ksat, ...)

Arguments

aquifer_type	"confined" or "unconfined"
Ksat	Saturated hydraulic conductivity
...	Optional parameters including boundaries, recharge_params, and aquifer thickness

Value

This function returns an S3 "aquifer" object that behaves as a list and contains the named items:

aquifer_type: "confined" or "unconfined"
Ksat: Saturated hydraulic conductivity
h0: Hydraulic head of the undisturbed aquifer (i.e., no pumping)
z0: Thickness of the aquifer (applies only to confined aquifers)
bounds: Boundaries of the aquifer, if any (generated using define_boundaries())

Examples

(aquifer <- define_aquifer("confined",1e-4))
#> # aquifer_type: confined
#> # Ksat: 1e-04
#> # h0: --
#> # z0: --
#> # bounds: 
#>     No boundaries
#> # recharge (undisturbed water table): 
#>     No recharge zones
(aquifer <- define_aquifer("confined",1e-4,h0=100,z0=10))
#> # aquifer_type: confined
#> # Ksat: 1e-04
#> # h0: 100
#> # z0: 10
#> # bounds: 
#>     No boundaries
#> # recharge (undisturbed water table): 
#>     No recharge zones

bounds_df1 <- data.frame(bound_type=c("CH","NF","NF","NF"),
  x1=c(0,10,13,1),y1=c(0,10,9,-1),x2=c(10,13,1,0),y2=c(10,9,-1,0))
aquifer_confined <- define_aquifer("confined",1e-3,bounds=bounds_df1,h0=100,z0=10)
print(aquifer_confined)
#> # aquifer_type: confined
#> # Ksat: 0.001
#> # h0: 100
#> # z0: 10
#> # bounds: 
#> # A tibble: 4 x 8
#>     bID bound_type      m      b     x1     y1     x2     y2
#>   <int> <chr>       <dbl>  <dbl>  <dbl>  <dbl>  <dbl>  <dbl>
#> 1     1 CH          0.909  0.455 10.8   10.3   -0.201  0.271
#> 2     2 NF         -1.1   22.2   10.8   10.3   12.2    8.73 
#> 3     3 NF          0.909 -2.36  12.2    8.73   1.20  -1.27 
#> 4     4 NF         -1.1    0.05  -0.201  0.271  1.20  -1.27 
#> # recharge (undisturbed water table): 
#>     No recharge zones

bounds_df2 <- data.frame(bound_type=c("CH","CH","NF","NF"),
  x1=c(0,0,10,10),y1=c(0,10,10,0),x2=c(0,10,10,0),y2=c(10,10,0,0))
aquifer_confined <- define_aquifer("unconfined",1e-3,bounds=bounds_df2,h0=100)
aquifer_confined
#> # aquifer_type: unconfined
#> # Ksat: 0.001
#> # h0: 100
#> # bounds: 
#> # A tibble: 4 x 8
#>     bID bound_type     m     b    x1    y1    x2    y2
#>   <int> <chr>      <dbl> <dbl> <dbl> <dbl> <dbl> <dbl>
#> 1     1 CH           Inf     0     0    10     0     0
#> 2     2 CH             0    10     0    10    10    10
#> 3     3 NF           Inf    10    10    10    10     0
#> 4     4 NF             0     0     0     0    10     0
#> # recharge (undisturbed water table): 
#>     No recharge zones

recharge_params <- list(recharge_type="F",recharge_vector=c(0,0,3,3),flow=1,x0=3,y0=3)
aquifer <- define_aquifer("confined",1e-3,h0=50,z0=10,recharge=recharge_params)