📑 Table of Contents (Tap to view sections)
1. Flash Drum Thermodynamic Principles
A Flash Drum (also known as an equilibrium flash vessel or vapor-liquid separator) is an essential unit operation designed to separate a multicomponent fluid mixture into coexisting vapor and liquid phases following a thermodynamic state transition. In industrial chemical processing and oil & gas production, this transition is almost universally achieved by adiabatic throttling across a control valve (isenthalpic expansion) or through thermal duty addition in a feed preheater.
Key industrial applications include:
- Upstream Production Separators: Multi-stage stabilization of crude oil and rich condensate to reduce Reid Vapor Pressure (RVP) for safe storage and pipeline transport.
- Refrigeration Chiller Economizers: Interstage flash economization in propane or ethylene refrigeration loops to boost compressor thermodynamic efficiency.
- Petrochemical Depressurization: Flashing reactor effluent streams to quickly disengage unreacted volatile monomers (e.g., ethylene, propylene) from heavy polymer solutions.
Flash drum sizing combines Rachford-Rice Vapor-Liquid Equilibrium (VLE) thermodynamics with API Specification 12J, GPSA Section 7 gravity separation dynamics, and API RP 14E erosional velocity criteria.
2. Vapor-Liquid Equilibrium (VLE) & Rachford-Rice Derivation
For an $N$-component mixture with feed molar composition $z_i$ ($i = 1, 2, \dots, N$), flash separation at specified temperature $T$ and pressure $P$ produces an equilibrium vapor stream of molar flow $V$ (composition $y_i$) and a liquid stream of molar flow $L$ (composition $x_i$).
2.1 Equilibrium Ratio ($K$-Value)
Thermodynamic phase equilibrium requires equality of chemical potentials, or equivalently, equal fugacities in both phases ($f_i^V = f_i^L$). This defines the equilibrium distribution coefficient ($K_i$):
$$K_i = \frac{y_i}{x_i} = \frac{\phi_i^L(T, P, \mathbf{x})}{\phi_i^V(T, P, \mathbf{y})}$$For ideal mixtures at low pressures, Raoult-Dalton's Law applies ($K_i = P_i^{sat}(T) / P$). For non-ideal hydrocarbon systems at high pressures, $K$-values must be evaluated using cubic equations of state (such as Peng-Robinson or Soave-Redlich-Kwong).
2.2 Material Balance & Derivation of Rachford-Rice Equation
Define the vapor molar flash fraction as:
$$\psi = \frac{V}{F} \implies \frac{L}{F} = 1 - \psi$$where $F$ is the total feed molar flow. Overall and component material balances yield:
$$F z_i = V y_i + L x_i = F \psi y_i + F (1 - \psi) x_i$$ $$z_i = \psi y_i + (1 - \psi) x_i$$Substituting $y_i = K_i x_i$ gives the liquid phase mole fraction:
$$x_i = \frac{z_i}{1 + \psi(K_i - 1)}$$Similarly, the vapor phase mole fraction is:
$$y_i = \frac{z_i K_i}{1 + \psi(K_i - 1)}$$Since the mole fractions in each phase must independently sum to unity ($\sum x_i = 1$ and $\sum y_i = 1$), their difference must be identically zero:
$$\sum_{i=1}^N (y_i - x_i) = 0$$Substituting $y_i$ and $x_i$ yields the celebrated Rachford-Rice Objective Function:
$$f(\psi) = \sum_{i=1}^N \frac{z_i (K_i - 1)}{1 + \psi(K_i - 1)} = 0$$2.3 Numerical Convergence Properties
The Rachford-Rice equation is mathematically unique because its derivative with respect to $\psi$ is strictly negative monotonic across the physical domain ($\psi \in [0, 1]$):
$$f'(\psi) = -\sum_{i=1}^N \frac{z_i (K_i - 1)^2}{[1 + \psi(K_i - 1)]^2} < 0$$Because $f'(\psi) < 0$ everywhere, $f(\psi)$ has exactly one unique real root between 0 (bubble point boundary) and 1 (dew point boundary). The root is calculated rapidly using the standard Newton-Raphson iteration:
$$\psi^{(k+1)} = \psi^{(k)} - \frac{f(\psi^{(k)})}{f'(\psi^{(k)})}$$Convergence to $10^{-8}$ absolute tolerance is typically achieved in fewer than 5 iterations.
3. Flash Drum Physical Sizing Methodology
Once the thermodynamic equilibrium split ($\psi$, phase compositions, mass flows, and fluid densities $\rho_v, \rho_l$) is determined, the physical vessel is sized using gravity separation principles.
3.1 Maximum Vapor Velocity (Souders-Brown)
The maximum allowable upward vapor velocity to avoid liquid droplet carryover is:
$$v_{max} = K_{SB} \sqrt{\frac{\rho_l - \rho_v}{\rho_v}}$$For vertical flash drums:
- With Wire Mesh Demister Pad: $K_{SB} = 0.09\text{--}0.107\,\text{m/s}$ ($0.30\text{--}0.35\,\text{ft/s}$), derated for pressure per GPSA Section 7.
- Without Demister (Bare Vessel): $K_{SB} = 0.046\text{--}0.061\,\text{m/s}$ ($0.15\text{--}0.20\,\text{ft/s}$).
3.2 Vessel Diameter Calculation
Setting the design vapor velocity to $v_{design} = 0.80\text{--}0.85 \times v_{max}$:
$$A_{vessel} = \frac{Q_v}{v_{design}} = \frac{W_v / \rho_v}{v_{design}} \implies D_{vessel} = \sqrt{\frac{4 A_{vessel}}{\pi}}$$3.3 Liquid Holdup & Surge Capacity (API 12J)
The liquid surge volume $V_L$ must provide sufficient residence time for process control stability, automated level dumping, and upstream/downstream operational decoupling:
| Downstream Destination | Minimum Surge Time (NLL to LLL) | Recommended Total Holdup |
|---|---|---|
| Storage Tank (Atmospheric) | $3\text{--}5\,\text{minutes}$ | $5\text{--}10\,\text{minutes}$ |
| Fired Heater / Reboiler Feed | $8\text{--}10\,\text{minutes}$ | $15\text{--}20\,\text{minutes}$ |
| Fractionation Column Feed | $5\text{--}8\,\text{minutes}$ | $10\text{--}15\,\text{minutes}$ |
4. Two-Phase Inlet Nozzle Sizing & API RP 14E
Downstream of the flash control valve, the fluid is a turbulent two-phase mixture. Sizing the flash drum inlet nozzle requires evaluating the API RP 14E erosional velocity:
$$v_e = \frac{C}{\sqrt{\rho_m}}$$where $\rho_m$ is the homogeneous mixture density:
$$\rho_m = \left[\frac{x_m}{\rho_v} + \frac{1 - x_m}{\rho_l}\right]^{-1}$$and $x_m$ is the mass vapor fraction. The empirical empirical factor $C = 100$ for continuous service in solids-free service, and $C = 125$ for intermittent service.
The actual inlet nozzle velocity $v_{in}$ should satisfy:
$$v_{in} \le 0.80 \times v_e$$Additionally, to avoid shattering liquid droplets at the inlet, the momentum flux must be kept within limits: $\rho_m v_{in}^2 \le 1400\,\text{Pa}$ (without baffle) or $\le 4000\text{--}6000\,\text{Pa}$ (with inlet half-pipe / Schoepentoeter).
5. Step-by-Step Worked Numerical Example
A rich hydrocarbon condensate stream is flashed across a pressure reducing valve into a vertical flash drum:
- Total feed molar flow rate: $F = 1{,}000\,\text{kmol/h}$
- Flash Drum Conditions: $T = 40^\circ\text{C}$ ($313.15\,\text{K}$), $P = 8.0\,\text{bar a}$
- Liquid product destination: Fractionation deethanizer ($t_{hold} = 6\,\text{minutes}$)
- Vessel equipped with a standard stainless steel demister pad ($K_0 = 0.107\,\text{m/s}$)
Feed Composition & Equilibrium $K$-Values (Peng-Robinson EOS):
| Component | Feed Mole Frac ($z_i$) | $K$-Value ($K_i$) | $K_i - 1$ | MW (kg/kmol) |
|---|---|---|---|---|
| Methane ($C_1$) | $0.25$ | $8.20$ | $+7.20$ | $16.04$ |
| Ethane ($C_2$) | $0.15$ | $1.85$ | $+0.85$ | $30.07$ |
| Propane ($C_3$) | $0.20$ | $0.65$ | $-0.35$ | $44.10$ |
| n-Butane ($nC_4$) | $0.20$ | $0.22$ | $-0.78$ | $58.12$ |
| n-Pentane ($nC_5$) | $0.20$ | $0.065$ | $-0.935$ | $72.15$ |
Step 1: Solve Rachford-Rice Objective Function
Solving $f(\psi) = 0$ via Newton-Raphson:
At iteration 1 ($\psi^{(0)} = 0.350$): $f(0.350) = +0.0214$
At iteration 2 ($\psi^{(1)} = 0.3625$): $f(0.3625) = -0.0003$
Converged vapor flash fraction: $\psi = 0.3621$ (36.21% molar vapor fraction).
Step 2: Calculate Phase Compositions & Mass Flows
- Vapor molar flow: $V = 1{,}000 \times 0.3621 = 362.1\,\text{kmol/h}$
- Liquid molar flow: $L = 1{,}000 \times (1 - 0.3621) = 637.9\,\text{kmol/h}$
Evaluating phase compositions ($y_i, x_i$):
- $y_{C1} = 0.25 \times 8.20 / [1 + 0.3621(7.20)] = 2.050 / 3.6071 = 0.5683$
- $y_{C2} = 0.15 \times 1.85 / [1 + 0.3621(0.85)] = 0.2775 / 1.3078 = 0.2122$
- $y_{C3} = 0.20 \times 0.65 / [1 + 0.3621(-0.35)] = 0.1300 / 0.8733 = 0.1489$
- $y_{nC4} = 0.20 \times 0.22 / [1 + 0.3621(-0.78)] = 0.0440 / 0.7176 = 0.0613$
- $y_{nC5} = 0.20 \times 0.065 / [1 + 0.3621(-0.935)] = 0.0130 / 0.6614 = 0.0197$
Sum of vapor mole fractions: $\sum y_i = 1.0104 \approx 1.000$.
Vapor molecular weight: $MW_v = \sum y_i MW_i = 25.43\,\text{kg/kmol}$.
Liquid molecular weight: $MW_l = 61.22\,\text{kg/kmol}$.
Total vapor mass flow: $W_v = 362.1 \times 25.43 = 9{,}208\,\text{kg/h} = 2.558\,\text{kg/s}$.
Total liquid mass flow: $W_l = 637.9 \times 61.22 = 39{,}052\,\text{kg/h} = 10.848\,\text{kg/s}$.
Step 3: Physical Properties & Vessel Sizing
At $P = 8.0\,\text{bar a}$, $T = 40^\circ\text{C}$ ($Z \approx 0.92$):
$$\rho_v = \frac{P \cdot MW_v}{Z R T} = \frac{800 \times 25.43}{0.92 \times 8.314 \times 313.15} = 8.49\,\text{kg/m}^3$$Liquid density from EOS: $\rho_l = 585\,\text{kg/m}^3$.
Derated Souders-Brown $K$-factor ($P = 8\,\text{bar a} = 116\,\text{psia}$):
$$K_{derated} = 0.107 \times [1 - 0.0001(116 - 100)] = 0.1068\,\text{m/s}$$Maximum allowable gas velocity:
$$v_{max} = 0.1068 \times \sqrt{\frac{585 - 8.49}{8.49}} = 0.1068 \times \sqrt{67.90} = 0.880\,\text{m/s}$$Design gas velocity ($85\%$ of $v_{max}$): $v_{design} = 0.85 \times 0.880 = 0.748\,\text{m/s}$.
Vapor volumetric flow rate:
$$Q_v = \frac{W_v}{\rho_v} = \frac{2.558\,\text{kg/s}}{8.49\,\text{kg/m}^3} = 0.3013\,\text{m}^3\text{/s}$$Required internal diameter ($D$):
$$A_{vessel} = \frac{0.3013}{0.748} = 0.4028\,\text{m}^2 \implies D = \sqrt{\frac{4 \times 0.4028}{\pi}} = 0.716\,\text{m}$$Select standard nominal shell: $D = 750\,\text{mm}$ ($30\,\text{inches}$). Actual area: $A_{act} = 0.4418\,\text{m}^2$.
Step 4: Liquid Surge Height
Liquid volumetric flow rate:
$$Q_l = \frac{39{,}052\,\text{kg/h}}{585\,\text{kg/m}^3} = 66.76\,\text{m}^3\text{/h} = 0.01854\,\text{m}^3\text{/s}$$Required surge volume for 6 minutes holdup:
$$V_{surge} = 66.76 \times \left(\frac{6}{60}\right) = 6.676\,\text{m}^3$$Liquid surge height:
$$H_{surge} = \frac{V_{surge}}{A_{act}} = \frac{6.676}{0.4418} = 15.11\,\text{m}\quad (\text{Vessel would be excessively tall!})$$Because the liquid volume dominates over the vapor load, a vertical vessel of 750 mm diameter would require a ridiculous $L/D$ ratio of $> 22$. The optimal engineering solution is to switch to a horizontal flash drum ($D = 1.8\,\text{m}$, $L = 5.4\,\text{m}$, $L/D = 3.0$), or increase the vertical diameter to $1600\,\text{mm}$ to provide the required liquid surge height within $3.3\,\text{m}$.
6. Common Operational Failure Modes
- Level Control Oscillations (Foaming): When hydrocarbons contain dissolved asphaltenes, amine traces, or corrosion inhibitors, rapid pressure drop creates micro-bubbles and persistent foam. Foam carryover into the vapor line can trigger compressor trips. Antifoam injection or vane demisters are required.
- Hydrate Formation Across the Throttle Valve: Flashing wet natural gas causes Joule-Thomson cooling, easily dropping temperatures into the gas clathrate hydrate envelope. Continuous MEG/methanol injection upstream of the valve is standard design practice.
- Choked Flow at Inlet Nozzle: If the inlet pipe is undersized, critical sonic velocity occurs at the vessel entrance nozzle, generating severe acoustic fatigue ($\text{dB} > 115$) and droplet shattering.
7. ChemProCal Integration
Model flash equilibrium and vessel sizing using ChemProCal calculation tools:
- Flash Drum Sizing Tool: Complete thermodynamic flash and vessel geometry calculations.
- Gas Intelligence Engine (EOS): Calculate multicomponent Peng-Robinson $K$-values and fugacities.
- Vertical Separator Sizing Tool: Detailed nozzle momentum and height sizing.
Flash Drum Sizing
Apply this methodology directly in the ChemProCal calculator.
Open Calculator →Live Separator Souders-Brown Capacity Estimator
Adjust parameters below to test the methodology equations in real time before running full simulations: