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

Key Takeaways

  • In the two-phase vapor dome, thermodynamic state is defined by pressure (or saturation temperature) and vapor quality x=mgmf+mgx = \frac{m_g}{m_f + m_g}; any property yy is computed as y=yf+xyfgy = y_f + x y_{fg}.

  • Superheat is the sensible temperature difference above the boiling point at a given pressure (ΔTsh=Tactual−Tsat(P)\Delta T_{sh} = T_{actual} - T_{sat}(P)); subcooling is the sensible temperature reduction below the boiling point (ΔTsc=Tsat(P)−Tactual\Delta T_{sc} = T_{sat}(P) - T_{actual}).

  • Liquid water and refrigerants in the compressed liquid region can be approximated using saturated liquid properties at the local fluid temperature: h(T,P)≈hf(T)+vf(T)[P−Psat(T)]≈hf(T)h(T, P) \approx h_f(T) + v_f(T)[P - P_{sat}(T)] \approx h_f(T).

  • The ideal gas law Pv=RTP v = R T applies accurately to gases well above their critical temperature or at very low pressures (Pr=P/Pc≪1P_r = P/P_c \ll 1); enthalpy changes for ideal gases depend strictly on temperature: Δh=cpΔT\Delta h = c_p \Delta T.

  • Real gas deviations are quantified by the compressibility factor Z=PvRTZ = \frac{P v}{R T}; when Z≠1.0Z \ne 1.0, generalized compressibility charts based on reduced coordinates (Pr,TrP_r, T_r) must be used.

Last updated: August 2026

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−vP-v, T−vT-v, T−sT-s, and P−hP-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<Tsat(P)T < T_{sat}(P)), or whose pressure is above the saturation pressure at the given temperature (P>Psat(T)P > P_{sat}(T)).
  2. Saturated Liquid (ff): Liquid on the verge of vaporizing (x=0.0x = 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.00 < x < 1.0). Temperature and pressure are mutually dependent (P=Psat(T)P = P_{sat}(T)).
  4. Saturated Vapor (gg): Vapor on the verge of condensing (x=1.0x = 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>Tsat(P)T > T_{sat}(P)), or whose pressure is below the saturation pressure at the given temperature (P<Psat(T)P < P_{sat}(T)).
  6. Supercritical Fluid: State at pressures and temperatures exceeding the Critical Point (P>Pc,T>TcP > 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 (TsatT_{sat}), and vice versa. An additional intensive property—such as vapor quality (xx)—is required to fix the state.

Vapor Quality (xx)

x=mgmtotal=mgmf+mg(0≤x≤1.0)x = \frac{m_g}{m_{total}} = \frac{m_g}{m_f + m_g} \quad (0 \le x \le 1.0)

Where:

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

Two-Phase Mixture Property Formula

Any specific property yy (v,u,h,sv, u, h, s) of a saturated two-phase mixture is computed from saturated liquid (yfy_f) and saturated vapor (ygy_g) values:

y=(1−x)yf+xyg=yf+x(yg−yf)=yf+xyfgy = (1 - x) y_f + x y_g = y_f + x (y_g - y_f) = y_f + x y_{fg} Specific Volume: v=vf+xvfg=vf+x(vg−vf)\text{Specific Volume: } v = v_f + x v_{fg} = v_f + x (v_g - v_f) Specific Enthalpy: h=hf+xhfg=hf+x(hg−hf)\text{Specific Enthalpy: } h = h_f + x h_{fg} = h_f + x (h_g - h_f) Specific Entropy: s=sf+xsfg=sf+x(sg−sf)\text{Specific Entropy: } s = s_f + x s_{fg} = s_f + x (s_g - s_f) Specific Internal Energy: u=uf+xufg=uf+x(ug−uf)\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 hh:

x=h−hfhfg=h−hfhg−hfx = \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 (ΔTsh\Delta T_{sh})

  • Definition: The temperature difference between the actual vapor temperature and the saturation (boiling) temperature at the local pressure: ΔTsh=Tactual−Tsat(P)\Delta T_{sh} = T_{actual} - T_{sat}(P)
  • HVAC/R Importance: Evaporator superheat (8∘F−15∘F8^\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 (ΔTsc\Delta T_{sc})

  • Definition: The temperature reduction of a liquid below its saturation (condensing) temperature at the local pressure: ΔTsc=Tsat(P)−Tactual\Delta T_{sc} = T_{sat}(P) - T_{actual}
  • HVAC/R Importance: Condenser subcooling (10∘F−15∘F10^\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 (qe=h1−h4q_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 psiaP < 500 \text{ psia}) are evaluated using saturated liquid properties at the local fluid temperature:

v(T,P)≈vf(T)v(T, P) \approx v_f(T) ν(T,P)≈uf(T)\nu(T, P) \approx u_f(T) s(T,P)≈sf(T)s(T, P) \approx s_f(T) h(T,P)≈hf(T)+vf(T)[P−Psat(T)]≈hf(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≥500 psiaP \ge 500 \text{ psia}), the flow work correction term vf(T)[P−Psat(T)]v_f(T)[P - P_{sat}(T)] becomes significant (1 to 5 Btu/lbm1 \text{ to } 5 \text{ Btu/lb}_m) and must be added:

h(T,P)=hf(T)+vf(T)[P−Psat(T)]×(144 in2/ft2778.17 ft⋅lbf/Btu)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≫TcT \gg T_c) and pressures far below their critical pressure (P≪PcP \ll P_c) follow the Ideal Gas Law:

Pv=RT  ⟹  PV=mRT=nRˉTP v = R T \quad \implies \quad P V = m R T = n \bar{R} T

Where:

  • Rˉ=1545.35 ft⋅lbf/(lbmol⋅∘R)=8.31447 kJ/(kmol⋅K)=Universal Gas Constant\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=RˉM=Specific Gas Constant of the gas (ft⋅lbf/(lbm⋅∘R) or kJ/(kg⋅K))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 lbm/lbmolM = 28.966 \text{ lb}_m/\text{lbmol}): Rair=53.352 ft⋅lbf/(lbm⋅∘R)=0.06855 Btu/(lbm⋅∘R)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)u = u(T), h=h(T)h = h(T)):

Δu=u2−u1=∫T1T2cv(T) dT≈cv(T2−T1)\Delta u = u_2 - u_1 = \int_{T_1}^{T_2} c_v(T) \, dT \approx c_v (T_2 - T_1) Δh=h2−h1=∫T1T2cp(T) dT≈cp(T2−T1)\Delta h = h_2 - h_1 = \int_{T_1}^{T_2} c_p(T) \, dT \approx c_p (T_2 - T_1) Mayer’s Relation: cp−cv=R\text{Mayer's Relation: } c_p - c_v = R Specific Heat Ratio: k=cpcv  ⟹  cp=kRk−1,cv=Rk−1\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 (ZZ)

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 (ZZ) quantifies this deviation:

Z=PvRT=vactualvidealZ = \frac{P v}{R T} = \frac{v_{actual}}{v_{ideal}}
  • When Z=1.0Z = 1.0: the fluid behaves as an exact ideal gas.
  • When Z<1.0Z < 1.0 (dominant intermolecular attractions) or Z>1.0Z > 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 (PrP_r) and Reduced Temperature (TrT_r) exhibit approximately the same compressibility factor ZZ on Nelson-Obert Generalized Compressibility Charts:

Pr=PPc,Tr=TTc,vr′=vactualRTc/PcP_r = \frac{P}{P_c}, \quad T_r = \frac{T}{T_c}, \quad v'_r = \frac{v_{actual}}{R T_c / P_c}

Where PcP_c and TcT_c are the substance critical pressure and critical temperature in absolute units.

Loading diagram...
Thermodynamic Phase Regions & Vapor Dome on P-h Coordinates

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 P1=100.0 psiaP_1 = 100.0\text{ psia} and T1=400.0∘FT_1 = 400.0^\circ\text{F} at a steady mass flow rate of m˙=1,500 lbm/hr\dot{m} = 1,500\text{ lbm/hr}. The steam condenses and cools as it passes through the coil, exiting as subcooled liquid condensate at P2=90.0 psiaP_2 = 90.0\text{ psia} and T2=180.0∘FT_2 = 180.0^\circ\text{F}.

From NCEES steam property tables:

  • At 100.0 psia100.0\text{ psia}: Tsat=327.86∘FT_{sat} = 327.86^\circ\text{F}, hf=298.61 Btu/lbmh_f = 298.61\text{ Btu/lbm}, hfg=889.2 Btu/lbmh_{fg} = 889.2\text{ Btu/lbm}, hg=1187.8 Btu/lbmh_g = 1187.8\text{ Btu/lbm}.
  • Superheated Steam at 100.0 psia100.0\text{ psia} and 400.0∘F400.0^\circ\text{F}: h1=1227.5 Btu/lbmh_1 = 1227.5\text{ Btu/lbm}, v1=4.934 ft3/lbmv_1 = 4.934\text{ ft}^3/\text{lbm}, s1=1.6517 Btu/(lbm⋅∘R)s_1 = 1.6517\text{ Btu/(lbm}\cdot^\circ\text{R)}.
  • At 90.0 psia90.0\text{ psia}: Tsat=320.27∘FT_{sat} = 320.27^\circ\text{F}, hf=290.73 Btu/lbmh_f = 290.73\text{ Btu/lbm}, hfg=895.0 Btu/lbmh_{fg} = 895.0\text{ Btu/lbm}, hg=1185.7 Btu/lbmh_g = 1185.7\text{ Btu/lbm}.
  • Saturated Liquid Water at 180.0∘F180.0^\circ\text{F}: Psat=7.515 psiaP_{sat} = 7.515\text{ psia}, hf=147.99 Btu/lbmh_f = 147.99\text{ Btu/lbm}, vf=0.01651 ft3/lbmv_f = 0.01651\text{ ft}^3/\text{lbm}.

Calculate:

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

Solution Steps

Step 1: Determine inlet degrees of superheat At P1=100.0 psiaP_1 = 100.0\text{ psia}, Tsat,1=327.86∘FT_{sat,1} = 327.86^\circ\text{F}.

ΔTsh,1=T1−Tsat,1=400.0∘F−327.86∘F=72.14∘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 P2=90.0 psiaP_2 = 90.0\text{ psia}, Tsat,2=320.27∘FT_{sat,2} = 320.27^\circ\text{F}. The exit temperature is T2=180.0∘FT_2 = 180.0^\circ\text{F}.

ΔTsc,2=Tsat,2−T2=320.27∘F−180.0∘F=140.27∘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 (h2h_2) Using the incompressible liquid approximation evaluated at local temperature T2=180.0∘FT_2 = 180.0^\circ\text{F} with the flow work correction term:

h2=hf(180∘F)+vf(180∘F)×[P2−Psat(180∘F)]×(144778.17)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) h2=147.99+0.01651×[90.0−7.515]×0.18505h_2 = 147.99 + 0.01651 \times [90.0 - 7.515] \times 0.18505 h2=147.99+0.01651×82.485×0.18505=147.99+0.252=148.24 Btu/lbmh_2 = 147.99 + 0.01651 \times 82.485 \times 0.18505 = 147.99 + 0.252 = 148.24\text{ Btu/lbm}

(Note: Saturated liquid approximation h2≈hf(180∘F)=147.99 Btu/lbmh_2 \approx h_f(180^\circ\text{F}) = 147.99\text{ Btu/lbm} is accurate within 0.17%0.17\%). We use h2=148.24 Btu/lbmh_2 = 148.24\text{ Btu/lbm}.

Step 4: Calculate total heat transfer duty Q˙coil\dot{Q}_{coil} Applying the steady-state open system energy balance:

Q˙coil=m˙(h1−h2)\dot{Q}_{coil} = \dot{m} (h_1 - h_2) Q˙coil=1,500 lbm/hr×(1227.5−148.24 Btu/lbm)\dot{Q}_{coil} = 1,500\text{ lbm/hr} \times (1227.5 - 148.24\text{ Btu/lbm}) Q˙coil=1,500×1079.26=1,618,890 Btu/hr\dot{Q}_{coil} = 1,500 \times 1079.26 = 1,618,890\text{ Btu/hr}

Converting to MBH\text{MBH} (1 MBH=1,000 Btu/hr1\text{ MBH} = 1,000\text{ Btu/hr}):

Q˙coil=1,618,8901,000=1618.9 MBH\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 (hfg≈895 Btu/lbmh_{fg} \approx 895\text{ Btu/lbm}).

Test Your Knowledge

A steam boiler operates at a drum pressure of 150 psia. A sample of wet steam extracted from the boiler drum has a specific enthalpy of 1,020.0 Btu/lbm. Using steam tables at 150 psia where h_f = 330.75 Btu/lbm, h_fg = 863.60 Btu/lbm, and h_g = 1,194.35 Btu/lbm, what is the dryness quality (x) of this steam and its specific volume if v_f = 0.01809 ft^3/lbm and v_g = 3.015 ft^3/lbm?

A

Quality x = 0.798 (79.8%); Specific volume v = 2.410 ft^3/lbm

B

Quality x = 0.854 (85.4%); Specific volume v = 2.576 ft^3/lbm

C

Quality x = 0.742 (74.2%); Specific volume v = 2.240 ft^3/lbm

D

Quality x = 0.912 (91.2%); Specific volume v = 2.752 ft^3/lbm

Test Your Knowledge

A refrigeration technician measures R-410A suction line conditions at the evaporator outlet: pressure gauge reads 118.0 psig (barometric pressure is 14.7 psia, so absolute pressure is 132.7 psia) and pipe surface temperature is 52.0°F. If the saturation temperature of R-410A at 132.7 psia is 45.0°F, what is the operating evaporator superheat, and what does this measurement verify about fluid entering the compressor?

A

Superheat = 14.7°F; verifies two-phase refrigerant mixture is entering the compressor.

B

Superheat = 7.0°F; verifies 100% dry vapor is entering the compressor, protecting against liquid slugging.

C

Subcooling = 7.0°F; verifies solid liquid refrigerant is entering the compressor.

D

Superheat = 0.0°F; verifies saturated liquid refrigerant is entering the compressor.

Test Your Knowledge

Liquid water at 120°F is pumped into a boiler economizer at a discharge pressure of 800 psia. At 120°F, saturated water properties are: h_f = 88.00 Btu/lbm, v_f = 0.01620 ft^3/lbm, and P_sat = 1.693 psia. What is the estimated specific enthalpy of this compressed liquid water including the Poynting flow work correction term?

A

88.00 Btu/lbm

B

98.40 Btu/lbm

C

90.39 Btu/lbm

D

104.22 Btu/lbm

Test Your Knowledge

A high-pressure nitrogen accumulator vessel stores gas at 3,000 psia and -50°F. For nitrogen, critical properties are P_c = 492 psia and T_c = 227.1°R (-232.55°F). Why must an engineer use a generalized compressibility chart rather than the ideal gas law (P v = R T) for sizing calculations at this state point?

A

Because nitrogen undergoes spontaneous condensation into liquid at any pressure exceeding 1,000 psia.

B

Because specific heats c_p and c_v become equal to zero near the critical point.

C

Because nitrogen dissociates into atomic nitrogen radicals under cryogenic pressures exceeding 2,000 psia.

D

Because the high reduced pressure (P_r = 6.10) and low reduced temperature (T_r = 1.80) produce a compressibility factor Z significantly different from 1.0, where ideal gas equations introduce substantial error.

Sections you finish are checked off in the contents.