Pumping System Head Estimation
Pumping System Head Estimation
Pumping System
Head Estimation
These spreadsheets highlight losses,
uncertainties, and value of power used,
James Anthony helping you optimize pumping system
Lockwood Greene
designs and save on operating costs.
T
O DESIGN A PUMPING SYSTEM WELL, YOU suction friction head hfs consumed on the way to the pump.
need to accurately estimate the total head, available
suction total head, driver power, and energy cost. hts = hus – hfs (3)
You can estimate these quantities and their uncertainties
for pumping systems with as many as 12 piping runs using a The discharge total head is the discharge friction head hfd
spreadsheet I’ve written for centrifugal pumps or a compan- consumed on the way to the destination, plus the destination
ion spreadsheet for positive displacement pumps, both of usable head hud.
which are available at www.cepmagazine.org (1). The spread-
htd = hfd + hud (4)
sheets handle flows in all Reynolds number regimes — turbu-
lent, intermittently turbulent, and laminar — for fluids for The usable head is sum of the elevation head hue, the pres-
which the viscosity and density can be taken as constant sure head hup, and the velocity head huv.
along a given piping run.
In this article I will describe the formulas in the hu = hue + hup + huv (5)
spreadsheets, give tips on choosing input values, and
review a sample calculation. The usable head is also the destination usable head hud
minus the source usable head hus, and its components are
Formulas defined similarly. Here, the subscript “d” stands for “destina-
First I’ll cover the formulas used to calculate total head, tion” and the subscript “s” stands for “source.”
usable head, friction head, uncertainties, the piping complexi- hu = hud – hus (6)
ty factor, available suction total head, power required, and
hue = hued – hues (7)
energy cost.
The total head ht that a pump provides is the usable head hup = hupd – hups (8)
hu plus the friction head hf. huv = huvd – huvs (9)
ht = hu + hf (1) The source usable head is the sum of the source elevation
head, the source pressure head, and the source velocity head.
The total head is also the discharge total head htd minus
the suction total head hts. hus = hues + hups + huvs (10)
ht = htd – hts (2)
The source elevation head equals the source elevation zs,
The suction total head is the source usable head hus minus the which is the elevation of the source fluid above the eye of
For an exit:
mula for intermittently turbulent flows and laminar flows (8).
2
For Re > 5,000: ⎛ ⎞
⎜ q ⎟
⎜ 2⎟
⎛ ⎛ ∆ r − 5.042 × ⎞⎞
–2
⎜ π⎛ d⎞ ⎟
⎜ ⎝ ⎠ ⎟
⎜ ⎜ 3.7065 ⎟⎟ ⎝ 2 ⎠
Re h ffx = (31)
f = 0.25⎜ –2 log⎜ ⎟⎟ (25) 2g
⎜ log⎛ ∆ r 5.8506 ⎞ ⎟ ⎟
1.1098
⎜ +
⎜ ⎜ ⎜⎝ 2.8257 Re 0.8981 ⎟⎠ ⎟ ⎟
⎝ ⎝ ⎠⎠ The increaser head or reducer head of a given reducer is
estimated from the inlet Fanning friction factor f1, inlet and
outlet diameters d1 and d2, inlet velocity v1, and inlet
For Re ≤ 5,000:
Reynolds number Re1 (10).
For Re1 > 4,000:
⎡⎛ 8 ⎞
12
f =2× +
⎢⎣⎝ Re ⎠ 2
⎛ ⎛ d ⎞ 2⎞ v 2
h ffi = (1 + 3.2 f1 )⎜1 − ⎜ 1 ⎟ ⎟ 1
1
3
⎤ (32)
⎝ ⎝ d2 ⎠ ⎠ 2 g
– 12
⎛⎛ ⎞
16
⎛ ⎞⎞ 2
⎥
⎜⎜ ⎜ 1 ⎟⎟ ⎛ 37, 530 ⎞ ⎟
16
⎥
⎜ ⎜ 2.457 ln ⎜ 7 +
⎝ Re ⎠ ⎟
0.9 ⎟⎟ ⎥ (26)
⎜⎜ ⎜⎛ ⎞ + 0.27 ∆ r ⎟⎟ ⎟ ⎥ For Re1 ≤ 4,000:
⎝⎝ ⎝ ⎝ Re ⎠ ⎠⎠ ⎠ ⎥⎦
⎛ ⎛ d ⎞ 4⎞ v 2
h ffi = 2⎜1 − ⎜ 1 ⎟ ⎟ 1 (33)
The relative equivalent sand roughness is calculated from the ⎝ ⎝ d2 ⎠ ⎠ 2 g
absolute equivalent sand roughness ε and the diameter.
For all Re1:
∆r = ε/d (27)
60 ⎞ ⎛ ⎛ d1 ⎞ ⎞v2
4
The Reynolds number is calculated from the density, diame- ⎛
h ffr = ⎜ 0.1 + ⎟ ⎜ ⎜ ⎟ − 1⎟ 1 (34)
ter, velocity, and viscosity µ. ⎝ Re 1 ⎠ ⎝ ⎝ d2 ⎠ ⎠ 2g
Re = ρdv/µ (28) The heads for other fittings and for valves other than throt-
tling control valves are estimated from Re, d, and v using the
The fitting and valve heads are the entrance and exit heads constants km, ki, and kd (11), the constants k1 and kinf (9), or
hffn and hffx, increaser and reducer heads hffi and hffr, and heads constant k values (12–14).
for other fittings and for valves other than throttling control
valves hffo. ⎛ ⎛ ⎞⎞
The entrance head is estimated from the Reynolds num- ⎜k ⎜ kd ⎟ ⎟ v 2
ber and the velocity (9). The exit head — defined logically h ffo = ⎜ m + ki ⎜1 + 0.3 ⎟ ⎟ (35)
⎜ Re ⎜ ⎛ d ⎞ ⎟ ⎟ 2g
here as the velocity head lost to friction when the exiting ⎜ ⎜ ⎝ in. ⎠ ⎟⎠ ⎟⎠
⎝ ⎝
flow discharges above the liquid surface or enters a vessel
below the liquid surface — is calculated from the flow and
the exit pipe diameter. ⎛ ⎛ ⎞⎞
For a flush entrance: ⎜ k1 ⎜ 1 ⎟ ⎟ v2
h ffo =⎜ + kinf ⎜1 + ⎟⎟ (36)
⎜⎜ Re ⎜ ⎛ d ⎞ ⎟ ⎟ 2g
v2 ⎝ ⎝ ⎝ in. ⎠ ⎠ ⎟⎠
h ffn = ⎛ + 0.50⎞
160
(29)
⎝ Re ⎠ 2g
The brake horsepower Pb is the pumping power divided might be just 70% of the values used for pump rating. So
by the pump efficiency ηp. the total head that the system will actually need in order
to operate might be around 90% of the rating estimate,
Pb = Pp/ηp (42) and might stay that way for years — maybe for the life-
time of the system. Besides that, the pump will very like-
The energy cost present value e is calculated from esti- ly deliver a little more head than the rated head. In a
mates of the fraction of time the pumping system is run- throttled system, the extra head will need to be taken up
ning at design conditions t, brake horsepower, driver or by the throttling valve. In an unthrottled system, the flow
motor efficiency ηd, utility rate r, real interest rate i, and will need to increase also. Such excesses in head and
design life n (18). flow waste energy.
You can provide for uncertainty, wear, and future
requirements easily enough without using a flow safety
Pb (1 + i ) n − 1
e=t r ( 43) factor. The uncertainty is often provided for, to some
ηd i(1 + i ) n extent, by allowances for corrosion in heat exchangers
and in piping, which often provide more conservative
When comparing different design options, it’s best to com- safety margins than turn out to be needed for a given sys-
pare their total costs, including both capital and operating tem. Wear is greatest for slow, low-flow, high-head
costs. Since the capital costs are present values, it makes pumps of soft construction that handle abrasive flows;
sense to convert the operating costs to present values. The for systems susceptible to wear, wear can be managed by
operating cost that is largest and that can be estimated the choosing pumps of suitable construction or by restoring
most accurately is the energy cost. You can convert the internal clearances, if needed over time (19). Future
energy cost to a present value using the equation above by capacity increases can be provided for by installing larg-
entering the long-term stock market real return as the real er impellers. When you make reasonable plans to provide
interest rate, entering the present utility rate, and entering for future needs, and you match present designs to pres-
the presently-allowed depreciation life, or by entering bet- ent needs, you’ll get guaranteed savings up front.
ter estimates if available. The present value of the operat- You’ll need to run the highest temperature case rou-
ing cost for different design options tells you how much tinely. Temperature-induced changes in the density and
capital cost you are justified in spending for different viscosity do not change the total head much, but temper-
design options. ature-induced changes in the vapor pressure can affect
the available suction total head substantially.
Choosing input values Increased densities at lower temperatures decrease the
You should choose the flow for the rating case so that pressure head and friction head, but increase the power.
it’s the maximum flow needed in the near term, with no Although the total head usually doesn’t change much, the
safety factor. power does — it increases almost directly with the density.
If the rated flow exceeds the near-term need, the oper- Increased viscosities at lower temperatures increase the
ating cost rises quickly. Increasing the flow by 10% friction head, but so little that for low-viscosity flows,
increases the friction head by 21% in turbulent flow, and order-of-magnitude estimates of viscosity are normally all
by 10% even in laminar flow. that are needed. Consider the system in the Figure at start-
If the rated flow is too high, the reliability falls also, up, taking the pump’s head as approximately constant. For
because wear increases and pump internal recirculation this low-viscosity system, reducing the temperature by
increases. Wear increases because the pump produces 190°C increases the viscosity 250 times but reduces the
more head and, in unthrottled systems, produces more flow only by 1/2. For high-viscosity systems with viscosi-
flow. Recirculation increases because the flows at the ties exceeding 100 cP, though, reducing the temperature
lower-flow operating points become lower fractions of increases the viscosity faster, and as a result, more accurate
the rated flow. Recirculation is more of a problem for estimates of viscosity are needed.
pumps that have higher heads per stage, because it begins Increased vapor pressures at higher temperatures
closer to the rated flow and builds more quickly to dam- decrease the available suction total head significantly.
aging levels. For both water and heptane, the increase in the vapor
When a pumping system is operated, the source level pressure head with temperature is on the order of 2% per
might be 5 ft higher than the value used for pump rating, degree F. Because of this, you need to be careful to iden-
and the like-new heat exchanger and straight pipe heads tify the maximum temperature.
Heater
Ele
va
tio i
n= ps
48
ft 13
Suction Drum ft
13.5 .5
= 15
psi
10 tion
ft va
Ele
30
n. ft
4i
10
ft
n.
4i Ele
ft ia. va
24 .d tio
f t -in 4-i n=
10 18 n. 12
ft Orifice Sc
3 h. ft
40
Discharge Vessel
Ele ft
ter
va
tio = 12
He
a n ion 40
=2
va
t h.
ft Sc
Ele -in
.
t 6
5f
i
ps ft
5.2 15
P= 4-in. Control Valve
Ele
va
tio
n=
1.5 on
ft cti ft
Su 1.5
n.
4-i ion
=
t
va
Ele
3-in. Discharge
Sample calculation one 4-in. × 3-in. increaser, five gate valves, 20 standard
In the Figure, a centrifugal pump transfers gas oil radius elbows, and a 4.026-in.-dia. exit pipe where the
from a suction drum to another vessel. Table 2 shows the velocity head is lost. The discharge line is a little simpler
inputs and the results. than normal piping.
I took the rated flow to be 250 gal/min, with no safety fac- The rated total head is 83 ±3 ft. The usable head is
tor. The fluid temperature is 555°F; ρ is 64.87 lb/ft3; µ is 0.6 40% of the total head and the friction head is 60% of the
cP; zs is 14 ft; ps is 13 psig. I took the source diameter ds to be total head. The available suction total head is 76 ft. At
100 in., which is much larger than the pump suction line. All the spreadsheet-estimated pump efficiency of 66%, the
the piping is flanged schedule 40 carbon steel. I took its brake horsepower is 8.3 hp. If the pump is onstream all
absolute equivalent sand roughness ε to be 0.2 mm, which is but 11 days a year (t = 97%), and if the motor efficiency
the value listed in the spreadsheet for corroded steel pipe car- is 90%, the utility rate is 4.45 ¢/kWh, the design life is
rying oil. The suction line diameter is 6 in. and its length is 39 9.5 yr, and the real interest rate is 6%/yr, then per Eq. 43,
ft. The line has one flush entrance, five standard radius the present value of the lifetime energy cost is $19,000.
elbows, one gate valve, one strainer, and one 6-in. × 4-in.
reducer. I took the strainer’s number of velocity heads k from Conclusion
Ref. 14 to be 0.88. The suction line is of normal complexity. The methods described here will help you be sure of the
On the discharge side, zd is 46.5 ft and pd is 13.5 psig. magnitude of the head losses due to corrosion and be sure of
I took the destination diameter dd to be 100 in., which the magnitude of the uncertainties in your head estimates, so
again is much larger than the pipeline. The orifice plate, you won’t have to use safety factors you don’t need.
control valve, and heat exchanger pressure drops are Since you can quickly estimate operating energy
1.52 psi, 5.2 psi, and 10 psi, respectively. The discharge costs, you’ll find it easier to evaluate variable speed
line diameter is 4 in. and its length is 177 ft. The line has drives that can save up to 80% of energy costs, or
† Underlines prompt users to enter values, such as fitting quantities, or to check values suggested by the spreadsheet. Calculated values of pumping system heads are shown at the
right in the “Suction” and “Discharge” columns and are summarized in the “Overall” column. The spreadsheets include additional header lines, valve types, and optional and explanatory
worksheets that are not shown in this table.
Liquids Handling
Literature Cited
1. Anthony, J. M., “Pumping System Head Estimate,” v. 1.0, the Ninth International Pump Users Symposium,” Turbomachinery
02CEP_P.xls and 02CEP_PD.xls, available via Laboratory, Dept. of Mechanical Engineering, Texas A&M Univ.,
www.cepmagazine.org (2005). College Station, pp. 51–60 (1992); available at
2. Perry, R. H., et al., Perry’s Chemical Engineers’ Handbook, 7th ed., www.gouldspumps.com/ download_files/Technews/desconsi.pdf.
McGraw-Hill, New York, p. 10-22 (1997). 24. Raju, K. S. N., and J. Chand, “Consider the Plate Heat
3. Perry, R. H., et al., p. 6-7. Exchanger,” Chem. Eng., 87 (16), pp. 133–144 (Aug. 11, 1980).
4. Streeter, V. L., “The Kinetic Energy and Momentum Correction 25. Jenssen, S. K., “Assessment of Heat Exchanger Data,” in Cichelli,
Factors for Pipes and Open Channels of Great Width,” Civil M., “Heat Transfer – Storrs,” Chem. Eng. Prog. Symposium Series,
Engineering, 12 (4), pp. 212–213 (1942). 56 (30), AIChE, New York, pp. 195–201 (1960).
5. Perry, R. H., et al., p. 6-10. 26. Baumann, H. D., “Control Valve Primer: A User’s Guide,”
6. Perry, R. H., et al., p. 6-17. Instrument Soc. of America, Research Triangle Park, NC (1991).
7. Chen, N. H., “An Explicit Equation for Friction Factor in Pipe,” 27. Lefebvre, A. H., “Atomization and Sprays,” Hemisphere, New York, pp.
Ind. Eng. Chem. Fundam., 18 (3), pp. 296–297 (1979). 157, 204 (Merrington and Richardson correlation) (1989).
8. Churchill, S. W., “Friction-Factor Equation Spans All Fluid-Flow 28. Budris, A. R., “Designing Pump Piping,” Chem. Processing, 65 (8),
Regimes,” Chem. Eng., 84 (24), pp. 91–92 (Nov. 7, 1977). pp. 33–37 (Aug. 2002); available at www.chemicalprocessing.com/
9. Hooper, W. B., “The Two-K Method Predicts Head Losses in Pipe Web_First/cp.nsf/ArticleID/CBOH-5CZJW3.
Fittings,” Chem. Eng., 88 (17), pp. 96–100 (Aug. 24, 1981). 29. Branan, C., “Rules of Thumb for Chemical Engineers,” Gulf
10. Hooper, W. B., “Calculate Head Loss Caused by Change in Pipe Publishing, Houston, TX, p. 80 (1994).
Size,” Chem. Eng., 95 (16), pp. 89–92 (Nov. 7, 1988). 30. National Electrical Manufacturers Assn., “Energy Management
11. Darby, R., “Correlate Pressure Drops Through Fittings,” Chem. Guide for Selection and Use of Fixed Frequency Medium AC
Eng., 108 (4), pp. 127–130 (April 2001); 108 (12), pp. 8–10 Squirrel-Cage Polyphase Induction Motors,” NEMA Standards
(Nov. 2001). Publication MG 10-2001, National Electrical Manufacturers
12. Perry, R. H., et al., p. 6-18. Association, Rosslyn, Virginia, p. 9 (2001); available at
13. Crane Valve, “Flow of Fluids Through Valves, Fittings, and Pipe,” www.nema.org/DocUploads/57D634E6-FF3C-4492-
Technical Paper No. 410, pp. A-27 to A-29 (1988). B3D73108C1E4F02D/MG10.pdf.
14. Hydraulic Institute, “Hydraulic Institute Engineering Data 31. U.S. Census Bureau, “Statistical Abstract of the United States:
Book,” 2nd ed., Hydraulic Institute, Cleveland, OH, pp. 25, 47, 2003,” 123 ed., U.S. Census Bureau, Washington, DC, p. 572 (2002);
and 77–79 (1990). available at www.census.gov/prod/2003pubs/02statab/energy.pdf.
15. Brown, G. S., “How to Predict Pressure Drop before Designing the 32. Siegel, J., “The Rise in Stock Valuations and Future Equity Returns,”
Piping,” Chem. Eng., 94 (4), pp. 85–86 (Mar. 16, 1987). The Journal of Management Consulting, 5 (1), pp. 9–19, Denver, CO
16. Hooper, W. B., “Predict Fittings for Piping Systems,” Chem. Eng., (June/July 2002); available at www.jeremysiegel.com/
89 (10), pp. 127–129 (May 17, 1982). view_article.asp?p=233.
17. Perry, R. H., et al., p. 10-23. 33. Internal Revenue Service, “How To Depreciate Property,”
18. Perry, R. H., et al., p. 9-13. Publication 946, Cal. No. 13081F, pp. 93–102 (2002); available at
19. Karassik, I. J., “Centrifugal Pumps and System Hydraulics,” Chem. www.irs.gov/pub/irs-pdf/p946.pdf.
Eng., 89 (10), pp. 84–106 (Oct. 4, 1982). 34. Kern, R., “How To Size Piping for Pump-Discharge Conditions,”
20. Chen, C. C., “Cope with Dissolved Gases in Pump Calculations,” Chem. Eng., 82 (11), pp. 113–120 (May 26, 1975).
Chem. Eng., 100 (10), pp. 106–112 (Oct. 1993). 35. Xenergy, Inc., “United States Industrial Electric Motor Systems
21. Goulds Pumps, “Goulds Pump Manual, 6th Ed.,” Goulds Pumps, Market Opportunities Assessment,” U.S. Dept. of Energy (DOE),
Seneca Falls, New York, pp. 607–611 (1995); available via Washington, DC, pp. 13–17, 56–59 (1998); available at
http://www.gouldspumps.com/cattechnews.ihtml?pid=76&last- www.oit.doe.gov/bestpractices/ pdfs/mtrmkt.pdf.
catid=46&step=4. 36. Lawrence Berkeley National Laboratory and Resource
22. Doolin, J. H., “Centrifugal Pumps and Entrained-Air Problems,” Dynamics Corp., “Improving Pumping System Performance:
Chem. Eng., 70 (1), pp. 103–106 (Jan. 7, 1963). A Sourcebook for Industry,” DOE, Office of Industrial
23. Cappellino, C. A., et al., “Design Considerations and Application Technologies, Motor Challenge Program, Washington, DC (1999);
Guidelines for Pumping Liquids with Entrained Gas using Open available at www.oit.doe.gov/bestpractices/pdfs/pump.pdf.
Impeller Centrifugal Pumps,” in Bailey, J. C., Ed., “Proceedings of