Distillation Column Calculator
Estimate the number of theoretical stages using McCabe-Thiele method for binary distillation. Enter values for instant results with step-by-step formulas.
Reviewed by Daniel Agrici, Founder & Lead Developer
Formula
Nmin = ln[(xD(1-xB))/((1-xD)xB)] / ln(alpha)
The Fenske equation gives minimum stages at total reflux. Combined with the Underwood equation for minimum reflux ratio and the Gilliland correlation, the actual number of theoretical stages is estimated. Actual trays = theoretical stages / tray efficiency.
Worked Examples
Example 1: Benzene-Toluene Binary Separation
Problem:Separate a 50/50 benzene-toluene feed (1000 mol/hr) to 95% benzene distillate and 5% benzene bottoms. Relative volatility = 2.5, reflux ratio = 1.5, saturated liquid feed, 70% tray efficiency.
Solution:Material balance: D = 1000 x (0.50 - 0.05)/(0.95 - 0.05) = 500 mol/hr, B = 500 mol/hr\n\nFenske (minimum stages): Nmin = ln[(0.95/0.05)(0.95/0.05)] / ln(2.5)\n= ln(361) / 0.916 = 5.89 / 0.916 = 6.4 stages\n\nUnderwood (minimum reflux): Rmin = (1/1.5)(0.95/0.50 - 2.5 x 0.05/0.50)\n= 0.667 x (1.90 - 0.25) = 1.10\n\nGilliland: X = (1.5 - 1.10)/(1.5 + 1) = 0.16\nY = 1 - exp[...] = 0.528\nN = (0.528 + 6.4 x 0.528 + 6.4)/(1 - 0.528) = 10.5 / 0.472 = 22.2\n\nActual trays = 22.2 / 0.70 = 32 trays
Result:D = 500 mol/hr | Nmin = 6.4 | Rmin = 1.10 | N = 22 theoretical | 32 actual trays
Example 2: High Purity Methanol-Water Separation
Problem:Separate methanol-water (40% methanol feed, 500 mol/hr) to 99% methanol distillate and 1% bottoms. Alpha = 3.0, R = 2.0, q = 1.0, efficiency = 65%.
Solution:D = 500 x (0.40 - 0.01)/(0.99 - 0.01) = 199.0 mol/hr\nB = 500 - 199 = 301.0 mol/hr\n\nNmin = ln[(0.99/0.01)(0.99/0.01)] / ln(3.0)\n= ln(9801) / 1.099 = 9.19 / 1.099 = 8.4 stages\n\nRmin = (1/2.0)(0.99/0.40 - 3.0 x 0.01/0.60)\n= 0.50 x (2.475 - 0.05) = 1.213\n\nX = (2.0 - 1.213)/(2.0 + 1) = 0.262\nN = apply Gilliland: approx 17.8 theoretical\nActual = 17.8 / 0.65 = 28 trays
Result:D = 199 mol/hr | Nmin = 8.4 | N = 17.8 theoretical | 28 actual trays
Frequently Asked Questions
What is distillation and how does a distillation column work?
Distillation is a separation process that exploits differences in the boiling points (volatilities) of components in a liquid mixture. A distillation column provides multiple stages of contact between a rising vapor and a descending liquid, with each stage enriching the vapor in the more volatile (lighter) component. Feed enters the column at an intermediate point, dividing it into a rectifying section (above the feed) and a stripping section (below the feed). Vapor rises through the column, contacting liquid on each tray where mass transfer occurs. At the top, vapor is condensed and part is returned as reflux while the rest is withdrawn as distillate product. At the bottom, liquid is partially vaporized in the reboiler and the remaining liquid is withdrawn as bottoms product. Each theoretical stage achieves an equilibrium separation, and the number of stages determines the degree of separation achievable between the components.
What is the McCabe-Thiele method for binary distillation design?
The McCabe-Thiele method is a graphical technique for determining the number of theoretical stages required for binary distillation. It uses operating lines and the equilibrium curve on an x-y diagram (where x and y are liquid and vapor mole fractions of the light component). The rectifying section operating line runs from (xD, xD) with slope R/(R+1), where R is the reflux ratio. The stripping section operating line runs from (xB, xB) to the intersection with the rectifying line at the feed point. Stages are stepped off between the operating lines and the equilibrium curve, starting from the distillate composition. Each step represents one theoretical stage. The equilibrium curve for an ideal binary system is y = alpha*x / (1 + (alpha-1)*x), where alpha is the relative volatility. The method assumes constant molar overflow (equal molar flow rates throughout each section) and is remarkably accurate for preliminary design despite its simplicity.
How does the feed condition (q value) affect distillation column performance?
The feed condition parameter q represents the thermal state of the feed and significantly affects the column internal flows and the number of stages required. A q value of 1.0 means saturated liquid feed (all liquid at its bubble point), 0.0 means saturated vapor (all vapor at its dew point), values between 0 and 1 represent partially vaporized feed, q > 1 means subcooled liquid, and q < 0 means superheated vapor. The q value determines the slope of the feed line (q-line) on the McCabe-Thiele diagram: slope = q/(q-1). A saturated liquid feed (q=1) gives a vertical feed line, while saturated vapor (q=0) gives a horizontal line. Subcooled liquid feeds increase the internal liquid flow in the stripping section, potentially improving separation but requiring more reboiler duty. Vaporized feeds reduce internal liquid flow. The optimal feed condition minimizes total energy input and usually corresponds to a partially preheated feed that matches the column internal temperature at the feed tray.
How do you determine the optimal feed stage location in a distillation column?
The optimal feed stage location minimizes the total number of stages required for the desired separation. In the McCabe-Thiele method, the feed stage is where the operating line transitions from the rectifying section to the stripping section, which occurs at the intersection of the two operating lines. The Kirkbride equation provides an analytical estimate: log(Nr/Ns) = 0.206 x log[(xF/(1-xF)) x ((1-xB)/xD)^2 x (B/D)], where Nr and Ns are the number of rectifying and stripping stages, B and D are bottoms and distillate flow rates. Feeding too high (above optimal) wastes stages in the stripping section by not utilizing the rectifying section fully. Feeding too low wastes rectifying stages. In practice, columns often include multiple feed nozzles at different elevations to allow optimization during operation. Rigorous simulation using software like Aspen Plus, ChemCAD, or HYSYS determines the exact optimal feed location by minimizing reboiler duty or total annual cost.
References
Reviewed by Daniel Agrici, Founder & Lead Developer · Editorial policy