step_transport_ogata_banks: step_transport_ogata_banks

In jkennel/hydrorecipes: Hydrogeology steps

step_transport_ogata_banks

Description

Ogata, A., Banks, R.B., 1961. A solution of the differential equation of longitudinal dispersion in porous media. U. S. Geol. Surv. Prof. Pap. 411-A. 1-D, infinite source, uniform flow, constant parameters, decay, retardation

To have values match the excel sheet https://www.civil.uwaterloo.ca/jrcraig/pdf/OgataBanks.xlsm the decay coefficient needs to be scaled by the retardation coefficient.

Care must be taken so that input values do not lead to NaN. -Need to fix this.

1-D infinite source uniform flow constant parameters no decay no retardation

Usage

step_transport_ogata_banks(
  .rec,
  time,
  distance,
  concentration_initial = 1,
  velocity = 0.1,
  diffusion = 0.1,
  retardation = 1,
  decay = 0,
  role = "predictor",
  ...
)

Arguments

.rec	the R6 recipe object.
time	vector time
distance	vector x position
concentration_initial	double concentration
velocity	double velocity
diffusion	double diffusion coefficient
retardation	double retardation coefficient
decay	double decay coefficient
role	character - the name of the role
...	additional arguments

Value

Ogata-Banks solution for time and distance pairs

References

Ogata, A., Banks, R.B., 1961. A solution of the differential equation of longitudinal dispersion in porous media. U. S. Geol. Surv. Prof. Pap. 411-A. 1-D, infinite source, uniform flow, constant parameters, decay, retardation

See Also

Other transport: step_transport_fractures_heat(), step_transport_fractures_solute()

Examples

formula <- as.formula(y~x)
rows <- 100

dat <- data.frame(expand.grid(as.numeric(1:rows), as.numeric(1:10)))
names(dat) <- c('x', 'y')
frec1 = recipe(formula = formula, data = dat) |>
  step_transport_ogata_banks(time = x, distance = y) |>
  plate("dt")

jkennel/hydrorecipes documentation built on Dec. 24, 2024, 5:38 p.m.