Mass balance calculations determine batch distillate composition by tracking what leaves the reboiler, what accumulates in the receiver, and how the remaining liquid composition changes over time. For a batch still with no feed, the reboiler holdup (B) decreases as distillate is withdrawn, while the distillate accumulator holdup (H_a) increases. The instantaneous distillate composition (x_D) is obtained from the dynamic component balance and the vapor-liquid equilibrium behavior of the changing reboiler liquid.

The key relationship is conservation: the light component removed from the reboiler must equal the light component accumulated in the distillate, while the instantaneous vapor composition changes as the reboiler becomes depleted.

How the Overall Holdup Changes

Reboiler holdup decreases

Let (B(t)) be the liquid holdup in the reboiler and (L_D) or (D) be the distillate withdrawal rate. With no feed and negligible losses,

[ \frac{dB}{dt}=-D ]

Therefore,

[ B(t)=B_0-\int_0^t D,dt ]

For a constant distillate withdrawal rate,

[ B(t)=B_0-Dt ]

The reboiler holdup continues to decline until the planned endpoint, minimum operating inventory, or target bottoms composition is reached.

Distillate accumulator holdup increases

The receiver or distillate accumulator holdup (H_a) increases at the rate of product collection:

[ \frac{dH_a}{dt}=D ]

For constant withdrawal,

[ H_a(t)=\int_0^t D,dt=Dt ]

When vapor is condensed and reflux is returned internally to the column, only the net distillate withdrawal contributes to (H_a). Reflux changes the separation behavior but does not accumulate in the product receiver.

Total mass balance provides a check

For an idealized batch system containing only the reboiler and distillate receiver,

[ B_0=B(t)+H_a(t) ]

This overall balance is also a practical operating check. A significant discrepancy indicates measurement error, unaccounted liquid holdup, leakage, sampling losses, or an incorrect flow-rate assumption.

How Component Balances Determine Composition

The reboiler component balance

For the more volatile component, let (x_B) be its mole fraction in the reboiler liquid and (x_D) the mole fraction in the withdrawn distillate. The component balance is

[ \frac{d(Bx_B)}{dt}=-D x_D ]

The left-hand side represents the rate at which the volatile component inventory in the reboiler changes. The right-hand side represents its removal in the distillate.

Applying the product rule and using (dB/dt=-D) gives

[ B\frac{dx_B}{dt}=D(x_B-x_D) ]

[ \frac{dx_B}{dt}=\frac{D}{B}(x_B-x_D) ]

Because the distillate is normally richer in the volatile component than the reboiler liquid, (x_D>x_B), so (dx_B/dt) is negative. The volatile component concentration in the reboiler therefore declines during operation.

The accumulated distillate composition

The average composition of all material collected in the accumulator, (x_a), follows from the integrated component balance:

[ H_a x_a=B_0x_{B0}-B_1x_{B1} ]

Thus,

[ x_a=\frac{B_0x_{B0}-B_1x_{B1}}{H_a} ]

Here, (B_0) and (x_{B0}) are the initial reboiler holdup and composition, while (B_1) and (x_{B1}) are their values at the selected operating time.

This is an average composition over the collected batch, not the composition of the most recent distillate. The two values differ whenever (x_D) changes with time.

The instantaneous distillate composition

The instantaneous (x_D) is found from the current reboiler state and column operating conditions. At each time step, the process model typically:

Determines the current (B) and (x_B).
Calculates the bubble-point temperature at the operating pressure.
Uses vapor-liquid equilibrium to determine the vapor composition leaving the reboiler or column top.
Applies the column mass-transfer or equilibrium-stage model at the selected reflux ratio.
Identifies the resulting overhead composition as the instantaneous distillate composition (x_D(t)).

For a simple equilibrium relationship,

[ y_i=\frac{x_i P_i^{\text{sat}}(T)}{P} ]

with the temperature selected to satisfy the bubble-point condition. Antoine equations or another vapor-pressure correlation can provide (P_i^{\text{sat}}(T)).

In a real pilot column, the top distillate composition may not equal the vapor composition directly above the reboiler. The trays, packing, reflux, pressure, efficiency, and column holdup determine how the changing bottoms composition is translated into the overhead composition.

Why Distillate Composition Changes During the Run

The volatile component is progressively depleted

At the beginning of a batch, the reboiler contains its highest concentration of the more volatile component. The vapor and overhead product are therefore relatively rich in that component.

As operation continues, the volatile component is preferentially removed. The reboiler liquid becomes richer in the less volatile component, causing the equilibrium vapor composition and usually the distillate purity to decline.

Reflux controls the degree of enrichment

At a higher reflux ratio, more condensed overhead liquid returns to the column instead of leaving as product. This provides additional vapor-liquid contacting and generally increases the instantaneous distillate purity.

The trade-off is a lower net product withdrawal rate for a given boil-up rate, so a higher reflux ratio can increase batch duration and energy consumption.

Column holdup creates a dynamic delay

Liquid retained on trays, packing, piping, or in the condenser does not respond instantaneously to a change in reboiler composition. These holdups contain material with their own compositions and must be included when a detailed transient model is required.

For a section with holdup (H), the total balance has the general form

[ \frac{dH}{dt}=\text{inflows}-\text{outflows} ]

The component balance is

[ \frac{d(Hx_i)}{dt}

\sum \text{component inflows}

\sum \text{component outflows} ]

These equations explain why a pilot plant can continue producing a relatively concentrated distillate for some time after the reboiler composition has begun to change.

How the Calculations Are Used Operationally

Establishing a composition profile

The dynamic balances are solved over successive time intervals. At each interval, the calculation updates the reboiler holdup, reboiler composition, equilibrium temperature, vapor composition, overhead composition, and accumulated product composition.

The result is a profile such as

[ B(t),\quad x_B(t),\quad x_D(t),\quad H_a(t),\quad x_a(t) ]

This profile shows both the current product quality and the total material collected.

Selecting main-cut and off-cut times

A main cut is collected while the instantaneous distillate composition remains within the required specification. When (x_D(t)) falls below that specification, the operator can switch receivers or divert the product to an off-cut vessel.

The switch time should be based on the predicted or measured instantaneous composition, not only on the average composition in the receiver. A receiver can still have acceptable average purity even after the latest distillate has fallen below specification.

Estimating the required batch duration

Suppose the process must reduce the volatile-component mole fraction in the reboiler from (x_{B0}) to a target (x_{B1}). The required time is obtained by integrating the composition balance:

[ \frac{dx_B}{dt}=\frac{D}{B}(x_B-x_D) ]

Because both (B) and (x_D) can vary with time, the equation is generally solved numerically. The result gives the operating time, remaining reboiler inventory, collected distillate quantity, and average product composition at the target state.

Comparing model predictions with measurements

Operators can compare calculated and measured values of:

Reboiler holdup.
Distillate receiver holdup.
Reboiler composition.
Instantaneous distillate composition.
Average receiver composition.
Batch duration and product yield.

This comparison can reveal incorrect flow measurements, unmodeled column holdup, poor equilibrium assumptions, heat losses, pressure changes, or tray and packing inefficiency.

Understanding the Trade-offs

Higher purity versus production rate

Increasing reflux generally improves separation and raises distillate purity. However, it reduces the fraction of condensed overhead that is withdrawn as product, increasing the time needed to collect a specified amount.

Simple balances versus rigorous models

A reboiler-and-receiver balance is often sufficient for estimating total product quantity and average composition. It is less reliable for predicting the timing of composition transitions when tray, packing, condenser, or piping holdups are significant.

Constant-rate assumptions versus actual operation

The simplified relationship (H_a=Dt) assumes a constant net distillate rate. Actual pilot plants may experience changing boil-up, reflux, pressure, condensation, or control-valve behavior, in which case the measured or modeled (D(t)) must be integrated.

Instantaneous composition versus cumulative composition

Confusing (x_D) with (x_a) is a common operational error. (x_D) describes the product leaving at a particular moment, whereas (x_a) describes everything collected up to that moment.

Numerical accuracy and physical data

The dynamic calculation is only as reliable as its vapor-liquid equilibrium data, pressure measurement, flow-rate data, and initial composition. Numerical integration methods such as Runge-Kutta or RK45 improve solution accuracy, but they cannot correct inaccurate physical assumptions or measurements.

Making the Right Choice for Your Goal

The balances should be selected according to the decision the pilot plant must support.

If your primary focus is predicting distillate composition: Use a dynamic reboiler balance coupled with vapor-liquid equilibrium and a column model that includes the selected reflux ratio.
If your primary focus is tracking holdup: Integrate the total balances for the reboiler and accumulator, and include tray, packing, condenser, and piping holdups when their inventories are material.
If your primary focus is setting receiver-switch times: Track the instantaneous (x_D(t)), compare it with the product specification, and do not rely solely on the receiver average (x_a).
If your primary focus is estimating product yield: Use (H_a=B_0-B) together with (H_a x_a=B_0x_{B0}-Bx_B), while accounting for any retained column inventory.
If your primary focus is optimizing reflux: Compare purity, recovery, batch duration, and energy use over the predicted dynamic composition profile.

A correctly formulated dynamic mass balance turns changing pilot-plant measurements into a defensible prediction of composition, holdup, product yield, and cut-off time.

Summary Table:

Parameter	Calculation or relationship	Operational use
Reboiler holdup	`dB/dt = -D`	Tracks remaining batch inventory
Accumulator holdup	`dHₐ/dt = D`	Measures collected distillate
Reboiler composition	`dxB/dt = (D/B)(xB - xD)`	Predicts volatile-component depletion
Instantaneous distillate	Determined from VLE, reflux, pressure, and column behavior	Sets current product quality
Average distillate	`xₐ = (B₀xB₀ - B₁xB₁)/Hₐ`	Evaluates cumulative batch purity
Column holdup	`d(Hxᵢ)/dt = component inflows - outflows`	Predicts composition delays and cut transitions

Need reliable batch distillation data for teaching, process development, or production? Contact LABPARK to discuss Educational and Vocational Unit Operations Pilot Plants for chemical engineering, bioprocess and biotech, and environmental and water treatment applications. LABPARK helps universities, research institutes, and enterprises evaluate dynamic mass balances, reflux strategies, holdup effects, product quality, and receiver-switch timing with practical pilot-scale systems tailored to your process goals.