Plot the Nukiyama boiling curve in real time. Automatically identify free convection, nucleate boiling, transition, and film boiling regimes. Compute CHF and Leidenfrost point with the Rohsenow & Zuber correlations.
Film boiling (Bromley): $h_{film}=0.62\left[\dfrac{k_v^3\rho_v(\rho_l-\rho_v)g h_{fg}}{\mu_v D\Delta T_e}\right]^{1/4}$
What is Boiling Heat Transfer?
🙋
What exactly is a "boiling curve," and why is it so important for engineers?
🎓
Basically, it's a graph that shows how the heat flux (the energy flow per area) changes as you make a surface hotter than the liquid's boiling point. It's crucial because it reveals four distinct boiling regimes. For instance, in electronics cooling, you want to operate in the efficient nucleate boiling zone, but accidentally cross into film boiling and your component can overheat catastrophically. Try moving the "Surface Excess Temp. ΔT" slider in the simulator above to see the curve and these regimes appear instantly.
🙋
Wait, really? So the curve isn't just a straight line? What's happening in that "nucleate boiling" peak you mentioned?
🎓
Right, it's highly nonlinear! In nucleate boiling, bubbles form vigorously at nucleation sites on the surface, which stirs the fluid and transfers a huge amount of heat. The peak of that curve is called the Critical Heat Flux (CHF) – it's the maximum cooling power you can get. In practice, this is the limit for a safe nuclear reactor core or a high-performance computer chip. The simulator uses the Rohsenow correlation to calculate this part of the curve. Try changing the "Fluid" from water to a refrigerant like R-134a and see how much lower the peak heat flux becomes.
🙋
That makes sense. So what's the deal with the "Surface–Fluid Pair" parameter? It just says "C_sf"...
🎓
Great question! $C_{sf}$ is an empirical constant that captures how easily bubbles form on a specific material-fluid combination. A polished copper surface in water has a different $C_{sf}$ than a stainless steel surface in refrigerant. It directly affects the temperature needed to start boiling. A common case is selecting a coated surface in an immersion cooling system to optimize bubble formation. When you change the "Surface–Fluid Pair" dropdown in the tool, you're changing this constant and seeing how the entire nucleate boiling curve shifts.
Physical Model & Key Equations
The Rohsenow correlation models the heat flux during the nucleate boiling regime, where bubble dynamics dominate heat transfer. It relates the heat flux to the excess temperature ($\Delta T_e = T_{surface}- T_{sat}$).
Where:
$q''$ = Heat flux (W/m²)
$\mu_l$ = Liquid dynamic viscosity (Pa·s)
$h_{fg}$ = Latent heat of vaporization (J/kg)
$g$ = Gravitational acceleration (m/s²)
$\rho_l, \rho_v$ = Liquid and vapor densities (kg/m³)
$\sigma$ = Surface tension (N/m)
$c_{pl}$ = Liquid specific heat (J/kg·K)
$\Delta T_e$ = Excess temperature (K)
$C_{sf}$ = Surface-fluid coefficient (from experiments)
$\mathrm{Pr}$ = Liquid Prandtl number
$n$ = Experimental exponent (often 1 for water, 1.7 for other fluids)
The Critical Heat Flux (CHF) marks the peak of the boiling curve and the transition to dangerous film boiling. A widely used correlation is the Zuber relation, which models the hydrodynamic instability limiting bubble departure.
$$q''_{max}= 0.149 h_{fg}\rho_v \left[\frac{\sigma g (\rho_l-\rho_v)}{\rho_v^2}\right]^{1/4}$$
This equation shows that CHF depends heavily on fluid properties like latent heat ($h_{fg}$), density difference ($\rho_l-\rho_v$), and surface tension ($\sigma$). The physical meaning is that at a certain vapor generation rate, bubbles coalesce into a continuous insulating vapor film, causing the heat flux to drop sharply even as temperature rises—a condition known as "burnout."
Frequently Asked Questions
They are determined based on the value of superheat (ΔT_e). Nucleate boiling is defined as low superheat, transition boiling as the range between the critical heat flux (CHF) and the Leidenfrost point, and film boiling as superheat beyond that. The regions are color-coded on the graph for visual confirmation.
C_sf and n are experimental constants that depend on the combination of heating surface material and fluid. Default values (e.g., C_sf=0.013, n=1.0 for water-copper) are set, but they can be manually changed by referring to literature values. For improved accuracy, fitting with actual measurement data is recommended.
CHF is calculated using Zuber's equation (q''_CHF=0.131ρ_g^{1/2}h_fg[gσ(ρ_l-ρ_g)]^{1/4}), and the Leidenfrost point is calculated using an empirical formula based on the saturation temperature. The results are displayed as markers on the graph and can be referenced as superheat thresholds.
This tool uses representative values at the saturation temperature. For systems with strong temperature dependence (e.g., high superheat), errors may occur, so please input values near the saturation temperature as needed. For higher-precision analysis, consider incorporating a temperature-dependent model separately.
Real-World Applications
Nuclear Reactor Safety (Burnout Analysis): The most critical application is preventing the "departure from nucleate boiling" (DNB) in reactor fuel rods. Engineers use CHF calculations to set absolute power limits, ensuring the cladding temperature never enters the film boiling regime, which could lead to meltdown.
High-Power Electronics Immersion Cooling: Servers for AI and cryptocurrency mining are submerged in dielectric fluids. Designers use boiling curves to select fluids and surface finishes that maximize heat flux in the nucleate regime, allowing compact, silent cooling of chips exceeding 1000 W/cm².
Industrial Boiler & Evaporator Design: In steam generation systems, the boiling curve dictates the required heating surface area and tube material. Operating near the CHF maximizes efficiency but requires precise control of pressure and temperature to avoid scaling or tube failure.
Metal Quenching (Heat Treatment): When hot steel is plunged into oil or water, it passes through all boiling regimes rapidly. Controlling the cooling rate by manipulating the boiling process (e.g., via agitation) is essential to achieve the desired material hardness and prevent cracking.
Common Misconceptions and Points to Note
When starting to use this tool, there are several common misconceptions, especially among those learning CAE on the job. First and foremost, the Rohsenow correlation is not a universal solution. This equation provides a guideline for the nucleate boiling regime, and relying solely on it for actual equipment design is risky. For instance, if the flow passage is narrow or the heating surface is inclined, significant discrepancies can arise between calculated and measured values. Second, choose the experimental constant Csf with care. While the tool allows you to select typical combinations, actual surface roughness and fouling (scaling) vary immensely. It's not uncommon for a heat exchanger designed using a value for polished copper (Csf~0.013) to underperform predictions due to slight changes in surface condition during manufacturing. The third pitfall is avoiding setting the CHF as a target value without a safety margin. Operating right at the peak of the boiling curve can lead to a transition to film boiling with just minor condition fluctuations (e.g., system pressure changes), causing a rapid temperature rise (burnout). In practice, it's standard to incorporate a safety factor of at least 1.5 to 2.0 for the CHF (actual heat flux ≤ CHF / safety factor).