The Shortcut Distillation Column Calculator solves preliminary design and rating problems for binary distillation columns. It utilizes the classical Fenske-Underwood-Gilliland (FUG) shortcut framework combined with Kirkbride's feed tray location formula and Souders-Brown vapor capacity hydraulics to size columns on a mass basis (kg/h \text{kg/h} kg/h ) and molar basis (kmol/h \text{kmol/h} kmol/h ).
Plant distillation columns are typically operated on a mass flow basis (F mass F_{\text{mass}} F mass in kg/h \text{kg/h} kg/h ).
Given:
Feed mass flow rate: F mass F_{\text{mass}} F mass (kg/h \text{kg/h} kg/h )
Light Key component Molecular Weight: M L M_L M L (g/mol \text{g/mol} g/mol or kg/kmol \text{kg/kmol} kg/kmol )
Heavy Key component Molecular Weight: M H M_H M H (g/mol \text{g/mol} g/mol or kg/kmol \text{kg/kmol} kg/kmol )
z F = w F M L w F M L + 1 − w F M H z_F = \frac{\frac{w_F}{M_L}}{\frac{w_F}{M_L} + \frac{1 - w_F}{M_H}} z F = M L w F + M H 1 − w F M L w F
x D = w D M L w D M L + 1 − w D M H x_D = \frac{\frac{w_D}{M_L}}{\frac{w_D}{M_L} + \frac{1 - w_D}{M_H}} x D = M L w D + M H 1 − w D M L w D
x B = w B M L w B M L + 1 − w B M H x_B = \frac{\frac{w_B}{M_L}}{\frac{w_B}{M_L} + \frac{1 - w_B}{M_H}} x B = M L w B + M H 1 − w B M L w B
M F , avg = z F M L + ( 1 − z F ) M H [ kg/kmol ] M_{F, \text{avg}} = z_F M_L + (1 - z_F) M_H \quad [\text{kg/kmol}] M F , avg = z F M L + ( 1 − z F ) M H [ kg/kmol ]
M D , avg = x D M L + ( 1 − x D ) M H [ kg/kmol ] M_{D, \text{avg}} = x_D M_L + (1 - x_D) M_H \quad [\text{kg/kmol}] M D , avg = x D M L + ( 1 − x D ) M H [ kg/kmol ]
M B , avg = x B M L + ( 1 − x B ) M H [ kg/kmol ] M_{B, \text{avg}} = x_B M_L + (1 - x_B) M_H \quad [\text{kg/kmol}] M B , avg = x B M L + ( 1 − x B ) M H [ kg/kmol ]
F molar = F mass M F , avg [ kmol/h ] F_{\text{molar}} = \frac{F_{\text{mass}}}{M_{F, \text{avg}}} \quad [\text{kmol/h}] F molar = M F , avg F mass [ kmol/h ]
D molar = F molar × ( z F − x B x D − x B ) [ kmol/h ] D_{\text{molar}} = F_{\text{molar}} \times \left(\frac{z_F - x_B}{x_D - x_B}\right) \quad [\text{kmol/h}] D molar = F molar × ( x D − x B z F − x B ) [ kmol/h ]
B molar = F molar − D molar [ kmol/h ] B_{\text{molar}} = F_{\text{molar}} - D_{\text{molar}} \quad [\text{kmol/h}] B molar = F molar − D molar [ kmol/h ]
D mass = D molar × M D , avg [ kg/h ] D_{\text{mass}} = D_{\text{molar}} \times M_{D, \text{avg}} \quad [\text{kg/h}] D mass = D molar × M D , avg [ kg/h ]
B mass = B molar × M B , avg [ kg/h ] B_{\text{mass}} = B_{\text{molar}} \times M_{B, \text{avg}} \quad [\text{kg/h}] B mass = B molar × M B , avg [ kg/h ]
Fenske Minimum Stages (N min N_{\text{min}} N min at Total Reflux) :
N min = ln [ ( x D 1 − x D ) × ( 1 − x B x B ) ] ln ( α ) N_{\text{min}} = \frac{\ln\left[ \left(\frac{x_D}{1 - x_D}\right) \times \left(\frac{1 - x_B}{x_B}\right) \right]}{\ln(\alpha)} N min = ln ( α ) ln [ ( 1 − x D x D ) × ( x B 1 − x B ) ]
Underwood Minimum Reflux Ratio (R min R_{\text{min}} R min ) : For saturated liquid feed (q = 1.0 q = 1.0 q = 1.0 ):
R min = 1 α − 1 [ x D z F − α ( 1 − x D ) 1 − z F ] R_{\text{min}} = \frac{1}{\alpha - 1} \left[ \frac{x_D}{z_F} - \frac{\alpha (1 - x_D)}{1 - z_F} \right] R min = α − 1 1 [ z F x D − 1 − z F α ( 1 − x D ) ]
Operating Reflux Ratio (R actual R_{\text{actual}} R actual ) :
R actual = R min × f R ( typically f R = 1.20 to 1.50 ) R_{\text{actual}} = R_{\text{min}} \times f_R \quad (\text{typically } f_R = 1.20 \text{ to } 1.50) R actual = R min × f R ( typically f R = 1.20 to 1.50 )
Gilliland Stage Correlation (N theoretical N_{\text{theoretical}} N theoretical ) :
X G = R actual − R min R actual + 1 X_G = \frac{R_{\text{actual}} - R_{\text{min}}}{R_{\text{actual}} + 1} X G = R actual + 1 R actual − R min
Y G = 0.75 × ( 1 − X G 0.5668 ) Y_G = 0.75 \times (1 - X_G^{0.5668}) Y G = 0.75 × ( 1 − X G 0.5668 )
N theoretical = N min + Y G 1 − Y G N_{\text{theoretical}} = \frac{N_{\text{min}} + Y_G}{1 - Y_G} N theoretical = 1 − Y G N min + Y G
Actual Physical Trays :
N actual = ⌈ N theoretical η stage / 100 ⌉ N_{\text{actual}} = \left\lceil \frac{N_{\text{theoretical}}}{\eta_{\text{stage}} / 100} \right\rceil N actual = ⌈ η stage /100 N theoretical ⌉
Kirkbride Optimum Feed Stage (N feed N_{\text{feed}} N feed from top) :
log ( N rect N strip ) = 0.206 log [ ( B molar D molar ) ( 1 − z F z F ) ( x B 1 − x D ) 2 ] \log\left(\frac{N_{\text{rect}}}{N_{\text{strip}}}\right) = 0.206 \log\left[ \left(\frac{B_{\text{molar}}}{D_{\text{molar}}}\right) \left(\frac{1 - z_F}{z_F}\right) \left(\frac{x_B}{1 - x_D}\right)^2 \right] log ( N strip N rect ) = 0.206 log [ ( D molar B molar ) ( z F 1 − z F ) ( 1 − x D x B ) 2 ]
Overhead Vapor Flow :
V top, molar = D molar × ( R actual + 1 ) [ kmol/h ] V_{\text{top, molar}} = D_{\text{molar}} \times (R_{\text{actual}} + 1) \quad [\text{kmol/h}] V top, molar = D molar × ( R actual + 1 ) [ kmol/h ]
V top, mass = V top, molar × M D , avg [ kg/h ] V_{\text{top, mass}} = V_{\text{top, molar}} \times M_{D, \text{avg}} \quad [\text{kg/h}] V top, mass = V top, molar × M D , avg [ kg/h ]
Souders-Brown Maximum Flooding Velocity (u flood u_{\text{flood}} u flood ) :
u flood = C SB ρ L − ρ V ρ V [ m/s ] u_{\text{flood}} = C_{\text{SB}} \sqrt{\frac{\rho_L - \rho_V}{\rho_V}} \quad [\text{m/s}] u flood = C SB ρ V ρ L − ρ V [ m/s ]
u design = 0.75 × u flood u_{\text{design}} = 0.75 \times u_{\text{flood}} u design = 0.75 × u flood
Column Diameter (D col D_{\text{col}} D col ) :
A col = Q vapor u design ⟹ D col = 4 A col π [ meters ⟹ mm ] A_{\text{col}} = \frac{Q_{\text{vapor}}}{u_{\text{design}}} \quad \implies \quad D_{\text{col}} = \sqrt{\frac{4 A_{\text{col}}}{\pi}} \quad [\text{meters} \implies \text{mm}] A col = u design Q vapor ⟹ D col = π 4 A col [ meters ⟹ mm ]
Methanol / Water (α = 2.80 \alpha = 2.80 α = 2.80 )
Ethanol / Water (α = 2.20 \alpha = 2.20 α = 2.20 )
Acetone / Water (α = 3.50 \alpha = 3.50 α = 3.50 )
Isopropanol (IPA) / Water (α = 2.10 \alpha = 2.10 α = 2.10 )
Benzene / Toluene (α = 2.40 \alpha = 2.40 α = 2.40 )
Ethyl Acetate / Ethanol (α = 1.85 \alpha = 1.85 α = 1.85 )