get_recharge_undisturbed_potential: Get undisturbed hydraulic potential due to recharge

In gopalpenny/anem: Analytical element models for simple aquifers

Description

Get undisturbed hydraulic potential due to recharge. This is the water table potential (head or discharge potential) without any well pumping.

Usage

get_recharge_undisturbed_potential(loc, aquifer)

Arguments

loc	Location for evaluating Vector location as c(x,y) or data.frame containing x, y columns
aquifer	Aquifer containing containing h0 and recharge$params (see Details)

Details

Note that this function ignores any boundaries, meaning the user must ensure that the recharge does not conflict with boundaries. For example, if uniform flow is specified, it will pass through any boundary, including no-flow boundaries.

The aquifer and aquifer$recharge must be correctly specified using define_aquifer and/or define_recharge.

Value

Returns hydraulic potential at the location(s) using helper functions for each recharge_type / aquifer_type. If aquifer_type=="confined", hydraulic head is returned. If it is "unconfined", discharge potential is returned. Note that these functions are not exported to the package namespace:

• 1. Constant flow, "F": get_recharge_flow_potential()

• 2. Head boundaries, "H": not yet implemented

• 3. Recharge divide, "D": get_recharge_divide_potential()

Examples

## Flow - confined aquifer
recharge_params <- list(recharge_type="F",recharge_vector=c(0,0,1,1),flow=1,x0=0,y0=0)
aquifer <- define_aquifer("confined",1,h0=0,z0=1,recharge=recharge_params)
get_recharge_undisturbed_potential(c(1/sqrt(2),1/sqrt(2)), aquifer)

loc <- expand.grid(x=-2:2,y=-2:2)
loc$h <- get_recharge_undisturbed_potential(loc, aquifer)

library(ggplot2)
ggplot(loc) + geom_raster(aes(x,y,fill=h)) + scale_fill_gradient2()

recharge_params <- list(recharge_type="F",
  recharge_vector=c(-1,-5,0,0),flow=1,x0=0,y0=0)
aquifer <- define_aquifer("confined",1,h0=10,z0=1,recharge=recharge_params)
loc <- expand.grid(x=-100:100,y=-100:100)
loc$h <- get_recharge_undisturbed_potential(loc, aquifer)
ggplot(loc) + geom_raster(aes(x,y,fill=h)) + scale_fill_gradient2()

gopalpenny/anem documentation built on Dec. 20, 2020, 5:27 a.m.