Inline Pumped Storage Hydropower towards Smart and Flexible Energy Recovery in Water Networks

: Energy and climate change are thoroughly linked since fossil energy generation highly a ﬀ ects the environment, and climate change inﬂuences the renewable energy generation capacity. Hence, this study gives a new contribution to the energy generation in water infrastructures by means of an inline pumped-storage hydro (IPSH) solution. The selection of the equipment is the ﬁrst step towards good results. The energy generation through decentralized micro-hydropower facilities can o ﬀ er a good solution since they are independent of the hydrologic cycle associated with climate change. The current study presents the methodology and analyses to use water level di ﬀ erence between water tanks or reservoirs in a base pumping system (BPS) to transform it into the concept of a pump-storage hydropower solution. The investigation was developed based on an experimental facility and numerical simulations using WaterGEMS in the optimization of the system operation and for the selection of the characteristic curves, both for the pump and turbine modes. The model simulation of the integrated system was calibrated, and the conceptual IPSH that can be installed was then investigated. The achieved energy for di ﬀ erent technical scale systems was estimated using proper dimensional analysis applied to di ﬀ erent scaled hydraulic circuits, as well as for hydropower response.


Introduction
The increasing need for energy in current societies is inducing more emissions of carbon dioxide to the atmosphere worsening the climate change issues. For that reason, the use of renewable energies has received serious attention in recent years. Ensuring a clean environment and a sustainable development of renewable energy sources widely and globally are appointed as future targets (in UN 2030 Agenda [1]). Although there are great interests in wind and solar as green energy sources, the hydropower should not be overlooked with huge given proof. Currently, hydropower is considered as one of the most flexible and preferred sources to produce electricity [2] and for renewable integration. Therefore, the idea of power production using water based on its available flow energy can contribute to the reduction in significant environmental impacts [3,4].
The basic principle of hydropower is driving a turbine by using the power of water through two common configurations: with or without reservoirs. In hydropower with a reservoir, the water can be stored and is able to generate a considerable amount of energy depending on its capacity, while in hydropower without a reservoir, it produces less, operating preferentially with a constant flow, such as water trunk mains or transmission lines [5]. For that reason, among the diverse use of water sense, this study tried to examine the idea of creating the required head by adding storage tanks to the system. Hence, the current study introduces an inline pumped-storage hydropower (IPSH) solution based on experimental tests, numerical simulations and parametric analysis. Among several discussed practical solutions, IPSH can add more flexibility to pumping systems by providing higher head difference enabling the reduction in energy consumption through reusing the pumped water in gravity branches. Since it will not ask for major changes and large additional investments, using a by-pass to the main pumping system can offer a profitable hydro-energy solution. Among the energy generated, the water power potential energy is more flexible, adaptive and feasible to combine with other renewable sources to feed the pumping system as a hybrid solution. Based on that, this study aims to (a) define the electromechanical selection rules; (b) present a conceptual idea of energy generation by using storage tanks in WDN; (c) analyze results from an experimental small-scaled system by means of numerical simulations to calculate energy generation in a conceptual energy recovery prototype; (d) use suitable dimensional analysis (for hydraulic system and turbomachinery) to predict energy output in different systems' scales. In general, this study presents a low-cost energy prediction using a pumped-storage hydropower solution through experimental measurements, calibration of the numerical model and dimensional analyses.

Characteristic Curves and Operational Point
Pump characteristic curves describe the relationship between the flow rate and the pump head for a specific pump type ( Figure 1). Other important information is graphs for different impeller diameters, the net positive suction head (NPSH), the efficiency and power curves. In the case of any pump, its designation determines its specific nominal discharge impeller diameter. The pump's efficiency throughout its characteristic curve should not drift too much from the best efficiency point (BEP). The motors whose pole number is associated with the rotational speed value also have their own efficiencies to be considered.
Water 2020, 12, x FOR PEER REVIEW 3 of 19 Hence, the current study introduces an inline pumped-storage hydropower (IPSH) solution based on experimental tests, numerical simulations and parametric analysis. Among several discussed practical solutions, IPSH can add more flexibility to pumping systems by providing higher head difference enabling the reduction in energy consumption through reusing the pumped water in gravity branches. Since it will not ask for major changes and large additional investments, using a by-pass to the main pumping system can offer a profitable hydro-energy solution. Among the energy generated, the water power potential energy is more flexible, adaptive and feasible to combine with other renewable sources to feed the pumping system as a hybrid solution. Based on that, this study aims to (a) define the electromechanical selection rules; (b) present a conceptual idea of energy generation by using storage tanks in WDN; (c) analyze results from an experimental small-scaled system by means of numerical simulations to calculate energy generation in a conceptual energy recovery prototype; (d) use suitable dimensional analysis (for hydraulic system and turbomachinery) to predict energy output in different systems' scales. In general, this study presents a low-cost energy prediction using a pumped-storage hydropower solution through experimental measurements, calibration of the numerical model and dimensional analyses.

Characteristic Curves and Operational Point
Pump characteristic curves describe the relationship between the flow rate and the pump head for a specific pump type ( Figure 1). Other important information is graphs for different impeller diameters, the net positive suction head (NPSH), the efficiency and power curves. In the case of any pump, its designation determines its specific nominal discharge impeller diameter. The pump's efficiency throughout its characteristic curve should not drift too much from the best efficiency point (BEP). The motors whose pole number is associated with the rotational speed value also have their own efficiencies to be considered.   The operation point in each pump curve is dependent upon the characteristics of the system in which it is operating. The system head curve is the head equation or the relationship between flow and hydraulic losses in the hydraulic system. Figure 2 shows the pump's operating point (Q 1 , H 1 ) which can change with the differential water level, closure flow control valve or the rotational speed of the pump. The operation point in each pump curve is dependent upon the characteristics of the system in which it is operating. The system head curve is the head equation or the relationship between flow and hydraulic losses in the hydraulic system. Figure 2 shows the pump's operating point ( 1 1 , Q H ) which can change with the differential water level, closure flow control valve or the rotational speed of the pump.
When two or more pumps were installed in parallel, the increasing of flow rates were obtained ( Figure 2). The operating point ( 2 2 , Q H ) represented a higher volumetric flow rate than a single pump for a consequence of greater system head loss, and the volumetric flow rate was lower than twice the flow rate achieved by using a single pump ( 2,1 1 Q Q < ). All of these changes can influence the system efficiency. The affinity laws expressed in Equation (1) represent the mathematical relationship between the rotational speed ( n ), flow rate ( Q ), head ( H ) and pump power ( P ) for the same impeller diameter. The pump specific speed is given by Equation (2).
Variable speed drive (VSD) in pumps induces smooth speed variations of the rotating shaft, directly proportional to the flow, translating into significant pump power variations, which can increase the efficiency of the pumping operation when compared to pumps equipped with fixed speed drive (FSD). Nevertheless, this also affects the pump head, rendering it inoperable below the point where it will not cross the system curve. As possible operating points come closer the operation of the pump becomes unstable if transient regimes occur, inducing the flow variation. Additionally, the pump curve can have a shut-off head inferior to the system curve, meaning that the pump is not able to start at that particular speed.

Selection of a Pump
In a water network system, the total water supply needed by the population must be met by the operation of one or more pumps, whereby the flow rate must be superior to the average daily demand, considering pumps with a flow up to double that value. Regarding the total head, the pumps should comprise a range which takes into account the increase in roughness of the pipes over time such as through Hazen-Williams (H-W) roughness coefficients. Based on the available pump curve, the pump needs to be suitable for its water supply system (WSS) or it will probably operate with reduced efficiency or even with flow instabilities. To have a real idea with a pump efficiency at a BEP of 72% and a motor with 92% efficiency, the total efficiency at that point is 66%. Since the When two or more pumps were installed in parallel, the increasing of flow rates were obtained ( Figure 2). The operating point (Q 2 , H 2 ) represented a higher volumetric flow rate than a single pump for a consequence of greater system head loss, and the volumetric flow rate was lower than twice the flow rate achieved by using a single pump (Q 2,1 < Q 1 ). All of these changes can influence the system efficiency.
The affinity laws expressed in Equation (1) represent the mathematical relationship between the rotational speed (n), flow rate (Q), head (H) and pump power (P) for the same impeller diameter. The pump specific speed is given by Equation (2).
Variable speed drive (VSD) in pumps induces smooth speed variations of the rotating shaft, directly proportional to the flow, translating into significant pump power variations, which can increase the efficiency of the pumping operation when compared to pumps equipped with fixed speed drive (FSD). Nevertheless, this also affects the pump head, rendering it inoperable below the point where it will not cross the system curve. As possible operating points come closer the operation of the pump becomes unstable if transient regimes occur, inducing the flow variation. Additionally, the pump curve can have a shut-off head inferior to the system curve, meaning that the pump is not able to start at that particular speed.

Selection of a Pump
In a water network system, the total water supply needed by the population must be met by the operation of one or more pumps, whereby the flow rate must be superior to the average daily demand, considering pumps with a flow up to double that value. Regarding the total head, the pumps should comprise a range which takes into account the increase in roughness of the pipes over time such as through Hazen-Williams (H-W) roughness coefficients. Based on the available pump curve, the pump needs to be suitable for its water supply system (WSS) or it will probably operate with reduced efficiency or even with flow instabilities. To have a real idea with a pump efficiency at a BEP of 72% and a motor with 92% efficiency, the total efficiency at that point is 66%. Since the operating point is far from the BEP, the operational costs resulting from additional energy consumption can be quite significant.
Then, in a pump system design, several pumps which fitted the considered flow and head ranges can be selected, such as the practical application presented in Figure 3.
operating point is far from the BEP, the operational costs resulting from additional energy consumption can be quite significant.
Then, in a pump system design, several pumps which fitted the considered flow and head ranges can be selected, such as the practical application presented in Figure 3.
1. Pump 1 is not appropriate since it would operate with flow rates not recommended by the manufacturer; 2. Pump 2 has a flow rate near to the average daily demand, increasing its probability of becoming obsolete if the demand is intensified or if the flow is reduced due to a pipe roughness increase over time. The maximum efficiency is also inferior to the pump 3 and, if equipped with a VSD, the possible speed range is minor given its inferior heads. It is also relevant to analyze the speed range in which a pump can operate, where the solution space in the optimization process by the hydraulic simulator WaterGEMS does not include unfeasible solutions, which is determined by comparison to the system curve. The selection has to avoid system instability considering all the former considerations on the pump selection.

Pump as Turbine Curves
When a pump works in the turbine zone, the motor will operate as a generator. During pump operation, the discharge, Q , is a function of the rotating speed, n, and the pumping head, H , whereas the alteration of the speed will depend upon the torque of the motor, T .
For normal turbine operation, the rotating speed ( n ) and discharge (Q ) were negatives and the head ( H ) and torque (T ) were positives ( Figure 4).

1.
Pump 1 is not appropriate since it would operate with flow rates not recommended by the manufacturer; 2.
Pump 2 has a flow rate near to the average daily demand, increasing its probability of becoming obsolete if the demand is intensified or if the flow is reduced due to a pipe roughness increase over time. The maximum efficiency is also inferior to the pump 3 and, if equipped with a VSD, the possible speed range is minor given its inferior heads.
It is also relevant to analyze the speed range in which a pump can operate, where the solution space in the optimization process by the hydraulic simulator WaterGEMS does not include unfeasible solutions, which is determined by comparison to the system curve. The selection has to avoid system instability considering all the former considerations on the pump selection.

Pump as Turbine Curves
When a pump works in the turbine zone, the motor will operate as a generator. During pump operation, the discharge, Q, is a function of the rotating speed, n, and the pumping head, H, whereas the alteration of the speed will depend upon the torque of the motor, T.
For normal turbine operation, the rotating speed (n) and discharge (Q) were negatives and the head (H) and torque (T) were positives ( Figure 4). operating point is far from the BEP, the operational costs resulting from additional energy consumption can be quite significant. Then, in a pump system design, several pumps which fitted the considered flow and head ranges can be selected, such as the practical application presented in Figure 3.
1. Pump 1 is not appropriate since it would operate with flow rates not recommended by the manufacturer; 2. Pump 2 has a flow rate near to the average daily demand, increasing its probability of becoming obsolete if the demand is intensified or if the flow is reduced due to a pipe roughness increase over time. The maximum efficiency is also inferior to the pump 3 and, if equipped with a VSD, the possible speed range is minor given its inferior heads. It is also relevant to analyze the speed range in which a pump can operate, where the solution space in the optimization process by the hydraulic simulator WaterGEMS does not include unfeasible solutions, which is determined by comparison to the system curve. The selection has to avoid system instability considering all the former considerations on the pump selection.

Pump as Turbine Curves
When a pump works in the turbine zone, the motor will operate as a generator. During pump operation, the discharge, Q , is a function of the rotating speed, n, and the pumping head, H , whereas the alteration of the speed will depend upon the torque of the motor, T .
For normal turbine operation, the rotating speed ( n ) and discharge (Q ) were negatives and the head ( H ) and torque (T ) were positives ( Figure 4).  Based on pumps operating as turbines, the turbine characterization was developed by the following next steps: (i) different parameters must be defined based on an established database of synthetic PATs (a large number of different PATs) through characteristic curves that are shown in Figure 5 as well as changing the specific speed (n sT ) defined according to Equation (3): where N R is the rated rotational speed in rpm, P R is the rated power in kW, H R is the rated head in m, since R means the rated conditions or the PAT design point for the best efficiency condition; (ii) when the specific speed is defined, n sT the values of head number (ψ) for each flow rate number (ϕ) can be estimated. When the specific speed was defined and the n s was not known, the values of head number (ψ int ) and efficiency (η int ) for each discharge number value (ϕ int ) were estimated by linear interpolation. When the non-dimensional number was defined, for each diameter and rotational speed (N) the head and efficiency curves were determined by Equations (4)- (6): Water 2020, 12, x FOR PEER REVIEW 6 of 19 Based on pumps operating as turbines, the turbine characterization was developed by the following next steps: (i) different parameters must be defined based on an established database of synthetic PATs (a large number of different PATs) through characteristic curves that are shown in Figure 5 as well as changing the specific speed ( sT n ) defined according to Equation (3): where R N is the rated rotational speed in rpm, R P is the rated power in kW, R H is the rated head in m, since R means the rated conditions or the PAT design point for the best efficiency condition; (ii) when the specific speed is defined, sT n the values of head number (ψ ) for each flow rate number ( ϕ ) can be estimated. When the specific speed was defined and the ns was not known, the values of head number ( int ψ ) and efficiency ( int η ) for each discharge number value ( int ϕ ) were estimated by linear interpolation. When the non-dimensional number was defined, for each diameter and rotational speed ( N ) the head and efficiency curves were determined by Equations (4)- (6): Hence, the correspondent net head and flow rate for the turbine mode will be also influenced by the system curve of the hydraulic system as represented in Figure 5.  Hence, the correspondent net head and flow rate for the turbine mode will be also influenced by the system curve of the hydraulic system as represented in Figure 5.

Methodology
This research used a hydraulic numerical simulation model of an energy recovery system to evaluate the conceptual idea of an IPSH solution and also to estimate the potential energy available in a WSS. This approach used an experimental apparatus of a pumping system to collect the required data of pressure, flow rate, efficiency and rotational speed of a pump for a certain controlling flow. These data were exploited to establish a numerical model using a WaterGEMS simulation tool and also to calibrate the model for different operating conditions. When the numerical model for the pumping system was validated to reproduce the measured data, it was upgraded by adding a by-pass branch to create the desired energy recovery system. The numerical results from the model were used to estimate the energy output of a possible energy recovery solution. This approach has been depicted in the current section by, first, discussing the experimental system and presenting measurement data.
The experimental system is known as the base pumping system (BPS) since it is the base of the future energy recovery system. Then a discussion about the calibration of the numerical model will be provided. This section will be closed by introducing the energy recovery system known as inline pumped-storage hydropower (IPSH) through numerical modelling in a WaterGEMS environment.

System Configuration
The experimental facility BPS is located at the laboratory of hydraulic (LH), Instituto Superior Técnico, Universidade de Lisboa and consists of several components (

Methodology
This research used a hydraulic numerical simulation model of an energy recovery system to evaluate the conceptual idea of an IPSH solution and also to estimate the potential energy available in a WSS. This approach used an experimental apparatus of a pumping system to collect the required data of pressure, flow rate, efficiency and rotational speed of a pump for a certain controlling flow. These data were exploited to establish a numerical model using a WaterGEMS simulation tool and also to calibrate the model for different operating conditions. When the numerical model for the pumping system was validated to reproduce the measured data, it was upgraded by adding a bypass branch to create the desired energy recovery system. The numerical results from the model were used to estimate the energy output of a possible energy recovery solution. This approach has been depicted in the current section by, first, discussing the experimental system and presenting measurement data. The experimental system is known as the base pumping system (BPS) since it is the base of the future energy recovery system. Then a discussion about the calibration of the numerical model will be provided. This section will be closed by introducing the energy recovery system known as inline pumped-storage hydropower (IPSH) through numerical modelling in a WaterGEMS environment.

System Configuration
The experimental facility BPS is located at the laboratory of hydraulic (LH), Instituto Superior Técnico, Universidade de Lisboa and consists of several components (   The measuring range of the mentioned sensors is presented in Table 1. A three-phase motor activated the pump, while the interface control panel provided the possibility of adjustment and measurement of the rotating speed and also the transmitted mechanic torque. A free surface tank with a capacity of 85 dm 3 existed at downstream of the operating system. A valve located at downstream made it possible to induce flow variations by maneuvering and applying local head loss of the valve. An electromagnetic flowmeter (SC-1) and two pressure transducers (SP-1 and SP-2) were used to record the flow rate and pressure data, respectively. Different tests were carried out for different conditions depending on the different opening percentages of the VR-2 valve and also the pump rotational speeds, as presented in Table 2. The data measurements included the flow rate (Q), head (H), rotational speed (N), pump upstream and downstream pressures, the efficiency of the pump (η) and hydraulic and mechanic powers P h and P M , respectively). During the pump operation, the flow rate was a function of the rotational speed and the pumping head [35]. The flow rate varied from zero, for a fully closed valve, to 61.98 L/min for a rotational speed of 2950 rpm, as shown in Figure 7. Additionally, the minimum and maximum measured heads were 1.8 and 16.92 m for rotational speeds of 1600 and 2950 rpm, respectively. Hence, Figure 7 presents characteristic curves for different rotational speeds, flow rate and head, covering the range of operation for the pumped storage system and efficiency variation for N = 2950 rpm. The measuring range of the mentioned sensors is presented in Table 1. A three-phase motor activated the pump, while the interface control panel provided the possibility of adjustment and measurement of the rotating speed and also the transmitted mechanic torque. 1600 to 2950 0 to 100 A free surface tank with a capacity of 85 dm 3 existed at downstream of the operating system. A valve located at downstream made it possible to induce flow variations by maneuvering and applying local head loss of the valve. An electromagnetic flowmeter (SC-1) and two pressure transducers (SP-1 and SP-2) were used to record the flow rate and pressure data, respectively. Different tests were carried out for different conditions depending on the different opening percentages of the VR-2 valve and also the pump rotational speeds, as presented in Table 2. The data measurements included the flow rate (Q ), head ( H ), rotational speed ( N ), pump upstream and downstream pressures, the efficiency of the pump (η ) and hydraulic and mechanic powers h P and M P , respectively). During the pump operation, the flow rate was a function of the rotational speed and the pumping head [35]. The flow rate varied from zero, for a fully closed valve, to 61.98 L/min for a rotational speed of 2950 rpm, as shown in Figure 7. Additionally, the minimum and maximum measured heads were 1.8 and 16.92 m for rotational speeds of 1600 and 2950 rpm, respectively. Hence, Figure 7 presents characteristic curves for different rotational speeds, flow rate and head, covering the range of operation for the pumped storage system and efficiency variation for N = 2950 rpm.

Model Calibration
To perform the system analyses, it was necessary to calibrate the numerical model based on the measured data. The WaterGEMS software from Bentley was used for numerical simulation providing an optimized simulation tool and a user-friendly environment for water distribution networks. WaterGEMS calculated the hydraulic head and pressure at every node along with the flow rate, the flow velocity and the head loss in each pipe branch and as well as the hydraulic head using the gradient algorithm based on the EPANET solver. The water system shown in Figure 8 was built by including one tank, one centrifugal pump and two flow control valves similar to the experimental setup. The experimental setup (shown in Figure 6) is called BPS (base pumping system) in this paper.
Water 2020, 12, x FOR PEER REVIEW 9 of 19

Model Calibration
To perform the system analyses, it was necessary to calibrate the numerical model based on the measured data. The WaterGEMS software from Bentley was used for numerical simulation providing an optimized simulation tool and a user-friendly environment for water distribution networks. WaterGEMS calculated the hydraulic head and pressure at every node along with the flow rate, the flow velocity and the head loss in each pipe branch and as well as the hydraulic head using the gradient algorithm based on the EPANET solver. The water system shown in Figure 8 was built by including one tank, one centrifugal pump and two flow control valves similar to the experimental setup. The experimental setup (shown in Figure 6) is called BPS (base pumping system) in this paper. The calibration process involved the optimization of BPS parameters correspondent to actual measured conditions [35,36]. The BPS was calibrated considering the characteristic curves of the pump, for several pump rotational speeds and valve opening percentages (TCV-2).
Hence, in the simulation process, a variable speed pump (as a VSD) was introduced into WaterGEMS through proper pump curves. The rated pump characteristic curve was defined for the maximum rotational speed, i.e., 2950 rpm, and for each different rotational speed a relative speed factor (RSF) was defined as a coefficient of this maximum rotational speed. A throttle control valve was also defined to induce different flow rates into the system by the partial opening of the valve. In summary, each test was simulated by defining proper RSF and then changing opening flow percentages as in Table 2. This process was repeated for different RSFs until all the rotational speeds were simulated. During this simulation process, the valve discharge coefficients and other associated losses were optimized and calibrated. Results obtained from the simulation by WaterGEMS presented in Figure 9a show good relative accordance with the measured data, associated with scale effects, offering the root mean square errors (RMSEs) shown in Figure 9b with an average of around 0.72. The numerical model calibration based on experimental data guaranteed reliability to follow simulations that take place to assess the behavior of a new adaptation system for energy recovery. A traditional pumping system is composed of different elements that correspond to energy consumption and head losses. However, interesting potential usually exists in pipe branches of pumping systems for energy recovery which can supply energy to treatment plants, electric data base measurement and control devices, and in general, reduce costs of energy in water networks. The calibration process involved the optimization of BPS parameters correspondent to actual measured conditions [35,36]. The BPS was calibrated considering the characteristic curves of the pump, for several pump rotational speeds and valve opening percentages (TCV-2).
Hence, in the simulation process, a variable speed pump (as a VSD) was introduced into WaterGEMS through proper pump curves. The rated pump characteristic curve was defined for the maximum rotational speed, i.e., 2950 rpm, and for each different rotational speed a relative speed factor (RSF) was defined as a coefficient of this maximum rotational speed. A throttle control valve was also defined to induce different flow rates into the system by the partial opening of the valve. In summary, each test was simulated by defining proper RSF and then changing opening flow percentages as in Table 2. This process was repeated for different RSFs until all the rotational speeds were simulated. During this simulation process, the valve discharge coefficients and other associated losses were optimized and calibrated. Results obtained from the simulation by WaterGEMS presented in Figure 9a show good relative accordance with the measured data, associated with scale effects, offering the root mean square errors (RMSEs) shown in Figure 9b with an average of around 0.72. The numerical model calibration based on experimental data guaranteed reliability to follow simulations that take place to assess the behavior of a new adaptation system for energy recovery.

Pumping System Operation
The best efficiency point (  (Table 3). The characteristic curves in the dimensionless form ( Figure 10) were constructed based on a rated condition associated with the best efficiency point. The dimensionless curves provided a tool to transfer information to other equivalent systems. A traditional pumping system is composed of different elements that correspond to energy consumption and head losses. However, interesting potential usually exists in pipe branches of pumping systems for energy recovery which can supply energy to treatment plants, electric data base measurement and control devices, and in general, reduce costs of energy in water networks.

Pumping System Operation
The best efficiency point (Q R , H R , η R ) of the pump was characterized by a flow rate of 35 L/min, with a head of 13.5 m and an efficiency of 75% for the rotational speed of 2950 rpm and specific speed 10.05 (Table 3). The characteristic curves in the dimensionless form ( Figure 10) were constructed based on a rated condition associated with the best efficiency point. The dimensionless curves provided a tool to transfer information to other equivalent systems. The best efficiency point of the pump was selected to present the head and efficiency variations with a flow rate, for the rated rotational speed of 2950 rpm, in Figure 11.

Inline Pumped-Storage Hydropower (IPSH)
In this stage, the numerical model was ready to be upgraded to an energy recovery system. The improvement was performed by adding a by-pass line, as shown in Figure 12. In order to examine the IPSH solution, a free surface tank (T-2) was installed in the by-pass branch that was added to the base pumping system (BPS) with the capacity of 0.032 m 3 . The by-pass line was equipped with a throttle control valve (TCV) working in an open or closed position to include or isolate the by-pass line. The T-2 was located at a lower elevation to generate a gravity flow from T-1 to T-2. The new IPSH system was considered as a loop system to use the previously assessed characteristics attained in the experimental loop system. It is worth mentioning that the idea is not limited to loop systems but can be adapted to a real system with direct flow condition. In a direct flow system, based on the available head at T-1 and downstream demand, the by-pass line can be activated to use the available head difference for energy generation. Hence, a turbine was considered in the by-pass line to generate energy from the gravity flow. Since WaterGEMS does not include a built-in turbine element, the general purpose valve (GPV) was used for this purpose by defining the flow-head loss curve correspondent to the turbine characteristic curves. The best efficiency point of the pump was selected to present the head and efficiency variations with a flow rate, for the rated rotational speed of 2950 rpm, in Figure 11. The best efficiency point of the pump was selected to present the head and efficiency variations with a flow rate, for the rated rotational speed of 2950 rpm, in Figure 11.

Inline Pumped-Storage Hydropower (IPSH)
In this stage, the numerical model was ready to be upgraded to an energy recovery system. The improvement was performed by adding a by-pass line, as shown in Figure 12. In order to examine the IPSH solution, a free surface tank (T-2) was installed in the by-pass branch that was added to the base pumping system (BPS) with the capacity of 0.032 m 3 . The by-pass line was equipped with a throttle control valve (TCV) working in an open or closed position to include or isolate the by-pass line. The T-2 was located at a lower elevation to generate a gravity flow from T-1 to T-2. The new IPSH system was considered as a loop system to use the previously assessed characteristics attained in the experimental loop system. It is worth mentioning that the idea is not limited to loop systems but can be adapted to a real system with direct flow condition. In a direct flow system, based on the available head at T-1 and downstream demand, the by-pass line can be activated to use the available head difference for energy generation. Hence, a turbine was considered in the by-pass line to generate energy from the gravity flow. Since WaterGEMS does not include a built-in turbine element, the general purpose valve (GPV) was used for this purpose by defining the flow-head loss curve correspondent to the turbine characteristic curves.

Inline Pumped-Storage Hydropower (IPSH)
In this stage, the numerical model was ready to be upgraded to an energy recovery system. The improvement was performed by adding a by-pass line, as shown in Figure 12. In order to examine the IPSH solution, a free surface tank (T-2) was installed in the by-pass branch that was added to the base pumping system (BPS) with the capacity of 0.032 m 3 . The by-pass line was equipped with a throttle control valve (TCV) working in an open or closed position to include or isolate the by-pass line.
The T-2 was located at a lower elevation to generate a gravity flow from T-1 to T-2. The new IPSH system was considered as a loop system to use the previously assessed characteristics attained in the experimental loop system. It is worth mentioning that the idea is not limited to loop systems but can be adapted to a real system with direct flow condition. In a direct flow system, based on the available head at T-1 and downstream demand, the by-pass line can be activated to use the available head difference for energy generation. Hence, a turbine was considered in the by-pass line to generate energy from the gravity flow. Since WaterGEMS does not include a built-in turbine element, the general purpose valve (GPV) was used for this purpose by defining the flow-head loss curve correspondent to the turbine characteristic curves. The numerical simulation was used to assess the main variations in the IPSH system. The characteristic curves in Figure 13 were used to adjust the flow rate in the main pipe branch by changing the rotational speed of the pump and the TCV-2 opening based on the flow rate in gravity by-pass line for energy recovery. The simulations were carried out for an extended period of 24 h. Two scenarios were considered, i.e., identical and variable flow rates in the pipe system. If the flow rates in different branches of the system were equal, it led to a steady flow regime resulting in a constant water level in T-1 and T-2 ( Figure 13). In this case, a gravity flow rate of 29.11 L/min existed from T-1 to T-2. To maintain this flow rate in the main pipe, a pump rotational speed of 2600 rpm and a TCV-2 opening of 72.5% were found to be appropriate based on the experimental measurements. In other words, the characteristic curves gave the ability to accurately adjust the pump working condition and valve opening percentage in order to establish a flow rate equal to the gravity flow between T-1 and T-2. Based on that, the water levels in T-1 and T-2 for the identical flow rate scenario over 24 h remained constant, as 0.50 and 0.18 m, respectively.  The numerical simulation was used to assess the main variations in the IPSH system. The characteristic curves in Figure 13 were used to adjust the flow rate in the main pipe branch by changing the rotational speed of the pump and the TCV-2 opening based on the flow rate in gravity by-pass line for energy recovery. The simulations were carried out for an extended period of 24 h. Two scenarios were considered, i.e., identical and variable flow rates in the pipe system. If the flow rates in different branches of the system were equal, it led to a steady flow regime resulting in a constant water level in T-1 and T-2 ( Figure 13). In this case, a gravity flow rate of 29.11 L/min existed from T-1 to T-2. To maintain this flow rate in the main pipe, a pump rotational speed of 2600 rpm and a TCV-2 opening of 72.5% were found to be appropriate based on the experimental measurements. In other words, the characteristic curves gave the ability to accurately adjust the pump working condition and valve opening percentage in order to establish a flow rate equal to the gravity flow between T-1 and T-2. Based on that, the water levels in T-1 and T-2 for the identical flow rate scenario over 24 h remained constant, as 0.50 and 0.18 m, respectively. The numerical simulation was used to assess the main variations in the IPSH system. The characteristic curves in Figure 13 were used to adjust the flow rate in the main pipe branch by changing the rotational speed of the pump and the TCV-2 opening based on the flow rate in gravity by-pass line for energy recovery. The simulations were carried out for an extended period of 24 h. Two scenarios were considered, i.e., identical and variable flow rates in the pipe system. If the flow rates in different branches of the system were equal, it led to a steady flow regime resulting in a constant water level in T-1 and T-2 ( Figure 13). In this case, a gravity flow rate of 29.11 L/min existed from T-1 to T-2. To maintain this flow rate in the main pipe, a pump rotational speed of 2600 rpm and a TCV-2 opening of 72.5% were found to be appropriate based on the experimental measurements. In other words, the characteristic curves gave the ability to accurately adjust the pump working condition and valve opening percentage in order to establish a flow rate equal to the gravity flow between T-1 and T-2. Based on that, the water levels in T-1 and T-2 for the identical flow rate scenario over 24 h remained constant, as 0.50 and 0.18 m, respectively.  In a water supply system, there was a pattern of flow demand along each 24 h. This pattern had a typical representation of each system characterization depending on the flow consumption used for water network design. Therefore, in this research, the same procedure was adopted since the flow pattern varied along the time, between rush consumption or peak hours to fewer consumptions, normally associated with the night period-the one which was also used for leak detection, since the level of consumption attained the minimum. Hence, a suitable design project for this extended period of 24 h represents a more granted solution to face flow and water level variations, head losses, leakages occurrence associated with high pressure values, along the hydraulic system and machine operation adaptation, which fit both results as the operating point. This is a complex issue that requires an extended period and is typically used in the design of water systems. Under some operating conditions, the flow rate in different branches can be variable and unequal. Then, to evaluate this scenario, a flow rate pattern of 24 h was adapted, based on a typical demand configuration, as presented in Figure 14a. A TCV-2 valve closure pattern was then calibrated to induce the desired flow rate in the hydraulic system (Figure 14b). In a water supply system, there was a pattern of flow demand along each 24 h. This pattern had a typical representation of each system characterization depending on the flow consumption used for water network design. Therefore, in this research, the same procedure was adopted since the flow pattern varied along the time, between rush consumption or peak hours to fewer consumptions, normally associated with the night period-the one which was also used for leak detection, since the level of consumption attained the minimum. Hence, a suitable design project for this extended period of 24 h represents a more granted solution to face flow and water level variations, head losses, leakages occurrence associated with high pressure values, along the hydraulic system and machine operation adaptation, which fit both results as the operating point. This is a complex issue that requires an extended period and is typically used in the design of water systems. Under some operating conditions, the flow rate in different branches can be variable and unequal. Then, to evaluate this scenario, a flow rate pattern of 24 h was adapted, based on a typical demand configuration, as presented in Figure 14a. A TCV-2 valve closure pattern was then calibrated to induce the desired flow rate in the hydraulic system (Figure 14b).
Despite the previous case of having a balance of flow rate in the whole system and constant water levels in tanks, the water level in both tanks changed with time, as shown in Figure 15. The water level variations in T-1 and T-2 were correlated, decreasing in T-1, in turbine mode and increasing in T-2, in pump mode, in a controlled optimized way between the maximum and minimum limits for tank water levels ( Figure 15).  Despite the previous case of having a balance of flow rate in the whole system and constant water levels in tanks, the water level in both tanks changed with time, as shown in Figure 15. The water level variations in T-1 and T-2 were correlated, decreasing in T-1, in turbine mode and increasing in T-2, in pump mode, in a controlled optimized way between the maximum and minimum limits for tank water levels ( Figure 15). Water 2020, 12, x FOR PEER REVIEW 14 of 19 Figure 15. The water level in T-1 and T-2 based on the variable demand pattern.
The operating curve of the GPV valve (acting as a turbine) was calculated with the available gross head and total head losses to obtain the net head of the turbine, as presented in the dimensionless graph of Figure 13. Equation (7) was used to calculate turbine power: where is the power, is the specific weight of the water, is the flow rate, is the turbine head and is the efficiency. Equation (7) was exploited to calculate the power for two different mentioned scenarios of constant and variable flow rates, as presented in Figure 16. The energy production of the system for a fixed flow rate of 35 L/min was 1.24 kWh, while the variable scenario led to 1.81 kWh daily energy production.

Dimensional Analysis and Discussion
Experimental tests are used to predict the performance of a real device including the behavior under different operating conditions. In this case, the model was evaluated on a laboratory scale, and then the results were simulated for real conditions. The similarity theory requires principles of geometric, kinetic and dynamic parities between a model and a prototype. These affinity laws in different categories can be attained by expressing shape, size, velocity component and acting forces [37]. In this study, the scales of velocity and flow rate were calculated based on the Froude criterion, as analyzed in [38], from the length scale ( ) using the following equations: The operating curve of the GPV valve (acting as a turbine) was calculated with the available gross head and total head losses to obtain the net head of the turbine, as presented in the dimensionless graph of Figure 13. Equation (7) was used to calculate turbine power: where P is the power, γ is the specific weight of the water, Q is the flow rate, H is the turbine head and η is the efficiency. Equation (7) was exploited to calculate the power for two different mentioned scenarios of constant and variable flow rates, as presented in Figure 16. The energy production of the system for a fixed flow rate of 35 L/min was 1.24 kWh, while the variable scenario led to 1.81 kWh daily energy production. The operating curve of the GPV valve (acting as a turbine) was calculated with the available gross head and total head losses to obtain the net head of the turbine, as presented in the dimensionless graph of Figure 13. Equation (7) was used to calculate turbine power: where is the power, is the specific weight of the water, is the flow rate, is the turbine head and is the efficiency. Equation (7) was exploited to calculate the power for two different mentioned scenarios of constant and variable flow rates, as presented in Figure 16. The energy production of the system for a fixed flow rate of 35 L/min was 1.24 kWh, while the variable scenario led to 1.81 kWh daily energy production.

Dimensional Analysis and Discussion
Experimental tests are used to predict the performance of a real device including the behavior under different operating conditions. In this case, the model was evaluated on a laboratory scale, and then the results were simulated for real conditions. The similarity theory requires principles of geometric, kinetic and dynamic parities between a model and a prototype. These affinity laws in different categories can be attained by expressing shape, size, velocity component and acting forces [37]. In this study, the scales of velocity and flow rate were calculated based on the Froude criterion, as analyzed in [38], from the length scale ( ) using the following equations:

Dimensional Analysis and Discussion
Experimental tests are used to predict the performance of a real device including the behavior under different operating conditions. In this case, the model was evaluated on a laboratory scale, and then the results were simulated for real conditions. The similarity theory requires principles of geometric, kinetic and dynamic parities between a model and a prototype. These affinity laws in different categories can be attained by expressing shape, size, velocity component and acting forces [37].
In this study, the scales of velocity and flow rate were calculated based on the Froude criterion, as analyzed in [38], from the length scale (λ L ) using the following equations: Flow rate scale : Velocity scale : Classic similarity laws for pumps and turbines with the same impeller diameter state that the discharge is proportional to the rotational speed, while the head is proportional to the squared rotational speed as Equation (1). Then based on [39,40], the correlation between pump and turbine mode was proved to be: n S T (in m, m 3 /s) = 0.8793n S P (11) To evaluate the results of the model in a technical approach, for a loop system where Q T = Q P the length scales of 20 and 50 were considered. Additionally, for geometrically similar impellers operating at the same specific speed, the affinity laws are as follows:  (13) or separating by specific flow, head and power: The results are presented in Table 4. The power of the model was 0.052 kW but, by upscaling, it grew to values of 1532 and 41,657 kW, for scales of 20 and 50, respectively. Table 4. Scale-up parameters for a hydraulic system and turbomachine affinity characteristic parameters with n sT (in m, m 3 /s) = 8.8.

Hydraulic System
Turbine Impeller Turbine Affinity Laws

Conclusions
In a base pumping system (BPS), the characteristic curves in pump and turbine modes were defined and analyzed in order to obtain the best operating conditions in water networks. The research study included experimental analyses, hydraulic simulations and optimized conditions to better define the best operating point. The experimental results were exploited to calibrate a numerical hydraulic simulator model in the WaterGEMS environment, which uses optimization in searching for