🎓 Lesson 13
D5
Isothermal Flash Algorithm: Newton-Raphson Implementation
It’s a step-by-step math method that finds how much liquid and vapor form when a mixture is flashed at constant temperature and pressure.
🎯 Learning Objectives
- ✓ Calculate the vapor fraction (V/F) and equilibrium compositions for a 3-component hydrocarbon mixture using the Newton-Raphson isothermal flash algorithm
- ✓ Analyze convergence behavior by monitoring residual norms and Jacobian conditioning in successive iterations
- ✓ Apply thermodynamic consistency checks (e.g., ∑x_i = 1, ∑y_i = 1, K_i = φ_i^L/φ_i^V) to validate flash results
- ✓ Explain the role of the Rachford-Rice equation in eliminating explicit vapor/liquid split dependencies
- ✓ Design initial guesses for V/F and K-values that ensure robust convergence for near-critical mixtures
📖 Why This Matters
In mining and mineral processing, isothermal flash calculations underpin solvent extraction, acid leaching off-gas treatment, and volatile reagent recovery (e.g., HCl, NH₃, organic solvents). For example, flash drums in hydrometallurgical circuits must precisely separate vapors from pregnant leach solutions—getting the vapor fraction wrong risks corrosion, emissions violations, or loss of valuable reagents. The Newton-Raphson implementation ensures speed and reliability in process simulators used daily by plant engineers and designers.
📘 Core Principles
The isothermal flash problem begins with known T, P, and overall mole fractions z_i. At equilibrium, each component distributes between vapor (y_i) and liquid (x_i) phases according to K_i = y_i/x_i. Material balance gives y_i = (V/F)·K_i·x_i + (1 − V/F)·x_i → leading to the Rachford-Rice equation: ∑[z_i(K_i−1)/(1 + (V/F)(K_i−1))] = 0. This single nonlinear equation in V/F is solved via Newton-Raphson: V/F_{k+1} = V/F_k − f(V/F_k)/f′(V/F_k). For multicomponent systems, the full Newton-Raphson method extends to a vector form solving both V/F and x_i (or y_i) simultaneously using the Jacobian of residuals—enabling convergence even when K-values are uncertain or near unity.
📐 Rachford-Rice Equation & Newton-Raphson Update
The Rachford-Rice equation reduces the isothermal flash to one variable (V/F), assuming K-values are known (e.g., from Wilson or UNIFAC correlations). Its derivative provides the Newton-Raphson slope for rapid convergence. When K_i values are strongly composition-dependent, a full 2N−1 variable Newton scheme (solving for V/F and all x_i except one) is required.
💡 Worked Example
Problem: A 3-component mixture (z₁=0.4, z₂=0.3, z₃=0.3) at 60°C and 200 kPa has K-values: K₁=2.5, K₂=0.8, K₃=0.3. Find V/F using Newton-Raphson (tolerance = 1e−5).
1.
Step 1: Define f(V/F) = ∑[z_i(K_i−1)/(1 + (V/F)(K_i−1))]. Initial guess: V/F₀ = 0.5.
2.
Step 2: Compute f(0.5) = 0.4(1.5)/(1+0.5×1.5) + 0.3(−0.2)/(1+0.5×−0.2) + 0.3(−0.7)/(1+0.5×−0.7) = 0.4615 − 0.0667 − 0.2857 = 0.1091.
3.
Step 3: Compute f′(V/F) = −∑[z_i(K_i−1)²/(1 + (V/F)(K_i−1))²] = −[0.4(2.25)/1.5625 + 0.3(0.04)/0.81 + 0.3(0.49)/0.6025] = −[0.576 + 0.0148 + 0.2439] = −0.8347.
4.
Step 4: Update: V/F₁ = 0.5 − (0.1091)/(−0.8347) = 0.6306. Repeat until |f| < 1e−5 → converges at V/F = 0.6482.
5.
Step 5: Compute x_i = z_i / [1 + (V/F)(K_i−1)]: x₁ = 0.4 / [1 + 0.6482×1.5] = 0.206; x₂ = 0.3 / [1 − 0.6482×0.2] = 0.343; x₃ = 0.3 / [1 − 0.6482×0.7] = 0.451. Sum = 1.000 ✓.
Answer:
The converged vapor fraction is V/F = 0.648, with liquid composition x = [0.206, 0.343, 0.451], satisfying ∑x_i = 1.000 within rounding error.
🏗️ Real-World Application
At the Tenke Fungurume copper-cobalt mine (DRC), sulfuric acid regeneration units use vacuum flash drums to recover H₂SO₄ vapor from spent electrolyte. An isothermal flash algorithm—implemented in Aspen HYSYS with Newton-Raphson core—models the 7-component (H₂O, H₂SO₄, Cu²⁺, Co²⁺, Fe³⁺, Al³⁺, Mn²⁺) aqueous system at 80°C and 15 kPa. Accurate V/F prediction ensures acid concentration in the bottoms stays >93 wt% to prevent corrosion in downstream concentrators—deviations >0.5% V/F cause measurable stainless-steel pitting rates per ASTM G140.
✏️ Student Exercise
Given a binary mixture of methanol (z₁ = 0.6) and water (z₂ = 0.4) at 80°C and 100 kPa, with K₁ = 2.1 and K₂ = 0.45: (a) Write the Rachford-Rice equation explicitly; (b) Perform two Newton-Raphson iterations starting from V/F₀ = 0.4; (c) Report final V/F and x₁, y₁ after convergence to 1e−4 tolerance. Verify ∑x_i = 1 and ∑y_i = 1.