Phase Change Heat Transfer: Boiling Curves, Critical Heat Flux, and Condensation Modes
When a liquid turns to vapor (boiling) or vapor turns back to liquid (condensation), it absorbs or releases large amounts of heat — and how that happens depends on temperature, pressure, and surface conditions.
⚠️ Why It Matters
📘 Definition
Phase change heat transfer refers to the exchange of thermal energy during first-order phase transitions—primarily boiling and condensation—governed by latent enthalpy, interfacial dynamics, and hydrodynamic stability. It is characterized by non-linear heat flux–temperature difference relationships, critical thresholds (e.g., CHF), and distinct regimes (e.g., nucleate, transition, film boiling; dropwise vs. film condensation). Accurate prediction requires coupling thermodynamics, fluid mechanics, and surface wettability physics.
🎨 Concept Diagram
AI-generated illustration for visual understanding
💡 Engineering Insight
CHF isn’t a material property—it’s a system response. A 5% reduction in inlet subcooling or 2% increase in mass flux can shift CHF by ±12%, yet most field failures trace to unmodeled inlet plenum asymmetry or localized fouling—not the correlation itself. Always anchor design margins to *measured* CHF under representative hydraulic and thermal boundary conditions—not textbook curves.
📖 Detailed Explanation
Beyond CHF lies transition boiling—a unstable, low-efficiency zone where vapor blankets the surface intermittently—followed by stable film boiling where vapor fully insulates the heater. In forced convection (e.g., inside pipes), flow patterns dominate: bubbly → slug → annular → mist flow, each with distinct CHF mechanisms (e.g., liquid film dryout in annular flow vs. bubble crowding in low-velocity bubbly flow).
Advanced treatment requires resolving interfacial physics: bubble nucleation depends on cavity radius distribution and contact angle hysteresis; CHF prediction now incorporates wall superheat gradients, microlayer evaporation, and vapor recoil momentum; modern condensation modeling uses statistical droplet growth/merging algorithms and dynamic contact angle models validated against high-speed IR thermography—going far beyond classical Nusselt film theory.
🔄 Engineering Workflow
📋 Decision Guide
| Rock/Field Condition | Recommended Design Action |
|---|---|
| High-pressure water system (>10 MPa) with tight thermal margins | Use subcooled nucleate boiling design with CHF margin ≥15%; specify surface-enhanced tubes (e.g., sintered metal powder) and enforce minimum ΔTsub ≥ 15 K |
| Low-surface-tension fluid (e.g., refrigerant R-134a) in compact heat exchanger | Select microfin or herringbone-tube geometry; operate in high-quality annular flow regime; avoid CHF-prone low-mass-flux zones (<200 kg/m²·s) |
| Condensing steam on vertical tube bundle with fouling risk | Apply hydrophobic nano-coating (θ ≈ 105°) to promote dropwise condensation; include periodic steam-jet cleaning ports; size drain lines for intermittent slug flow |
📊 Key Properties & Parameters
Critical Heat Flux (CHF)
0.1–2.5 MW/m² (subcooled water at 7–15 MPa, typical PWR/BWR conditions)Maximum heat flux sustainable before departure from nucleate boiling (DNB) or dryout, beyond which heat transfer deteriorates catastrophically.
Defines absolute safety limit for core power density and dictates minimum mass velocity and subcooling requirements.
Boiling Number (Bo)
1×10⁻⁵ to 5×10⁻³ (vertical upflow in industrial boilers and nuclear steam generators)Dimensionless ratio of convective heat flux to latent heat transport capacity: Bo = q'' / (G h_fg x), where x is quality.
Predicts onset of flow instability and dryout location in two-phase evaporators and fuel assemblies.
Surface Roughness (Ra)
0.1–10 µm (machined copper, polished stainless steel, enhanced surfaces like micro-structured or porous coatings)Arithmetic average deviation of surface profile from its mean line, governing nucleation site density in pool boiling.
Directly controls nucleate boiling heat transfer coefficient—rougher surfaces increase CHF up to an optimal Ra (~1–3 µm), then degrade due to bubble coalescence.
Contact Angle (θ)
10°–120° (hydrophilic: θ < 90°; hydrophobic: θ > 90°; e.g., bare Cu ≈ 65°, silanized Cu ≈ 110°)Angle formed between liquid–vapor interface and solid surface, quantifying wettability and influencing condensation mode.
Determines whether condensation proceeds as high-efficiency dropwise (θ > 90°) or low-efficiency filmwise (θ < 70°), impacting heat transfer coefficients by 3–8×.
📐 Key Formulas
Kutateladze CHF Correlation (Pool Boiling)
q''_CHF = C σ g (ρ_l − ρ_v)^{0.5} h_fgEmpirical correlation for maximum heat flux in saturated pool boiling
| Symbol | Name | Unit | Description |
|---|---|---|---|
| q''_CHF | Critical Heat Flux | W/m² | Maximum heat flux in saturated pool boiling |
| C | Empirical Constant | dimensionless | Correlation-specific constant dependent on fluid and surface conditions |
| σ | Surface Tension | N/m | Interfacial tension between liquid and vapor phases |
| g | Gravitational Acceleration | m/s² | Acceleration due to gravity |
| ρ_l | Liquid Density | kg/m³ | Density of the saturated liquid phase |
| ρ_v | Vapor Density | kg/m³ | Density of the saturated vapor phase |
| h_fg | Latent Heat of Vaporization | J/kg | Enthalpy required to vaporize unit mass of liquid at saturation conditions |
Chen Boiling Heat Transfer Coefficient
h = h_conv + h_nuc = F·h_{conv, single-phase} + S·h_{nuc, pool}Semi-empirical model combining convective and nucleate components for flow boiling
| Symbol | Name | Unit | Description |
|---|---|---|---|
| h | Boiling heat transfer coefficient | W/(m²·K) | Total heat transfer coefficient for flow boiling |
| h_conv | Convective heat transfer coefficient | W/(m²·K) | Component due to forced convection |
| h_nuc | Nucleate boiling heat transfer coefficient | W/(m²·K) | Component due to nucleate boiling |
| F | Convective enhancement factor | dimensionless | Multiplier for single-phase convective coefficient |
| h_{conv, single-phase} | Single-phase convective heat transfer coefficient | W/(m²·K) | Convective coefficient without boiling effects |
| S | Nucleate boiling enhancement factor | dimensionless | Multiplier for pool nucleate boiling coefficient |
| h_{nuc, pool} | Pool nucleate boiling heat transfer coefficient | W/(m²·K) | Nucleate boiling coefficient under pool boiling conditions |
🏭 Engineering Example
Vogtle Electric Generating Plant Unit 3 (Georgia, USA)
N/A — nuclear steam generator application🏗️ Applications
- Pressurized Water Reactor (PWR) fuel assembly thermal-hydraulics
- Waste heat recovery from data center immersion coolants
- Cryogenic propellant management in upper-stage rocket tanks
📋 Real Project Case
Air-Cooled Condenser Retrofit for 600 MW Coal Power Plant
Retrofit of legacy water-cooled condenser at Midwest US plant