Abstract

The empirical logarithmic Colebrook equation for hydraulic resistance in pipes implicitly considers the unknown flow friction factor. Its explicit approximations, used to avoid iterative computations, should be accurate but also computationally efficient. We present a rational approximate procedure that completely avoids the use of transcendental functions, such as logarithm or non-integer power, which require execution of the additional number of floating-point operations in computer processor units. Instead of these, we use only rational expressions that are executed directly in the processor unit. The rational approximation was found using a combination of a Padé approximant and artificial intelligence (symbolic regression). Numerical experiments in Matlab using 2 million quasi-Monte Carlo samples indicate that the relative error of this new rational approximation does not exceed 0.866%. Moreover, these numerical experiments show that the novel rational approximation is approximately two times faster than the exact solution given by the Wright omega function.

Abstract

In this reply, we present updated approximations to the Colebrook equation for flow friction. The equations are equally computational simple, but with increased accuracy thanks to the optimization procedure, which was proposed by the discusser, Dr. Majid Niazkar. Our large-scale quasi-Monte Carlo verifications confirm that the here presented novel optimized numerical parameters further significantly increase accuracy of the estimated flow friction factor.

Abstract

The Colebrook equation is a popular model for estimating friction loss coefficients in water and gas pipes. The model is implicit in the unknown flow friction factor, f . To date, the captured flow friction factor, f , can be extracted from the logarithmic form analytically only in the term of the Lambert W-function. The purpose of this study is to find an accurate and computationally efficient solution based on the shifted Lambert W-function also known as the Wright ω-function. The Wright ω-function is more suitable because it overcomes the problem with the overflow error by switching the fast growing term, y=W(ex), of the Lambert W-function to series expansions that further can be easily evaluated in computers without causing overflow run-time errors. Although the Colebrook equation transformed through the Lambert W-function is identical to the original expression in terms of accuracy, a further evaluation of the Lambert W-function can be only approximate. Very accurate explicit approximations of the Colebrook equation that contain only one or two logarithms are shown. The final result is an accurate explicit approximation of the Colebrook equation with a relative error of no more than 0.0096%. The presented approximations are in a form suitable for everyday engineering use, and are both accurate and computationally efficient.

Abstract

This reply gives two corrections of typographical errors in respect to the commented article, and then provides few comments in respect to the discussion and one improved version of the approximation of the Colebrook equation for flow friction, based on the Wright ω-function. Finally, this reply gives an exact explicit version of the Colebrook equation expressed through the Wright ω-function, which does not introduce any additional errors in respect to the original equation. All mentioned approximations are computationally efficient and also very accurate. Results are verified using more than 2 million of Quasi Monte-Carlo samples.

Abstract

This paper provides a new unified formula for Newtonian fluids valid for all pipe flow regimes from laminar to fully rough turbulent flow.This includes laminar flow; the unstable sharp jump from laminar to turbulent flow; and all types of turbulent regimes, including the smooth turbulent regime, the partial non-fully developed turbulent regime, and the fully developed rough turbulent regime. The new unified formula follows the inflectional form of curves suggested in Nikuradse’s experiment rather than the monotonic shape proposed by Colebrook and White. The composition of the proposed unified formula uses switching functions and interchangeable formulas for the laminar, smooth turbulent, and fully rough turbulent flow regimes. Thus, the formulation presented below represents a coherent hydraulic model suitable for engineering use. This new flow friction model is more flexible than existing literature models and provides smooth and computationally cheap transitions between hydraulic regimes.

Abstract

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.

Abstract

This review paper gives Excel functions for highly precise Colebrook’s pipe flow friction approximations developed by users. All shown codes are implemented as User Defined Functions – UDFs written in Visual Basic for Applications – VBA, a common programming language for MS Excel spreadsheet solver. Accuracy of the friction factor computed using nine to date the most accurate explicit approximations is compared with the sufficiently accurate solution obtained through an iterative scheme which gives satisfying results after sufficient number of iterations. The codes are given for the presented approximations, for the used iterative scheme and for the Colebrook equation expressed through the Lambert W-function (including its cognate Wright ωfunction). The developed code for the principal branch of the Lambert W-function has additional and more general application for solving different problems from variety branches of engineering and physics. The approach from this review paper automates computational processes and speeds up manual tasks.