Abstract

Centrifugal pumps as turbine (PAT) are widely used in petrochemical and water conservancy industries. The research on the internal flow field and energy loss of PAT is of great significance to improve the performance and efficiency of PAT. In this paper, experimental and numerical simulation methods are used to study the energy loss and flow field. The results show that the numerical simulation method can accurately simulate the internal flow field of PAT. And the entropy generation theory is applied to visualize the internal energy loss of PAT through the comparison of total pressure loss and entropy generation. The highest energy loss among PAT components is the guide vane. The loss in the guide vane is mainly caused by the flow separation caused by the wake of the guide vane and the asymmetric structure of the volute. The losses in the impeller are mainly due to flow separation and wake. Besides, the unsteady simulation results show that rotor-stator interaction has a great influence on the gap between the impeller and the guide vane. The research results provide a reference for the design of the PAT. This study is beneficial to studying the dynamic and static interference and PAT vibration to improve the stability of the PAT.

1. Introduction

Centrifugal pumps as turbines (PAT) are widely used in steel, petrochemical, and water conservancy industries to recover energy from high-pressure liquids in the production process due to their low cost and high efficiency. PAT is in a long-term operating state, so it is of great significance to study the performance and efficiency of PAT [14].

At present, the methods of experiment and numerical simulation are used by most scholars to study the performance of PAT [510]. Bozorgi et al. [11] used numerical simulation to numerically simulate and verify the reverse mode industrial axial pump. The results show that the numerical simulation method can accurately predict the PAT performance and confirm that the axial flow pump in the reverse mode can operate normally in a small hydropower station. Su et al. [12] used the CFD method to numerically simulate the internal flow field of PAT and revealed the periodic law of flow velocity distribution in PAT. Wang et al. [13] used ANSYS and NX software to design a special forward-curved blade impeller. The comparison between the two impellers shows that the new impeller can significantly improve the efficiency of the PAT. Liu and Tan [14] conducted numerical simulation and experimental research on PAT with tip clearance. The results show that the increase of the blade tip clearance will cause the drop of PAT head and efficiency and reveal the vortex characteristics of PAT. Liu et al. [15] used the SST k-ω turbulence model to numerically simulate PAT. A new method to determine the best efficiency point (BEP) of PAT was proposed, and it is found that the high-efficiency area of PAT is relatively wide, which exists in large flow conditions. Tang et al. [16, 17] used machine learning methods to conduct in-depth research on fault diagnosis of rotating machinery. Gao et al. [18] discussed the numerical characterization of the influence of flow on the pressure and velocity distribution of PAT. Analyzing the internal flow phenomenon of PAT under different operating conditions, the results show that the unsteady turbulent flow in PAT will cause the vortex to have a greater impact on the flow separation [19].

Recently, the research on the internal flow of hydraulic machinery from the perspective of energy has attracted more and more attention from scholars [2024]. Entropy generation theory has been applied in pipelines and pumps, but its application in PAT is also meaningful. In this paper, the entropy generation theory is used to analyze the energy loss characteristics in PAT. The SST k-ω turbulence model is used to simulate the unsteady result of PAT and verify it by experiment. Based on the numerical simulation results verified by experiments, the energy loss distribution and formation mechanism inside the PAT are analyzed, and rotor-stator interaction phenomena in the PAT are analyzed.

2. Entropy Production Theory

The internal flow structure of PAT is complicated, and it is difficult to obtain detailed internal energy loss distribution using conventional methods. The entropy generation theory can be used to analyze the internal energy loss in detail. Entropy is one of the parameters that characterize the state of matter in thermodynamics; the balance equation (1) is the entropy in a single-phase incompressible fluid [25]:where s represents the specific entropy, Φ is the viscous dissipation term; and ΦΘ is the dissipation term caused by the heat transfer. The entropy production analysis is based on the second law of thermodynamics, combined with heat transfer theory and fluid mechanics to calculate the actual process dissipation of energy loss. Among them, the two rightmost items are entropy production items.

The temperature is set to a constant temperature in this paper, the entropy production item on the right side is processed by the Reynolds time-averaged process, and the local entropy production term expression is obtained as the two following equations:where represents the entropy production by average velocity (W/m3) and is the entropy production by the pulsation velocity (W/m3). is the effective dynamic viscosity, which can be calculated according towhere μt is the turbulence dynamic viscosity (N·s/m2).

In the k-ε turbulent model, the local entropy production rate caused by the pulsation velocity is given in the following equation [26]:

The energy loss in the wall region (the area between the wall and the first layer of the grid) also occupies an important part. Since the entropy production in the wall could not be predicted precisely by these equations, the local entropy production rate for the wall is calculated using where is the wall shear stress (Pa) and is the velocity vector of the center of the first layer of the wall (m/s).

The total entropy produced by the mainstream area can be obtained by integrating the following equations:where is the total entropy produced by the time-averaged velocity in the main flow area (W) and is the total entropy production caused by the pulsation velocity in the main flow area (W).

So, the total entropy production of the PAT is shown in the following equation:

3. Numerical Model and Methods

3.1. PAT Model

To analyze the internal flow characteristics and energy loss mechanism of PAT, a numerical simulation of a low specific speed single-stage centrifugal pump was carried out. The main parameters of the centrifugal pump are shown in Table 1. The PAT is composed of five main parts: inlet volute, guide vane, impeller, and export. The computational domain grids are all structured grids generated by commercial software ANSYS ICEM 18.0. Figure 1 shows the computational domain of the PAT model. To ensure accurate calculation of the boundary area, the boundary layer grid is encrypted and the y+ value is controlled to be less than 200, as shown in Figure 2.

ANSYS CFX 18.0 was used for the steady and unsteady numerical simulation of PAT. The static pressure was set as the inlet boundary condition, the outlet boundary condition was set as the mass flow rate, all walls were set as nonslip wall surfaces, the frozen rotor is adopted at the rotor-stator interfaces, and the convergence criterion, RMS (root mean square), is set as 10−5. The turbulence model adopted the SST k-ω turbulence model which was widely used in turbomachinery [27]. For unsteady numerical simulation, the time step was set to the time it takes for the impeller to turn 1°, which was 5.747 × 10−5 s. And, the steady simulation result is adopted to the initial value of the unsteady simulation.

Four grids with 4.7 × 106, 6.57 × 106, 8.05 × 106, and 9.4 × 106 were selected to validate the independence of the grid. Each set of grids operates under 0.8Qd, 1.2Qd, and 1.6Qd conditions. The dimensionless pressure drop coefficient is calculated by the following equation:where and are the outlet and inlet pressures (Pa) and is the circumferential velocity of impeller tail edge (m/s). Table 2 shows the calculation result for different calculation grid numbers. The pressure drop coefficient calculated by the first three groups of grids has a significant increase, and the pressure drop coefficient calculated by the third and fourth groups of grids has little change. It shows that when the number of grids is greater than 8.05 × 106, the increase in the number of grids has little effect on the calculation results. So, to fully consider the computing resources to ensure the accuracy of the calculation, the fourth set of grids is selected for the numerical simulation of PAT.

3.2. Experimental Verification

The flow rate is measured by the Yokogawa AE215 flowmeter and the measurement error is ±0.5%. The pressure is measured by the Yokogawa EJA510A pressure sensor, and the measurement error is ±0.075%. JN338 type torque sensor is used to measure the torque, and its measurement error is ±0.2%. The unsteady flow field in PAT is studied and analyzed by numerical simulation. Figure 3 shows the test rig of model PAT and the comparison between numerical simulation results and experimental results. It can be seen that the numerical simulation results have the same trend as the experimental results, with the maximum error not exceeding 5%. Figure 4 shows the comparison between the numerical simulation and the experiment of the pressure pulsation at the P3 and P4 monitoring points at 1.0Qd. As can be seen from the figure, due to the effect of rotor-stator interaction, the pressure pulsation coefficient shows periodic fluctuations, and the time-domain curve distribution trend of numerical simulation of pressure pulsation is similar to the experiment. The time-domain distribution of the pressure pulsation at the monitoring point P3 agrees well with the experimental value. The time-domain distribution of monitoring point P3 is quite different from the experimental value. And the pressure pulsation coefficients at the monitoring points P3 and P4 are higher than the test values. However, the amplitude of the pressure pulsation in the numerical simulation is higher than the experimental value at both the leaf frequency and the harmonic frequency.

In summary, the monitoring positions can be accurately arranged during the numerical simulation. Moreover, factors such as gap leakage and disc friction loss were not considered in the numerical simulation process, so the numerical simulation pressure pulsation result was higher than the experimental value. But the distribution trend of numerical simulation results agrees well with the experiment. Therefore, the meshing, boundary conditions, and selection of turbulence models of each overcurrent component in the numerical method are reliable.

3.3. Entropy Generation Theory Verification

The total energy loss in PAT can be accurately obtained by the method of pressure drop. The entropy production theory can be used to analyze the internal loss of the PAT and its mechanism. The total entropy production can be calculated in equation (9) and the mechanical energy loss of the PAT can be calculated by the following equation:

The runner shaft power can be calculated in the following equation:

So the accuracy of entropy production theory can be verified by these two equations, calculated as follows:

The left side of the above equation is defined as Ploss, and the accuracy of entropy generation theory can be verified by comparing Ploss and Spro.

Figure 5 shows the energy loss of each component calculated by two different methods and the error of the results calculated by two methods. The energy loss calculated by the two methods tends to decrease first and then increase. The energy loss reaches the minimum value when Q = 40 m3/h. This is because the design operating condition of the pump used as PAT is Q = 40 m3/h. Under this working condition, the performance is the best and the internal flow is the most stable. At low flow rates, many vortices were generated in the PAT to block the flow passage, which caused higher losses. At the flow rates, the turbulence in the PAT was severe due to the high velocity, resulting in a high energy loss. The error between the entropy generation theory and the pressure drop method is less than 5% at Q = 48. The error can reach up to 10% when the flow rate is small. This is because the flow is stable when the operating conditions are designed, and the small flow and large flow are complicated. The entropy generation theory still has errors in some complex flow situations. In general, the entropy generation theory can be used to predict the energy loss inside the PAT.

4. Results and Discussion

Figure 6 shows the identification of the blades and flow passages, among which the impeller blades are RB1, RB2, ..., RB6. The flow channel between the blades RB1 and RB2 is RP1. The naming rules in the guide vanes are similar to the impeller.

Figure 7 shows the total entropy production and wall loss under different flow conditions. The wall loss increases linearly with the increase in the flow rate, and the entropy production decreases first and then increases with the increase in the flow rate and reaches the minimum value at Qd. This is because the fluid flow in the PAT is the best and the energy loss is the lowest at Qd. Besides, the loss of entropy generation inside the PAT is higher than the loss near the wall.

Figure 8 shows the energy loss obtained by the pressure drop method among the three main components of the PAT over one impeller rotation period. The energy loss of each component shows obvious periodic fluctuations. There are 6 obvious wave crests and troughs on all three parts, which are the same as the number of impeller blades. And the loss of the three components reaches the crest and trough at the same time. Comparing the amplitude of the energy loss fluctuation of the three components, it is found that the amplitude of the loss in the impeller is the highest, and the amplitude of the loss in the volute is the lowest. Also, there are 5 smaller wave crests between the two wave crests of the impeller and the guide vane component, the number of which is the same as the number of the guide vane blades. The fluctuation between the two main wave crests of the volute component is not very obvious. This shows that the internal loss of PAT is affected by rotor-stator interaction, the impeller is most affected by rotor-stator interaction, and the volute is weakly affected by the interaction. Besides, the impeller is strongly affected by the rotor-stator interaction of the impeller-guide vane. The guide vanes are also affected by the impeller-volute and impeller-guide vane rotor-stator interaction. The volute is strongly affected by the impeller-volute rotor-stator interaction. The interaction of the impeller-volute is mainly the interaction between the impeller and the volute tongue.

Figure 9 shows the entropy generation rate at three different moments in the volute at 48 m3/h. The choices of three time periods are shown in the black vertical lines in Figure 8. High energy loss occurs near the volute tongue. After the fluid passes through the tongue, energy loss begins to occur on the wall of the spiral section of the volute and increases in the circumferential direction. Also, energy loss also occurs at the exit of the volute near the leading edge of the guide vane. Compared with the velocity and streamline distribution in the volute (see Figure 10), the flow field at the volute tongue is complicated. Since the trailing edge of the guide vane blade of the centrifugal pump was not rounded, the fluid will impact the guide vane vertically. Part of the fluid directly bounces back and impacts with the fluid in the spiral section of the volute. The two consecutive impacts complicate the flow field near the tongue and form two low-velocity zones and high-velocity gradients. This causes high energy loss at the volute tongue. In the spiral section of the volute, the wall of the volute forces the fluid to change the direction of velocity; in addition, the influence of the boundary layer resulted in a large flow velocity change and a high-velocity gradient near the volute wall, which caused high losses. The reduction of the cross-sectional area of the volute leads to a certain increase in the velocity gradient of the cross section of the volute, which caused increased losses. When the volute section is too small, the impact between the fluid and the guide vane blades will also become stronger, so obvious energy loss occurs at the guide vane blade DB2. There was no obvious change in the energy loss everywhere in the volute at different times, which shows that the rotor-stator interaction will not have a significant impact on the local loss in the volute.

Figure 11 depicts the distribution of entropy production on the guide vane at different times. The energy loss of the flow passage near the suction surface of the guide vane blade close to the tongue is significantly higher than that of other flow passages. Higher energy losses are also produced near the trailing edges of the 5 blades at the exit of the guide vane. Also, the energy loss near the wall of the guide vane blade DB2 is significantly higher than other blades. Due to the complex flow at the separation tongue, the liquid flow angle near the suction surface of the guide vane blade 1 is too large. The guide vane velocity distribution graph is shown in Figure 12. As a result, flow separation occurs at the leading edge of the suction force of the blade, and a wide range of low-speed areas appear, which affects the flow of the entire flow passage and causes a high energy loss on the suction surface of the blade DB1. The energy loss to the trailing edge of the blade is due to the round shape of the trailing edge of the blade, which easily forms a wake flow, which causes a high loss. The reduction of the flow cross-section of the upstream volute causes a certain difference between the angle of the fluid entering the flow channel and the angle of the blade, which results in a strong impact of the fluid on the blade and high energy loss.

The entropy production rate at the exit of the guide vane changes at different times, and the wake loss changes significantly. The energy loss near the volute has a little change within three moments. This is because when the impeller blades rotate to the vicinity of the wake, the transmission of the wake vortex will be blocked. When the wake hits the blade, the loss in the guide vane and impeller will increase. The blade inlet is far away from the impeller, so the rotation of the impeller has little effect on the energy loss of the guide vane inlet. This shows that the influence of the guide vane by the impeller-guide vane rotor-stator interaction is mainly caused by affecting the wake of the guide vane. Also, the loss at the guide vane channel DP1 has changed significantly at different times. This is because the fluid in the guide vane channel DP1 causes a strong impact on the impeller blades due to the influence of the volute tongue. This is the main reason why the guide vane is affected by the rotor-stator interaction of the impeller and the volute.

Figure 13 shows the energy loss distribution of different impeller spans at T2 and T3 under Q = 48 m3/h. The losses in the impeller are concentrated near the entrance and exit of the impeller. And the energy loss near the shroud and the hub is significantly higher than the loss in the middle section. Thus, the impeller inlet is strongly affected by the rotor-stator interaction of the impeller-guide vane. In the streamline diagram (as shown in Figure 14), it was found that a large range of vortices appeared near the working surface of the impeller blades, causing nearly half of the impeller flow passage to be blocked, therefore, high energy loss near the impeller inlet. The loss at the exit of the impeller is caused by the blade wake flow. When it is close to the front and rear cover, the loss caused by the increased speed gradient is greater than the loss in the middle section. When the exit of the blade is close to the front cover, a wide range of low-speed zone appears, and the velocity distribution in the other two sections is more uniform. Therefore, the impeller outlet loss is the highest near the front cover. Although vortexes appear on the pressure surface, the loss of the pressure surface is lower than the energy loss of the suction surface due to the low flow velocity.

Comparing the loss in the impeller at different times, it can be seen that the position where the loss changes greatly is near the impeller inlet. The energy loss at the impeller inlet at T2 is significantly higher than the energy loss at the impeller inlet at T3. And the energy loss of the impeller flow path increases when passing through the guide vane flow path DP1. This shows that the flow in the impeller is affected by static and dynamic interference, is greatly affected by rotor-stator interaction of the impeller-guide vane, and is less affected by the rotor-stator interaction.

5. Conclusion

This paper uses the entropy production theory to obtain the detailed situation of the loss caused by PAT during operation and analyzes the phenomenon of rotor-stator interaction in PAT. Find the following conclusions:(1)The largest energy loss is the guide vanes, followed by the impellers, and the loss of the volute is relatively small in PAT. The asymmetric structure of the volute also affects the hydraulic loss of downstream components.(2)The loss in the volute was mainly composed of friction loss and impact loss at the tongue. The loss in the guide vanes was due to the flow separation in the flow channel, the impact of the leading edge of the blade, and the blade wake. The loss in the impeller was mainly due to flow separation and wake flow.(3)Rotor-stator interaction has the greatest influence on the vicinity of the interface between the impeller and the guide vane. The impeller is mainly affected by the rotor-stator interaction of the impeller and the guide vane, and the guide vane is affected by the rotor-stator interaction of the impeller and the volute. Rotor-stator interaction is one of the main causes of PAT vibration. This study is beneficial to study the dynamic and static interference and PAT vibration to improve the stability of the PAT.

Nomenclature

BEP:Best efficiency point
CFD:Computation fluid dynamics
PAT:Pumps as turbines
s:Specific entropy
Φ:Viscous dissipation term
ΦΘ:Dissipation term
T:Temperature (K)
:Inlet pressure (Pa)
:Outlet pressure (Pa)
Ploss:Energy loss in the runner (W)
:The entropy production by average velocity (W/m3)
:The entropy production by the pulsation velocity (W/m3)
Spro:Total entropy production (W)
:The entropy production caused by the pulsation velocity (W)
:The entropy produced by the time-averaged velocity (W)
Spro,w:The local entropy production rate for the wall
:Velocity vector (m/s)
ρ:Density (kg/m3)
ε:Turbulent dissipation rate (m2/s3)
Wsft:Runner shaft power (W)
μeff:Effective dynamic viscosity (N·s/m2)
μt:Turbulence dynamic viscosity (N·s/m2)
:Wall shear stress (Pa)
φ:Pressure drop coefficient
Qd:Flow rate (m3/s)
H:Head (m)
n:Rotating speed (rev/min)
D1:Outer diameter of impeller (mm)
b1:Impeller blade width (mm)
z1:Number of impeller blades
D2:Diameter of guide vane inlet (mm)
D3:Diameter of guide vane outlet (mm)
b2:Guide vane blade width (mm)
z2:Number of guide vanes
D4:Diameter of volute base circle (mm)
u2:Circumferential velocity of impeller tail edge (m/s).

Data Availability

The data used to support the findings of this study are available from the corresponding author upon request.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the Scientific Research Program funded by the Shaanxi Provincial Education Department (Grant no. 20JK0957), the Natural Science Basic Research Plan in Shaanxi Province of China (Grant no. 2020JM-645), and Xijing University Special Research Foundation (Grant no. XJ17T09).