Boiling Curve Simulator Back
Heat Transfer Simulator

Boiling Curve Simulator — Nukiyama Curve and Pool Boiling

As you raise or lower the heating, the operating point traverses the Nukiyama curve in real time. Under heat-flux control, crossing CHF jumps the point to film boiling (burnout); under temperature control it follows the full S-curve — a clear view of the dangerous hysteresis.

Parameters
Control mode

Heat-flux control = burnout when q'' exceeds CHF (jump to film boiling). Temperature control = traverse the full S-curve including the transition regime.

Heating target (input level)
%

During the animation the operating point moves along the curve toward this target. 0% = near room temperature, 100% = maximum input.

Nucleate boiling coefficient C
Critical heat flux CHF
MW/m²
Emissivity ε (film boiling)
Presets

Water at 1 atm (T_sat = 100°C) is fixed. Natural convection h_nc ≈ 1000 W/(m²·K) and film boiling h_fb ≈ 30 W/(m²·K) are assumed. The Zuber CHF is about 1.1 MW/m² for water at 1 atm.

Results (live)
Boiling regime
Wall superheat ΔT_e
Heat flux q''
Surface temp T_w
Critical heat flux CHF
Leidenfrost point q_min
Margin to CHF
Curve traversal and surface state

Top = boiling curve (log q'' vs ΔT_e). Yellow dot = current operating point / red dashed = CHF and Leidenfrost / orange arrow = burnout jump. Bottom = heating surface (bubbles → vapor film).

Theory & Key Formulas

Pool boiling heat flux $q''$ exhibits four distinct regimes as a function of wall superheat $\Delta T_e = T_w - T_\text{sat}$, as first mapped by Nukiyama (1934).

Natural convection ($\Delta T_e < 5$ K). $h_{nc}$ is the natural convection coefficient:

$$q'' = h_{nc}\,\Delta T_e$$

Nucleate boiling ($5 \leq \Delta T_e < 30$ K). Simplified Rohsenow form:

$$q'' = C\,\Delta T_e^{3}$$

Critical heat flux (Zuber, ≈ 1.1 MW/m² for water at 1 atm):

$$q''_\text{CHF} \approx 0.131\,\rho_v^{1/2}\,h_{fg}\,[\sigma g(\rho_l-\rho_v)]^{1/4}$$

Film boiling (Bromley convection plus radiation):

$$q'' = h_{fb}\,\Delta T_e + \varepsilon\sigma_{SB}(T_w^4 - T_\text{sat}^4)/2$$

Under heat-flux control, once $q''_\text{input} > q''_\text{CHF}$ the surface cannot remain in nucleate boiling and jumps discontinuously to the film branch, so the surface temperature $T_w$ soars by hundreds to a thousand degrees (burnout = hysteresis).

What is the Boiling Curve Simulator

🙋
I always thought boiling was just water making bubbles, but this curve has four different regions. What makes them so different from each other?
🎓
Great catch. Back in 1934, Professor Nukiyama in Japan showed experimentally that boiling actually has four distinct faces. Roughly: as you slowly raise the surface temperature, the liquid first moves by quiet natural convection, then begins vigorous bubbling (nucleate boiling), then a strange unstable regime where heat transfer gets worse (transition), and finally film boiling where vapor blankets the surface. Press "Play" and the yellow operating point actually travels along the curve, while the lower panel shows bubbles turning into a vapor film. The plot uses log axes because the heat flux varies by orders of magnitude.
🙋
Wait — the nucleate region shoots up steeply, then suddenly the curve goes back down. What is actually happening there physically?
🎓
In nucleate boiling, vapor bubbles nucleate from tiny surface cavities and vigorously stir the liquid, so heat transfer ramps up as roughly $q\propto\Delta T^3$. But once superheat exceeds about 30 K, bubbles merge into vapor patches that cling to the wall. That is transition boiling. Vapor is a thermal insulator, so even as wall temperature rises, the net heat flux falls. Run the "Temperature-controlled full curve" preset and you will see the operating point crest the CHF peak and then slide back down the descending branch.
🙋
There is a red line labeled CHF. What happens if I cross it?
🎓
That is the dreaded burnout. When you control the heat flux — like an electric heater — and push above CHF, the system cannot stay in nucleate boiling. Following the orange arrow, it jumps instantly to the film-boiling branch. That is the hysteresis. Run the "Ramp up to burnout" preset: the moment the point reaches the nucleate peak, the surface temperature (the T_w card) soars by hundreds of degrees. That is what melts boiler tubes or fuel cladding. It is why we design with CHF margins of 1.5 to 2.
🙋
So heat-flux control and temperature control behave differently? Does temperature control avoid the jump?
🎓
Exactly — that is the crux. If you can control the wall temperature directly (steam reheaters, solar receivers), raising superheat lets you traverse the transition regime as a smooth descending branch, so no jump occurs and you ease into film boiling. Under heat-flux control the transition regime is unstable and cannot exist — so the point has no choice but to leap to the film branch at CHF. Same curve, but what you hold fixed decides whether you burn out. On the right, the Leidenfrost point is the gateway to film boiling — the same state where droplets dance on a hot pan. The insulating vapor film makes heat transfer plummet.

Frequently Asked Questions

When heat flux exceeds the critical heat flux (CHF), the heating surface becomes blanketed by a vapor film, heat transfer collapses, and surface temperature jumps by hundreds of degrees, often destroying the wall. It is the dominant failure mode for heat-flux-controlled systems like electric heaters, boiler tubes, and electronic chip cooling. Designers typically maintain a safety margin of 1.5 to 2 against CHF.
In BWRs and in PWR accident cooling (ECCS), avoiding the Departure from Nucleate Boiling (DNB) on fuel rod surfaces is the central safety concern. Once film boiling sets in, the heat transfer coefficient drops by orders of magnitude and cladding temperature soars, risking fuel damage. Safety analysis requires the DNB Ratio (DNBR) to stay above about 1.3 under all transients.
When a droplet contacts a very hot surface, the bottom vaporizes instantly and the droplet rides on its own vapor layer. This corresponds to the Leidenfrost point (q_min) on the boiling curve; beyond it lies the stable film boiling regime. The familiar dancing of water on a hot skillet is the same phenomenon. The insulating vapor film causes a dramatic drop in heat transfer.
Pool boiling occurs in stagnant liquid heated only by the surface; it is described by the simple Nukiyama curve. Forced convection boiling occurs in flowing liquid and depends additionally on flow velocity, subcooling and channel geometry, making it far more complex. Industrial steam generators and reactor fuel bundles operate in forced convection boiling, but understanding the four pool boiling regimes is the essential starting point.

Real-world Applications

Thermal and nuclear power boilers: Power stations boil water at high pressure to produce steam, and a tube that crosses CHF will rupture instantly. Designers evaluate CHF with correlations such as the Groeneveld look-up tables and maintain safety margins of at least 1.3 to 2.0. Understanding the boiling curve is the foundation for safe power-plant operation.

High-density electronics cooling: CPUs, GPUs and power semiconductors already dissipate over 100 W/cm², beyond what air cooling can handle. Immersion cooling and jet-impingement cooling exploit the high heat transfer coefficients of nucleate boiling (CHF reaches the MW/m² range). Data centers and EV inverters are early adopters of this technology.

Steel quenching and heat treatment: When hot steel is quenched into water or oil, the surface initially enters film boiling and cools slowly. As temperature drops past the Leidenfrost point, the regime shifts back to nucleate boiling and rapid cooling begins. The dramatic change in cooling rate controls phase transformations, so quench-oil selection and agitation are essentially design choices about which boiling regime to traverse.

Cryogenics and aerospace: Superconducting magnets and rocket engines using liquid helium or liquid nitrogen face boiling at the cold-fluid side of warm walls. Such systems often start in film boiling, which is inefficient. Designers use porous or treated surfaces to promote an early transition to nucleate boiling and improve cool-down rates.

Common Misconceptions and Pitfalls

The most common misconception is to assume that hotter surface always means better heat transfer. In reality, beyond the CHF peak of nucleate boiling, raising wall temperature actually reduces the heat flux in the transition regime. Run the "Temperature-controlled full curve" preset: once the point passes the peak, the heat flux card decreases instead of increasing. This counter-intuitive behavior is the defining feature of boiling heat transfer — and the root cause of many design accidents.

The second common error is treating CHF as an absolute physical limit. CHF is the maximum heat flux at which nucleate boiling can be maintained; crossing it does not always destroy the surface. In heat-flux-controlled systems (electric heaters) you jump immediately to film boiling and burnout occurs. In wall-temperature-controlled systems (reheaters, solar receivers) you traverse the transition regime smoothly and may avoid catastrophic failure. Switch the control mode in this simulator and you will see that the same curve either burns out or does not. Always confirm whether your real system is flux-controlled or temperature-controlled before applying the results.

Finally, remember that this simulator uses a simplified model for water at 1 atm. The coefficients C, CHF and h_fb are representative textbook values. Accurate engineering predictions require evaluating Rohsenow's surface coefficient $C_{sf}$, Zuber's CHF correlation, forced-convection corrections, and channel-geometry effects individually. CHF varies sharply with pressure (rising at high pressure, falling at vacuum). Real plant designs always rely on validated correlations such as the Groeneveld look-up table or subchannel codes (VIPRE, COBRA).

How to Use

  1. Choose the control mode (heat-flux or temperature). Heat-flux control burns out when CHF is crossed; temperature control traverses the full S-curve continuously
  2. Press "Play" and the operating point moves along the boiling curve toward the heating target, while the lower panel shows the bubble-to-vapor-film state simultaneously
  3. Adjust the critical heat flux CHF (MW/m²) for your material and pressure; about 1.1 MW/m² is representative for water at 1 atm (Zuber correlation)
  4. Vary the nucleate coefficient and emissivity to reshape the curve, and watch the live numbers (ΔT_e, q'', surface temperature, margin to CHF)

Worked Example

For water at 1 atm with CHF = 1.1 MW/m² and nucleate coefficient C = 100, nucleate boiling reaches its peak (CHF) at a superheat of ΔT_e ≈ 22 K. Under heat-flux control, running "Ramp up to burnout" makes the operating point jump to the film branch the instant it hits ΔT_e ≈ 22 K, and the surface temperature T_w soars from around 100 °C to the 1000 °C range. Under temperature control over the same curve, the point crests the ΔT_e = 22 K peak and continues smoothly down the transition regime: the surface gets hotter yet the heat flux q'' falls from about 1.1 MW/m². This non-monotonicity, and whether hysteresis (the jump) occurs depending on the control mode, is the heart of the boiling curve.

Practical Notes

  1. Csf values: 0.013 for stainless steel, 0.011 for copper, 0.006 for polished surfaces—use surface-specific data to avoid non-physical predictions
  2. CHF is fluid- and pressure-dependent; water at 1 atm reaches 1200–1300 kW/m², while cryogenic fluids and high-pressure coolants differ significantly
  3. Film boiling (ΔT > 120 K for water) exhibits unstable heat transfer; maintain margin >15% to CHF in design to prevent cycling and burnout
  4. Radiative contribution becomes dominant only above 200 K superheat in atmospheric pool boiling; verify emissivity input for oxidized or coated surfaces

🎬 Watch it in motion

Phase Transitions of Matter Explained | Melting, Boiling and Sublimation Visualized
Phase Transitions of Matter Explained | Melting, Boiling and Sublimation Visualized