Boiling Heat Transfer Analysis

Category: Coupled Analysis – Thermo-Fluid | Revised 2026-10-01
CAE visualization for boiling heat transfer theory - technical simulation diagram
Boiling Heat Transfer Analysis

Theory: the four regions of the boiling curve

Overview

🙋

I hear boiling carries heat very well — how much better? And I've heard that too much boiling is dangerous.

🎓

Bubbles carry away latent heat as they form and stir the liquid as they leave, so nucleate boiling of water reaches 10,000–60,000 W/m²K — over a thousand times natural convection in air. But raise the heat flux and at some point bubbles merge into a vapour film blanketing the surface, and heat transfer drops sharply. That's the critical heat flux (CHF); staying below it is a basic design premise for reactors and electronics cooling. Let's follow the boiling curve with numbers.

The four regions

RegionSuperheat (water, 1 atm)Features
Natural convectionBelow about 5 KNo bubbles
Nucleate boilingAbout 5 to 30 KBest heat transfer
Transition boilingAbout 30 to 120 KHeat flux falls as temperature rises
Film boilingAbove about 120 KHeat passes through a vapour film

Nucleate boiling and CHF correlations

$$ q'' = \mu_l h_{fg}\left[\frac{g(\rho_l-\rho_v)}{\sigma}\right]^{1/2}\left(\frac{c_{p,l}\,\Delta T_e}{C_{sf}\,h_{fg}\,Pr_l^{\,n}}\right)^{3} $$
$$ q''_{max} = 0.149\,h_{fg}\,\rho_v\left[\frac{\sigma g(\rho_l-\rho_v)}{\rho_v^2}\right]^{1/4},\qquad q''_{min} = 0.09\,\rho_v h_{fg}\left[\frac{g\sigma(\rho_l-\rho_v)}{(\rho_l+\rho_v)^2}\right]^{1/4} $$

Top is Rohsenow's nucleate boiling correlation, with heat flux proportional to the cube of superheat $\Delta T_e$; $C_{sf}$ depends on the liquid–surface pair (0.013 for water on copper). Below are Zuber's critical and minimum (Leidenfrost) heat fluxes, independent of surface material.

Coffee Break Trivia

Nukiyama's boiling curve

The boiling curve was first measured by Shiro Nukiyama at Tohoku Imperial University in 1934. He heated platinum and nichrome wires in water with electric current and studied heat flux against wire temperature. As current (heat flux) rose, the wire temperature suddenly jumped at a certain value, glowing red or melting — the discovery of the critical heat flux. Nukiyama reasoned that temperature control would reveal the region in between, and transition boiling was later confirmed. Textbooks worldwide still call it the Nukiyama curve.

Worked example: water at 1 atm

Nucleate boiling heat flux

SuperheatHeat fluxHeat transfer coefficient
5 K0.017 MW/m²3,400 W/m²K
10 K0.14 MW/m²13,700 W/m²K
15 K0.46 MW/m²30,800 W/m²K
20 K1.10 MW/m²54,800 W/m²K

Critical and minimum heat flux

ItemValue
Critical heat flux (Zuber)1.26 MW/m²
Superheat at CHF (Rohsenow)20.9 K
Heat transfer coefficient at CHF60,000 W/m²K
Minimum heat flux (Leidenfrost point)18.9 kW/m²

Burnout under heat-flux control

Suppose a 1 mm wire is electrically heated and moves into film boiling while held at the critical 1.26 MW/m². Finding the superheat at which vapour-film conduction (Bromley) plus radiation (emissivity 0.8) carries that flux gives about 1,700 K. Film-side coefficients are roughly 360 W/m²K convective and 470 radiative — under a hundredth of nucleate boiling's 60,000.

🙋

From 21 K to 1,700 K — just from a slight increase in heat flux?

🎓

With imposed heat flux, only two points on the boiling curve carry that flux: one in nucleate and one in film boiling. Exceed CHF even slightly and there's no nucleate point left, so the system jumps to the film-boiling point at the same flux. The vapour film acts like insulation, so the temperature difference grows by orders of magnitude. Copper melts at 1,085 °C, so a copper wire would melt. That's why reactor fuel rods and high-power electronics keep margin to CHF. With imposed surface temperature (e.g. steam heating), transition boiling can be followed and no jump occurs.

Analysis and design tips

  1. First draw the boiling curve from correlations and locate the operating point.
  2. For heat-flux-driven devices (heaters, fuel rods, chips), always evaluate margin to CHF.
  3. Csf changes with finish and fouling; calibrate with measurements if available.
  4. With flow (forced-convection boiling), use flow-boiling correlations, not pool boiling ones.
  5. When resolving bubbles in CFD, verify interface tracking and wall nucleation models.
Coffee Break Trivia

“The vapour film in quenching”

When hot steel is quenched in water, the surface is at first wrapped in a vapour film and cools surprisingly slowly (film boiling). As temperature falls past the Leidenfrost point, the film collapses into nucleate boiling and cooling suddenly accelerates. If this switch happens at different times across the part, cooling becomes uneven, causing distortion and quench cracks. Agitating the water or using additives that destabilise the film are ways to shorten and even out film boiling.

Common mistakes

Mistakes and fixes

MistakeEffectFix
Constant heat transfer coefficientLarge error with superheatVary along the boiling curve
Ignoring CHF marginBurnoutEvaluate margin
Unchecked CsfHeat flux off by 2×Calibrate
Pool correlation for flow boilingWrong CHFFlow-boiling correlations
Neglecting radiation in film boilingWall temperature overestimatedAdd radiation
🙋

I'd like to learn related topics.

🎓
Related Simulators

Experience the theory firsthand with the interactive simulator for this field

All Simulators

Related fields

Structural AnalysisElectromagnetic AnalysisThermal Analysis
Rate this article
Thank you for your feedback!
Helpful
More details
Report error
Helpful
0
More details
0
Report error
0