Abstract

This work analyzes the causes of the slip phenomenon in the impeller on the basis of the internal flow mechanism. Detailed optical measurements of the flow inside the rotation passages of a five-bladed centrifugal pump impeller are obtained through particle image velocimetry (PIV). On the basis of experimental data, the deviation coefficient of slip velocity is proposed and then revised according to the slip factor calculation formula of Stechkin. Results show that, at the same rotation speed, the slip factor increases with the flow rate and reaches the maximum value at 1.0 QBEP flow rate. At different rotation speeds, the slip factor increases with the rotation speed and shows a relatively large variation range. Moreover, a revised slip factor formula is proposed. The modified model is suitable for the correction of slip factor at part-load flow rates and serves as a guide for the hydraulic performance design and prediction of centrifugal pumps.

1. Introduction

The slip factor is an important parameter in the hydraulic design of a centrifugal pump, and it can be obtained experimentally and theoretically. However, obtaining a relatively accurate solution for this phenomenon is difficult due to various factors affecting pump parameters. Over the past years, researchers have studied this phenomenon and put forward corresponding calculation formulas. Based on the 2D potential flow theory, Busemann [1] proposed a theoretical analytical equation for the slip factor of a backward curved blade. Stodola et al. [2, 3] considered that the slip phenomenon is mainly caused by the vortex in the channel. Based on the theory of velocity circulation and vortex strength, they also developed a formula to solve the slip factor. Weisner [4] conducted a comprehensive arrangement and statistical analysis on the experimental results of 65 impellers and obtained the empirical calculation formula for the slip factor. Balje [5] provided a simplified formula that is only related to the number of blades and the ratio of inlet and outlet radius. In Stanitz’s [6] formula, the blade number, outlet axial velocity, and blade outlet angle are conducive to improving the accuracy. Pfleiderer [7] first analyzed the differences in the relative velocity and then the theoretical head of impeller between finite and infinite blades. Stechkin [8] has improved the Pfleiderer formula by taking the empirical coefficient from the Pfleiderer calculation formula as constant for easy application. Fujie [9] proposed the theoretical analytical formula of slip factor based on the 2D potential flow theory; Backstrom [10] put forward the hypothesis of single relative eddy (SRE), which considers that the slip factor is caused by a relative vortex within the entire impeller channel, instead of a single vortex in a single blade channel.

Previous studies on slip factor mainly considered the influence of impeller geometry [1113]. However, this phenomenon is also affected by the flow rate [1416], rotation speed, and fluid viscosity [17, 18]. At present, the main research methods for slip factor include theoretical analysis, numerical simulation [1923], and experimental research [24, 25]. Qiu et al. [26] proposed a unified slip model for axial, radial, and mixed flow impellers. The core assumption of this model is that the slip velocity at the impeller outlet is mainly derived from the blade load near the impeller outlet and the relative eddy current in the impeller channel. When blade rotation occurs at the impeller outlet, the flow rate becomes an important variable that affects the slip factor. By using the Euler equation in turbomachinery and related eddy current theory, Ji et al. [27] derived a correlation for the slip factor of radial and mixed flow impeller. This formula considers the effect of geometric parameters and flow rate to accurately describe the changing trend of slip factor with the flow rate. Caridad et al. [28] conducted a numerical simulation on a centrifugal pump and used the results to develop a method that predicts the pump head. Li [29] researched the influence of flow rate and viscosity on the slip factor of heavy oil centrifugal pump by using numerical simulation and laser doppler velocimetry (LDV) technology. The results show that the change of slip factor largely depends on the flow rate and is minimally affected by the liquid viscosity. By using probe technology, Memardezfouli et al. [30] studied the outlet slip of five industrial pumps at different flow rates. They found that the experimental values correspond with the theoretical values under the design pump conditions but show significant differences under off-design conditions. The correlation formula of slip factor under off-design condition is derived by defining the flow distortion coefficient.

Owing to the large deviation of slip factor at part-load flow rates, most of the above articles consider the influence of flow rate on slip factor and put forward a new calculation model for the corresponding slip factor. However, only a few works considered the influence of slip velocity on slip factor. The present work attempts to solve these problems using the particle image velocimetry (PIV) technique, which offers information on instantaneous spatial flow structures, and the revised calculation formula of Stechkin [8]. In addition to traditional geometric parameters, such as the number of blades and radius ratio, this paper considers slip velocity, an important influencing element of the slip factor. On the basis of the experimental data, the slip velocity deviation coefficient is defined, and a new formula for the slip factor is proposed. This paper has a certain guiding importance for the hydraulic performance design and prediction of centrifugal pumps.

2. Experiment Platform and Measurement Schemes

2.1. Slip Phenomenon

When the theoretical impeller has an infinite number of blades, the relative velocity of the fluid is always tangent to the blade surface, the outlet flow angle is equal to the geometrical blade angle, and the fluid has no flow slip quantity. When the actual impeller is finite, the relative velocity of the fluid is no longer tangent to the blade surface, and the outlet flow angle is less than the geometrical blade angle. Hence, the fluid produces a flow slip quantity , which leads to flow slip as shown in Figure 1.

According to the different definitions of slip factor, it can be divided into two forms. The first definition is commonly used in the United States:

The famous empirical formulas in this definition include Stodola [2] and Weisner [4]. The second definition is commonly used in Europe:

The famous empirical formulas in this definition include Busemann [1] and Stechkin [8]. In the present paper, the common definition of slip factor in Europe is adopted, and the local slip factor in the impeller passage is defined as follows:

The average slip factor is defined as follows:

2.2. Measurement Schemes

The PIV system of this experiment is produced by the TSI company of the United States and is set as follows. The time delay between two laser pulses is 130 μs. The CCD camera with a spatial resolution of 2048 px × 2048 px is used for image acquisition, and its spatial resolution is 102.38 μm/px. The thickness of the laser source is 1 mm. Solid alumina particles with a density of 1050 kg/m3 and diameters ranging from 20 to 60 μm are selected as seeding particles. For convenience, the impeller channels are marked as 1 to 3. The flow data of the No. 2 channel are selected for analysis. As shown in Figure 2(a), , , and u are defined as the relative, absolute, and circumferential velocities, respectively.

The middle span Sm of the impeller is selected for the measurement. As shown in Figure 2(b), Sm is located in the middle of the front and rear cover plates of the impeller. The impeller is calibrated by using the pixels and the actual physical size of the circle located on the impeller and the diameter is 30 mm. The shooting speeds are divided into five speeds (600, 800, 1000, 1200, and 1400 r/min). On the basis of the performance curve of the experimental pump, PIV data are obtained at six flow rates (0.2, 0.4, 0.6, 0.8, 1.0, and 1.2 QBEP), and 300 double pictures are captured under operating conditions. Li et al. [31] studied the influence of flow patterns on the hydraulic performance and energy conversion characteristics in a centrifugal pump. The results demonstrated that the velocity at a sample size of n = 300 is well converged. For each flow rate, a total of 123  123 vectors were obtained. The experimental data are then processed as follows. First, the absolute velocity is derived by the internal software based on the following equation:

The phase-averaged relative velocity is the difference between the two figures, namely, the absolute velocity and the local circumferential impeller speed at each point, as shown in the following equation:

The velocity data of each calculation point of the channel arc are obtained through interpolation, which is based on least squares method. According to the definition of the local slip factor, the slip factor values of each calculation point of the channel arc are calculated. In the slip factor processing program, the blades of the suction side and pressure side of the No. 2 channel are divided into 44 calculation points, and the connecting line between the suction surface and pressure surface formed 44 channel arcs. Each channel arc is split into 21 slip factor calculation points as shown in Figure 2(c). A total of 924 slip factor calculation points can be obtained in the whole flow channel. Thus, the slip factor processing program can acquire the distribution of slip factors in the entire channel. Given that the slip factor is mainly the data at the outlet of the impeller, the 21 slip factor calculation points of the channel arc at the impeller outlet are averaged to obtain the slip factor value of the whole impeller.

2.3. Test Pump and Test Rig

The test pump and test rig have been introduced in the study of Li et al. [31]. The five-blade impeller and rear cover of the experimental pump are composed of plexiglass, which provides an optical channel for PIV measurement. For a reduced laser reflection, the nontransparent area is blackened, and the surface of the transparent impeller is polished. The main parameters of the impeller are shown in Table 1. The centrifugal pump is a special volute-free structure used to obtain the flow data in the entire impeller channel. The test pump and test rig are shown in Figure 3. A control valve is set at the outlet of the test pump to control the flow rate, and the electromagnetic flowmeter is used to measure the flow rate. The inlet and outlet pressure levels of the centrifugal pump are measured by a pressure transmitter. Speed, torque, and shaft power are calculated by the frequency conversion control cabinet. The main equipment includes the following: frequency conversion control cabinet, three-phase asynchronous motor, torque meter, valve, electromagnetic flowmeter, pressure sensor, water storage tank, laser, camera, and computer.

2.4. Performance Measurements

Figure 4 presents the performance curves of the centrifugal pump at five rotation speeds. The Reynolds number for the pump at various rotation speeds is as follows: 1.25 × 106 at 600 r/min, 1.67 × 106 at 800 r/min, 2.11 × 106 at 1000 r/min, 2.51 × 106 at 1200 r/min, and 2.92 × 106 at 1400 r/min. υ is the kinematic viscosity of the fluid, and Ω is the impeller angular velocity. On the basis of the performance curves of the test pump, the flow rates at best effect point QBEP and head HBEP at various rotation speeds are as follows: 1.0 m3/h and 1.0 m at 600 r/min, 1.7 m3/h and 1.5 m at 800 r/min, 1.8 m3/h and 2.67 m at 1000 r/min, 2.3 m3/h and 3.7 m at 1200 r/min, and 2.5 m3/h and 5.05 m at 1400 r/min.

3. Experimental Results and Discussion

3.1. Experimental Value of Slip Factor

Data on the slip factor of five rotation speeds are shown in Figure 5. The slip factor increases with the flow rate at the same rotation speed. When the flow rate is 1.0 QBEP, the slip factor reaches the maximum value. This trend is consistent with [30]. Under the same rotation speed, the variation range of the slip factor is relatively small. In addition, the slip factor increases with the rotation speed at different levels and shows a relatively large variation range.

3.2. Influence of Flow Rate on Slip Factor

The influence of flow rates on slip factor is analyzed when the rotation speed is 600 r/min. The distribution of the slip factor is depicted in Figure 5. The slip factor increases with the flow rate because its distribution in the channel is nonuniform (Figure 6). At the outlet of the channel, the value of the slip factor near the suction side is larger than that near the pressure side. With the decrease in the flow rate, the slip phenomenon near the outlet of the pressure side tends to be strengthened. At overload flow rate, the slip factor greater than 1 is distributed in the major region of the flow channel.

Stodola [2] believes that the slip phenomenon is mainly caused by the vortex, and the fluid flow in the impeller can be divided into two forms: the inertial flow through the fixed impeller and the vortex motion. Therefore, the phase-averaged relative velocity and flow streamlines are obtained and displayed in Figure 7 to explain the distribution of the slip factor at the outlet of the channel. At an overload flow rate of 1.2 QBEP, the streamlines follow the blade curvature without flow separation, and a low-velocity region develops at the blade outlet. This occurrence has been previously reported [32, 33]. The distribution of velocity is uniform at the outlet of the impeller and has minimal effects on the slip factor. This finding explains the value of slip factor greater than 1 being distributed in the major region of the flow channel. As a result, the average slip factor at the outlet of the channel reaches the maximum. At nominal flow rate and part-load flow rate of 0.8 QBEP, a clockwise vortex occurs on the suction side of the outlet flow channel. The rotation direction of the vortex is the same as that of the impeller. After the inertial flow through the fixed impeller and the vortex motion are combined, the circumferential component of the absolute velocity is approximately equal to the . According to the definition of the slip factor, slip did not occur. Therefore, the slip factor near the outlet of the suction side is relatively large. However, no vortex can be found at the pressure side of the outlet channel. The direction of streamline movement is opposite to the rotation of the impeller, thus reducing the average slip factor at the outlet. At part-load flow rates of 0.6, 0.4, and 0.2 QBEP, a clockwise vortex occurs at the outlet of the suction side of the channel. Moreover, a counterclockwise vortex is found on the outlet of the pressure side of the channel. The rotation direction of this vortex is opposite to that of the impeller. After the inertial flow through the fixed impeller and the vortex motion are combined, the circumferential component of the absolute velocity is smaller than . According to the definition of the slip factor, the slip phenomenon occurs in the flow channel. Thus, the slip factor near the outlet of the pressure side is relatively small. With the decrease in the flow rate, the vortex intensity increases gradually and leads to aggravate slip. This finding explains the slip factor near the outlet of the pressure side decreasing with the flow rate. This phenomenon is applicable to the decrease in the average slip factor at the outlet under three flow rates (0.6, 0.4, and 0.2 QBEP). As the turbulence is random, the flow fields in different channels are somewhat different.

3.3. Influence of Rotation Speed on Slip Factor

The influence of rotation speeds (600, 800, 1000, 1200, and 1400 r/min) on slip factor is analyzed at a part-load flow rate of 0.2 QBEP and nominal flow rate. At a part-load flow rate of 0.2 QBEP, the slip factor distribution in the impeller channel is plotted in Figure 8. The slip factor near the outlet of the blade suction side is relatively large, whereas that near the outlet of the blade pressure side is relatively small. With the increase in the rotation speeds, the slip area at the outlet of the blade pressure side presents a decreasing trend.

The phase-averaged relative velocity and flow streamlines are obtained and exhibited in Figure 9. Two vortexes with opposite rotation directions can be found in the channel at five rotation speeds. Each channel is filled with a low-velocity region, resulting in nonuniform velocity distribution at the outlet of the impeller. The slip velocity at the outlet of the impeller is illustrated in Figure 10. The circumferential component of absolute velocity at infinite blade number is shown in Figure 11. In Figures 10 and 11, L = 1 represents the suction side, and L = 21 represents the pressure side. The channel arc at the outlet is divided into 21 calculation points of slip factor.

Figures 10 and 11 show that with the increase in the rotation speeds, the slip velocity at the outlet of the channel also increases gradually, but the variation range is relatively small. The circumferential component of the absolute velocity also increases with the rotation speed when the number of blades is infinite, and the variation range is relatively wide. In accordance with equation (2) of the definition of slip factor, the ratio of decreases with the increase in the rotation speeds. Therefore, the slip factor increases with the rotation speeds.

The slip factor distribution in the impeller channel for five rotation speeds at a nominal flow rate is acquired and shown in Figure 12. The phase-averaged relative velocity and flow streamlines are obtained and exhibited in Figure 13. The slip factor is relatively large due to a clockwise vortex on the pressure side of the outlet at the rotation speed of 600 r/min. For the other four rotation speeds at the nominal flow rate, the streamline is distributed along with the blade curvature. No vortex structure, which only exists in the local low-velocity region, can be observed. Therefore, the slip factor is quite large and relatively uniform in the channel. Hence, the slip phenomenon is not evident. In addition, the velocity distribution at the outlet of the channel is uniform.

The slip velocity at the outlet of the impeller is illustrated in Figure 14, and the circumferential component of absolute velocity under infinite blades is revealed in Figure 15. The slip velocity increases with the rotation speeds, but its variation range is relatively small. When the number of blades is infinite, the circumferential component of the absolute velocity also increases with the rotation speeds, and its variation range is relatively large. This finding is similar to the distribution rule for the five rotation speeds at a part-load flow rate of 0.2 QBEP. In accordance with equation (2) of the definition of slip factor, when the rotation speed increases, the ratio of decreases, and the slip factor increases.

4. Empirical Formula and Revision of Slip Factor

Stechkin’s [8] empirical formula is used to revise the definition of slip factor as follows:

Owing to its importance, the deviation of slip factor at part-load flow rate is strivingly corrected in this work. Memardezfouli [30] transformed the real-flow model to a uniform flow one by defining the flow distortion coefficient AN, as the ratio of to the same component for uniform flow, . In the present paper, the slip factor is found to be related to the change of . Hence, the slip velocity deviation coefficient is defined as follows:where is the slip speed at the nominal flow rate and is the slip speed at the overload and part-load flow rates. At nominal flow rate, the slip velocity deviation coefficient is equal to 1. The modified slip factor is obtained by incorporating the slip velocity deviation coefficient into the Stechkin empirical formula:

According to the experimental data, the slip velocity deviation coefficients under different flow rates are obtained and shown in Table 2. The impeller C has been verified in the paper of Memardezfouli [30]. For impeller C, the number of blades (Z) is 6, and the ratio of inner and outer diameter (D1/D2) is 0.46. Figure 16 shows the correction results of the slip factor changing with the flow coefficient. Owing to the large deviation of the slip factor at part-load flow rates, the σskk modified model with a good correction effect at part-load flow rates is proposed. The difference between the experimental and theoretical slip factor is negligent at the nominal flow rate but is important at the part-load flow rates. Therefore, the σskk modified mode is highly suitable for correcting the slip factor at part-load flow rates.

5. Conclusion

(1)PIV is used to measure the velocity field of a five-blade centrifugal pump and obtain the slip factor at five rotation speeds. The results show that the slip factor increases with the flow rates at the same rotation speed and increases with the rotation speeds at varying rotation speeds. The changing trend is evident.(2)The influence of flow rate on slip factor is analyzed when the rotation speed is 600 r/min. When the flow rate is less than 1.2 QBEP, a clockwise vortex occurs near the suction side at the outlet of the blade, resulting in a large slip factor in this region. When the flow rate is less than 0.8 QBEP, a clockwise vortex can be observed near the pressure side at the outlet of the blade, resulting in a relatively small slip factor in this region. The existence of these vortexes affects the average slip factor. Therefore, the slip factor increases with the flow rate and reaches the maximum value at the nominal flow rate.(3)Comparison of the data from five rotation speeds at a nominal flow rate and part-load flow rate of 0.2 QBEP revealed that the slip velocity increases with the rotation speeds, but its variation range is relatively small. When the number of blades is infinite, the circumferential component of the absolute velocity also increases with the rotation speeds, and its variation range is relatively large. This phenomenon leads to the slip factor increasing with the rotation speeds.(4)On the basis of the experimental data, the deviation coefficient of slip velocity is proposed, and the Stechkin slip factor calculation formula is modified. The results show that the σskk modified model is suitable for the correction of slip factor at part-load flow rates. This work has a certain guiding importance for the hydraulic performance design and prediction of centrifugal pumps.

Nomenclature

D:Diameter (mm)
r:Radius (mm)
b:Blade height (mm)
Z:Number of blades
∞:Theoretical impeller
Q:Flow rate, m3/h
QBEP:Flow rate at best effect point (m3/h)
H:Pump head (m)
σ′:Slip factor, American definition
σ:Slip factor, European definition
:Relative velocity (m/s)
:Absolute velocity (m/s)
u:Circumferential velocity
:Circumferential component of absolute velocity (m/s)
:Radial component of relative velocity (m/s)
:Circumferential component of relative velocity (m/s)
Δ:Slip velocity (m/s)
N:Number of calculation points of slip factor
α:Absolute flow angle (°)
β:Relative flow angle (°)
x:Horizontal axis (mm)
y:Vertical axis (mm)
φ:Impeller discharge flow coefficient, dimensionless = c m2/ u2
cm:Meridional velocity
j:Slip factor calculation point
Subscripts
1:Impeller inlet
2:Impeller outlet.

Data Availability

The experimental data used to support the findings of this study are included within the article.

Conflicts of Interest

The authors declare that there are no conflicts of interest regarding the publication of this paper.

Acknowledgments

This work was supported by the National Natural Science Foundation of China (Grant numbers: U2006221 and 51806197), the Natural Science Foundation of Zhejiang Province (Grant number: LY20E060006), and Key Research and Development Program of Zhejiang Province (Grant number: 2021C05006). The supports are gratefully acknowledged.