Darcy-Weisbach equation

From Wikipedia, the free encyclopedia

The Darcy-Weisbach equation is an important and widely used phenomenological equation in hydraulics. It relates the head loss or pressure loss due to friction along a given length of pipe to the average velocity of the fluid flow.

Contents

[edit] Head loss form

Head loss can be calculated with

h_f = f \cdot \frac{L}{D} \cdot \frac{V^2}{2g}

where

  • hf is the head loss due to friction;
  • L is the length of the pipe;
  • D is the hydraulic diameter of the pipe (for a pipe of circular section, this equals the internal diameter of the pipe);
  • V is the average velocity of the fluid flow, equal to the volumetric flow rate per unit cross-sectional wetted area;
  • g is the local acceleration due to gravity;
  • f is a dimensionless coefficient called the Darcy friction factor.

[edit] Pressure loss form

Given that the head loss hf expresses the pressure loss Δp as the height of a column of fluid,

Δp = ρghf

where ρ is the density of the fluid, the Darcy-Weisbach equation can also be written in terms of pressure loss:

\Delta p = \lambda \cdot \frac{L}{D} \cdot \frac{\rho V^2}{2}

where the pressure loss due to friction Δp is a function of:

  • the ratio of the length to diameter of the pipe, L/D;
  • the density of the fluid, ρ;
  • the average velocity of the flow, V, as defined above;
  • a (dimensionless) coefficient of laminar, or turbulent flow, λ

Since the pressure loss equation can be derived from the head loss equation by multiplying each side by ρ and g, clearly λ equals f. The use of different symbols for the same numerical coefficient depending on whether head loss or pressure is considered is a historical accident due to different conventions being used by different communities of scientists and engineers.

[edit] Darcy friction factor

See also Darcy friction factor formulae

The friction factor f or flow coefficient λ is not a constant and depends on the parameters of the pipe and the velocity of the fluid flow, but it is known to high accuracy within certain flow regimes. It may be evaluated for given conditions by the use of various empirical or theoretical relations, or it may be obtained from published charts. These charts are often referred to as Moody diagrams, after L. F. Moody, and hence the factor itself is sometimes called the Moody friction factor. It is also sometimes called the Blasius friction factor, after the approximate formula he proposed.

For laminar (slow) flows, it is a consequence of Poiseuille's Law that λ equals 64/Re, where Re is the Reynolds Number calculated substituting for the characteristic length the hydraulic diameter of the pipe, which equals the inside diameter for circular pipe geometries.

For turbulent flow, methods for finding the friction factor f include using a diagram such as the Moody chart; or solving equations such as the Colebrook-White equation, or the Swamee-Jain equation. While the diagram and Colebrook-White equation are iterative methods, the Swamee-Jain equation allows f to be found directly for full flow in a circular pipe.

[edit] Confusion with the Fanning friction factor

The Darcy-Weisbach friction factor is 4 times larger than the Fanning friction factor, so attention must be paid to note which one of these is meant in any "friction factor" chart or equation being used. Of the two, the Darcy-Weisbach factor is more commonly used by civil engineers, and the Fanning factor by chemical and mechanical engineers, but care should be taken to identify the correct factor regardless of the source of the chart or formula.

Most charts or tables indicate the type of friction factor, or at least provide the formula for the friction factor with laminar flow. If the formula is f=16/Re, it's the Fanning factor, and if the formula is f=64/Re, it's the Darcy-Weisbach factor.

Which friction factor is plotted in a Moody diagram may be determined by inspection if the publisher did not include the formula described above:

  1. Observe the value of the friction factor for laminar flow at a Reynolds number of 1000.
  2. If the value of the friction factor is 0.064, then the Darcy friction factor is plotted in the Moody diagram. Note that the nonzero digits in 0.064 are the numerator in the formula for the laminar Darcy friction factor: f = 64/Re.
  3. If the value of the friction factor is 0.016, then the Fanning friction factor is plotted in the Moody diagram. Note that the nonzero digits in 0.016 are the numerator in the formula for the laminar Fanning friction factor: f = 16/Re.

The procedure above is similar for any available Reynolds number that is an integral power of ten. It is not necessary to remember the value 1000 for this procedure – only that an integral power of ten is of interest for this purpose.

[edit] History

Historically this equation arose as a variant on the Prony equation; this variant was developed by Henry Darcy of France, and further refined into the form used today by Julius Weisbach of Saxony in 1845. Initially, data on the variation of f with velocity was lacking, so the Darcy-Weisbach equation was outperformed at first by the empirical Prony equation in many cases. In later years it was eschewed in many special-case situations in favor of a variety of empirical equations valid only for certain flow regimes, notably the Hazen-Williams equation or the Manning equation, most of which were significantly easier to use in calculations. However, since the advent of the calculator, ease of calculation is no longer a major issue, and so the Darcy-Weisbach equation's generality has made it the preferred one.

[edit] Derivation

The Darcy-Weisbach equation is a phenomenological formula obtainable by dimensional analysis.

Away from the ends of the pipe, the characteristics of the flow are independent of the position along the pipe. The key quantities are then the pressure drop along the pipe per unit length, Δp/L, and the volumetric flow rate. The flow rate can be converted to an average velocity V by dividing by the wetted area of the flow (which equals the cross-sectional area of the pipe if the pipe is full of fluid).

Pressure has dimensions of energy per unit volume and, so, the pressure drop between two points must be proportional to (1/2)ρV2, which has the same dimensions as it resembles (see below) the expression for the kinetic energy per unit volume. We also know that it must be proportional to the length of the pipe between the two points L as the pressure drop per unit length is a constant. To turn that into a dimensionless quantity we can divide by the hydraulic diameter of the pipe, D, which is also constant along the pipe. Therefore,

\Delta p \propto \frac{L}{D} \times \frac{1}{2}\rho V^2

The proportionality coefficient is the dimensionless "Darcy friction factor" or "flow coefficient". This dimensionless coefficient will be a combination of geometric factors such as π, the Reynolds number and (outside the laminar regime) the relative roughness of the pipe (the ratio of the roughness height to the hydraulic diameter).

Note that (1/2)ρV2 is not the kinetic energy of the fluid per unit volume, for the following reasons. Even in the case of laminar flow, where all the flow lines are parallel to the length of the pipe, the velocity of the fluid on the inner surface of the pipe is zero due to viscosity, and the velocity in the center of the pipe must therefore be larger than the average velocity obtained by dividing the volumetric flow rate by the wet area. The average kinetic energy then involves the mean-square velocity, which always exceeds the square of the mean velocity. In the case of turbulent flow, the fluid acquires random velocity components in all directions, including perpendicular to the length of the pipe, and thus turbulence contributes to the kinetic energy per unit volume but not to the average lengthwise velocity of the fluid.

[edit] References

  • Moody diagram (GIF; ~25k)
  • The History of the Darcy-Weisbach Equation
  • De Nevers (1970) Fluid Mechanics. Addison-Wesley ISBN 0-201-01497-1
  • Shah R.K. and A.L. London, 1978, 'Laminar Flow Forced Convection in Ducts', Supplement 1 in Advances in Heat Transfer, Academic, NY, 1978.
  • Rohsenhow, W.M., J.P. Hartnett and E.N. Ganic, 1985, 'Handbook of Heat Transfer Fundamentals, 2nd Ed., McGraw-Hill Book Company.