Abstract

In engineering, the highest operating head of the pumping station is usually controlled to be slightly lower than the lowest saddle bottom head of the axial-flow pump. However, in the practical operation, it is found that the highest operating head of the pumping station is obviously lower than the saddle bottom head of the pump device, which leads to the reduction of the operating range of the pumping station. To investigate the difference of lowest saddle bottom head between axial flow pump and axial flow pump device and apply it correctly, the energy performance tests of the TJ04-ZL-06 hydraulic model and its corresponding pump device were carried out to obtain the external curves, and numerical simulation was carried out to analyze and compare the internal flow field and pressure distribution. The results show that when the flow rate decreases, the first saddle-shaped region of the axial-flow pump and the saddle-shaped region of the pump device are caused by the decrease of the lift coefficient due to the increase of the attack angle between flow and blade. When the flow rate is less than 0.32Qd, the influence range of backflow in the inlet pipe is large, which leads to the high-pressure zone near the wall of the inlet pressure measurement section during the pump performance test, and hence the second saddle-shaped region of the axial-flow pump is essentially a measurement illusion. It is suggested that the inlet pressure measurement section should be set at least 4Dp away from the inlet flange of the impeller when testing the performance of the axial-flow pump under the condition of small flow rate, and the first saddle bottom head of the axial-flow pump or the saddle bottom head of the corresponding pump device can be considered as the control value of the highest head of the pumping station.

1. Introduction

Axial-flow pumps play significant roles in irrigation, urban water supply, and flood control [13]. In the design of pumping stations, to maximize the pumping station’s efficiency and save energy, the pump must be equipped with an inlet and outlet passage to form a pump device. At present, in addition to a large number of studies on the operation of the axial-flow pump and axial-flow pump device under optimal conditions [47], some scholars have studied the operation and harm in the saddle-shaped region. It is found that there is vibration and noise when the axial-flow pump works in saddle-shaped region, which reduces the life of the pump and affects the safety of the pumping station [810]. Therefore, when selecting a pump for an axial-flow pumping station, the highest head of the pumping station must be lower than the lowest saddle bottom head of the axial-flow pump [11]. However, it is found in practical engineering that the highest head of the pumping station is much lower than the saddle bottom head of the pump device, which leads to the reduction of the operating range of the pumping station. The Q-H curves of a large number of axial-flow pumps and their corresponding pump devices are compared, and it is found that there are two saddle-shaped regions of the axial-flow pump [12], while there is only one saddle-shaped region of the pump device [13]. Typical Q-H curves of the axial-flow pump and axial-flow pump device are shown in Figure 1. In addition, the first saddle bottom head of the pump is close to the saddle bottom head of the corresponding pump device and is obviously higher than the second saddle bottom head.

At present, there are some research studies on the operation of the axial-flow pump and pump device in saddle-shaped region, but there is no study on the difference between the saddle-shaped region in Q-H curves of the axial-flow pump and its corresponding pump device. Yang et al. [12] conducted a model test of the axial-flow pump, found that there are two saddle-shaped regions in the head curve, and measured the flow field with PIV technology. The results showed that when the flow rate decreases to 0.55Qd, backflow starts to appear at the impeller inlet, and water flows out from the rim of the impeller. Zheng et al. [13] carried out a CFD analysis of the axial-flow pump device and found that a saddle-shaped region exists at 50%∼65% Qd. The axial-flow pump cannot operate stably under this condition. Goltz et al. [14] found that the impeller starts to stall when an inflection point occurs in the head curve by carrying out a test of hydraulic characteristics of a pump and captured some characteristics of the initial stall. Cheng et al. [15] conducted a CFD analysis of the pump device. The results showed that under small flow rate working condition, there is a large range of spiral backflow in the inlet passage. Wang et al. [16] found that the flow pattern in the impeller is well in the optimal condition. When the flow rate decreases to the saddle-shaped region, there is backflow at the impeller inlet.

Model test and numerical simulation are used in this study to investigate the cause of the difference between the saddle-shaped region of the axial-flow pump and its corresponding pump device. This study can provide suggestions for the axial-flow pump selection of large- and medium-sized pumping stations and the performance test of the axial-flow pump; the inlet pressure measurement section should be set at least 4Dp away from the inlet flange of impeller when testing the performance of axial flow pump under the condition of small flow rate; and the first saddle bottom head of axial flow pump or the saddle bottom head of corresponding pump device can be considered as the control value of the highest head of pumping station. This study has theoretical value and application value.

2. Model Test

2.1. Research Object

The inlet passage is the connection between the forebay and the impeller chamber. The purpose is to make the flow uniform and smooth and provide a good inlet condition for the impeller. The outlet passage is commonly used to recover kinetic energy [17]. The object of this study is an axial-flow pump (TJ04-ZL-06) and its corresponding shaft tubular pump device. The axial-flow pump has 3 blades and 6 guide vanes. The corresponding device consists of the aforementioned axial-flow pump, shaft inlet passage, and straight outlet passage.

The design requirements of the axial-flow pump are presented as follows: the impeller diameter D = 300 mm and the rotation speed n = 1450 r/min. The tip clearance of model pump is within 0.2 mm.

2.2. Experimental Setup

According to standards (e.g., SL140-2006, ISO 9906:2012, and ANSI/HI 14.6-2011), for the axial-flow pump model test, the inlet and outlet pressure measurement sections are set at twice the pipe diameter of the inlet and outlet flanges of the axial-flow pump, respectively, and the pipe diameter Dp = 350 mm. For the pump device model test, the inlet pressure measurement section is arranged at the inlet section of the shaft inlet passage, whereas the outlet pressure measurement section is arranged at the outlet of the straight outlet passage. Four measurement points are arranged on each measurement section; for the circular section, points are set at the ends of two vertical diameters, whereas for the rectangular section, points are set at the midpoints of four edges. The model diagram of the experiment is shown in Figure 2. The location of the pressure measurement sections is shown in Figure 3.

2.3. Analysis of the Test

To ensure the comparability of the results of the two model tests, both model tests were conducted on the general hydraulic model test bench of the China Water Resources Beifang Co., Ltd. The general hydraulic test bench is composed of a hydraulic circulation device, power control device, and data acquisition and computer analysis device. The tests were completed under the requirements of Chinese standards. The information about the test instruments is shown in Table 1. The tests were completed under the requirements of Chinese standards. The performance of the axial-flow pump and pump device at different blade mounting angles of −4°, −2°, 0°, and +2° was tested. The performance curves at different blade mounting angles were obtained from more than 15 test points.

The uncertainty of the head measurement system is less than ±0.1%, the uncertainty of the torque and rotation speed measurement system is less than ±0.1%, and the uncertainty of the flow rate measurement system is less than ±0.2%; according to the above uncertainty, the uncertainty of the test bench system is calculated to be ±0.25%. In this study, the data of the axial-flow pump device with blade angle of −4° under optimal conditions were collected 10 times, and the efficiency of each time was calculated, and then the average efficiency and standard deviation were calculated. Student’s t-distribution is adopted with a confidence of 95%, and the random uncertainty of the efficiency test is ±0.0728%. According to the system uncertainty and random uncertainty, the comprehensive uncertainty of the efficiency test is ±0.26%. Table 2 shows the efficiency data of 10 repetitions.

Figure 4 shows the model test results. The performance curves at different angles exhibit basically the same trend. Figure 5 compares the Q-H curves at different blade mounting angles. It can be found that there are two saddle-shaped regions of the pump, while there is only one saddle-shaped region of the corresponding pump device, and this is consistent with the results in [12, 13]. For the convenience of discussion, the saddle-shaped region with larger flow rate is called the first saddle-shaped region of the axial-flow pump, and the other with the smaller flow rate is called the second saddle-shaped region. The lowest head of the first saddle-shaped region is called the first saddle bottom head. The lowest head of the second saddle-shaped region is called the second saddle bottom head. The saddle bottom head of the pump and its corresponding pump device at different angles was statistically analyzed and is shown in Tables 3 and 4. It can be found that at different blade angles, the first saddle bottom head of the pump is close to the saddle bottom head of the corresponding pump device, while the second saddle bottom head is obviously lower than the saddle bottom head of the corresponding pump device.

As the trend of performance curves at different blade mounting angles is basically the same, in order to describe in detail, two curves with blade mounting angle of 0° were compared, and the x axis was replaced by the ratio of flow rate to design flow rate (Qd). The comparison is shown in Figure 6. Under the design condition, the efficiency of pump measured by the model test is 85.74%, and the head measured by the model test is 5.20 m. When the axial-flow pump is operating at 0.6∼1.2Qd, the head increases with decreasing flow rate, and the head exhibits one-to-one correspondence with the flow rate, so this region is called the stable operating region. When the flow rate is less than 0.6Qd, the curve fluctuates, forming two saddle-shaped regions. In this region, one head corresponds to multiple flow rate values, and thus this region is called the unstable operating region of the pump.

According to Figure 6, under the design condition, the efficiency of the pump device measured by the model test is 81.77%, and the model test head is 4.82 m. The head of the pump device head is 0.38 m lower than that of the pump, and the efficiency of the pump device is 3.97% lower than that of the pump. This is because under the design condition, the hydraulic loss of inlet and outlet passage is greater than that of inlet and outlet pipe, so the head and efficiency of the pump device will be lower. When the axial-flow pump device is operating at 0.6∼1.2Qd, the head increases with decreasing flow rate, and the curve exhibits a quadratic relationship. In this operation region, the head exhibits one-to-one correspondence with the flow rate, so this region is called the stable operation region. When the flow rate is less than 0.6Qd, the head decreases in the range of 0.45∼0.6Qd and then increases. In this region, the curve fluctuates, forming one saddle-shaped region, and the head corresponds to multiple flow rate values. Thus, this region is called the unstable operating region of the pump device. Comparing the Q-H curves of the pump and corresponding pump device, it is found that the two curves are basically the same in the stable operation region. In the unstable operation region, it is found that there are two saddle-shaped regions of the axial-flow pump, while there is only one saddle-shaped region of corresponding pump device.

3. Numerical Calculations

3.1. Computational Model

To research the difference of saddle-shaped region of Q-H curves, numerical calculations of the axial-flow pump with a blade angle of 0° and the axial-flow pump device with a blade angle of 0° were completed. The size of the calculation domain is the same as the model test, and the arrangement of pressure measurement points is also consistent with the model test. Figures 7 and 8 show the three-dimensional models of each component in the calculation domain.

3.2. Turbulence Model and Boundary Conditions

In recent years, the numerical simulation method has become a common method for the study of rotating machinery, such as the axial-flow pump [1822], centrifugal pump [2326], multistage multiphase pump [2729], pump as turbine [30], and so on. In this numerical simulation, the time-averaged Navier–Stokes equation is adopted to describe flow [31, 32]. Because of the advantage of the RNG k-ε turbulence model in solving the flow field of rotating machinery, the model is used in this paper [3337]. In this study, the inlet boundary condition was a standard atmospheric pressure, and the outlet was set as the mass flow outlet and all solid surfaces were set to wall. Water temperature is 25°C, regardless of heat transfer. The equations are shown below:where Gk represents the turbulent kinetic energy generated by the average velocity gradient; C1ε = 1.42 and C2ε = 1.68; μeff is the effective viscosity of turbulence; k is the turbulence energy; and ε is the dissipation rate of the turbulent kinetic energy.

3.3. Grid Division and Irrelevance Verification

According to the design parameters, inlet pipe, shaft inlet passage, water-guide cone, outlet pipe, and straight outlet passage are modeled, and the hexahedral structured grid is divided in ICEM software. The impeller and guide vane are modeled according to coordinate points and divided into structural grids in TurboGrid software. A large number of research studies show that the value of y+ and the setting of tip clearance have a great influence on the calculation accuracy [3840]. The tip clearance was 0.2 mm in this calculation. The y+ value of the impeller wall was 18. In this study, six sets of grids were divided for the calculation domain, and the efficiency under design condition was used as an index to judge grid independence. Figure 9 shows the verification results. The grid division is shown in Figure 10.

Figure 9(a) shows that the efficiency of the pump is almost constant when the grid number is more than 4.92 million. Figure 9(b) shows that the efficiency of the pump device is almost constant when the grid number is more than 6.29 million. After comprehensive consideration of calculation accuracy and computing resources, 4.92 million and 6.29 million were finally selected.

4. Results and Discussion

4.1. Numerical Simulation Results

To verify the accuracy of simulation, the experimental results were compared with the numerical simulation results. In this study, the head was calculated according to the Bernoulli equation. The efficiency was predicted according to the velocity, pressure, and torque on the impeller. Comparison between numerical and experimental results is shown in Figure 11.

The calculation formula of the head is expressed as follows:where and represent the total pressure of outlet and inlet pressure measurement section, respectively, Pa.

The calculation formula of hydraulic efficiency is expressed as follows:where is torque, N·m.

According to Figure 11(a), in the design condition, the experimental efficiency value is 85.74%, the efficiency predicted by the CFD is 85.62%, the absolute error of efficiency is 0.12%, and the relative error of efficiency is 0.14%. The experimental head value is 5.20 m, the head predicted by the CFD is 5.05 m, the absolute error of the head is 0.15 m, and the relative error of the head is 2.9%. According to Figure 11(b), in the design condition, the experimental efficiency value is 81.77%, the efficiency predicted by the CFD is 81.58%, the absolute error of efficiency is 0.19%, and the relative error of efficiency is 0.23%. The experimental head value is 4.82 m, the head predicted by the CFD is 4.61 m, the absolute error of the head is 0.21 m, and the relative error of the head is 4.3%. The error between numerical simulation and test is very small in both design and non-design conditions, indicating that the CFD method is reliable.

4.2. Internal Flow of the Impeller

This study established different blade spans; span represents the blade span, span=0 represents the hub position, and span=1 represents the rim positioned. Figure 12 shows the internal flow characteristics for the blade spans under different working conditions. According to Figure 12, when the flow rate is 1.0Qd∼0.8Qd, the flow pattern is smooth and ordered, the flow moves along the blade airfoil, and there is no backflow in the impeller. When the flow rate is 0.6Qd, there is a backflow zone near the impeller rim, but the backflow zone is small, whereas swirl appears near the impeller hub. When the flow is 0.32Qd, the flow pattern is disordered, and the backflow zone near the impeller hub and rim increases greatly. When the flow rate continues to decrease to 0.12Qd, the backflow zone in the impeller further increases, whereas the backflow zone near the impeller rim occupies almost the whole back of the blade. It was also found in reference [16] that there was backflow at the inlet of the impeller under the condition of small flow rate.

According to Figures 11 and 12, when the flow rate is about 0.6Qd, the head begins to decrease, the flow pattern of the impeller is disordered, and backflow phenomenon is observed. This is due to the impact angle of water flow and blade increasing with the decrease of flow rate. When the impact angle increases to a certain extent, the lift coefficient decreases. Therefore, the first saddle-shaped region of the axial-flow pump and the saddle-shaped region of the corresponding pump device are the external characteristic responses of the pump’s internal flow characteristics.

4.3. Internal Flow of Inlet Pipe and Shaft Inlet Passage

The backflow at the inlet section of the impeller enters the inlet straight pipe or the inlet passage, thus affecting the internal flow in the inlet pipe or the inlet passage. For further study, the internal flow of inlet pipe and shaft inlet passage are removed for comparison. Figure 13 shows the internal flow of inlet pipe. Figure 14 shows the axial velocity distribution of the inlet pressure measurement section of the axial-flow pump test. Figure 15 shows the internal flow of shaft inlet passage.

According to Figures 13 and 14, when the axial-flow pump is operating at 0.8∼1.0Qd, the flow in the inlet pipe is uniform and smooth without backflow. When the flow decreases to 0.6Qd, the backflow occurs at the impeller inlet near the rim, but it is small. When the flow rate is 0.32Qd, the backflow zone in the inlet pipe increases greatly along the radial and axial directions. The axial velocity distribution of the inlet pressure measurement section is uniform, and the axial velocity value is basically the same when the flow rate is greater than 0.6Qd. However, due to the influence of backflow, the axial velocity in the inner and outer sides of the section is inconsistent, and the axial velocity decreases from the center to the outer, which may lead to inconsistent pressure distribution in the section when the flow rate is less than 0.32Qd. Zheng et al. [15] also found the same backflow phenomenon but did not determine the influence range of the backflow in combination with the position of the inlet pressure measurement section.

According to Figure 15, When the axial-flow pump device is operating at 0.8∼1.0Qd, the flow in the inlet passage is uniform and smooth without backflow. When the flow decreases to 0.6Qd, the backflow phenomenon is observed. With the decrease of flow rate, the range of backflow zone increases along the axial and radial directions. Since the inlet pressure measurement section is set far away from the inlet of the impeller, it will not be affected by the backflow in all working conditions.

4.4. Internal Pressure Distribution

According to the above analysis, for the axial-flow pump test, the flow field of the inlet pressure measurement section will be affected by the backflow at low flow rate, resulting in an uneven axial velocity distribution in the section. The pressure distribution of inlet pressure measurement section under different working conditions was compared. For further study, the pressure distribution of different inlet pressure measuring sections was compared and the head was calculated according to the pressure of different sections. The position of inlet pressure measurement section A is the position of the inlet pressure measurement section in the axial-flow pump model test, and the inlet pressure measurement section B is 4Dp from the inlet flange of the impeller chamber. The location of different sections is shown in Figure 16. The pressure distributions of different sections are shown in Figures 17 and 18, respectively. Different Q-H curves are shown in Figure 19.

When the flow rate is 1.0Qd∼0.6Qd, the pressure distribution in inlet pressure measurement section A is uniform, and the pressure near the outside and inside of the inlet pipe is almost the same. When the flow is lower than 0.32Qd, the pressure is not uniform, and a high-pressure zone is formed in the rim of inlet pipe. With the decrease of the flow rate, the pressure in the high-pressure zone becomes higher. However, the pressure distribution on the inlet pressure measurement section B is uniform under all working conditions, and the pressure values at each point on the section are basically the same.

As shown in Figure 19, when the inlet pressure measurement section A is set at the same position as that in the model test, two saddle-shaped regions appear in the Q-H curve. However, when the inlet pressure measurement section is set at section B, the Q-H curve has only one saddle-shaped region. This demonstrates that when the flow rate is less than 0.32Qd, the influence range of backflow in the inlet pipe is large, which leads to the high-pressure zone near the wall of the inlet pressure measurement section, and hence the second saddle-shaped region of the axial-flow pump is essentially a measurement illusion.

5. Conclusions

In this study, the energy performance tests of the TJ04-ZL-06 hydraulic model and its corresponding pump device were carried out to obtain the external curves, and numerical simulation was carried out. The conclusions are as follows:(1)There are two saddle-shaped regions of the axial-flow pump, while there is only one saddle-shaped region of the corresponding device. The first saddle bottom head of the axial-flow pump is close to the saddle bottom head of the corresponding pump device, while the second saddle bottom head is obviously lower than the saddle bottom head of the corresponding pump device.(2)The first saddle-shaped region of the axial-flow pump and the saddle-shaped region of the corresponding pump device are caused by the decrease of the lift coefficient due to the increase of the attack angle between flow and blade. The second saddle-shaped region of the axial-flow pump is essentially a measurement illusion under small flow rate.(3)It is suggested that the inlet pressure measurement section should be set at least 4Dp away from the inlet flange of the impeller when testing the performance of the axial-flow pump with flow rate less than 0.32Qd, which avoids the measurement illusion under small flow rate conditions.(4)When selecting a pump for an axial-flow pumping station, the first saddle bottom head of the axial-flow pump or the saddle bottom head of the corresponding pump device should be taken as the control value of the highest head of the pumping station. The method makes up for the defect of the traditional method of choosing the control value and avoids the reduction of the operation range of the pumping station.

Abbreviations

D:Impeller diameter, mm
Dp:Pipe diameter, mm
Q:Flow, L/s
Qd:Design flow, L/s
n:Rotation speed, r/min
:Inflow velocity, m/s
p:Flow field pressure, Pa
:Fluid density, kg/m3
:Fluid dynamic viscosity, m2/s
Gk:The generation of turbulent kinetic energy
αk:Prandtl numbers for k
αε:Prandtl numbers for ε
C1ε, C1ε:Empirical constants
H:Head, m
:Total pressure of the outlet pressure measurement section, Pa
:Total pressure of the inlet pressure measurement section, Pa
:Local acceleration of gravity, m/s2
:Torque, N·m.

Data Availability

The data in this paper were obtained by physical experiment and numerical simulation, and 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 regarding the publication of this paper.

Authors’ Contributions

Dongtao Ji, Weigang Lu, and Lei Xu curated the data. Linguang Lu conducted formal analysis. Dongtao Ji wrote the original draft. Wei Shi, Jun Liu, and Ye Zhu reviewed and edited the article.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (grant no. 51779215) and Water Conservancy Science and Technology Project of Jiangsu Province (grant no. 2019015).