1. Introduction
Water hammer — also referred to as hydraulic transient or pressure surge — occurs in pressurised pipeline systems whenever the rate of flow changes rapidly. The kinetic energy of the moving fluid column converts instantaneously into pressure energy, generating a wave that propagates at the acoustic wave speed. Left unmanaged, this surge pressure can exceed the static working pressure by several multiples, leading to pipe rupture, joint failure, pump damage, or structural collapse.
In practice, the hydraulic engineer has relatively few tools to actively mitigate water hammer once pipeline geometry and material are fixed. The control valve, however, offers a proactive design lever: by extending closure time, the designer can dramatically reduce peak surge pressure without additional infrastructure. This article examines that relationship in full mathematical and practical detail, and links the analytical approach directly to industry-standard numerical simulation using Bentley OpenFlows HAMMER.
2. Physics of Water Hammer
When a control valve begins to close, the cross-sectional area available to flow decreases. The fluid immediately upstream decelerates, generating a pressure rise that propagates upstream as a compression wave at the acoustic (pressure wave) speed a. The full cycle of wave travel — from valve to the upstream reservoir and back — defines the critical period. If closure completes within this period, the reflected relief wave arrives too late, and the full instantaneous pressure rise is realised.
3. The Joukowsky Equation — Instantaneous Closure
For an instantaneous and complete valve closure, the maximum pressure rise is given by the Joukowsky equation:
a = Acoustic wave speed (m/s) | ΔV = Velocity change = V₀ for full closure
The acoustic wave speed depends on both fluid compressibility and pipe-wall elasticity (Korteweg formula):
| Pipe Material | E (GPa) | Typical a (m/s) | Typical Application |
|---|---|---|---|
| Rigid (theoretical) | ∞ | 1,450 | Reference case |
| Steel (Carbon) | 206 | 900–1,200 | Transmission mains |
| Ductile Iron | 170 | 1,000–1,250 | Distribution, WTPs |
| GRP / FRP | 20–40 | 400–700 | Gravity mains |
| HDPE (PE100) | 0.8–1.0 | 200–400 | Water supply, effluent |
| PVC-U | 2.7–3.5 | 250–450 | Low-pressure distribution |
4. Critical Time Tc — The Design Threshold
Rapid Closure
Full Joukowsky pressure always develops — regardless of how much faster the valve closes. This is the worst-case design scenario.
Slow (Controlled) Closure
Peak pressure is attenuated below the Joukowsky value. Longer closure relative to Tc = lower surge. This is the target design condition.
5. Slow-Closure Attenuation Formula
- Compute acoustic wave speed for pipe material and geometry: a = √[ K/ρ / (1 + (K/E)(D/e)) ]
- Calculate the pipeline critical period: Tc = 2L / a
- Compute the Joukowsky upper-bound pressure rise: ΔPJ = ρ · a · V₀
- Compare closure time against Tc: ΔP = ΔPJ · min(1, Tc / tv)
- Check maximum system pressure against pipe PN rating: Pmax = Pstatic + ΔP ≤ 1.5 × PN
6. Interactive Surge Pressure Calculator
Adjust the pipeline parameters below to instantly compute the critical time, Joukowsky pressure, and the effect of different closure times.
ΔP vs. Closure Time for current system. Dashed line = Tc. Red dot = current tv.
Idealised pressure-time trace at the valve. Blue = static pressure. Red = surge envelope. Uses current calculator parameters.
Joukowsky pressure (bar) for five pipe materials at the same V₀ and L from the slider above.
7. Numerical Example — Multi-Scenario Comparison
The following example is representative of a steel water transmission main. Base parameters: pipeline length L = 2,000 m, diameter DN 900 mm, wave speed a = 1,000 m/s (steel), flow velocity V₀ = 1.5 m/s, static pressure P₀ = 8.0 bar, pipe rating PN 16.
| Closure Time | tv/Tc | Regime | Factor | ΔP (bar) | Pmax (bar) | PN Needed | Status |
|---|---|---|---|---|---|---|---|
| 1.0 s | 0.25 | Rapid | 1.00 | 15.0 | 23.0 | PN 25 | Unacceptable |
| 2.0 s | 0.50 | Rapid | 1.00 | 15.0 | 23.0 | PN 25 | Unacceptable |
| 4.0 s = Tc | 1.00 | Boundary | 1.00 | 15.0 | 23.0 | PN 25 | Unacceptable |
| 8.0 s | 2.00 | Slow | 0.50 | 7.5 | 15.5 | PN 16 | Marginal |
| 20.0 s | 5.00 | Slow | 0.20 | 3.0 | 11.0 | PN 12.5 | Acceptable |
| 40.0 s | 10.0 | Slow | 0.10 | 1.5 | 9.5 | PN 10 | Good |
| 80.0 s | 20.0 | Slow | 0.05 | 0.75 | 8.75 | PN 10 | Excellent |
| 120.0 s | 30.0 | Slow | 0.033 | 0.50 | 8.50 | PN 10 | Excellent |
8. Bentley OpenFlows HAMMER — Modelling Valve Closure
While the analytical formulas provide invaluable first-pass design estimates, final design validation for large-scale water infrastructure requires full numerical simulation using the Method of Characteristics (MOC). Bentley OpenFlows HAMMER is the industry-standard software for hydraulic transient analysis, referenced in major project specifications including those of NWC (National Water Company, Saudi Arabia).
8.1 The Method of Characteristics (MOC)
HAMMER solves the two governing partial differential equations of unsteady pipe flow simultaneously:
These are transformed into characteristic equations along the C+ and C− characteristic lines, enabling the propagation of pressure waves through the system to be tracked at each time step Δt = Δx/a. This is fundamentally more accurate than simplified analytical formulas for complex systems with branching, multiple valves, pumps, and surge protection devices.
8.2 Setting Up Valve Closure in HAMMER
8.3 Key HAMMER Parameters
| Parameter | Effect on Result | Notes |
|---|---|---|
| Wave Speed (a) | Directly scales both ΔPJ and Tc | Auto-computed from material/geometry. Verify against Korteweg formula for non-standard materials. |
| Closure Time (tv) | Primary control parameter | Must be specified as total travel time. For two-speed actuators, use the custom closure curve instead. |
| Valve Characteristic Curve | Determines how Cv changes with stem position over time | Built-in options: Gate, Globe, Butterfly, Ball, Check. Custom import via CSV. This curve — not just tv — determines the effective flow-reduction profile. |
| Time Step (Δt) | Numerical accuracy of wave propagation | Must satisfy Courant condition: Δt = Δx/a. HAMMER enforces this automatically — do not override. |
| Friction Model | Attenuates pressure waves over time | Options: Steady, Quasi-steady, or Unsteady (Vardy-Brown). For preliminary design, steady friction is acceptable. |
| Max / Min Pressure Limits | Design validation output | Minimum pressure below vapour pressure (≈ −10 m) indicates column separation — activate DVCM (Discrete Vapour Cavity Model) in HAMMER. |
8.4 When Do Analytical and HAMMER Results Diverge?
| Condition | Analytical Formula | HAMMER (MOC) |
|---|---|---|
| Simple series pipe, rapid closure | Accurate (within 5%) | Accurate — use as verification check |
| Simple series pipe, slow closure | Approximate (valid up to ~5 × Tc) | Accurate for any closure time |
| Non-linear valve characteristic | Not applicable — assumes linear flow reduction | Accurately models any Cv curve |
| Pipeline with high points | Cannot predict column separation | Models column separation with DVCM |
| Multi-branch networks | Not applicable | Handles all network topologies |
| Multiple transient events | Cannot superimpose | Simulates all events simultaneously |
| Surge protection devices | Cannot model device interaction | Full modelling of all device types |
8.5 Reading the HAMMER Pressure Envelope
- Maximum pressure (red envelope) must not exceed 1.5 × PN rating of the pipe at any point along the pipeline. Identify the location of peak pressure — typically immediately upstream of the valve, but reflections can amplify pressure at intermediate high points.
- Minimum pressure (blue envelope) must remain above −10 m (vapour pressure at 20°C) at all points. Pressure falling below this at a high point indicates column separation — potentially more destructive than the positive surge wave when the columns rejoin.
- Pressure class transitions along the pipeline (e.g., PN 16 to PN 10 sections) must be identified on the profile overlay. Ensure the maximum pressure envelope does not exceed the PN rating of each pipe section independently.
- Time-history plots at the valve node and at critical intermediate nodes confirm the duration of the surge event and whether pressure oscillations are adequately damped by friction before the next transient event.
9. Practical Design Recommendations
| Parameter | Recommendation | Basis |
|---|---|---|
| Minimum closure time | Target tv ≥ 5 × Tc as starting point | At 5 × Tc, ΔP = 20% of Joukowsky — manageable within standard PN margin |
| Closure time from PN rating | tv = ΔPJ · Tc / (PNallowable − P₀) | Directly solves for minimum closure time without over-engineering |
| Two-speed actuators | Fast initial (0→80% travel) then slow final (20% over ≥ 5 × Tc) | Most flow restriction occurs near full closure — only the last 20% of stroke is hydraulically critical |
| Valve characteristic curve | Always use manufacturer Cv curve in HAMMER — not assumed linear | Butterfly valves are highly non-linear: ~70% of flow passes through the last 30% of travel |
| Power failure scenario | Evaluate spring-return closure time separately from motorised closure time | Fail-safe spring actuators typically close 3–5× faster than motorised operation |
| MOC simulation trigger | Require HAMMER for L > 1,000 m or any system with high points, branches, or surge devices | Analytical formulas lose accuracy; only MOC resolves reflections, column separation, and device interaction |
| Maximum pressure limit | Pmax ≤ 1.5 × PN (transient allowance per EN 805 / AWWA M11) | Temporary surge pressure permitted up to 1.5× PN for infrequent events |
| Minimum pressure limit | Pmin > −10 m at all pipeline high points | Column separation below vapour pressure causes rejoining impact waves exceeding the positive surge in severity |
10. Summary and Conclusions
- The Joukowsky equation (ΔP = ρ·a·ΔV) defines the absolute upper limit of surge pressure for any closure equal to or faster than the critical period.
- The critical time Tc = 2L/a is the decisive design threshold. Closure time must be compared against Tc, not evaluated in absolute seconds alone.
- For slow closures (tv > Tc), surge pressure is inversely proportional to closure time — doubling the closure time halves the surge pressure.
- Pipe material has a profound effect: steel pipes generate surge pressures 3–4× higher than HDPE under identical hydraulic conditions, making closure time specification especially critical for metallic transmission mains.
- Bentley OpenFlows HAMMER (MOC-based) is the appropriate tool for final design validation — it captures non-linear valve characteristics, pipeline branches, surge device interaction, and column separation.
- The recommended workflow: analytical formulas for preliminary design → interactive verification → HAMMER MOC for final design validation and regulatory submission.
References
- Joukowsky, N. (1898). Über den hydraulischen Stoss in Wasserleitungsröhren. Mémoires de l'Académie Impériale des Sciences de St.-Pétersbourg, 9(5). [ASME Translation, 1904]
- Streeter, V. L., & Wylie, E. B. (1993). Fluid Mechanics (9th ed.). McGraw-Hill.
- Chaudhry, M. H. (2014). Applied Hydraulic Transients (3rd ed.). Springer. ISBN 978-1-4614-8538-4.
- American Water Works Association (2004). AWWA Manual M11 — Steel Pipe: A Guide for Design and Installation (4th ed.). AWWA, Denver.
- European Committee for Standardisation (2000). EN 805:2000 — Water Supply: Requirements for Systems Outside Buildings. BSI, London.
- Thorley, A. R. D. (2004). Fluid Transients in Pipeline Systems (2nd ed.). Professional Engineering Publishing, London.
- Fisher Controls International / Emerson Process Management (2005). Control Valve Handbook (4th ed.). Emerson.
- Bentley Systems, Incorporated (2023). OpenFlows HAMMER User Guide — Transient Analysis Reference. Bentley Documentation Portal.
- Plastic Pipe Institute (2012). Handbook of PE Pipe (2nd ed.). PPI, Irving TX.
- National Water Company — Saudi Arabia / NWC (2020 ed.). Technical Specifications for Water Transmission Networks. NWC, Riyadh.