Advanced iterative procedures for solving the implicit Colebrook equation for fluid flow friction

Abstract

The empirical Colebrook equation from 1939 is still accepted as an informal standard way to calculate the friction factor of turbulent flows (4000 < Re < 10(8)) through pipes with roughness between negligible relative roughness (epsilon/D -> 0) to very rough (up to epsilon/D = 0.05). The Colebrook equation includes the flow friction factor lambda in an implicit logarithmic form, lambda being a function of the Reynolds number Re and the relative roughness of inner pipe surface epsilon/D: lambda = f(lambda, Re, epsilon/D). To evaluate the error introduced by the many available explicit approximations to the Colebrook equation, lambda approximate to f(Re, epsilon/D), it is necessary to determinate the value of the friction factor lambda from the Colebrook equation as accurately as possible. The most accurate way to achieve that is by using some kind of the iterative method. The most used iterative approach is the simple fixed-point method, which requires up to 10 iterations to achieve a good level of accuracy. The simple fixed-point method does not require derivatives of the Colebrook function, while the most of the other presented methods in this paper do require. The methods based on the accelerated Householder's approach (3rd order, 2nd order: Halley's and Schroder's method, and 1st order: Newton-Raphson) require few iterations less, while the three-point iterative methods require only 1 to 3 iterations to achieve the same level of accuracy. The paper also discusses strategies for finding the derivatives of the Colebrook function in symbolic form, for avoiding the use of the derivatives (secant method), and for choosing an optimal starting point for the iterative procedure. The Householder approach to the Colebrook' equations expressed through the Lambert W-function is also analyzed. Finally, it is presented one approximation to the Colebrook equation with an error of no more than 0.0617%.