Hazen–Williams equation

The Hazen–Williams equation is an empirical relationship which relates the flow of water in a pipe with the physical properties of the pipe and the pressure drop caused by friction. It is used in the design of water pipe systems[1] such as fire sprinkler systems,[2] water supply networks, and irrigation systems. It is named after Allen Hazen and Gardner Stewart Williams.

The Hazen–Williams equation has the advantage that the coefficient C is not a function of the Reynolds number, but it has the disadvantage that it is only valid for water. Also, it does not account for the temperature or viscosity of the water.[3]

General form

Henri Pitot discovered that the velocity of a fluid was proportional to the square root of its head in the early 18th century. It takes energy to push a fluid through a pipe, and Antoine de Chézy discovered that the head loss was proportional to the velocity squared.[4] Consequently, the Chézy formula relates hydraulic slope S (head loss per unit length) to the fluid velocity V and hydraulic radius R:

V=C\sqrt{RS}=C\, R^{0.5}\, S^{0.5}

The variable C expresses the proportionality, but the value of C is not a constant. In 1838 and 1839, Gotthilf Hagen and Jean Léonard Marie Poiseuille independently determined a head loss equation for laminar flow, the Hagen–Poiseuille equation. Around 1845, Julius Weisbach and Henry Darcy developed the Darcy–Weisbach equation.[5]

The Darcy-Weisbach equation was difficult to use because the friction factor was difficult to estimate.[6] In 1906, Hazen and Williams provided an empirical formula that was easy to use. The general form of the equation relates the mean velocity of water in a pipe with the geometric properties of the pipe and slope of the energy line.

V = k\, C\, R^{0.63}\, S^{0.54}

where:

The equation is similar to the Chézy formula but the exponents have been adjusted to better fit data from typical engineering situations. A result of adjusting the exponents is that the value of C appears more like a constant over a wide range of the other parameters.[7]

The conversion factor k was chosen so that the values for C were the same as in the Chézy formula for the typical hydraulic slope of S=0.001.[8] The value of k is 0.001−0.04.[9]

Typical C factors used in design, which take into account some increase in roughness as pipe ages are as follows:[10]

Material C Factor low C Factor high Reference
Asbestos-cement 140 140 -
Cast iron new 130 130 [10]
Cast iron 10 years 107 113 [10]
Cast iron 20 years 89 100 [10]
Cast iron 30 years 75 90 [10]
Cast iron 40 years 64 83 [10]
Cement-Mortar Lined Ductile Iron Pipe 140 140 -
Concrete 100 140 [10]
Copper 130 140 [10]
Steel 90 110 -
Galvanized iron 120 120 [10]
Polyethylene 140 140 [10]
Polyvinyl chloride (PVC) 150 150 [10]
Fibre-reinforced plastic (FRP) 150 150 [10]

Pipe equation

The general form can be specialized for full pipe flows. Taking the general form

V = k\, C\, R^{0.63}\, S^{0.54}

and exponentiating each side by 1/0.54 gives (rounding exponents to 2 decimals)

V^{1.85} = k^{1.85}\, C^{1.85}\, R^{1.17}\, S

Rearranging gives

S = {V^{1.85} \over k^{1.85}\, C^{1.85}\, R^{1.17}}

The flow rate Q = V A, so

S = {V^{1.85} A^{1.85}\over k^{1.85}\, C^{1.85}\, R^{1.17}\, A^{1.85}} = {Q^{1.85}\over k^{1.85}\, C^{1.85}\, R^{1.17}\, A^{1.85}}

The hydraulic radius R (which is different from the geometric radius r) for a full pipe of geometric diameter d is d/4; the pipe's cross sectional area A is \pi d^2 / 4, so

S =  {4^{1.17}\, 4^{1.85}\,Q^{1.85}\over \pi^{1.85}\,k^{1.85}\, C^{1.85}\, d^{1.17}\, d^{3.70}}
=  {4^{3.02}\,Q^{1.85}\over \pi^{1.85}\,k^{1.85}\, C^{1.85}\, d^{4.87}}
=  { 4^{3.02} \over \pi^{1.85}\,k^{1.85}} {Q^{1.85}\over C^{1.85}\, d^{4.87}}
=  { 7.916 \over k^{1.85}} {Q^{1.85}\over C^{1.85}\, d^{4.87}}

U.S. customary units (Imperial)

When used to calculate the pressure drop using the US customary units system, the equation is:[11]

S_{\mathrm{psi\ per\ foot}} = \frac{P_d}{L} = \frac{4.52\ Q^{1.85}}{C^{1.85}\ d^{4.87}}

where:

Note: Caution with U S Customary Units is advised. The equation for head loss in pipes, also referred to as slope, S, expressed in "feet per foot of length" vs. in 'psi per foot of length' as described above, with the inside pipe diameter, d, being entered in feet vs. inches, and the flow rate, Q, being entered in cubic feet per second, cfs, vs. gallons per minute, gpm, appears very similar. However, the constant is 4.75 vs. the 4.52 constant as shown above in the formula as arranged by NFPA for sprinkler system design. The exponents and the Hazen-Williams "C" values are unchanged.

SI units

When used to calculate the head loss with the International System of Units, the equation becomes:[12]

S = \frac{h_f}{L} = \frac{10.67\ Q^{1.85}}{C^{1.85}\ d^{4.87}}

where:

Note: pressure drop can be computed from head loss as hf × the unit weight of water (e.g., 9810 N/m3 at 4 deg C)

See also

References

  1. "Hazen–Williams Formula". Retrieved 2008-12-06.
  2. "Hazen–Williams equation in fire protection systems". Canute LLP. 27 January 2009. Archived from the original on 2013-04-06. Retrieved 2009-01-27.
  3. Brater, Ernest F.; King, Horace W.; Lindell, James E.; Wei, C. Y. (1996). "6". Handbook of Hydraulics (Seventh ed.). New York: McGraw Hill. p. 6.29. ISBN 0-07-007247-7.
  4. Walski, Thomas M. (March 2006), "A history of water distribution", Journal of the American Water Works Association (American Water Works Association) 98 (3): 110–121, p. 112.
  5. Walski 2006, p. 112
  6. Walski 2006, p. 113
  7. Williams & Hazen 1914, p. 1, stating "Exponents can be selected, however, representing approximate average conditions, so that the value of c for a given condition of surface will vary so little as to be practically constant."
  8. Williams & Hazen 1914, p. 1
  9. Williams & Hazen 1914, pp. 1–2
  10. 1 2 3 4 5 6 7 8 9 10 11 12 Hazen-Williams Coefficients, Engineering ToolBox, retrieved 7 October 2012
  11. 2007 version of NFPA 13: Standard for the Installation of Sprinkler Systems, page 13-213, eqn 22.4.2.1
  12. "Comparison of Pipe Flow Equations and Head Losses in Fittings" (PDF). Retrieved 2008-12-06.

External links

This article is issued from Wikipedia - version of the Monday, March 21, 2016. The text is available under the Creative Commons Attribution/Share Alike but additional terms may apply for the media files.