Gas Pipeline Sizing: Friction & Compressibility

In chemical processing plants, oil and gas gathering networks, natural gas transmission pipelines, and refinery gas treating units, sizing gas piping is fundamentally more complex than sizing liquid lines. While liquids maintain nearly constant density regardless of pressure drop, gases expand continuously as they traverse a conduit. A drop in pressure produces an immediate decrease in gas density ($\rho$), which forces the volumetric flow rate to expand and the linear velocity to accelerate toward the pipe exit.

Applying standard incompressible equations (such as the standard liquid Darcy-Weisbach or Hazen-Williams formulas) to gas lines when the pressure drop exceeds 10% of inlet pressure leads to catastrophic under-sizing. At high velocities, gas lines risk acoustic-induced vibration (AIV), intense noise exceeding OSHA thresholds, localized pipe erosion per API RP 14E, and thermal expansion shockwaves. At the extreme limit, the expanding gas accelerates until it reaches the local speed of sound ($Ma = 1.0$), resulting in choked flow where downstream pressure reductions cannot increase throughput.

This engineering reference covers the full mathematical theory of compressible gas flow in conduits: real gas equations of state ($Z$-factor), sonic velocity and Mach number criteria, the complete Isothermal Compressible Flow equation with kinetic acceleration, classical transmission formulas (Weymouth, Panhandle A & B), acoustic screening criteria ($\rho v^2$), and multi-schedule discrete optimization using the free ChemProCal Gas Line Sizer.

1. Why Gas Piping Differs Fundamentally from Incompressible Liquids

To understand gas line sizing, one must appreciate the coupled interaction between pressure, density, and velocity. Consider a gas flowing through a pipe of constant cross-sectional area $A$. By the continuity equation:

Because frictional shear against the pipe wall dissipates mechanical energy, the pressure must decrease along the length ($P_2 < P_1$). By the real gas equation of state:

As pressure falls, density $\rho$ must decrease proportionally. For mass flow to remain constant across a uniform cross-section, the linear velocity must increase:

This acceleration creates an additional pressure drop over and above pure wall friction: part of the pressure gradient must physically accelerate the expanding fluid mass.

2. Thermodynamics of Compressible Flow: Speed of Sound & Mach Number

The acoustic behavior of a gas establishes the upper physical speed limit for any piping system.

Speed of Sound in Real Gases ($c$)

The speed of sound ($c$) is the velocity at which infinitesimal pressure disturbances propagate through a compressible medium. For an ideal or real gas undergoing isentropic acoustic compression:

Where:

• $\gamma = c_p / c_v$ = Ratio of specific heats (isentropic expansion coefficient, typically $1.30 - 1.40$ for diatomic gases and dry natural gas, $1.15 - 1.25$ for heavy hydrocarbons)
• $P$ = Absolute static pressure ($\text{Pa}$)
• $\rho$ = Gas density ($\text{kg/m}^3$)
• $Z$ = Gas compressibility factor (dimensionless)
• $R_{univ}$ = Universal gas constant ($8,314.46\text{ J/(kmol}\cdot\text{K)}$)
• $M$ = Molecular weight ($\text{kg/kmol}$, e.g., $16.04\text{ kg/kmol}$ for methane, $28.96\text{ kg/kmol}$ for air)
• $T$ = Absolute temperature ($\text{K} = ^\circ\text{C} + 273.15$)

The Mach Number ($Ma$) and Flow Regimes

The Mach number ($Ma$) is the ratio of local fluid velocity to the local speed of sound:

Compressible pipe flow is classified into four distinct operational regimes:

3. The Isothermal Compressible Flow Equation

In process piping and long transmission lines, gas flowing through an uninsulated or underground pipe rapidly exchanges heat with the pipe wall and surroundings. For lines longer than roughly 50 to 100 meters, Isothermal Flow ($T \approx \text{constant}$) is the universally accepted standard model (as formalized in Crane Technical Paper 410, Section 1-8).

Derivation from Momentum and Energy Balances

Beginning with the differential one-dimensional momentum balance for compressible flow with wall friction in a conduit of constant cross-sectional area:

Expressing wall shear stress via the Darcy friction factor $\tau_w = \frac{f_D \rho v^2}{8}$, defining the mass flux $G = \frac{\dot{m}}{A} = \rho v$, and integrating between inlet ($P_1$) and outlet ($P_2$) at constant temperature $T$ yields the complete Fundamental Isothermal Flow Equation:

Where:

• $P_1, P_2$ = Absolute inlet and outlet pressures ($\text{Pa}$)
• $G = \frac{\dot{m}}{\pi D^2 / 4}$ = Mass flux or mass velocity ($\text{kg/(s}\cdot\text{m}^2\text{)}$)
• $\rho_1 = \frac{P_1 M}{Z R T}$ = Gas density evaluated at inlet conditions ($\text{kg/m}^3$)
• $f_D$ = Darcy friction factor evaluated via Colebrook-White or Swamee-Jain at the average Reynolds number
• $L / D$ = Pipe length to internal diameter ratio
• $\sum K$ = Total minor loss resistance coefficient for fittings and valves
• $2 \ln(P_1 / P_2)$ = Dimensionless kinetic energy acceleration term

Numerical Solution for Downstream Pressure ($P_2$)

Because $P_2$ appears both algebraically on the left side ($P_2^2$) and inside the natural logarithm on the right side ($\ln(P_1/P_2)$), this equation is nonlinear and implicit. ChemProCal’s engine defines the objective residual function:

The engine solves this root-finding problem using `scipy.optimize.fsolve`. If no physical solution exists where $P_2 > 0$ and $Ma_{exit} \le 1.0$, the engine immediately detects that the pipe is acoustically choked.

4. Fanno Flow & Sonic Choking Mechanics

To understand the physical limits of gas throughput, consider Fanno Flow: adiabatic flow with friction in a constant-area duct. As friction dissipates energy along the duct, the thermodynamic state of the gas moves along the Fanno line toward the point of maximum entropy—which corresponds precisely to Mach 1.0 (Sonic Velocity).

Critical Pressure Ratio ($P^* / P_0$)

When a gas accelerates to Mach 1.0 at an orifice or nozzle throat, the critical pressure ratio ($P^* / P_0$) is governed by the isentropic relationship:

For air and diatomic gases ($\gamma = 1.40$), $\frac{P^*}{P_0} = 0.528$. For natural gas ($\gamma \approx 1.30$), $\frac{P^*}{P_0} = 0.546$. If the receiver pressure downstream drops below this critical threshold, flow cannot accelerate any further inside the pipe. Instead, the remaining pressure difference expands through an under-expanded supersonic jet shock cell outside the pipe discharge nozzle.

5. Long-Distance Pipeline Equations (Weymouth, Panhandle A & B)

In cross-country gas transmission and field gathering networks spanning miles rather than meters, industry standards rely on specialized empirical integrations of the compressible flow equation. These formulations eliminate the need to solve friction factors iteratively by directly embedding empirical Reynolds-roughness relationships.

The General Flow Equation

All pipeline transmission equations derive from the steady-state isothermal integral:

Where $Q_b$ is flow rate at base conditions ($P_b, T_b$), $G_g$ is gas specific gravity ($M_{gas} / 28.96$), and $\Delta Z$ is elevation change.

6. Industry Design Criteria & Velocity Envelopes (API 14E & NORSOK P-002)

Sizing gas lines is a balance between keeping lines small enough to minimize steel tonnage and large enough to respect three operational boundaries: Erosional Velocity, Acoustic Vibration, and Permissible Pressure Drop.

1. API Recommended Practice 14E Erosional Velocity Ceiling

The American Petroleum Institute standard API RP 14E provides the universal criterion used worldwide in offshore and production piping:

Where $\rho$ is the gas density in $\text{kg/m}^3$ at operating conditions. In metric units, the constant is:

• Continuous Service (Solids-Free Gas): $C = 122$ ($100\text{ in US customary}$) $\implies v_e = \frac{122}{\sqrt{\rho}}$
• Intermittent Service (Relief / Flare / Blowdown): $C = 149$ ($122\text{ in US customary}$) $\implies v_e = \frac{149}{\sqrt{\rho}}$
• Corrosion-Resistant Alloys (Duplex / Super Duplex / 316L SS): $C = 183 - 244$ ($150 - 200\text{ in US customary}$) where sand is monitored.

2. NORSOK Standard P-002 Velocity & Pressure Drop Envelopes

NORSOK P-002 provides prescriptive sizing rules for process gas lines:

3. Kinetic Energy Flux ($\rho v^2$) & Acoustic Induced Vibration (AIV)

To prevent severe aerodynamic noise and fatigue failure of small-bore piping connections (thermowells, sample points, drain branches), industry standards (including Energy Institute guidelines and ISO 15649) screen the dynamic kinetic energy density ($\rho v^2$ in $\text{Pa}$):

7. The ChemProCal Smart Pipe Gas Optimization Algorithm

Manual calculation of compressible gas hydraulics requires guessing an exit pressure, computing the average density, looking up pipe schedules, solving the friction factor, evaluating the log term, and re-iterating until convergence.

The ChemProCal Gas Line Sizer automates this workflow:

1. Flexible Input Flow Modes: Accepts process throughput specified as Mass Flow ($\text{kg/h}$ or $\text{lb/h}$), Actual Volumetric Flow ($\text{m}^3\text{/h}$ or $\text{ACFM}$), or Standard Volumetric Flow ($\text{Nm}^3\text{/h}$ or $\text{MMSCFD}$).
2. ASME B36.10M / B36.19M Geometry Integration: Queries the `fluids` engineering core to retrieve exact internal diameters and thicknesses across all schedules (5 to XXS).
3. Nonlinear Numerical Solver (`fsolve`): Solves the full Isothermal Flow Equation with the $2 \ln(P_1/P_2)$ acceleration term for downstream pressure $P_2$ and exit velocity $v_{exit}$.
4. Acoustic & Choking Evaluation: Evaluates local speed of sound ($c = \sqrt{\gamma P_1 / \rho_1}$), computes inlet Mach ($Ma_{in}$) and exit Mach ($Ma_{exit}$), and flags choked lines where $Ma_{exit} \ge 1.0$.
5. Automated Optimizer Grid: Ranks standard nominal pipe sizes from NPS 1/2" to 24", testing each against:
   • $v \le v_{max}$ (Velocity constraint)
   • $v \le v_e = 122 / \sqrt{\rho}$ (API 14E erosional ceiling)
   • $\Delta P \le \Delta P_{max}$ (Pressure drop constraint)
   • $Ma_{exit} < 1.0$ (Non-choked flow)
6. Preferred Size Tagging: Highlights the most economical standard NPS pipe that satisfies all constraints without over-sizing.

8. Comprehensive Worked Engineering Example (Hand Calculation vs ChemProCal)

Let us solve a practical offshore natural gas export line design problem to demonstrate the calculation procedure.

Problem Statement: Natural Gas Gathering Header

An offshore production separator discharges dry natural gas into a gathering pipeline leading to a compression platform:

• Mass Flow Rate ($\dot{m}$): $25,000\text{ kg/h}$ ($6.944\text{ kg/s}$)
• Gas Properties: Molecular weight $M = 18.2\text{ kg/kmol}$ ($SG = 0.628$), ratio of specific heats $\gamma = 1.30$, dynamic viscosity $\mu = 0.012\text{ cP} = 1.2 \times 10^{-5}\text{ Pa}\cdot\text{s}$
• Operating Temperature ($T$): $35.0^\circ\text{C} = 308.15\text{ K}$
• Inlet Pressure ($P_1$): $35.0\text{ bar-a} = 3,500,000\text{ Pa}$
• Compressibility Factor ($Z$): $0.90$ at operating $P$ and $T$
• Pipeline Length ($L$): $250.0\text{ meters}$
• Pipe Material: Carbon steel, Schedule 40 ($\epsilon = 0.0457\text{ mm}$)
• Total Minor Loss Coefficient ($\sum K$): $2.50$ (elbows and isolation valves)
• Design Criteria: $v_{max} = 18.0\text{ m/s}$, $\Delta P_{max} = 1.50\text{ bar}$

Step 1: Baseline Thermodynamic Properties

Inlet gas density:

Local speed of sound at inlet:

API 14E erosional velocity limit:

Step 2: Evaluating Candidate Pipe Sizes

Option A: NPS 3" Schedule 40 ($D = 0.0779\text{ m}$)

Area: $A = \frac{\pi (0.0779)^2}{4} = 0.004766\text{ m}^2$

Inlet velocity: $v_1 = \frac{\dot{m}}{\rho_1 A} = \frac{6.944}{27.62 \times 0.004766} = \mathbf{52.75\text{ m/s}}$

Inlet Mach: $Ma_1 = \frac{52.75}{405.9} = 0.130$.

Evaluation: $v_1 = 52.75\text{ m/s}$ vastly exceeds $v_{max} = 18\text{ m/s}$ and exceeds the erosional ceiling ($23.21\text{ m/s}$). At this massive mass flux ($G = 1,457\text{ kg/m}^2\text{s}$), solving the isothermal equation reveals that the required pressure drop exceeds inlet pressure—the 3" pipe is heavily choked!

Option B: NPS 4" Schedule 40 ($D = 0.1023\text{ m}$)

Area: $A = \frac{\pi (0.1023)^2}{4} = 0.008219\text{ m}^2$

Mass flux: $G = \frac{6.944}{0.008219} = 844.9\text{ kg/(s}\cdot\text{m}^2\text{)}$

Inlet velocity: $v_1 = \frac{844.9}{27.62} = \mathbf{30.59\text{ m/s}}$ ($> 18.0\text{ m/s}$, exceeds allowable criteria and exceeds $v_e = 23.21\text{ m/s}$).

Reynolds number: $Re = \frac{27.62 \times 30.59 \times 0.1023}{1.2 \times 10^{-5}} = 7,205,000$. Friction factor $f_D = 0.0168$.
Solving the isothermal equation yields: $P_2 = 25.8\text{ bar-a} \implies \Delta P = 9.2\text{ bar} \gg 1.50\text{ bar}$ allowable limit.

Option C: NPS 6" Schedule 40 ($D = 0.1541\text{ m}$)

Area: $A = \frac{\pi (0.1541)^2}{4} = 0.01865\text{ m}^2$

Mass flux: $G = \frac{6.944}{0.01865} = 372.3\text{ kg/(s}\cdot\text{m}^2\text{)}$

Inlet velocity: $v_1 = \frac{372.3}{27.62} = \mathbf{13.48\text{ m/s}}$

Inlet Mach number: $Ma_1 = \frac{13.48}{405.9} = \mathbf{0.0332}$ (Well within subsonic regime!).

Velocity evaluation: $13.48\text{ m/s} \le 18.0\text{ m/s}$ (Pass!) and $13.48\text{ m/s} < 23.21\text{ m/s}$ (Erosional check pass!).

Reynolds number:

Relative roughness: $\frac{\epsilon}{D} = \frac{0.0000457}{0.1541} = 0.000297$. Friction factor from Swamee-Jain: $f_D = 0.0154$.

Equivalent friction factor terms:

Substitute into the Isothermal Equation:

Solving iteratively for $P_2$:

Exit velocity: $v_{exit} = 13.48 \times \left(\frac{35.00}{34.30}\right) = \mathbf{13.75\text{ m/s}}$. Exit Mach: $Ma_{exit} = 0.034$.

Verdict on NPS 6": 🟢 PREFERRED SIZE — Velocity ($13.48\text{ m/s}$) and pressure drop ($0.70\text{ bar}$) both satisfy criteria with an economical carbon steel specification.

Optimizer Comparison Grid from ChemProCal Gas Sizer

9. Flow-Induced Vibration (FIV) & Acoustic Fatigue (AIV)

In high-velocity gas lines and downstream of pressure let-down valves (such as Joule-Thomson refrigeration valves or turbine bypasses), acoustic vibrations can generate severe pipe failure within hours if not screened early.

Carucci-Mueller Acoustic Screening

The Carucci-Mueller methodology calculates the total Sound Power Level ($PWL$) generated inside the fluid stream:

Where $W$ is mass flow rate ($\text{lb/hr}$ or $\text{kg/s}$) and $\Delta P$ is pressure drop across the restriction.

10. Software Comparison: ChemProCal Gas Sizer vs. Generic Tools

11. Frequently Asked Questions (FAQ)

Why does gas velocity increase toward the outlet of a pipe?

As a gas flows through a pipe, friction dissipates mechanical energy, causing pressure to drop ($P_2 < P_1$). Because gases are compressible, a lower pressure means the molecules spread farther apart, reducing gas density ($\rho = PM/ZRT$). Since the total mass flow rate ($\dot{m} = \rho A v$) must remain constant by the law of conservation of mass, the gas must accelerate to higher velocities toward the outlet.

What happens if a gas line is undersized?

An undersized gas line causes excessive pressure drop that starves downstream equipment, creates extreme erosive shear on pipe walls per API 14E, generates deafening aerodynamic noise, induces acoustic fatigue on pipe branches, and can trigger sonic choking ($Ma = 1.0$) where throughput is capped regardless of how low downstream pressure falls.

How does temperature affect gas pressure drop compared to liquids?

In liquids, higher temperature decreases viscosity and reduces pressure drop. In gases, higher temperature has the opposite effect:

1. Gas dynamic viscosity ($\mu$) increases with temperature because molecular thermal agitation increases momentum exchange.
2. Gas density ($\rho$) decreases with temperature at constant pressure, increasing volumetric flow rate and velocity ($v = \dot{m} / \rho A$).

Therefore, heating a gas increases frictional pressure drop, whereas heating a liquid decreases pressure drop.

Can the Mach number exceed 1.0 inside a straight pipe?

No. In a constant-cross-section pipe with friction (Fanno flow), subsonic flow can only accelerate up to a maximum of Mach 1.0 at the exit plane. Supersonic flow ($Ma > 1.0$) inside a conduit can only be achieved by accelerating gas through a converging-diverging (de Laval) nozzle. If you attempt to push more mass through a choked pipe, the upstream pressure simply builds up until the inlet density increases to accommodate the flow.

What is the difference between Actual Flow and Standard Flow?

Actual Flow Rate (ACFM or $\text{m}^3\text{/h}$) is the physical volume occupied by the gas at its exact operating pressure and temperature inside the pipe. Standard Flow Rate (SCFM, MMSCFD, or $\text{Nm}^3\text{/h}$) is a mass-equivalent volumetric flow referenced to standardized base conditions (typically $1.01325\text{ bar-a}$ and $15^\circ\text{C}$ or $60^\circ\text{F}$). Standard flow remains constant throughout a process regardless of pressure, whereas actual flow changes at every point along the pipeline.