step_transport_fractures_heat: step_transport_fractures_heat

In jkennel/hydrorecipes: Hydrogeology steps

step_transport_fractures_heat

Description

Sudicky and Frind 1982 solution adapted for heat. Two parallel fractures.

Usage

step_transport_fractures_heat(
  .rec,
  time,
  distance_fracture,
  distance_matrix,
  temperature_influent = 15,
  time_influent = 0,
  temperature_initial = 10,
  fracture_aperture = 2e-04,
  fracture_spacing = 1,
  velocity = 0.1/86400,
  thermal_conductivity_water = 0.615,
  thermal_conductivity_solids = 3.4,
  specific_heat_water = 4192,
  specific_heat_solids = 908,
  density_water = 1,
  density_solids = 2.5,
  porosity = 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)
temperature_influent	vector temperature history (t_in)
time_influent	vector time of influent values (t_in)
temperature_initial	double temperature (t_0)
fracture_aperture	double fracture aperture (2b)
fracture_spacing	double fracture aperture (2B)
velocity	double water velocity in fracture (v)
thermal_conductivity_water	double water thermal conductivity (λ_f)
thermal_conductivity_solids	double solids thermal conductivity (λ_s)
specific_heat_water	double specific heat of water
specific_heat_solids	double specific heat of solid particles
density_water	double density of the water (ρ_w)
density_solids	double density of the solid particles (ρ_s)
porosity	double matrix porosity (θ)
n_terms	integer the number of laplace terms
role	character - the name of the role
...	additional arguments

Value

an updated recipe

See Also

Other transport: step_transport_fractures_solute(), step_transport_ogata_banks()

Examples

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

dat <- setDT(expand.grid(10^(3:8),
                         seq(0.0, 100, 1),
                         c(0.0, 0.05)))

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


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

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