Slug Flow Regime Map (Taitel-Dukler 1976) — Fundamentals

Two-phase flow pattern prediction, F/K/T dimensionless groups, severe slugging, and how regime drives equipment design.

1. Why regime matters

Gas-liquid two-phase flow can take fundamentally different forms in the same pipe depending on the ratio of phases and the geometry. Each form has different pressure drop, holdup, heat transfer, erosion behavior, and operational risk. Slug flow is the most challenging — a liquid plug repeatedly accelerates through the pipe at near-mixture velocity, slamming into bends, fittings, and downstream vessels.

Knowing the regime up front lets you size the right downstream vessel (slug catcher vs. separator vs. scrubber), pick the right correlations for pressure drop and holdup (Beggs-Brill, Hagedorn-Brown, OLGA), and identify operational risks (severe slugging in a riser system, liquid carryover into a dry-gas compressor).

2. The five regimes

RegimeAppearanceDriver
Stratified SmoothGas on top, liquid on bottom; flat interfaceLow v_SG, low v_SL — gravity dominates
Stratified WavyGas shear creates surface wavesIncreased v_SG over stratified-smooth
Intermittent (Slug / Plug)Liquid slugs alternate with elongated gas bubblesWave growth bridges the pipe → slug
AnnularLiquid film on wall; gas core with entrained dropletsHigh v_SG strips liquid into core
Dispersed BubbleSmall gas bubbles in continuous liquidHigh v_SL — turbulence dominates buoyancy

3. Taitel-Dukler F/K/T groups

Taitel & Dukler (1976) showed that regime boundaries collapse onto a single chart when expressed in three dimensionless groups (plus the Lockhart-Martinelli X):

F = √(ρG / (ρL − ρG)) · vSG / √(D · g · cos θ)

F is a modified Froude number for the gas phase. The gravity term in the denominator makes it dimensionally a velocity ratio: the gas superficial velocity vs the gravitational wave-speed scale. F > ~0.5 says gas-phase inertia is overcoming gravity → stratified flow can't be maintained.

K = F · √(ReSL)

K extends F with a liquid-Reynolds-number factor that captures viscous wave amplification. The K = 20 boundary separates smooth stratified from wavy stratified.

T = √[ (dP/dL)L,sup / ( (ρL − ρG) · g · cos θ ) ]

T is the ratio of liquid turbulence (frictional pressure gradient) to buoyancy. T > ~1.5 means turbulence breaks bubbles into a fine dispersion → dispersed-bubble flow.

Together with the Lockhart-Martinelli X:

X2 = (dP/dL)L,sup / (dP/dL)G,sup

which compares the single-phase frictional pressure gradients each phase would have if it flowed alone in the pipe — X > 1.6 means liquid friction dominates → intermittent, X < 1.6 means gas friction dominates → annular.

4. Regime boundaries

The full Taitel-Dukler procedure requires iteratively solving for the equilibrium liquid level hL/D in stratified flow, then checking four stability criteria (one per regime transition). For screening, the simplified boundaries:

BoundaryApproximate criterion
Stratified → WavyK > 20
Stratified → IntermittentF > 0.5 (with X-dependent refinement)
Intermittent → AnnularX < 1.6
Intermittent → Dispersed BubbleT > 1.5
For final design, use a full code. OLGA, LedaFlow, and PIPESIM implement the full Taitel-Dukler procedure plus newer mechanistic models. This screening calc tells you which regime to expect; the commercial code tells you the pressure drop and holdup once you're committed.

5. Severe slugging

A special case worth its own screen: when stratified flow in a horizontal/downhill pipeline meets a vertical riser. Liquid accumulates at the base of the riser, eventually plugs the riser cross-section, and the gas pushes the entire column up at once. Then gas blows through, pressure crashes, liquid falls back, and the cycle repeats — slugs orders of magnitude longer than ordinary intermittent slugs.

Schmidt & Brill (1980) parameter for severe-slugging risk:

ΠSS = ρL · g · vSL · Lpipe / (P · vSG)

ΠSS > 1 indicates severe-slugging risk. Mitigations:

  • Backpressure on the riser top (operating cost — burns compressor HP).
  • Gas lift at the riser base (operating cost — sources of lift gas).
  • Slug-control valve (Texaco Slug Control or equivalent) at the riser top.
  • Slug catcher with 5–10× normal slug volume capacity to absorb the transient.

6. Design implications by regime

RegimeEquipment / design action
Stratified SmoothStandard separator at end; minimal upset risk; pig-frequency based on liquid line packing.
Stratified WavyAdd demister to receiving vessel; check for mist carryover into compressor suction.
Intermittent (Slug)Slug catcher sized for liquid-slug volume; downstream ESDV trip logic on slug catcher high-high level.
AnnularMist eliminator on receiving vessel; erosion check per API RP 14E if particulates; consider liquid drain at low points.
Dispersed BubbleTreat as pseudo-single-phase liquid; pressure drop dominated by liquid friction; minimal vessel design changes.

7. References

  • Taitel, Y. & Dukler, A. E. (1976). "A model for predicting flow regime transitions in horizontal and near horizontal gas-liquid flow." AIChE J. 22(1), 47–55.
  • Mandhane, J. M., Gregory, G. A., & Aziz, K. (1974). "A flow pattern map for gas-liquid flow in horizontal pipes." Int. J. Multiphase Flow. 1, 537–553.
  • Schmidt, Z., Brill, J. P., & Beggs, H. D. (1980). "Choking can eliminate severe pipeline slugging." Oil & Gas J. Nov 1980.
  • Beggs, H. D. & Brill, J. P. (1973). "A study of two-phase flow in inclined pipes." JPT, May 1973.
  • Zhang, H. Q., Sarica, C. (2003). "Unified Modeling of Gas-Liquid Pipe Flow." SPE 95749.
  • Brill, J. P. & Mukherjee, H. (1999). Multiphase Flow in Wells, SPE Monograph 17.
  • Shoham, O. (2006). Mechanistic Modeling of Gas-Liquid Two-Phase Flow in Pipes, SPE.

Map your flow regime

→ Open Regime Map Calculator