2.3 Gas & Vapor Property Evaluation: Steam Tables, Superheat & Subcooling

2.3 Gas & Vapor Property Evaluation: Steam Tables, Superheat & Subcooling

Accurate evaluation of thermodynamic fluid properties is fundamental to sizing HVAC equipment, steam distribution networks, and refrigeration circuits. Whether determining the enthalpy drop across a steam turbine, the required subcooling before a thermostatic expansion valve, or the density of moist air at high elevation, mechanical engineers must navigate saturated, superheated, subcooled, and ideal/real gas states with complete precision.

---

1. Pure Substance Phase Change & Property Diagrams

A pure substance has a fixed, homogeneous chemical composition throughout (e.g., water, nitrogen, pure refrigerants). Phase equilibrium is represented on thermodynamic coordinate diagrams including $P-v$, $T-v$, $T-s$, and $P-h$ (pressure-enthalpy) diagrams.

Pressure (P) ^             Critical Point (CP)
             |                 .---.
             |   Subcooled    /     \   Superheated
             |    Liquid     /       \    Vapor
             |    Region    / Two-    \   Region
             |             /  Phase    \
             |   Saturated/   Vapor     \Saturated
             |     Liquid/     Dome      \Vapor
             |       Line/                \Line
             +----------+------------------+--------> Enthalpy (h)

The Thermodynamic Phase Regions

1. Subcooled (Compressed) Liquid: Liquid whose temperature is below the saturation temperature at the given pressure ($T < T_{sat}(P)$), or whose pressure is above the saturation pressure at the given temperature ($P > P_{sat}(T)$).
2. Saturated Liquid ($f$): Liquid on the verge of vaporizing ($x = 0.0$). Any heat addition initiates phase change at constant temperature and pressure.
3. Saturated Liquid-Vapor Mixture (Two-Phase Region): Coexisting liquid and vapor phases in thermodynamic equilibrium ($0 < x < 1.0$). Temperature and pressure are mutually dependent ($P = P_{sat}(T)$).
4. Saturated Vapor ($g$): Vapor on the verge of condensing ($x = 1.0$). Any heat removal initiates condensation at constant temperature and pressure.
5. Superheated Vapor: Vapor whose temperature is above the saturation temperature at the given pressure ($T > T_{sat}(P)$), or whose pressure is below the saturation pressure at the given temperature ($P < P_{sat}(T)$).
6. Supercritical Fluid: State at pressures and temperatures exceeding the Critical Point ($P > P_c, T > T_c$), where liquid and vapor phases merge into a single homogeneous fluid with no distinct phase transition boundary.

---

2. Vapor Quality & Two-Phase Property Evaluation

Inside the saturation dome, temperature and pressure are not independent properties. Specifying pressure automatically fixes saturation temperature ($T_{sat}$), and vice versa. An additional intensive property—such as vapor quality ($x$)—is required to fix the state.

Vapor Quality ($x$)

$x = \frac{m_g}{m_{total}} = \frac{m_g}{m_f + m_g} \quad (0 \le x \le 1.0)$

Where:

• $m_g = \text{mass of saturated vapor } (\text{lb}_m \text{ or kg})$
• $m_f = \text{mass of saturated liquid } (\text{lb}_m \text{ or kg})$
• For saturated liquid: $x = 0.0$
• For saturated vapor: $x = 1.0$
• Quality is undefined in the single-phase compressed liquid and superheated vapor regions.

Two-Phase Mixture Property Formula

Any specific property $y$ ($v, u, h, s$) of a saturated two-phase mixture is computed from saturated liquid ($y_f$) and saturated vapor ($y_g$) values:

$y = (1 - x) y_f + x y_g = y_f + x (y_g - y_f) = y_f + x y_{fg}$

$\text{Specific Volume: } v = v_f + x v_{fg} = v_f + x (v_g - v_f)$

$\text{Specific Enthalpy: } h = h_f + x h_{fg} = h_f + x (h_g - h_f)$

$\text{Specific Entropy: } s = s_f + x s_{fg} = s_f + x (s_g - s_f)$

$\text{Specific Internal Energy: } u = u_f + x u_{fg} = u_f + x (u_g - u_f)$

To determine quality from a known mixture enthalpy $h$:

$x = \frac{h - h_f}{h_{fg}} = \frac{h - h_f}{h_g - h_f}$

---

3. Superheat & Subcooling in HVAC/R Systems

In single-phase regions, temperature and pressure are independent properties. The degree of departure from the saturation boundary is quantified by superheat and subcooling:

+-----------------------------------------------------------------------------------------+
| SUPERHEAT AND SUBCOOLING DEFINITIONS                                                    |
+-----------------------------------------------------------------------------------------+
| Superheat:   Delta_T_sh = T_actual - T_sat(P)    [Evaluated in Superheated Vapor Region]|
| Subcooling:  Delta_T_sc = T_sat(P) - T_actual    [Evaluated in Compressed Liquid Region]|
+-----------------------------------------------------------------------------------------+

Superheat ($\Delta T_{sh}$)

• Definition: The temperature difference between the actual vapor temperature and the saturation (boiling) temperature at the local pressure: $\Delta T_{sh} = T_{actual} - T_{sat}(P)$
• HVAC/R Importance: Evaporator superheat ($8^\circ\text{F} - 15^\circ\text{F}$) ensures that 100% dry vapor enters the compressor suction port, preventing catastrophic liquid droplet slugging and oil dilution.

Subcooling ($\Delta T_{sc}$)

• Definition: The temperature reduction of a liquid below its saturation (condensing) temperature at the local pressure: $\Delta T_{sc} = T_{sat}(P) - T_{actual}$
• HVAC/R Importance: Condenser subcooling ($10^\circ\text{F} - 15^\circ\text{F}$) ensures that solid liquid refrigerant enters the thermostatic expansion valve (TXV), preventing premature flash gas formation in liquid lines and expanding net evaporator cooling capacity ($q_e = h_1 - h_4$).

---

4. Compressed Liquid Region & Incompressible Approximations

Water and liquid refrigerants are nearly incompressible fluids. Because compressed liquid property tables are tabulated only at very high pressures, properties in standard HVAC conditions ($P < 500 \text{ psia}$) are evaluated using saturated liquid properties at the local fluid temperature:

$v(T, P) \approx v_f(T)$

$u(T, P) \approx u_f(T)$

$s(T, P) \approx s_f(T)$

$h(T, P) \approx h_f(T) + v_f(T) [P - P_{sat}(T)] \approx h_f(T)$

For high-pressure boiler feed pumps ($P \ge 500 \text{ psia}$), the flow work correction term $v_f(T)[P - P_{sat}(T)]$ becomes significant ($1 \text{ to } 5 \text{ Btu/lb}_m$) and must be added:

$h(T, P) = h_f(T) + v_f(T) [P - P_{sat}(T)] \times \left(\frac{144 \text{ in}^2/\text{ft}^2}{778.17 \text{ ft}\cdot\text{lb}_f/\text{Btu}}\right)$

---

5. Ideal Gas Equations of State & Gas Mixtures

Gases at temperatures well above their critical temperature ($T \gg T_c$) and pressures far below their critical pressure ($P \ll P_c$) follow the Ideal Gas Law:

$P v = R T \quad \implies \quad P V = m R T = n \bar{R} T$

Where:

• $\bar{R} = 1545.35 \text{ ft}\cdot\text{lb}_f/(\text{lbmol}\cdot^\circ\text{R}) = 8.31447 \text{ kJ}/(\text{kmol}\cdot\text{K}) = \text{Universal Gas Constant}$
• $R = \frac{\bar{R}}{M} = \text{Specific Gas Constant of the gas } (\text{ft}\cdot\text{lb}_f/(\text{lb}_m\cdot^\circ\text{R}) \text{ or } \text{kJ}/(\text{kg}\cdot\text{K}))$
• For standard dry air ($M = 28.966 \text{ lb}m/\text{lbmol}$): $R{air} = 53.352 \text{ ft}\cdot\text{lb}_f/(\text{lb}_m\cdot^\circ\text{R}) = 0.06855 \text{ Btu}/(\text{lb}_m\cdot^\circ\text{R})$

Specific Heat Relationships for Ideal Gases

For ideal gases, internal energy and enthalpy depend strictly on temperature ($u = u(T)$, $h = h(T)$):

$\Delta u = u_2 - u_1 = \int_{T_1}^{T_2} c_v(T) \, dT \approx c_v (T_2 - T_1)$

$\Delta h = h_2 - h_1 = \int_{T_1}^{T_2} c_p(T) \, dT \approx c_p (T_2 - T_1)$

$\text{Mayer's Relation: } c_p - c_v = R$

$\text{Specific Heat Ratio: } k = \frac{c_p}{c_v} \implies c_p = \frac{k R}{k - 1}, \quad c_v = \frac{R}{k - 1}$

---

6. Real Gas Behavior & Compressibility Factor ($Z$)

When gases approach the saturation curve or operate at high pressures and cryogenic temperatures, intermolecular attractive and repulsive forces cause significant deviations from ideal behavior. The Compressibility Factor ($Z$) quantifies this deviation:

$Z = \frac{P v}{R T} = \frac{v_{actual}}{v_{ideal}}$

• When $Z = 1.0$: the fluid behaves as an exact ideal gas.
• When $Z < 1.0$ (dominant intermolecular attractions) or $Z > 1.0$ (molecular volume exclusions), real gas charts or equations of state (Peng-Robinson, Redlich-Kwong, REFPROP) must be utilized.

Principle of Corresponding States

Gases at the same Reduced Pressure ($P_r$) and Reduced Temperature ($T_r$) exhibit approximately the same compressibility factor $Z$ on Nelson-Obert Generalized Compressibility Charts:

$P_r = \frac{P}{P_c}, \quad T_r = \frac{T}{T_c}, \quad v'_r = \frac{v_{actual}}{R T_c / P_c}$

Where $P_c$ and $T_c$ are the substance critical pressure and critical temperature in absolute units.

7. Step-by-Step Worked Example: Steam Heating Coil State Analysis

Problem Statement

A central heating coil in an industrial HVAC air handler receives superheated steam at $P_1 = 100.0\text{ psia}$ and $T_1 = 400.0^\circ\text{F}$ at a steady mass flow rate of $\dot{m} = 1,500\text{ lbm/hr}$. The steam condenses and cools as it passes through the coil, exiting as subcooled liquid condensate at $P_2 = 90.0\text{ psia}$ and $T_2 = 180.0^\circ\text{F}$.

From NCEES steam property tables:

• At $100.0\text{ psia}$: $T_{sat} = 327.86^\circ\text{F}$, $h_f = 298.61\text{ Btu/lbm}$, $h_{fg} = 889.2\text{ Btu/lbm}$, $h_g = 1187.8\text{ Btu/lbm}$.
• Superheated Steam at $100.0\text{ psia}$ and $400.0^\circ\text{F}$: $h_1 = 1227.5\text{ Btu/lbm}$, $v_1 = 4.934\text{ ft}^3/\text{lbm}$, $s_1 = 1.6517\text{ Btu/(lbm}\cdot^\circ\text{R)}$.
• At $90.0\text{ psia}$: $T_{sat} = 320.27^\circ\text{F}$, $h_f = 290.73\text{ Btu/lbm}$, $h_{fg} = 895.0\text{ Btu/lbm}$, $h_g = 1185.7\text{ Btu/lbm}$.
• Saturated Liquid Water at $180.0^\circ\text{F}$: $P_{sat} = 7.515\text{ psia}$, $h_f = 147.99\text{ Btu/lbm}$, $v_f = 0.01651\text{ ft}^3/\text{lbm}$.

Calculate:

1. The degrees of superheat at the coil inlet ($\Delta T_{sh,1}$).
2. The degrees of subcooling at the coil exit ($\Delta T_{sc,2}$).
3. The exit specific enthalpy ($h_2$) using compressed liquid evaluation.
4. The total rate of thermal heat transfer released to the building supply air stream ($\dot{Q}_{coil}$) in $\text{Btu/hr}$ and $\text{MBH}$.

Solution Steps

Step 1: Determine inlet degrees of superheat At $P_1 = 100.0\text{ psia}$, $T_{sat,1} = 327.86^\circ\text{F}$. $\Delta T_{sh,1} = T_1 - T_{sat,1} = 400.0^\circ\text{F} - 327.86^\circ\text{F} = 72.14^\circ\text{F}$

Step 2: Determine exit degrees of subcooling At $P_2 = 90.0\text{ psia}$, $T_{sat,2} = 320.27^\circ\text{F}$. The exit temperature is $T_2 = 180.0^\circ\text{F}$. $\Delta T_{sc,2} = T_{sat,2} - T_2 = 320.27^\circ\text{F} - 180.0^\circ\text{F} = 140.27^\circ\text{F}$

Step 3: Evaluate exit specific enthalpy ($h_2$) Using the incompressible liquid approximation evaluated at local temperature $T_2 = 180.0^\circ\text{F}$ with the flow work correction term: $h_2 = h_f(180^\circ\text{F}) + v_f(180^\circ\text{F}) \times [P_2 - P_{sat}(180^\circ\text{F})] \times \left(\frac{144}{778.17}\right)$ $h_2 = 147.99 + 0.01651 \times [90.0 - 7.515] \times 0.18505$ $h_2 = 147.99 + 0.01651 \times 82.485 \times 0.18505 = 147.99 + 0.252 = 148.24\text{ Btu/lbm}$ (Note: Saturated liquid approximation $h_2 \approx h_f(180^\circ\text{F}) = 147.99\text{ Btu/lbm}$ is accurate within $0.17%$). We use $h_2 = 148.24\text{ Btu/lbm}$.

Step 4: Calculate total heat transfer duty $\dot{Q}_{coil}$ Applying the steady-state open system energy balance: $\dot{Q}_{coil} = \dot{m} (h_1 - h_2)$ $\dot{Q}_{coil} = 1,500\text{ lbm/hr} \times (1227.5 - 148.24\text{ Btu/lbm})$ $\dot{Q}_{coil} = 1,500 \times 1079.26 = 1,618,890\text{ Btu/hr}$

Converting to $\text{MBH}$ ($1\text{ MBH} = 1,000\text{ Btu/hr}$): $\dot{Q}_{coil} = \frac{1,618,890}{1,000} = 1618.9\text{ MBH}$

This validates that steam heating coils yield exceptionally high thermal heat duty primarily due to latent heat release during phase change ($h_{fg} \approx 895\text{ Btu/lbm}$).