step_transport_fractures_solute: step_transport_fractures_solute

In jkennel/hydrorecipes: Hydrogeology steps

step_transport_fractures_solute

Description

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 <- setDT(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("dt")

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