Two-Phase Flow: Predicting Regimes, Slugging & Line Sizing Hydraulics

1. Physics of Concurrent Gas-Liquid Flow

Concurrent flow of gas and liquid in a closed conduit is among the most challenging fluid mechanics problems encountered in petroleum engineering and chemical processing. Common industrial examples include:

• Upstream wellhead flowlines and subsea tiebacks carrying crude oil, produced water, and associated gas.
• Boiler and fired heater tubes experiencing nucleate boiling and convective vaporization.
• Distillation column overhead lines leading to partial condensers.
• Thermodynamic relief headers discharging two-phase flashing liquids during plant emergencies.

Single-phase hydraulic equations (such as Darcy-Weisbach) fail completely in two-phase systems because of three interconnected physical realities:

1. Phase Slip & Holdup: Gas, having vastly lower density and viscosity than liquid, travels at a substantially higher in-situ velocity ($v_g > v_l$). The ratio of their velocities is the slip ratio ($S = v_g / v_l \ge 1.0$). Consequently, liquid accumulates inside the pipe, making the in-situ liquid volume fraction (liquid holdup, $H_l$) substantially greater than the input volume fraction (no-slip liquid fraction, $\lambda_l$).
2. Deformable Phase Interface: Unlike a solid boundary, the interface between gas and liquid continuously deforms, generating interfacial waves, shear ripples, droplets, and liquid bridges.
3. Interfacial Shear Stress: High-speed gas transfers momentum to the liquid through interfacial friction, generating complex pressure gradient profiles that can be orders of magnitude higher than single-phase flow.

2. Two-Phase Flow Regimes in Pipes

Depending on relative phase flow rates, fluid properties, pipe diameter, and inclination angle, two-phase flow organizes into distinct spatial geometries known as flow regimes:

2.1 Horizontal & Slightly Inclined Pipes

2.2 Vertical Upflow Regimes

In vertical risers, buoyancy acts in the flow direction. Regimes progress from Bubble flow $\to$ Slug flow (Taylor bubbles occupying $\approx 90\%$ of pipe area) $\to$ Churn flow (highly chaotic, oscillatory churning) $\to$ Annular flow.

3. Slugging Mechanics: Hydrodynamic vs. Severe Terrain Slugging

Slugs originate from two distinct physical mechanisms that require completely different mitigation strategies:

3.1 Hydrodynamic Slugging

In horizontal lines, the Kelvin-Helmholtz instability causes surface waves to grow when the Bernoulli suction pressure drop above a wave crest overcomes gravitational restoring forces. Hydrodynamic slugs are short ($10\text{--}40 \times D$) and frequent (frequencies of $0.1\text{--}1.0\,\text{Hz}$). They are managed using standard slug catchers and robust pipe supports.

3.2 Severe Terrain-Induced Slugging (Riser Slugging)

Severe slugging occurs in offshore pipelines with descending sea-floor terrain terminating in a vertical riser to a production platform:

1. Slug Formation: Liquid accumulates at the riser base (low point), forming a hydrostatic liquid column that blocks gas passage.
2. Slug Growth: The liquid level rises up the entire height of the vertical riser while upstream pipeline pressure builds.
3. Blowout: Once upstream gas pressure exceeds the hydrostatic head of the full riser ($\Delta P_{gas} > \rho_l g H_{riser}$), the entire liquid column is violently ejected into the topside separator at velocities up to 20–30 m/s.
4. Liquid Fallback: Following gas depressurization, residual liquid falls back down the riser, and the cycle repeats with periods of 10 to 60 minutes.

Riser slugging can completely flood topside separation trains and trip platforms. Mitigation requires topside choking (increasing backpressure to stabilize the gas-liquid interface) or automated active slug controllers.

4. Multiphase Pressure Drop Modeling

4.1 Homogeneous No-Slip Model

The simplest model treats the two-phase mixture as a single pseudo-fluid with homogeneous properties:

The homogeneous model is accurate only for dispersed bubble or high-pressure mist flow ($P > 70\,\text{bar}$). At lower pressures, it severely underestimates liquid holdup and pressure drop.

4.2 Lockhart-Martinelli Separated Flow Model

The classic separated flow method introduces the Martinelli parameter ($\chi$), defined as the square root of the ratio of single-phase liquid to gas pressure gradients:

The two-phase frictional pressure drop is then calculated using the two-phase multiplier ($\phi_l^2$):

Chisholm formulated the analytical correlation for $\phi_l^2$:

where $C$ depends on the flow regimes of the individual phases ($C = 20$ for turbulent liquid / turbulent gas).

5. API RP 14E Erosional Velocity & Momentum Limits

Two-phase piping must be sized between two strict hydraulic boundaries:

1. Minimum Velocity (Slugging Prevention): Velocity must be high enough to sweep liquid continuously along the pipe, preventing massive liquid dropout in low terrain ($v_m \ge 3\text{--}5\,\text{m/s}$).
2. Maximum Velocity (Erosional Limit): Velocity must not exceed the API RP 14E erosional threshold to prevent pipe wall erosion-corrosion: $$v_e = \frac{C}{\sqrt{\rho_m}}$$ where $C = 100$ for continuous service in solids-free carbon steel, $C = 125$ for intermittent relief lines, and $C = 150\text{--}200$ for corrosion-resistant alloys (duplex stainless steel, Inconel).

Additionally, the mixture momentum flux ($\rho_m v_m^2$) must be evaluated for pipe support design:

6. Step-by-Step Worked Engineering Example

Evaluate a subsea production flowline transporting wet gas and condensate:

• Gas flow rate: $Q_g = 6.0\,\text{MMSCFD} = 7{,}080\,\text{Sm}^3\text{/h} = 1.967\,\text{Sm}^3\text{/s}$
• Condensate flow rate: $Q_l = 2{,}500\,\text{BPD} = 16.56\,\text{m}^3\text{/h} = 0.00460\,\text{m}^3\text{/s}$
• Operating pressure: $P = 40.0\,\text{bar a}$ ($4000\,\text{kPa}$)
• Operating temperature: $T = 35^\circ\text{C} = 308.15\,\text{K}$
• Gas properties at operating conditions: $\rho_g = 32.5\,\text{kg/m}^3$, $\mu_g = 1.30 \times 10^{-5}\,\text{Pa}\cdot\text{s}$
• Condensate properties: $\rho_l = 780\,\text{kg/m}^3$, $\mu_l = 1.20 \times 10^{-3}\,\text{Pa}\cdot\text{s}$
• Trial pipe: $6^{\prime\prime}$ Sch 80 carbon steel ($D = 146.3\,\text{mm} = 0.1463\,\text{m}$, Area $A = 0.01681\,\text{m}^2$)

Step 1: Calculate In-Situ Volumetric Flow Rates

Gas volumetric flow rate at operating conditions ($Z \approx 0.86$):

Liquid volumetric flow rate: $Q_l = 0.00460\,\text{m}^3\text{/s}$.

Total mixture flow rate: $Q_m = 0.0458 + 0.00460 = 0.0504\,\text{m}^3\text{/s}$.

Step 2: Superficial & Mixture Velocities

No-slip liquid input volume fraction:

Step 3: Calculate Mixture Density & API 14E Erosional Velocity

Homogeneous mixture density:

API RP 14E Erosional Velocity ($C = 100$ for continuous carbon steel):

Checking operating velocity margin:

Velocity is safely below the erosional threshold.

Step 4: Check Momentum Flux ($\rho_m v_m^2$)

This is well below the $5{,}000\,\text{Pa}$ threshold, indicating minimal pipe vibration and low cyclic stress on pipe bends.

Step 5: Determine Flow Regime (Mandhane / Beggs-Brill Map)

With superficial velocities $v_{sg} = 2.73\,\text{m/s}$ and $v_{sl} = 0.27\,\text{m/s}$ in a 6-inch horizontal pipe:

• The flow regime lies squarely in the Stratified Wavy to gentle Elongated Bubble transition boundary.
• Because $v_m \approx 3.0\,\text{m/s}$, liquid is continuously swept along the pipe invert without building into destructive high-energy slugs.
• Verdict: The 6-inch Sch 80 pipe is hydraulically well-sized.

7. ChemProCal Integration

Model complex multiphase piping with ChemProCal's specialized calculators:

• Two-Phase Line Sizing Tool: Comprehensive Beggs-Brill and Lockhart-Martinelli hydraulic modeling.
• Separator Utilities: Size slug catchers to handle upstream pipeline liquid inventory.
• Knock-Out Drum Calculator: Design the downstream reception vessel to separate multiphase fluids.