1. The Joukowsky Surge
When flowing liquid is stopped suddenly, its momentum converts to a pressure spike that travels as a wave — "water hammer." For instantaneous flow stoppage the peak pressure rise is the Joukowsky equation:
where ρ is the fluid density, a the pressure-wave speed, and ΔV the change in flow velocity. This is the instantaneous single-wave surge — the bound for a short or essentially frictionless line. It is not a universal upper limit: in a long, high-friction liquid pipeline, line packing drives the peak pressure at the closing valve above the single-wave Joukowsky value, because the pressure keeps rising behind the wavefront as flow is still being arrested upstream (see §4). A 2 m/s liquid line with a ≈ 1,200 m/s and ρ ≈ 800 kg/m³ gives ΔP ≈ 1.9 MPa (~280 psi) from a single fast closure — and the packed peak in a long line can exceed this.
2. Wave Speed
The wave speed is set by the fluid's compressibility and the pipe's elasticity (the pipe "gives" slightly, slowing the wave):
where K is the fluid bulk modulus, E the pipe Young's modulus, D/t the diameter-to-wall ratio, and c₁ a restraint factor. Stiffer pipe (lower D/t, higher E) → faster wave → larger Joukowsky surge. Entrained gas dramatically lowers K and a, which both reduces surge and complicates the analysis.
3. Critical Closure Time
The full Joukowsky surge only develops if the valve closes faster than the wave can travel to the far end and back — the critical (pipe-period) time:
Closing in less than tc ("rapid closure") gives the full ΔP; closing slower ("gradual closure") gives a reduced surge roughly proportional to tc/tclose. This is exactly why valve closure rates and pump-trip ramps are engineered — extending the effective closure time is the primary surge-mitigation lever, alongside relief/surge tanks and slow-opening check valves.
4. The Method of Characteristics
Real surge events — pump trips, valve sequences, reflections off tees and elevation changes — need the full transient. The method of characteristics (MOC) solves the coupled continuity and momentum partial differential equations along characteristic lines (dx/dt = ±a), marching pressure and flow forward in time on a space-time grid. It captures wave reflection and superposition, line friction, and boundary conditions (valves, pumps, tanks) that a single Joukowsky number cannot. Crucially, in a long, friction-dominated liquid line MOC reproduces line packing: because the steady friction gradient is recovered as flow decelerates, the head at the closing valve continues to climb behind the initial wavefront and the peak can rise above the single-wave Joukowsky ΔP. This is the Wylie & Streeter result, and it is exactly why the surge calculator reports a separate line-packing term that can push the computed peak past ρ·a·ΔV.
5. Line Pack (Inventory)
Line pack is the mass of gas stored in a pipeline at operating pressure — the buffer that lets a gas system absorb a mismatch between supply and demand. It is computed from first principles as the real-gas inventory of the pipe volume:
where V is the internal pipe volume, M the gas molar mass, T the flowing temperature, and Z the real-gas compressibility factor evaluated by the AGA Report No. 8 / GERG equation of state at the local pressure and temperature. Because pressure varies along the line, an accurate line-pack uses the integrated pressure profile rather than a single average. There is no "AGA Report No. 17" line-pack method — line pack is a thermodynamic inventory calculation, and the relevant AGA reference is Report No. 8 for the Z-factor; API RP 1130 covers leak detection, not line pack. Line pack underpins pack/draft operations, nomination balancing, and the pressure-decay timing of a shutdown.
6. Code Surge Limits
For hazardous-liquid pipelines, ASME B31.4 limits the surge pressure such that the sum of steady operating pressure and surge does not exceed 110% of the internal design (MAOP) at any point — so the surge analysis is a code-compliance calculation, not just an operability check. Gas transmission transients are governed by ASME B31.8. The shut-in/static-head case must also be checked: in liquids the basis is 1.10 × max(steady, shut-in static head).
7. References
- Joukowsky, N. (1898/1900) — water-hammer theory; Wylie & Streeter, Fluid Transients in Systems (MOC).
- AGA Report No. 8 (AGA-8 / GERG) — Compressibility Factor of Natural Gas (Z for line-pack inventory).
- ASME B31.4 — Pipeline Transportation Systems for Liquids and Slurries (surge ≤ 110% design).
- ASME B31.8 — Gas Transmission and Distribution Piping Systems.
- API RP 1130 — Computational Pipeline Monitoring for Liquids (leak detection; not line pack).
Ready to use the calculator?
→ Launch Surge Analysis Calculator