When the power fails, a pump does not stop at once. It coasts down on the energy in its rotating parts, and while it coasts it keeps the column moving. Make that rundown long enough and the downsurge never reaches vapour. On our 12 km reference main the motor and pump alone, at about 25 kg·m², protect nothing. Sixteen times that, 400 kg·m² in total, keeps every point on the line at or above +9.3 m (+7.3 m under the least favourable torque law tried), with no upsurge above the steady 85 m. The price comes at every start: with half of rated torque available for acceleration, the set takes 28.1 s to reach speed instead of 1.8.

1 · Rotating energy as surge protection

The load case that usually sizes the protection is a power failure that trips every pump at once (see Surge Scenarios in Pumping Stations). The reference main for this series is a 12 km, DN800 ductile iron K9 line [1] on a flat profile. It carries 2,520 m³/h (1.39 m/s) from an HGL of 85.0 m at the pump to 44.8 m at the reservoir, with a wave speed of 1,050 m/s. The pipe is PN16 (136 m allowable) [2], and the design minimum for this series is +3.0 m. Unprotected, with the pump modelled as stopping at once, the downsurge reaches vapour (−9.8 m) along the line, and the peak when the cavities collapse is 165–190 m. That is above PN16 whichever cavity model you believe.

A hydropneumatic vessel, a one-way surge tank or air valves put water or air back into the line. A flywheel goes after the cause instead: the column is being stopped too quickly. Until the first reflection returns from the reservoir, each cubic metre per second the pump stops delivering drops the head at the pump by

\[ \Delta H \;=\; \frac{a}{gA}\,\Delta Q \;=\; \frac{1050}{9.81 \times 0.5027}\,\Delta Q \;=\; 212.9\ \text{m per m}^3\text{/s} \]

So losing all 0.70 m³/s within 2L/a = 22.9 s would call for a 149.1 m drop, from a head of 85 m that can fall only 94.8 m before it reaches vapour [3]. A pump with enough rotating mass is still delivering when the reflections come back, and the head never reaches vapour. Nothing is added to the line, there is no gas to maintain and nothing to switch on: the protection acts on every trip [4].

What a flywheel does, and what it does not It slows the rate at which the flow falls, which tackles the downsurge at its source. It does nothing for an upsurge caused elsewhere, it cannot help a high point that separates for its own reasons, and it makes every normal start slower and harder on the motor. In our practice, that last point decides more flywheel designs than the hydraulics does.

2 · The rundown equation: stored energy, τ and GD²

Once the motor loses its supply, the only thing driving the pump is the kinetic energy of everything on the shaft: impeller, rotor, coupling and flywheel [3][5]:

\[ I\,\frac{d\omega}{dt} \;=\; -\,T, \qquad T \;=\; \frac{\rho g\,Q\,H_p}{\eta\,\omega} \;+\; T_{loss} \]

Here \(I\) is the total inertia on the shaft (kg·m²), \(\omega\) the speed, \(Q\) and \(H_p\) the flow and head the pump is still producing, \(\eta\) its efficiency and \(T_{loss}\) the windage and friction torque. Pump head falls with the square of speed. On the model pump the zero-flow head is \(100\,w^2\) metres above a 5 m suction, where \(w\) is the speed ratio. So once the speed has dropped by about 11 % (\(5 + 100\,w^2 = 85\) at \(w\) = 0.89), the pump cannot hold 85 m even at zero flow. A flywheel does not keep the pressure up. It keeps the flow going while the pressure falls.

\[ E_k \;=\; \tfrac12\,I\,\omega_0^2, \qquad \tau \;=\; \frac{E_k}{P_0}, \qquad GD^2 \;=\; 4I \]

\(E_k\) is the energy stored at running speed \(\omega_0\), and \(\tau\) is how long that energy could supply the rated power \(P_0\). At a constant rated torque \(T_0 = P_0/\omega_0\), the set would stop in \(I\omega_0/T_0 = 2\tau\). In reality the torque falls with speed and the rundown tails off more slowly, but \(\tau\) still puts sets of any speed and power on one scale [6]. Parmakian's classical charts for power failure at a pump rest on the same idea, a parameter that sets the rotor's stored energy against the pipeline's wave travel time [7].

Many data sheets give GD², which in kg·m² is numerically \(4I\). US data sheets give WR² in lb·ft², which is \(I\) itself; multiply by 0.0421 to get kg·m². A mix-up here is a factor of four in the input that matters most.

3 · The model, and six inertias on the reference main

The numbers come from a method-of-characteristics model with a vapour cavity model, cross-checked with a gas cavity model and an independent second code [3][8]. The pump is a screening model. None of this is Bentley HAMMER output: a project has to be run in HAMMER, or an equivalent program, on the real profile with the real pump data.

HAMMER describes the pump by four-quadrant characteristics [5][10], which set the torque during the rundown. The simple torque law used here is the main uncertainty, and section 5 puts a number on it. There is also a quirk: once its head rise has fallen to zero, the slowing model pump still passes water from its 5 m suction but takes no hydraulic torque. The bare set's speed therefore levels off at about 0.2 instead of reaching zero, and as \(I\) falls towards zero the model does not become exactly an instantaneous stop. Read that part of the speed trace as a feature of the torque law, not a prediction.

What each inertia stores: 685 kW at 1,480 rpm, 2L/a = 22.9 s, start time with 50 % of rated torque available for acceleration (section 8)
Total I (kg·m²)GD² (kg·m²)Ek (MJ)τ (s)τ / (2L/a)Start time (s)
25 (bare set)1000.300.440.021.8
1004001.201.750.087.0
2008002.403.510.1514.0
4001,6004.807.010.3128.1
8003,2009.6114.030.6156.1
1,6006,40019.2228.051.23112.2
What each inertia does to the line: all pumps trip, no other protection, pressure heads above the pipe
Total I (kg·m²)Lowest head on the line (m)WhereLowest head at the pump (m)Highest head on the line (m)Speed ratio at 5 / 10 / 22.9 s
25−9.8 (vapour)7.4–11.3 km+2.5not quoted (cavitating)0.22 / 0.22 / 0.20
100−8.77.8 km−1.388.2 (indicative)0.40 / 0.27 / 0.25
200−1.55.6 km+5.185.00.58 / 0.40 / 0.28
400+9.33.0 km+10.885.00.73 / 0.57 / 0.37
800+17.8at the pump+17.885.00.85 / 0.73 / 0.54
1,600+26.0at the pump+26.085.00.92 / 0.85 / 0.70
Why one maximum is not quoted With 25 kg·m² cavities form and collapse between 7.4 and 11.3 km. The vapour cavity run shows nothing above the steady 85 m, but collapse peaks depend on how the cavity is represented [8] (across this series the gas cavity model puts them between 18 % lower and 10 % higher than the vapour cavity model), and this case was not rerun with the gas cavity model, so no maximum is quoted. The 165–190 m of section 1 does not apply either: it belongs to the pump stopping at once. The design answer is to keep the line out of vapour rather than argue about the peak. Every case from 200 kg·m² up stays clear of vapour under every torque law tried, so those results are quoted to 0.1 m.

At 100 kg·m² the line stays only 1.1 m clear of vapour, with a small upsurge to 88.2 m at the pump. Under the least favourable torque law tried (section 5) it reaches vapour, so treat both figures as indicative. With a short rundown the lowest point is well down the line; with a long one it moves back to the pump.

4 · Interactive: the rundown

Pick an inertia and watch the pump slow down and the head at the pump respond. The grey traces are the bare motor and pump.

Pump speed and pump-end pressure head after a power failure
Method-of-characteristics model with a vapour cavity model, 12 km DN800, a = 1,050 m/s, 685 kW at 1,480 rpm, simplified homologous rundown with an ideal non-return valve. Traces are thinned to 260 points for the page; the readouts come from the full run.
Motor rotor, impeller, coupling and flywheel together.
The vertical line is 2L/a = 22.9 s, when the first reflection returns from the reservoir.
Inertia I
400 kg·m²
GD²
1,600 kg·m²
Stored energy
4.80 MJ
τ = E/P
7.01 s
Lowest head on line
+9.3 m, at 3.0 km
Verdict
meets +3.0 m
Highest head on line
85.0 m

At 400 kg·m² the speed ratio is 0.73 at 5 s, 0.57 at 10 s and 0.37 at 2L/a. The pump-end head bottoms out at +10.8 m just as the first reflection returns, the line never drops below +9.3 m, and nothing rises above 85 m. Switch to 25 kg·m²: the speed is down to 0.22 within 5 s (where this model holds it; see section 3) and the pump-end head is below 10 m within two seconds. The pump end holds at a few metres only because the slowing pump still passes water from its suction; meanwhile the far end of the line goes to vapour. Now try 200 kg·m². The pump-end trace never falls below +5 m, yet the line fails at 5.6 km. Judge the design on the minimum envelope along the whole line, not on the history at the pump.

5 · Worked example: how much inertia the 12 km main needs

This is the sequence we follow on a project, run here on the reference main.

Step 1: the duty and the set

\[ P_0 = \frac{\rho g Q H}{\eta} = \frac{998 \times 9.81 \times 0.70 \times 80}{0.80} = 685\ \text{kW}, \qquad \omega_0 = \frac{2\pi \times 1480}{60} = 155.0\ \text{rad/s}, \qquad T_0 = \frac{P_0}{\omega_0} = 4{,}420\ \text{N m} \]

Step 2: test the bare set

With 25 kg·m², \(E_k\) = ½ × 25 × 155.0² = 0.30 MJ and \(\tau\) = 0.44 s, a ratio of 0.02 to 2L/a. The speed is down to 0.22 within five seconds, and the line reaches vapour over about 4 km, from 7.4 to 11.3 km. The motor and pump alone protect nothing.

Step 3: double the inertia until the whole line passes

The response is far from linear. At 100 kg·m² the lowest head is −8.7 m at 7.8 km. At 200 the pump reads +5.1 m, but the line reaches −1.5 m at 5.6 km, with about 7 km of it below +3.0 m. At 400 the whole line holds at +9.3 m or better, with the lowest point at 3.0 km, and nothing rises above 85 m. That is 16 times the bare set. The true threshold lies somewhere between 200 and 400; we design at the tested value that passes with margin.

Step 4: check that the answer survives the torque law

The weakest assumption is the torque law, so the audit reran the same line with other plausible torque characteristics. The tables above use the law that gives the highest minima of those tried. The minima moved by 2–3 m: about −4 to −1.5 m at 200 kg·m² (failing under every law tried), +7.3 to +9.3 m at 400 (passing under every law tried) and +24 to +26 m at 1,600. The conclusion is robust; the decimals are not. At worst, 400 kg·m² still clears +3.0 m by more than 4 m.

Step 5: check the pressure side, and compare

The maximum is 85.0 m: no upsurge at all, and far inside 136 m. The reference vessel used across this site (20 m³ shell, 3.5 m³ gas, DN400 differential connection) gives +4.3 m and 119.2 m on the same line with the pump stopped at once, so the comparison is not strictly like for like; even at +7.3 m the flywheel keeps more margin on the downsurge. Cost, space and operation are compared in Choosing surge protection on one pipeline.

Step 6: turn the inertia into steel

The flywheel must add 375 kg·m², which stores 4.50 MJ. For a solid steel disc \(I = \tfrac12\rho\pi t\,r^4\), so at \(t\) = 120 mm:

\[ r = \left(\frac{2I}{\rho\pi t}\right)^{1/4} = \left(\frac{2 \times 375}{7850\,\pi \times 0.12}\right)^{1/4} = 0.7095\ \text{m} \;\rightarrow\; 0.71\ \text{m} \]

At 0.71 m the disc is 1.42 m across and weighs 1,492 kg. It adds 376 kg·m², for a total of 401 kg·m², and its rim runs at 110.0 m/s.

Step 7: price the start

With 50 % of rated torque available for acceleration (section 8), the start time rises from 1.8 s to 28.1 s. That figure goes to the motor manufacturer before anything else is fixed.

τ against 2L/a: a screening indicator, not a rule

Protection arrived at τ ≈ 7 s, about a third of 2L/a (0.31); at 0.15 it failed. Do not turn that into a rule. The ratio a line needs depends on its profile, friction, pump curve, torque characteristic and criterion, and a knee can separate whatever the ratio. Use it only to decide whether a flywheel is worth modelling at all [6][7].

6 · Interactive: inertia against minimum pressure

These are the six runs: minima on the line and at the pump against total inertia, with τ on the right-hand axis, and below them the minimum envelope along the main for the selected case. Set your own criterion.

How much inertia is enough on the 12 km main
Six model runs, 685 kW at 1,480 rpm, all pumps tripping. Upper: line and pump-end minima against total inertia (log scale). Lower: minimum pressure head along the line for the selected inertia, with the steady HGL; the thin grey lines are the other five inertias. Length below the criterion is counted on the 100 m computational grid.
Highlights the case above and draws its envelope below.
+3.0 m is the reference criterion. Owners and pipe classes differ.
Selected I
400 kg·m²
τ / (2L/a)
0.31
Lowest on line
+9.3 m
Where
at 3.0 km
Lowest at pump
+10.8 m
Line below criterion
0.0 km
Verdict
meets +3.0 m
Smallest tested I that passes
400 kg·m²

At +3.0 m the smallest tested inertia that passes is 400 kg·m², with τ/(2L/a) = 0.31. Select 200 and the two minima split apart: +5.1 m at the pump but −1.5 m at 5.6 km, with about 7 km of line below the criterion. Lower the criterion to 0 m and 400 is still needed; at −1.5 m, 200 passes. Raise it to +10 m and the answer doubles to 800. Treat any pass by less than 2–3 m with suspicion, since that is the spread from the torque law alone. Up to 100 kg·m² the line is at or near vapour. From there each doubling of inertia raises the line minimum by 7–11 m (8.2 m from 800 to 1,600), but it also doubles the start time, so every extra metre of margin costs more seconds at every start.

7 · Designing the flywheel

For a solid steel disc (\(\rho\) = 7,850 kg/m³):

\[ m = \rho\pi r^2 t, \qquad I = \tfrac12 m r^2 = \tfrac12\rho\pi t\,r^4, \qquad \frac{E_k}{m} = \frac{(\omega_0 r)^2}{4} = \frac{v_{rim}^2}{4} \]

Inertia goes with the fourth power of radius. Doubling the thickness doubles both inertia and mass, while adding 19 % to the radius doubles the inertia for only 41 % more mass. The energy per kilogram depends only on rim speed. Rim speed sets the stress in the disc, the hub and the shaft fixing, so it is what limits the diameter, and that limit comes from the manufacturer, not from a rule of thumb.

Solid steel discs on the reference set: 1,480 rpm, 685 kW, bare set 25 kg·m², start with 50 % of rated torque available for acceleration
Disc r × tMass (kg)Disc I (kg·m²)Disc GD² (kg·m²)Rim speed (m/s)Total I (kg·m²)τ (s)Start (s)
0.60 m × 100 mm88816063993.01853.2413.0
0.71 m × 120 mm (design)1,4923761,504110.04017.0328.1
0.80 m × 120 mm1,8946062,424124.063111.0644.3
1.00 m × 150 mm3,6991,8507,398155.01,87532.87131.5

In our practice, these items decide whether the flywheel gets built:

8 · The motor start penalty

Whatever the flywheel gives back on a trip, the motor must put in on every start:

\[ t_{acc} \;=\; \frac{I\,\omega_0}{f\,T_0} \;=\; \frac{2\tau}{f} \]

Here \(f\) is the average fraction of rated torque left over to accelerate the rotor: motor torque minus the pump's load torque. Taking \(f\) = 0.5 is our screening assumption, not a motor constant. At \(f\) = 0.5 the start takes 4τ: 1.8 s for the bare set, 7.0 s at 100 kg·m², 28.1 s at 400 and 56.1 s at 800.

9 · Interactive: flywheel and motor calculator

Size a disc for your own set. The chart plots τ and start time against total inertia, with your set marked on both lines.

Flywheel disc, stored energy and motor start time
Solid steel disc, \(\rho\) = 7,850 kg/m³: \(m = \rho\pi r^2 t\), \(I = \tfrac12 m r^2\). \(E_k = \tfrac12 I\omega_0^2\), \(\tau = E_k/P\), start time \(= I\omega_0/(f\,T_0)\) with \(T_0 = P/\omega_0\). The green line is the τ at which the reference main was protected (7.0 s at 1,480 rpm and 685 kW), shown as a screening marker for that main only.
Inertia goes with r⁴, rim speed with r.
Inertia and mass both in proportion.
From the data sheets: GD² ÷ 4, or WR² (lb·ft²) × 0.0421.
Sets both τ and the rated torque.
Stored energy goes with speed squared.
Average motor torque minus pump load torque, as % of rated. 50 % is our screening assumption.
Disc mass
1,894 kg
Disc I
606 kg·m²
Total I
631 kg·m²
Total GD²
2,524 kg·m²
Stored energy
7.58 MJ
τ = E/P
11.06 s
Start time
44.3 s
Rim speed
124.0 m/s

The default disc, 0.8 m × 120 mm, weighs 1,894 kg and adds 606 kg·m². That makes 631 kg·m² in total, 7.58 MJ stored, τ = 11.06 s, a 44.3 s start and a rim speed of 124.0 m/s, which is more than the main needs. Drag the radius to 0.71 m: the disc drops to 1,492 kg and 376 kg·m², the total is the design 401 kg·m², and the start takes 28.1 s. Raise the torque available to accelerate to 100 % and the start halves to 14.1 s. Switch to 2,960 rpm and the same disc stores four times the energy, but its rim speed doubles to 248.0 m/s, well beyond the discs of section 7; get the permissible rim speed from the manufacturer.

10 · When a flywheel is the wrong answer

Flywheels belong on short and medium-length mains, where a modest addition of inertia buys a rundown long enough to matter over the wave period [5][6][11]. Outside that range the answer is usually something else.

A surge relief valve covers only the pressure side. The full comparison is in Choosing surge protection on one pipeline.

11 · Setting it up in Bentley HAMMER

This is how we check a flywheel in HAMMER [10]; the general transient workflow and common pitfalls are covered separately. Field names differ slightly between HAMMER versions.

  1. Build the line and check the steady state. Model a Reservoir for the wet well, the Pump, Pipes with Junctions along the route, and a delivery Reservoir (see boundary conditions). Set wave speeds with the Wave Speed Calculator (see Wave speed), and confirm the steady HGL (85.0 m and 44.8 m here) first.
  2. Define the pump. Enter the duty, pump curve, efficiency and speed (1,480 rpm). Select the 4-quadrant characteristic curves from the specific speed closest to the real impeller, or ask the manufacturer for the complete characteristics; the rundown torque is where the answer is sensitive.
  3. Enter the inertia (pump and motor) as the total on the shaft: rotor, impeller, coupling and flywheel. Check whether each value is \(I\), GD² or WR², and whether the pump value includes the water in the impeller. Before data sheets exist, estimate the inertia from published correlations with power and speed [6], and record the estimate as an assumption.
  4. Set the pump trip (shut down) at the start of the run, with every pump in the station tripping together.
  5. Model the check valve on the pump, first as ideal and then with its real closure time/delay. A valve that closes slowly passes reverse flow before it seats and can then slam; check the reverse velocity at closure with the real closure time, even though the flywheel delays the reversal.
  6. Set up each inertia as its own alternative and scenario (the six inertias of the tables, 25 to 1,600 kg·m²), changing nothing else, so the envelopes can be compared and the threshold bracketed.
  7. Set the transient run options. Make the run duration several times 2L/a (150 s here). Accept the time step computed from the shortest pipe and the wave speeds, and check that the wave speed adjustment tolerance has not moved any wave speed materially. Turn on vapour pressure / column separation, and keep the friction method the same in every run.
  8. Read the time histories in the Transient Results Viewer: head and pump speed at the pump. A set that stops in about a second when you expected a slow rundown usually has its inertia in the wrong unit.
  9. Read the profile, not the pump. Plot the profile (path) with maximum and minimum head envelopes. Check the lowest head and its chainage against +3.0 m and the maximum against 136 m, and use the animation to see where the minimum forms.
  10. Test the sensitivity. Rerun the chosen inertia with neighbouring specific-speed characteristics, the lowest credible inertia and the wave speed range. If the minimum moves by more than your margin, the design is not yet robust.
  11. Close the loop on site. In our practice a recorded trip after commissioning is the best check of the modelled inertia; see SCADA records in transient analysis.

12 · Design checklist

Surge protection design series
  1. Wave speed: the number that sets the surge
  2. The differential orifice: empty freely, refill slowly
  3. Bladder, diaphragm or air-over-water vessel
  4. One-way surge tanks at the knee
  5. Surge relief valves: what a valve at the pump can protect
  6. Pump inertia and the flywheel
  7. Choosing surge protection on one pipeline
The sizing method itself is in Sizing the Hydropneumatic Surge Vessel.

References & standards

  1. ISO 2531 Ductile iron pipes, fittings, accessories and their joints for water applications — the DN800 K9 ductile iron pipe of the reference main.
  2. EN 805 Water supply — Requirements for systems and components outside buildings — design and allowable pressures, and the surge allowance the maximum is checked against.
  3. Wylie, E.B. & Streeter, V.L. Fluid Transients in Systems. Prentice Hall, 1993 — the Joukowsky relation, the method of characteristics, pump boundary conditions with rotor inertia, and the discrete vapour cavity model.
  4. Boulos, P.F., Karney, B.W., Wood, D.J. & Lingireddy, S. “Hydraulic transient guidelines for protecting water distribution systems.” Journal AWWA, 97(5), 2005 — overview of transient control strategies, including added pump inertia.
  5. Chaudhry, M.H. Applied Hydraulic Transients, 3rd ed. Springer, 2014 — pump rundown equations, complete (four-quadrant) pump characteristics, and flywheels as a control measure.
  6. Thorley, A.R.D. Fluid Transients in Pipeline Systems, 2nd ed. Professional Engineering Publishing, 2004 — pump trip and rundown, estimating pump and motor inertia, and where flywheels are and are not practical.
  7. Parmakian, J. Waterhammer Analysis. Dover, 1963 — classical charts for power failure at a pump, built on pipeline and pump inertia parameters.
  8. Bergant, A., Simpson, A.R. & Tijsseling, A.S. “Water hammer with column separation: a historical review.” Journal of Fluids and Structures, 22(2), 2006 — column separation, vapour and gas cavity models, and why collapse peaks are model-sensitive.
  9. Larock, B.E., Jeppson, R.W. & Watters, G.Z. Hydraulics of Pipeline Systems. CRC Press, 2000 — homologous pump relations and transient pump boundary conditions.
  10. Bentley Systems. OpenFlows HAMMER product documentation and help — pump inertia, pump trip, 4-quadrant characteristic curves, check valve closure, transient run options and results viewing.
  11. Stephenson, D. Pipeline Design for Water Engineers, 3rd ed. Elsevier, 1989 — water hammer protection of pumping lines, and the practical limits of flywheels on long mains.
#SurgeAnalysis #WaterHammer #HydraulicTransients #PumpTrip #PowerFailure #PumpInertia #Flywheel #GD2 #MomentOfInertia #PumpRundown #ColumnSeparation #VapourCavity #MethodOfCharacteristics #BentleyHAMMER #OpenFlowsHAMMER #SurgeProtection #TransmissionMain #DuctileIron #PumpStation #PumpStationDesign #MotorStarting #VariableSpeedDrive #SoftStarter #FourQuadrantCharacteristics #CheckValve #TransientAnalysis #PipelineDesign #HydraulicDesign #WaterInfrastructure #WaterEngineering