R/reynolds_valve.R

In ratral/hyd4gpv: General Purpose Valve functions

#' Valve Reynolds number
#' @param flow flow in m³/h
#' @param kv Flow coefficient in m³/h
#' @param dn diameter in millimeter (mm).
#' @param d1 Inlet diameter reducer only in millimeter (mm).
#' @param d2 Outlet diameter increase only in millimeter (mm).
#' @param temp temperature is in Celsius.
#' @param fl Liquid pressure recovery factor of a control valve without attached fittings
#' @param fd Valve style modifier (1.0 fuer Ventile mit V-fourmigen Drosselquerschnitt)
#' @return valve Reynolds number
#' @export
#'
reynolds_valve <- function( flow, kv, dn, d1, d2, temp = 20, fl, fd){

  Fp <- fp(kv, dn, d1, d2)  # Fp Piping geometry factor

  n4 <- 7.07e4
  n2 <- 0.0016

  kvisc <- kinematic_viscosity(temp) # in mm2/s

  reynolds <- (n4*fd*flow)/(kvisc*sqrt(Fp*fl*kv))*
    ((Fp^2*fl^2*kv^2)/(n2*dn^4)+1)^(1/4)

  return(reynolds)

}