9.3 Boiling Regimes, Condensation Phenomena, and Fouling Resistances

Key Takeaways

  • The overall thermal resistance network sums five series resistances: 1/(U_o*A_o) = 1/(h_i*A_i) + R_f,i/A_i + ln(r_o/r_i)/(2*pi*k*L) + R_f,o/A_o + 1/(h_o*A_o); fouling factors degrade the design overall coefficient via 1/U_dirty = 1/U_clean + R_f,total.
  • The classical pool boiling curve comprises four distinct regimes characterized by excess temperature Delta_T_e = T_s - T_sat: natural convection (Delta_T_e < 5°C), nucleate boiling (5°C < Delta_T_e < 30°C, yielding peak heat transfer), transition boiling (unstable, negative slope dq''/dDelta_T_e < 0), and stable film boiling (Delta_T_e > 120°C).
  • Critical Heat Flux (CHF / Departure from Nucleate Boiling DNB) defines the upper limit of nucleate boiling, predicted by Zuber's equation: q''_max = 0.149 * rho_v * h_fg * [sigma * g * (rho_l - rho_v) / rho_v²]^(1/4); exceeding CHF in a heat-flux-controlled system triggers catastrophic wall temperature runaway ('burnout') into the film boiling regime.
  • Nusselt film condensation theory demonstrates that horizontal tubes achieve significantly higher heat transfer coefficients than vertical tubes of identical length: h_horiz / h_vert = 0.773 * (L / D_o)^(1/4), making horizontal condenser bundles vastly superior when L / D_o >> 1.
  • Dropwise condensation produces heat transfer coefficients 5 to 10 times higher than filmwise condensation (30,000-100,000 W/(m²*K) vs 2,000-10,000 W/(m²*K)) because droplets continually coalesce and roll off, exposing bare metal; however, industrial designs conservatively assume filmwise condensation due to the difficulty of permanently maintaining non-wetting surface promoters.
Last updated: September 2026

9.3 Boiling Regimes, Condensation Phenomena, and Fouling Resistances

Phase-change operations—specifically boiling and condensation—exhibit heat transfer coefficients ($h$) that are orders of magnitude larger than single-phase liquid or gas convection. However, these operations are constrained by complex hydrodynamic boundaries, such as the Critical Heat Flux (CHF) in reboilers and condensate film thickening in condensers. Furthermore, real process equipment inevitably accumulates insulating deposits over time, requiring rigorous application of fouling resistances ($R_f$) in overall heat transfer calculations.


1. Overall Heat Transfer Resistance & Fouling Factors

Heat transfer across a cylindrical tube wall involves five distinct thermal resistances acting in series:

  1. Convective resistance of the tube-side boundary layer.
  2. Conductive resistance of the internal fouling deposit (scale, biological growth, coking).
  3. Conductive resistance of the metallic tube wall.
  4. Conductive resistance of the external shell-side fouling layer.
  5. Convective resistance of the shell-side boundary layer.
  Tube Fluid       Inner Scale        Tube Wall       Outer Scale      Shell Fluid
   (T_c)            (R_f,i)             (k_w)           (R_f,o)          (T_h)
     |                 |                  |                |               |
  ==[1/(h_i*A_i)]=====[R_f,i/A_i]=====[ln(r_o/r_i)]======[R_f,o/A_o]=====[1/(h_o*A_o)]==
                                       (2*pi*k_w*L)

Formulations Based on Outside Surface Area ($A_o$)

The total thermal resistance is expressed as:

1UoAo=1hiAi+Rf,iAi+ln(ro/ri)2πkwL+Rf,oAo+1hoAo\frac{1}{U_o A_o} = \frac{1}{h_i A_i} + \frac{R_{f,i}}{A_i} + \frac{\ln(r_o / r_i)}{2 \pi k_w L} + \frac{R_{f,o}}{A_o} + \frac{1}{h_o A_o}

Multiplying through by outside area $A_o$ ($A_o / A_i = r_o / r_i = d_o / d_i$) yields the design equation for the outside overall coefficient ($U_o$):

1Uo=1ho+Rf,o+roln(ro/ri)kw+Rf,i(rori)+1hi(rori)\frac{1}{U_o} = \frac{1}{h_o} + R_{f,o} + \frac{r_o \ln(r_o / r_i)}{k_w} + R_{f,i} \left( \frac{r_o}{r_i} \right) + \frac{1}{h_i} \left( \frac{r_o}{r_i} \right)

Where:

  • $U_o$ = overall heat transfer coefficient referenced to outside tube area ($\text{W/(m}^2\cdot\text{K)}$).
  • $h_i, h_o$ = inner and outer convective heat transfer coefficients ($\text{W/(m}^2\cdot\text{K)}$).
  • $R_{f,i}, R_{f,o}$ = inner and outer fouling resistances ($\text{m}^2\cdot\text{K/W}$ or $\text{hr}\cdot\text{ft}^2\cdot^\circ\text{F/Btu}$).
  • $r_i, r_o$ = inner and outer tube radii ($d_i/2, d_o/2$).
  • $k_w$ = thermal conductivity of the tube wall material ($\text{W/(m}\cdot\text{K)}$).
  • $L$ = total tube length ($\text{m}$).

Clean vs. Dirty Overall Heat Transfer Coefficients

In process design, the clean overall heat transfer coefficient ($U_{clean}$) represents newly fabricated, un-fouled equipment ($R_f = 0$):

1Uclean=1ho+roln(ro/ri)kw+1hi(rori)\frac{1}{U_{clean}} = \frac{1}{h_o} + \frac{r_o \ln(r_o / r_i)}{k_w} + \frac{1}{h_i} \left( \frac{r_o}{r_i} \right)

The dirty (design) overall coefficient ($U_{dirty}$) incorporates the total combined fouling resistance ($R_{f,total}$):

1Udirty=1Uclean+Rf,total    Udirty=Uclean1+UcleanRf,total\frac{1}{U_{dirty}} = \frac{1}{U_{clean}} + R_{f,total} \implies U_{dirty} = \frac{U_{clean}}{1 + U_{clean} \cdot R_{f,total}}

Where $R_{f,total} = R_{f,o} + R_{f,i}(r_o / r_i)$.

Required Overdesign Surface Margin

Because fouling reduces the heat transfer coefficient from $U_{clean}$ to $U_{dirty}$, the required surface area must be expanded to guarantee performance throughout the operating cycle between turnarounds:

AdesignAclean=UcleanUdirty=1+UcleanRf,total\frac{A_{design}}{A_{clean}} = \frac{U_{clean}}{U_{dirty}} = 1 + U_{clean} \cdot R_{f,total}

Overdesign %=(AdesignAcleanAclean)×100%=UcleanRf,total×100%\text{Overdesign \%} = \left( \frac{A_{design} - A_{clean}}{A_{clean}} \right) \times 100\% = U_{clean} \cdot R_{f,total} \times 100\%

Standard TEMA Fouling Resistances

Process Fluid ServiceMetric Fouling Resistance $R_f$ ($\text{m}^2\cdot\text{K/W}$)US Customary $R_f$ ($\text{hr}\cdot\text{ft}^2\cdot^\circ\text{F/Btu}$)
Treated Boiler Feedwater ($<50^\circ\text{C}$)$0.00009 - 0.00018$$0.0005 - 0.0010$
Treated Cooling Tower Water ($<50^\circ\text{C}$)$0.00018 - 0.00035$$0.0010 - 0.0020$
River Water / Untreated Surface Water$0.00035 - 0.00053$$0.0020 - 0.0030$
Clean Light Hydrocarbons (Naphtha, LPG)$0.00018$$0.0010$
Heavy Gas Oil / Fuel Oil$0.00053 - 0.00088$$0.0030 - 0.0050$
Crude Oil Bottoms / Vacuum Resid$0.00088 - 0.00176$$0.0050 - 0.0100$
Dry Gas (Air, Nitrogen, Natural Gas)$0.00018 - 0.00035$$0.0010 - 0.0020$

2. The Pool Boiling Curve and Regimes

Boiling occurs when a solid heating surface is maintained at a temperature higher than the saturation temperature of the surrounding liquid. The thermodynamic driving force is the excess temperature ($\Delta T_e$):

ΔTeTsTsat\Delta T_e \equiv T_s - T_{sat}

Plotting heat flux ($q'' = Q/A$) versus $\Delta T_e$ on logarithmic coordinates reveals the classical Nukiyama Pool Boiling Curve, which spans four distinct heat transfer regimes:

  log(q'') ^                      Point C: Critical Heat Flux (CHF)
           |                          /\                      Regime IV:
           |                         /  \    Regime III:     Film Boiling
           |                        /    \   Transition     /------------
           |         Regime II:    /      \    Boiling     /
           |          Nucleate    /        \              /  Point D:
           |           Boiling   /          \------------/  Leidenfrost Point
           |                    / 
           |   Regime I:       /
           |    Natural       /
           |   Convection    /
           +----------------------------------------------------> log(Delta_T_e)
           0                 5°C       30°C         120°C

Regime I: Natural / Free Convection ($\Delta T_e \lesssim 5^\circ\text{C}$)

Liquid adjacent to the heated surface warms up and rises due to buoyancy. No vapor bubbles nucleate because the liquid superheat is insufficient to overcome the surface tension forces required to form a bubble nucleus. Heat transfer is governed purely by single-phase natural convection correlations ($q'' \propto \Delta T_e^{5/4}$).

Regime II: Nucleate Boiling ($5^\circ\text{C} \lesssim \Delta T_e \lesssim 30^\circ\text{C}$)

Vapor bubbles nucleate at microscopic pits and cavities on the heated surface. This regime contains two sub-regions:

  1. Isolated Bubbles ($5^\circ\text{C} < \Delta T_e < 10^\circ\text{C}$): Individual bubbles form, expand, and detach, inducing intense micro-convective mixing in the adjacent thermal boundary layer.
  2. Vapor Columns and Slugs ($10^\circ\text{C} < \Delta T_e < 30^\circ\text{C}$): Bubble nucleation frequency becomes so rapid that adjacent bubbles coalesce into continuous vapor jets. Heat flux increases rapidly with excess temperature: $q'' \propto \Delta T_e^3$.

The heat flux in nucleate boiling is predicted by the Rohsenow Correlation:

q=μlhfg[g(ρlρv)σ]1/2[Cp,l(TsTsat)CsfhfgPrln]3q'' = \mu_l h_{fg} \left[ \frac{g (\rho_l - \rho_v)}{\sigma} \right]^{1/2} \left[ \frac{C_{p,l} (T_s - T_{sat})}{C_{sf} h_{fg} Pr_l^n} \right]^3

Where $C_{sf}$ is an empirical surface-fluid constant reflecting surface wettability (e.g., $0.013$ for water-copper; $0.006$ for water-brass), and $n = 1.0$ for water ($1.7$ for organic liquids).

The Critical Heat Flux (CHF / Burnout Point)

The peak of the nucleate boiling curve is the Critical Heat Flux ($q''_{max}$), also known as the Departure from Nucleate Boiling (DNB) or the Burnout Point. At this point, vapor generation is so intense that escaping vapor jets prevent liquid from returning to re-wet the heating surface.

Zuber developed the hydrodynamic stability equation for pool boiling on horizontal surfaces:

qmax=0.149ρvhfg[σg(ρlρv)ρv2]1/4q''_{max} = 0.149 \cdot \rho_v h_{fg} \left[ \frac{\sigma g (\rho_l - \rho_v)}{\rho_v^2} \right]^{1/4}

For water at atmospheric pressure, $q''_{max} \approx 1.1\text{ to }1.3\text{ MW/m}^2$ ($350,000\text{ to }400,000\text{ Btu/(hr}\cdot\text{ft}^2)$).

Regime III: Transition Boiling ($30^\circ\text{C} \lesssim \Delta T_e \lesssim 120^\circ\text{C}$)

Also known as unstable film boiling or partial boiling. A continuous vapor film intermittently forms and collapses over portions of the heating surface. Because vapor has a thermal conductivity an order of magnitude lower than liquid ($k_v \ll k_l$), the vapor patches act as an insulating blanket.

Crucially, the transition boiling curve has a negative slope ($dq'' / d\Delta T_e < 0$). Increasing the surface temperature expands the vapor blanket faster than convective transport can rise, causing net heat flux to paradoxically decrease!

[!CAUTION] Temperature-Controlled vs. Heat-Flux-Controlled Boiling:

  • In temperature-controlled systems (such as reboilers heated by condensing steam), the wall temperature is fixed by steam saturation pressure. The system can safely traverse the transition regime and find a stable operating point.
  • In heat-flux-controlled systems (such as fired heaters, combustion furnaces, and electric immersion heaters), heat flux $q''$ is fixed by fuel firing rate. If heat flux is increased beyond $q''_{max}$ (CHF), the operating point cannot traverse the negative slope. Instead, it undergoes a catastrophic horizontal jump across to the film boiling curve (Point D), where the wall temperature spikes by several hundred degrees ($> 1,000^\circ\text{C}$). This thermal runaway melts or ruptures the heater tubes—a failure mode termed burnout.

Regime IV: Film Boiling & The Leidenfrost Point ($\Delta T_e > 120^\circ\text{C}$)

The minimum heat flux ($q''_{min}$) in film boiling is called the Leidenfrost Point:

qmin=0.09ρvhfg[σg(ρlρv)(ρl+ρv)2]1/4q''_{min} = 0.09 \cdot \rho_v h_{fg} \left[ \frac{\sigma g (\rho_l - \rho_v)}{(\rho_l + \rho_v)^2} \right]^{1/4}

Beyond the Leidenfrost point, the heating surface is completely covered by a stable, continuous vapor blanket. Heat must be transferred across this insulating blanket by conduction and radiation:

h=hconv(hconvh)1/3+hradhconv+34hradh = h_{conv} \left( \frac{h_{conv}}{h} \right)^{1/3} + h_{rad} \approx h_{conv} + \frac{3}{4} h_{rad}

Where the radiation coefficient is $h_{rad} = \epsilon_{em} \sigma_{SB} (T_s^4 - T_{sat}^4) / (T_s - T_{sat})$.


3. Condensation Phenomena & Nusselt Film Theory

Condensation occurs when a saturated or superheated vapor contacts a solid surface maintained below its saturation temperature ($T_s < T_{sat}$).

Filmwise vs. Dropwise Condensation

  1. Filmwise Condensation: The liquid condensate completely wets the surface, forming a continuous liquid film that flows downward under gravity. This liquid film acts as a conductive thermal barrier between the vapor and the cold wall. Filmwise heat transfer coefficients typically range from $2,000$ to $10,000\text{ W/(m}^2\cdot\text{K)}$.
  2. Dropwise Condensation: The liquid condensate does not wet the surface (contact angle $\theta > 90^\circ$). Instead, microscopic droplets nucleate at surface cavities, grow by coalescence, and roll down the plate, continually exposing bare metal surface directly to vapor. Dropwise condensation achieves heat transfer coefficients of $30,000$ to $100,000\text{ W/(m}^2\cdot\text{K)}$—a $5$- to $10$-fold enhancement over filmwise condensation!

[!NOTE] Industrial Condenser Design Practice:
Although dropwise condensation provides extraordinary heat transfer rates, it requires hydrophobic chemical promoter coatings (e.g., fluorocarbons, noble metal plating, oleic acid) that degrade and wash away after weeks of process exposure. Because sustaining dropwise condensation in chemical plants is unreliable, all industrial heat exchangers are conservatively sized assuming filmwise condensation.

Nusselt Laminar Film Condensation Theory (1916)

Wilhelm Nusselt formulated the classical analytical model for laminar film condensation on a vertical flat plate or tube based on six core physical assumptions:

  1. The condensate film flows under pure laminar conditions ($Re_f < 30$).
  2. Vapor is pure, saturated, and stationary (zero vapor shear stress at the interface).
  3. Heat transfer across the condensate film occurs purely by steady-state 1D conduction.
  4. Subcooling of the liquid film is negligible ($T$ decreases linearly across the film from $T_{sat}$ at the liquid-vapor interface to $T_s$ at the solid wall).
  5. Liquid acceleration and convective momentum terms in the Navier-Stokes equations are negligible.
  6. Thermophysical properties are evaluated at the film reference temperature: $T_f = (T_{sat} + T_s) / 2$.

From the balance of gravity and viscous shear forces ($d\tau/dy = -\rho_l g$), the condensate film thickness $\delta(x)$ grows along the vertical drainage length $x$ as:

δ(x)=[4klμl(TsatTs)xgρl(ρlρv)hfg]1/4x1/4\delta(x) = \left[ \frac{4 k_l \mu_l (T_{sat} - T_s) x}{g \rho_l (\rho_l - \rho_v) h_{fg}'} \right]^{1/4} \propto x^{1/4}

Where $h_{fg}'$ is the modified latent heat of vaporization, corrected for liquid film subcooling (Rohsenow correction):

hfg=hfg+0.68Cp,l(TsatTs)h_{fg}' = h_{fg} + 0.68 \cdot C_{p,l} (T_{sat} - T_s)

Because local conduction resistance is $\delta(x) / k_l$, the local heat transfer coefficient is $h_x = k_l / \delta(x) \propto x^{-1/4}$. Integrating over a vertical surface of length $L$ yields the average vertical condensation coefficient:

hˉvert=0.943[gρl(ρlρv)kl3hfgμlL(TsatTs)]1/4\bar{h}_{vert} = 0.943 \left[ \frac{g \rho_l (\rho_l - \rho_v) k_l^3 h_{fg}'}{\mu_l L (T_{sat} - T_s)} \right]^{1/4}

Nusselt Condensation on a Horizontal Tube

For condensation on the outer surface of a single horizontal tube of diameter $d_o$, the condensate film drains around the circumference, traveling a characteristic drainage distance of only $\pi d_o / 2$ before dripping off the bottom. Nusselt's derivation yields:

hˉhoriz=0.729[gρl(ρlρv)kl3hfgμldo(TsatTs)]1/4\bar{h}_{horiz} = 0.729 \left[ \frac{g \rho_l (\rho_l - \rho_v) k_l^3 h_{fg}'}{\mu_l d_o (T_{sat} - T_s)} \right]^{1/4}

Vertical vs. Horizontal Tube Comparison

Dividing the horizontal expression by the vertical expression for identical fluid properties, $\Delta T$, and tube dimensions ($d_o, L$):

hˉhorizhˉvert=0.7290.943(Ldo)1/4=0.773(Ldo)1/4\frac{\bar{h}_{horiz}}{\bar{h}_{vert}} = \frac{0.729}{0.943} \left( \frac{L}{d_o} \right)^{1/4} = 0.773 \left( \frac{L}{d_o} \right)^{1/4}

For a standard industrial tube where $L = 4.0\text{ m}$ and $d_o = 0.0254\text{ m}$ ($1.0\text{ in}$): Ldo=4.00.0254=157.48\frac{L}{d_o} = \frac{4.0}{0.0254} = 157.48 (Ldo)1/4=(157.48)0.25=3.541\left( \frac{L}{d_o} \right)^{1/4} = (157.48)^{0.25} = 3.541 hˉhorizhˉvert=0.773×3.541=2.737\frac{\bar{h}_{horiz}}{\bar{h}_{vert}} = 0.773 \times 3.541 = \mathbf{2.737}

[!IMPORTANT] Why Process Condensers Are Arranged Horizontally:
A horizontal tube provides nearly $3$ times higher heat transfer coefficient than the same tube placed vertically! On a vertical tube, condensate drains down the entire length ($4.0\text{ m}$), accumulating into a thick, thermally insulating liquid blanket at the bottom. On a horizontal tube, condensate only travels around the narrow outer perimeter before falling away, keeping the average boundary film extraordinarily thin.

Condensate Inundation in Horizontal Tube Bundles

In a multi-row horizontal condenser, condensate dripping from upper tubes lands on lower tubes, thickening their liquid films. For a vertical column of $N$ tubes, Nusselt derived the average bundle coefficient as:

hˉbundle=hˉsingleN1/4\bar{h}_{bundle} = \bar{h}_{single} \cdot N^{-1/4}

Kern proposed a less conservative empirical exponent to account for turbulence and splashing caused by dripping droplets:

hˉbundle,Kern=hˉsingleN1/6\bar{h}_{bundle, Kern} = \bar{h}_{single} \cdot N^{-1/6}


4. Summary Table of Boiling and Condensation Regimes & Equations

Phenomenon / RegimeGoverning MechanismCharacteristic $\Delta T_e$ RangeGoverning Correlation / Analytical ModelKey Engineering Design Rule
Natural Convection BoilingSingle-phase buoyant liquid circulation$\Delta T_e < 5^\circ\text{C}$$q'' = h_{nc} \Delta T_e \propto \Delta T_e^{5/4}$Superheat insufficient to nucleate bubbles
Nucleate BoilingBubble nucleation, growth, and micro-convection$5^\circ\text{C} < \Delta T_e < 30^\circ\text{C}$Rohsenow: $q'' \propto \Delta T_e^3$Optimal design regime; highest heat flux per $m^2$
Critical Heat Flux (CHF)Vapor columns choke liquid replenishmentPeak at $\Delta T_e \approx 30^\circ\text{C}$Zuber: $q''{max} = 0.149 \rho_v h{fg} [\sigma g (\Delta\rho)/\rho_v^2]^{1/4}$Maximum safe thermal threshold in fired reboilers
Transition BoilingIntermittent, unstable vapor blanket formation$30^\circ\text{C} < \Delta T_e < 120^\circ\text{C}$Negative slope: $dq''/d\Delta T_e < 0$Unstable in flux-controlled heaters (burnout danger)
Film BoilingStable continuous vapor cushion$\Delta T_e > 120^\circ\text{C}$$h = h_{conv} + \frac{3}{4} h_{rad}$; $q''_{min}$ via LeidenfrostExtremely high wall temperatures ($T_s > 800^\circ\text{C}$)
Vertical CondensationLaminar draining liquid filmN/A ($T_s < T_{sat}$)$\bar{h}{vert} = 0.943 [g \rho_l (\Delta\rho) k_l^3 h{fg}' / (\mu_l L \Delta T)]^{1/4}$Film thickens down tube; $h \propto L^{-1/4}$
Horizontal CondensationPeripheral drainage around circumferenceN/A ($T_s < T_{sat}$)$\bar{h}{horiz} = 0.729 [g \rho_l (\Delta\rho) k_l^3 h{fg}' / (\mu_l d_o \Delta T)]^{1/4}$Superior performance: $h_{horiz} / h_{vert} \approx 2.5-3.0$

5. Comprehensive Worked Numerical Example: Surface Condenser Sizing with Fouling

Problem Statement

A chemical processing facility must design a horizontal shell-and-tube surface condenser to condense $15,000\text{ kg/h}$ ($4.167\text{ kg/s}$) of pure saturated organic process vapor at $T_{sat} = 60.0^\circ\text{C}$ on the shell side.

Condensate & Vapor Properties at $60.0^\circ\text{C}$:

  • Latent heat of vaporization: $h_{fg} = 2,358\text{ kJ/kg} = 2.358 \times 10^6\text{ J/kg}$
  • Liquid condensate density: $\rho_l = 985\text{ kg/m}^3$
  • Vapor density: $\rho_v \approx 0.13\text{ kg/m}^3$ (so $\rho_l - \rho_v \approx \rho_l$)
  • Condensate liquid thermal conductivity: $k_l = 0.650\text{ W/(m}\cdot\text{K)}$
  • Condensate liquid viscosity: $\mu_l = 4.80 \times 10^{-4}\text{ Pa}\cdot\text{s}$

Tube Bundle Specifications:

  • Admiralty brass tubes: $k_w = 110.0\text{ W/(m}\cdot\text{K)}$
  • Tube outside diameter: $d_o = 25.4\text{ mm} = 0.0254\text{ m}$
  • Tube wall thickness: $t = 1.65\text{ mm} = 0.00165\text{ m}$
  • Tube inside diameter: $d_i = d_o - 2t = 25.4 - 3.30 = 22.1\text{ mm} = 0.0221\text{ m}$
  • Tube length: $L = 4.50\text{ m}$
  • Average vertical tube count in a vertical column: $N = 9\text{ tubes}$

Operating & Utility Conditions:

  • Cooling water enters tubes at $t_1 = 20.0^\circ\text{C}$ and exits at $t_2 = 35.0^\circ\text{C}$.
  • Tube-side convective heat transfer coefficient: $h_i = 4,800.0\text{ W/(m}^2\cdot\text{K)}$
  • Tube-side cooling water fouling resistance: $R_{f,i} = 0.00025\text{ m}^2\cdot\text{K/W}$
  • Shell-side organic fouling resistance: $R_{f,o} = 0.00010\text{ m}^2\cdot\text{K/W}$
  • Average outer tube wall surface temperature: $T_s \approx 42.0^\circ\text{C}$ (giving $\Delta T = T_{sat} - T_s = 60.0 - 42.0 = 18.0^\circ\text{C}$)

Calculate:

  1. The total condensation heat duty ($Q$) in $\text{kW}$.
  2. The average shell-side condensation coefficient for a single horizontal tube ($\bar{h}_{single}$) using Nusselt theory, and the corrected bundle coefficient ($\bar{h}_o$).
  3. The clean overall heat transfer coefficient ($U_{clean}$) based on outside area.
  4. The design fouled overall heat transfer coefficient ($U_{dirty}$) and total fouling resistance ($R_{f,total}$).
  5. The required outside heat transfer area ($A_o$) using $\Delta T_{lm}$.
  6. The required number of tubes ($N_t$) and the overdesign margin specified for fouling.

Step 1: Condensation Heat Duty

Q=m˙hfg=(4.1667 kg/s)×(2,358 kJ/kg)=9,825.0 kW=9.825×106 WQ = \dot{m} \cdot h_{fg} = (4.1667\text{ kg/s}) \times (2,358\text{ kJ/kg}) = \mathbf{9,825.0}\text{ kW} = \mathbf{9.825 \times 10^6}\text{ W}


Step 2: Nusselt Shell-Side Condensation Coefficient

For a single horizontal tube with $d_o = 0.0254\text{ m}$ and $\Delta T = 18.0\text{ K}$:

hˉsingle=0.729[gρl(ρlρv)kl3hfgμldo(TsatTs)]1/4\bar{h}_{single} = 0.729 \left[ \frac{g \rho_l (\rho_l - \rho_v) k_l^3 h_{fg}}{\mu_l d_o (T_{sat} - T_s)} \right]^{1/4}

Evaluate the bracketed term: Numerator=(9.81 m/s2)×(985 kg/m3)2×(0.650 W/(mK))3×(2.358×106 J/kg)\text{Numerator} = (9.81\text{ m/s}^2) \times (985\text{ kg/m}^3)^2 \times (0.650\text{ W/(m}\cdot\text{K)})^3 \times (2.358 \times 10^6\text{ J/kg}) Numerator=9.81×(970,225)×(0.274625)×(2.358×106)\text{Numerator} = 9.81 \times (970,225) \times (0.274625) \times (2.358 \times 10^6) Numerator=9.81×970,225×647,566=6.1633×1012\text{Numerator} = 9.81 \times 970,225 \times 647,566 = 6.1633 \times 10^{12}

Denominator=μldoΔT=(4.80×104 Pas)×(0.0254 m)×(18.0 K)\text{Denominator} = \mu_l \cdot d_o \cdot \Delta T = (4.80 \times 10^{-4}\text{ Pa}\cdot\text{s}) \times (0.0254\text{ m}) \times (18.0\text{ K}) Denominator=2.1946×104\text{Denominator} = 2.1946 \times 10^{-4}

Ratio=6.1633×10122.1946×104=2.8084×1016\text{Ratio} = \frac{6.1633 \times 10^{12}}{2.1946 \times 10^{-4}} = 2.8084 \times 10^{16}

Taking the fourth root: (2.8084×1016)1/4=12,943.4(2.8084 \times 10^{16})^{1/4} = 12,943.4

Single tube coefficient: hˉsingle=0.729×12,943.4=9,435.7 W/(m2K)\bar{h}_{single} = 0.729 \times 12,943.4 = \mathbf{9,435.7}\text{ W/(m}^2\cdot\text{K)}

Correct for tube bundle inundation across $N = 9$ vertical rows (using Kern's $N^{-1/6}$ correction): hˉo=hˉsingleN1/6=9,435.7×(9)1/6=9,435.7×0.6934=6,542.7 W/(m2K)\bar{h}_o = \bar{h}_{single} \cdot N^{-1/6} = 9,435.7 \times (9)^{-1/6} = 9,435.7 \times 0.6934 = \mathbf{6,542.7}\text{ W/(m}^2\cdot\text{K)}


Step 3: Clean Overall Coefficient ($U_{clean}$)

Radii and area ratios:

  • $r_o = d_o / 2 = 0.0254 / 2 = 0.0127\text{ m}$
  • $r_i = d_i / 2 = 0.0221 / 2 = 0.01105\text{ m}$
  • Ratio $r_o / r_i = 0.0127 / 0.01105 = 1.1493$

Evaluate individual clean thermal resistances referenced to $A_o$:

  1. Shell-side condensation resistance: Ro=1ho=16,542.7=0.0001528 m2K/WR_o = \frac{1}{h_o} = \frac{1}{6,542.7} = \mathbf{0.0001528}\text{ m}^2\cdot\text{K/W}
  2. Tube wall conduction resistance: Rwall=roln(ro/ri)kw=(0.0127)×ln(1.1493)110.0=0.0127×0.13915110.0=0.0000161 m2K/WR_{wall} = \frac{r_o \ln(r_o / r_i)}{k_w} = \frac{(0.0127) \times \ln(1.1493)}{110.0} = \frac{0.0127 \times 0.13915}{110.0} = \mathbf{0.0000161}\text{ m}^2\cdot\text{K/W}
  3. Tube-side convective resistance: Ri=1hi(rori)=14,800.0×1.1493=0.0002394 m2K/WR_i = \frac{1}{h_i} \left( \frac{r_o}{r_i} \right) = \frac{1}{4,800.0} \times 1.1493 = \mathbf{0.0002394}\text{ m}^2\cdot\text{K/W}

Total clean resistance: 1Uclean=Ro+Rwall+Ri=0.0001528+0.0000161+0.0002394=0.0004083 m2K/W\frac{1}{U_{clean}} = R_o + R_{wall} + R_i = 0.0001528 + 0.0000161 + 0.0002394 = \mathbf{0.0004083}\text{ m}^2\cdot\text{K/W}

Uclean=10.0004083=2,449.2 W/(m2K)U_{clean} = \frac{1}{0.0004083} = \mathbf{2,449.2}\text{ W/(m}^2\cdot\text{K)}


Step 4: Dirty Overall Coefficient ($U_{dirty}$)

Fouling resistances referenced to outside area:

  • Shell-side fouling: $R_{f,o} = 0.00010\text{ m}^2\cdot\text{K/W}$
  • Tube-side fouling: $R_{f,i}(r_o / r_i) = 0.00025 \times 1.1493 = 0.0002873\text{ m}^2\cdot\text{K/W}$

Total fouling resistance: Rf,total=0.00010+0.0002873=0.0003873 m2K/WR_{f,total} = 0.00010 + 0.0002873 = \mathbf{0.0003873}\text{ m}^2\cdot\text{K/W}

Total dirty thermal resistance: 1Udirty=1Uclean+Rf,total=0.0004083+0.0003873=0.0007956 m2K/W\frac{1}{U_{dirty}} = \frac{1}{U_{clean}} + R_{f,total} = 0.0004083 + 0.0003873 = \mathbf{0.0007956}\text{ m}^2\cdot\text{K/W}

Udirty=10.0007956=1,256.9 W/(m2K)U_{dirty} = \frac{1}{0.0007956} = \mathbf{1,256.9}\text{ W/(m}^2\cdot\text{K)}

Notice that fouling cuts the overall heat transfer coefficient nearly in half (from $2,449\text{ W/(m}^2\cdot\text{K)}$ down to $1,257\text{ W/(m}^2\cdot\text{K)}$)!


Step 5: Temperature Driving Force and Required Area

Because the shell side is isothermal ($T_{sat} = 60.0^\circ\text{C}$), $C_r = 0$ and $F = 1.00$:

  • $\Delta T_1 = T_{sat} - t_1 = 60.0 - 20.0 = 40.0^\circ\text{C}$
  • $\Delta T_2 = T_{sat} - t_2 = 60.0 - 35.0 = 25.0^\circ\text{C}$

ΔTlm=40.025.0ln(40.0/25.0)=15.0ln(1.6000)=15.00.47000=31.915C\Delta T_{lm} = \frac{40.0 - 25.0}{\ln(40.0 / 25.0)} = \frac{15.0}{\ln(1.6000)} = \frac{15.0}{0.47000} = \mathbf{31.915}^\circ\text{C}

Required outside design heat transfer area: Ao=QUdirtyΔTlm=9.825×106 W(1,256.9 W/(m2K))×(31.915 K)=9,825,00040,114=244.92 m2A_o = \frac{Q}{U_{dirty} \cdot \Delta T_{lm}} = \frac{9.825 \times 10^6\text{ W}}{(1,256.9\text{ W/(m}^2\cdot\text{K)}) \times (31.915\text{ K})} = \frac{9,825,000}{40,114} = \mathbf{244.92}\text{ m}^2


Step 6: Tube Count and Overdesign Percentage

Outside surface area of a single tube: Asingle=πdoL=π×(0.0254 m)×(4.50 m)=0.35908 m2A_{\text{single}} = \pi \cdot d_o \cdot L = \pi \times (0.0254\text{ m}) \times (4.50\text{ m}) = \mathbf{0.35908}\text{ m}^2

Total tubes required: Nt=AoAsingle=244.92 m20.35908 m2/tube=682.08683 tubesN_t = \frac{A_o}{A_{\text{single}}} = \frac{244.92\text{ m}^2}{0.35908\text{ m}^2/\text{tube}} = 682.08 \to \mathbf{683}\text{ tubes}

Clean surface area required: Aclean=QUcleanΔTlm=9,825,0002,449.2×31.915=9,825,00078,166=125.69 m2A_{clean} = \frac{Q}{U_{clean} \cdot \Delta T_{lm}} = \frac{9,825,000}{2,449.2 \times 31.915} = \frac{9,825,000}{78,166} = \mathbf{125.69}\text{ m}^2

Fouling overdesign margin: Overdesign %=(AoAcleanAclean)×100%=(244.92125.69125.69)×100%=+94.86%\text{Overdesign \%} = \left( \frac{A_o - A_{clean}}{A_{clean}} \right) \times 100\% = \left( \frac{244.92 - 125.69}{125.69} \right) \times 100\% = \mathbf{+94.86}\%

(To accommodate long-term fouling, the condenser requires almost double ($+94.9%$) the surface area of a brand-new clean unit).


6. Critical PE Exam Traps & Pitfalls

Trap 1: Omitting the Area Ratio ($r_o / r_i$) on Inside Fouling Resistances
When calculating $U_o$ based on outside area, internal resistances ($1/h_i$ and $R_{f,i}$) must be multiplied by $(A_o / A_i) = (r_o / r_i) = (d_o / d_i)$. Omitting this geometric correction factor introduces a $10%$ to $25%$ error on thick-walled tubes, leading to significant undersizing.

Trap 2: Assuming Transition Boiling is Stable in Fired Reboilers
On qualitative questions, candidates often confuse temperature-controlled boiling (steam reboilers) with flux-controlled boiling (fired heaters). Never design a fired heater or electric reboiler to operate near or past Critical Heat Flux (CHF). The negative slope ($dq''/d\Delta T_e < 0$) causes an uncontrollable temperature excursion into the film boiling regime, rupturing tube metallurgy.

Trap 3: Sizing Condensers Assuming Dropwise Condensation
Although dropwise condensation delivers heat transfer coefficients up to $100,000\text{ W/(m}^2\cdot\text{K)}$, you must never use dropwise correlations for industrial design unless the prompt explicitly specifies non-degrading hydrophobic promoter treatment. Industrial designs are always sized using Nusselt filmwise theory.

Test Your Knowledge

A newly installed shell-and-tube reboiler operating with clean surfaces exhibits an overall heat transfer coefficient of U_clean = 850 W/(m²K). Process specifications require accounting for a tube-side cooling water fouling resistance of R_f,i = 0.00025 m²K/W and a shell-side boiling hydrocarbon fouling resistance of R_f,o = 0.00020 m²*K/W (both already referenced to the outside tube area). What is the design fouled overall coefficient (U_dirty) and what percentage increase in heat transfer area is required to accommodate this fouling?

A
B
C
D
Test Your Knowledge

In a high-heat-flux fired reboiler tube where heat flux q'' is externally imposed by combustion gas radiation, what physical mechanism causes the negative slope (dq''/dDelta_T_e < 0) in the transition boiling regime, and what catastrophic event occurs if the heat flux exceeds the Critical Heat Flux (CHF)?

A
B
C
D
Test Your Knowledge

A chemical process condenser must condense saturated organic vapor at atmospheric pressure. The design engineer evaluates two bundle configurations using the same 1-inch OD (d_o = 0.0254 m) tubes of length L = 3.60 m: Configuration A places the tubes vertically, while Configuration B places the tubes horizontally in a single-row arrangement. According to Nusselt film condensation theory, what is the ratio of the average condensation heat transfer coefficient of the horizontal tube to that of the vertical tube (h_horiz / h_vert)?

A
B
C
D