Next Article in Journal
Extended k-Gamma and k-Beta Functions of Matrix Arguments
Previous Article in Journal
New Approach for a Weibull Distribution under the Progressive Type-II Censoring Scheme
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Improvement of Internal Flow Performance of a Centrifugal Pump-As-Turbine (PAT) by Impeller Geometric Optimization

1
National Research Center of Pumps, Jiangsu University, Zhenjiang 212013, China
2
Institute of Fluid Engineering Equipment, JITRI, Jiangsu University, Zhenjiang 212013, China
3
School of Chemistry, Physics and Mechanical Engineering, Queensland University of Technology, Brisbane 4001, QLD, Australia
*
Authors to whom correspondence should be addressed.
Mathematics 2020, 8(10), 1714; https://doi.org/10.3390/math8101714
Submission received: 5 September 2020 / Revised: 27 September 2020 / Accepted: 1 October 2020 / Published: 5 October 2020

Abstract

:
Rotor-stator interaction (RSI) in the centrifugal pump-as-turbine (PAT) is a significant source of high amplitude of the pressure pulsation and the flow-induced vibration, which is detrimental to the stable operation of PAT. It is therefore imperative to analyze the rotor-stator interaction, which can subsequently be used as a guideline for reducing the output of PAT noise, vibration and cavitation. In addition, it is important for a PAT to have a wide operating range preferably at maximum efficiency. In order to broaden the operating range, this work proposes a multi-condition optimization scheme based on numerical simulations to improve the performance of a centrifugal PAT. In this paper, the optimization of PAT impeller design variables (b2, β1, β2 and z) was investigated to shed light upon its influence on the output efficiency and its internal flow characteristics. Thus, the aim of the study is to examine the unsteady pressure pulsation distributions within the PAT flow zones as a result of the impeller geometric optimization. The numerical results of the baseline model are validated by the experimental test for numerical accuracy of the PAT. The optimized efficiencies based on three operating conditions (1.0Qd, 1.2Qd, and 1.4Qd) were maximally increased by 13.1%, 8.67% and 10.62%, respectively. The numerical results show that for the distribution of PAT pressure pulsations, the RSI is the main controlling factor where the dominant frequencies were the blade passing frequency (BPF) and its harmonics. In addition, among the three selected optimum cases, the optimized case C model exhibited the highest level of pressure pulsation amplitudes, while optimized case B reported the lowest level of pressure pulsation.

1. Introduction

The need for electrical power today is on the rise, owing to the population and economic growth. Life without electricity today is practically impossible and unpleasant because it has now become a part of our everyday life. However, getting electricity to remote areas poses a problem due to the topographical nature of these areas. Yet, most remote areas have streams and rivers with high and low head that can be difficult to generate electricity at a cheaper rate. To curb this problem, the ideal method is to apply the pump-as-turbine (PAT) technology which has the tendency to generate electricity at a cheaper rate and be more environmentally friendly. In contrast to other types of small hydraulic devices, PAT has the advantages of being small, relatively cheap and easy to produce and maintain. Due to it numerous advantages, the PAT has been widely used in many fields of engineering, including oil and chemical industries, small- or micro-hydropower plants, et cetera. [1]. Although the use of PAT is associated with substantial benefits, there are several significant disadvantages that should not be ignored. A major drawback of PAT is the lack of guide vanes and flow control devices which limit the operating range and the ability to control the operating point efficiently [2,3]. A wide range of operations is of great significance for a PAT to perform at optimum efficiency by saving energy [4,5]. Another issue relating to PAT is the difficulty to predict the actual performance of PAT mode performance [6,7]. The efficiency of PAT compared with pump mode operation is also significantly lower at off-design conditions [8].
The current PAT research is mainly directed towards developing precise prediction of the performance of turbine mode [9] and optimizing the design [10]. Very few researches have focused on the unsteady pressure field analysis of PAT [11]. The transmission of fluid machinery is unsteady due to the impulses of the rotating blades with respect to the stator [12]. The interaction between rotor and stator produces strong pressure pulsations which affects not only the hydraulic performance but also causes excessive noise and vibration [13,14]. This in turn jeopardizes the longevity of the components. Hence, to improve the overall performance and reliability, it is very important to investigate and analyze the mechanism of the pressure pulsations in PAT.
Several research studies have investigated the unstable flow characteristics in hydraulic pumps, and the following are among them. Pei et al. [15] numerically researched a radial pump and predicted the pressure pulsations field within the entire flow channel. Numerical investigations were carried out by Spence and Aaral-Teixeira concerning the effect of four geometrical parameters on pressure pulsation in the centrifugal pump: blade tip clearance, vane configuration, hemlock radial gap and sidewall clearance [16]. As a result, increasing the distance between the rotor and the stator was found to be an effective and reasonable means to reduce the amplitude of pressure pulsation. In the pump mode, the pump-turbine flow characteristics were examined by Yao et al. [17]. The results showed that in the tube draft, runner and diffuser complex flows, such as backflow or vortex, lead to the large pulsations of pressure. Aiming at a double suction centrifugal pump, Yao et al. [18] again studied time-frequency characteristics of pressure pulsation using a fast Fourier transform (FFT) algorithm and time-frequency methods. It was found that under off-design conditions, some unsteady flow phenomena occurring inside pumps have considerable effect on the pressure pulsation. Furthermore, the authors investigated the pressure pulsations caused by rotor-stator interactions in axially split centrifugal pumps under design and off-design operating conditions. They showed that the pulsating behavior of the blade passing-through frequency was most prominent near the volute tongue, and the pressure wave propagates in both radial and circumferential directions. Tan et al. [19] showed that the volute tongue region of a centrifugal pump is relatively uniform in the design operational condition of unsteady pressure and velocity pulsations. The effect of the volute tongue on the pulsations of PAT was investigated by Yang et al. [20]. The result indicated that the amplitude of the pressure pulsation decreases as the radial clearance increases. Barrio et al. [21] also performed a computational analysis using computational fluid dynamics (CFD) codes on the unstable PAT flow. The study found that the flow only matches the geometry of the impeller in reverse mode under design conditions, and recirculating fluid regions form at part load (close to the blade discharge) and at overload (near the suction side). From the above literature, it can be established that most studies through numerical methods on the pressure pulsations were carried out in pump mode with very few in PAT mode.
As pumps are not designed to operate as turbines, their geometrical parameters are determined without taking the PAT mode operation into account [22]. Therefore, to improve the operational reliability and efficiency, the hydraulic design of the PAT has to be optimized to suit both operational modes. The traditional hydraulic design of the PAT however deeply depends on the designers’ experience, which restricts its development [23]. As such, to improve the efficiency and operational reliability of the PAT, there is the need to apply advanced optimization algorithms so that optimization would not solely depend on the designers’ experience. Several researchers have proposed modifications to optimize the performance of the hydraulic characteristics. These modifications include trimming the impeller blades, adjusting the blade number, adding splitter blades, installing guide vanes and rounding impeller leading edges, et cetera. [24,25,26].
To save computational resources and obtain the optimal parameter combination accurately within a certain range, a surrogate model and an optimization algorithm must be applied. The application of optimization algorithm with the surrogate model (which can build a function between design variable and optimization objective) can improve pump performance efficiently. This method has been widely applied in centrifugal pumps optimization. Nourbakhsh et al. [27] used a particle swarm optimization method and non-dominated sorting genetic algorithm (NSGA) II algorithm with an artificial neural network (ANN) model to find the Pareto front of two conflict objectives of centrifugal pumps: efficiency and NPSHr. Derakhshan et al. [28] also used a global optimization method based on the ANN and the artificial bee colony algorithm along with a validated 3-D Navier–Stokes flow solver to redesign the impeller geometry and improve the performance of a pump.
Above all, the primary aim of the study is to optimize the geometrical parameters of the PAT impeller (b2, β1, β2 and z) at a wide range of flow conditions to examine its influence on hydraulic performance. Secondly, the aim is to examine, analyze and predict the unsteady pressure pulsation distributions within the PAT flow zones under different operating conditions. This study proposes an optimization procedure combining the design of experiment (DOE), surrogate models and Pareto-based genetic algorithm (PBGA) with numerical simulation to improve hydraulic efficiencies and operational reliability of a PAT under 1.0Qd, 1.2Qd and 1.4Qd (Qd is the design flow rate) flow conditions. The surrogate model is applied to construct a nonlinear function between the objective functions and the design variables at design and off-design conditions. Moreover, the simulation accuracy is verified by comparison with the experimental data for the baseline PAT model, and the flow field analysis of the original design and Pareto-optimal design are compared to illustrate the performance improvement of PAT.

2. Numerical and Experimental Setup

2.1. PAT Modelling

This study selects a baseline centrifugal pump which operates both in pump and turbine modes in National Research Centre of Pumps. The detailed design parameters of the pump are as follows: 25 m3/h flow rate; 8.5 m head; 1450 r/min rotating speed and 81.5 specific speed. Figure 1 shows the computational flow domain of the pump model working as a turbine. As can be seen, the flow enters from the volute, then drives the impeller to rotate and finally exits from the outlet pipe.

2.2. Computational Setup

In the study, ANSYS ICEM software is applied to create hexahedral structured grids for all components of the numerical model. The impeller component is periodically meshed which is set as the rotational domain. The hub, shroud, back and front chamber parts are meshed separately, while the volute and the discharge pipe are put together as one domain. The structural grids for the three domain parts are shown in Figure 2. Particularly, the meshes close to the interfaces between the impeller and the radial volute are densified to increase the complexity of the flow fields in that region due to its flow exchange characteristics. To consider the effect of the number of grids on the calculation accuracy, the grids were defined in the near-wall flow regions and the mesh independence check was performed with the grid numbers. The total number of grids was approximately 2.9 million [29]. After the analysis, the mesh of 2.9 × 106 grids is selected for the CFD simulation study of the baseline pump model since it reaches a good balance of the accuracy and computational cost.
For the computational setup of the simulation, the same setting as in reference [15] was preferred. For the turbulence model, the standard k-ω STT turbulence model as shown in Equations (1)–(4) is employed. The impeller operates at 1450 r/min, and 1000 iterations are calculated in the steady state. For the unsteady simulation, eight complete revolutions are computed while the time step of 1° rotational degree is set.
D D t ( ρ ω ) = ( ρ D ω ω ) + ρ γ G ν 2 3 ρ γ ω ( u ) ρ β ω 2 ρ ( F 1 1 ) C D k ω + S ω
D D t ( ρ k ) = ( ρ D k k ) + ρ G 2 3 ρ k ( u ) ρ β * ω k + S k
where, in the equation, the expression of turbulent viscosity coefficient ν t is expressed by:
ν t = a 1 k max ( a 1 ω b 1 F 23 S )
α k 1 = 0.85 , α k 2 = 1 , α ω 1 = 0.5 , α ω 1 = 0.856 , β 1 = 0.075 , β 2 = 0.0828 γ 1 = 5 / 9 , γ 2 = 0.44 , β * = 0.09 , a 1 = 0.31 , b 1 = 1 , c 1 = 10

2.3. Experimental Setup

A complete schematic of the experimental set up is shown in Figure 3. The experimental set up with pressure gauge, valves, underground tank and pipes was installed. The baseline pump model introduced above is used as the PAT. The flow rate is measured using an electromagnetic flow meter installed in the main pipe line. The inlet pressure is measured by an electronic pressure gauge and outlet pressure is measured by a static pressure sensor with the measuring range up to 2 MPa. To measure the output power, a synchronous generator is coupled to the PAT. The torque is acquired with a torque sensor and an auxiliary pump is also used to supply the water at desired head and flow rate. The auxiliary pump supply water is from the underground tank to the inlet of the pump tested. The real PAT test rig operation is shown in Figure 4. Due to the chaos of the real test rig background, the labels of different components in Figure 4 are not marked which can be easily identified by comparison with Figure 3.

3. Optimization Process

As shown in Figure 5, the optimization process is divided into four main steps. First, the optimization objective is selected, and the design variables and variables range are determined by a sensitivity test. Second, the sampling space is obtained by design of experiment (DOE). In this space, the objective values according to every design variable are calculated by ANSYS CFX software. Third, the surrogate model is constructed and the accuracy of the surrogate model is verified by error analysis. Finally, the Pareto-optimal solutions are found by running the optimization algorithm.

3.1. Optimization Objective

To broaden the operating range of the PAT efficiency area, the efficiencies based on design point, 1.0Qd, best efficiency point (BEP), 1.2Qd and overload point, 1.4Qd are selected as the optimization objectives. The PAT efficiency is defined as:
η = max ( P ρ g Q H )
where ρ and g are the flow density and acceleration of gravity, respectively. Q, H and P are the flow rate, head and mechanical power, respectively.

3.2. Design Variables and Surrogate Models

Due to the design experience and impeller geometrical parameters that have great impact on the efficiency and head [30], the blade outlet width b2, blade exit angles β1 and β2 on the hub and shroud and blade number z are chosen as the input variables, and their limits used as the optimization constraint are shown in Table 1 [31].
In order to establish the relationship between the input parameter and the output properties of PAT, the artificial neural network (ANN) is used for the surrogate modeling [32]. Figure 6 shows the general structure of the ANN and it has been adopted to fit the PAT efficiencies at the part load, BEP and overload flow conditions with four design variables of the impeller. The steps for the basic principles of the ANN were carried out in our previous study [31]. As shown in the Figure 6, the ANN includes input layer, hidden layer and output layer. The hidden layer and output layer are connected by a weight and activation function. As a result, the trained predictive model can accurately predict the PAT efficiency.

3.3. Multi-Objective Optimization

A multi-objective optimization problem is different from a single objective optimization problem because of the requirement to optimize several objectives. For multi-objective problems with several conflicting objectives, there is usually an optimal solution trade-off. The advantage of a PBGA is that the method for determining multiple objectives is reduced to a single fitness measure by the creation of a number of fronts that are then sorted according to non-domination. The mathematical steps of the multi-objective optimization calculations are done with reference to [31]. The detailed setting parameters of the PBGA are listed in Table 2. The population size and number of generations will directly affect the calculation cost of the optimization process and the quality of the optimization results. Therefore, in order to avoid inbreeding and poor convergence in the optimization process, the population number and the number of generations are determined. In addition, crossover fraction and function tolerance are set to 0.85 and 10, respectively, to ensure population diversity and better convergence.

4. Results and Discussion

4.1. Experimental Validation

The optimization of a PAT relies on the numerical simulation method which starts from the baseline PAT model. Hence, the accuracy of CFD simulation is critical to the reliability of the optimization results. In this section, the experimental test result of the baseline PAT model is presented before introducing the optimization outcomes. The hydraulic PAT characteristics curves of the CFD and the experimental test are plotted at different flow conditions. The numerical calculations in this research are carried out in line with the experimental test results for validation. Figure 7 shows characteristics of the efficiency and head curves. Here, the efficiencies of the PAT increase to reach the best efficiency point (BEP) at flow rate of 30 m3/h and begin to decline afterwards under both CFD and experimental measurements as established in Figure 7a. At the best efficiency point, the CFD and the experimental setup recorded hydraulic efficiency of 75.35% and 74.15%, respectively. The efficiencies of the CFD data compared with their respective experimental data increase about 1.6%, respectively. Thus, the CFD results record a slightly higher efficiency than the experiments and this can account for the flow leakages and other high turbulent effects that exist in the flow under these operating points. This can also be attributed to the mechanical losses ignored during the numerical calculations.
Figure 7b shows the pressure head characteristics curve. Generally, it can be observed that with increasing flow rate, the head of the CFD and experimental test increased. In turbine mode the head increased steadily as the flow rate increased from part-load to overload operating condition.

4.2. Analysis of Sensitivity Test

The sensitivity test analysis on design variables of the geometry is very important to optimization procedures. These tests are applied to establish the level of influence between the design variables and the objective function. In this study, Isight software is employed for the sensitivity analysis. Figure 8 shows the effect of the main geometry parameters of the impeller on PAT efficiency under flow conditions (1.0Qd, 1.2Qd and 1.4Qd). As shown in Figure 8, the selected geometric parameters which are the blade exit angles, blade inlet width and blade number (β2, β1, b2 and z) have minimal impact on PAT efficiency under design point 1.0Qd. However, the biggest impact on the PAT efficiency was established under overload flow condition 1.4Qd.
Figure 8 presents a correlation between the input and outputs values. The correlations values are ranked from zero (0) to positive one (+1). If the value is closer to (+1) it indicates a high correlation between the input and output, and if the value is closer to (0) it indicates a weak correlation. The negative and positive signs shows whether the inputs and the outputs variables are negatively or positively correlated. At the design operating point 1.0Qd, the design variable β2 has the highest correlation value of 0.52, followed by the blade width b2 which has a negative correlation value of −0.17. The blade number had the least correlation. At the best efficiency point (BEP), β2 recorded a positive correlation value of 0.62, followed by the blade number z which recorded a positive correlation value of 0.14. The blade width b2 had the least correlation value of 0.2. Similarly, in the overload operation 1.4Qd, the design variable β2 recorded the highest correlation value of 0.58, followed by the blade width b2 with a positive correlation value of 0.36. The blade angle β1 also recorded the least correlation value of 0.07 as established in Figure 8. It has been established that the four main geometry parameters of the impeller (β2, β1, b2 and z) plays an important role under 1.0Qd, 1.2Qd and 1.4Qd for PAT efficiency. Hence, these four geometry parameters were selected as design variables.
Figure 9 shows a bar graph that ranks the design variables (inputs) over the objective function (outputs). The significance of these graphs is to establish which of the design variables have the most influence on the objective function. The blue and red colors in Figure 10 indicate a positive and negative relationship between the input and the output variables, respectively. It can be seen that at design point 1.0Qd that the order of influence of the parameters on the PAT efficiency are, β2 > z > β1 > b2; at the BEP 1.2Qd the order is β2 > β1 > b2 > z; at the overload operating condition 1.4Qd the order is β2 > b2 > β1 > z. It implies that β2 has a dominant influence on the single objective of maximizing the PAT efficiency.

4.3. Analysis of Multi-Objective Optimization

The prediction functions between design variables and objective functions were built by using an ANN model. The Pareto-based genetic algorithm (PBGA) was used to obtain 1000 groups of Pareto-optimal solutions as shown in Figure 10. Three representative points A, B and C were selected to verify the accuracy of PBGA prediction. The comparison of the efficiencies for the original baseline model and the three optimized cases at the design point, best efficiency point (BEP) and overload flow conditions is presented in Table 3. Evidently, the optimized cases have a great improvement in the PAT efficiency compared to the baseline model.
Figure 11 presents a comparison of the PAT model and the optimized cases A, B and C at different flow conditions. Table 4 indicates the improvement of the efficiencies for different optimized PAT cases with respect to the baseline model. At design point, optimized case B performed better with an increase of 13.1% in efficiency, while optimized cases A and C increased by 10.1% and 3.34%, respectively. At the best efficiency point (BEP), optimized case B increases by 8.67% in efficiency as compared to cases A and C which increase by 7.56% and 0.76%, respectively. At overload condition, optimized case B’s increase in efficiency is 10.62% while the optimized cases A and C’s increases are 9.32% and 2.74%, respectively.

4.4. Analysis of Internal Flow Characteristics

To investigate how the impeller design variables (β2, β1, b2 and z) impact the PAT internal flow field, the unsteady CFD simulations of the optimized PAT models are carried out at the flow rate of 25 m3/h and rotational speed of 1450 r/min. The resultant optimized PAT flow characteristics are presented and analyzed.

4.4.1. Static Pressure Distribution

(a) Volute Flow Passage
For the present optimized cases A, B and C, the turbulent pressure distribution contours under different flow conditions in the PAT volute are presented in Figure 12. It can be observed that the static pressure within the flow domain declines from the volute inlet to the outlet of the PAT. A relatively high pressure region is observed at the inlet of the volute in all the optimized cases. As the flow progresses through the volute flow channel, the static pressure decreases due to an increase in velocity as a result of the gradual reduction of the cross-sectional area of the volute. The pressure distribution in the volute flow channel is not periodic due to the asymmetric shape of the volute and the interaction between the impeller and the volute.
The flow at the volute tongue region exhibits higher pressure gradients, for all three optimized impeller cases. In the regions close to the volute tongue under design point, BEP and the overload exhibited high pressure fluctuation intensity in all the optimized impeller cases A, B and C. It can also be observed that as the flow condition increases from design point to overload, the magnitude of the pressure fluctuation intensity increases in unison. It is worth noting that pressure distribution within the volute casing is dependent on flow rate. Under all the flow conditions, the BEP (1.2Qd) exhibits a slightly smoother flow in all the optimized cases.
However, by comparing pressure fluctuation intensity in the optimized cases A, B, and C; optimized case C, exhibited high-pressure fluctuation intensity under all flow conditions, indicating high operational instability. This can be attributed to poor parameter combination between the volute and the optimized impeller and rotor-stator matching.
(b) Impeller Flow Passage
The pressure fluctuation intensity distribution within the impeller flow zones are shown in Figure 13. The pressure fluctuation intensity shows a continuous decrease from the blade leading edge (BLE) to the blade trailing edge (BTE) of the impeller. It is observed that static pressure on the pressure side (PS) of the blade is more intense than the suction side (SS). Thus, the impeller is pushed to rotate by the pressure difference between the surfaces. The pressure difference from the pressure side to the suction side of the impeller blade decreases from the blade leading edge (BLE) to blade trailing edge (BTE).
The minimum value of the static pressure inside the impeller is located at the trailing edge of the blades on the suction side. The decreasing static pressure from the BLE to the BTE of the impeller is due to the continuous transmission of energy from the fluid to the impeller blades. The pressure fluctuation intensity within the optimized cases increased as the flow conditions increased from design point (1.0Qd) to overload point (1.4Qd). Comparatively, optimized case C generated the higher pressure fluctuation intensity indicating likelihood of high vibration and noise; this phenomenon can be attributed to bad parameter combination of the design variables or impeller-volute matching. However, optimized case A exhibited moderate pressure fluctuation intensity compared with optimized case B, indicating that the design variables for optimized case A have a high tendency of improving PAT performance.
(c) Outlet Pipe Flow Passage
Pressure fluctuation distribution contours at the outlet pipe of each optimized case have been presented in Figure 14. It can be observed that the pressure fluctuation increases radially and decreases axially from the outlet pipe entrance. This is mainly because the flow at the impeller outlet shows a higher tangential flow velocity component than the radial, which ends up retarding the axial flow towards the downstream flow zones. However, as the flow gets away from the impeller outlet region, the impeller rotational speed influence on the flow field decreases. These occurrences weaken the tangential flow velocity and enhance the axial flow velocity. Optimized cases A and B exhibited very low pressure fluctuations at the outlet pipe as compared to optimized case C.

4.4.2. Time Domain History of Pressure Fluctuation Coefficient

To investigate the pressure fluctuation characteristics within the PAT flow zones quantitatively, five pressure monitoring points were placed in the volute and the impeller flow channel. As shown in Figure 15, two pressure monitoring points were placed in the outlet pipe. The recorded pressure fluctuation data were processed using fast Fourier transform (FFT). The pressure fluctuation coefficient is plotted at the design point, BEP and overload flow conditions at different monitoring points. For time dependent history, a non-dimensional coefficient, Cp, is employed with the expression below.
C p = ( P m o n P r e f ) 0.5 ρ U 2 2
where, the referred pressure P r e f is the average static pressure for a complete rotation period of the impeller and P m o n is the pressure at the monitoring points. In the figure, “V” denotes the volute cross-sectional area, which reduces from “V1” to “V5” as depicted in Figure 15. “Imp” denotes the impeller low zone. “Imp1”, “Imp2” and “Imp3” cover the impeller inlet width. “Imp4” is situated at the mid-section of the impeller while “Imp5” monitors the impeller outlet zone. “D” monitors the outlet pipe flow channel.
(a) Volute Flow Passage
Figure 16 demonstrates the time domain history of pressure fluctuation coefficient, Cp, in the volute flow channel for all the optimized impeller cases A, B and C. The time-domain curves represent the periodic variations in pressure coefficient amplitude harmonics with respect to time step. It should be noted that the time step here is used to represent impeller position. During the impeller rotation of 360°, eight peaks and eight valleys were seen, which represent the eight blades of the impeller passing within a shaft revolution. When the impeller rotates, the blades sweep over the monitoring points, and they alternate between the high and low pressure zone.
In the volute flow passage, the pressure fluctuation at monitoring points “V1”–“V5” had quite similar trends in all the optimized cases A, B and C. Within the volute flow passage, monitoring points “V1” and “V5” recorded the highest amplitude. These high pressure fluctuation occurrences may be due to the geometrical position of the monitoring points, which are situated around the volute tongue region. It has been observed that as the flow condition increase from design point (1.0Qd), best efficiency point (BEP) (1.2Qd) and to overload point (1.4Qd), the peak-to-peak pressure fluctuation amplitude appears to increase along.
Optimized case C recorded a higher amplitude of Cp = 15 and above, followed by optimized cases A and B with amplitude values of Cp = 13 and Cp = 12 under flow condition of 1.4Qd, respectively. Comparatively, monitoring points “V4” and “V3” recorded the weakest amplitude in all the optimized cases with the exception of optimized case C. Above all, the pressure fluctuation in the optimized case C is the strongest. Thus, optimized case C has the tendency to induce vibration and noise.
(b) Impeller Flow Passage
The time-domain pressure fluctuation for one impeller revolution at the five monitoring points (Imp1, Imp2, Imp3, Imp4 and Imp5) within the impeller flow channel is depicted in Figure 17. During impeller rotation of 360°, eight peaks and eight valleys are seen, which represent the eight blades of the impeller passing within a shaft revolution. This confirms that the pressure fluctuation intensity in the impeller flow channel depends on the interaction between the rotating impeller and the stationary volute. The average value of the pressure fluctuation intensity decreases rapidly along the impeller flow channel from the BLE to the BTL. This occurrence is due to the energy transfer from the high pressure fluid to the mechanical shaft.
For all the optimized cases A, B and C, the pressure fluctuation at monitoring point “Imp3” recorded the most intensive pressure fluctuation among the five monitoring points located within the impeller flow zone. The high-pressure fluctuation of these regions may be due to the geometrical position of the monitoring points. Thus, at the volute outlet, where the flow still possesses a high percentage of its momentum hits the blade leading edge (BLE) surface, this impact results in different complex flow structure formation.
Optimized case C recorded the highest amplitude of Cp = 15–16 with a flow condition of 1.4Qd. Optimized cases A and B had a quite similar trend, with optimized case A, recording the least amplitude of Cp = 11.2. A similar trend was also observed at monitoring point “Imp2” which is located at the mid-section of the blade leading edge (BLE). Monitoring point “Imp5”, which is located at the impeller trailing edge, rendered the lowest pressure coefficients, which are lower than Cp = 8, especially at the design flow condition, in all the optimized cases. This makes the impeller eye highly susceptible to cavitation. It can be seen that the amplitude of pressure pulsation increases with increasing flow conditions. Among all the optimized cases, optimized case C exhibited the most intense pressure fluctuation compared to optimized cases A and B.
(c) Outlet Pipe Flow Passage
Figure 18 shows the pressure fluctuation in the outlet pipe. In one complete revolution; the flow at “D1” and “D2” in all three optimized cases fluctuates eight times, which equals the blade number. This confirms that the pressure fluctuation in the outlet pipe is due to the interaction between the rotating impeller and the stationary volute. Pressure fluctuation distribution at monitoring points “D1” and “D2” are selected to study the pressure fluctuation intensity in the discharge pipe. The value of the pressure fluctuation increases from the center of the outlet pipe to its near wall as indicated in Figure 16.
The similarity of pressure fluctuation trends in all the investigated optimized cases A, B and C under all the flow conditions were expected to exhibit the lowest pressure fluctuation as compared to the volute and impeller monitoring points. It is evidential that pressure fluctuation reduces as the flow moves away from the impeller eye region as depicted by the two monitoring points “D1” and “D2”. At monitoring point “D1”, optimized case C recorded the lowest pressure fluctuation value of Cp = 5.8, under all the flow conditions. However, optimized cases A and B exhibited a similar trend with Cp = 6.1 under all the flow conditions. Optimized cases A, B and C at monitoring point “D2” exhibited a similar trend of low-pressure fluctuation under all the flow conditions.

4.4.3. Frequency Domain History

(a) Volute Flow Passage
Figure 19 shows fast Fourier transform (FFT)-based frequency domain pressure pulsations at points “V1”–“V5” for the three optimized cases. Note that the impeller rotational speed is 1450 r/min, which indicates that the impeller rotating frequency fn is 24.167 Hz, leading to the blade passing frequency (BPF) of 193.33 Hz as there are eight impeller blades. The frequencies excited by pressure fluctuation happened at 193.33 Hz (8 × fn), 386.6 Hz (16 × fn) and 773.3 Hz (32 × fn).
For the three optimized cases, the blade passing frequency (BPF) was generally found to be the main frequency. The pressure fluctuation amplitudes continuously increase from monitoring point “V3” all the way to “V1”, with the exception of “V2”, in all the optimized impeller cases. Monitoring point “V5” recorded the highest pressure coefficient amplitude on the blade passing frequency harmonics. It recorded a Cp value of about 0.96 with optimized case A under 1.4Qd flow condition, followed by optimized case B with a Cp value of 0.84 under flow condition of 1.4Qd. However, the monitoring point “V2” and “V3” recorded the weakest pressure coefficient amplitude. At monitoring point “V2” of optimized case C, a Cp value of 0.18 was attained under flow condition of 1.4Qd, followed by optimized case B, with a Cp value of 0.17 under flow condition of 1.4Qd.
(b) Impeller Flow Passage
Figure 20 shows the properties in the impeller flow passage. The pressure fluctuation amplitudes continue to decrease from the blade leading edge (BLE) to the blade trailing edge (BTE) as depicted in the figure. The monitoring point “Imp4” recorded the peak pressure coefficient amplitude of the blade passing frequency harmonics. At monitoring point “Imp4”, optimized case B, recorded a high Cp value of 0.5 under 1.4Qd flow condition, followed by optimized cases A and B underflow condition of 1.4Qd and 1.2Qd, respectively, recording a Cp value of 0.42.
Monitoring point “Imp1”, denotes the suction side (SS) of the impeller blade. Optimized case B recorded a high Cp value of 0.4, followed by optimized case A, with a Cp value of 0.27 under flow condition of 1.0Qd, respectively. Monitoring point, “Imp2” denoting the mid-section of the blade inlet width (b2) recorded a Cp value of 0.43 with optimized case A, under 1.2Qd flow condition, followed by optimized case B, with a Cp value of 0.38 under 1.4Qd flow condition. Monitoring point “Imp3”, representing the pressure side (PS) of the impeller blade, recorded a high Cp value of 0.44 with optimized case A under 1.2Qd flow condition, followed by optimized case B, under flow condition of 1.4Qd recorded a Cp value of 0.4. Monitoring point “Imp5”, denoting the impeller eye, recorded a comparatively lower Cp value of 0.16 with optimized case C under flow condition of 1.0Qd and 1.4Qd, respectively.
(c) Outlet Pipe Flow Passage
Figure 21 shows the properties in the monitoring points of the flow passage of the outlet pipe “D1” and “D2”. The pressure fluctuation decreases as the flow moves away from the impeller eye region. At monitoring point “D1”, optimized case C recorded high Cp values of 0.16 and 0.14 underflow condition 1.4Qd and 1.0Qd, respectively. Similarly, at monitoring point “D2”, optimized case C recorded Cp values of 0.11 and 0.10 under flow condition 1.4Qd and 1.2Qd. It can also be noticed that for optimized case C, outlet pipe pressure pulsation was the highest, whereas optimized cases A and B exhibited the lowest pressure pulsation amplitudes. It is established that the pressure fluctuation within the outlet pipe is the spread of the pressure fluctuation caused by the rotor-stator interaction within the volute and the impeller.

5. Conclusions

The selected centrifugal PAT is distinguished by the problem of low operating range and high pressure pulsations. Thus, it remains a major challenge to expand the operating range and reduce the pressure pulsation during operation of the PAT. In order to maximize the PAT efficiency under multi working conditions, the design variables of the PAT impeller blade have been optimized in this paper to broaden its operational range. Furthermore, the intensity of the pressure pulsation was examined with numerical methods throughout the entire PAT flow passage. Below are some of the key findings and recommendations that can be taken into account during the PAT design optimization to boost PAT performance and operating reliability.
  • In this paper, the three-objective optimization problem is successfully solved and 3D Pareto solutions are achieved. There is a big performance improvement in the internal flow characteristic for all three optimized cases. Among which, case B exhibits the best performance under all flow conditions. At the design point (1.0Qd), best efficiency point (1.2Qd) and overload point (1.4Qd) flow conditions, the efficiency of case B increases by 13.1%, 8.67% and 10.62%, respectively. This implies that larger blade exit angles on the hub (β2) and smaller blade exit angles on the shroud (β1) have a combinatorial effect on PAT efficiency improvement.
  • The analysis of the time and frequency domain of the unstable pressure distribution within PAT shows that the high frequency pulsation is reduced by design variables of optimized cases A and B. However, case C exhibits the highest pressure pulsations and the lowest efficiency, while case B presents the lowest level of pressure pulsations and the highest efficiency. Therefore, high frequency pulsation directly affects PAT noise and stability, resulting in increased energy loss and reduced efficiency.
  • The pressure pulsation main frequencies are usually found at the impeller rotation frequency and also at different monitoring points within the entire PAT flow passage. The rotor-stator interaction (RSI) therefore constitutes the main factor affecting the characteristics of the pressure pulsations within the PAT flow passage. In the three PAT components (volute, impeller and outlet pipes), the blade passing frequency (BPF) is found to be the main frequency.
  • An investigation of the unsteady pressure field indicates that the unsteady pressure pulsation could spread downstream. Therefore, the two forms of pressure pulsations are the impeller and the outlet pipe flow passage. The low-frequency pressure pulsation is also found to decrease rapidly along the PAT flow channel. Thus, the unstable pressure pulsation inside the outlet pipe results mainly from the distribution of the unsteady pressure within the volute.
  • The RSI produces the highest pressure pulsation distribution within the volute tongue region for the three optimized cases. Thus, volute tongue modification is a better approach to reducing the intensity of the pressure pulsations in the volute and the entire PAT.
  • For future research directions, it is recommended that more accurate and also computationally heavy CFD simulation models such as detached eddy simulation [33] should be applied to the numerical prediction of PAT performance. Moreover, it is worthwhile trying to convert the multi-objective optimization problem into a single-objective one [34], so that the optimization process can be simplified and hence increase the likelihood of reaching a global optimal solution.

Author Contributions

Methodology, J.X.; Writing—original draft, S.N.A.; Writing—review & editing, L.W., W.L. and R.L. All authors have read and agreed to the published version of the manuscript.

Funding

The authors would like to express their great gratitude for the financial support from the independent project of Institute of Fluid Engineering Equipment, Jiangsu University. The funding from the Natural Science Foundation of Jiangsu Province (Grant No. BK20180879), the High-level Talent Research Foundation of Jiangsu University (Grant No. 18JDG012) and the Construction of Dominant Disciplines in Colleges and Universities and Universities in Jiangsu (PAPD) are equally acknowledged.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Nautiyal, H.; Kumar, A. Reverse running pumps analytical, experimental and computational study: A review. Renew. Sustain. Energy Rev. 2010, 14, 2059–2067. [Google Scholar] [CrossRef]
  2. Luo, B.; Wang, C.; Xia, Y.; Ye, J.; Yang, X. Numerical simulation of flow-induced vibration of double-suction centrifugal pump as turbine. J. Drain. Irrig. Mach. Eng. 2019, 37, 313–318. [Google Scholar]
  3. Liu, Y.; Yang, S.; Kong, F.; Wang, T.; Chen, K. Numerical simulation of hydraulic turbines with forward-curved blades. J. Drain. Irrig. Mach. Eng. 2018, 36, 21–27. [Google Scholar]
  4. Wei, H.; Xuefeng, L.I.; Min, S.U.; Rennian, L.I.; Hao, C. Pressure fluctuation of solid-liquid flow in stator and rotor cascades of pump as turbine. J. Drain. Irrig. Mach. Eng. 2018, 36, 99–103. [Google Scholar]
  5. Kang, Y.; Liu, M.; Wang, X.; Tan, Y.; Liu, J. Numerical simulation of pressure distribution in bristle seal for turbomachinery. J. Drain. Irrig. Mach. Eng. 2018, 36, 420–425. [Google Scholar]
  6. Shi, G.; Liu, X.; Wei, W.; Liu, Y. Pressure pulsation characteristics in pump as hydraulic turbine with guide vanes. J. Drain. Irrig. Mach. Eng. 2017, 35, 6–12. [Google Scholar]
  7. Shi, F.; Yang, J.; Wang, X. Numerical analysis on transient characteristics of hydraulic turbine under variable working conditions. J. Drain. Irrig. Mach. Eng. 2017, 35, 200–206. [Google Scholar]
  8. Williams, A. The turbine performance of centrifugal pumps: A comparison of prediction methods. Proc. Inst. Mech. Eng. Part A J. Power Energy 1994, 208, 59–66. [Google Scholar] [CrossRef]
  9. Singh, P.; Nestmann, F. An optimization routine on a prediction and selection model for the turbine operation of centrifugal pumps. Exp. Therm. Fluid Sci. 2010, 34, 152–164. [Google Scholar] [CrossRef]
  10. Sun-Sheng, Y.; Fan-Yu, K.; Jian-Hui, F.; Ling, X. Numerical research on effects of splitter blades to the influence of pump as turbine. Int. J. Rotating Mach. 2012, 2012, 1–9. [Google Scholar] [CrossRef] [Green Version]
  11. Morros, C.S.; Oro, J.M.F.; Díaz, K.M.A. Numerical modelling and flow analysis of a centrifugal pump running as a turbine: Unsteady flow structures and its effects on the global performance. Int. J. Numer. Methods Fluids 2011, 65, 542–562. [Google Scholar] [CrossRef]
  12. Barrio, R.; Parrondo, J.; Blanco, E. Numerical analysis of the unsteady flow in the near-tongue region in a volute-type centrifugal pump for different operating points. Comput. Fluids 2010, 39, 859–870. [Google Scholar] [CrossRef]
  13. Spence, R.; Purdom, T. Numerical prediction of transient loadings on multistage pump impellers. Adv. CFD Fluid Mach. Des. 2003, 1–10. [Google Scholar]
  14. Yang, S.-S.; Liu, H.-L.; Kong, F.-Y.; Xia, B.; Tan, L.-W. Effects of the radial gap between impeller tips and volute tongue influencing the performance and pressure pulsations of pump as turbine. J. Fluids Eng. 2014, 136, 054501. [Google Scholar] [CrossRef]
  15. Pei, J.; Wang, W.; Yuan, S.; Zhang, J. Optimization on the impeller of a low-specific-speed centrifugal pump for hydraulic performance improvement. Chin. J. Mech. Eng. 2016, 29, 992–1002. [Google Scholar] [CrossRef]
  16. Spence, R.; Teixeira, J. A CFD parametric study of geometrical variations on the pressure pulsations and performance characteristics of a centrifugal pump. Comput. Fluids 2009, 38, 1243–1257. [Google Scholar] [CrossRef] [Green Version]
  17. Yao, Y.; Xiao, Y.X.; Zhu, W.; An, S.H.; Wang, Z.W. Numerical prediction of the pressure fluctuations on small discharge condition of a pump-turbine at pump mode. IOP Conf. Ser. Mater. Sci. Eng. 2015, 72, 032029. [Google Scholar] [CrossRef] [Green Version]
  18. Yao, Z.; Wang, F.; Qu, L.; Xiao, R.; He, C.; Wang, M. Experimental investigation of time-frequency characteristics of pressure fluctuations in a double-suction centrifugal pump. J. Fluids Eng. 2011, 133, 101303. [Google Scholar] [CrossRef]
  19. Tan, L.; Wang, Y.C.; Cao, S.L.; Zhu, B.S. Characteristics of unsteady flow around the tongue region in a centrifugal pump. Trans. Beijing Inst. Technol. 2014, 34, 670–675. [Google Scholar]
  20. Yang, S.; Kong, F.; Chen, B. Research on pump volute design method using CFD. Int. J. Rotating Mach. 2011, 2011, 1–7. [Google Scholar] [CrossRef]
  21. Barrio, R.; Fernández, J.; Parrondo, J.; Blanco, E. Performance prediction of a centrifugal pump working in direct and reverse mode using computational fluid dynamics. In Proceedings of the International Conference on Renewable Energies and Power Quality, Granada, Spain, 23–25 March 2010; pp. 23–25. [Google Scholar]
  22. 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]
  23. Shah, S.; Jain, S.V.; Patel, R.N.; Lakhera, V.J. CFD for centrifugal pumps: A review of the state-of-the-art. Procedia Eng. 2013, 51, 715–720. [Google Scholar] [CrossRef] [Green Version]
  24. Yang, S.-S.; Kong, F.-Y.; Qu, X.-Y.; Jiang, W.-M. Influence of blade number on the performance and pressure pulsations in a pump used as a turbine. J. Fluids Eng. 2012, 134, 124503. [Google Scholar] [CrossRef]
  25. Yang, S.-S.; Kong, F.-Y.; Jiang, W.-M.; Qu, X.-Y. Effects of impeller trimming influencing pump as turbine. Comput. Fluids 2012, 67, 72–78. [Google Scholar] [CrossRef]
  26. Patel, V.; Jain, S.; Motwani, K.; Patel, R.N. Numerical optimization of guide vanes and reducer in pump running in turbine mode. Procedia Eng. 2013, 51, 797–802. [Google Scholar] [CrossRef]
  27. Nourbakhsh, A.; Safikhani, H.; Derakhshan, S. The comparison of multi-objective particle swarm optimization and NSGA II algorithm: Applications in centrifugal pumps. Eng. Optim. 2011, 43, 1095–1113. [Google Scholar] [CrossRef]
  28. Derakhshan, S.; Pourmahdavi, M.; Abdolahnejad, E.; Reihani, A.; Ojaghi, A. Numerical shape optimization of a centrifugal pump impeller using artificial bee colony algorithm. Comput. Fluids 2013, 81, 145–151. [Google Scholar] [CrossRef]
  29. Asomani, S.N.; Yuan, J.; Wang, L.; Appiah, D.; Adu-Poku, K.A. The impact of surrogate models on the multi-objective optimization of Pump-As-Turbine (PAT). Energies 2020, 13, 2271. [Google Scholar] [CrossRef]
  30. Shouqi, Y. Low Specific Speed Centrifugal Pump Theory and Design [M]; Mechanical Industry Press: Beijing, China, 1997. [Google Scholar]
  31. McKay, M.D.; Beckman, R.J.; Conover, W.J. Comparison of three methods for selecting values of input variables in the analysis of output from a computer code. Technometrics 1979, 21, 239–245. [Google Scholar]
  32. McCulloch, W.; Pitts, W.J.M.B. A logical calculus of the ideas imminent in nervous activity. Bull. Math. Biophys. 1943, 5, 115–133. [Google Scholar] [CrossRef]
  33. Yuan, J.; Chen, Y.; Wang, L.; Fu, Y.; Zhou, Y.; Xu, J.; Lu, R. Dynamic analysis of cavitation tip vortex of pump-jet propeller based on DES. Appl. Sci. 2020, 10, 5998. [Google Scholar] [CrossRef]
  34. Wang, L.; Zuo, M.J.; Xu, J.; Zhou, Y.; Tan, A.C. Optimizing wind farm layout by addressing energy-variance trade-off: A single-objective optimization approach. Energy 2019, 189, 116149. [Google Scholar] [CrossRef]
Figure 1. Computational flow domain of the numerical model.
Figure 1. Computational flow domain of the numerical model.
Mathematics 08 01714 g001
Figure 2. Meshes for the different computational flow domains. (a) Hub and shroud, (b) impeller blades, and (c) volute with tongue enlargement.
Figure 2. Meshes for the different computational flow domains. (a) Hub and shroud, (b) impeller blades, and (c) volute with tongue enlargement.
Mathematics 08 01714 g002
Figure 3. Schematic layout of the pump-as-turbine (PAT) experimental setup: (a) motor, (b) pump inlet pipe, (c) pump outlet pipe, (d) bypass valve, (e) bypass pipe, (f) electromagnetic flow rate, (g) generator, (h) PAT, (i) inlet and outlet pressure transmit and (j) PAT outlet pipe.
Figure 3. Schematic layout of the pump-as-turbine (PAT) experimental setup: (a) motor, (b) pump inlet pipe, (c) pump outlet pipe, (d) bypass valve, (e) bypass pipe, (f) electromagnetic flow rate, (g) generator, (h) PAT, (i) inlet and outlet pressure transmit and (j) PAT outlet pipe.
Mathematics 08 01714 g003
Figure 4. Completed PAT setup for validating the baseline model.
Figure 4. Completed PAT setup for validating the baseline model.
Mathematics 08 01714 g004
Figure 5. Optimization process of PAT.
Figure 5. Optimization process of PAT.
Mathematics 08 01714 g005
Figure 6. Artificial neural network (ANN) structure for PAT efficiency objective.
Figure 6. Artificial neural network (ANN) structure for PAT efficiency objective.
Mathematics 08 01714 g006
Figure 7. Performance characteristic curve of the tested PAT and computational fluid dynamics (CFD): (a) efficiency (b) head.
Figure 7. Performance characteristic curve of the tested PAT and computational fluid dynamics (CFD): (a) efficiency (b) head.
Mathematics 08 01714 g007
Figure 8. Correlation of the input and output variables.
Figure 8. Correlation of the input and output variables.
Mathematics 08 01714 g008
Figure 9. Ranks of the design variables over the PAT efficiency at (a) 1.0Qd, (b) 1.2Qd and (c) 1.4Qd.
Figure 9. Ranks of the design variables over the PAT efficiency at (a) 1.0Qd, (b) 1.2Qd and (c) 1.4Qd.
Mathematics 08 01714 g009
Figure 10. Pareto-frontiers from ANN.
Figure 10. Pareto-frontiers from ANN.
Mathematics 08 01714 g010
Figure 11. Comparison of PAT efficiencies for different case scenarios.
Figure 11. Comparison of PAT efficiencies for different case scenarios.
Mathematics 08 01714 g011
Figure 12. Pressure distribution within the volute of the optimized cases A, B and C under flow conditions of: (a) 1.0Qd (b) 1.2Qd and (c) 1.4Qd.
Figure 12. Pressure distribution within the volute of the optimized cases A, B and C under flow conditions of: (a) 1.0Qd (b) 1.2Qd and (c) 1.4Qd.
Mathematics 08 01714 g012
Figure 13. Pressure distribution within the impeller of optimized cases A, B and C under flow conditions of: (a) 1.0Qd (b) 1.2Qd and (c) 1.4Qd.
Figure 13. Pressure distribution within the impeller of optimized cases A, B and C under flow conditions of: (a) 1.0Qd (b) 1.2Qd and (c) 1.4Qd.
Mathematics 08 01714 g013
Figure 14. Pressure distribution within the outlet pipe of optimized cases A, B and C under flow conditions of: (a) 1.0Qd (b) 1.2Qd and (c) 1.4Qd.
Figure 14. Pressure distribution within the outlet pipe of optimized cases A, B and C under flow conditions of: (a) 1.0Qd (b) 1.2Qd and (c) 1.4Qd.
Mathematics 08 01714 g014
Figure 15. Monitoring points in the impeller, volute and outlet pipe channel.
Figure 15. Monitoring points in the impeller, volute and outlet pipe channel.
Mathematics 08 01714 g015
Figure 16. Time-domain of pressure fluctuations at different monitoring points within the volute of the three optimized impeller cases at the locations of: (a) V1, (b) V2, (c) V3, (d) V4 and (e) V5.
Figure 16. Time-domain of pressure fluctuations at different monitoring points within the volute of the three optimized impeller cases at the locations of: (a) V1, (b) V2, (c) V3, (d) V4 and (e) V5.
Mathematics 08 01714 g016
Figure 17. Time-domain of pressure fluctuations at different monitoring points within the impeller of the three optimized impeller cases at the locations of: (a) Imp1, (b) Imp2, (c) Imp3, (d) Imp4 and (e) Imp5.
Figure 17. Time-domain of pressure fluctuations at different monitoring points within the impeller of the three optimized impeller cases at the locations of: (a) Imp1, (b) Imp2, (c) Imp3, (d) Imp4 and (e) Imp5.
Mathematics 08 01714 g017
Figure 18. Time-domain of pressure fluctuations at different monitoring points within the outlet pipe of the three optimized impeller cases at the locations of: (a) D1 and (b) D2.
Figure 18. Time-domain of pressure fluctuations at different monitoring points within the outlet pipe of the three optimized impeller cases at the locations of: (a) D1 and (b) D2.
Mathematics 08 01714 g018
Figure 19. Frequency domain history of pressure fluctuation in volute flow passage at the locations of: (a) V1, (b) V2, (c) V3, (d) V4 and (e) V5.
Figure 19. Frequency domain history of pressure fluctuation in volute flow passage at the locations of: (a) V1, (b) V2, (c) V3, (d) V4 and (e) V5.
Mathematics 08 01714 g019
Figure 20. Frequency domain history of pressure fluctuation in impeller flow passage at the locations of: (a) Imp1, (b) Imp2, (c) Imp3, (d) Imp4 and (e) Imp5.
Figure 20. Frequency domain history of pressure fluctuation in impeller flow passage at the locations of: (a) Imp1, (b) Imp2, (c) Imp3, (d) Imp4 and (e) Imp5.
Mathematics 08 01714 g020
Figure 21. Frequency domain history of pressure fluctuation in outlet pipe flow passage at the locations of: (a) D1 and (b) D2.
Figure 21. Frequency domain history of pressure fluctuation in outlet pipe flow passage at the locations of: (a) D1 and (b) D2.
Mathematics 08 01714 g021
Table 1. Ranges of design parameters.
Table 1. Ranges of design parameters.
Design VariablesLowerOriginalUpper
b2/mm0.0100.0120.015
β2203035
β1203035
z468
Table 2. Setting parameters of the Pareto-based genetic algorithm (PBGA).
Table 2. Setting parameters of the Pareto-based genetic algorithm (PBGA).
Design VariablesLower
Population size100
Number of generations1000
Pareto-front population0.8
Crossover fraction0.85
Function tolerance10
Table 3. Comparison of baseline model and the optimized cases.
Table 3. Comparison of baseline model and the optimized cases.
Casesb2beta2beta1z1.0Qd1.2Qd1.4Qd
Baseline0.0123030673.918478.845478.5914
A0.01433.5223.49881.373484.811685.9235
B0.01323.0628.09883.581585.695786.9401
C0.01323.8829.26876.389679.448480.741
Table 4. Improvement of optimized cases with respect to the baseline model.
Table 4. Improvement of optimized cases with respect to the baseline model.
CasesDesign Point (1.0Qd)BEP (1.2Qd)Overload Point (1.4Qd)
A10.1%7.56%9.32%
B13.1%8.67%10.62%
C3.34%0.76%2.74%

Share and Cite

MDPI and ACS Style

Xu, J.; Wang, L.; Ntiri Asomani, S.; Luo, W.; Lu, R. Improvement of Internal Flow Performance of a Centrifugal Pump-As-Turbine (PAT) by Impeller Geometric Optimization. Mathematics 2020, 8, 1714. https://doi.org/10.3390/math8101714

AMA Style

Xu J, Wang L, Ntiri Asomani S, Luo W, Lu R. Improvement of Internal Flow Performance of a Centrifugal Pump-As-Turbine (PAT) by Impeller Geometric Optimization. Mathematics. 2020; 8(10):1714. https://doi.org/10.3390/math8101714

Chicago/Turabian Style

Xu, Jian, Longyan Wang, Stephen Ntiri Asomani, Wei Luo, and Rong Lu. 2020. "Improvement of Internal Flow Performance of a Centrifugal Pump-As-Turbine (PAT) by Impeller Geometric Optimization" Mathematics 8, no. 10: 1714. https://doi.org/10.3390/math8101714

APA Style

Xu, J., Wang, L., Ntiri Asomani, S., Luo, W., & Lu, R. (2020). Improvement of Internal Flow Performance of a Centrifugal Pump-As-Turbine (PAT) by Impeller Geometric Optimization. Mathematics, 8(10), 1714. https://doi.org/10.3390/math8101714

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