R/cavtation_regime.R

In ratral/hyd4gpv: General Purpose Valve functions

#' Cavitation Level (Regime)
#'
#' @param x valve position
#' @param b steepness
#' @param d upper value
#' @param e the effective dose
#' @param fls liquid pressure recovery full open (max between fl and Flp/Fp)
#' @param sigma_value cavitation index (Reference downstream pressure P2)
#'
#' @return regime
#' @export
#'
cavtation_regime <- function(x, b, d, e, fls, sigma_value) {
  limit_1 <- sigma_i(x, b, d, e, fls)
  limit_2 <- sigma_c(x, b, d, e, fls)
  limit_3 <- sigma_mv(x, b, d, e, fls)

  regime <- dplyr::case_when(
    sigma_value > limit_1 ~ "Regime 1: Cavitation free",
    sigma_value < limit_1 & sigma_value > limit_2 ~ "Regime 2: Incipient cavitation.",
    sigma_value < limit_2 & sigma_value > limit_3 ~ "Regime 3: Constant cavitation.",
    sigma_value < limit_3 ~ "Regime 4: Maximum vibration cavitation."
  )

  return(regime)
}