Gas Processing

Distillation & Fractionation Design

The shortcut FUG method that sizes a fractionation column โ€” minimum stages, minimum reflux, the reflux/stage trade-off โ€” plus the tray hydraulics that set the column's diameter.

1. The FUG shortcut method

Before rigorous tray-by-tray simulation, columns are scoped with the Fenske-Underwood-Gilliland (FUG) shortcut: Fenske gives the minimum number of stages at total reflux, Underwood gives the minimum reflux ratio at infinite stages, and Gilliland correlates the real column operating between those two limits. It is a screening tool โ€” fast, transparent, and good for first sizing โ€” but it assumes constant relative volatility and constant molar overflow.

2. Fenske โ€” minimum stages

Nmin = ln[ (xLK/xHK)D ยท (xHK/xLK)B ] / ln(ฮฑLK,HK)

where LK/HK are the light and heavy key components, the subscripts D/B the distillate and bottoms, and ฮฑ the relative volatility between the keys (geometric-mean across the column). Nmin is the stage count at total reflux โ€” the theoretical floor.

3. Underwood โ€” minimum reflux

Underwood solves for the root ฮธ (lying between the key components' volatilities) of:

ฮฃi ฮฑiยทzi,F / (ฮฑi โˆ’ ฮธ) = 1 โˆ’ q,   then   Rmin + 1 = ฮฃi ฮฑiยทxi,D / (ฮฑi โˆ’ ฮธ)

where zi,F are feed compositions, xi,D distillate compositions, and q the feed thermal condition. Caveat: for a true multicomponent separation the first equation can have multiple roots and distributing non-key components; collapsing it to a single binary key-only root (and clamping it) can under-predict Rmin and therefore under-size the column. Treat the binary form as a screening approximation and confirm wide or non-sharp splits in a rigorous simulator.

4. Gilliland โ€” actual stages

Real columns run at R = (1.1โ€“1.5)ยทRmin. The Gilliland correlation ties the operating reflux to the actual stage count:

Y = (N โˆ’ Nmin)/(N + 1),   X = (R โˆ’ Rmin)/(R + 1),   Y โ‰ˆ 1 โˆ’ exp[ (1+54.4X)/(11+117.2X) ยท (Xโˆ’1)/โˆšX ]

(Molokanov form). Kirkbride's equation then estimates the feed-stage location. The reflux choice is the central economic trade-off: more reflux โ†’ fewer trays (smaller vessel) but more reboiler/condenser duty and a larger diameter (higher vapour load).

5. Condensate stabilization

A condensate stabilizer is a fractionation column (often reboiled, sometimes cold-feed) that strips light ends (Cโ‚โ€“Cโ‚„, Hโ‚‚S) from raw condensate to hit a Reid Vapor Pressure (RVP) or true-vapor-pressure spec so the liquid can be stored/transported safely without flashing. It is sized by the FUG/hydraulics methods above, with the bottoms RVP as the separation target and overhead recovery balancing product loss against vapour-handling.

6. Tray hydraulics

Stage count comes from FUG; column diameter comes from tray hydraulics. The vapour capacity uses the same Souders-Brown form as separators, with a system/flooding factor:

UN,flood = CSB ยท โˆš( (ฯL โˆ’ ฯV)/ฯV ),   design at ~80โ€“85% of flood

where CSB (the capacity factor) depends on tray spacing and the flow-parameter FLV = (L/V)ยทโˆš(ฯV/ฯL) (Fair correlation). The operating window is bounded above by jet flooding and downcomer flooding/backup, and below by weeping/dumping; weir loading and downcomer residence time must also be checked. A tray that passes the stage count but fails hydraulics either floods or weeps.

7. References

  • GPSA Engineering Data Book (14th Ed) โ€” Section 19, Fractionation & Absorption (FUG method, tray hydraulics, Fair correlation).
  • Fenske (1932); Underwood (1948); Gilliland (1940); Molokanov et al. โ€” Gilliland correlation fit; Kirkbride โ€” feed-stage location.
  • Fair, J.R. โ€” tray flooding / capacity-factor correlation.
  • King, C.J., Separation Processes; Kister, H.Z., Distillation Design.

Frequently Asked Questions

What is the Fenske-Underwood-Gilliland (FUG) shortcut method?

The FUG shortcut scopes a fractionation column before rigorous simulation: Fenske gives the minimum number of stages at total reflux, Underwood gives the minimum reflux ratio at infinite stages, and Gilliland correlates the real column operating between those two limits. It assumes constant relative volatility and constant molar overflow.

Why can the binary form of the Underwood equation under-size a column?

For a true multicomponent separation, the Underwood feed equation can have multiple roots with distributing non-key components. Collapsing it to a single binary key-only root and clamping it can under-predict the minimum reflux and therefore under-size the column. Confirm wide or non-sharp splits in a rigorous simulator.

What is a condensate stabilizer?

A condensate stabilizer is a fractionation column that strips light ends (C1-C4, H2S) from raw condensate to meet a Reid Vapor Pressure (RVP) or true-vapor-pressure spec, so the liquid can be stored and transported safely without flashing.

What sets a distillation column's diameter?

Stage count comes from the FUG method, but column diameter comes from tray hydraulics. Vapor capacity uses the Souders-Brown form with a flooding factor; the operating window is bounded above by jet flooding and downcomer flooding/backup, and below by weeping and dumping. Design is typically at about 80-85% of flood.