Abstract

In many situations ranging from geophysics to chemical engineering, turbulent drag moves particle clouds. I will present and compare various numerical approaches. On the one hand the mean velocity profile above ground is systematically constructed subtracting momentum loss; on the other hand the intrinsic spatio-temporal fluctuations are imposed from empirical distributions on point-like fluid particles. Various applications are explored. One is saltation, i.e. Aeolian transport of sand, discovering that the onset of particles flux exhibits a first order transition with hysteresis. The inclusion of mid-air grain collisions is found to increase the flux considerably due to the formation of a floating “soft bed” that screens energy-rich grains (saltons) from hitting the ground. Solving the fluid motion with the Lattice Boltzmann Method the effect of particle-particle collisions on preferential concentration is also investigated. Another application is powder mixing in a channel due to turbulent fluctuations. Following A.M. Reynolds (2003), a stochastic differential equation is solved for the motion of fluid particles that are attached to real particles. The dependence of the observed diffusive behaviour on Reynolds and Stokes number is monitored. Finally, also spatial correlations in the velocity field are imposed by a Heisenberg-type Hamiltonian. 

Abstract

This work is about the development of a general algorithm for the numerical solution of flow equations: the Navier-Stokes set. This set of differential equations models the time dependent behavior of fluids. It is formed by continuity, linear momentum and an energy transport equations. The algorithm here described is a general one since it can handle equally a great variety of problems, ranging from incompressible to compressible flows, viscous to inviscid, stationary and transient, all of them phisically modeled by the same set of differential equations.

In the present work, a quest for a general algorithm is described, following one of many possible ways to tackle the problem. In general, this is done extending methods either from compressible to incompressible flows or from incompressible to compressible ones.

Abstract

We present a formulation for analysis of turbulent incompressible flows using a stabilized finite element method (FEM) based on the finite calculus (FIC) procedure. The stabilization terms introduced by the FIC approach allow to solve a wide range of fluid flow problems at different Reynolds numbers, including turbulent flows, without the need of a turbulence model. Examples of application of the FIC/FEM formulation to the analysis of 2D and 3D incompressible flows at large Reynolds numbers exhibiting turbulence features are presented.

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

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.

Abstract

The Colebrook equation is implicitly given in respect to the unknown flow friction factor λ; λ=ζ(Re, ε∗, λ)which cannot be expressed explicitly in exact way without simplifications and use of approximate calculus. A common approach to solve it is through the Newton–Raphson iterative procedure or through the fixed-point iterative procedure. Both require in some cases, up to seven iterations. On the other hand, numerous more powerful iterative methods such as threeor two-point methods, etc. are available. The purpose is to choose optimal iterative method in order to solve the implicit Colebrook equation for flow friction accurately using the least possible number of iterations. The methods are thoroughly tested and those which require the least possible number of iterations to reach the accurate solution are identified. The most powerful three-point methods require, in the worst case, only two iterations to reach the final solution. The recommended representatives are Sharma–Guha–Gupta, Sharma–Sharma, Sharma–Arora, Džuni´c–Petkovi´c–Petkovi´c; Bi–Ren–Wu, Chun–Neta based on Kung–Traub, Neta, and the Jain method based on the Steffensen scheme. The recommended iterative methods can reach the final accurate solution with the least possible number of iterations. The approach is hybrid between the iterative procedure and one-step explicit approximations and can be used in engineering design for initial rough, but also for final fine calculations.

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 ε*; λ= ξ (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

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.