Distillation Column Design — McCabe-Thiele Method

🧪 How a Distillation Column Works

Distillation separates a liquid mixture into its components by exploiting differences in volatility. The more volatile component (lower boiling point) concentrates in the vapour phase, while the less volatile component concentrates in the liquid phase. By contacting vapour and liquid repeatedly across multiple trays or packing sections, the separation is progressively enriched.

// COLUMN SECTIONS

🔵 CONDENSER

Distillate (D) + Reflux (L) withdrawn here

RECTIFYING SECTION (above feed)

Enriches more volatile component in vapour

⟶ FEED TRAY (F, z_F)

Optimal location minimises stages

STRIPPING SECTION (below feed)

Strips volatile component from liquid

🔴 REBOILER

Bottoms (B) withdrawn here

Key Design Variables

Symbol	Parameter	Defined By
F, z_F	Feed flowrate and composition	Process specification
D, x_D	Distillate flowrate and composition	Product specification
B, x_B	Bottoms flowrate and composition	Product specification
R = L/D	Reflux ratio	Designer choice (R > R_min)
N	Number of theoretical stages	McCabe-Thiele construction
E_o	Overall tray efficiency	Correlation or vendor data

⚖️ Overall Material Balance

Before drawing the McCabe-Thiele diagram, the overall material balance must be solved to find distillate and bottoms flowrates:

// Overall Material Balance

F = D + B

F = Feed | D = Distillate | B = Bottoms (all in mol/h or kmol/h)

// Component Balance (more volatile component)

F × z_F = D × x_D + B × x_B

z_F = Feed mole fraction | x_D = Distillate mole fraction | x_B = Bottoms mole fraction

// Solving for D and B

D = F × (z_F − x_B) / (x_D − x_B)

B = F − D

✅

Always verify: Check that D × x_D + B × x_B = F × z_F within ±0.1%. Any discrepancy means an arithmetic error before you proceed to the McCabe-Thiele construction.

📈 VLE Equilibrium Curve & Relative Volatility

The McCabe-Thiele method requires the vapour-liquid equilibrium (VLE) curve for the binary system. This plots vapour composition (y) vs liquid composition (x) of the more volatile component at the column operating pressure.

Relative Volatility (α)

For ideal or near-ideal mixtures, the equilibrium relationship is described by the relative volatility α:

// Relative Volatility

α = (y/x) / ((1−y)/(1−x)) = p_A* / p_B*

α = relative volatility of A relative to B | p_A*, p_B* = vapour pressures of A and B at the same temperature | α > 1 always (A is more volatile)

// Equilibrium y from x (constant α)

y* = α×x / (1 + (α−1)×x)

This is the equilibrium curve equation used to construct the VLE curve on the McCabe-Thiele diagram when α is approximately constant.

Typical Relative Volatility Values

Binary System	α (approx.)	Separation Difficulty
Ethanol – Water	1.5 – 2.0	Moderate (azeotrope forms at 89.4 mol%)
Benzene – Toluene	2.3 – 2.5	Relatively easy — nearly ideal
n-Heptane – n-Octane	1.8 – 2.2	Moderate — paraffinic system
Methanol – Water	3.0 – 4.0	Relatively easy
Propane – n-Butane	4.0 – 5.0	Easy — LPG fractionation
α ≈ 1.0 – 1.1	~1.05	Very difficult — needs many stages or extractive distillation

⚠️

Non-ideal systems: For systems that deviate significantly from Raoult's Law (e.g. ethanol-water, acetone-water), α is not constant and the equilibrium curve must be constructed from experimental VLE data or thermodynamic models (NRTL, UNIQUAC, UNIFAC). The McCabe-Thiele method still applies — only the equilibrium curve changes shape.

📐 Operating Lines

The McCabe-Thiele diagram has two straight operating lines — one for the rectifying section (above feed) and one for the stripping section (below feed). These lines represent the material balance on each tray.

Rectifying Section Operating Line (ROL)

// Rectifying Operating Line

y = (R/(R+1)) × x + x_D/(R+1)

R = L/D = reflux ratio | x_D = distillate composition | Slope = R/(R+1) | y-intercept = x_D/(R+1)

Stripping Section Operating Line (SOL)

// Stripping Operating Line

y = ((L + qF)/(V − (1−q)F)) × x − (B × x_B)/(V − (1−q)F)

Simplified: passes through (x_B, x_B) on the 45° diagonal and intersects the ROL at the q-line. In practice, draw it from (x_B, x_B) to the intersection of ROL and q-line.

🎯

Practical construction tip: Draw the ROL first (from x_D on the 45° line, with slope R/(R+1)). Then draw the q-line. The SOL starts at (x_B, x_B) and passes through the intersection of the ROL and q-line. You never need to calculate the SOL equation explicitly.

🔀 The q-Line — Feed Condition

The q-line (or feed line) represents the thermal condition of the feed. It passes through (z_F, z_F) on the 45° diagonal and has a slope that depends on the feed quality q:

// q-Line Equation

y = (q/(q−1)) × x − z_F/(q−1)

z_F = Feed mole fraction | q = Feed thermal condition (see below) | Slope = q/(q−1)

Feed Condition (q) — All Five Cases

q > 1

Subcooled Liquid

Feed is below its bubble point. q-line slopes steeply upward (left of vertical).

q = 1

Saturated Liquid

Feed is at bubble point. q-line is vertical. Most common feed condition.

0 < q < 1

Partially Vaporised

Feed is a vapour-liquid mixture. q-line slopes negatively downward.

q = 0

Saturated Vapour

Feed is at dew point. q-line is horizontal (slope = 0).

q < 0

Superheated Vapour

Feed is above dew point. q-line slopes positively upward (right of vertical).

// q Value for Subcooled Liquid Feed

q = 1 + (Cp_L × (T_bubble − T_feed)) / λ_feed

Cp_L = liquid heat capacity (kJ/kmol·K) | T_bubble = bubble point of feed | T_feed = actual feed temperature | λ_feed = molar latent heat of vaporisation (kJ/kmol)

🪜 Stepping Off Theoretical Stages

Once the equilibrium curve, ROL, SOL, and q-line are drawn on the x-y diagram, the number of theoretical stages is counted by "stepping off" stages between the equilibrium curve and operating lines — starting at x_D and working down to x_B.

1

Start at point (x_D, x_D) on the 45° diagonal

This represents the top of the column — distillate composition on both axes.
2

Draw a horizontal line to the equilibrium curve

Move left horizontally from (x_D, x_D) until you hit the equilibrium curve. This represents one theoretical stage (tray 1 — the condenser).
3

Draw a vertical line down to the operating line

From the equilibrium curve point, drop vertically to the ROL. This gives the liquid composition leaving that tray.
4

Switch to SOL at the feed tray

When the stepping crosses the intersection of ROL and q-line, switch from stepping off the ROL to stepping off the SOL. The tray where this switch occurs is the optimal feed tray.
5

Continue until x_B is reached or passed

Count each horizontal step as one theoretical stage. Stop when the staircase reaches or goes below x_B. Total count = N theoretical stages (including reboiler).

🎯

The reboiler counts as one theoretical stage. So the number of actual trays required = N_theoretical − 1 (before applying tray efficiency). Switch to the SOL as early as possible — this gives the minimum number of stages.

♻️ Reflux Ratio — Minimum, Actual & Optimal

The reflux ratio R = L/D is the most important operating variable in distillation. It directly controls the tradeoff between capital cost (number of trays) and operating cost (energy).

Minimum Reflux Ratio (R_min)

At minimum reflux, the operating lines pinch the equilibrium curve — requiring an infinite number of stages. R_min is found graphically as the slope of the ROL that just touches (pinches) the equilibrium curve at the feed intersection:

// R_min from Pinch Point

R_min/(R_min+1) = (x_D − y_pinch) / (x_D − x_pinch)

y_pinch, x_pinch = equilibrium curve coordinates at the pinch point (intersection of q-line with equilibrium curve for most cases)

// Underwood Equation (for constant α, saturated liquid feed)

R_min = (1/(α−1)) × (x_D/z_F − α×(1−x_D)/(1−z_F))

Valid for ideal binary systems with saturated liquid feed (q=1) and constant relative volatility α

Actual Reflux Ratio

// Design Reflux Ratio

R = 1.2 to 1.5 × R_min (typical design rule)

R = 1.2 × R_min: fewer trays, higher energy cost. R = 1.5 × R_min: more trays, lower energy. Optimum usually at R = 1.2–1.3 × R_min for most columns.

Effect of Reflux Ratio on Design

Reflux Ratio R	Number of Stages	Energy Consumption	Column Diameter
R = R_min	Infinite	Minimum	Smallest
R = 1.2 × R_min	High but finite	Low	Moderate
R = 1.3 × R_min	Moderate	Moderate	Moderate — optimal zone
R = 1.5 × R_min	Lower	Higher	Larger
R = Total reflux	Minimum stages	Maximum	Largest

⚙️ Tray Efficiency & Column Sizing

Theoretical stages assume perfect vapour-liquid equilibrium on every tray. Real trays never achieve this — tray efficiency accounts for the actual separation achieved per real tray.

// Actual Number of Trays

N_actual = N_theoretical / E_o

E_o = Overall tray efficiency (typically 0.60–0.85 for sieve trays in most systems) | N_theoretical = stages from McCabe-Thiele (excluding reboiler)

Tray Types and Typical Efficiencies

⬡

Sieve Trays

Perforated plate. Simple, low cost, easy to clean. E_o = 60–80%. Most widely used.

🔘

Bubble Cap Trays

Best turndown ratio, no weeping. E_o = 70–85%. Higher cost, older technology.

🎯

Valve Trays

Moving valves — good turndown and efficiency. E_o = 70–85%. Standard in modern columns.

🌀

Structured Packing

Low pressure drop, high efficiency, no stage count — uses HETP instead. Ideal for vacuum and low-flow duties.

O'Connell Correlation (E_o Estimation)

// O'Connell Correlation

E_o = 51 − 32.5 × log(μ_L × α)

μ_L = liquid viscosity at average column temperature (cP) | α = relative volatility | E_o in % | Valid range: E_o = 40–90%

Column Diameter — Souders-Brown

// Column Diameter from Vapour Flowrate

u_flood = C_sb × √[(ρ_L − ρ_V) / ρ_V]

A_col = V_vol / (u_flood × f)

D_col = √(4 × A_col / π)

C_sb = Souders-Brown coefficient (0.06–0.12 m/s for sieve trays) | f = fraction of flooding (design at 70–80%) | V_vol = vapour volumetric flowrate (m³/s)

🧮 Complete Worked Example — Benzene/Toluene

System: Separate a benzene-toluene mixture. Feed = 100 kmol/h at bubble point (q = 1). z_F = 0.45 (benzene). x_D = 0.95. x_B = 0.05. α = 2.45. Design at R = 1.3 × R_min.

Step 1 — Overall Material Balance

Calculation	Result
D = F(z_F − x_B)/(x_D − x_B) = 100(0.45−0.05)/(0.95−0.05)	D = 44.4 kmol/h
B = F − D = 100 − 44.4	B = 55.6 kmol/h
Check: D×x_D + B×x_B = 44.4×0.95 + 55.6×0.05	= 44.9 ≈ 45 ✅

Step 2 — Minimum Reflux (Underwood, q=1)

// R_min

R_min = (1/(2.45−1)) × (0.95/0.45 − 2.45×0.05/0.55)

R_min = (1/1.45) × (2.111 − 0.222) = 0.690 × 1.889 = 1.30

Step 3 — Actual Reflux Ratio

// Operating Reflux

R = 1.3 × R_min = 1.3 × 1.30 = 1.69

ROL slope = R/(R+1) = 1.69/2.69 = 0.628

ROL y-intercept = x_D/(R+1) = 0.95/2.69 = 0.353

Step 4 — q-Line (Saturated Liquid, q=1)

// q-Line

q = 1 → vertical line at x = z_F = 0.45

Pinch at x=0.45 on equilibrium curve: y* = 2.45×0.45/(1+1.45×0.45) = 1.1025/1.6525 = 0.667

Step 5 — Theoretical Stages (McCabe-Thiele)

Stage	x (liquid leaving)	y* (equilibrium)	Section
Condenser	0.950	0.950	Start (45° line)
Stage 1	0.881	0.950	Rectifying
Stage 2	0.779	0.881	Rectifying
Stage 3	0.635	0.779	Rectifying
Stage 4 (Feed)	0.450	0.635	Feed tray — switch to SOL
Stage 5	0.280	0.450	Stripping
Stage 6	0.130	0.280	Stripping
Reboiler (Stage 7)	0.050	0.130	Reboiler = 1 stage

Step 6 — Actual Trays (Tray Efficiency)

// Actual Tray Count

N_theoretical = 7 (including reboiler) → 6 trays above reboiler

E_o = 75% (benzene-toluene, sieve trays)

N_actual = 6 / 0.75 = 8 actual trays + reboiler

Feed tray = Tray 4 theoretical → actual feed nozzle at approximately tray 5–6 from top

Summary

Parameter	Result	Unit
Distillate (D)	44.4	kmol/h
Bottoms (B)	55.6	kmol/h
Minimum reflux (R_min)	1.30	—
Design reflux (R)	1.69	—
Theoretical stages (N)	7	(incl. reboiler)
Overall tray efficiency (E_o)	75%	—
Actual trays required	8	trays + reboiler
Optimal feed tray	Stage 4 (tray 5–6 from top)	—

Distillation Column Design — McCabe-Thiele Method Step by Step

🧪 How a Distillation Column Works

// COLUMN SECTIONS

Key Design Variables

⚖️ Overall Material Balance

📈 VLE Equilibrium Curve & Relative Volatility

Relative Volatility (α)

Typical Relative Volatility Values

📐 Operating Lines

Rectifying Section Operating Line (ROL)

Stripping Section Operating Line (SOL)

🔀 The q-Line — Feed Condition

Feed Condition (q) — All Five Cases

🪜 Stepping Off Theoretical Stages

Start at point (x_D, x_D) on the 45° diagonal

Draw a horizontal line to the equilibrium curve

Draw a vertical line down to the operating line

Switch to SOL at the feed tray

Continue until x_B is reached or passed

♻️ Reflux Ratio — Minimum, Actual & Optimal

Minimum Reflux Ratio (R_min)

Actual Reflux Ratio

Effect of Reflux Ratio on Design

⚙️ Tray Efficiency & Column Sizing

Tray Types and Typical Efficiencies

Sieve Trays

Bubble Cap Trays

Valve Trays

Structured Packing

O'Connell Correlation (E_o Estimation)

Column Diameter — Souders-Brown

🧮 Complete Worked Example — Benzene/Toluene

Step 1 — Overall Material Balance

Step 2 — Minimum Reflux (Underwood, q=1)

Step 3 — Actual Reflux Ratio

Step 4 — q-Line (Saturated Liquid, q=1)

Step 5 — Theoretical Stages (McCabe-Thiele)

Step 6 — Actual Trays (Tray Efficiency)

Summary

Practice More Distillation Problems — GATE Style