Skip to contents

step_transport_fractures_solute

step_transport_fractures_solute.Rd

Sudicky and Frind 1982 solution. Two parallel fractures

Usage

step_transport_fractures_solute(
  .rec,
  time,
  distance_fracture,
  distance_matrix,
  concentration_influent = 1,
  time_influent = 0,
  concentration_initial = 0,
  fracture_aperture = 2e-04,
  fracture_spacing = 1,
  velocity = 0.1/86400,
  dispersivity_longitudinal = 0.1,
  diffusion = 1e-09,
  sorption_fracture = 0,
  sorption_matrix = 0,
  decay = 1e+15,
  density_bulk = 2.5,
  porosity = 0.1,
  tortuosity = 0.1,
  n_terms = 30L,
  role = "predictor",
  ...
)

Arguments

.rec

the R6 recipe object.

time

vector elapsed time (t)

distance_fracture

vector distance along fracture (z)

distance_matrix

vector distance into matrix (x)

concentration_influent

vector concentration history (c_in)

time_influent

vector concentration history (t_in)

concentration_initial

double concentration (c_0)

fracture_aperture

double fracture aperture (2b)

fracture_spacing

double fracture aperture (2B)

velocity

double water velocity in fracture (v)

dispersivity_longitudinal

double longitudinal dispersivity (α_l)

diffusion

double free-water diffusion coefficient (D*)

sorption_fracture

double fracture distribution coefficient (K_f)

sorption_matrix

double matrix distribution coefficient (K_m)

decay

double radioactive half-life for solute (λ)

density_bulk

double dry bulk density (ρ_b)

porosity

double porosity (θ)

tortuosity

double tortuosity (τ)

n_terms

integer number of terms for laplace inversion

role

character - the name of the role

...

additional arguments

Value

an updated recipe

References

Sudicky, E.A., Frind, E.O., Contaminant transport in fractured porous media: Analytical solutions for a system of parallel fractures, December 1982 https://doi.org/10.1029/WR018i006p01634

See also

Other transport: step_transport_fractures_heat(), step_transport_ogata_banks()

Examples

formula <- as.formula(~time+z+x)

dat <- as.data.frame(expand.grid(10^(3:8),
                         seq(0.0, 10, 1),
                         c(0.0)))

names(dat) <- c("time", "z", "x")

frec1 = recipe(formula = formula, data = dat) |>
  step_transport_fractures_solute(time = time,
                                  distance_fracture = z,
                                  distance_matrix = x) |>
  plate()