step_transport_fractures_heat
step_transport_fractures_heat.RdSudicky 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
See also
Other transport:
step_transport_fractures_solute(),
step_transport_ogata_banks()
Examples
formula <- as.formula(~time+z+x)
dat <- as.data.frame(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()