step_transport_ogata_banks

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")