R/q_max.R

In ratral/hyd4gpv: General Purpose Valve functions

#' @title the maximum flow through the valve
#' @param p1 Gauge upstream pressure (bar)
#' @param fl liquid pressure recovery factor
#' @param fr Reynolds number factor
#' @param kv Flow coefficient value in (m3/h).
#' @param dn valve diameter (mm).
#' @param d1 downstream pipe diameter (mm).
#' @param d2 upstream pipe diameter (mm).
#' @param masl meters above sea level (m).
#' @param temp The temperature is in Celsius.
#' @return q_max (m3/h)
#'
#' @export
#'
q_max <- function(p1, fl, fr, kv, dn, d1, d2, masl = 0, temp = 20){
  p1 = p1 + atm_pressure(masl)
  flp_value <-  flp(kv, fl, dn, d1, d2)
  ff_value <- ff(temp)
  r_density <- water_density(temp)/water_density(20)
  flow <- kv * flp_value * fr * sqrt((p1 - ff_value * vapour_pressure(temp))/r_density)
  return(flow)
}