R/permeability_to_hydraulic_conductivity.R

Defines functions permeability_to_hydraulic_conductivity

Documented in permeability_to_hydraulic_conductivity

#' Convert permeability to hydraulic conductivity
#' 
#' Converts an intrinsic permeability value to a value for hydraulic conductivity.
#' 
#' @param k Intrinsic permeability of the medium (in square metres)
#' @param mu Dynamic viscosity of the fluid; defaults to 1.308E-3 Pa*s
#' @param rho Density of the fluid; defaults to 999.7 kg/m³
#' @param g Gravitational acceleration; defaults to 9.80665 m/s²
#' @param t Temperature; defaults to 10°C
#' 
#' @return Hydraulic conductivity value
#' @details \code{mu} and \code{rho} are automatically determined in function of t in a temperature range between 0 and 100 °C, unless the values are explicitely specified.
#' @seealso \code{\link{hydraulic_conductivity_to_permeability}}
#' @export
permeability_to_hydraulic_conductivity <- function(k,mu=1.308E-3,rho=999.7,g=9.80665,t=10)
{
  if(t!=10 & mu==1.308E-3 & rho==999.7)
  {
    tVector <- c(0,5,10,20,30,40,50,60,70,80,90,100)
    rhoVector <- c(999.8,1000,999.7,998.2,995.7,992.2,988.1,983.2,977.8,971.8,965.3,958.4)
    muVector <- c(1.787,1.519,1.307,1.002,0.798,0.653,0.547,0.467,0.404,0.355,0.315,0.282)*10^(-3)
    mu <- approx(x=tVector,y=muVector,xout=t)$y
    rho <- approx(x=tVector,y=rhoVector,xout=t)$y
  }
  return(k*(rho*g)/mu)
}
rogiersbart/groundwaterr documentation built on May 27, 2019, 1:10 p.m.