Batch Distillation with Reflux

Reference ID: MET-5590 | Process Engineering Reference Sheets Calculation Guide

Introduction & Context

Batch distillation with reflux is a fundamental unit operation in process engineering, characterized by unsteady-state rectification. Unlike continuous distillation, the composition of the pot charge and the distillate changes over time as the more volatile component is depleted. This process is critical in industries such as fine chemicals, pharmaceuticals, and craft spirits production, where high-purity separation is required from a finite batch of raw material.

The calculation of batch distillation performance involves balancing the rectifying operating line against the vapor‑liquid equilibrium (VLE) of the mixture. By applying an effective reflux ratio control strategy, engineers can predict the evolution of distillate purity and the remaining mass in the reboiler, allowing for the precise management of product cuts (heads, hearts, and tails).

Methodology & Formulas

The calculation relies on the McCabe-Thiele graphical method adapted for unsteady-state conditions, integrated with the Rayleigh equation to account for the depletion of the pot charge.

The vapor-liquid equilibrium is modeled using the relative volatility (\(\alpha\)):

\[ y = \frac{\alpha \cdot x}{1 + (\alpha - 1) \cdot x} \]

The rectifying operating line is defined by the reflux ratio (\(R\)) and the distillate composition (\(x_D\)):

\[ y = \left( \frac{R}{R + 1} \right) \cdot x + \left( \frac{x_D}{R + 1} \right) \]

To determine the number of theoretical stages (\(N_{\text{theoretical}}\)), the system iterates through the equilibrium curve and operating line until the pot composition (\(x_W\)) is reached. The total stages include the reboiler:

\[ N_{\text{theoretical}} = N_{\text{trays}} + N_{\text{reboiler}} \]

The actual number of physical trays (\(N_{\text{actual}}\)) is adjusted by the overall plate efficiency (\(\eta\)), applied to the trays only (the reboiler is treated as a perfect equilibrium stage):

\[ N_{\text{actual, trays}} = \frac{N_{\text{trays}}}{\eta} \]

The final mass remaining in the pot (\(W_{\text{final}}\)) is calculated by integrating the Rayleigh equation over the change in pot composition from the initial charge (\(W_{\text{initial}}\)) to the final state:

\[ \ln\left( \frac{W_{\text{initial}}}{W_{\text{final}}} \right) = \int_{x_{W,\text{final}}}^{x_F} \frac{1}{x_D - x_W} \, dx_W \]

Parameter	Constraint/Regime
Reflux Ratio (\(R\))	\(1.5 \leq R \leq 6.0\)
Distillate Purity (\(x_D\))	\(x_D < 0.956\) (Azeotropic Limit)
Composition Gradient	\(x_F < x_D\)
Pot Composition (\(x_W\))	\(0 < x_W < 1\)

Feed composition – initial mass or mole fractions of key components.
Desired product purity – target mass or mole fraction for the overhead distillate (hearts cut).
Batch size – total mass or volume of the pot charge.
Reflux ratio – proportion of condensed overhead returned to the column (typically 1.5–6.0).
Relative volatility – measure of separability between the key components.
Column internals – tray type, spacing, and efficiency (typically 60–85%).
Operating pressure – influences boiling points, relative volatility, and reflux duties.

For a given number of trays, compute the instantaneous distillate composition \(x_D\) as a function of the current pot composition \(x_W\) using McCabe-Thiele stepping down from an assumed \(x_D\).
Perform a Rayleigh integration over the target pot composition range to find the total distillate mass and its average composition.
Repeat the analysis for several candidate reflux ratios (e.g., \(R = 2.0, 3.0, 4.0, 5.0\)).
Plot average product purity versus batch time (or energy consumption) for each ratio.
Select the reflux ratio that meets the purity target with the shortest feasible batch time while keeping energy consumption acceptable.
Validate the chosen ratio on a small-scale trial before full-scale implementation.

Foaming or entrainment – install anti-foam agents and ensure proper tray spacing.
Temperature control drift – use calibrated thermocouples and implement PID control loops.
Reflux flow fluctuations – maintain a stable pump and incorporate flow-rate sensors with feedback control.
Product carry-over – adjust the cut-off point and consider a short-duration purge after the main draw.
Excessive energy use – optimize reflux ratio and consider heat integration with other plant streams.

Specify the target distillate composition \(x_{D,\text{target}}\) and the final pot composition \(x_{W,\text{final}}\) that define the end-of-batch separation requirement.
For a candidate number of trays \(N_{\text{trays}}\), start from the reboiler liquid at \(x_{W,\text{final}}\) and compute the equilibrium vapor.
Step up through the operating line (with the chosen \(R\) and the assumed \(x_D\)) and the VLE curve for \(N_{\text{trays}}\) stages.
Compare the computed top tray vapor composition to the assumed \(x_D\); iterate the assumed \(x_D\) until they match (a single-variable root-finding problem).
If the converged \(x_D\) meets or exceeds \(x_{D,\text{target}}\), the tray count is sufficient; otherwise, increment \(N_{\text{trays}}\) and repeat.
Adjust the theoretical tray count for overall tray efficiency (typically 60–85%) to obtain the physical stage count.

Batch Distillation of Ethanol-Water: Tray Sizing and Rayleigh Integration

A batch column is used to concentrate ethanol from a fermented wash. The pot is initially charged with 1000 kg of wash at 8% w/w ethanol. The column operates at a constant reflux ratio of 4.0. The target distillate composition for the hearts cut is 25% w/w, and the operation continues until the pot composition falls to 4% w/w. The system exhibits a relative volatility \(\alpha = 2.5\) (assumed constant for the range of compositions), and the overall tray efficiency is 60%. The reboiler acts as one equilibrium stage with near-perfect efficiency.

Knowns

Pot charge: \(W_0 = 1000.0\ \text{kg}\)
Initial pot composition: \(x_F = 0.08\) (8% w/w ethanol)
Target distillate composition: \(x_{D,\text{target}} = 0.25\) (25% w/w)
Reflux ratio: \(R = 4.0\) (dimensionless)
Final pot composition: \(x_{W,\text{final}} = 0.04\) (4% w/w)
Relative volatility: \(\alpha = 2.5\)
Overall plate efficiency: \(\eta = 0.60\)

Step-by-Step Calculation

Determine the number of theoretical trays required. Using the iterative tray-counting function with the inputs \(R = 4.0\), \(\alpha = 2.5\), \(x_{W,\text{final}} = 0.04\), and \(x_{D,\text{target}} = 0.25\), the number of theoretical trays (excluding the reboiler) is found to be \(N_{\text{trays}} = 3\). Including the reboiler, the total theoretical stages are: \[ N_{\text{theoretical}} = N_{\text{trays}} + N_{\text{reboiler}} = 3 + 1 = 4. \]
Calculate the actual number of trays. The tray efficiency is applied only to the trays (not the reboiler, which is taken as a perfect equilibrium stage): \[ N_{\text{actual, trays}} = \frac{N_{\text{trays}}}{\eta} = \frac{3}{0.6} = 5.0. \] Including the reboiler as one actual stage gives a total of \(5.0 + 1 = 6.0\) actual stages. The trays are rounded up to the nearest integer: \[ N_{\text{trays, installed}} = 6\ \text{actual trays.} \] (Note: The calculation treats the reboiler separately; applying efficiency to the theoretical stages including reboiler would yield \(4/0.6 = 6.667\) total stages and \(6.667-1 = 5.667\) trays, also rounding to 6 trays.)
Perform the Rayleigh equation integration. The Rayleigh equation for a batch distillation with a variable overhead composition \(x_D\) is: \[ \ln\left(\frac{W_0}{W_{\text{final}}}\right) = \int_{x_{W,\text{final}}}^{x_F} \frac{1}{x_D - x_W}\,dx_W. \] The integral is evaluated numerically over \(N = 100\) steps: \[ \Delta x = \frac{x_F - x_{W,\text{final}}}{N} = \frac{0.08 - 0.04}{100} = 0.0004. \] At each step, the corresponding instantaneous distillate composition \(x_D\) is computed using the function find_xD_for_still with \(N_{\text{trays}} = 3\), \(R = 4.0\), and \(\alpha = 2.5\). Applying the trapezoidal rule yields: \[ I = \int_{0.04}^{0.08} \frac{1}{x_D - x_W}\,dx_W = 0.125285. \]
Compute the final pot mass. From the integrated Rayleigh equation: \[ W_{\text{final}} = W_0 \cdot e^{-I} = 1000.0 \cdot e^{-0.125285} = 882.246\ \text{kg}. \]
Determine the total distillate collected. \[ D_{\text{total}} = W_0 - W_{\text{final}} = 1000.0 - 882.246 = 117.754\ \text{kg}. \]
State the operating line parameters at the target distillate composition. The rectifying operating line for a batch column (no feed stage) is: \[ y = \frac{R}{R+1}x + \frac{x_{D,\text{target}}}{R+1}. \] Substituting the values: \[ \text{Slope} = \frac{R}{R+1} = \frac{4}{4+1} = 0.8, \] \[ \text{Intercept} = \frac{x_{D,\text{target}}}{R+1} = \frac{0.25}{4+1} = 0.05. \]

Final Answer

Theoretical stages (including reboiler): 4
Theoretical trays (excluding reboiler): 3
Actual trays required (at 60% efficiency): 6 (rounded up from 5.0 theoretical trays ÷ 0.60)
Final pot mass after operation: 882.246 kg
Total distillate collected: 117.754 kg
Rayleigh integral value: 0.125285
Operating line slope at target \(x_D\): 0.8; intercept: 0.05

"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