What the calculator does
Water and steam properties follow the industrial formulation IAPWS-IF97 (International Association for the Properties of Water and Steam, Revised Release R7-97(2012)), the standard behind modern steam tables and Mollier charts. The coefficients were transferred from the official IAPWS text, and the code is checked against the verification tables of the formulation itself: 72 property values (18 each in tables 5, 15, 33 and 42 for regions 1, 2, 3 and 5) and 6 saturation points (tables 35 and 36) match in all 9 significant digits.
IAPWS-IF97 regions
- Region 1: liquid water from 0 to 350 °C above the saturation pressure, up to 100 MPa. Gibbs free energy g(p, T).
- Region 2: steam from 0 to 800 °C below the saturation pressure (and below the B23 boundary at 350…590 °C). Gibbs free energy with ideal-gas and residual parts.
- Region 3: the near-critical zone (350…590 °C, high pressures). Helmholtz free energy f(ρ, T); density from p and T is found iteratively.
- Region 4: the saturation line from 0 °C to the critical point (373.946 °C, 22.064 MPa), equations p_s(T) and T_s(p). Wet steam is a mixture of saturated liquid and vapor with quality x.
- Region 5: steam from 800 to 2000 °C up to 50 MPa.
Backward problems
The IF97 basic equations use (p, T) or (ρ, T). For (p, h) and (p, s) the temperature is found by iteration (regula falsi with bisection); for (h, s) the pressure is found too, since at constant s the enthalpy grows monotonically with pressure. IF97 also provides fast backward equations T(p, h) and T(p, s), but they differ from the basic equations by hundredths of a kelvin; iteration gives an answer consistent with the basic equations.
Processes on the h-s chart
- Turbine expansion. The ideal process is the isentrope 1–2s: s₂s = s₁, enthalpy drop H₀ = h₁ − h₂s. With isentropic efficiency η: h₂ = h₁ − η·H₀, work l = h₁ − h₂, power N = G·l.
- Throttling (valve): h₂ = h₁, a horizontal line on the h-s chart; pressure falls, entropy rises, superheated steam usually cools.
Example
Steam at 10 MPa, 500 °C: h₁ = 3375.1 kJ/kg, s₁ = 6.5993 kJ/(kg·K). Expansion to 10 kPa: the isentropic end is wet, x₂s = 0.793, h₂s = 2089.6 kJ/kg, H₀ = 1285.4 kJ/kg. With η = 0.85: h₂ = 2282.5 kJ/kg, x₂ = 0.874, l = 1092.6 kJ/kg; at G = 100 kg/s the power is about 109 MW.
Limitations
Accuracy is that of IF97 itself (a fraction of a percent or better for most properties). Metastable states (subcooled vapor, superheated liquid) and ice are not covered. Very close to the critical point (within fractions of a kelvin of 373.946 °C) the two-phase values v′, v″ come from the region 3 equation and are less accurate. Viscosity, thermal conductivity and surface tension are separate IAPWS formulations and are not included.
Sources: IAPWS R7-97(2012), Revised Release on the IAPWS Industrial Formulation 1997 for the Thermodynamic Properties of Water and Steam, iapws.org; W. Wagner, A. Kruse, Properties of Water and Steam, Springer, 1998.