Flash Distillation Material Balance

Reference ID: MET-6D48 | Process Engineering Reference Sheets Calculation Guide

Introduction & Context

Flash distillation, also known as equilibrium flash vaporization, is a fundamental unit operation in process engineering used to separate a multi-component mixture into vapor and liquid phases. The process involves taking a pressurized feed stream and reducing its pressure (or increasing its temperature) to induce partial vaporization. The resulting mixture is then separated in a flash drum, where the vapor and liquid phases reach thermodynamic equilibrium.

This calculation is critical for designing separation equipment, sizing downstream piping, and determining the thermal duty required for heat exchangers, as detailed in the flash distillation energy balance. It is widely used in oil and gas processing, petrochemical refining, and chemical manufacturing to estimate the split of components between phases based on feed composition and operating conditions.

Methodology & Formulas

The flash distillation calculation relies on the simultaneous solution of total mass balances, component mass balances, and thermodynamic equilibrium relationships; for a visual approach, see the graphical solution for flash distillation. Given a feed rate F, feed composition zᵢ, and a vaporization fraction f, the system is solved as follows:

1. Total Mass Balance:

\[ F = V + L \]

Where the vapor flow rate V is determined by the vaporization fraction f:

\[ V = f \cdot F \]

\[ L = F \cdot (1 - f) \]

2. Component Mass Balance:

For any component i, the total amount in the feed must equal the sum of the amounts in the vapor and liquid phases:

\[ F \cdot z_{i} = V \cdot y_{i} + L \cdot x_{i} \]

3. Equilibrium Relationship:

Assuming constant relative volatility α for a binary system, the vapor phase composition y_i is linked to the liquid phase composition x_i:

\[ y_{i} = \frac{\alpha \cdot x_{i}}{1 + (\alpha - 1) \cdot x_{i}} \]

4. Quadratic Solution:

By substituting the equilibrium relationship into the component mass balance, we derive a quadratic equation in the form a · x_i² + b · x_i + c = 0, where:

\[ a = L \cdot (\alpha - 1) \]

\[ b = (V \cdot \alpha) + L - (F \cdot z_{i} \cdot (\alpha - 1)) \]

\[ c = -F \cdot z_{i} \]

The liquid composition x_i is then solved using the quadratic formula:

\[ x_{i} = \frac{-b + \sqrt{b^{2} - 4 \cdot a \cdot c}}{2 \cdot a} \]

Condition	Criteria
Vaporization Fraction	\( 0 < f < 1 \)
Relative Volatility	\( \alpha > 1 \) (for the more volatile component)
Composition Summation	\( \sum x_{i} = 1.0 \) and \( \sum y_{i} = 1.0 \)
Phase Regime	\( T_{\mathrm{bubble}} \leq T_{\mathrm{flash}} \leq T_{\mathrm{dew}} \)

The overall material balance is established by assuming a steady-state operation where the total molar feed rate equals the sum of the vapor and liquid product molar flow rates. The fundamental equation is:

\( F = V + L \)
Where \( F \) represents the total feed rate, \( V \) represents the vapor product rate, and \( L \) represents the liquid product rate.

For any specific component \( i \), the mass or molar balance must account for the distribution of that component between the two exiting phases. The component balance is expressed as:

\( F \cdot z_{i} = V \cdot y_{i} + L \cdot x_{i} \)
In this equation, \( z_{i} \) is the mole fraction in the feed, \( y_{i} \) is the mole fraction in the vapor phase, and \( x_{i} \) is the mole fraction in the liquid phase.

The equilibrium constant, denoted as \( K_{i} \), relates the vapor and liquid compositions at a given temperature and pressure. To solve the material balance, you must substitute \( y_{i} \) with \( K_{i} \cdot x_{i} \):

\( F \cdot z_{i} = V \cdot (K_{i} \cdot x_{i}) + L \cdot x_{i} \)
This allows process engineers to solve for \( x_{i} \) in terms of the feed composition and the vapor-liquid ratio.

The Rachford-Rice equation is a specialized form of the material balance used to determine the vapor fraction (\( V/F \)) without needing to know the individual phase compositions beforehand. It is defined as:

\[ \sum_{i} \frac{z_{i} \cdot (K_{i} - 1)}{1 + (V/F) \cdot (K_{i} - 1)} = 0 \]
This equation is solved iteratively to find the value of \( V/F \) that satisfies the equilibrium constraints.

Worked Example: Binary Flash Distillation of Benzene/Toluene

Scenario: A binary mixture of benzene and toluene is fed to a single-stage flash separator operating at steady-state. The vaporization fraction is known, and the system is ideal with constant relative volatility. Determine the product flow rates and compositions.

Knowns

Feed rate, \(F = 1000.0\) kmol/h
Feed mole fraction of benzene, \(z_{B} = 0.5\)
Vaporization fraction, \(f = 0.4\)
Operating pressure, \(P_{\mathrm{op}} = 1.0\) atm
Relative volatility of benzene to toluene, \(\alpha = 2.4\)
Equilibrium relationship: \(y_{B} = \dfrac{\alpha \cdot x_{B}}{1 + (\alpha - 1) \cdot x_{B}}\)

Step-by-Step Calculation

Determine total stream flow rates
Vapor flow rate: \(V = f \cdot F = 0.4 \cdot 1000.0 = 400.0\) kmol/h
Liquid flow rate: \(L = F - V = 1000.0 - 400.0 = 600.0\) kmol/h
Set up the component material balance
Overall benzene balance: \(F \cdot z_{B} = V \cdot y_{B} + L \cdot x_{B}\)
Substitute known values: \(1000.0 \cdot 0.5 = 400.0 \cdot y_{B} + 600.0 \cdot x_{B}\)
Equation A: \(500.0 = 400.0 \cdot y_{B} + 600.0 \cdot x_{B}\)
Introduce the equilibrium relationship
From the given VLE: \(y_{B} = \dfrac{2.4 \cdot x_{B}}{1 + (2.4 - 1) \cdot x_{B}} = \dfrac{2.4 \cdot x_{B}}{1 + 1.4 \cdot x_{B}}\) (Equation B)
Solve the coupled system
The component balance combined with VLE yields a quadratic equation in \(x_{B}\):
Let \(A = \alpha - 1 = 1.4\)
Quadratic coefficients:
\(a = L \cdot A = 600.0 \cdot 1.4 = 840.0\)
\(b = V \cdot \alpha + L - F \cdot z_{B} \cdot A = 400.0 \cdot 2.4 + 600.0 - 1000.0 \cdot 0.5 \cdot 1.4 = 960.0 + 600.0 - 700.0 = 860.0\)
\(c = -F \cdot z_{B} = -500.0\)
Discriminant: \(\Delta = b^{2} - 4 \cdot a \cdot c = 860.0^{2} - 4 \cdot 840.0 \cdot (-500.0) = 2419600.0\)
Positive root: \(x_{B} = \dfrac{-b + \sqrt{\Delta}}{2 \cdot a} = \dfrac{-860.0 + \sqrt{2419600.0}}{2 \cdot 840.0} = 0.414\)
Then from Equation B: \(y_{B} = \dfrac{2.4 \cdot 0.414}{1 + 1.4 \cdot 0.414} = 0.629\)
Compute remaining compositions
Liquid: \(x_{T} = 1 - x_{B} = 1 - 0.414 = 0.586\)
Vapor: \(y_{T} = 1 - y_{B} = 1 - 0.629 = 0.371\)
Verify the benzene balance
Benzene in: \(F \cdot z_{B} = 1000.0 \cdot 0.5 = 500.0\) kmol/h
Benzene out: \(V \cdot y_{B} + L \cdot x_{B} = 400.0 \cdot 0.629 + 600.0 \cdot 0.414 = 251.6 + 248.4 = 500.0\) kmol/h
The balance closes.

Final Results

Vapor stream: \(V = 400.0\) kmol/h, \(y_{B} = 0.629\), \(y_{T} = 0.371\)
Liquid stream: \(L = 600.0\) kmol/h, \(x_{B} = 0.414\), \(x_{T} = 0.586\)
Benzene balance error: \(0.0\) kmol/h

"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