Jump to content

Darcy friction factor formulae

From Wikipedia, the free encyclopedia
(Redirected from Colebrook equation)

In fluid dynamics, the Darcy friction factor formulae are equations that allow the calculation of the Darcy friction factor, a dimensionless quantity used in the Darcy–Weisbach equation, for the description of friction losses in pipe flow as well as open-channel flow.

The Darcy friction factor is also known as the Darcy–Weisbach friction factor, resistance coefficient or simply friction factor; by definition it is four times larger than the Fanning friction factor.[1]

Notation

[edit]

In this article, the following conventions and definitions are to be understood:

  • The Reynolds number Re is taken to be Re = V D / ν, where V is the mean velocity of fluid flow, D is the pipe diameter, and where ν is the kinematic viscosity μ / ρ, with μ the fluid's Dynamic viscosity, and ρ the fluid's density.
  • The pipe's relative roughness ε / D, where ε is the pipe's effective roughness height and D the pipe (inside) diameter.
  • f stands for the Darcy friction factor. Its value depends on the flow's Reynolds number Re and on the pipe's relative roughness ε / D.
  • The log function is understood to be base-10 (as is customary in engineering fields): if x = log(y), then y = 10x.
  • The ln function is understood to be base-e: if x = ln(y), then y = ex.

Flow regime

[edit]

Which friction factor formula may be applicable depends upon the type of flow that exists:

  • Laminar flow
  • Transition between laminar and turbulent flow
  • Fully turbulent flow in smooth conduits
  • Fully turbulent flow in rough conduits
  • Free surface flow.

Transition flow

[edit]

Transition (neither fully laminar nor fully turbulent) flow occurs in the range of Reynolds numbers between 2300 and 4000. The value of the Darcy friction factor is subject to large uncertainties in this flow regime.

Turbulent flow in smooth conduits

[edit]

The Blasius correlation is the simplest equation for computing the Darcy friction factor. Because the Blasius correlation has no term for pipe roughness, it is valid only to smooth pipes. However, the Blasius correlation is sometimes used in rough pipes because of its simplicity. The Blasius correlation is valid up to the Reynolds number 100000.

Turbulent flow in rough conduits

[edit]

The Darcy friction factor for fully turbulent flow (Reynolds number greater than 4000) in rough conduits can be modeled by the Colebrook–White equation.

Free surface flow

[edit]

The last formula in the Colebrook equation section of this article is for free surface flow. The approximations elsewhere in this article are not applicable for this type of flow.

Choosing a formula

[edit]

Before choosing a formula it is worth knowing that in the paper on the Moody chart, Moody stated the accuracy is about ±5% for smooth pipes and ±10% for rough pipes. If more than one formula is applicable in the flow regime under consideration, the choice of formula may be influenced by one or more of the following:

  • Required accuracy
  • Speed of computation required
  • Available computational technology:
    • calculator (minimize keystrokes)
    • spreadsheet (single-cell formula)
    • programming/scripting language (subroutine).

Colebrook–White equation

[edit]

The phenomenological Colebrook–White equation (or Colebrook equation) expresses the Darcy friction factor f as a function of Reynolds number Re and pipe relative roughness ε / Dh, fitting the data of experimental studies of turbulent flow in smooth and rough pipes.[2][3] The equation can be used to (iteratively) solve for the Darcy–Weisbach friction factor f.

For a conduit flowing completely full of fluid at Reynolds numbers greater than 4000, it is expressed as:

or

where:

  • Hydraulic diameter, (m, ft) – For fluid-filled, circular conduits, = D = inside diameter
  • Hydraulic radius, (m, ft) – For fluid-filled, circular conduits, = D/4 = (inside diameter)/4

Note: Some sources use a constant of 3.71 in the denominator for the roughness term in the first equation above.[4]

Solving

[edit]

The Colebrook equation is usually solved numerically due to its implicit nature. Recently, the Lambert W function has been employed to obtain explicit reformulation of the Colebrook equation.[5][6][7]

or

will get:

then:

Expanded forms

[edit]

Additional, mathematically equivalent forms of the Colebrook equation are:

where:
1.7384... = 2 log (2 × 3.7) = 2 log (7.4)
18.574 = 2.51 × 3.7 × 2

and

or
where:
1.1364... = 1.7384... − 2 log (2) = 2 log (7.4) − 2 log (2) = 2 log (3.7)
9.287 = 18.574 / 2 = 2.51 × 3.7.

The additional equivalent forms above assume that the constants 3.7 and 2.51 in the formula at the top of this section are exact. The constants are probably values which were rounded by Colebrook during his curve fitting; but they are effectively treated as exact when comparing (to several decimal places) results from explicit formulae (such as those found elsewhere in this article) to the friction factor computed via Colebrook's implicit equation.

Equations similar to the additional forms above (with the constants rounded to fewer decimal places, or perhaps shifted slightly to minimize overall rounding errors) may be found in various references. It may be helpful to note that they are essentially the same equation.

Free surface flow

[edit]

Another form of the Colebrook-White equation exists for free surfaces. Such a condition may exist in a pipe that is flowing partially full of fluid. For free surface flow:

The above equation is valid only for turbulent flow. Another approach for estimating f in free surface flows, which is valid under all the flow regimes (laminar, transition and turbulent) is the following:[8]

where a is:

and b is:

where Reh is Reynolds number where h is the characteristic hydraulic length (hydraulic radius for 1D flows or water depth for 2D flows) and Rh is the hydraulic radius (for 1D flows) or the water depth (for 2D flows). The Lambert W function can be calculated as follows:

Approximations of the Colebrook equation

[edit]

Haaland equation

[edit]

The Haaland equation was proposed in 1983 by Professor S.E. Haaland of the Norwegian Institute of Technology.[9] It is used to solve directly for the Darcy–Weisbach friction factor f for a full-flowing circular pipe. It is an approximation of the implicit Colebrook–White equation, but the discrepancy from experimental data is well within the accuracy of the data.

The Haaland equation[10] is expressed:

Swamee–Jain equation

[edit]

The Swamee–Jain equation is used to solve directly for the Darcy–Weisbach friction factor f for a full-flowing circular pipe. It is an approximation of the implicit Colebrook–White equation.[11]

Serghides's solution

[edit]

Serghides's solution is used to solve directly for the Darcy–Weisbach friction factor f for a full-flowing circular pipe. It is an approximation of the implicit Colebrook–White equation. It was derived using Steffensen's method.[12]

The solution involves calculating three intermediate values and then substituting those values into a final equation.

The equation was found to match the Colebrook–White equation within 0.0023% for a test set with a 70-point matrix consisting of ten relative roughness values (in the range 0.00004 to 0.05) by seven Reynolds numbers (2500 to 108).

Goudar–Sonnad equation

[edit]

Goudar equation is the most accurate approximation to solve directly for the Darcy–Weisbach friction factor f for a full-flowing circular pipe. It is an approximation of the implicit Colebrook–White equation. Equation has the following form[13]

Brkić solution

[edit]

Brkić shows one approximation of the Colebrook equation based on the Lambert W-function[14]

The equation was found to match the Colebrook–White equation within 3.15%.

Brkić-Praks solution

[edit]

Brkić and Praks show one approximation of the Colebrook equation based on the Wright -function, a cognate of the Lambert W-function[15]

, , , and

The equation was found to match the Colebrook–White equation within 0.0497%.

Praks-Brkić solution

[edit]

Praks and Brkić show one approximation of the Colebrook equation based on the Wright -function, a cognate of the Lambert W-function[16]

, , , and

The equation was found to match the Colebrook–White equation within 0.0012%.

Niazkar's solution

[edit]

Since Serghides's solution was found to be one of the most accurate approximation of the implicit Colebrook–White equation, Niazkar modified the Serghides's solution to solve directly for the Darcy–Weisbach friction factor f for a full-flowing circular pipe.[17]

Niazkar's solution is shown in the following:

Niazkar's solution was found to be the most accurate correlation based on a comparative analysis conducted in the literature among 42 different explicit equations for estimating Colebrook friction factor.[17]

Blasius correlations

[edit]

Early approximations for smooth pipes[18] by Paul Richard Heinrich Blasius in terms of the Darcy–Weisbach friction factor are given in one article of 1913:[19]

.

Johann Nikuradse in 1932 proposed that this corresponds to a power law correlation for the fluid velocity profile.[20]

Mishra and Gupta in 1979 proposed a correction for curved or helically coiled tubes, taking into account the equivalent curve radius, Rc:[21]

,

with,

where f is a function of:

  • Pipe diameter, D (m, ft)
  • Curve radius, R (m, ft)
  • Helicoidal pitch, H (m, ft)
  • Reynolds number, Re (dimensionless)

valid for:

  • Retr < Re < 105
  • 6.7 < 2Rc/D < 346.0
  • 0 < H/D < 25.4

Swamee equation

[edit]

The Swamee equation is used to solve directly for the Darcy–Weisbach friction factor (f) for a full-flowing circular pipe for all flow regimes (laminar, transitional, turbulent). It is an exact solution for the Hagen–Poiseuille equation in the laminar flow regime and an approximation of the implicit Colebrook–White equation in the turbulent regime with a maximum deviation of less than 2.38% over the specified range. Additionally, it provides a smooth transition between the laminar and turbulent regimes to be valid as a full-range equation, 0 < Re < 108. [22]

Table of Approximations

[edit]

The following table lists historical approximations to the Colebrook–White relation[23] for pressure-driven flow. Churchill equation[24] (1977) is the only equation that can be evaluated for very slow flow (Reynolds number < 1), but the Cheng (2008),[25] and Bellos et al. (2018)[8] equations also return an approximately correct value for friction factor in the laminar flow region (Reynolds number < 2300). All of the others are for transitional and turbulent flow only.

Table of Colebrook equation approximations
Equation Author Year Range Ref

Moody 1947

where
Wood 1966

Eck 1973

Swamee and Jain 1976

Churchill 1973

Jain 1976

where
Churchill 1977 All flow regimes

Chen 1979

Round 1980

Barr 1981

or

Zigrang and Sylvester 1982

Haaland[10] 1983

or

where
Serghides 1984

if then and if then

Tsal 1989 [26]

Manadilli 1997

Romeo, Royo, Monzon 2002

where:
Goudar, Sonnad 2006

where:
Vatankhah, Kouchakzadeh 2008

where
Buzzelli 2008

where


Cheng 2008 All flow regimes [25]

Avci, Kargoz 2009

Evangelides, Papaevangelou, Tzimopoulos 2010

Fang 2011

,

Brkić 2011

where
S.Alashkar 2012

where

Bellos, Nalbantis, Tsakiris 2018 All flow regimes [8][27]

where

Niazkar 2019 [17]
Tkachenko, Mileikovskyi 2020 Deviation 5.36 %,

[28]

where

Tkachenko, Mileikovskyi 2020 Deviation 0.00072 %,

[28]

References

[edit]
  1. ^ Manning, Francis S.; Thompson, Richard E. (1991). Oilfield Processing of Petroleum. Vol. 1: Natural Gas. PennWell Books. ISBN 978-0-87814-343-6., 420 pages. See page 293.
  2. ^ Colebrook, C. F.; White, C. M. (1937). "Experiments with Fluid Friction in Roughened Pipes". Proceedings of the Royal Society of London. Series A, Mathematical and Physical Sciences. 161 (906): 367–381. Bibcode:1937RSPSA.161..367C. doi:10.1098/rspa.1937.0150. Often erroneously cited as the source of the Colebrook-White equation. This is partly because Colebrook (in a footnote in his 1939 paper) acknowledges his debt to White for suggesting the mathematical method by which the smooth and rough pipe correlations could be combined.
  3. ^ Colebrook, C F (1939). "Turbulent Flow in Pipes, with Particular Reference to the Transition Region Between the Smooth and Rough Pipe Laws". Journal of the Institution of Civil Engineers. 11 (4): 133–156. doi:10.1680/ijoti.1939.13150. ISSN 0368-2455.
  4. ^ VDI Gesellschaft (2010). VDI Heat Atlas. Springer. ISBN 978-3-540-77876-9.
  5. ^ More, A. A. (2006). "Analytical solutions for the Colebrook and White equation and for pressure drop in ideal gas flow in pipes". Chemical Engineering Science. 61 (16): 5515–5519. Bibcode:2006ChEnS..61.5515M. doi:10.1016/j.ces.2006.04.003.
  6. ^ Brkić, D. (2012). "Lambert W Function in Hydraulic Problems" (PDF). Mathematica Balkanica. 26 (3–4): 285–292.
  7. ^ Keady, G. (1998). "Colebrook-White Formula for Pipe Flows". Journal of Hydraulic Engineering. 124 (1): 96–97. CiteSeerX 10.1.1.1027.8918. doi:10.1061/(ASCE)0733-9429(1998)124:1(96).
  8. ^ a b c Bellos, Vasilis; Nalbantis, Ioannis; Tsakiris, George (December 2018). "Friction Modeling of Flood Flow Simulations". Journal of Hydraulic Engineering. 144 (12): 04018073. doi:10.1061/(asce)hy.1943-7900.0001540. ISSN 0733-9429.
  9. ^ Haaland, SE (1983). "Simple and Explicit Formulas for the Friction Factor in Turbulent Flow". Journal of Fluids Engineering. 105 (1): 89–90. doi:10.1115/1.3240948.
  10. ^ a b Massey, Bernard Stanford (1989). Mechanics of fluids. Chapman & Hall. ISBN 978-0-412-34280-6.
  11. ^ Swamee, P.K.; Jain, A.K. (1976). "Explicit equations for pipe-flow problems". Journal of the Hydraulics Division. 102 (5): 657–664. doi:10.1061/JYCEAJ.0004542.
  12. ^ T.K, Serghides (1984). "Estimate friction factor accurately". Chemical Engineering Journal. 91 (5): 63–64. ISSN 0009-2460.
  13. ^ Goudar, C. T; Sonnad, J. R. (2008). "Comparison of the iterative approximations of the Colebrook-White equation: Here's a review of other formulas and a mathematically exact formulation that is valid over the entire range of Re values". Hydrocarbon Processing. 87 (8).
  14. ^ Brkić, Dejan (2011). "An Explicit Approximation of Colebrook's equation for fluid flow friction factor" (PDF). Petroleum Science and Technology. 29 (15): 1596–1602. doi:10.1080/10916461003620453. S2CID 97080106.
  15. ^ Brkić, Dejan; Praks, Pavel (2019). "Accurate and Efficient Explicit Approximations of the Colebrook Flow Friction Equation Based on the Wright ω-Function". Mathematics. 7 (1): 34. doi:10.3390/math7010034. hdl:10084/134214.
  16. ^ Praks, Pavel; Brkić, Dejan (2020). "Review of new flow friction equations: Constructing Colebrook's explicit correlations accurately". Revista Internacional de Métodos Numéricos para Cálculo y Diseño en Ingeniería. 36 (3). arXiv:2005.07021. doi:10.23967/j.rimni.2020.09.001.
  17. ^ a b c Majid, Niazkar (2019). "Revisiting the Estimation of Colebrook Friction Factor: A Comparison between Artificial Intelligence Models and C-W based Explicit Equations". KSCE Journal of Civil Engineering. 23 (10): 4311–4326. doi:10.1007/s12205-019-2217-1. S2CID 203040860.
  18. ^ Massey, B. S. (2006). Mechanics of Fluids (8th ed.). Taylor & Francis. p. 254 eq 7.5. ISBN 978-0-415-36205-4.
  19. ^ Trinh, Khanh Tuoc (2010), On the Blasius correlation for friction factors, arXiv:1007.2466, Bibcode:2010arXiv1007.2466T
  20. ^ Nikuradse, Johann (1932). "Gesetzmässigkeiten der Turbulenten Stromung in Glatten Rohren". VDI Forschungsheft. 359 B (3). Verein Deutscher Ingenieure: 1–36.
  21. ^ Bejan, Adrian; Kraus, Allan D. (2003). Heat Transfer Handbook. John Wiley & Sons. ISBN 978-0-471-39015-2.
  22. ^ Swamee, P.K. (1993). "Design of a Submarine Oil Pipeline". Journal of Transportation Engineering. 119 (1). doi:10.1061/(ASCE)0733-947X(1993)119:1(159).
  23. ^ Brkić, Dejan (March 2012). "Determining Friction Factors in Turbulent Pipe Flow". Chemical Engineering. Beograd: 34–39.(subscription required)
  24. ^ Churchill, S.W. (November 7, 1977). "Friction-factor equation spans all fluid-flow regimes". Chemical Engineering: 91–92.
  25. ^ a b Cheng, Nian-Sheng (September 2008). "Formulas for Friction Factor in Transitional Regimes". Journal of Hydraulic Engineering. 134 (9): 1357–1362. doi:10.1061/(asce)0733-9429(2008)134:9(1357). hdl:10220/7647. ISSN 0733-9429.
  26. ^ Zeyu, Zhang; Junrui, Chai; Zhanbin, Li; Zengguang, Xu; Peng, Li (2020-06-01). "Approximations of the Darcy–Weisbach friction factor in a vertical pipe with full flow regime". Water Supply. 20 (4): 1321–1333. doi:10.2166/ws.2020.048. ISSN 1606-9749.
  27. ^ Bellos, Vasilis; Nalbantis, Ioannis; Tsakiris, George (2020-10-01). "Erratum for "Friction Modeling of Flood Flow Simulations" by Vasilis Bellos, Ioannis Nalbantis, and George Tsakiris". Journal of Hydraulic Engineering. 146 (10): 08220005. doi:10.1061/(ASCE)HY.1943-7900.0001802. ISSN 1943-7900.
  28. ^ a b Mileikovskyi, Viktor; Tkachenko, Tetiana (2020-08-17). "Precise Explicit Approximations of the Colebrook-White Equation for Engineering Systems". Proceedings of EcoComfort 2020. Lecture Notes in Civil Engineering. Vol. 100. pp. 303–310. doi:10.1007/978-3-030-57340-9_37. ISBN 978-3-030-57339-3. ISSN 2366-2557. S2CID 224859478.(subscription required)

Further reading

[edit]
[edit]