Batch Distillation Product Collection Strategy

Reference ID: MET-4D87 | Process Engineering Reference Sheets Calculation Guide

Introduction & Context

The Batch Distillation Product Collection Strategy is a fundamental process engineering model used to determine optimal cut points during the separation of multi-component liquid mixtures. In industries such as beverage production (e.g., whisky distillation) or fine chemical synthesis, the quality of the final product is highly dependent on the precise separation of volatile components—categorized as heads, hearts, and tails—based on their respective vapor pressures and concentrations.

This calculation is critical for process control, as it allows engineers to predict the transition times between collection fractions by modeling the convective mass transfer from the liquid-vapor interface. By applying mass balance and convective flux equations, operators can automate the diversion of distillate streams to ensure product consistency and purity.

Methodology & Formulas

The model assumes that vapor-phase resistance controls the mass transfer rate and that the bulk vapor concentration is negligible (i.e., perfect sweep or rapid removal of evolved vapor). The system dynamics are governed by the following physical relationships:

The convective mass flux for each component i is defined by the mass transfer coefficient, the surface area, and the interface concentration:

\[ \dot{m}_{i} = h_{m} \cdot A \cdot C_{s,i} \]

The interface concentration is derived from the liquid mole fraction and the vapor pressure of the component at the boiling temperature, using Raoult's law and the ideal gas law (with vapor pressure \(P_{\mathrm{vap},i}\) expressed in Pa):

\[ C_{s,i} = \frac{M_{i} \cdot P_{\mathrm{vap},i}}{R \cdot T} \cdot x_{i} \]

The liquid mole fraction is calculated based on the instantaneous mass of each component in the pot:

\[ x_{i} = \frac{m_{i} / M_{i}}{\sum_{j} (m_{j} / M_{j})} \]

The mass transfer coefficient is determined by the Sherwood number, which accounts for the flow regime at the interface:

\[ h_{m} = \frac{Sh \cdot D_{AB}}{L} \]

The instantaneous vapor composition, which dictates the cut point logic, is the ratio of the individual component flux to the total mass flux:

\[ y_{i} = \frac{\dot{m}_{i}}{\sum_{j} \dot{m}_{j}} \]

Regime	Condition	Sherwood Number Correlation
Laminar	\( Re < 5 \cdot 10^{5} \)	\( Sh = 0.664 \cdot Re^{0.5} \cdot Sc^{1/3} \)
Turbulent	\( Re \geq 5 \cdot 10^{5} \)	\( Sh = 0.037 \cdot Re^{0.8} \cdot Sc^{1/3} \)

The cut points are controlled by monitoring the vapor-phase mass-fraction ratios against defined threshold values. Recommended thresholds for a typical potable-spirits separation are:

Cut Point	Logic Criteria	Typical Threshold
Heads Cut End	\( \displaystyle \frac{y_{\mathrm{heads}}}{y_{\mathrm{hearts}}} < \varepsilon_{H} \)	\( \varepsilon_{H} = 0.05 \)
Tails Cut Start	\( \displaystyle \frac{y_{\mathrm{tails}}}{y_{\mathrm{hearts}}} > \varepsilon_{T} \)	\( \varepsilon_{T} = 0.02 \)

Note: The cut criteria must be evaluated continuously throughout the batch simulation, not only at the final time step. The cumulative distillate volume corresponding to the first time step at which the heads ratio falls below \(\varepsilon_{H}\) is recorded as \(V_{\mathrm{heads}}\); similarly, the tails cut volume \(V_{\mathrm{hearts}}\) is recorded when the tails ratio first exceeds \(\varepsilon_{T}\).

To ensure product purity and minimize waste, process engineers should monitor the following indicators, which directly relate to the vapor-phase mass transfer model:

Real-time refractive index measurements of the distillate stream, correlated to the calculated vapor mass fraction \(y_{i}\).
Head temperature stability relative to the vapor-liquid equilibrium curve for the target component.
Gas chromatography analysis of intermediate samples taken at the product condenser outlet to validate the modeled component fluxes \(\dot{m}_{i}\).
Calculated mass balance thresholds based on the initial charge composition; compare the instantaneous ratio \(y_{\mathrm{heads}} / y_{\mathrm{hearts}}\) against the defined threshold \(\varepsilon_{H}\) (typically 0.05) and \(y_{\mathrm{tails}} / y_{\mathrm{hearts}}\) against \(\varepsilon_{T}\) (typically 0.02).

Managing the transition requires a controlled adjustment of the reflux ratio to sharpen the separation. The modeled vapor composition \(y_{i}\) guides this decision. Follow these steps:

Increase the reflux ratio as the head temperature approaches the target boiling point of the hearts component.
Continuously evaluate \(y_{\mathrm{heads}} / y_{\mathrm{hearts}}\) from the mass transfer model; isolate the fore-run receiver once this ratio drops below the specification limit \(\varepsilon_{H}\).
Verify the purity of the initial distillate volume before diverting flow to the main product tank.
Maintain a steady boil-up rate to prevent pressure surges during the valve switch, ensuring the vapor velocity \(V\) remains stable and the Sherwood number correlation stays within the laminar regime.

To maximize recovery of the hearts fraction while preventing tails contamination, consider implementing these strategies, informed by the mass transfer model:

Monitor \(y_{\mathrm{tails}} / y_{\mathrm{hearts}}\) continuously; terminate hearts collection as soon as this ratio exceeds \(\varepsilon_{T}\) to avoid off-flavors in the product.
Utilize a vacuum system to lower the boiling point of the remaining heavy components, which reduces their vapor pressure \(P_{\mathrm{vap,T}}\) and delays the onset of the tails cut.
Perform a final bottoms analysis to determine if further stripping is economically viable based on the residual component masses \(m_{i}\).
Implement a nitrogen sparge to assist in the removal of residual product from the reboiler by increasing the effective mass transfer coefficient \(h_{m}\).
Ensure the condenser duty is sufficient to capture low-boiling fractions during the final stages of the batch, maintaining the assumption of negligible bulk vapor concentration.

Batch Distillation Product Collection: Worked Example

Scenario: A 1000 L whisky wash is charged to a batch pot distillation column. The objective is to determine the cumulative distillate volumes at which the heads and hearts cuts end, using the vapor-phase mass transfer model. The system operates at constant temperature with the vapor bulk concentration assumed zero (perfect sweep).

Knowns:

Operating temperature, \(T = 363.0\) K
Vapor-liquid interfacial area, \(A = 1.5\) m²
Vapor velocity, \(V = 5.0\) m/s
Characteristic length, \(L = 1.0\) m
Schmidt number, \(Sc = 2.0\)
Diffusivity, \(D_{AB} = 1.0 \times 10^{-5}\) m²/s
Vapor density, \(\rho_{\mathrm{vap}} = 1.2\) kg/m³
Vapor dynamic viscosity, \(\mu_{\mathrm{vap}} = 1.8 \times 10^{-5}\) Pa·s
Distillate density, \(\rho_{\mathrm{dist}} = 800.0\) kg/m³
Universal gas constant, \(R = 8314.0\) J/(kmol·K)
Time step, \(\Delta t = 1.0\) s
Heads threshold, \(\varepsilon_{H} = 0.05\)
Tails threshold, \(\varepsilon_{T} = 0.02\)
Congeners (Heads/Hearts/Tails):
Heads (Ethyl Acetate proxy): Molar mass \(M_{\mathrm{H}} = 88.1\) kg/kmol, vapor pressure \(P_{\mathrm{vap,H}} = 200.0\) kPa, initial mass \(m_{0,\mathrm{H}} = 20.0\) kg
Hearts (Ethanol): \(M_{\mathrm{E}} = 46.0\) kg/kmol, \(P_{\mathrm{vap,E}} = 160.0\) kPa, \(m_{0,\mathrm{E}} = 100.0\) kg
Tails (Isoamyl Alcohol proxy): \(M_{\mathrm{T}} = 88.1\) kg/kmol, \(P_{\mathrm{vap,T}} = 30.0\) kPa, \(m_{0,\mathrm{T}} = 1.0\) kg

Step-by-Step Calculation:

Compute Reynolds number: \[ Re = \frac{\rho_{\mathrm{vap}} \cdot V \cdot L}{\mu_{\mathrm{vap}}} = \frac{1.2 \cdot 5.0 \cdot 1.0}{1.8 \times 10^{-5}} = 3.333 \times 10^{5} \] Since \(Re = 3.333 \times 10^{5} < 5 \times 10^{5}\), the flow is laminar.
Compute Sherwood number using laminar correlation: \[ Sh = 0.664 \cdot Re^{1/2} \cdot Sc^{1/3} = 0.664 \cdot (3.333 \times 10^{5})^{0.5} \cdot (2.0)^{1/3} = 483.0 \]
Compute mass transfer coefficient: \[ h_{m} = \frac{Sh \cdot D_{AB}}{L} = \frac{483.0 \cdot (1.0 \times 10^{-5})}{1.0} = 4.83 \times 10^{-3} \text{ m/s} \]
Initial liquid mole fractions (time \(t = 0\) s):
First, compute the initial kilomoles of each component: \[ n_{\mathrm{H}}(0) = \frac{20.0}{88.1} = 0.2270 \text{ kmol},\quad n_{\mathrm{E}}(0) = \frac{100.0}{46.0} = 2.1739 \text{ kmol},\quad n_{\mathrm{T}}(0) = \frac{1.0}{88.1} = 0.01135 \text{ kmol} \] Total kilomoles: \(n_{\mathrm{tot}} = 2.4123\) kmol. \[ x_{\mathrm{H}}(0) = \frac{0.2270}{2.4123} = 0.0941,\quad x_{\mathrm{E}}(0) = \frac{2.1739}{2.4123} = 0.9012,\quad x_{\mathrm{T}}(0) = \frac{0.01135}{2.4123} = 0.0047 \] (Verification: \(0.0941 + 0.9012 + 0.0047 = 1.0000\) ✓)
Simulate mass transfer over 3600 s: For each time step, compute the saturation concentration at the interface using Raoult's law and the ideal gas law (with \(P_{\mathrm{vap}}\) converted from kPa to Pa): \[ C_{s,i}(t) = \frac{M_{i} \cdot (P_{\mathrm{vap},i} \cdot 1000)}{R \cdot T} \cdot x_{i}(t) \] Then determine the evaporation flux: \[ \dot{m}_{i}(t) = h_{m} \cdot A \cdot C_{s,i}(t) \] Update liquid masses (forward Euler integration): \[ m_{i}(t+\Delta t) = m_{i}(t) - \dot{m}_{i}(t) \cdot \Delta t \] At each time step, monitor the vapor mass-fraction ratios to detect cut-point crossings. After 3600 s of simulation with \(\Delta t = 1.0\) s, the remaining masses are: \[ m_{\mathrm{H}} \approx 7.23 \text{ kg},\quad m_{\mathrm{E}} \approx 44.3 \text{ kg},\quad m_{\mathrm{T}} \approx 0.86 \text{ kg} \]
Compute cumulative distillate volume: Total mass evaporated: \[ \Delta m_{\mathrm{total}} = (20.0 - 7.23) + (100.0 - 44.3) + (1.0 - 0.86) = 12.77 + 55.7 + 0.14 = 68.6 \text{ kg} \] \[ V_{\mathrm{total}} = \frac{\Delta m_{\mathrm{total}}}{\rho_{\mathrm{dist}}} = \frac{68.6}{800.0} \approx 0.086 \text{ m}^{3} \]
Evaluate cut point thresholds at the final time step:
Vapor mole fractions (mass-flux fractions) at \(t = 3600\) s: \[ y_{\mathrm{H}} \approx 0.169,\quad y_{\mathrm{E}} \approx 0.828,\quad y_{\mathrm{T}} \approx 0.003 \] Heads-to-hearts ratio: \[ \frac{y_{\mathrm{H}}}{y_{\mathrm{E}}} \approx \frac{0.169}{0.828} \approx 0.204 > \varepsilon_{H}\;(0.05) \] Tails-to-hearts ratio: \[ \frac{y_{\mathrm{T}}}{y_{\mathrm{E}}} \approx \frac{0.003}{0.828} \approx 0.0036 < \varepsilon_{T}\;(0.02) \]
Determine cut status:
Throughout the 3600 s simulation, the heads-to-hearts ratio \(y_{\mathrm{H}}/y_{\mathrm{E}}\) started at approximately 0.25 and decreased monotonically, but never fell below the threshold \(\varepsilon_{H} = 0.05\). Therefore, the heads cut has not yet ended; all distillate collected to this point is heads product (\(V_{\mathrm{heads}} = 0.086\) m³).

The tails-to-hearts ratio \(y_{\mathrm{T}}/y_{\mathrm{E}}\) remained below \(\varepsilon_{T} = 0.02\) throughout the simulation, indicating that the tails cut has not yet started (\(V_{\mathrm{hearts}} = 0.0\) m³). The simulation must be extended further to capture both the end of the heads cut and the onset of the tails cut.

Final Answer: After 3600 s of distillation, the cumulative volume collected is 0.086 m³. The heads cut is still active because the ratio of heads to hearts in the vapor (\(\approx 0.204\)) exceeds the threshold of 0.05. The heads cut volume collected to date is 0.086 m³; the hearts cut volume is 0.0 m³ (not yet started). To complete the separation, the simulation should be extended until \(y_{\mathrm{H}}/y_{\mathrm{E}} < 0.05\) at which point the cumulative volume at that instant defines \(V_{\mathrm{heads}}\), and collection of hearts can commence.

"Un projet n'est jamais trop grand s'il est bien conçu." — André Citroën

"La difficulté attire l'homme de caractère, car c'est en l'étreignant qu'il se réalise." — Charles de Gaulle