Next Article in Journal
Religious Participation: Does It Matter for Sustainable Culture and Entertainment Consumption?
Next Article in Special Issue
Hydrostatic Pressure Wheel for Regulation of Open Channel Networks and for the Energy Supply of Isolated Sites
Previous Article in Journal
The Influence of Tilia tomentosa Moench on Plant Species Diversity and Composition in Mesophilic Forests of Western Romania–A Potential Tree Species for Warming Forests in Central Europe?
Previous Article in Special Issue
Inner Damping of Water in Conduit of Hydraulic Power Plant
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Numerical Investigation into the PAT Hydrodynamic Response to Impeller Rotational Speed Variation

1
College of Water Conservancy and Hydropower Engineering, Hohai University, Nanjing 210098, China
2
College of Energy and Electrical Engineering, Hohai University, Nanjing 210098, China
3
College of Agricultural Science and Engineering, Hohai University, Nanjing 210098, China
4
College of Petroleum Engineering, Liaoning Shihua University, Fushun 113001, China
*
Author to whom correspondence should be addressed.
Sustainability 2021, 13(14), 7998; https://doi.org/10.3390/su13147998
Submission received: 21 June 2021 / Revised: 11 July 2021 / Accepted: 12 July 2021 / Published: 17 July 2021

Abstract

:
The utilization of pump as turbines (PATs) within water distribution systems for energy regulation and hydroelectricity generation purposes has increasingly attracted the energy field players’ attention. However, its power production efficiency still faces difficulties due to PAT’s lack of flow control ability in such dynamic systems. This has eventually led to the introduction of the so-called “variable operating strategy” or VOS, where the impeller rotational speed may be controlled to satisfy the system-required flow conditions. Taking from these grounds, this study numerically investigates PAT eventual flow structures formation mechanism, especially when subjected to varying impeller rotational speed. CFD-backed numerical simulations were conducted on PAT flow under four operating conditions (1.00 QBEP, 0.82 QBEP, 0.74 QBEP, and 0.55 QBEP), considering five impeller rotational speeds (110 rpm, 130 rpm, 150 rpm, 170 rpm, and 190 rpm). Study results have shown that both PAT’s flow and pressure fields deteriorate with the machine influx decrease, where the impeller rotational speed increase is found to alleviate PAT pressure pulsation levels under high-flow operating conditions, while it worsens them under part-load conditions. This study’s results add value to a thorough understanding of PAT flow dynamics, which, in a long run, contributes to the solution of the so-far existent technical issues.

1. Introduction

The use of fossil fuels as main energy sources goes with a number of shortfalls, environmental disruption, and possible reserves for depletion, among others [1]. Fossil fuels are linked to the widely condemned greenhouse gas emission, i.e., gases which are believed to be at the source of climate change and other detrimental effects [2]. Of course, the non-renewable nature of fossil fuels itself poses a threat of resources exhaustion in the near or far future [3]. With the international will on reducing the CO2 emissions, renewable energy sources (REN) have been advocated, and countries have taken action to not only involve RENs in their power systems, but to continue increasing RENs annual productions [4,5]. Among the so far available RENs, hydropower is one of the most mature REN technologies that dominates the market with a multitude of hydrosites already expanding the whole world over [6,7]. Hydropower technology is linked to different advantages when compared to other REN technologies, including their huge energy storage ability and operational flexibility, among others. Depending on different aspects, the targeted production and site-available potential among others, hydropower technology exploits natural waters from all sorts, such as seas, lakes, rivers, and estuaries. This also somehow leads to the classification of hydrosites as large- or small-scale hydrosites [8]. While large hydropower plants are generally known to feed large-size power grids such as the national or even international grids, small hydropower is mostly used under off-grid operational mode, contributing to the electrification of rural areas which not only provides electric lighting for households but also contributes to other economic activities such as the agricultural irrigation [9,10,11].
With the continued research on hydropower development, water supply and irrigation systems have been found to be another alternative for hydropower production on small scales [12]. In these systems, water flows continuously, thus presenting a stable hydropower potential. Besides, system excess energy in terms of pressure is usually wasted through pressure-reducing valves (PRVs) to avoid water leakages [13]. The installation of hydraulic turbines in these networks not only replaces the formerly utilized valves in regulating the system’s energy levels but also uses that same energy to produce usable electric power [14,15,16,17].
Now, as also mentioned by different investigators, for small-scale hydropower (SHP) projects, the initial installation budget is the most expensive as the electromechanical equipment takes up a considerably high percentage of the total investments, in the range from 35% to 40% or even beyond [18,19]. In line with this, the utilization of normal hydraulic pumps in reverse mode, as an economic alternative to SHP installation budget, has recently attracted the attention of many REN field players [20,21]. In addition to its cost-effectiveness, pump as turbines (PATs) present different advantages over comparable mini-scale hydroturbines, including the easy availability of the machine itself and spare parts, structural simplicity, easy maneuverability, short delivery time, and a long life span [22,23].
Nevertheless, PAT technology has also shown different shortfalls. Among others, the so far controversial conversion of pump-mode performance characteristics to turbine mode is the most researched, where the so far available analytical conversion methods present errors in the range of ±20% [24]. Scientific research on this matter still goes on as of today. For more about PAT turbine mode performance predictions, review articles by Jain et al. [24], Elbatran et al. [25], Adu et al. [26], and Ntiri et al. [27] are recommended.
Another problem about PAT technology resides in its structural design that does not necessarily meet the energy production requirements. This aspect itself plays a big role in the above-mentioned difficulty in PAT’s turbine mode performance prediction. To be more specific, the pump impeller blades are curved in the opposite direction of its rotation (backward). This blade design becomes forward-curved when the flow direction changes in turbine mode, thus serving the source to large flow separations and associated hydraulic losses as a result of flow impingement on the blade’s pressure side in the vicinities of its leading edge. In addition, hydraulic pumps are not necessarily provided with flow control mechanism. Therefore, the machine is unable to regulate the flow rate or the water flow incidence angle to satisfy the demanded operating conditions. This also explains its inability to handle off-design operating conditions. Therefore, PATs have been found to optimally function in a tiny range of flows around the best efficiency point (BEP), while they operate poorly under off-design conditions [28]. This constitutes one of the biggest shortfalls of PAT technology and has attracted the attention of many researchers, mainly targeting the understanding of PAT flow dynamics and subsequently attempting to widen its optimum operating range. Among the recently published findings, one would cite a study conducted by Štefan et al. [29] on PAT internal flow field under both pump and turbine modes. In total agreement with the literature, the turbine mode was found to attain its BEP under a higher head and flow rate, as compared to the pumping mode (41% and 27%, respectively), while the efficiency was slightly lower (ηT = 0.96 ηP). Moreover, a steep drop on the PAT efficiency curve has been noticed for deep part-load operating conditions, supposedly owing to the intensification of swirling flows within the PAT’s outlet pipe. Based on the fact that, for fluid machinery, adverse flow phenomena mostly occurring under unstable conditions led to increased energy loss rates, Lin et al. [30] has used the entropy dissipation theory to study PAT operating characteristics under off-design operating conditions. As expected, unsteady flow phenomena such as the backflows, shock flows, and vortical flows, have inflicted higher entropy and energy dissipation, where flow zones with the highest energy losses have been pointed out as the vicinities of the blade leading edge, the cavity zone (impeller’s back and front chambers), and volute throat. In the same respect, Guan et al. [31] and Morabito et al. [32] found that the position and geometry of the volute throat may influence the machine energy loss characteristics and its efficiency.
More studies have also been conducted to try and decrease PAT hydraulic losses through different technics, geometry modification, and optimization, among others. To say the least, Ghorani et al. [33,34], after finding that irreversible hydraulic losses associated with flow shock at the blade leading edge and other secondary flow phenomena such as backflow and vortices take place within PAT flow zones, have opted to optimizing the PAT impeller geometry, leading to substantial reduction in the above-stated unsteady flow phenomena and a corresponding increase in the machine hydraulic efficiency for different operating conditions (3.8% and 5.72% for QBEP and 1.3 QBEP flow conditions, respectively). In the same respect, Wang et al. [35] mentioned that one of most cost-effective ways to improve PAT’s turbine mode performance would be by redesigning its blade geometry without touching other components. Therefore, they designed a special impeller with forward-swept blades and compared its performance with the original back-swept bladed one. The special impeller, owing to the associated considerable reduction in impeller hydraulic and shock losses, exhibited a flatter efficiency-flow rate curve (η-Q curve), thus promising a larger range of optimum operating conditions around the BEP (1.5% efficiency variation for flow conditions within the range from 0.9 QBEP to 1.2 QBEP). Moreover, a 7.9% increase in the maximum attainable efficiency was achieved, thus confirming the special impeller’s superiority over the original one. In the same respect, Sengpanich et al. [36] proposed “an impulse pump-as-turbine” using an injector as a flow control mechanism. This technique also assured a continuously high efficiency for a wide range of flows.
Coming back to the application of PAT technology in water distribution networks (WDNs), different studies have so far been carried out [25,37,38,39,40], which concern the above-discussed PAT technology-associated problems. Consequently, the current technology readiness level (TRL) for PAT technology is estimated at TRL−4 [41,42]. The utilization of PATs in WDNs must consider the variability of pressure and flow rate as imposed by the system, due to the instantaneity in water demands [43,44,45]. In this respect, Carraveta et al. [46,47,48] have developed a so-called “variable operating strategy” or VOS, with the aim of efficiently using PAT technology to regulate the WDN energy. This strategy consists of two distinct regulation methods, namely hydraulic regulation (HR) and electric regulation (ER), though they can also be combined to give new regulation features. The HR method consists of the usage of PAT in a series-parallel circuit with PRVs where the system energy regulation is mainly achieved by just acting on the utilized PRVs. The valve in series with PAT is used to dissipate the surplus energy leading to a lower system flow as compared to PAT’s BEP flow (QS < QT), while the valve in parallel with PAT bypasses a quantity of flow thus shifting the system operating point to higher flows (QS > QT). As for the ER method, an electric inverter is used to regulate the impeller rotational speed (IRS) where, for instance, one would increase the IRS to shift the system operating point to lower flows, or decrease it to obtain the system operating point to higher flow conditions. In this method, however, owing to a comparatively narrow range of system operating conditions between the runaway (T = 0 Nm) and resistance (RRS = 0 rpm) characteristic curves, high IRS values may lead to a system operating point where the runner rotational speed is close to the runaway value, possibly causing PAT’s poor performance. To avoid this situation, a series-PRV can be used to shift the IRS to lower values, while respecting the system-required operating point. VOS functionality discussions and applications have been presented through a number of studies by Fecarrota et al. [49] and Stefanizzi et al. [41], among others.
Specifically, the utilization of variable-speed PATs in WDNs as a means to control the flow while keeping a good system operational efficiency has so far been researched by different authors, including Delgado et al. [50], Lima et al. [51], and Tahani et al. [52], among others. In most of these studies, the main focus has been the optimization of selection, position, and utilization of variable-speed PATs, targeting the improved WDN energy regulation and efficient energy recovery. On the other hand, a few more studies have been conducted on the effect of PAT speed variation, but considering its application to normal hydropower production within small-scale hydroelectric plants, where the machine external performance characteristics have been the focus [53,54]. What is also worth mentioning is a study on the effect of PAT rotational speed variation on tip leakage vortex (TLV) dynamics under pump mode, as conducted by Han et al. [55], where the TLV was found to worsen with an increase in PAT impeller rotational speed. To the authors’ best knowledge; literature is so far quite limited on the impact of impeller rotational speed variation on PAT flow dynamics under turbine operating mode.
Therefore, by targeting a thorough understanding of PAT flow dynamics under these circumstances, especially considering PAT’s off-design operating conditions, the present article presents a numerical investigation on the formation mechanism of PAT flow structures, and the effect of impeller rotational speed variation on the same. CFD-backed numerical simulations have been conducted on PAT’s 3D turbulent flow under optimal condition and three part-load conditions (1.00 QBEP, 0.82 QBEP, 0.74 QBEP, and 0.55 QBEP), where five impeller rotating speeds, namely 110 rpm, 130 rpm, 150 rpm, 170 rpm, and 190 rpm, were considered.
The investigated model is a mixed flow PAT with a three-bladed impeller and a flow control system of seven guide vanes. This article is composed of three main sections. Section 1 provides an extensive background on PAT technology and the associated techno-scientific issues, as well as the so far attained level of knowledge. Section 2 describes both the investigated case and utilized methods. Section 3 provides a thorough discussion on the present study’s findings, which finally led to a number of concluding remarks.

2. Materials and Methods

In order to achieve this study’s objective, a roadmap was designated and followed, as shown by the flowchart in Figure 1. The PAT prototype was first experimentally tested under various operating conditions to finally obtain the pump as turbine’s external performance characteristics in terms of Head (H) and Efficiency (η).
Then, a PAT geometrical model was built, from which a predominantly hexahedral computational grid was later generated. The resultant grid was taken to ANSYS CFX-Solver (developed by Ansys, Inc., an American company based in Canonsburg, Pennsylvania) for PAT flow simulation under various operating conditions, as also mentioned in the last sections of the article.
During simulation, the experimental test-extracted data were used as system boundary conditions, and the utilized numerical scheme validation was later performed by comparing experimental and numerical simulation results. Finally, farther study results in terms of machine flow and pressure field characteristics with gradually changing machine influx (Q) and impeller rational speed (n) were extracted, and deeply analyzed, leading to a number of concluding remarks.

2.1. PAT Geometry

The research object of this study is a mixed-flow pump as turbine (PAT) composed of four main components, namely the inlet pipe (IP), the impeller, the guide vanes (GV) zone, and the outlet pipe (OP). The impeller had 3 blades (ZR = 3) while the guide vanes zone was composed of 7 guide vanes (ZG = 7). Both the inlet and outlet pipes were split type tubes. The model’s specific speed (NQ) was 163.9 min−1. The impeller diameter was 320 mm, while its outlet diameter was 300 mm. The guide vanes inlet diameter was 350 mm, and there was no existent blade tip clearance. A corresponding computational domain was built and simulated using Unigraphix NX8.0 (owned by Siemens PLM, an American company based in Plano, TX, USA) and Ansys CFX softwares, respectively. Figure 2 shows the utilized PAT computational domain and its components.

2.2. Grid Generation

Grid generation constitutes one of very important phases of the numerical simulation process. The generated grid’s quality is a very crucial parameter, as it plays a fundamental role on different aspects, numerical solution convergence, and accuracy within numerical simulation results, among others. Therefore, great attention is always paid to this phase. Depending on different factors such as the geometry complexity or the nature of simulated system’s operations, grid characteristics can be varied at different flow zones.
For the case of an axial-flow PAT, which is this study’s case, flow separations are likely to take place in rotating parts and the adjacent flow zones. Therefore, for the sake of accurately capturing the flow dynamics at those zones, a finer grid was created around the impeller blades and guide vanes, leading to a global y+ value less than 50. Ansys ICEM software was used to generate this grid, where the hexahedral mesh type was utilized for all of the components of the computational domain. In order to avoid the possibility of numerical simulation results being influenced by the utilized grid number, a grid independence test was carried out where six different grid numbers ranging from 4 to 8 million were tested. Having been simulated under the same boundary conditions, different grid numbers have led to different values of PAT’s hydraulic head. Table 1 displays the head variation with the increasing grid number from 4 to 8 million elements. It can be seen that the head increased with the grid number until 7.2 million, after which the head value approximately stabilized. Therefore, in line with the available computational resources, a grid number of 7.2 million was chosen for further transient simulations. Details of the selected grid are presented in Table 2.

2.3. Numerical Simulation Scheme

The water flow that passes through a PAT’s full computational domain is three dimensional. This study considers the working fluid as the water at 25 °C. It is considered incompressible, isothermal, and is governed by the Reynolds-averaged Navier-Stokes equations, otherwise known as RANS. The RANS-associated continuity and momentum equations are shown in Equations (1) and (2).
Continuity equation:
x i ( ρ u i ) = 0
Momentum equation:
t ( ρ u j ) + x j ( ρ u j u i ) = p x j + τ i j x j + S i j
In these equations, ui and uj stand for velocity components, τij stand for shear stress, Sij for the source term, p for pressure, and ρ for density. In this study, the shear stress transport (SST) computational model was selected for the closure of the RANS equations [56]. This model combines the advantages of the renowned k-ε and k-ω models [57,58], thus giving it superiority over both of them. This is due to its ability to deliver high-accuracy predictions for flow dynamics in the wall vicinities and flow zones away from it, especially when solving for study cases involving adverse pressure gradients and flow separations [59]. Equations (3)–(5) are the mathematical expressions of SST’s turbulent kinetic energy, its specific dissipation rate, and the kinematic eddy viscosity, respectively.
k t + u j k x j = P k β k ω + x j ν + σ k ν T k x j
ω t + u j ω x j = α S 2 β ω 2 + x j ν + σ ω ν T ω x j + 2 1 F 1 σ ω , 2 1 ω k x i ω x i
ν T = a 1 k max a 1 ω , S F 2
The 3D numerical simulation of PAT flow under different operating conditions is conducted using the CFD’s commercial code Ansys CFX 18.0. Simulations were set in such a way that the water flow entered the computational domain at the inlet zone of the inlet pipe (see Figure 1), then passed through the guide vanes inter-spaces before whirling through the impeller inter-blade channels, and flowed down to the outlet pipe to finally exit the computational domain through the outlet pipe’s outlet zone.
For the sake of a quick solution convergence, steady state numerical simulations were initially run, where their results files served the initial states for the followed transient simulations. A constant mass flow (kg s−1) was imposed at the inlet of the computational domain as the inlet boundary condition, while a gauge pressure (Pa) was selected as the outlet boundary conditions. Different interfaces were utilized between components of the computational domain, depending on the location and simulation type. The frozen rotor interface was used at both ends of the impeller under steady state numerical simulations, while it was changed to transient rotor-stator interface for transient numerical simulations. The general grid interface (GGI) was used between the stationary components for both simulation types. All the remaining solid walls were given a non-slip boundary condition.
The length of a single simulation session was set to ten impeller revolutions, where the selected timestep (∆t) was equivalent to the impeller’s one degree rotation (1°), which therefore means that a single impeller revolution was equivalent to 360 × ∆t. Five internal loops were selected for each timestep. The solution convergence criteria was chosen as the root mean square of residuals inferior to 1 × 10−4. For the sake of checking the trustworthiness of generated simulation results by the above-presented numerical scheme, experimental and numerical simulation results in terms of machine external performance characteristics were compared.
Experiments were conducted on a PAT test rig, as shown in Figure 3. In this figure, the test rig components are shown on the right side, and their respective positions on a fully assembled PAT system are shown by letters A, B, and C. The test rig had the highest head (H) and flow rate (Q) of 100 m and 1000 L/s, respectively. The adoptable impeller speed range was 0 to 1500 rpm, which could accommodate impeller diameters in the range from 250 to 460 mm.
As shown in Figure 4, a global error of less than 8% and 6.5% for head and efficiency predictions was recorded, respectively. The overestimation of experimental results by the utilized numerical scheme is believed to take source for different reasons, including the computational geometry simplification, among others. Nevertheless, the utilized numerical scheme showed the ability to generate reasonably accurate results, thus proving its trustworthiness.

3. Results

3.1. Flow Field Characteristics Evolution within the Impeller and Guide Vanes

Though the experimental tests were conducted on a wide range of flows, this study has only considered four operating conditions, namely 1 QBEP, 0.82 QBEP, 0.74 QBEP, and 0.55 QBEP, considering a constant impeller rotational speed of 150 rpm. This was carried out with a target of investigating the eventual changes in machine flow dynamics for a gradually decreasing machine influx. To do so, three consecutive and equidistant turbo-surfaces were selected in the streamwise direction from the impeller inlet to the outlet.
Considering the impeller inlet and outlet positions as 0 and 1, respectively, the 1st plane (S1) was at 0.2, the 2nd (S2) at 0.5, and the 3rd (S3) at 0.8. Each of the three planes expanded from the impeller hub to the shroud. During the investigation of local flow parameter, changes took place as the machine flow rate decreased, and the above-mentioned surfaces were taken as the reference flow zones within the impeller. In this respect, an attempt was completed to showcase some of these changes through both Figure 5 and Figure 6.
Figure 5 shows the pressure distribution mode and associated flow pattern within the impeller flow zone under 0.55 QBEP. In this figure, from left to right, are plane S3 (impeller outlet vicinal zone), plane S2 (impeller most central zone), and plane S1 (impeller inlet vicinal zone).
It can be seen that the flow region in the vicinities of the impeller inlet (Plane S1) was marked by considerably low-pressure zones on the blade pressure side and high-pressure zones on the blade’s suction side. This, as confirmed by Figure 5b, is linked to the occurred flow separations on the blade pressure side around the blade leading edge under this condition. However, low-pressure zones and corresponding flow separations weakened in the downstream direction. The evolution of local pressure distribution mode, and associated flow structures as the machine influx decreased, are presented in Figure 6. In this figure, the selected plane is S2 which is, as already mentioned, located in the impeller’s most central zone (blade mid-length zone). As presented in Figure 6 left side, high-pressure zones were located on the blade pressure surface (BPS) in the vicinal flow zones of the impeller shroud. On the other hand, low-pressure zones were located on each blade’s suction side (BSS), in the shroud’s vicinal flow zones.
With the machine influx decrease from 1.00 QBEP to 0.55 QBEP, local pressure on the BSS was found to gradually increase, while the BPS pressure correspondingly decreased. The right side of Figure 6 presents the evolution of velocity-colored flow streams on the blade’s two surfaces (BPS and BSS). Keeping in mind that the plane in the left hand of Figure 6 was located in the impeller’s most central zone, the velocity-colored streams at the same zone showed high and comparatively low speed flows on the BSS and BPS, respectively. This explains the presence of low- and high-pressure zones on respective surfaces, as presented in Figure 6 left side. In addition, the level of flow unsteadiness was found to gradually worsen with the machine influx decrease. Under 0.82 QBEP conditions, there was an emergence of flow separations within the vicinal flow zones of the blade leading edge (BLE), especially on the BPS.
Under 0.74 QBEP conditions, the area occupied by BLE flow separations widened up, especially in the shroud vicinities, and extended towards downstream, eventually reaching the blade’s mid-length zone. In addition, another vortex flow zone appeared at the blade trailing edge (BTE) within the hub region and vicinities. Under this condition, the flow on the BSS also became slightly disturbed, especially in the hug region towards the BTE. Under 0.55 QBEP conditions, the BLE flow separation became even wider, while the BTE one almost did not grow. However, flow vortices on the BSS grew far larger, both at the BLE’s hub region and the BTE’s whole zone from hub to shroud, probably causing the blockage of the flow passage in the blade’s vicinal zone. It is therefore obvious that the impeller flow field becomes more and more disturbed with the decreasing machine influx, where flow separations in form of shock and vortical flows emerge and eventually grow. It was found that these secondary flow structures are more likely to take place at the impeller blade’s leading edge and flow zones in the vicinities of the impeller shroud, respectively.
As for the flow characteristics in the guide vanes region, Figure 7 shows the eventual changes in local flow structures as the machine influx decreased. The left side shows the velocity-colored streams at different flows. It can be seen that, as is the case with the impeller flow zone, the guide vanes flow pattern becomes messy when the flow rate decreases, and where vortical structures are emerge and eventually grow at certain flow zones, especially the vaneless space zone and the vane pressure side’s shroud vicinities. For high-flow conditions (1.00 QBEP and 0.82 QBEP), flow vortices in shroud vicinities occupy a tiny area (D), while flow wakes (E) formed at the vane’s trailing edge expand to the vaneless space.
Under 0.74 QBEP, there was a new tubular vortex that emerged in the vicinity of the guide vane trailing edge, expanding spanwise from hub to the shroud side (B). In addition, another vortex structure appeared near the guide vane’s outlet zone, within the guide vane shroud region (C). Under deep part-load condition (0.55 QBEP), type C vortex grew even larger, occupying the flow zones in the layer below, while type B grew bigger and slightly shifted its location towards upstream zones. In addition, a new and large vortical structure (A) developed from the guide vane’s leading edge on its suction side, and travelled downstream obliquely, finally touching the next vane’s pressure side in the vicinities of its trailing edge. The Type A vortex occurred mostly within the hub region.
It can therefore be briefly concluded that the guide vane flow pattern deteriorates with the decreasing flow where most of vortical structures form in the close proximity of the vaneless space.

3.2. Effect of Impeller Rotational Speed Variation

In line with this study’s main objectives, an effort was also made to investigate the effect of impeller rotational speed variation on the above discussed PAT flow dynamics under part-load conditions. In this respect, 0.74 QBEP was selected as the tested part-load flow condition where five impeller rotational speed ranging from 110 to 190 rpm were considered. For the sake of investigating the occurred changes in impeller flow field characteristics as its rotational speed increased, a mid-span plane (Span = 0.5) was taken as a reference flow zone.
Figure 8 displays the evolution of local pressure distribution within the impeller, as the impeller rotational speed gradually increased. It is generally found that low-pressure zones are located on BLE’s pressure side and BSS towards the trailing edge. With a continuously increasing impeller rotational speed, low-pressure zones grew increasingly larger at both mentioned zones, leading to a 190 rpm situation where they can then obviously extend to the draft tube flow zone (adjacent flow zone). In addition, high-pressure zones are found to be on the BLE’s suction side and continuously enlarge within the increase in impeller rotational speed.
In Figure 9, velocity-colored flow streamlines and correspondent turbulence kinetic energy distribution on the BPS are presented. Having considered a part-load condition (0.74 QBEP), the above-discussed flow separations were already present for the lowest of five tested speed values (110 rpm). Those were the BLE’s pressure side shock flow, its extension downstream in shroud vicinities, and the BTE’s hub vortex. With the increasing impeller rotational speed, the above-cited vortical structures strengthened and occupied even wide areas. This aspect is confirmed through Figure 9 right side where the evolution of BPS turbulence kinetic energy with the increasing speed is presented. It is shown that a high-BPS turbulence kinetic energy zone was located at the BLE close to the shroud under 110 rpm condition.
With the increasing rotational speed, this zone gradually widened up in both the spanwise and streamwise directions, reaching the blade mid-length zone under 190 rpm condition. It is therefore obvious that, for a PAT already operating under part-load conditions, the increase in impeller rotational speed worsens the level of flow instability.

3.3. Pressure Pulsation Characteristics

In order to study the evolution of pressure pulsation characteristics and its distribution mode for the above-discussed operating conditions, different flow zones were selected, at which pressure-monitoring points were positioned. To be more specific, nine pressure monitors were positioned on the mid-span plane of the inter-blade flow zone, as shown in Figure 10a. In addition, considering three consecutive span-wise planes from hub to the shroud, 18 monitors were positioned on the BPS, where each of the planes held six equidistant monitors expanding from the blade’s leading edge to the trailing edge as shown in Figure 10b.
Six monitors, BP16 to BP11, were located within the shroud’s vicinal zone. The other six, BP26 to BP21, were positioned on the mid-span plane, and the last six were located in the hub’s vicinal zone. Finally, 24 equidistant monitors were positioned around the whole vaneless space between the impeller and the guide vanes as shown in Figure 10c.
Figure 11 shows the pressure pulsation characteristics in the vaneless space zone for the four investigated flow conditions. The VS pressure pulsation spectra were found to be dominated by 7 fn component for both BEP and 0.82 QBEP conditions (fn: impeller rotational frequency).
Considering the impeller rotational speed of 150 rpm, this component represented the effect of rotor-stator interaction (RSI), known as the guide vane passing frequency (GPF = ZG × fn). Under BEP conditions, the VS pressure pulsation spectrum was marked by GPF as the 1st dominant frequency, followed by a low-frequency component (LFC) 0.17 fn. Moreover, more LFCs, which are multiples of fn, emerged between these two components (0.17 fn and GPF), adding to GPF harmonics which are considered high-frequency components. When the machine influx decreased from 1.00 QBEP to 0.82 QBEP, the above-described components persisted, but their amplitudes considerably increased, reaching to pulsation amplitude values four times higher.
Considerably remarkable changes in VS pressure pulsation characteristics were again recorded when the machine flow conditions reached 0.74 QBEP. Under this condition, unlike the first two conditions, the dominant frequency component shifted to LFC class, where 5.5 fn serves as the 1st dominant frequency component, followed by 0.64 fn as the 2nd dominant frequency. The RSI-born components are still present but with comparatively small amplitudes.
It is also worth mentioning that, compared to previous flow conditions, there was a global increase in pressure pulsation amplitudes, and a multitude of more LFCs emerged. This corresponds to increased local flow unsteadiness (LFU), as also discussed in the above sections. Under 0.55 QBEP condition, pulsation characteristics from the precedent flow condition persisted where, however, pulsation amplitudes continued to grow, and the 1st dominant frequency shifted to 4.32 fn, which is still classified as the LFC, while 0.66 fn served the 2nd dominant frequency.
In brief, the vaneless space pressure pulsation characteristics were found to be consistent with PAT flow dynamics evolution as the machine influx decreased. Under the first two high flow conditions, PAT flow experienced low-level of unsteadiness. Therefore, correspondent vaneless space pressure pulsation spectra were dominated by RSI-born frequency components.
Moreover, the number of LFCs was still low under these two conditions, which means a low level of local flow instability. On the other hand, the last two low flow conditions were marked by the dominance of LFCs and the increase in their numbers.
Though the RSI-born components were still present, their amplitudes were comparatively very low. This coincides with the increase in local flow unsteadiness under these conditions, where the impact of RSI on VS pressure pulsation characteristics became negligible, as compared to that inflicted by local flow instabilities.
In the same respect, an attempt was made to investigate the evolution of pressure pulsation levels in the vicinities of the impeller blade pressure surface for different flow conditions, as shown in Figure 12. This was performed using the relative pressure pulsation coefficient CP, for which the mathematical expression is presented in Equation (6).
C P = ( P i max P i min ) ρ g H
In this formula, Pimax, and Pimin, ρ, g, and H stand for the individual point’s maximum and minimum pressure pulsation amplitudes, fluid density, gravitational acceleration, and head, respectively. Figure 12 considers BPS monitoring points as shown in Figure 10b. In this figure, pressure pulsation amplitudes are found to decrease in the streamwise direction from the blade’s leading edge (vaneless space vicinities) towards its trailing edge (impeller outlet zone vicinities).
As also discussed in the above sections, it is also shown that the first two high flow conditions (1.00 QBEP and 0.82 QBEP) exhibited the lowest of pressure pulsation levels, while the last two (0.74 QBEP and 0.55 QBEP) were marked by high-pressure pulsation levels. Generally speaking, BPS pressure pulsation levels worsened with the decrease in machine flow conditions.
In the next section, an attempt was made to investigate the evolution in pressure pulsation distribution modes at different PAT flow zones, with the increasing impeller rotational speed, under two operating conditions, namely 1.00 QBEP and 0.74 QBEP. Figure 13 shows the pressure pulsation distribution mode within the vaneless space zone for different speeds, considering the above-mentioned flow conditions.
In this figure, as also explained in the above sections, the VS pressure pulsation level under optimum operating conditions were far lower compared to the 0.74 QBEP case. In addition, Figure 13 shows that the impeller rotational speed inflicts opposite reactions from both conditions’ VS pressure pulsation characteristics. It is shown that, under optimum operating conditions, the increase in impeller rotational speed leads to correspondingly decreasing vaneless pressure pulsation levels.
On the other hand, under 0.74 QBEP conditions, speed increases from 110 to 130 rpm led to the increase in pulsation levels, while the situation was complex for the rest of speeds. The pressure pulsation distribution mode for the rest of speeds (150, 170, and 190 rpm) was quite similar and their amplitudes quite randomly varied, but within short ranges.
Nevertheless, one would globally assume the worsening of pressure pulsation levels with the increase in impeller rotational speed under this condition, which would also be supported by correspondent flow field characteristics evolution, as presented in both Figure 8 and Figure 9.
In Figure 14, the effect of impeller rotational speed on pressure pulsation distribution mode along PAT’s full flow passage is displayed. In the direction of PAT water flow from the computational domain’s inlet zone towards the outlet (see Figure 1), different monitoring points were successively positioned along the full flow passage as follows: inlet pipe (OP4, OP33, OP32, and OP31), inter-guide vane channel (GV24, GV23, GV22, and GV21), vaneless space zone (VS), impeller inter-blade channel (IB32, IB22, and IB12), and outlet pipe (IP12, IP22, IB32, and IP4).
Under optimum operating conditions, pressure pulsation amplitudes within the guide vanes flow zone were globally found to be the highest, with the exception of 130 rpm case where the outlet pipe presented the highest of pressure pulsations. Under this condition, it is again confirmed that the increase in impeller rotational speed led to a global weakening of pressure pulsations along PAT’s full flow passage.
Under the 0.74 QBEP condition, however, the lowest of pressure pulsations were recorded for 110 rpm cases, followed by 130 rpm for the majority of flow zones. The situation is complex for the rest of speed cases, where, however, unlike the case optimum flow condition case, the highest pressure pulsation level was recorded within the vaneless space.
Therefore, one would globally conclude that, for a PAT under optimum operating conditions, the impeller rotational speed increase would improve the pressure pulsation levels along the machine at full flow passage, while causing an opposite effect for the case of part-load operating conditions.

4. Conclusions

This study numerically investigated the evolution of flow dynamics in a mixed-flow PAT under turbine operating mode, and the effect of impeller rotational speed on the same. Four different operating conditions ranging from optimum to part-load flow conditions (1.00 QBEP, 0.82 QBEP, 0.74 QBEP, and 0.55 QBEP) were selected, and five impeller rotational speeds ranging from 110 to 190 rpm were tested. After a thorough discussion and deep analysis of study results, a number of concluding remarks were drawn as follows:
With the machine influx decrease from optimum to 0.82 QBEP condition, PAT flow pattern slightly deteriorated. A further decrease to lower flow conditions inflicted the emergence of very large flow instabilities where, of the four investigated conditions, 0.55 QBEP presented the least stable flow structures. Moreover, the vaneless space zone between the impeller and guide vanes was found to experience the highest flow unsteadiness under part-load conditions.
Vaneless space pressure pulsation characteristics depended on the machine operating conditions, where pulsation amplitudes continuously increased as the machine influx decreased. For high flow conditions, vaneless space pressure pulsation spectrum was dominated by RSI-born frequency components, while the LFU-born frequency components dominated for low flow conditions.
The impeller rotational speed was found to differently impact the machine flow and pressure fields characteristics, depending on the considered operating condition. Under optimum operating conditions, the impeller rotational increase led to a global decrease in pressure pulsations at a majority of flow zones. On the other hand, the increase in impeller rotational speed worsened the pressure pulsation level at different flow zones, under part-load conditions.
It is therefore obvious that the impeller rotational speed exercises a considerable impact of the machine flow and pressure field characteristics, which under certain conditions may end up affecting the overall machine operational efficiency. This, for instance, sends an alert to electric regulation (ER) users within WDNs, where careful impeller rotational speed change is recommended. Further, CFD-supported parametric studies are still needed to reveal the flow dynamics of PATs under turbine mode, and to compare and establish a relationship with the pump mode. This would lead to a complete understanding of PAT flow dynamics, and probably contribute to an improvement in terms of PAT performance prediction accuracy, for the so far existent pump-to-turbine performance conversion expressions.

Author Contributions

Conceptualization, M.B., K.K. and H.-X.C.; Methodology, M.B., H.-X.C., W.-T.S. and X.-F.G.; Software, M.B., K.K., X.-F.G. and J.N.; Formal analysis, M.B., K.K., H.-X.C., Y.Z. and D.-Q.Z.; Writing-original draft, M.B., X.-F.G. and J.N.; Writing—review & editing, Y.Z., D.-Q.Z., H.-X.C. and K.K.; Resources, K.K., Y.Z. and D.-Q.Z.; Supervision, Y.Z., D.-Q.Z., and K.K.; Funding acquisition, M.B., K.K. and H.-X.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Postdoctoral Research Fund of Jiangsu Province (2020Z152), National Natural Science Foundation of China (52009033; 52006053), Natural Science Foundation of Jiangsu Province (BK20200509; BK20200508), and the Fundamental Research Funds for Central Universities (B210202066).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Acknowledgments

The support from the College of Water Conservancy and Hydropower Engineering at Hohai University, China is gratefully acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

Nomenclature

PATPump as Turbine
RANSReynolds Averaged Nevier-Stokes equations
BEPBest Efficiency Point
SSTShear Stress Transport turbulence model
IRSImpeller Rotational Speed
GGIGeneral Grid Interface
VSVaneless Space
LFCLow Frequency Component
CFDComputational Fluid Dynamics
TRLTechnology Readiness Level
HRHydraulic Regulation
LFULocal Flow Unsteadiness
GPFGuide vane Passing Frequency
RSIRotor-Stator Interactions
ERElectric Regulation
BLEBlade Leading Edge
BPSBlade Pressure Surface
BTEBlade Trailing Edge
BSSBlade Suction Surface
GVGuide Vane
WDNWater Distribution Networks
GSSGuide vane Suction Surface
IECInternational Electrotechnical Commission
GPSGuide vane Pressure Surface
IPInlet Pipe
OPOutlet Pipe
VOSVariable Operating Strategy
PRVPressure Relieve Valve
y+Dimensional height from the wall
fnRunner rotational frequency
QDischarge
nRunner rotational speed
HHead
gGravitational acceleration
ηEfficiency
PStatic pressure
ρDensity
ZRNumber of runner blades
ZGNumber of guide vanes

References

  1. Algieri, A.; Zema, D.A.; Nicotra, A.; Zimbone, S.M. Potential energy exploitation in collective irrigation systems using pumps as turbines: A case study in Calabria (Southern Italy). J. Clean. Prod. 2020, 257, 120538. [Google Scholar] [CrossRef]
  2. Kandi, A.; Moghimi, M.; Tahani, M.; Derakhshan, S. Optimization of pump selection for running as turbine and performance analysis within the regulation schemes. Energy 2021, 217, 119402. [Google Scholar] [CrossRef]
  3. Al-Shetwi, A.Q.; Hannan, M.A.; Jern, K.P.; Mansur, M.; Mahlia, T.M.I. Grid-connected renewable energy sources: Review of the recent integration requirements and control methods. J. Clean. Prod. 2020, 253, 119831. [Google Scholar] [CrossRef]
  4. Namahoro, J.P.; Wu, Q.; Xiao, H.; Zhou, N. The asymmetric nexus of renewable energy consumption and economic growth: New evidence from Rwanda. Renew. Energy 2021, 174, 336–346. [Google Scholar] [CrossRef]
  5. Namahoro, J.P.; Wu, Q.; Xiao, H.; Zhou, N. The Impact of Renewable Energy, Economic and Population Growth on CO2 Emissions in the East African Region: Evidence from Common Correlated Effect Means Group and Asymmetric Analysis. Energies 2021, 14, 312. [Google Scholar] [CrossRef]
  6. Hamududu, B.; Killingtveit, A. Assessing Climate Change Impacts on Global Hydropower. Energies 2012, 5, 305–322. [Google Scholar] [CrossRef] [Green Version]
  7. Hamududu, B.H.; Killingtveit, Å. Hydropower Production in Future Climate Scenarios: The Case for Kwanza River, Angola. Energies 2016, 9, 363. [Google Scholar] [CrossRef] [Green Version]
  8. Haidar, A.M.A.; Senan, M.F.M.; Noman, A.; Radman, T. Utilization of pico hydro generation in domestic and commercial loads. Renew. Sustain. Energy Rev. 2012, 16, 518–524. [Google Scholar] [CrossRef]
  9. Ismail, M.A.; Othman, A.K.; Islam, S.; Zen, H. End Suction Centrifugal Pump Operating in Turbine Mode for Microhydro Applications. Adv. Mech. Eng. 2014, 6, 139868. [Google Scholar] [CrossRef]
  10. Motwani, K.H.; Jain, S.V.; Patel, R.N. Cost Analysis of Pump as Turbine for Pico Hydropower Plants–A Case Study. Procedia Eng. 2013, 51, 721–726. [Google Scholar] [CrossRef] [Green Version]
  11. Ranjitkar, G.; Huang, J.; Tung, T. In Proceedings of the Application of Micro-Hydropower Technology for Remote Regions, EIC Climate Change Technology. Ottawa, ON, Canada, 10–12 May 2006; pp. 1–10. [Google Scholar]
  12. Renzi, M.; Rudolf, P.; Štefan, D.; Nigro, A.; Rossi, M. Installation of an axial Pump-as-Turbine (PaT) in a wastewater sewer of an oil refinery: A case study. Appl. Energy 2019, 250, 665–676. [Google Scholar] [CrossRef]
  13. Moazeni, F.; Khazaei, J. Optimal energy management of water-energy networks via optimal placement of pumps-as-turbines and demand response through water storage tanks. Appl. Energy 2021, 283, 116335. [Google Scholar] [CrossRef]
  14. Du, J.; Shen, Z.; Yang, H. Study on the effects of runner geometries on the performance of inline cross-flow turbine used in water pipelines. Sustain. Energy Technol. Assess. 2020, 40, 100762. [Google Scholar] [CrossRef]
  15. Du, J.; Shen, Z.; Yang, H. Effects of different block designs on the performance of inline cross-flow turbines in urban water mains. Appl. Energy 2018, 228, 97–107. [Google Scholar] [CrossRef]
  16. Jiyun, D.; Zhicheng, S.; Hongxing, Y. Numerical study on the impact of runner inlet arc angle on the performance of inline cross-flow turbine used in urban water mains. Energy 2018, 158, 228–237. [Google Scholar] [CrossRef]
  17. Jiyun, D.; Zhicheng, S.; Hongxing, Y. Study on the effects of blades outer angle on the performance of inline cross-flow turbines. Energy Procedia 2019, 158, 1039–1045. [Google Scholar] [CrossRef]
  18. Hosseini, S.M.H.; Forouzbakhsh, F.; Rahimpoor, M. Determination of the optimal installation capacity of small hydro-power plants through the use of technical, economic and reliability indices. Energy Policy 2005, 33, 1948–1956. [Google Scholar] [CrossRef] [Green Version]
  19. Mishra, S.; Singal, S.K.; Khatod, D.K. Optimal installation of small hydropower plant—A review. Renew. Sustain. Energy Rev. 2011, 15, 3862–3869. [Google Scholar] [CrossRef]
  20. Williams, A.A. Pumps as turbines for low cost micro hydro power. Renew. Energy 1996, 9, 1227–1234. [Google Scholar] [CrossRef]
  21. Ramos, H.; Borga, A. Pumps as turbines: An unconventional solution to energy production. Urban Water 1999, 1, 261–263. [Google Scholar] [CrossRef]
  22. Mercier, T.; Hardy, C.; Van Tichelen, P.; Olivier, M.; De Jaeger, E. Control of variable-speed pumps used as turbines for flexible grid-connected power generation. Electr. Power Syst. Res. 2019, 176, 105962. [Google Scholar] [CrossRef]
  23. Kusakana, K. A survey of innovative technologies increasing the viability of micro-hydropower as a cost effective rural electrification option in South Africa. Renew. Sustain. Energy Rev. 2014, 37, 370–379. [Google Scholar] [CrossRef]
  24. Jain, S.V.; Patel, R.N. Investigations on pump running in turbine mode: A review of the state-of-the-art. Renew. Sustain. Energy Rev. 2014, 30, 841–868. [Google Scholar] [CrossRef]
  25. Patelis, M.; Kanakoudis, V.; Gonelas, K. Pressure Management and Energy Recovery Capabilities Using PATs. Procedia Eng. 2016, 162, 503–510. [Google Scholar] [CrossRef]
  26. Adu, D.; Jianguo, D.; Darko, R.O.; Boamah, K.B.; Emmanuel, A. Investigating the state of renewable energy and concept of pump as turbine for energy generation development. Energy Rep. 2020, 6, 60–66. [Google Scholar] [CrossRef]
  27. Asomani, S.N.; Yuan, J.; Wang, L.; Appiah, D.; Zhang, F. Geometrical effects on performance and inner flow characteristics of a pump-as-turbine: A review. Adv. Mech. Eng. 2014, 12. [Google Scholar] [CrossRef]
  28. Kaunda, C.S.; Kimambo, C.Z.; Nielsen, T.K. A technical discussion on microhydropower technology and its turbines. Renew. Sustain. Energy Rev. 2014, 35, 445–459. [Google Scholar] [CrossRef]
  29. Štefan, D.; Rossi, M.; Hudec, M.; Rudolf, P.; Nigro, A.; Renzi, M. Study of the internal flow field in a pump-as-turbine (PaT): Numerical investigation, overall performance prediction model and velocity vector analysis. Renew. Energy 2020, 156, 158–172. [Google Scholar] [CrossRef]
  30. Lin, T.; Li, X.; Zhu, Z.; Xie, J.; Li, Y.; Yang, H. Application of enstrophy dissipation to analyze energy loss in a centrifugal pump as turbine. Renew. Energy 2021, 163, 41–55. [Google Scholar] [CrossRef]
  31. Hongyu, G.; Wei, J.; Yuchuan, W.; Hui, T.; Ting, L.; Diyi, C. Numerical simulation and experimental investigation on the influence of the clocking effect on the hydraulic performance of the centrifugal pump as turbine. Renew. Energy 2021, 168, 21–30. [Google Scholar] [CrossRef]
  32. Morabito, A.; Vagnoni, E.; Di Matteo, M.; Hendrick, P. Numerical investigation on the volute cutwater for pumps running in turbine mode. Renew. Energy 2021, 175, 807–824. [Google Scholar] [CrossRef]
  33. Ghorani, M.M.; Sotoude Haghighi, M.H.; Maleki, A.; Riasi, A. A numerical study on mechanisms of energy dissipation in a pump as turbine (PAT) using entropy generation theory. Renew. Energy 2020, 162, 1036–1053. [Google Scholar] [CrossRef]
  34. Ghorani, M.M.; Sotoude Haghighi, M.H.; Riasi, A. Entropy generation minimization of a pump running in reverse mode based on surrogate models and NSGA-II. Int. Commun. Heat Mass Transf. 2020, 118, 104898. [Google Scholar] [CrossRef]
  35. Wang, T.; Wang, C.; Kong, F.; Gou, Q.; Yang, S. Theoretical, experimental, and numerical study of special impeller used in turbine mode of centrifugal pump as turbine. Energy 2017, 130, 473–485. [Google Scholar] [CrossRef]
  36. Sengpanich, K.; Bohez, E.L.J.; Thongkruer, P.; Sakulphan, K. New mode to operate centrifugal pump as impulse turbine. Renew. Energy 2019, 140, 983–993. [Google Scholar] [CrossRef]
  37. Patelis, M.; Kanakoudis, V.; Gonelas, K. Combining pressure management and energy recovery benefits in a water distribution system installing PATs. J. Water Supply Res. Technol. Aqua 2017, 66, 520–527. [Google Scholar] [CrossRef] [Green Version]
  38. Cimorelli, L.; D’Aniello, A.; Cozzolino, L.; Pianese, D. Leakage reduction in WDNs through optimal setting of PATs with a derivative-free optimizer. J. Hydroinform. 2020, 22, 713–724. [Google Scholar] [CrossRef]
  39. Postacchini, M.; Darvini, G.; Finizio, F.; Pelagalli, L.; Soldini, L.; Di Giuseppe, E. Hydropower Generation Through Pump as Turbine: Experimental Study and Potential Application to Small-Scale WDN. Water 2020, 12, 958. [Google Scholar] [CrossRef] [Green Version]
  40. Pérez-Sánchez, M.; Sánchez-Romero, F.J.; Ramos, H.M.; López-Jiménez, P.A. Improved Planning of Energy Recovery in Water Systems Using a New Analytic Approach to PAT Performance Curves. Water 2020, 12, 468. [Google Scholar] [CrossRef] [Green Version]
  41. Stefanizzi, M.; Capurso, T.; Balacco, G.; Binetti, M.; Camporeale, S.M.; Torresi, M. Selection, control and techno-economic feasibility of Pumps as Turbines in Water Distribution Networks. Renew. Energy 2020, 162, 1292–1306. [Google Scholar] [CrossRef]
  42. Kougias, I.; Aggidis, G.; Avellan, F.; Deniz, S.; Lundin, U.; Moro, A.; Muntean, S.; Novara, D.; Pérez-Díaz, J.I.; Quaranta, E.; et al. Analysis of emerging technologies in the hydropower sector. Renew. Sustain. Energy Rev. 2019, 113, 109257. [Google Scholar] [CrossRef]
  43. Lima, G.M.; Luvizotto, E.; Brentan, B.M. Selection and location of Pumps as Turbines substituting pressure reducing valves. Renew. Energy 2017, 109, 392–405. [Google Scholar] [CrossRef]
  44. Polák, M. The Influence of Changing Hydropower Potential on Performance Parameters of Pumps in Turbine Mode. Energies 2019, 12, 2103. [Google Scholar] [CrossRef] [Green Version]
  45. Spedaletti, S.; Rossi, M.; Comodi, G.; Salvi, D.; Renzi, M. Energy recovery in gravity adduction pipelines of a water supply system (WSS) for urban areas using Pumps-as-Turbines (PaTs). Sustain. Energy Technol. Assess. 2021, 45, 101040. [Google Scholar]
  46. Carravetta, A.; Del Giudice, G.; Fecarotta, O.; Ramos, H.M. Pump as Turbine (PAT) Design in Water Distribution Network by System Effectiveness. Water 2013, 5, 1211–1225. [Google Scholar] [CrossRef] [Green Version]
  47. Carravetta, A.; Del Giudice, G.; Fecarotta, O.; Ramos, H.M. Energy Production in Water Distribution Networks: A PAT Design Strategy. Water Resour. Manag. 2012, 26, 3947–3959. [Google Scholar] [CrossRef]
  48. Carravetta, A.; del Giudice, G.; Fecarotta, O.; Ramos, H. PAT Design Strategy for Energy Recovery in Water Distribution Networks by Electrical Regulation. Energies 2013, 6, 411–424. [Google Scholar] [CrossRef] [Green Version]
  49. Fecarotta, O.; Carravetta, A.; Ramos, H.M.; Martino, R. An improved affinity model to enhance variable operating strategy for pumps used as turbines. J. Hydraul. Res. 2016, 54, 332–341. [Google Scholar] [CrossRef]
  50. Delgado, J.; Ferreira, J.P.; Covas, D.I.C.; Avellan, F. Variable speed operation of centrifugal pumps running as turbines. Experimental investigation. Renew. Energy 2019, 142, 437–450. [Google Scholar] [CrossRef]
  51. Lima Gustavo, M.; Luvizotto, E.; Brentan Bruno, M.; Ramos Helena, M. Leakage Control and Energy Recovery Using Variable Speed Pumps as Turbines. J. Water Resour. Plan. Manag. 2018, 144, 04017077. [Google Scholar] [CrossRef]
  52. Tahani, M.; Kandi, A.; Moghimi, M.; Houreh, S.D. Rotational speed variation assessment of centrifugal pump-as-turbine as an energy utilization device under water distribution network condition. Energy 2020, 213, 118502. [Google Scholar] [CrossRef]
  53. Jain, S.V.; Swarnkar, A.; Motwani, K.H.; Patel, R.N. Effects of impeller diameter and rotational speed on performance of pump running in turbine mode. Energy Convers. Manag. 2015, 89, 808–824. [Google Scholar] [CrossRef]
  54. Fernández, J.; Blanco, E.; Parrondo, J.; Stickland, M.T.; Scanlon, T.J. Performance of a centrifugal pump running in inverse mode. Proc. Inst. Mech. Eng. Part A J. Power Energy 2004, 218, 265–271. [Google Scholar] [CrossRef] [Green Version]
  55. Han, Y.; Tan, L. Influence of rotating speed on tip leakage vortex in a mixed flow pump as turbine at pump mode. Renew. Energy 2020, 162, 144–150. [Google Scholar] [CrossRef]
  56. Menter, F.R. Two-Equation Eddy-Viscosity Turbulence Models for Engineering Applications. AIAA J. 1994, 32, 1598–1605. [Google Scholar] [CrossRef] [Green Version]
  57. Wilcox, D.C. Reassessment of the scale-determining equation for advanced turbulence models. AIAA J. 1988, 26, 1299–1310. [Google Scholar] [CrossRef]
  58. Jones, W.P.; Launder, B.E. The prediction of laminarization with a two-equation model of turbulence. Int. J. Heat Mass Transf. 1972, 15, 301–314. [Google Scholar] [CrossRef]
  59. Menter, F. Zonal Two Equation k-w Turbulence Models For Aerodynamic Flows. In Proceedings of the 23rd Fluid Dynamics, Plasmadynamics, and Lasers Conference, American Institute of Aeronautics and Astronautics, Orlando, FL, USA, 6–9 July 1993. [Google Scholar]
Figure 1. Study roadmap.
Figure 1. Study roadmap.
Sustainability 13 07998 g001
Figure 2. PAT computational domain and components.
Figure 2. PAT computational domain and components.
Sustainability 13 07998 g002
Figure 3. PAT Experimental test rig and components (a) Inlet pipe (b) Shaft casing (c) Outlet pipe.
Figure 3. PAT Experimental test rig and components (a) Inlet pipe (b) Shaft casing (c) Outlet pipe.
Sustainability 13 07998 g003
Figure 4. Comparison of test and numerical results.
Figure 4. Comparison of test and numerical results.
Sustainability 13 07998 g004
Figure 5. Impeller flow field characteristics at different planes in the streamwise direction for 0.55 QBEP condition (a) pressure distribution contours (b) velocity curl-colored flow streamlines.
Figure 5. Impeller flow field characteristics at different planes in the streamwise direction for 0.55 QBEP condition (a) pressure distribution contours (b) velocity curl-colored flow streamlines.
Sustainability 13 07998 g005
Figure 6. Evolution of Impeller flow field characteristics (a) 1.00 QBEP (b) 0.82 QBEP (c) 0.74 QBEP and (d) 0.55 QBEP.
Figure 6. Evolution of Impeller flow field characteristics (a) 1.00 QBEP (b) 0.82 QBEP (c) 0.74 QBEP and (d) 0.55 QBEP.
Sustainability 13 07998 g006aSustainability 13 07998 g006b
Figure 7. Guide vane flow structures evolution for different operating conditions (a) 1 QBEP, (b) 0.82 QBEP, (c) 0.74 QBEP, and (d) 0.55 QBEP.
Figure 7. Guide vane flow structures evolution for different operating conditions (a) 1 QBEP, (b) 0.82 QBEP, (c) 0.74 QBEP, and (d) 0.55 QBEP.
Sustainability 13 07998 g007
Figure 8. Local pressure variation with the impeller rotating speed under 0.74 QBEP flow condition.
Figure 8. Local pressure variation with the impeller rotating speed under 0.74 QBEP flow condition.
Sustainability 13 07998 g008
Figure 9. Variation in flow patterns under 0.74 QBEP operating condition (a) BPS velocity-colored flow streamlines (b) BPS turbulence kinetic energy contours.
Figure 9. Variation in flow patterns under 0.74 QBEP operating condition (a) BPS velocity-colored flow streamlines (b) BPS turbulence kinetic energy contours.
Sustainability 13 07998 g009aSustainability 13 07998 g009b
Figure 10. Pressure monitoring point locations. (a) Inter-blade channel (b) Blade pressure side and (c) Vaneless space.
Figure 10. Pressure monitoring point locations. (a) Inter-blade channel (b) Blade pressure side and (c) Vaneless space.
Sustainability 13 07998 g010
Figure 11. Vaneless space pressure pulsation spectra under different flow conditions (a) 1.00 QBEP (b) 0.82 QBEP (c) 0.74 QBEP and (d) 0.55 QBEP.
Figure 11. Vaneless space pressure pulsation spectra under different flow conditions (a) 1.00 QBEP (b) 0.82 QBEP (c) 0.74 QBEP and (d) 0.55 QBEP.
Sustainability 13 07998 g011
Figure 12. Pressure pulsation extent in BPS vicinities at three span-wise planes under different flow conditions (a) Shroud vicinal zone (b) Mid-span zone and (c) Hub vicinal zone.
Figure 12. Pressure pulsation extent in BPS vicinities at three span-wise planes under different flow conditions (a) Shroud vicinal zone (b) Mid-span zone and (c) Hub vicinal zone.
Sustainability 13 07998 g012
Figure 13. Vaneless space pressure pulsation distribution mode for different impeller rotational speeds (a) 1.00 QBEP (b) 0.74 QBEP.
Figure 13. Vaneless space pressure pulsation distribution mode for different impeller rotational speeds (a) 1.00 QBEP (b) 0.74 QBEP.
Sustainability 13 07998 g013
Figure 14. Pressure pulsation distribution mode along PAT flow passage for different impeller rotational speeds (a) 1.00 QBEP (b) 0.74 QBEP.
Figure 14. Pressure pulsation distribution mode along PAT flow passage for different impeller rotational speeds (a) 1.00 QBEP (b) 0.74 QBEP.
Sustainability 13 07998 g014
Table 1. Hydraulic head variation with different grid numbers.
Table 1. Hydraulic head variation with different grid numbers.
Testing ParameterGrid No. 1 (3.9 M)Grid No. 2 (4.63 M)Grid No. 3 (5.84 M)Grid No. 4 (7.19 M)Grid No.5 (7.8 M)Grid No.6 (8.1 M)
Head (m)6.37.538.269.29.219.19
Table 2. Selected grid details for each PAT component.
Table 2. Selected grid details for each PAT component.
Grid DetailsInlet PipeImpellerGuide VanesOutlet Pipe
Grid typeHexahedralTetrahedralHexahedralHexahedral
Grid number (Million)1.22.512.171.31
Orthogonal Quality (0–1)0.450.320.40.4
Skewness (0–1)0.450.30.40.35
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Binama, M.; Kan, K.; Chen, H.-X.; Zheng, Y.; Zhou, D.-Q.; Su, W.-T.; Ge, X.-F.; Ndayizigiye, J. A Numerical Investigation into the PAT Hydrodynamic Response to Impeller Rotational Speed Variation. Sustainability 2021, 13, 7998. https://doi.org/10.3390/su13147998

AMA Style

Binama M, Kan K, Chen H-X, Zheng Y, Zhou D-Q, Su W-T, Ge X-F, Ndayizigiye J. A Numerical Investigation into the PAT Hydrodynamic Response to Impeller Rotational Speed Variation. Sustainability. 2021; 13(14):7998. https://doi.org/10.3390/su13147998

Chicago/Turabian Style

Binama, Maxime, Kan Kan, Hui-Xiang Chen, Yuan Zheng, Da-Qing Zhou, Wen-Tao Su, Xin-Feng Ge, and Janvier Ndayizigiye. 2021. "A Numerical Investigation into the PAT Hydrodynamic Response to Impeller Rotational Speed Variation" Sustainability 13, no. 14: 7998. https://doi.org/10.3390/su13147998

APA Style

Binama, M., Kan, K., Chen, H. -X., Zheng, Y., Zhou, D. -Q., Su, W. -T., Ge, X. -F., & Ndayizigiye, J. (2021). A Numerical Investigation into the PAT Hydrodynamic Response to Impeller Rotational Speed Variation. Sustainability, 13(14), 7998. https://doi.org/10.3390/su13147998

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop