heavy17 {iapws} | R Documentation |
IAPWS-17 Formulation in Single Phase Regions
Description
Compute thermodynamic and transport properties of heavy water and steam using the IAPWS formulation 2017.
Usage
heavy17(what, p, t, rho, state = NULL)
Arguments
what |
a character vector listing the output properties. See Details for available properties. |
p |
a numeric vector giving the pressure values in MPa. |
t |
a numeric vector giving the temperature values in K. |
rho |
a numeric vector giving the density values in kg/m |
state |
|
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.
A valid pair of p
, rho
and t
is needed:
either (p, t)
or (rho, t)
.
Value
A numeric matrix of dimension c(n, length(what))
with n
the maximum length of either rho
, p
, t
or state
.
Note
Computing several properties in a single call may be more efficient than separately because most of the computation time is shared.
Author(s)
Jonathan Debove
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.
International Association for the Properties of Water and Steam, IAPWS R17-20, Release on the IAPWS Formulation 2020 for the Viscosity of Heavy Water.
International Association for the Properties of Water and Steam, IAPWS R18-21, Release on the IAPWS Formulation 2021 for the Thermal Conductivity of Heavy Water.
Examples
# IAPWS-17 Tab. 7
M <- 20.027508 # g/mol
tab <- heavy17(c("p", "cv", "w", "s"),
t = c(300, 300, 300,
500, 500, 500, 500, 500,
643.8,
800, 800),
rho = c(0.55126e2, 0.6e2, 0.65e2,
0.5e-1, 0.5, 0.4626e2, 0.5e2, 0.6e2,
0.2e2,
0.1e-1, 0.25) * M)
tab[, c("cv", "s")] <- tab[, c("cv", "s")] / M
format(tab, scientific = TRUE, digits = 9)
# Viscosity Tab. 3
tab <- heavy17("eta", t = c(298.15, 298.15, 298.15, 373.15,
775.0, 775.0, 775.0),
rho = c(0, 1105, 1130, 1064, 1, 100, 400))
format(tab, digits = 8)
# Viscosity Tab. 4
tab <- heavy17("eta", t = 644.101, rho = c(145, 245, 295, 345, 395, 445))
format(tab, digits = 8)
# Thermal conductivity Tab. 3
tab <- heavy17("lambda", t = c(298.15, 298.15, 298.15, 825),
rho = c(0, 1104.5, 1200, 0))
format(tab, digits = 6)
# Thermal conductivity Tab. 4
tab <- heavy17("lambda", t = 644.10,
rho = c(1, 106, 256, 306, 356, 406, 456, 750))
format(tab, digits = 6)