Calculate thermodynamic properties of water following the IAPWS-95 formulation.

1 | ```
IAPWS95(property, T = 298.15, rho = 1000)
``` |

`property` |
character, name(s) of property(s) to calculate |

`T` |
numeric, temperature (K) |

`rho` |
numeric, density (kg m |

`IAPWS95`

provides an implementation of the IAPWS-95 formulation for properties (including pressure) calculated as a function of temperature and density.

The `IAPWS95`

function returns values of thermodynamic properties in specific units (per gram).
The IAPWS-95 formulation follows the triple point convention used in engineering (values of internal energy and entropy are taken to be zero at the triple point).

For `IAPWS95`

the upper temperature limit of validity is 1000 *°*C, but extrapolation to much higher temperatures is possible (Wagner and Pruss, 2002).
Valid pressures are from the greater of zero bar or the melting pressure at temperature to 10000 bar (with the provision for extrapolation to more extreme conditions).
The function does not check these limits and will attempt calculations for any range of input parameters, but may return `NA`

for properties that fail to be calculated at given temperatures and pressures and/or produce warnings or even errors when problems are encountered.

A data frame the number of rows of which corresponds to the number of input temperature, pressure and/or density values.

Wagner, W. and Pruss, A. (2002) The IAPWS formulation 1995 for the thermodynamic properties of ordinary water substance for general and scientific use. *J. Phys. Chem. Ref. Data* **31**, 387–535. http://dx.doi.org/10.1063/1.1461829

`util.water`

for properties along the saturation curve (`WP02.auxiliary`

) and calculation of density from pressure and temperature (`rho.IAPWS92`

).
`water.IAPWS95`

is a wrapper around `IAPWS95`

and the utility functions, which converts the specific units to molar quantities, and is used in higher-level functions (`water`

).

1 2 | ```
## calculate pressure for given temperature, density
P <- as.numeric(IAPWS95("P", T=500, rho=838.0235))
``` |

Questions? Problems? Suggestions? Tweet to @rdrrHQ or email at ian@mutexlabs.com.

Please suggest features or report bugs with the GitHub issue tracker.

All documentation is copyright its authors; we didn't write any of that.