dyn_visc_water: Absolute or Dynamic Viscosity for Liquid Water

In iemisc: Irucka Embry's Miscellaneous Functions

Absolute or Dynamic Viscosity for Liquid Water

Description

This function solves for the absolute or dynamic viscosity of water using only the temperature of the water in either units of degrees Celsius, degrees Fahrenheit, or Kelvin.

Usage

dyn_visc_water(
  Temp,
  units = c("SI", "Eng", "Absolute"),
  Eng_units = c("slug/ft/s", "lbf*s/ft^2")
)

Arguments

Temp	numeric vector that contains the temperature (degrees Celsius, degrees Fahrenheit, or Kelvin)
units	character vector that contains the system of units options are SI for International System of Units, Eng for English units (United States Customary System in the United States and Imperial Units in the United Kingdom), or Absolute for Absolute Units
Eng_units	character vector that contains the unit for the dynamic viscosity of water in the English system options are slug/ft/s or lbf*s/ft^2

Details

The simplified equation is expressed as

\\mu_s = \frac{1}{a + bT + cT^2 + dT^3}

with

a = 557.82468

b = 19.408782

c = 0.1360459

d = -3.1160832 * 10 ^ -4

\mu_s

Water Absolute or Dynamic Viscosity (kg/m*s, slug/ft/s, or lbf*s/ft^2)

T

the water temperature, degrees Celsius

Value

the absolute or dynamic viscosity as a numeric vector. The units are not returned.

Author(s)

Irucka Embry

References

C. O. Popiel & J. Wojtkowiak (1998). "Simple Formulas for Thermophysical Properties of Liquid Water for Heat Transfer Calculations (from 0C to 150C)". Heat Transfer Engineering, 19:3, 87-101, article from ResearchGate: https://www.researchgate.net/publication/239243539_Simple_Formulas_for_Thermophysical_Properties_of_Liquid_Water_for_Heat_Transfer_Calculations_from_0C_to_150C.

Examples

# Example 1 (Compare to the tabulated values in the Reference paper)

install.load::load_package("iemisc", "data.table", "round")

Temp <- c(0, 0.01, 3.86, seq(5, 95, by = 5), 99.974, seq(100, 150, by = 5))

dynamic_viscosity <- data.table("Temperature (degrees C)" = Temp,
"mu (* 10 ^ 6, kg / m*s)" = round_r3(dyn_visc_water(Temp, units = "SI")
* 10^6, d = 1))
dynamic_viscosity





# Example 2 - Example from the hydraulics package

library(iemisc)

mu <- hydraulics::dvisc(T = 55, units = "Eng"); mu

mu2 <- dyn_visc_water(Temp = 55, units = "Eng", Eng_units = "lbf*s/ft^2"); mu2

iemisc documentation built on Sept. 25, 2023, 5:09 p.m.