Everything, Everywhere
Verified Specification | Standardized Formulas | Instant Precision
Secure & Private (Zero Data Retention) Free Access • No Sign-Up

Distillation Column FUG Shortcut Sizing

Fenske Stages, Underwood Rmin, Gilliland N, Kirkbride Feed Tray & Column Diameter

Units:
Binary / Key System Preset
Relative Volatility (α_LK/HK)
Total Feed Rate (F)
Feed LK Mole Fraction (z_LK)
Distillate LK Purity (x_D,LK)
Bottoms LK Loss (x_B,LK)
Feed Thermal Condition (q)
Operating Reflux Factor (R / Rmin)
Column Top Pressure
Overall Tray Efficiency (E_o)
Tray Spacing
% of Flooding Design
Actual Theoretical Stages (N)
21.8
Minimum Stages Nmin: 9.3
Operating Reflux Ratio (R)
1.54
Underwood Rmin: 1.23
Optimum Feed Tray (NF)
Tray 11
Rect: 10 | Strip: 11
Column Inner Diameter
5.8 ft
80% Flood (Fair Csb)

Distillation Column Mass Balance & Sizing Summary

Distillate & Bottoms Rates:
D: 448.5 | B: 551.5 lbmol/h
Distillate Split: 44.8% of Feed
Actual Trays Required (Eo = 75%):
29 Actual Trays
Tangent-to-Tangent Height: ~68 ft
Condenser & Reboiler Traffic:
Vapor V: 1,139 lbmol/h
Reflux Liquid L: 691 lbmol/h
Gilliland Reflux vs Stages Operating Curve N vs R/Rmin
Fractionation Column Structural Schematic Feed Tray & Section Heights

Fatal Traps & Distillation Column Engineering Pitfalls

Trap 1: Close-Boiling Relative Volatility Pinch Point (α ≤ 1.15) Causing Stage Explosion

In close-boiling systems like propane/propylene (α = 1.12 to 1.15) or ethylbenzene/styrene (α = 1.30), junior engineers frequently underestimate how fast stages escalate. The denominator in Fenske is ln(α). As α approaches 1.0, ln(α) approaches zero, causing N_min to explode from 10 stages up to 90 to 180 stages. Furthermore, Underwood R_min skyrockets, demanding enormous vapor reboil and gigantic column diameters. Attempting to design a tight close-boiling splitter with standard rules without considering heat-pumped distillation or extractive dividing-wall columns (DWC) results in non-economic monsters that fail financial justification.

Trap 2: Ignoring Feed Tray Mismatch and Composition Shock

Introducing feed at an incorrect tray level introduces severe thermodynamic irreversibility and composition pinching. If a feed with 45% LK is introduced on a tray where internal liquid has already concentrated to 70% LK, the incoming feed severely dilutes the rectified liquid. This forces the rectifying section to re-separate material that was already refined, demanding excessive reflux and destroying separation efficiency. The Kirkbride equation must be rigorously checked, and industrial columns should always be fabricated with at least two alternate feed nozzles (one tray above and one tray below the design feed tray).

Trap 3: Jet Flooding and Downcomer Choking from Foaming Hydrocarbons

Fair's entrainment flooding correlation assumes clean, non-foaming systems. In amine treaters, sour water strippers, or heavy oil vacuum towers containing surface-active contaminants, stable foams develop on tray decks. Froth backs up into downcomers, choking liquid drainage and causing liquid to carry over into the overhead vapor line (jet flooding) at vapor velocities well below 60% of Fair's theoretical flood limit. Designers must apply an API 521 foaming system derating factor (typically F_foam = 0.70 to 0.85) to the allowable capacity factor C_sb.

Trap 4: Tray Weeping at Turndown Leading to Total Fractionation Loss

While maximum column diameter is governed by flooding at peak rates, minimum diameter is governed by weeping at low rates. Sieve trays rely on upward vapor kinetic energy (F-factor = u_hole * sqrt(rho_V)) to support the liquid pool on the active tray deck. If plant throughput is turned down to 50% or 60% of design, vapor velocity through tray holes drops below the weeping threshold. Liquid dumps directly through tray perforations without contacting vapor, causing tray efficiency to collapse from 80% to below 20% and dumping light components into bottoms. Sieve trays require blanking strips or replacement with movable valve trays for broad turndown flexibility.

Trap 5: Column Top-to-Bottom Temperature and Pressure Gradient Invalidation

The standard FUG shortcut assumes constant relative volatility (α) throughout the column. In reality, column bottoms operate at higher pressure (due to cumulative tray pressure drop of 3 to 5 mmHg per tray) and much higher temperature than the top condenser. Because relative volatility declines as temperature rises, α at the reboiler can be 20% to 30% lower than α at the condenser. Sizing a column using top-condenser α alone under-estimates required stages and reboiler duty by 25%+. Engineers must use the geometric mean relative volatility: α_mean = sqrt(α_top * α_bottom).

Comprehensive Fenske-Underwood-Gilliland Mathematical Derivations

Multicomponent fractionation shortcuts evaluate equilibrium stage requirements and vapor-liquid capacity using rigorous thermodynamic equations:

1. Overall Component Mass Balance

Overall Molar Balance: F = D + B
Light Key Balance: F * z_LK = D * x_D_LK + B * x_B_LK
Distillate Rate: D = F * (z_LK - x_B_LK) / (x_D_LK - x_B_LK)
Bottoms Rate: B = F - D

2. Minimum Stages at Total Reflux (Fenske Equation)

Separation Factor: S = (x_D_LK / (1 - x_D_LK)) * ((1 - x_B_LK) / x_B_LK)
Fenske Minimum Stages: N_min = ln(S) / ln(alpha_LK_HK)

3. Minimum Reflux Ratio (Underwood Equations)

Underwood Root theta (1.0 < theta < alpha):
  1 - q = (alpha * z_LK) / (alpha - theta) + (1 * (1 - z_LK)) / (1 - theta)
Minimum Reflux Ratio R_min:
  R_min + 1 = (alpha * x_D_LK) / (alpha - theta) + (1 * (1 - x_D_LK)) / (1 - theta)

4. Actual Theoretical Stages (Gilliland Correlation)

Operating Reflux: R = R_factor * R_min
Gilliland Parameter: X = (R - R_min) / (R + 1)
Stage Fraction: Y = 1 - exp( [ (1 + 54.4*X) / (11 + 117.2*X) ] * [ (X - 1) / sqrt(X) ] )
Actual Theoretical Stages: N = (N_min + Y) / (1 - Y)

5. Optimum Feed Tray (Kirkbride Equation) & Column Diameter

Ratio of Stages: ln(NR / NS) = 0.206 * ln[ (z_HK / z_LK) * (x_B_LK / (1 - x_D_LK))^2 * (B / D) ]
Fair Flooding Velocity: u_flood = C_sb * sqrt((rho_L - rho_V) / rho_V)
Column Internal Diameter: D_col = sqrt( (4 * V_m3_s) / (pi * f_flood * u_flood) )

Frequently Asked Questions

What is the Fenske-Underwood-Gilliland (FUG) shortcut distillation method? +
What are the Light Key (LK) and Heavy Key (HK) components in multicomponent fractionation? +
How does the feed thermal condition factor (q) impact column vapor and liquid traffic? +
How is the optimum economic operating reflux ratio selected? +
How does Fair's entrainment flooding correlation determine column diameter? +
Sponsored Utility
While You're Here
Sponsored Recommendations
Advertisement