This article should examine ideas of how a chatbot can support decision making processes by col-lecting general opinion about the certain topics. The ChatGPT is tested in this article in respect of natural gas security of supply of the Republic of Serbia in the context of the current political situation in the world. The author first describes its own stance on the topic and then the ChatGPT is asked for its opinion (which most probably collects and compiles phrases from English speaking sources to build a text in a smart way). The topic is ort only about security of gas supply of Serbia but also political and highly sensitive, so it is examined if the ChatGPT can be useful tool to reach a balance between less on more favorable options in the process of policy making.

Here is shown method for the hydraulic solution of a looped gas-pipeline networks. Calculation of presented network is done according to principles of Hardy Cross method. The optimization was carried out by iteration of the pipes diameters, node consumptions are known and flow velocities through pipes have to stand below certain values. Accent is on determination of appropriate friction factor, and on selection of representative equation for natural gas flow under presented conditions in the network. Inappropriate usage of friction factor, equally as inappropriate usage of gas flow equation can lead to inaccurate final results. Here is shown new facts in comparison to previous calculation of gas distribution network in Kragujevac, Serbia which is done in 1994. After the implementation, measurements in situ have performed, and real measured values deviate from calculated. Causes for these errors are investigated, and improved and more accurate procedure is shown.

Paper presents the energy policy of the Republic of Serbia with special attention to the energy situation on the government controlled territory. South Serbian autonomous province of Kosovo and Metohija is under UN jurisdiction since the 1999 according to UNSC Resolution 1244. Renewable energy sources are rarely used in Serbia with exception of energy from hydropower plants, but in this sector priorities in geothermal and energy coming from biomass recently increased. In natural gas sector, Serbia has the deal with Russia for construction of South Stream gas-line through Serbia and for construction of the first underground storage in depleted gas reservoir in Banatski Dvor. In 2008, Serbia also sold 51% of the government founded petroleum industry – NIS which has exclusive monopoly for exploitation of crude oil. Serbian government has complete monopoly in electric power sector. Electric power infrastructure became technologically obsolete, and operative efficiency is at very low level. Serbia has not yet decided whether Serbian Electric Power Industry – EPS will be privatized. District heating sector mostly natural gas fuelled is highly inefficient and it is in jurisdiction of local municipalities but also has social component dictated by central government.

Accent is on determination of appropriate friction factor of the pipes and on selection of the representative equation for water or natural gas flow which is valuable for existing conditions in the looped network of pipelines. Note that in a municipal gas pipeline, natural gas can be treated as incompressible fluid (liquid) i.e. as water or oil. Even under this circumstance, calculation of water pipelines cannot be literary copied and applied for calculation of gas pipelines. Inappropriate friction factor, equally as e.g. inappropriate usage of water flow equations for calculation of gas networks can lead to inaccurate final results. Few iterative methods for determining the optimal hydraulic solution of water- and gas- looped pipeline networks, such as, Hardy Cross, modified Hardy Cross, node-loop method, node and M.M. Andrijashev method, will be shown. Speed of convergence will be compared and discussed using a simple network with three loops.

Today, two very efficient methods for calculation of flow distribution per branches of a looped gas pipeline are available. Most common is improved Hardy Cross method, while the second one is so-called unified node-loop method. For gas pipeline, gas flow rate through a pipe can be determined using Colebrook equation modified by AGA (American Gas Association) for calculation of friction factor accompanied with Darcy-Weisbach equation for pressure drop and second approach is using Renouard equation adopted for gas pipeline calculation. For the development of Renouard equation for gas pipelines some additional thermodynamic properties are involved in comparisons with Colebrook and Darcy-Weisbach model. These differences will be explained. Both equations, the Colebrook’s (accompanied with Darcy-Weisbach scheme) and Renouard’s will be used for calculation of flow through the pipes of one gas pipeline with eight closed loops which are formed by pipes. Consequently four different cases will be examined because the network is calculated using improved Hardy Cross method and unified node-loop method. Some remarks on optimization in this area of engineering also will be mentioned.

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.

Using only a limited number of computationally expensive functions, we show a way how to construct accurate and computationally efficient approximations of the Colebrook equation for flow friction. The presented approximations are based on the asymptotic series expansion of the Wright ω-function and symbolic regression. The results are verified with 8 million of Quasi-Monte Carlo points covering the domain of interest for engineers. In comparison with the built-in “wrightOmega” feature of Matlab R2016a, the herein introduced related approximations of the Wright ω-function significantly accelerate the computation. With only two logarithms and several basic arithmetic operations used, the presented approximations are not only computationally efficient but also extremely accurate. The maximal relative error of the most promising approximation which is given in the form suitable for engineers’ use is limited to 0.0012%, while for a little bit more complex variant is limited to 0.000024%.

Original and improved version of the Hardy Cross iterative method with related modifications are today widely used for calculation of fluid flow through conduits in loops-like distribution networks of pipes with known node fluid consumptions. Fluid in these networks is usually natural gas for distribution in the municipalities, water in waterworks or hot water in district heating system, air in the case of ventilation systems in buildings or mines, etc. Since, the resistances in these networks depend of flow, problem is not linear like in electrical circuits, and iterative procedure must be used. In both version of the Hardy Cross method, in original and in the improved one, initial results of calculation in iteration procedure is not flow, but rather the correction of flow. Unfortunately, these corrections should be added to or subtracted from a flow calculated in previous iteration according to complicate algebraic rules. After the here presented node-loop method, final results in each of the iterations is flow directly rather than flow correction. In that way complex algebraic scheme for sign of flow correction is avoided, while the final results still remain unchanged. Numbers of required iterations for the same results are comparable with the improved Hardy Cross method.

Hardy Cross originally proposed a method for analysis of flow in networks of conduits or conductors in 1936. His method was the first really useful engineering method in the field of pipe network calculation. Only electrical analogs of hydraulic networks were used before the Hardy Cross method. A problem with the flow resistance versus the electrical resistance makes these electrical analog methods obsolete. The method by Hardy Cross is taught extensively at faculties and it still remains an important tool for analysis of looped pipe systems. Engineers today mostly use a modified Hardy Cross method which threats the whole looped network of pipes simultaneously (use of these methods without computers is practically impossible). A method from the Russian practice published during 1930s, which is similar to the Hardy Cross method, is described, too. Some notes from the life of Hardy Cross are also shown. Finally, an improved version of the Hardy Cross method, which significantly reduces number of iterations, is presented and discussed. Also we tested multi-point iterative methods which can be used as substitution for the Newton-Raphson approach used by Hardy Cross, but this approach didn’t reduce number of required iterations to reach the final balanced solution. Although, many new models have been developed since the time of Hardy Cross, main purpose of this paper is to illustrate the very beginning of modeling of gas and water pipe networks or ventilation systems.

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 . To date, the captured flow friction factor can be extracted from the logarithmic form analytically only in the term of the Lambert -function. The purpose of this study is to find an accurate and computationally efficient solution based on the shifted Lambert -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 of the Lambert -function to the series expansions that further can be easily evaluated in computers without causing overflow run-time errors. Although the Colebrook equation transformed through the Lambert -function is identical to the original expression in term of accuracy, a further evaluation of the Lambert -function can be only approximate. Very accurate explicit approximations of the Colebrook equation that contains only one or two logarithms are shown. The final result is an accurate explicit approximation of the Colebrook equation with the relative error of no more than 0.0096%. The presented approximations are in the form suitable for everyday engineering use, they are both accurate and computationally efficient.

This paper gives a new unified formula for the Newtonian fluids valid for all pipe flow regimes from laminar to the fully rough turbulent. It includes laminar, unstable sharp jump from laminar to turbulent, and all types of the turbulent regimes: smooth turbulent regime, partial non-fully developed turbulent and fully developed rough turbulent regime. The formula follows the inflectional form of curves as suggested in Nikuradse’s experiment rather than monotonic shape proposed by Colebrook and White. The composition of the proposed unified formula consists of switching functions and of the interchangeable formulas for laminar, smooth turbulent and fully rough turbulent flow. The proposed switching functions provide a smooth and a computationally cheap transition among hydraulic regimes. Thus, the here presented formulation represents a coherent hydraulic model suitable for engineering use. The model is compared to existing literature models, and shows smooth and computationally cheap transitions among hydraulic regimes.

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, ε/D with the output unknown 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 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 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 λ, 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é approximation, practically the same error is obtained also using only one logarithm.

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. Common approach to solve it is through the Newton-Raphson iterative procedure or through the fixed-point iterative procedure. Both requires in some case even eight iterations. On the other hand numerous more powerful iterative methods such as three-or 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 worst case only two iterations to reach final solution. The recommended representatives are Sharma-Guha-Gupta, Sharma-Sharma, Sharma-Arora, Džunić-Petković-Petković; Bi-Ren-Wu, Chun-Neta based on Kung-Traub, Neta, and Jain method based on Steffensen scheme. The recommended iterative methods can reach the final accurate solution with the least possible number of iterations. The approach is hybrid between iterative procedure and one-step explicit approximations and can be used in engineering design for initial rough, but also for final fine calculations.

Empirical Colebrook equation from 1939 is still accepted as an informal standard to calculate friction factor during the turbulent flow (4000 < Re < 10⁸) through pipes from smooth with almost negligible relative roughness (ε/D→0) to the very rough (up to ε/D = 0.05) inner surface. The Colebrook equation contains flow friction factor λ in implicit logarithmic form where it is, aside of itself; λ, 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 many available explicit approximations to the Colebrook equation, λ ≈ f(Re, ε/D), it is necessary to determinate value of the friction factor λ 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’s approach (3^rd order, 2^nd order: Halley’s and Schröder’s method and 1^st 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’s approach to the Colebrook’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.

The eighty years 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 and the known relative roughness of a pipe inner surface ; . It is based on logarithmic law in the form that captures the unknown flow friction factor in a way from which 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 -call in total for the whole procedure (expensive -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.

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.