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
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 [...]
Even a relatively simple equation such as Colebrook’s offers a lot of possibilities to students to increase their computational skills. The Colebrook’s equation is implicit in the flow friction factor and, therefore, it needs to be solved iteratively or using explicit approximations, which need to be developed using different approaches. Various procedures can be used for iterative methods, such as single the fixed-point iterative method, Newton–Raphson, and other types of multi-point iterative methods, iterative methods in a combination with Padé polynomials, special functions such as Lambert W, artificial intelligence such as neural networks, etc. In addition, to develop explicit approximations or to improve their accuracy, regression analysis, genetic algorithms, and curve fitting techniques can be used too. In this learning numerical exercise, a few numerical examples will be shown along with the explanation of the estimated pedagogical impact for university students. Students can see what the difference is between the classical vs. floating-point algebra used in computers.
Abstract
Even a relatively simple equation such as Colebrook’s offers a lot of possibilities to students to increase their computational skills. The Colebrook’s equation is implicit in the flow friction factor and, therefore, it needs to be solved iteratively or using explicit [...]
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 R e and the known relative roughness of a pipe inner surface ε*; λ= ξ (R e, ε*, λ). 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 l o g-call in total for the whole procedure (expensive l o g-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
Abstract
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 R e and the known relative roughness of a pipe inner surface ε*; [...]