R/step_transport_fractures_heat.R

In jkennel/hydrorecipes: Hydrogeology steps

#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
#
# Sudicky and Frind 1982 Parallel Fractures Adapted for Heat Step --------------
#
#^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^^
StepTransportFracturesHeat <- R6Class(

  classname = 'step_transport_fractures_heat',
  inherit = Step,

  public = list(

    # step specific variables
    time = NULL,
    distance_fracture = NULL,
    distance_matrix = NULL,
    temperature_influent = NULL,
    time_influent = NULL,
    temperature_initial = NULL,
    fracture_aperture = NULL,
    fracture_spacing = NULL,
    velocity = NULL,
    thermal_conductivity_water = NULL,
    thermal_conductivity_solids = NULL,
    specific_heat_water = NULL,
    specific_heat_solids = NULL,
    density_water = NULL,
    density_solids = NULL,
    porosity = NULL,
    n_terms = NULL,

    initialize = function(time,
                          distance_fracture,
                          distance_matrix,
                          temperature_influent = 15.0,
                          time_influent = 0.0,
                          temperature_initial = 10,
                          fracture_aperture = 2e-4,
                          fracture_spacing = 1.0,
                          velocity = 0.1 / 86400.0,
                          thermal_conductivity_water = 0.615,
                          thermal_conductivity_solids = 3.4,
                          specific_heat_water = 4192,
                          specific_heat_solids = 908,
                          density_water = 1.0,
                          density_solids = 2.5,
                          porosity = 0.1,
                          n_terms = 30L,
                          role = "predictor",
                          ...) {

      # get function parameters to pass to parent
      time <- deparse(substitute(time))
      distance_fracture <- deparse(substitute(distance_fracture))
      distance_matrix <- deparse(substitute(distance_matrix))
      env_list <- get_function_arguments()
      env_list$step_name <- 'step_transport_fractures_heat'
      env_list$type <- 'add'
      super$initialize(terms = c(as.symbol(time),
                                 as.symbol(distance_fracture),
                                 as.symbol(distance_matrix)),
                       env_list,
                       ...)

      # step specific values
      self$time = time
      self$distance_fracture = distance_fracture
      self$distance_matrix = distance_matrix
      self$temperature_influent = temperature_influent
      self$time_influent = time_influent
      self$temperature_initial = temperature_initial
      self$fracture_aperture = fracture_aperture
      self$fracture_spacing = fracture_spacing
      self$velocity = velocity
      self$thermal_conductivity_water = thermal_conductivity_water
      self$thermal_conductivity_solids = thermal_conductivity_solids
      self$specific_heat_water = specific_heat_water
      self$specific_heat_solids = specific_heat_solids
      self$density_water = density_water
      self$density_solids = density_solids
      self$porosity = porosity
      self$n_terms = n_terms

      self$columns <- c(time, distance_fracture, distance_matrix)

      invisible(self)
    },

    bake = function(new_data) {

      pfh <- list(parallel_fractures_heat(
        new_data[[1L]], # time
        new_data[[2L]], # z
        new_data[[3L]], # x
        self$temperature_influent,
        self$time_influent,
        self$temperature_initial,
        self$fracture_aperture,
        self$fracture_spacing,
        self$velocity,
        self$thermal_conductivity_water,
        self$thermal_conductivity_solids,
        self$specific_heat_water,
        self$specific_heat_solids,
        self$density_water,
        self$density_solids,
        self$porosity,
        self$n_terms))

      self$new_columns <- self$prefix

      names(pfh) <- self$new_columns

      self$result <- pfh
      self$result

    }
  )
)