## ABSTRACT

In improving the energy efficiency of water transport systems, two critical stages are involved: assessment (to understand the system's operation and identify potential energy savings) and auditing (to locate and break down the energy losses). Both stages are based on energy balances, which can be conducted using either the extended Bernoulli equation or the energy integral equation. Both equations can be applied, but depending on the system, data availability, and the kind of study to be performed, one is preferable over the other. This paper analyses, applies and compares both equations, with a particular focus on the less commonly used energy integral equation in the hydraulic field. This more general equation includes thermal and transient effects and it is more suitable for analyzing complex systems. In contrast, the extended Bernoulli equation, while simpler to apply, can lead to the loss of relevant information, such as the evaluation of the topographic energy. The main objective of this work is to bridge the gap between these two fundamental energy equations and recommend the most appropriate one for hydraulic problems. Real examples are presented to show their differences and validate our recommendations.

## HIGHLIGHTS

A comparison is made between the energy integral and Bernoulli's equation, in order to know which one is more appropriate when analyzing a system energetically.

The most suitable energy analysis for multi-topology systems, which are those with different operating layouts, is explained in depth.

A guide to analyze the energy of pressurized water systems is presented, with the ultimate goal to improve their efficiency.

## NOTATIONS

*E*_{e}electric energy consumed

*E*_{hp}hydraulic energy delivered by pumps (

*E*_{hp}*=**E*._{e}*η*_{e}.*η*)_{p}*E*_{N}natural energy supplied to the system

*E*_{T}total energy supplied to the system (=

*E*_{e}*+**E*)_{N}*E*_{ou}output energy useful (supplied to users)

*E*_{ol}output energy lost (in leaks)

*E*_{o}output energy (=

*E*_{ou}*+**E*)_{ol}*E*_{oum}minimum energy required by users

*E*_{ΔE}excess (Δ) of energy supplied to users (=

*E*_{ou}*−**E*)_{oum}*H*_{oi,m}minimum piezometric required

*h*_{p}pumping head

*I*_{e}energy intensity

real energy intensity

*P*_{t}power supplied by compensation tanks

*P*_{N}power gravitational (natural)

*P*_{Sl}total power system losses

*P*_{pl}power pumping station losses

*P*_{T}total power supplied

*P*_{eh}hydraulic power input (entrance) to the system

*P*_{e}total electric power input

*P*_{f}power dissipated by friction

*P*_{oh}output hydraulic power

*P*_{ohl}power lost in leaks

*P*_{ohu}useful delivered power

*P*_{ohum}useful (minimum required) delivered power

*P*_{ml}power lost (other miscellaneous losses)

*P*_{ΔP}power in excess (Δ) supplied

*p*_{m}minimum pressure

*Q*_{t}tank net flow

*u*internal energy

*z*elevation

*W*_{s}hydraulic energy delivered by pumps (

*W*_{s}*=**E*)_{hp}friction energy variation

- Δ
*E*_{G} gravitational energy variation

thermal energy variation

- Δ
*T* temperature variation

electric motor efficiency

pump efficiency

global efficiency

global efficiency considering the minimum standard pressure

## INTRODUCTION

The first question to address is why the extended Bernoulli equation is preferentially used in hydraulic flows whereas the energy integral equation is, surprisingly, rarely applied. The primary reason is that since most hydraulic flows are assumed to be incompressible and one-dimensional, the extended Bernoulli equation fits these problems perfectly. However, the energy integral equation (White 1973), which includes thermal and mechanical terms as well as various inputs and outputs, is more general and therefore not as immediately applicable. In fact, it can be used to analyze any problem because it performs a general energy balance over a control volume (CV) limited by a control surface (CS). Furthermore, as these two equations have historically different origins, they can only be compared after formulating restrictive hypotheses to simplify the energy integral equation to equate it to the extended Bernoulli equation, a generalization of the original equation formulated by Bernoulli in the 18th century (Rouse & Ince 1963).

*u*), the kinetic energy (

*v*

^{2}/2, where

*v*is the velocity) and the gravitational energy (

*gz,*where

*g*is the acceleration due to gravity and z the elevation). The most general form of the energy integral equation (

*three*-

*dimensional*) is:

This form was derived by separating the flow work through the CS (between the system and the external medium) and the shaft work, commonly associated with turbomachines. According to the sign convention for work (work done by the system is considered positive), pump work is negative and turbine work is positive. The other left-hand side term of the equation, the heat transfer rate, , represents the thermal energy exchange between the external medium and the system through the CS because of the existence of a temperature gradient (heat transfer to the system is considered positive). Other variables are, *ρ*, the fluid density; d*V*, the differential of volume inside the CV; *p*, the pressure and, the differential flow through the CS. This formulation has few constraints and is widely used in thermal fluid mechanics.

*γ*is the specific weight of the fluid. This expression is conceptually important. It shows the different forms that a fluid has to convert mechanical energy but has limited applicability. Therefore, it needs to be generalized. As a starting point, the Bernoulli equation is applied to a pipe, assuming its axis as the average streamline of the flow, and including the friction between the fluid and the walls of the pipe throughout the Darcy–Weisbach equation. The addition of shaft work, provided by a pump or removed by a turbine, (White 1973) yields the extended Bernoulli equation:where

*f*is the friction factor (dimensionless),

*L*the pipe length (m),

*h*the pumping head (m) and

_{p}*h*the turbine head (m). A more complete formulation, including inertial terms, is used to model rigid transients, but in the following energy analyses, it is not of interest because of the very different scale of time of these problems. Henceforth, we simply refer to Equation (3) as the Bernoulli equation.

_{t}Since the Bernoulli equation (Equation (3)) involves only mechanical energy, whereas the energy integral equation (Equation (1)) includes both mechanical and thermal energy, a crucial point in comparing the equations is related to the coupled nature of the physical problem, that is, with the possible conversion between thermal and mechanical energy, which is the clear boundary between thermal fluid mechanics and hydraulics.

In hydraulics, the special case of incompressible fluids, with both density and viscosity constants, is of considerable importance (Batchelor 1967; Panton 2013). The continuity and the momentum equations are simpler and, more importantly, decoupled from the energy equation (the momentum equation does not include the temperature). Consequently, the mass and momentum equations form a set of four scalar equations with four mechanical unknowns ( and *p*). Later, when necessary, the obtained velocity field can be substituted into the energy equation to determine the temperature field. On the other hand, for an incompressible fluid, thermal and mechanical energies are only tied by friction.

The integral energy equation (Equation (1)) can be applied to both incompressible and compressible flows. One example is the analysis of a water hammer (elastic conduit with *slightly compressible fluid*) from an energy perspective (Karney 1990), where velocity and pressure changes are related to account for water compressibility (elastic effects), while the density is assumed only a function of the pressure (which is captured via the bulk modulus of elasticity). Thus, the flow is considered essentially isothermal (due to temperature, there is a negligible change in flow properties) and the thermal and mechanical variables are decoupled.

*et al.*1990). They are a fixed CV with one inlet and one outlet section (1 and 2, respectively), being the flow normal to control surfaces and uniform properties over plain perpendicular cross-sections. Briefly, the flow must be one-dimensional, a strong constraint that makes the difference (Figure 1).

*q*is due to the conversion of mechanical energy into thermal energy that happens because of the viscosity of the fluid. Therefore, represents the mechanical energy loss per unit of weight, through a mass friction loss term. So, in an ideal fluid flow (no viscosity and, therefore, no friction losses), , and therefore, . In conclusion, the simplest Bernoulli equation (Equation (2)) applies, while both mechanical and thermal energies are conserved separately.

The above analysis confirms that, for an incompressible fluid flow, friction is the only link between thermal and mechanical energy. From a mechanical perspective, friction takes its toll converting useful mechanical energy into dissipated heat energy. If there is no heat transfer, the temperature of the fluid increases because of the friction, and, in the end, *u*_{2} becomes greater than *u*_{1}. On the other hand, when the flow is not adiabatic (the system is not well insulated), due to the temperature gradient, heat will flow to the surroundings. Last, for an isothermal flow (both temperature and *u,* constants), the system has a loss of heat, *q,* at the same rate that friction converts mechanical into thermal energy.

In short, the same problem can be analyzed from different perspectives using both equations. With different origins (one being integral and thermal, the other differential and mechanical) they follow opposite paths: generalization, in the case of the Bernoulli equation, and simplification, for the integral energy equation, where an initial global power balance integrated over time becomes an energy balance. This last versatile equation can be employed for energy diagnosis and auditing but requires more information than the simplest Bernoulli equation. The latter is recommended to diagnose and audit simple systems and, when applied repeatedly, it can solve more complex systems (such as hydraulic networks), albeit in snapshots.

To show the particularities of each equation, this paper is organized as follows. The next section is a short reminder, through the energy intensity concept, of how to apply the Bernoulli equation to simple and complex systems. The rest of the study is devoted to the integral energy equation. After some previous remarks, two real examples are discussed. First, a closed and hydraulic system, that shows the complementarity of both equations. Later, a multi-topology system, with multiple inputs and outputs, proved the versatility of the global energy equation. A final summary collects all these recommendations.

## ASSESSING AND AUDITING WATER TRANSPORT SYSTEMS WITH THE BERNOULLI EQUATION

*I*(kWh/m

_{e}^{3}) that relates water and energy. Since its dimensions (energy per unit of volume) are pressure dimensions as well, for an incompressible fluid, the energy intensity can be expressed in terms of head (

*H*) in m, as

*p*=

*γ H*. This relationship (with

*γ*= 9,810 N/m

^{3}, International System units) is a simple unit change conversion, resulting in:

For simple systems (one input and one output), these concepts have been well developed (Del Teso *et al.* 2023) and can be summarized as follows. On one side, the real energy intensity, *I _{er}*, is known from field data, including the energy consumed (billed kWh) and the measured pumped water volume (m

^{3}) at the delivering point. On the other side, the energy intensity objective,

*I*, is determined by setting target efficiency values. The difference

_{eo}*I*–

_{er}*I*indicates the margin of improvement.

_{eo}This estimation has a weak point. Bernoulli is a stationary equation and during the period considered, the system's conditions can change (e.g. pumping groundwater systems, with water table level variations along the period of analysis). Therefore, the difference *I _{er} – I_{eo}* provides a snapshot assessment, that is strictly valid for a given instant. To some extent, this drawback can be mitigated by using average values over the whole considered period.

Bernoulli's equation can also be used to assess complex systems (many output nodes and, sometimes, more than one input node). A crucial aspect is to identify the final node of the equation because the first one is always the source. This final node, named the critical node, is the one among all the systems that need more energy to achieve the required pressure. Applying the equation between the source and the critical node provides the unitary energy that the system needs, and this maximum value is extrapolated to the whole system. But again, additional care must be taken with the Bernoulli equation, assessing systems with variable time conditions along the considered period of time. In fact, the critical node can depend on the scenario and, even being the same node, the energetic requirements to be covered between both nodes can be different. Again, the energy used to assess the system is the average energy required over the period.

Another fact to consider is that, in the rest of the output nodes, the pressure, although higher than required, is not known. This drawback can be addressed by applying successively the Bernoulli equation between the source and all these output nodes, as many times as necessary up to cover the whole network. However, this is most of the time an awkward strategy due to the high number of nodes. Therefore, a complete physical picture of how the system works can be calculated easily with Bernoulli's equation, although at the expense of simplicity (Cabrera *et al.* 2021).

In conclusion, Bernoulli is the appropriate energy equation to assess and audit simple systems. Being one-dimensional, both stages (assessment and audit) are, in practice, covered simultaneously. With these results, the energetic system's behavior can be easily labeled (Cabrera *et al.* 2023; del Teso *et al.* 2023).

## THE ENERGY INTEGRAL EQUATION. PREVIOUS REMARKS

*ρ*= C) and the flow is assumed to be adiabatic (no thermal exchange between the system and the exterior,

*dQ/dt*= 0). In these conditions, Equation (1) becomes:

Thus, the energy integral equation differs from the Bernoulli equation in dimensionality (a CV versus a stream tube, Figures 1 and 2) and the capability to analyze transient flows. This general equation has been previously applied to diagnose (Cabrera *et al.* 2015) and audit (Cabrera *et al.* 2010) complex systems. In this study, it is extended to multi-topology systems.

The terms of Equation (7) are physically meaningful. Firstly, because the shaft work added (by pumps) or removed (by turbines) from the system balances the change with the time of the energy stored inside the CV, first integral (it is non-zero if there are tanks inside the CV, with level's variation), and second integral which is a flow term that represents a power global balance between the outgoing and incoming flows. Both integrals include the kinetic energy (negligible in our analysis), the internal energy (its change represents the power dissipated by friction, Equation (5)) and the gravitational head. The pressure head is the only difference between both integral terms.

Some remarks concerning this equation follow:

### Closed and open systems

Systems can be open or closed (with or without flow throughout the CS, respectively). When dealing with pressurized water transport, most of the systems to analyze are open. In closed systems, because there are no input nor output flows of energy, it does not make sense to define efficiency.

### Thermal implications of water transport systems

There is a growing interest in modeling the water temperature evolution along urban water networks (Agudelo-Vera *et al.* 2020) because its value impacts the water quality. In fact, some countries have set a maximum temperature threshold of 25 °C at the customer's tap (Blokker & Pieterse-Quirijns 2013). In the context of climate change, with increasingly hot summers, this issue becomes more challenging. Obviously, for other uses (as irrigation), this question has no interest.

In any case, thermal and dynamic problems are decoupled because water is an incompressible fluid. Although the main focus of EPANET (Rossman 2000) is on water quality issues, its solver ignores any potential thermal exchange between soil and water. Even its recent update (Rossman *et al.* 2020) has ignored this question. Conversely, the models proposed to determine the water temperature evolution (Burch & Christensen 2007; Blokker & Pieterse-Quirijns 2013) are, as expected, decoupled of the hydraulic equations.

In summary, the adiabatic hypothesis assumed in this work may not apply in some cases, but this has no influence on the mechanical energy analysis being performed. The hydraulic variables calculated, and therefore the temperature variation Δ*T* due to friction, will be the same no matter whether the heat transfer soil–water effect is included or not. The only changed value should be the final water temperature, important from the water quality side, but irrelevant for the purpose of this analysis.

### Gravity terms

Gravitational energy, which depends on the potential elevation change of the considered mass, has no absolute value. Therefore, it is a crucial issue to set correctly the system's reference level. This is the lowest active node (consumption or source). Its elevation is set to zero. If real elevations (over the sea level) are kept, the gravitational energy is oversized although the final energy balance will be satisfied because inlet and outlet flows are equal (*ρ**=**C*) and the gravitational terms linear. Both inlet and outlet gravitational terms are equally incremented. However, a tank inside the CV breaks the energy balance, because the changes in the energy stored in the tank depend on the square elevation of the level while inlet and outlet flows are no longer equal. Therefore, node elevations must all refer to the lowest node, an invariable reference when it is a consumption node. Nevertheless, if the lowest node is a tank, its elevation varies with the surface level but can be kept as a reference. This will be further analyzed in the following section.

## THE ENERGY INTEGRAL EQUATION APPLIED TO A CLOSED SYSTEM

*A*

_{1}and

*A*

_{2}) connected by a pipe (section

*A*). The flow begins when the valve is instantaneously opened. The level of reference for the elevation (gravitational energy) is the initial water level of the downstream tank. No heat transfer is considered in the analysis.

_{p}*et al.*1999) if:

*C*is the integration constant.

For the considered period, Equation (10) shows that in the closed system the loss of gravitational energy, Δ*E _{G}*, equalizes the positive internal energy variation (Δ

*E*

_{T}*=*Δ

*E*). Assuming a quasi-static flow (Rao & Bree 1977), without changes of energy inside the pipe in the period

_{f}*t,*water thermal changes came from the absorption of heat generated by friction between pipe and fluid.

*E*< 0) is equal to the gained internal energy. That gravitational energy loss (Figure 3) is:

_{G}The negative value confirms the loss of gravitational energy. The first term is the energy lost at the first tank, equal to the weight of the discharged volume multiplied by the change of elevation of the center of gravity (assuming a constant reference level). The second term is the energy gained by the lower tank during the filling with respect to the initial reference. Therefore, this term is the correction of the error introduced because of the variation of the reference level.

In short, for a variable reference elevation, the compensation term of Equation (16) restores the balance and therefore can be kept as a reference level. If it is a reservoir, the lowest level is constant and the correction is null (). In a reverse flow situation (pumping line), the lower tank is emptied while the upper one is filled. In that case, Equation (14) is the same, but the signs of the terms are interchanged.

*c*the specific heat (1 cal/°C gr), and ΔT the temperature variation. Therefore, in this part of the analysis, the only mechanical variable that influences ΔT is the volume moved, regardless of the time required, which depends on the system's characteristics.

_{e}Data of the system depicted in Figure 3 are: *D*_{1} = 40 m; *D*_{2} = 20 m;_{1} = 200 m;_{2} = 100 m; *ξ*_{10} = 6 m; *ξ*_{20} = 1 m; *D _{p}* = 100 mm and

*L*= 5,000 m. The volume transferred between both tanks is 1,256.64 m

^{3}. Two scenarios are considered: (a) valve completely open (no local losses;

*k*≈ 0) and (b) valve almost closed (local losses are relevant;

*k*= 500).

^{3}), the level in Tank 1 has decreased 1 m and the level in Tank 2 has increased by 4 m. Therefore, in the end, both tanks have the same water level (

*ξ*

_{1f}=

*ξ*

_{2f}= 5 m). Assuming a constant head loss equal to the average between the initial difference in tank levels (105 m) and the final one (100 m), the energy lost by friction (head loss times the discharged volume) results:

This result is equal to the loss of gravitational energy, Δ*E _{G}*, 350.99 kWh. This value has been obtained from Equation (14) with

*h*

_{o}= 105 m and

*h*= 104 m. From this value, the water temperature variation due to friction is Δ

_{f}*T*= 0.24 °C, Equation (17), a minuscule value, because of the high specific heat of water. In fact, the energy intensity to heat water by 1 °C is 1.16 kWh/m

^{3}(Equation (5)), which is equivalent to elevating 1 m

^{3}of water by 426 m. Thus, it is reasonable to ignore thermal effects due to hydraulic friction.

If the valve is open (*K* = 0), only the pipe dissipates energy but, if partially closed, the energy is lost at both elements (valve and pipe), although the integral energy equation is unable to break down their respective contribution. For this purpose, assuming a quasi-stationary flow (a valid assumption because ), it is necessary to model hydraulically, with EPANET, both elements. Results are:

(a) Valve open: Time required to transfer 1,256.64 m

^{3}, 32 h 40 min; initial velocity, 1.38 m/s; final velocity, 1.34 m/s; pipe average losses 102.5 m.(b) Valve partially closed: Time required to transfer 1,256.64 m

^{3}, 40 h 00 m; initial velocity, 1.13 m/s; final velocity, 1.10 m/s; pipe losses, at*t*= 0 h, 72.35 m and at*t*= 40 h 00 m, 68.90 m; average pipe friction losses, 70.63 m. Valve losses, at*t*= 0 h, 32.67 m and at*t*= 40 h 00 m, 31.06 m, average valve's losses, 31.87 m. Therefore, the valve is responsible for 31.09% of the total losses, because it is practically closed (*k*= 500).

From this simple but conceptually strong case study, it can be concluded that the initial level of the lowest tank, even being variable, can be set as the reference elevation node. Furthermore, it evidences that to complete the audit, it is necessary to solve the hydraulic system's behavior, showing that the discussed equations are also complementary to each other.

## THE ENERGY INTEGRAL EQUATION APPLIED TO MULTI-TOPOLOGY SYSTEMS

*n*and

_{e}, n_{o}, n_{p}*n*denote the number of entrances, outlets, pumping stations and tanks within the CV, respectively, while subscripts used are (

_{t}*e)*, entrance; (

*o)*, outlet; (

*p)*, pump and (

*t)*, tank. Last, the variables are:

*u,*internal energy;

*H,*piezometric hydraulic head;

*h,*pumping head;

*z,*elevation referred to the lowest active node;

*Q,*pipe flow and

*Q*tank net flow.

_{t},*P*, is the sum of the energies provided by sources (gravitational or natural power,

_{eh}*P*) and pumps minus (or plus) the tank's contribution, which is equal to the hydraulic output power,

_{N}*P*, plus the power dissipated by friction,

_{oh}*P*. Finally, the hydraulic output power is the sum of the useful delivered power (

_{f}*P*) plus the power lost through leaks,

_{ohu}*P*:

_{ohl}*L*, the hydraulic useful output power is:

_{j}A final concept must be considered: the minimum energy to be supplied (*H _{om,i}*), defined through the minimum pressure, because if the pressure exceeds the standard (minimum) value,

*p*, then

_{m}*H*

_{o,i}*>*

*H*and the system will deliver more energy than needed. That is, to some extent, an inefficiency that should be avoided.

_{om,i}### Assessment of multi-topology systems

The assessment of a specific topology within a multi-topology system without energy storage (no tanks inside the CV) can be conducted by comparing output with input powers, disregarding what occurs inside the CV. Within this approach for the diagnosis phase, the system can be assimilated to a black box, defined as ‘*a fictitious system representing a set of concrete systems onto which stimuli S impinge and out of which reactions R emerge. The constitution and structure of the box are altogether irrelevant to the approach under consideration, which is purely external or phenomenological*’ (Bunge 1963).

If the CV includes internal tanks, the flow *Q _{t}*, calculated from a differential quasi-static continuity relationship (Equation (20)), has to be known in order to perform the assessment. The additional inputs/outputs of these tanks (as water levels, z), although internal to the CV, are similar to the external inputs/outputs of the system.

Therefore, at the diagnostic stage, the system can be associated with a black box, with or without energy storage. Nevertheless, at the next stage, the audit, the whole inefficiency, , must be decoupled in (Equation (25), second term). To do that, each system's element must be hydraulically characterized to assess individually its work, complementing the integral energy equation. So, each possible topology must be diagnosed and audited. Nevertheless, it is convenient to carry out the analysis of each scenario to identify the most inefficient ones and to improve them.

#### Assessment of each of the topologies

*p*. It is calculated by replacing the hydraulic head outputs,

_{m}*H*by

_{o,i}*H*in Equation (24). The difference is the excess of delivered power,

_{om,i}*P*:

_{ΔP}#### Assessment of the whole system

Equations (27) and (29) provide, in power terms, the instantaneous efficiency (snapshots) for different topologies *L _{j}*. The whole system assessment over a period

*t*(month or a year) to calculate the average efficiency,

*η*, can be obtained by integrating Equation (27) (or Equation (29)) over time along all the scenarios that the system evolves through. However, it is easier and more direct to calculate

*η*through an energy balance, from:

Average piezometric heads (at the inlets/outlets of the system along t), based on the recorded pressures along

*t,*and the volumes measured at the inlets/outlets of the system for the same period t.The electric energy supplied to the system,

*E*, the sum of the electricity bills during period t._{ei}Over a long period of time, the contribution of the internal tanks is negligible (therefore, the third addend of the denominator, is zero).

### Audit of multi-topology systems

Once the system's energy losses have been determined with the assessment, it proceeds, when necessary (a poor performance result) to audit the system. This involves breaking down the inefficiencies. The audit will quantify all the power lost through friction,, power embedded in leaks,, excess of delivered power due to overpressure, and power inefficiency at the pumping station, . This breakdown is essential to identify what must be improved in the system.

#### Audit of each topology

*q*, and on the pipe's physical characteristics (diameters

_{k}*D*and lengths

_{k}*L*). Let

_{k}*n*

_{l}denote the total number of lines (with

*n*active for topology

_{l,j}*L*) and Δ

_{j}*h*the head loss of each active line. The power lost due to friction,

_{k}*P*, is

_{f,j}*n*, on the SC, internal nodes inside the CV,

_{oj}*n*, must be considered because, although there is no demand in these multi-topology systems, they can include leaks. Therefore,

_{i}*n*is the active internal nodes of the topology

_{i,j}*L*(with

_{j}*q*leak at node

_{ol,i}*i*, modeled as a pressure-demand), the power results:

This balance shows that the total power supplied, *P _{T,,j}*, is equal to the hydraulic input power , plus the pumping station power losses, , and equalizes the sum of the minimum useful output power, , the excess of the delivered power, , and the powers required to overcome internal inefficiencies (friction, and leaks, ), plus losses at the pumping station, . After the audit, the terms of Equation (39) are known.

#### Audit of the whole system

Formulating a global audit in terms of power in multi-topology systems is neither useful nor physically meaningful, but it is in terms of global energy as proposed in the assessment. It is better to analyze each topology with its inefficiencies, examine their causes, and amend them. The feasible solution is to improve, one by one, the efficiency of each system's topology, paying more attention to the most energy-consuming one.

## OPEN MULTI-TOPOLOGY SYSTEM: A REAL CASE STUDY

*n*), tanks (

_{e}and n_{o}*n*) and pumps (

_{t}*n*). The topology 2.2 (Figure 6) is not included in the analysis. as it is very simple. It consists of a decoupled simple system filling Tank 1, with no pumps working.

_{p}. | Topology 1 . | Topology 2 . | Topology 3 . | Topology 4 . | Topology 5 . | Topology 6 . |
---|---|---|---|---|---|---|

Pump 1 | ON | OFF | OFF | OFF | OFF | ON |

Pump 2 | ON | ON | ON | OFF | OFF | OFF |

Pump 3 | ON | OFF | ON | OFF | ON | ON |

. | Topology 1 . | Topology 2 . | Topology 3 . | Topology 4 . | Topology 5 . | Topology 6 . |
---|---|---|---|---|---|---|

Pump 1 | ON | OFF | OFF | OFF | OFF | ON |

Pump 2 | ON | ON | ON | OFF | OFF | OFF |

Pump 3 | ON | OFF | ON | OFF | ON | ON |

There are different energy sources (pumps or tanks) and different critical nodes for each topology. To assess each system with the Bernoulli equation is impractical because the initial and final nodes change permanently. However, the use of the energy integral equation is straightforward. Table 2 details the results of the diagnostic calculated with Equation (28). These results have been obtained from the inlet/outlet shown in Figure 6.

. | Topology 1 . | Topology 6 . | Topology 3 . | Topology 5 . | Topology 4 . | Topology 2.1 . |
---|---|---|---|---|---|---|

P (from wells) [kW] _{N} | 1.49 | 0 | 0 | 0 | 0 | 0 |

P (electric) [kW] _{Te} | 282.53 | 201.64 | 154.66 | 68.58 | 0 | 0 |

P (tanks) [kW] _{t} | 154.57 | 115.35 | 60.19 | 33.22 | −8.53 | −8.42 |

P = _{T}P_{N}+P_{Te}−P [kW] _{t} | 129.44 | 86.29 | 94.47 | 35.36 | 8.53 | 8.42 |

P [kW] _{ohu} | 28.99 | 52.27 | 20.14 | 10.43 | 7.02 | 7.77 |

P [kW] _{SI} | 100.45 | 34.02 | 74.34 | 24.93 | 1.51 | 0.65 |

Efficiency (η_{j}=P) _{ohu}/P_{T} | 0.22 | 0.61 | 0.22 | 0.30 | 0.82 | 0.92 |

. | Topology 1 . | Topology 6 . | Topology 3 . | Topology 5 . | Topology 4 . | Topology 2.1 . |
---|---|---|---|---|---|---|

P (from wells) [kW] _{N} | 1.49 | 0 | 0 | 0 | 0 | 0 |

P (electric) [kW] _{Te} | 282.53 | 201.64 | 154.66 | 68.58 | 0 | 0 |

P (tanks) [kW] _{t} | 154.57 | 115.35 | 60.19 | 33.22 | −8.53 | −8.42 |

P = _{T}P_{N}+P_{Te}−P [kW] _{t} | 129.44 | 86.29 | 94.47 | 35.36 | 8.53 | 8.42 |

P [kW] _{ohu} | 28.99 | 52.27 | 20.14 | 10.43 | 7.02 | 7.77 |

P [kW] _{SI} | 100.45 | 34.02 | 74.34 | 24.93 | 1.51 | 0.65 |

Efficiency (η_{j}=P) _{ohu}/P_{T} | 0.22 | 0.61 | 0.22 | 0.30 | 0.82 | 0.92 |

The water comes from the wells, which are the lowest nodes of the system. Therefore, only when both good pumps are ON (topology 1), is there a gravitational (or natural) input power contribution from well 1 (well 2 is the lowest). The total input power is the sum of the pump's power minus the power invested to fill the tanks, which only supply power to the system when they are emptying. This occurs in topologies 4 and 5. Summing up, for all topologies (except for L4 and L2.1), the pumping power exceeds the total input power because the internal tanks are being filled.

Efficiency is the ratio between the power output and power input (*P _{ohu}/P_{T}*) while their difference (

*P*) gives the power lost,

_{T}– P_{ohu}*P*, for each topology. The two most inefficient ones are L1 and L3 when Pumps 2 and 3 are simultaneously working.

_{Tl}The mathematical model is used to perform the audit (Equations (33) to (39)), calculating in which parts of the system the power is lost. Table 3 shows the results for topology (T1).

Breakdown Power Lost (P)
. _{Tl} | [kW] . |
---|---|

Losses at the pumping station (P) _{pl} | 74.57 |

Friction losses (P) _{f} | 13.58 |

Power lost by leaks (P) _{ohl} | 2.76 |

Power lost by filling the tanks from top (P) _{ml} | 9.52 |

Total Power Lost ( P)_{Tl} | 100.43 |

Breakdown Power Lost (P)
. _{Tl} | [kW] . |
---|---|

Losses at the pumping station (P) _{pl} | 74.57 |

Friction losses (P) _{f} | 13.58 |

Power lost by leaks (P) _{ohl} | 2.76 |

Power lost by filling the tanks from top (P) _{ml} | 9.52 |

Total Power Lost ( P)_{Tl} | 100.43 |

Similar tables for the other topologies will show their respective weak points. In this topology 1, major inefficiencies are located at the pumping stations.

A possible solution to enhance efficiency could be to remove tank 1 and pump 3. This would entail a direct pumping from well 2 to tank 2. Consequently, *P _{ml}* would disappear and less head would be needed in pump 2. Additionally, feeding the tanks from the bottom should be considered as an improvement, as it requires less head.

A final global assessment (in energy terms) provides relevant information. This value can be obtained from the different topologies' efficiencies, with an appropriated weighted average. However, it is often simpler to determine this value from annual inputs and outputs data. Table 4 presents volumes (elevated and delivered to users), annual energy consumption by the pumps, average pressure at the system outlets, elevation nodes, and minimum standard pressure.

. | Input . | Output . | ||
---|---|---|---|---|

Well 1 . | Well 2 + Pump 3 . | Industrial . | Town . | |

Volume [m^{3}] | 641,879 | 524,774 | 9,853 | 1,100,473 |

Electric energy consumed (billed) [kWh] | 403,085 | 338,885 | – | – |

Average output pressure [m] | – | – | 50 | 38 |

Minimum output pressure [m] | – | – | 30 | 30 |

Electric energy consumed, E _{e} | 741,970 | |||

Natural energy supplied, E _{N} | 3,498 | |||

Total energy injected, E _{T} | 745,468 | |||

Energy supplied to users, E _{ou} | 357,401 | |||

Minimum energy required by users, E _{oum} | 332,874 | |||

Average efficiency, | 0.48 | |||

Average efficiency (minimum pressure), | 0.45 |

. | Input . | Output . | ||
---|---|---|---|---|

Well 1 . | Well 2 + Pump 3 . | Industrial . | Town . | |

Volume [m^{3}] | 641,879 | 524,774 | 9,853 | 1,100,473 |

Electric energy consumed (billed) [kWh] | 403,085 | 338,885 | – | – |

Average output pressure [m] | – | – | 50 | 38 |

Minimum output pressure [m] | – | – | 30 | 30 |

Electric energy consumed, E _{e} | 741,970 | |||

Natural energy supplied, E _{N} | 3,498 | |||

Total energy injected, E _{T} | 745,468 | |||

Energy supplied to users, E _{ou} | 357,401 | |||

Minimum energy required by users, E _{oum} | 332,874 | |||

Average efficiency, | 0.48 | |||

Average efficiency (minimum pressure), | 0.45 |

As expected, the average efficiency calculated over a one-year period (0.45), is an intermediate value between the maximum efficiency (0.92 in T2.1) and the minimum value (0.22 in T1 and T3).

## DISCUSSION

Pressurized water transport systems can be divided into simple, complex (current networks) and multi-topology systems. All of them must be assessed in terms of their energy efficiency and, when it proceeds, audited. Previous papers have presented assessments and audits for simple (Cabrera *et al.* 2021; Cabrera *et al.* 2023) and complex systems (Cabrera *et al.* 2010; Cabrera *et al.* 2015) using either the Bernoulli equation or the integral energy equation.

The main objective of this study is to discuss which equation is most convenient to apply to a particular system for analyzing it from an energy perspective. Particular attention has been paid to multi-topology systems because they have not been analyzed previously. Each topology creates a different system that requires, at the assessment stage, an independent analysis.

Table 5 summarizes the recommended equation to be used in each analysis. Integral Energy (*P*) refers to the equation expressed in power terms, its original form, while an (*E)* indicates that the equation has been integrated over the analysis period to calculate the energy required or delivered.

. | Simple systems . | Complex systems (networks) . | Multi-topologic system . | |
---|---|---|---|---|

One topology . | Global . | |||

Assessment | Bernoulli | Integral energy (P or E) | Integral energy (P) | Integral energy (E) |

Audit | Bernoulli | Integral energy (E) | Integral energy (P) | – |

. | Simple systems . | Complex systems (networks) . | Multi-topologic system . | |
---|---|---|---|---|

One topology . | Global . | |||

Assessment | Bernoulli | Integral energy (P or E) | Integral energy (P) | Integral energy (E) |

Audit | Bernoulli | Integral energy (E) | Integral energy (P) | – |

Regardless of the energy equation used, the hydraulic behavior of the system (currently through its mathematic model) is needed to perform the audit. Otherwise, it is impossible to break down the losses.

## CONCLUSIONS

The transport of pressurized water consumes a large amount of energy that should be reduced within the current context of climate change and growing energy costs. Therefore, the first step is to diagnose the system understanding its operational mechanisms and potential areas for energy savings. If there is an important margin for improvement the system should be audited to identify where and how much energy is lost. Both steps can be performed using either the energy integral equation or the Bernoulli equation. While these equations differ in their origins, once adapted, both are suitable to perform these analyses. However, as each equation has its advantages and disadvantages (spatiality and simplicity, respectively), it is advisable to discern the cases for which each equation is best suited. Additionally, any water audit must complement the energy equation with the hydraulic behavior of the system, currently its mathematical model.

This study has placed particular emphasis on the integral energy equation, rarely used to solve hydraulic problems. In particular, it has been applied to assess and audit multi-topologic water transport systems. In addition, the application to a closed simple system has evidenced that the reference elevation level must be the lowest active node (source of consumption) even when the elevation of the lowest node is a tank with a variable level.

The type of system dictates the applicability of the equations. For simple and stationary systems, the Bernoulli equation is accurate and easy to use, making it the recommended equation. Both equations are useful for complex systems, as long as the topology does not change over time. The energy integral equation requires more data but provides more information than the Bernoulli equation. The use of the Bernoulli equation corresponds to a more direct and intuitive approach, but only provides a snapshot of the system, at the moment the equation is applied. This study has established recommendations on when to apply each equation depending on the system under analysis and the phase of the analysis.

## ACKNOWLEDGEMENTS

The authors acknowledge the staff of the water company FACSA for providing financial assistance, helpful advice, and the real case study presented in this work.

## DATA AVAILABILITY STATEMENT

All relevant data are included in the paper or its Supplementary Information.

## CONFLICT OF INTEREST

The authors declare there is no conflict.