# Shortcut Distillation Column Calculator Documentation

# 1. Objective

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}) and molar basis (kmol/h\text{kmol/h}).


# 2. Design Basis & Mass-Molar Conversion Methodology

# 2.1 Feed Mass Flow Rate & Composition Conversion

Plant distillation columns are typically operated on a mass flow basis (FmassF_{\text{mass}} in kg/h\text{kg/h}).

Given:

  • Feed mass flow rate: FmassF_{\text{mass}} (kg/h\text{kg/h})
  • Light Key component Molecular Weight: MLM_L (g/mol\text{g/mol} or kg/kmol\text{kg/kmol})
  • Heavy Key component Molecular Weight: MHM_H (g/mol\text{g/mol} or kg/kmol\text{kg/kmol})

# Converting Mass Fractions (ww) to Mole Fractions (z,xz, x):

zF=wFMLwFML+1−wFMHz_F = \frac{\frac{w_F}{M_L}}{\frac{w_F}{M_L} + \frac{1 - w_F}{M_H}}
xD=wDMLwDML+1−wDMHx_D = \frac{\frac{w_D}{M_L}}{\frac{w_D}{M_L} + \frac{1 - w_D}{M_H}}
xB=wBMLwBML+1−wBMHx_B = \frac{\frac{w_B}{M_L}}{\frac{w_B}{M_L} + \frac{1 - w_B}{M_H}}

# Average Stream Molecular Weights:

MF,avg=zFML+(1−zF)MH[kg/kmol]M_{F, \text{avg}} = z_F M_L + (1 - z_F) M_H \quad [\text{kg/kmol}]
MD,avg=xDML+(1−xD)MH[kg/kmol]M_{D, \text{avg}} = x_D M_L + (1 - x_D) M_H \quad [\text{kg/kmol}]
MB,avg=xBML+(1−xB)MH[kg/kmol]M_{B, \text{avg}} = x_B M_L + (1 - x_B) M_H \quad [\text{kg/kmol}]

# Molar & Mass Overall Material Balances:

Fmolar=FmassMF,avg[kmol/h]F_{\text{molar}} = \frac{F_{\text{mass}}}{M_{F, \text{avg}}} \quad [\text{kmol/h}]
Dmolar=Fmolar×(zF−xBxD−xB)[kmol/h]D_{\text{molar}} = F_{\text{molar}} \times \left(\frac{z_F - x_B}{x_D - x_B}\right) \quad [\text{kmol/h}]
Bmolar=Fmolar−Dmolar[kmol/h]B_{\text{molar}} = F_{\text{molar}} - D_{\text{molar}} \quad [\text{kmol/h}]
Dmass=Dmolar×MD,avg[kg/h]D_{\text{mass}} = D_{\text{molar}} \times M_{D, \text{avg}} \quad [\text{kg/h}]
Bmass=Bmolar×MB,avg[kg/h]B_{\text{mass}} = B_{\text{molar}} \times M_{B, \text{avg}} \quad [\text{kg/h}]

# 3. FUG Shortcut Method Equations

  1. Fenske Minimum Stages (NminN_{\text{min}} at Total Reflux):
Nmin=ln⁡[(xD1−xD)×(1−xBxB)]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)}
  1. Underwood Minimum Reflux Ratio (RminR_{\text{min}}):
    For saturated liquid feed (q=1.0q = 1.0):
Rmin=1α−1[xDzF−α(1−xD)1−zF]R_{\text{min}} = \frac{1}{\alpha - 1} \left[ \frac{x_D}{z_F} - \frac{\alpha (1 - x_D)}{1 - z_F} \right]
  1. Operating Reflux Ratio (RactualR_{\text{actual}}):
Ractual=Rmin×fR(typically fR=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)
  1. Gilliland Stage Correlation (NtheoreticalN_{\text{theoretical}}):
XG=Ractual−RminRactual+1X_G = \frac{R_{\text{actual}} - R_{\text{min}}}{R_{\text{actual}} + 1}
YG=0.75×(1−XG0.5668)Y_G = 0.75 \times (1 - X_G^{0.5668})
Ntheoretical=Nmin+YG1−YGN_{\text{theoretical}} = \frac{N_{\text{min}} + Y_G}{1 - Y_G}
  1. Actual Physical Trays:
Nactual=⌈Ntheoreticalηstage/100⌉N_{\text{actual}} = \left\lceil \frac{N_{\text{theoretical}}}{\eta_{\text{stage}} / 100} \right\rceil
  1. Kirkbride Optimum Feed Stage (NfeedN_{\text{feed}} from top):
log⁡(NrectNstrip)=0.206log⁡[(BmolarDmolar)(1−zFzF)(xB1−xD)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]

# 4. Column Hydraulic Sizing (Souders-Brown)

  1. Overhead Vapor Flow:
Vtop, molar=Dmolar×(Ractual+1)[kmol/h]V_{\text{top, molar}} = D_{\text{molar}} \times (R_{\text{actual}} + 1) \quad [\text{kmol/h}]
Vtop, mass=Vtop, molar×MD,avg[kg/h]V_{\text{top, mass}} = V_{\text{top, molar}} \times M_{D, \text{avg}} \quad [\text{kg/h}]
  1. Souders-Brown Maximum Flooding Velocity (ufloodu_{\text{flood}}):
uflood=CSBρL−ρVρV[m/s]u_{\text{flood}} = C_{\text{SB}} \sqrt{\frac{\rho_L - \rho_V}{\rho_V}} \quad [\text{m/s}]
udesign=0.75×ufloodu_{\text{design}} = 0.75 \times u_{\text{flood}}
  1. Column Diameter (DcolD_{\text{col}}):
Acol=Qvaporudesign  ⟹  Dcol=4Acolπ[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}]

# 5. Industrial Binary Solvent Systems Included

  • Methanol / Water (α=2.80\alpha = 2.80)
  • Ethanol / Water (α=2.20\alpha = 2.20)
  • Acetone / Water (α=3.50\alpha = 3.50)
  • Isopropanol (IPA) / Water (α=2.10\alpha = 2.10)
  • Benzene / Toluene (α=2.40\alpha = 2.40)
  • Ethyl Acetate / Ethanol (α=1.85\alpha = 1.85)