Friction Factor Diagrams for Pipe Flow
Friction Factor Diagrams for Pipe Flow
ARROW@TU Dublin
2011-10-03
Part of the Energy Systems Commons, Hydraulic Engineering Commons, and the Nanoscience and
Nanotechnology Commons
Recommended Citation
McGovern, J.: Technical Note: Friction Diagrams for Pipe Flow. Technological University Dublin, 2011.
This Working Paper is brought to you for free and open access by the School of Mechanical and Design
Engineering at ARROW@TU Dublin. It has been accepted for inclusion in Articles by an authorized administrator of
ARROW@TU Dublin. For more information, please contact [email protected], [email protected],
[email protected], [email protected].
Abstract
This technical note describes diagrams of friction factor for pipe flow
that have been prepared using, mainly, the equations that Lewis Moody
used to prepare his famous diagram in 1944. The preparation of the
new diagrams was prompted by the need for vector graphics versions
that could be used for teaching purposes and that could be distributed
freely to students and others under a Creative Commons Attribution-
Share-Alike license.
Using a structure very similar to that of Moody’s diagram, variants
with the Darcy friction factor and the Fanning friction factor have
been prepared. In addition, variants have been prepared that include
not only monotonic roughness curves, but also inflectional roughness
curves.
Keywords: friction factor, Moody diagram, vector graphic, Darcy
friction factor, Fanning friction factor, pipe flow, relative roughness,
Reynolds number, monotonic roughness, inflectional roughness.
Introduction
Since the mid 1940s, practicing engineers, engineering academics and stu-
dents of engineering have made use of a diagram of friction factor for pipe
flow of the type that was published by Lewis Moody in the Transactions
of A.S.M.E. in 1944 [1]. Textbooks in such areas as Fluid Mechanics, Hy-
draulics, Heat Transfer or Unit Operations commonly include re-drawn or
reproduced versions of Moody’s diagram. The diagram is semi-empirical,
based on some fundamental principles and the strong intuition of leading
researchers up to 1944. At every stage since then the diagram could be
regarded as a temporary solution, until sufficient further advance had been
made. Although considerable progress has been made, notably in computa-
tion, in the measurement of surface roughness and in the measurement and
1
Darcy Friction Factor for Pipe Flow Reynolds Number
103 1.4 2 3 4 5 6 104 1.4 2 3 4 5 6 105 1.4 2 3 4 5 6 106 1.4 2 3 4 5 6 107 1.4 2 3 4 5 6 108 1.4 2 3 4 5 6 109 ¥
0.08 0.08
0.06
0.055
0.05
0.07 0.07
0.045
0.04
0.06 0.035 0.06
0.03
0.025
Transition Zone
N
0.05 0.05
rV 2
0.02
2
0.018
0.045 0.016 0.045
= 4 tw J
0.014
0.04 0.012 0.04
0.01
0.035 0.008 0.035
Complete Turbulence Zone
N
0.006
rV 2
2
Friction Factor
0.024 0.024
f f 0.002
0.022 ¥ 0.001 5 0.022
=
1.
01
N
0.018
an 0.000 8
0.018
d
Friction Factor fD = hf J DL
0.000 6
Re
0.016 0.000 4 0.016
=
0.015 16 0.000 3 0.015
00
0.014 Critical He 0.014
D
0.000 2
Zone
0.013 L 0.013
0.012 Indicative Roughness Values 0.000 1 0.012
e @mmD
0.011
Smooth honed steel 0.00065
Sm 0.011
Drawn tubing: glass, brass, copper, lead, plastic 0.0015
o oth 0.000 05
0.01 0.01
Asphalted cast iron 0.12 Pip 0.000 03
0.009 Galvanized steel 0.15 e 0.000 02 0.009
Laminar Wood stave 0.18-0.91
Cast iron 0.26
0.008 Zone Concrete 0.3-3
0.000 01 0.008
fD = 64Re Heavy brush coat: asphalts, enamels, tars 0.45-0.6
0.000 005
General tuberculation 1-3 mm 0.6-1.9
0.007 0.007
Riveted steel 0.9-9 0.000 003
Severe tuberculation and incrustation 2.5-6.5 0.000 002
0.006 0.006
0.000 001
5. ´ 10-7
0.005 0.005
103 1.4 2 3 4 5 6 104 1.4 2 3 4 5 6 105 1.4 2 3 4 5 6 106 1.4 2 3 4 5 6 107 1.4 2 3 4 5 6 108 1.4 2 3 4 5 6 109 ¥
rVD VD
Jim McGovern, Dublin Institute of Technology Reynolds Number Re = m
= n
http://www.fun-engineering.net/pipefrictionfactor.html Version DM001 2011.10.02
2
fF = fD /4 (3)
The fact that the friction factor diagram in the form set out by Moody
still endures is due to the use of dimensionless quantities and an effective
consensus to continue to use the constants and underlying equations that
were used by Moody. A consistent skeleton structure is provided. Moody
and his contemporaries were under no illusions about its accuracy. Moody
stated ‘it must be recognized that any high degree of accuracy in determining
f is not to be expected’ and ‘fairly reasonable estimates of friction loss can
be made.’
The diagram is necessarily a simplification and a rough approximation:
the conditions it describes (fully developed, isothermal, incompressible, dis-
sipative, pseudo-steady-state flow) are never quite attained in practice. Con-
ditions within pipelines are inherently non-uniform over the flow cross sec-
tion and over a given length. Fluids are more or less compressible, rather
than incompressible, and have finite thermal conductivity. The flowing fluid
within pipelines is not isothermal in the radial direction because of the very
dissipation that is quantified by the Darcy-Weisbach equation. Frictional
dissipation and finite thermal conductivity together give rise to differences
in temperature—and density and viscosity—between different parts of the
flow. Inevitably too, any length of pipe for which the fully-developed flow
requirement is approximated must be connected to entry and exit systems
within which this is not the case. When the flow is turbulent, vortices form
and collapse relentlessly over time. In the flow of liquids, vapour pressure
may play a role. Nevertheless, the diagram is a very useful design tool. For
computational purposes, it is easily represented as a data set of discrete
points, or by the equations that define it along with appropriate solution
algorithms.
The preparation of new diagrams by the present author was prompted by
the need for vector graphics versions of the friction factor diagram that could
be used for teaching purposes and could be distributed freely to students
and others under a Creative Commons Attribution-Share-Alike license. The
same four zones that were labeled by Moody (laminar, critical, transition
and complete turbulence) are shown. Colour has been employed, although
the diagram is also usable in grayscale form.
3
usability at A4 or Letter size. Absolute size is of little relevance for use on a
screen or tablet. The diagrams can be used conveniently on any device that
can display a PDF document.
The mathematical curves have not been plotted directly. Rather, straight
lines have been drawn between precisely calculated data points of the un-
derlying curves. The underlying equations, such as the Colebrook White
[6] equation, Equation 4, were solved or evaluated with a precision of about
fifteen significant digits. The plotted points are sufficiently close together
that visual smoothness is achieved and, as the diagram is a vector graphic,
the file size is small.
!
1 /D 2.51
p = −2 log10 + p (4)
fD 3.7 Re fD
The horizontal axis in the diagram is for the Reynolds number Re, Equa-
tion 5, which is dimensionless. The range has been extended at the upper
end to 109 (Moody’s diagram went to 108 ) and then to infinity!
ρ VD VD
Re = = (5)
µ ν
The vertical axis of the diagram is for the friction factor. The Darcy and
the Fanning friction factors are given by Equations 6 and 7 respectively.
! ! !
L V2 L ρV 2 ρV 2
fD = hf / = ∆p/ = 4τw / (6)
D 2g D 2 2
! ! !
L V2 L ρV 2 ρV 2
fF = hf / 4 = ∆p/ 4 = τw / (7)
D 2g D 2 2
Friction factors were in use prior to the publication of Moody’s diagram.
For example, in 1930 Fred Scobey [7] compiled ‘Weisbach friction coeffi-
cients,’ fD , for riveted and analogous pipes under actual operating condi-
tions. These friction factors were determined from measurements of flow
rate and pressure drop or head loss in situ.
For turbulent flow in smooth pipes, the friction factor is calculated from
the von Kármán, Prandtl, Nikuradse expression, Equation 8.
!
1 2.51
p = −2 log10 p (8)
fD Re fD
4
where /D is the relative roughness of the pipe, being the surface absolute
roughness and D being the inside diameter of the pipe. Equations 8 and
9 are both included in the Colebrook White equation, Equation 4, which
achieves a blending of the two extremes to cover the full transitional and
complete turbulence zones of the diagram. Friction factor curves are shown
on the diagram for relative roughness values ranging from 1 × 10−7 to 0.06.
It can be noted that, in principle, the absolute roughness that is used is
an experimentally determined ‘equivalent sandgrain size’ parameter for the
particular pipe surface. It is referenced to tests carried out by J. Nikuradse
[8].
As part of the discussion, R.J.S. Pigott, who had published a friction factor
diagram with basically the same structure as Moody’s in 1933 [9], proposed
Equation 12 for the boundary curve, while Moody responded with a sugges-
tion that Equation 13 would approximate to a boundary at 1% above the
friction factor for infinite Reynolds number. This has been verified by the
present author. Therefore, Equation 13 also represents the boundary curve
shown in Figure 1.
3500
= Re (12)
/D
1600
= Re (13)
/D
Hence, for a maximum positive correction factor of less than 1% in using
the rough pipes expression (Equation 9) for the friction factor, the complete
turbulence zone can be defined by the inequality labeled Equation 14.
1600
Re > (14)
/D
5
This can be generalized, with excellent precision, for an arbitrary positive
maximum correction factor, η, as the expression labeled Equation 15.
16 × 10−Log10 (η)
Re > (15)
/D
where 2.5 × 10−7 ≤ η ≤ 0.02, 104 ≤ Re ≤ 109 and /D > 0. For example,
setting η = 0.01 yields Equation 13. A practical application of this is to
avoid solving the Colebrook White equation where it is unnecessary. The
correction can be ignored, where it is insignificant, or applied (up to a max-
imum value of 0.02) according to Equations 16 and 17. This correction can
be applied and has a value less than 0.01 in the entire ‘Complete Turbulence
Zone’ of Figure 1.
f = (1 + η)f∞ (16)
η = 10−Log10 (Re D /0.16) (17)
6
for from Equation 4, as in Equation 18. Alternatively, the friction factor
data points could be superimposed on a Moody diagram in order to read the
relative roughness value and then calculate the roughness by multiplying by
the diameter. This was done, for example, in reference [3].
√
!
−1/2 fD 2.51
= 3.7D 10 − p (18)
Re fD
7
important to take account of condition, corrosion or fowling. The summary
table for inclusion on the printed diagram is shown as Table 2.
Relative roughness values depend on the pipe diameter and, for rough
pipes, there could be some ambiguity about which ‘diameter’ to use—a
volume based average would seem the most appropriate, as in Equation 19.
‘Roughness’ would thus extend inside and outside the ‘diameter.’
s
4V
D= (19)
πL
8
Darcy Friction Factor for Pipe Flow Reynolds Number
103 1.4 2 3 4 5 6 104 1.4 2 3 4 5 6 105 1.4 2 3 4 5 6 106 1.4 2 3 4 5 6 107 1.4 2 3 4 5 6 108 1.4 2 3 4 5 6 109 ¥
0.08 0.08
0.06
0.055
0.05
0.07 0.07
0.045
0.04
0.06 0.035 0.06
0.03
0.025
Transition Zone
N
0.05 0.05
rV 2
0.02
2
0.018
0.045 0.016 0.045
= 4 tw J
0.014
0.04 0.012 0.04
0.01
0.035 0.008 0.035
Complete Turbulence Zone
N
0.006
rV 2
2
Friction Factor
0.024 0.024
f f 0.002
0.022 ¥ 0.001 5 0.022
=
1.
01
N
0.02 0.02
Sample 0.001
2g
V2
0.018
inflectional an 0.000 8
0.018
d
Friction Factor fD = hf J DL
roughness 0.000 6
curves Re
0.016 0.000 4 0.016
=
0.015 16 0.000 3 0.015
00
0.014 Critical He 0.014
D
0.000 2
Zone
0.013 L 0.013
0.012 Indicative Roughness Values 0.000 1 0.012
e @mmD
0.011
Smooth honed steel 0.00065
Sm 0.011
Drawn tubing: glass, brass, copper, lead, plastic 0.0015
o oth 0.000 05
0.01 0.01
Asphalted cast iron 0.12 Pip 0.000 03
0.009 Galvanized steel 0.15 e 0.000 02 0.009
Laminar Wood stave 0.18-0.91
Cast iron 0.26
0.008 Zone Concrete 0.3-3
0.000 01 0.008
fD = 64Re Heavy brush coat: asphalts, enamels, tars 0.45-0.6
0.000 005
General tuberculation 1-3 mm 0.6-1.9
0.007 0.007
Riveted steel 0.9-9 0.000 003
Severe tuberculation and incrustation 2.5-6.5 0.000 002
0.006 0.006
0.000 001
5. ´ 10-7
0.005 0.005
103 1.4 2 3 4 5 6 104 1.4 2 3 4 5 6 105 1.4 2 3 4 5 6 106 1.4 2 3 4 5 6 107 1.4 2 3 4 5 6 108 1.4 2 3 4 5 6 109 ¥
rVD VD
Jim McGovern, Dublin Institute of Technology Reynolds Number Re = m
= n
http://www.fun-engineering.net/pipefrictionfactor.html Version DMI001 2011.10.02
Figure 2: Diagram of friction factor for pipe flow, including sample inflec-
tional roughness curves.
9
of 3.66 micron. An equivalent sand grain roughness was not determined for
the smooth coated pipe, which had a measured rms of 1.41 micron.
There is still a scarcity of data on the absolute sand-grain-equivalent
roughness of commercial pipes. The values provided by Farshad et al. [20,
21, 22] are illustrative. An examination of the data available in [21] indicates
that flow measurements were made at a Reynolds number of about 6.5 × 105
and in all cases were within the transitional zone. Also, in [22] the authors
seem to have misinterpreted Moody’s chart (his Figure 2, p. 678) that related
relative roughness to pipe diameter for various absolute roughness values—
Moody’s log-log chart of /D versus D was simply based on the inverse
proportionality between relative roughness and pipe diameter.
Conclusions
The diagrams have been prepared by the author for students of engineering
and practicing engineers. He would welcome any suggestions for further
adjustments.
Acknowledgement
This work was supported in part by a National Digital Learning Resources
and Dublin Institute of Technology Award within the 2010–2011 academic
year.
Nomenclature
∆p pressure drop due to fluid friction
η correction factor
V volume
10
Re Reynolds number
µ absolute viscosity
ν kinematic viscosity
ρ density
D pipe diameter
f friction factor
L pipe length
11
Table 1: Indicative absolute roughness values
Ref. /[mm]
Finely honed smooth steel
(rms = 0.00015 mm, inflectional) 11 (2007) 0.00065
Drawn tubing (e.g. glass, brass,
copper, lead, plastic) 1 (1944) 0.0015
Butt-welded steel: new smooth pipe
or with centrifugally applied enamels 3 (1977) 0.009 - 0.06
Honed steel
(rms = 0.0025 mm, inflectional) 11 (2007) 0.01
(rms = 0.0037 mm, inflectional) 17 (2003) 0.021
Butt-welded steel: centrifugally
applied concrete linings 3 (1977) 0.045 - 0.15
Wrought iron, steel 1 (1944) 0.046
Butt-welded steel:
hot asphalt dipped 3 (1977) 0.06 - 0.15
Asphalted cast iron 1 (1944) 0.12
Galvanized iron 1 (1944) 0.15
Butt-welded steel: light rust 3 (1977) 0.15 - 0.35
Wood stave 1 (1944) 0.18 - 0.91
Cast iron 1 (1944) 0.26
Concrete 1 (1944) 0.3 - 3
Butt-welded steel: heavy brush coat,
asphalts, enamels and tars 3 (1977) 0.45 - 0.6
Riveted steel 1 (1944) 0.9 - 9
Butt-welded steel:
general tuberculation 1-3 mm 3 (1977) 1 - 1.85
severe tuberculation
and incrustation 3 (1977) 2.5 - 6.5
12
Table 2: Summary table of indicative absolute roughness values
/[mm]
Smooth honed steel 0.00065
Drawn tubing: glass, brass, copper, lead, plastic 0.0015
Asphalted cast iron 0.12
Galvanized steel 0.15
Wood stave 0.18 - 0.91
Cast iron 0.26
Concrete 0.3 - 3
Heavy brush coat: asphalts, enamels, tars 0.45 - 0.6
General tuberculation 1-3 mm 0.6 - 1.9
Riveted steel 0.9 - 9
Severe tuberculation and incrustation 2.5 - 6.5
References
[1] Lewis F. Moody. Friction factors for pipe flow. Transactions of the
A.S.M.E., pages 671–684, November 1944.
[2] Marco Lorenzini, Gian Luca Morini, and Sandro Salvigni. Lami-
nar, transitional and turbulent friction factors for gas flows in smooth
and rough microtubes. International Journal of Thermal Sciences,
49(2):248–255, 2010.
[5] Glenn O. Brown. The history of the Darcy-Weisbach equation for pipe
flow resistance. In A. Fredrich and J. Rogers, editors, Proceedings of
the 150th Anniversary Conference of ASCE, Washington D.C., pages
34–43. American Society of Civil Engineers, November 3-6 2002. Avail-
able at http://bioen.okstate.edu/Darcy/DarcyWeisbach/
HistoryoftheDarcyWeisbachEq.pdf, accessed 2 August 2011.
13
the transition region between the smooth and rough pipe laws. Journal
of the Institution of Civil Engineers, 11:133–156, 1938-1939.
[7] Fred C. Scobey. The flow of water in riveted steel and analogous pipes.
United States Department of Agriculture, January 1930. Available at
http://naldr.nal.usda.gov/, accessed 8 August 2011.
[9] R.J.S. Pigott. The flow of fluids in closed conduits. Mechanical Engi-
neering, 55:497–501,515, August 1933.
[10] B.J. McKeon, C.J. Swanson, M.V Zagarola, R.J. Donnelly, and A.J.
Smits. Friction factors for smooth pipe flow. Journal of Fluid Mechan-
ics, 511:41–44, 2004.
[11] J.J. Allen, M.A. Shockling, G.J. Kunkel, and A.J. Smits.
Turbulent flow in smooth and rough pipes. Philosophical
Transactions. Series A, Mathematical, physical, and engineer-
ing sciences, 365(1852):699–714, March 2007. Available at
http://rsta.royalsocietypublishing.org/cgi/content/
abstract/365/1852/699, accessed 22 August 2011.
[12] B.J. McKeon, M.V. Zagarola, and A.J. Smits. A new friction factor
relationship for fully developed pipe flow. Journal of Fluid Mechanics,
538:429–443, 2005.
[14] C.F. Colebrook and C.M. White. Experiments with fluid friction in
roughened pipes. In Proceedings of the Royal Society of London A,
Volume 161, pages 367–381. The Royal Society, February 1937.
14
[17] Elling Sletfjerding and Jon Steinar Gudmundsson. Friction factor di-
rectly from roughness measurements. Journal of Energy Resources
Technology, Transactions of the ASME, 125(2):126–130, 2003.
[18] J.J. Allen, M.A. Shockling, and A.J. Smits. Evaluation of a universal
transitional resistance diagram for pipes with honed surfaces. Physics
of Fluids, 17(12):121702, December 2005.
[19] M.A. Shockling, J.J. Allen, and A.J. Smits. Roughness effects in tur-
bulent pipe flow. Journal of Fluid Mechanics, 564(1):267–285, 2006.
[20] Fred Farshad, Herman Rieke, and James Garber. New developments in
surface roughness measurements, characterization, and modeling fluid
flow in pipe. Journal of Petroleum Science and Engineering, 29:139–
150, February 2001.
[21] Fred F. Farshad, H.H. Rieke, and Charlie Mauldin. Flow test validation
of direct measurement methods used to determine surface roughness in
pipes (OCTG). In SPE Western Regional/AAPG Pacific Section Joint
Meeting, pages 529–540, Anchorage, AK, United States, May 20–22
2002.
[22] F.F. Farshad and H.H. Rieke. Surface-roughness design values for mod-
ern pipes. SPE Drilling & Completion, pages 212–215, September 2006.
15