Number of Stages by the Gilliland Correlation
Estimates the number of theoretical stages of a distillation column with the Gilliland correlation, which links excess reflux to excess stages: with X = (R − R_min)/(R + 1) and Y = (N − N_min)/(N + 1), you get N = (Y + N_min)/(1 − Y). It is the third step of the FUG shortcut method, after N_min from the Fenske equation and R_min from Underwood, and it settles in one line the preliminary sizing that would otherwise need a McCabe-Thiele diagram or a simulator. Eduljee's analytical fit is adopted, Y = 0.75·(1 − X^0.5668), the usual form for hand calculation; the Molokanov correlation is more accurate at the extremes and gives a result a few percent different. Enter the operating reflux ratio, the minimum reflux ratio and the minimum number of stages.
Result
—
Gilliland: theoretical stages without a McCabe-Thiele plot
The FUG shortcut settles a distillation column in three moves: Fenske gives the minimum stage count at total reflux, Underwood gives the minimum reflux ratio, and Gilliland ties both into the stage count at real operating reflux. This is the calculation a process engineer runs during a feasibility study, ahead of any simulator, to land on column height, tray count and a capital cost estimate that survives the first review meeting. Without it, picking a reflux ratio happens blind — and reflux decides how many trays get bought and how much steam gets paid for every year.
The correlation runs on two dimensionless coordinates: X = (R − R_min)/(R + 1), the excess reflux, and Y = (N − N_min)/(N + 1), the excess stages. Eduljee's analytical fit appears here, Y = 0.75·(1 − X^0.5668), the form used for hand calculation since 1975; inverted, N = (Y + N_min)/(1 − Y). Defaults of R = 1.5, R_min = 1 and N_min = 5 give X = 0.2, Y = 0.4488 and N = 9.885 theoretical stages. Economic design reflux typically lands between 1.1 and 1.5 times the minimum, precisely where the correlation was fitted. Molokanov's alternative tracks the extremes better and would return 10.12 in this same case, about 2% higher.
Eduljee's fit carries a known flaw at the lower bound: as R closes on R_min, X heads to zero and Y saturates at 0.75, so N stalls at 4·N_min + 3 — with N_min = 5, a ceiling of 23 stages, where physics demands N grow without limit. Feed R = 1.001 against R_min = 1 and the screen prints 22.069, while Molokanov clears 380. Away from that edge, expect roughly 7% error, and rather more in systems whose relative volatility swings widely down the column. The N returned counts the reboiler as an equilibrium stage; a total condenser counts for nothing. Round up, then divide by overall tray efficiency.
Frequently asked questions
The defaults give N = 9.885. How many trays do I actually buy?
Why does the calculator demand R strictly greater than R_min?
Where do R_min and N_min themselves come from?
Related Tools
Optimum Reflux Ratio
Compute a distillation column's operating reflux ratio, R = factor · Rmin, multiplying the minimum reflux by a factor (typically 1.1 to 1.5). There is a classic economic trade-off: a low reflux (near minimum) requires many plates (more column investment); a high reflux requires fewer plates but much more energy in the reboiler and condenser (more operating cost). The optimum balances the two. Enter the minimum reflux and the factor.
Overall Plate Efficiency (O'Connell Correlation)
Estimates the overall plate efficiency of a distillation column by the O'Connell correlation, the step that converts theoretical stages into actual trays in preliminary design. The correlation uses a single combined variable, the product of liquid viscosity and mean relative volatility, and reads efficiency = 0.492 × (viscosity × relative volatility) raised to −0.245, with viscosity in centipoise. What governs it is the product μ·α: a low-viscosity liquid with relative volatility near 1 lands in the 70 to 80% band, and efficiency falls to 40% or less once the liquid is viscous or the relative volatility is high. Note that a difficult separation, with α close to 1, gives the highest efficiency per tray — that column is tall because of the theoretical stage count, not the efficiency, which there actually dampens the height. The analytical fit 0.492 × (μ·α)^−0.245 was adopted, the usual form for hand calculation; the Kessler and Wankat polynomial on the logarithm of the product differs by a few percentage points at the ends of the range. Enter the liquid viscosity at the mean column temperature and the mean relative volatility.
Actual Plates (Efficiency)
Compute the number of actual plates of a distillation column, N_actual = N_theoretical / (efficiency/100), from the number of theoretical (equilibrium) plates and the column's overall efficiency (%). Since no real plate reaches perfect equilibrium, more actual plates than theoretical are needed: a 50% efficiency doubles the plate count. It is the step that turns the theoretical design into the physical column. Enter the theoretical plates and the overall efficiency.
Minimum Reflux (Underwood)
Estimate the minimum reflux ratio of a binary distillation by Underwood's equation (saturated-liquid feed), Rmin = [xD/xF − α·(1−xD)/(1−xF)]/(α − 1), from the relative volatility (α) and the light-component mole fractions in the distillate (xD) and feed (xF). At minimum reflux, the column would need infinite plates; the operating reflux is a multiple of it (1.1–1.5×). It is a key number in column design. Enter α, xD and xF.
Absolute Permeability by the Timur Correlation
Estimates the absolute permeability of a sandstone from porosity and irreducible water saturation, when no core sample is available for laboratory testing. The Timur correlation (1968), fitted on 155 Alaskan sandstone samples, is k = 0.136 × porosity^4.4 ÷ irreducible saturation², with permeability in millidarcy. The result tells how well the reservoir can flow: below 1 mD the formation is considered tight, between 10 and 100 mD it is moderate and above 500 mD it is excellent. The convention that causes most errors here is the unit: the constants 0.136 and 4.4 were fitted with porosity and saturation in percent, not as fractions — entering 0.22 instead of 22 drops the result by several orders of magnitude, which is why both fields ask for percent. Enter the effective porosity and the irreducible water saturation.
Total Reduction Ratio
Calculate the total reduction ratio of a three-stage comminution circuit, RR_total = RR₁ × RR₂ × RR₃, multiplying the reduction ratios of each crusher/mill in series. The dimensionless result is the circuit's overall size reduction — from bench rock blocks (hundreds of mm) to fine particles (mm or µm). Since each stage has a limited reduction ratio (4 to 10 for crushers), large total reductions (100, 1000 or more) require several stages in series: primary, secondary, tertiary crushing and milling. Enter the reduction ratios of the three stages.
The results provided by this tool are for general informational and educational purposes only and do not constitute professional, financial, medical, legal, tax or accounting advice. Always confirm important decisions with a qualified professional and official sources.