Process Tools

Engineering Calculation Tools

Distillation Column

Multi-Stage Column Design

Calculator

Distillation Column Design Theory

The Fenske-Underwood-Gilliland (FUG) shortcut method for designing multi-component distillation columns with rigorous analytical equations.

Overview

The FUG method is a three-step shortcut procedure for preliminary distillation column design:

  1. Fenske Equation: Calculate minimum stages at total reflux
  2. Underwood Equations: Determine minimum reflux ratio
  3. Gilliland Correlation: Find actual stages for a given reflux ratio

Fenske Equation (Minimum Stages)

At total reflux (R → ∞), the minimum number of stages is:

Nmin = ln[(xD,LK/xD,HK) × (xB,HK/xB,LK)] / ln(αavg)

Where:

  • Nmin = Minimum number of equilibrium stages
  • xD,LK, xD,HK = Distillate mole fractions of light and heavy keys
  • xB,LK, xB,HK = Bottoms mole fractions of light and heavy keys
  • αavg = Average relative volatility between keys

Underwood Equations (Minimum Reflux)

The minimum reflux ratio requires solving two coupled equations:

First Equation (Find θ):

1 - q = Σ (αi × zi) / (αi - θ)

Second Equation (Find Rmin):

Rmin + 1 = Σ (αi × xD,i) / (αi - θ)

Where θ is the Underwood root (αHK < θ < αLK) and q is the feed quality.

Gilliland Correlation

Relates actual stages to minimum stages using dimensionless variables:

Define: X = (R - Rmin) / (R + 1)

Then: Y = (N - Nmin) / (N + 1)

Molokanov correlation (curve fit):

Y = 1 - exp[(1 + 54.4X) / (11 + 117.2X) × (X - 1) / √X]

Feed Quality (q)

The feed quality parameter q represents the thermal condition:

q ValueFeed Condition
q > 1Subcooled liquid
q = 1Saturated liquid (bubble point)
0 < q < 1Two-phase (partially vaporized)
q = 0Saturated vapor (dew point)
q < 0Superheated vapor

Design Procedure

  1. Specify feed composition, flow rate, and thermal condition
  2. Select light key (LK) and heavy key (HK) components
  3. Specify desired recoveries or product compositions
  4. Calculate average relative volatility
  5. Use Fenske equation to find Nmin
  6. Use Underwood equations to find Rmin
  7. Select operating reflux (typically 1.2-1.5 × Rmin)
  8. Use Gilliland correlation to find actual stages N
  9. Determine feed stage location using Kirkbride correlation

Kirkbride Equation (Feed Stage)

Estimates the optimal feed stage location:

NR/NS = [(zHK/zLK) × (xB,LK/xD,HK)² × (B/D)]0.206

Where NR is stages in rectifying section and NS is stages in stripping section.

References

  • Fenske, M.R. (1932). "Fractionation of Straight-Run Pennsylvania Gasoline". Ind. Eng. Chem.
  • Underwood, A.J.V. (1948). "Fractional Distillation of Multicomponent Mixtures". Chem. Eng. Prog.
  • Gilliland, E.R. (1940). "Multicomponent Rectification". Ind. Eng. Chem.
  • Seader, J.D. and Henley, E.J. (2006). Separation Process Principles. Wiley.