iapws (version 1.1)

heavy17_sat: IAPWS-17 Formulation in the Saturation Region

Description

Compute thermodynamic and transport properties of heavy water along the saturated line according to the IAPWS formulation 2017.

Usage

heavy17_sat(what, p, t)

Value

A numeric array of dimension c(n, length(what), 2L) with

n the length of either p or t. The last dimension indicate the physical state ("liquid" or "gas").

Arguments

what

a character vector listing the output properties. See Details for available properties.

t

a numeric vector giving the temperature values in K.

p

a numeric vector giving the pressure values in MPa.

Author

Jonathan Debove

Details

The available properties for what are:

• "f": the specific free energy in kJ/kg.

• "g": the specific Gibbs enthalpy in kJ/kg.

• "u": the specific internal energy in kJ/kg.

• "h": the specific enthalpy in kJ/kg.

• "s": the specific entropy in kJ/K/kg.

• "t": the temperature in K.

• "p": the pressure in MPa.

• "v": the specific volume in m\(^3\)/kg.

• "rho": the mass density in kg/m\(^3\).

• "cv": the specific isochoric heat capacity in kJ/K/kg.

• "cp": the specific isobaric heat capacity in kJ/K/kg.

• "alpha": the expansion coefficient in 1/K

• "beta": the pressure coefficient in MPa/K.

• "kappat": the isothermal compressibility in 1/MPa.

• "w": the speed of sound in m/s.

• "eta": the dynamic viscosity in \(\mu\)Pa.s.

• "lambda": the thermal conductivity in mW/K/m.

Only one of t or p is needed. If both arguments are given, p is ignored.

References

International Association for the Properties of Water and Steam, IAPWS R16-17(2018), Revised Release on the IAPWS Formulation 2017 for the Thermodynamic Properties of Heavy Water.

Examples

# IAPWS-17 Tab. 8
M <- 20.027508  # g/mol
tab <- heavy17_sat(c("p", "rho", "h", "s"), t = c(280, 450, 625))
format(tab, scientific = TRUE, digits = 9)
# \dontshow{
stopifnot(all.equal(tab[, "p", "gas", drop = FALSE],
          array(c(0.823054058e-3, 0.921212105, 0.172118129e2),
                c(3, 1, 1), dimnames = list(NULL, c("p"), c("gas"))),
          tolerance = 2e-9))
# }