PracticeVideosBlogFlashcards

9.3 Boiling Regimes, Condensation Phenomena, and Fouling Resistances

Key Takeaways

  • The overall thermal resistance network sums five series resistances: 1/(U_oA_o) = 1/(h_iA_i) + R_f,i/A_i + ln(r_o/r_i)/(2pikL) + R_f,o/A_o + 1/(h_oA_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 (hh) 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 (RfR_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 (AoA_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 AoA_o (Ao/Ai=ro/ri=do/diA_o / A_i = r_o / r_i = d_o / d_i) yields the design equation for the outside overall coefficient (UoU_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:

  • UoU_o = overall heat transfer coefficient referenced to outside tube area (W/(m2⋅K)\text{W/(m}^2\cdot\text{K)}).
  • hi,hoh_i, h_o = inner and outer convective heat transfer coefficients (W/(m2⋅K)\text{W/(m}^2\cdot\text{K)}).
  • Rf,i,Rf,oR_{f,i}, R_{f,o} = inner and outer fouling resistances (m2⋅K/W\text{m}^2\cdot\text{K/W} or hr⋅ft2⋅∘F/Btu\text{hr}\cdot\text{ft}^2\cdot^\circ\text{F/Btu}).
  • ri,ror_i, r_o = inner and outer tube radii (di/2,do/2d_i/2, d_o/2).
  • kwk_w = thermal conductivity of the tube wall material (W/(m⋅K)\text{W/(m}\cdot\text{K)}).
  • LL = total tube length (m\text{m}).

Clean vs. Dirty Overall Heat Transfer Coefficients

In process design, the clean overall heat transfer coefficient (UcleanU_{clean}) represents newly fabricated, un-fouled equipment (Rf=0R_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 (UdirtyU_{dirty}) incorporates the total combined fouling resistance (Rf,totalR_{f,total}):

1Udirty=1Uclean+Rf,total  ⟹  Udirty=Uclean1+Uclean⋅Rf,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 Rf,total=Rf,o+Rf,i(ro/ri)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 UcleanU_{clean} to UdirtyU_{dirty}, the required surface area must be expanded to guarantee performance throughout the operating cycle between turnarounds:

AdesignAclean=UcleanUdirty=1+Uclean⋅Rf,total\frac{A_{design}}{A_{clean}} = \frac{U_{clean}}{U_{dirty}} = 1 + U_{clean} \cdot R_{f,total} Overdesign %=(Adesign−AcleanAclean)×100%=Uclean⋅Rf,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 RfR_f (m2⋅K/W\text{m}^2\cdot\text{K/W})US Customary RfR_f (hr⋅ft2⋅∘F/Btu\text{hr}\cdot\text{ft}^2\cdot^\circ\text{F/Btu})
Treated Boiler Feedwater (<50∘C<50^\circ\text{C})0.00009−0.000180.00009 - 0.000180.0005−0.00100.0005 - 0.0010
Treated Cooling Tower Water (<50∘C<50^\circ\text{C})0.00018−0.000350.00018 - 0.000350.0010−0.00200.0010 - 0.0020
River Water / Untreated Surface Water0.00035−0.000530.00035 - 0.000530.0020−0.00300.0020 - 0.0030
Clean Light Hydrocarbons (Naphtha, LPG)0.000180.000180.00100.0010
Heavy Gas Oil / Fuel Oil0.00053−0.000880.00053 - 0.000880.0030−0.00500.0030 - 0.0050
Crude Oil Bottoms / Vacuum Resid0.00088−0.001760.00088 - 0.001760.0050−0.01000.0050 - 0.0100
Dry Gas (Air, Nitrogen, Natural Gas)0.00018−0.000350.00018 - 0.000350.0010−0.00200.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 (ΔTe\Delta T_e):

ΔTe≡Ts−Tsat\Delta T_e \equiv T_s - T_{sat}

Plotting heat flux (q′′=Q/Aq'' = Q/A) versus ΔTe\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 (ΔTe≲5∘C\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′′∝ΔTe5/4q'' \propto \Delta T_e^{5/4}).

Regime II: Nucleate Boiling (5∘C≲ΔTe≲30∘C5^\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∘C<ΔTe<10∘C5^\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∘C<ΔTe<30∘C10^\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′′∝ΔTe3q'' \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(Ts−Tsat)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 CsfC_{sf} is an empirical surface-fluid constant reflecting surface wettability (e.g., 0.0130.013 for water-copper; 0.0060.006 for water-brass), and n=1.0n = 1.0 for water (1.71.7 for organic liquids).

The Critical Heat Flux (CHF / Burnout Point)

The peak of the nucleate boiling curve is the Critical Heat Flux (qmax′′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, qmax′′≈1.1 to 1.3 MW/m2q''_{max} \approx 1.1\text{ to }1.3\text{ MW/m}^2 (350,000 to 400,000 Btu/(hr⋅ft2)350,000\text{ to }400,000\text{ Btu/(hr}\cdot\text{ft}^2)).

Regime III: Transition Boiling (30∘C≲ΔTe≲120∘C30^\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 (kv≪klk_v \ll k_l), the vapor patches act as an insulating blanket.

Crucially, the transition boiling curve has a negative slope (dq′′/dΔTe<0dq'' / 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′′q'' is fixed by fuel firing rate. If heat flux is increased beyond qmax′′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∘C> 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 (ΔTe>120∘C\Delta T_e > 120^\circ\text{C})

The minimum heat flux (qmin′′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+hrad≈hconv+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 hrad=ϵemσSB(Ts4−Tsat4)/(Ts−Tsat)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 (Ts<TsatT_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,0002,000 to 10,000 W/(m2⋅K)10,000\text{ W/(m}^2\cdot\text{K)}.
  2. Dropwise Condensation: The liquid condensate does not wet the surface (contact angle θ>90∘\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,00030,000 to 100,000 W/(m2⋅K)100,000\text{ W/(m}^2\cdot\text{K)}—a 55- to 1010-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 (Ref<30Re_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 (TT decreases linearly across the film from TsatT_{sat} at the liquid-vapor interface to TsT_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: Tf=(Tsat+Ts)/2T_f = (T_{sat} + T_s) / 2.

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

δ(x)=[4klμl(Tsat−Ts)xgρl(ρl−ρv)hfg′]1/4∝x1/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 hfg′h_{fg}' is the modified latent heat of vaporization, corrected for liquid film subcooling (Rohsenow correction):

hfg′=hfg+0.68⋅Cp,l(Tsat−Ts)h_{fg}' = h_{fg} + 0.68 \cdot C_{p,l} (T_{sat} - T_s)

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

hˉvert=0.943[gρl(ρl−ρv)kl3hfg′μlL(Tsat−Ts)]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 dod_o, the condensate film drains around the circumference, traveling a characteristic drainage distance of only πdo/2\pi d_o / 2 before dripping off the bottom. Nusselt's derivation yields:

hˉhoriz=0.729[gρl(ρl−ρv)kl3hfg′μldo(Tsat−Ts)]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, ΔT\Delta T, and tube dimensions (do,Ld_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 mL = 4.0\text{ m} and do=0.0254 md_o = 0.0254\text{ m} (1.0 in1.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 33 times higher heat transfer coefficient than the same tube placed vertically! On a vertical tube, condensate drains down the entire length (4.0 m4.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 NN tubes, Nusselt derived the average bundle coefficient as:

hˉbundle=hˉsingle⋅N−1/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ˉsingle⋅N−1/6\bar{h}_{bundle, Kern} = \bar{h}_{single} \cdot N^{-1/6}

4. Summary Table of Boiling and Condensation Regimes & Equations

Phenomenon / RegimeGoverning MechanismCharacteristic ΔTe\Delta T_e RangeGoverning Correlation / Analytical ModelKey Engineering Design Rule
Natural Convection BoilingSingle-phase buoyant liquid circulationΔTe<5∘C\Delta T_e < 5^\circ\text{C}q′′=hncΔTe∝ΔTe5/4q'' = h_{nc} \Delta T_e \propto \Delta T_e^{5/4}Superheat insufficient to nucleate bubbles
Nucleate BoilingBubble nucleation, growth, and micro-convection5∘C<ΔTe<30∘C5^\circ\text{C} < \Delta T_e < 30^\circ\text{C}Rohsenow: q′′∝ΔTe3q'' \propto \Delta T_e^3Optimal design regime; highest heat flux per m2m^2
Critical Heat Flux (CHF)Vapor columns choke liquid replenishmentPeak at ΔTe≈30∘C\Delta T_e \approx 30^\circ\text{C}Zuber: qmax′′=0.149ρvhfg[σg(Δρ)/ρv2]1/4q''_{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 formation30∘C<ΔTe<120∘C30^\circ\text{C} < \Delta T_e < 120^\circ\text{C}Negative slope: dq′′/dΔTe<0dq''/d\Delta T_e < 0Unstable in flux-controlled heaters (burnout danger)
Film BoilingStable continuous vapor cushionΔTe>120∘C\Delta T_e > 120^\circ\text{C}h=hconv+34hradh = h_{conv} + \frac{3}{4} h_{rad}; qmin′′q''_{min} via LeidenfrostExtremely high wall temperatures (Ts>800∘CT_s > 800^\circ\text{C})
Vertical CondensationLaminar draining liquid filmN/A (Ts<TsatT_s < T_{sat})hˉvert=0.943[gρl(Δρ)kl3hfg′/(μlLΔT)]1/4\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∝L−1/4h \propto L^{-1/4}
Horizontal CondensationPeripheral drainage around circumferenceN/A (Ts<TsatT_s < T_{sat})hˉhoriz=0.729[gρl(Δρ)kl3hfg′/(μldoΔT)]1/4\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: hhoriz/hvert≈2.5−3.0h_{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 kg/h15,000\text{ kg/h} (4.167 kg/s4.167\text{ kg/s}) of pure saturated organic process vapor at Tsat=60.0∘CT_{sat} = 60.0^\circ\text{C} on the shell side.

Condensate & Vapor Properties at 60.0∘C60.0^\circ\text{C}:

  • Latent heat of vaporization: hfg=2,358 kJ/kg=2.358×106 J/kgh_{fg} = 2,358\text{ kJ/kg} = 2.358 \times 10^6\text{ J/kg}
  • Liquid condensate density: ρl=985 kg/m3\rho_l = 985\text{ kg/m}^3
  • Vapor density: ρv≈0.13 kg/m3\rho_v \approx 0.13\text{ kg/m}^3 (so ρl−ρv≈ρl\rho_l - \rho_v \approx \rho_l)
  • Condensate liquid thermal conductivity: kl=0.650 W/(m⋅K)k_l = 0.650\text{ W/(m}\cdot\text{K)}
  • Condensate liquid viscosity: μl=4.80×10−4 Pa⋅s\mu_l = 4.80 \times 10^{-4}\text{ Pa}\cdot\text{s}

Tube Bundle Specifications:

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

Operating & Utility Conditions:

  • Cooling water enters tubes at t1=20.0∘Ct_1 = 20.0^\circ\text{C} and exits at t2=35.0∘Ct_2 = 35.0^\circ\text{C}.
  • Tube-side convective heat transfer coefficient: hi=4,800.0 W/(m2⋅K)h_i = 4,800.0\text{ W/(m}^2\cdot\text{K)}
  • Tube-side cooling water fouling resistance: Rf,i=0.00025 m2⋅K/WR_{f,i} = 0.00025\text{ m}^2\cdot\text{K/W}
  • Shell-side organic fouling resistance: Rf,o=0.00010 m2⋅K/WR_{f,o} = 0.00010\text{ m}^2\cdot\text{K/W}
  • Average outer tube wall surface temperature: Ts≈42.0∘CT_s \approx 42.0^\circ\text{C} (giving ΔT=Tsat−Ts=60.0−42.0=18.0∘C\Delta T = T_{sat} - T_s = 60.0 - 42.0 = 18.0^\circ\text{C})

Calculate:

  1. The total condensation heat duty (QQ) in kW\text{kW}.
  2. The average shell-side condensation coefficient for a single horizontal tube (hˉsingle\bar{h}_{single}) using Nusselt theory, and the corrected bundle coefficient (hˉo\bar{h}_o).
  3. The clean overall heat transfer coefficient (UcleanU_{clean}) based on outside area.
  4. The design fouled overall heat transfer coefficient (UdirtyU_{dirty}) and total fouling resistance (Rf,totalR_{f,total}).
  5. The required outside heat transfer area (AoA_o) using ΔTlm\Delta T_{lm}.
  6. The required number of tubes (NtN_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 do=0.0254 md_o = 0.0254\text{ m} and ΔT=18.0 K\Delta T = 18.0\text{ K}:

hˉsingle=0.729[gρl(ρl−ρv)kl3hfgμldo(Tsat−Ts)]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/(m⋅K))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=μl⋅do⋅ΔT=(4.80×10−4 Pa⋅s)×(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×10−4\text{Denominator} = 2.1946 \times 10^{-4} Ratio=6.1633×10122.1946×10−4=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/(m2⋅K)\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=9N = 9 vertical rows (using Kern's N−1/6N^{-1/6} correction):

hˉo=hˉsingle⋅N−1/6=9,435.7×(9)−1/6=9,435.7×0.6934=6,542.7 W/(m2⋅K)\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 (UcleanU_{clean})

Radii and area ratios:

  • ro=do/2=0.0254/2=0.0127 mr_o = d_o / 2 = 0.0254 / 2 = 0.0127\text{ m}
  • ri=di/2=0.0221/2=0.01105 mr_i = d_i / 2 = 0.0221 / 2 = 0.01105\text{ m}
  • Ratio ro/ri=0.0127/0.01105=1.1493r_o / r_i = 0.0127 / 0.01105 = 1.1493

Evaluate individual clean thermal resistances referenced to AoA_o:

  1. Shell-side condensation resistance: Ro=1ho=16,542.7=0.0001528 m2⋅K/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 m2⋅K/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 m2⋅K/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 m2⋅K/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/(m2⋅K)U_{clean} = \frac{1}{0.0004083} = \mathbf{2,449.2}\text{ W/(m}^2\cdot\text{K)}

Step 4: Dirty Overall Coefficient (UdirtyU_{dirty})

Fouling resistances referenced to outside area:

  • Shell-side fouling: Rf,o=0.00010 m2⋅K/WR_{f,o} = 0.00010\text{ m}^2\cdot\text{K/W}
  • Tube-side fouling: Rf,i(ro/ri)=0.00025×1.1493=0.0002873 m2⋅K/WR_{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 m2⋅K/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 m2⋅K/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/(m2⋅K)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 W/(m2⋅K)2,449\text{ W/(m}^2\cdot\text{K)} down to 1,257 W/(m2⋅K)1,257\text{ W/(m}^2\cdot\text{K)})!


Step 5: Temperature Driving Force and Required Area

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

  • ΔT1=Tsat−t1=60.0−20.0=40.0∘C\Delta T_1 = T_{sat} - t_1 = 60.0 - 20.0 = 40.0^\circ\text{C}
  • ΔT2=Tsat−t2=60.0−35.0=25.0∘C\Delta T_2 = T_{sat} - t_2 = 60.0 - 35.0 = 25.0^\circ\text{C}
ΔTlm=40.0−25.0ln⁡(40.0/25.0)=15.0ln⁡(1.6000)=15.00.47000=31.915∘C\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/(m2⋅K))×(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=π⋅do⋅L=π×(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.08→683 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 %=(Ao−AcleanAclean)×100%=(244.92−125.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%+94.9\%) the surface area of a brand-new clean unit).


6. Critical PE Exam Traps & Pitfalls

Trap 1: Omitting the Area Ratio (ro/rir_o / r_i) on Inside Fouling Resistances
When calculating UoU_o based on outside area, internal resistances (1/hi1/h_i and Rf,iR_{f,i}) must be multiplied by (Ao/Ai)=(ro/ri)=(do/di)(A_o / A_i) = (r_o / r_i) = (d_o / d_i). Omitting this geometric correction factor introduces a 10%10\% to 25%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ΔTe<0dq''/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 W/(m2⋅K)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

U_dirty = 425 W/(m²*K), Area increase = 100.0%

B

U_dirty = 735 W/(m²*K), Area increase = 15.6%

C

U_dirty = 512 W/(m²*K), Area increase = 66.0%

D

U_dirty = 615 W/(m²*K), Area increase = 38.3%

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

An unstable vapor blanket intermittently covers the heating surface, replacing high-conductivity liquid with low-conductivity vapor; exceeding CHF causes a catastrophic jump in wall temperature to the film boiling regime, potentially melting the tube wall (burnout)

B

Liquid molecules dissociate into ions that suppress vapor pressure; exceeding CHF causes sudden explosive cavitation that implodes the tube bundle

C

The boundary layer transitions from turbulent to laminar flow due to thermophoresis; exceeding CHF causes the fluid to freeze onto the heat transfer surface

D

Nucleate bubble density drops to zero due to surface de-wetting; exceeding CHF causes heat flux to drop to zero and flow to stagnate

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

0.77

B

1.00

C

2.67

D

5.34

Sections you finish are checked off in the contents.