scholarly journals Evolutionary Optimization of Colebrook’s Turbulent Flow Friction Approximations

Author(s):  
Dejan Brkić ◽  
Žarko Ćojbašić

Today, Colebrook’s equation is mostly accepted as an informal standard for modeling of turbulent flow in hydraulically smooth and rough pipes including transient zone in between. The empirical Colebrook’s equation relates the unknown flow friction factor (λ) with the known Reynolds number (R) and the known relative roughness of inner pipe surface (ε/D). It is implicit in unknown friction factor (λ). Implicit Colebrook’s equation cannot be rearranged to derive friction factor (λ) directly and therefore it can be solved only iteratively [λ=f(λ, R, ε/D)] or using its explicit approximations [λ≈f(R, ε/D)]. Of course, approximations carry in certain error compared with the iterative solution where the highest level of accuracy can be reached after enough number of iterations. The explicit approximations give a relatively good prediction of the friction factor (λ) and can reproduce accurately Colebrook’s equation and its Moody’s plot. Usually, more complex models of approximations are more accurate and vice versa. In this paper, numerical values of parameters in various existing approximations are changed (optimized) using genetic algorithms to reduce maximal relative error. After this improvement computational burden stays unchanged while accuracy of approximations increases in some of the cases very significantly.

2017 ◽  
Author(s):  
Dejan Brkić ◽  
Žarko Ćojbašić

This paper presents evolutionary optimization of explicit approximations of the empirical Colebrook’s equation that is used for the calculation of the turbulent friction factor (λ), i.e., for the calculation of turbulent hydraulic resistance in hydraulically smooth and rough pipes including the transient zone between them. The empirical Colebrook’s equation relates the unknown flow friction factor (λ) with the known Reynolds number (R) and the known relative roughness of the inner pipe surface (ε/D). It is implicit in the unknown friction factor (λ). The implicit Colebrook’s equation cannot be rearranged to derive the friction factor (λ) directly, and therefore, it can be solved only iteratively [λ = f(λ, R, ε/D)] or using its explicit approximations [λ ≈ f(R, ε/D)], which introduce certain error compared with the iterative solution. The optimization of explicit approximations of Colebrook’s equation is performed with the aim to improve their accuracy, and the proposed optimization strategy is demonstrated on a large number of explicit approximations published up to date where numerical values of the parameters in various existing approximations are changed (optimized) using genetic algorithms to reduce maximal relative error. After that improvement, the computational burden stays unchanged while the accuracy of approximations increases in some of the cases very significantly.


2017 ◽  
Author(s):  
Dejan Brkić

Empirical Colebrook equation implicit in unknown ow friction factor (λ) is an accepted standard for calculation of hydraulic resistance in hydraulically smooth and rough pipes. e Colebrook equation gives friction factor (λ) implicitly as a function of the Reynolds number (Re) and relative roughness (ε/D) of inner pipe surface; i.e. λ0=f(λ0, Re, ε/D). e paper presents a problem that requires iterative methods for the solution. In particular, the implicit method used for calculating the friction factor λ0 is an application of xed- point iterations. e type of problem discussed in this "in the classroom paper" is commonly encountered in uid dynamics, and this paper provides readers with the tools necessary to solve similar problems. Students’ task is to solve the equation using Excel where the procedure for that is explained in this “in the classroom” paper. Also, up to date numerous explicit approximations of the Colebrook equation are available where as an additional task for students can be evaluation of the error introduced by these explicit approximations λ≈f(Re, ε/D) compared with the iterative solution of implicit equation which can be treated as accurate.


2017 ◽  
Author(s):  
Dejan Brkić

Two new correlations of single-phase friction factor for turbulent pipe flow are shown in this paper. These two formulas are actually explicit approximations of iterative Colebrook's relation for calculation of flow friction factor. Calculated friction factors are valid for whole turbulent flow including hydraulically smooth and rough pipes with special attention on transient zone of turbulence between them. Hydraulically smooth regime of turbulence does not occur only in total absence of roughness of inner pipe surface, but also, four new relations for this theoretical regime are presented. Some recent formulas for turbulent flow friction calculation are also commented.


Water ◽  
2018 ◽  
Vol 10 (9) ◽  
pp. 1175 ◽  
Author(s):  
Pavel Praks ◽  
Dejan Brkić

Widely used in hydraulics, the Colebrook equation for flow friction relates implicitly to the input parameters; the Reynolds number, Re and the relative roughness of an inner pipe surface, ε/D with an unknown output parameter; the flow friction factor, λ; λ = f (λ, Re, ε/D). In this paper, a few explicit approximations to the Colebrook equation; λ ≈ f (Re, ε/D), are generated using the ability of artificial intelligence to make inner patterns to connect input and output parameters in an explicit way not knowing their nature or the physical law that connects them, but only knowing raw numbers, {Re, ε/D}→{λ}. The fact that the used genetic programming tool does not know the structure of the Colebrook equation, which is based on computationally expensive logarithmic law, is used to obtain a better structure of the approximations, which is less demanding for calculation but also enough accurate. All generated approximations have low computational cost because they contain a limited number of logarithmic forms used for normalization of input parameters or for acceleration, but they are also sufficiently accurate. The relative error regarding the friction factor λ, in in the best case is up to 0.13% with only two logarithmic forms used. As the second logarithm can be accurately approximated by the Padé approximation, practically the same error is obtained also using only one logarithm.


Energies ◽  
2018 ◽  
Vol 11 (7) ◽  
pp. 1825 ◽  
Author(s):  
Pavel Praks ◽  
Dejan Brkić

The 80 year-old empirical Colebrook function ξ, widely used as an informal standard for hydraulic resistance, relates implicitly the unknown flow friction factor λ, with the known Reynolds number Re and the known relative roughness of a pipe inner surface ε*; λ=ξ(Re,ε*,λ). It is based on logarithmic law in the form that captures the unknown flow friction factor λ in a way that it cannot be extracted analytically. As an alternative to the explicit approximations or to the iterative procedures that require at least a few evaluations of computationally expensive logarithmic function or non-integer powers, this paper offers an accurate and computationally cheap iterative algorithm based on Padé polynomials with only one log-call in total for the whole procedure (expensive log-calls are substituted with Padé polynomials in each iteration with the exception of the first). The proposed modification is computationally less demanding compared with the standard approaches of engineering practice, but does not influence the accuracy or the number of iterations required to reach the final balanced solution.


2018 ◽  
Vol 2018 ◽  
pp. 1-18 ◽  
Author(s):  
Pavel Praks ◽  
Dejan Brkić

The empirical Colebrook equation from 1939 is still accepted as an informal standard way to calculate the friction factor of turbulent flows (4000 < Re < 108) through pipes with roughness between negligible relative roughness (ε/D ⟶ 0) to very rough (up to ε/D = 0.05). The Colebrook equation includes the flow friction factor λ in an implicit logarithmic form, λ being a function of the Reynolds number Re and the relative roughness of inner pipe surface ε/D: λ = f(λ, Re, ε/D). To evaluate the error introduced by the many available explicit approximations to the Colebrook equation, λ ≈ f(Re, ε/D), it is necessary to determinate the value of the friction factor λ 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 Schröder’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%.


Author(s):  
Pavel Praks ◽  
Dejan Brkić

Widely used in hydraulics, the Colebrook equation for flow friction relates implicitly to the input parameters; the Reynolds number, Re and the relative roughness of inner pipe surface, &epsilon;/D with the output unknown parameter; the flow friction factor, &lambda;; &lambda;=f(&lambda;, Re, &epsilon;/D). In this paper, a few explicit approximations to the Colebrook equation; &lambda;&asymp;f(Re, &epsilon;/D), are generated using the ability of artificial intelligence to make inner patterns to connect input and output parameters in explicit way not knowing their nature or the physical law that connects them, but only knowing raw numbers, {Re, &epsilon;/D}&rarr;{&lambda;}. The fact that the used genetic programming tool does not know the structure of the Colebrook equation which is based on computationally expensive logarithmic law, is used to obtain better structure of the approximations which is less demanding for calculation but also enough accurate. All generated approximations are with low computational cost because they contain a limited number of logarithmic forms used although for normalization of input parameters or for acceleration, but they are also sufficiently accurate. The relative error regarding the friction factor &lambda;, in best case is up to 0.13% with only two logarithmic forms used. As the second logarithm can be accurately approximated by the Pad&eacute; approximation, practically the same error is obtained also using only one logarithm.


2017 ◽  
Author(s):  
Žarko Ćojbašić ◽  
Dejan Brkić

To date, the Colebrook equation is mostly accepted as an unofficial standard for calculation of the friction factor in turbulent flow through pipes. Unfortunately, the unknown friction factor in the Colebrook equation is given implicitly. Therefore, the implicit Colebrook equation has to be solved in an iterative procedure or using some of the appropriate explicit correlations proposed by many authors. Although the iterative solution is simple and very accurate, it can cause some problems during the calculation of looped network of pipes or similar systems of pipes. Therefore, explicit approximations are favorable in these cases. Up to date, the most accurate approximations have maximal relative error of no more than 0.14% compared to the very accurate iterative solution. Here two explicit approximations are presented, based on already existing models which are improved using genetic algorithms optimization. They are with the maximal relative error of no more than 0.0083% and 0.0026%.


2017 ◽  
Author(s):  
Dejan Brkić ◽  
Brkic Dejan

Maximal relative error of the explicit approximation to the Colebrook equation for flow friction presented in the discussed paper by Saeed Samadianfard [J. Pet. Sci. Eng. 92-93 (2012), 48-55; doi. 10.1016/j.petrol.2012.06.005] is investigated. Samadianfard claims that his approximation is very accurate with the maximal relative error of no more than 0.08152%. Here is shown that this error is about 7%. Related comments about the paper are also enclosed. ; JRC.F.3-Energy Security, Systems and Market


Author(s):  
Pavel Praks ◽  
Dejan Brkić

Empirical Colebrook equation from 1939 is still accepted as an informal standard to calculate friction factor during the turbulent flow (4000 &lt; Re &lt; 108) through pipes from smooth with almost negligible relative roughness (&epsilon;/D&rarr;0) to the very rough (up to &epsilon;/D = 0.05) inner surface. The Colebrook equation contains flow friction factor &lambda; in implicit logarithmic form where it is, aside of itself; &lambda;, 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 many available explicit approximations to the Colebrook equation, &lambda; &asymp; f(Re, &epsilon;/D), it is necessary to determinate value of the friction factor &lambda; from the Colebrook equation as accurate as possible. The most accurate way to achieve that is using some kind of iterative methods. Usually classical approach also known as simple fixed point method requires up to 8 iterations to achieve the high level of accuracy, but does not require derivatives of the Colebrook function as here presented accelerated Householder&rsquo;s approach (3rd order, 2nd order: Halley&rsquo;s and Schr&ouml;der&rsquo;s method and 1st order: Newton-Raphson) which needs only 3 to 7 iteration and three-point iterative methods which needs only 1 to 4 iteration to achieve the same high level of accuracy. Strategies how to find derivatives of the Colebrook function in symbolic form, how to avoid use of the derivatives (Secant method) and how to choose optimal starting point for the iterative procedure are shown. Householder&rsquo;s approach to the Colebrook&rsquo;s equations expressed through the Lambert W-function is also analyzed. One approximation to the Colebrook equation based on the analysis from the paper with the error of no more than 0.0617% is shown.


Sign in / Sign up

Export Citation Format

Share Document