1. Introduction
The vortex pump evolved from a centrifugal pump with a semi-open impeller, and it is a type of non-clogging pump with an impeller located completely within the impeller chamber behind the vortex chamber, as shown in
Figure 1. Its unique structure makes it suitable for conveying solid–liquid mixtures containing solid and long fibers [
1,
2,
3], and in particular, it shows excellent performance in the hydraulic transportation of gases containing manure, biogas, digesters, waste liquid, and agricultural products such as potato, etc., whilst solid media were not damaged. At present, if we compare the vortex pump with the centrifugal pump, the biggest disadvantage is the low efficiency of the vortex pump [
4,
5]. In recent years, researchers have carried out studies on pump design and numerical simulation to improve the performance of the vortex pump. Gerlach et al. analyzed the how variations in geometric parameters influence vortex pump characteristics, and investigated many articles considering parameters in both single-parameter and multi-parameter experiments [
6,
7]. Moreover, Gerlach et al. performed a numerical study of the internal flow field of a vortex pump, and made a comparison with the Hamel–Oseen vortex model. The study suggests that a vortex similar to the Hamel–Oseen vortex is only present at the strong part load operation [
8]. Through a numerical simulation and experimental analysis, Jiang and Zhu studied the influence of different types of hems and blades on the performance of the vortex pump [
9,
10]. Li carried out a numerical simulation and experiment analysis for the gas–liquid two-phase vortex pump [
11]. Because of the change in structure from a centrifugal pump to a vortex pump, the internal flow of the vortex pump is complicated, and the flow through the outlet of the pump is the intersection of the flow through the vortex chamber and the impeller. The 32WB8-12 vortex pump was built based on continuous research findings from the 1980s. In order to design and manufacture a vortex pump with good performance and establish a method of optimum design, we must master the internal flow rules of the pump. So, this paper carries out research work on a performance test and internal flow field simulation of the vortex pump. At present, the internal flow velocity field of the vortex pump can be measured by the particle image velocimetry (PIV) technique, but the pressure field of the high-speed rotating impeller cannot be measured without contact. Therefore, the numerical calculation of the flow field is an effective way to obtain the internal flow rules. In order to determine whether the calculation results are correct and can reflect the actual objective, we need scientific experimental verification. The volute chamber of the vortex pump is relatively wide, so it is possible to measure the flow field by probe contact. As long as verification of the numerical calculation results of the vortex chamber flow field through comparative analysis reaches a certain level of precision, it can be extended to deduce that the numerical calculation of the internal flow field of the whole vortex pump gives realistic results. Our paper is based on this idea.
The flow field inside the vortex pump was simulated by ANSYS® CFX software, and the performance curve was calculated based on the design, development, and performance test of the model pump. The flow field in the vortex chamber was measured by a five-hole probe, and the correctness of the numerical results was verified. After a comprehensive analysis, the internal flow rules of the vortex pump and the constraint relationship between the pump performance and the flow characteristics were clarified.
2. Design and Test Scheme of the Model Pump
For the study of the vortex pump, considering factors such as the scale of the test platform, this paper designed and developed the 32WB8-12 pump (Zhejiang University of Science and Technology, Hangzhou, Zhejiang, China) model. The hydraulic design and geometric parameters are shown in
Table 1, and its structure and impeller, pressurized water chamber parts, coordinate system, and the position of the measuring hole are shown in
Figure 1.
3. Model Pump Performance
The pump is directly connected to the motor, and the power of the pump shaft was measured by the electronic method. The no-load motor is in line with GB/T 12785-91 standards (Test Methods for Submersible Motor Pumps), the pump performance test was carried out according to GB/T 3216-2005 (Rotodynamic Pumps—Hydraulic Performance Acceptance Tests, Grades 1 and 2), and the frequency converter was used to adjust to the rated speed.
Figure 2 is the scheme of the experimental test bench. The technical data for instrumentation are provided in
Table 2.
Pump experimental performance curves are shown in
Figure 3 (
ηgr is the unit efficiency). The optimal operating conditions of the pump are: discharge (
qv) is 9.00 m
3/h, pump head (
H) is 12.88 m, maximum efficiency (
ηmax) is 53.13%, and specific speed
ns is 76.5. The pump
qv-
H curve is flat, higher than in the design condition. The
qv-
P curve rises faster than other curves. The discharge-critical net positive suction head (
qv-
NPSHc) curve is parabolic-diminishing and shows the opposite trend compared with the centrifugal pump’s
qv-
NPSHc curve. The critical
NPSH of small flow is larger, but the critical
NPSHc decreases with the increased discharge.
ηp: pump efficiency; ηgr: unit efficiency; H: pump head in m; P: shaft power in W; NPSHc: critical net positive suction head in m; qv: discharge in m3/h; qv-H: discharge-pump head; qv-P: discharge-pump shaft power in w; qv-ηp: discharge-pump efficiency; qv-ηgr: discharge–unit efficiency; qv-NPSHc: discharge-critical net positive suction head.
4. Numerical Calculation of Pump’s Internal Flow Field
4.1. Governing Equation
The flow field was calculated based on selecting the best condition of pump performance test. The internal flow field of the vortex pump is very complex, which is three-dimensional turbulent flow dominated by strong swirling flow. It is affected by many factors, such as curvature, rotation, and so on. Choosing the turbulence model for the calculation of the internal flow field is very important. The previous calculation used the standard
k-
ε two-equation model [
12,
13,
14], which is a local equilibrium model that cannot accurately predict turbulence with drastic changes in mean flow, such as separation flow and swirling flow. Later, the researchers applied the modified
k-
ε and the Reynolds stress equation model—which was derived from the anisotropic conditions—to the numerical calculation of the strong cyclonic flow field. The renormalization group
k-
ε turbulence model is called the RNG
k-
ε model. The coefficient in the RNG
k-
ε turbulence model is derived by using the RNG model theory. There is an additional item R in the
k-
ε equation, which represents the effect of the average strain rate on the turbulent dissipation rate
ε. The RNG
k-
ε turbulence model can simulate the flow separation and vortex well. Moreover, because of the complexity of the Reynolds stress model and the large amount of computational resources, we chose the RNG
k-
ε model to numerically calculate the three-dimensional flow in the vortex pump [
15].
There are three expressions of the governing equation:
- (1)
- (2)
- (3)
The
k and
ε equation
where
where
is unit tensor,
,
,
,
,
,
, and
.
The environment temperature is 25 °C. The medium is water, and the density is 997 kg/m3. The dynamic viscosity is 0.89 mPa·s. The quantities with over-bar (such as ) are mean values.
4.2. 3D Modeling and Grid Division
The 3D modeling software Pro/ENGINEER (Version 5.0, PTC, Needham, MA, USA, 2015) [
16,
17,
18,
19] was used to perform solid modeling of the vortex chamber of the vortex pump and impeller region. The computational domain consists of the impeller region and volute region, in which the volute region is the stationary region. The impeller region is the rotating region and the transmission coupling flow parameter between the impeller and volute region is processed by the frozen rotor method. The calculation of the impeller flow field is carried out in the relative coordinate system, while the vortex chamber flow field calculation is carried out in the absolute coordinate system and the integral calculation of the dynamic and static components is realized through the coupling calculations between the two parts. Because the interface grid between the computational domains is not exactly the same, the general grid interfaces (GGIs) are used to connect the parts. The computational grids are divided by ANSYS
® ICEM Tetra tools (Version 14.5, ANSYS, Canonsburg, PA, USA, 2012). The tetrahedral mesh is divided into the vortex chamber region and the impeller region, respectively, and the residual error is 10
−4. The number of mesh units is 616,369, and the number of nodes is 117,300, as shown in
Figure 4.
4.3. Boundary Conditions
4.3.1. Inlet Boundary Conditions
The inlet of the vortex pump calculation area is an axial inlet. The inlet axial velocity value is given at the inlet pipe according to the flow rate. Assuming that the inlet has no rotation, the tangential velocity and the radial velocity are 0, and the turbulence is fully developed and evenly distributed in the import section. A is the area of the inlet of the vortex pump. The inlet axial velocity
is given at the inlet pipe according to the flow rate
.
The turbulent kinetic energy and turbulent energy dissipation rate at the inlet are set according to the flow and inlet cross-sectional area. The formula is as follows:
is the mean flow velocity, and here
=
.
where
I is turbulence intensity
From the above formula,
is the root-mean-square of the turbulent velocity fluctuations,
is the is the mean velocity, and
is the inlet’s length,
is the Reynolds number based on
DH, and
DH represents the hydraulic diameter,
where
is 0.09, and
L is the associated length. For the full development of the turbulence, desirable
L is equal to the hydraulzic diameter. Here,
L is equal to the inlet round pipe diameter.
Specific boundary conditions: velocity vector (0 m/s, 0 m/s, 3.11 m/s), turbulent kinetic energy = 0.021 m2/s2, and turbulent kinetic energy dissipation = 0.22 m2/s3.
4.3.2. Outlet Boundary Conditions
Assuming that the outlet flow is fully developed for turbulence, the specific boundary conditions are as follows:
The outlet type is known as Outlet, and the average static pressure ps = 101,325 Pa.
4.3.3. Wall Boundary Conditions
In the non-slip boundary condition (No Slip) of the solid wall, the wall function method is used to calculate the turbulence in the vicinity of the wall. The influence of the wall roughness on the flow field is neglected in the calculation.
5. Numerical Results
5.1. The Flow Distribution of the Static Pressure Inside the Pump
Figure 5,
Figure 6,
Figure 7 and
Figure 8 are all based on the calculation result of CFX software (Version 14.5, ANSYS, Canonsburg, PA, USA, 2012). Based on the coordinate system in
Figure 1,
Figure 5 shows the static pressure distribution on the axial cross-section where
Z is −10 mm. This axial section is also the middle cross-section of the impeller blade. According to the static pressure value and its position (radial direction and
),
Figure 6 is generated.
Figure 6 shows the static pressure
ps variable curves along the 0°, 90°, 180°, and 270° radial directions.
Figure 5 and
Figure 6 show that the center of the impeller is the lowest static pressure area in this axial cross-section, which is negative pressure (vacuum). The pressure in the center of the impeller is also lower than the pressure in the center of the chamber, so the run-through flow from chamber to impeller can be formed. When the fluid flow enters into the region near the blade, centrifugal force applies work on fluid flow and the pressure is increasing. So, in the region near the impeller edge, the pressure is higher than other regions.
Figure 7 shows the static pressure distribution for the axial cross-section where
Z is 12.5 mm, which is also based on the coordinate system in
Figure 1. This axial section is the middle axial cross-section of the chamber.
Figure 8 is generated according to the static pressure value and its position (radial direction).
Figure 8 shows the static pressure
ps variable curves along the 0°, 90°, 180°, and 270° radial directions.
Figure 7 and
Figure 8 show that the center of the chamber is the lowest pressure area in this axial cross-section, which is negative pressure (vacuum). So, the fluid flow from inlet to chamber can be formed.
5.2. Flow Velocity Distribution in the Pump
Figure 9 and
Figure 10 are also the calculation result of CFX software. Based on the coordinate system in
Figure 1,
Figure 9 and
Figure 10 respectively show the absolute velocity
, the circumferential velocity
, the radial velocity
, and the axial velocity
of two different axial cross-sections. Absolute velocity
is the vector synthesis result of
.
Figure 9 shows the velocity variable on the axial cross-section, where
Z is −10 mm. This axial section is the middle cross-section of the impeller blade.
Figure 9a–d show the
variable curves along the 0°, 90°, 180°, and 270° radial directions on the axial cross-section where
Z is −10 mm.
Figure 10 shows the velocity variable on the axial cross-section, where
Z is 12.5 mm. This axial section is the middle axial cross-section of the chamber.
Figure 10a–d show how the
change along the 0°, 90°, 180°, and 270° radial directions on the axial cross-section, where
Z is 12.5 mm.
5.3. Numerical Simulation of the Pump Performance Curve
The discharge calculation software CFX was used to calculate the discharge at 2850 r/min. The calculated performance and measured curves are shown in
Figure 11.
Disk frictional loss and volume loss are difficult to determine, and the computational fluid dynamics (CFD) numerical simulation pumper efficiency only considers hydraulic efficiency. The calculation model is simplified, for instance, by neglecting the clearance between the impeller back shroud, water pressure chamber, and surface roughness, etc.
6. Experimental Verification
In this paper, five spherical probes were used to extend into the vortex chamber, and five points were measured. The location of the measuring points is shown in
Figure 12. The measuring system and principle of the probe are shown in
Figure 11. The flow field calculation, measurement data, and the curve contrast are shown in
Table 3 and
Figure 13. Curves in
Figure 13 are based on the coordinate system in
Figure 1.
It can be seen from the comparison and analysis that the numerical calculation is consistent with the experimental results, which can reflect the internal flow pattern of the vortex pump. The numerical simulation results are reliable, and can be applied to the qualitative and quantitative analysis of the internal flow field of the vortex pump.
7. Comprehensive Analysis of Research Results
Based on the CFD numerical calculation and experimental results, we can carry out a comprehensive analysis of the internal flow rules of the vortex pump
Due to the unique structure and working principle of the vortex pump (i.e., the use of a semi-open impeller), the impeller and vortex chambers are fully open and connected, and the vortex chamber and pump outlet are connected through a channel. Carrying on with the flow field analysis, we must consider the concept of potential strength; i.e., the flow field of gravity and centrifugal inertia forces. The unequal mass force of the impeller is an important reason for the asymmetric distribution of the impeller flow field and the degree of disorder.
In the center of the impeller, because the hub occupies the central space, in order to form a vacuum range in the 10 mm < r < 20 mm annular area (r represents the radius with the center of the impeller as the origin), the lowest point of static pressure is in the upper y = 20 mm position, up to −3.6 × 104 Pa, and it can be inferred that bubbles are first produced in this region, while the lower static pressure remained at −2.7 × 104 Pa or so, forming clear dropped gradient upward pressure. The pressure in the lower half of the whole impeller chamber and the vortex chamber is greater than the pressure at the corresponding position of the upper part, and the static pressure reaches a maximum of about 1.1 × 105 Pa at the outlet of the pump body and becomes the main body of the pump head. The vacuum area is greater than the pump inlet area (pump inlet d = 32 mm). The main reason is the gravity field mass force potential function W = gz; the potential energy at the lower position is greater than that at the top. In addition, the lower gravity is in accordance with the centrifugal inertial force, increasing the pressure energy in the region. The structure is mainly due to the influence of the pump’s outlet.
In the range of 20 mm < r < 48 mm, the static pressure of the impeller chamber and the vortex chamber grows linearly, which is influenced by the impeller and has a forced vortex property. In the range of 50 mm < r < 70 mm, the upward trend of static pressure is slowed down, which indicates that the imposition of the impeller is reduced and the steady growth of the static pressure is due to the increase in the volume of the vortex chamber and the partial kinetic energy. The hydraulic design of the radius of the volute chamber and the width of the channel is calculated by this principle, and the cross-sectional area of the flow passage is found in the determination of the optimum velocity in the flow passage.
In the range of r > 20 mm, the center speed of vu and the absolute speed of the impeller almost coincide with each other, and have greater values than the speed at the center of the vortex chamber. That is, the velocity of the volute chamber lags behind the rotation velocity of the impeller, and the vortex chamber and almost coincide. It can be seen that the width of the volute chamber should not be designed too wide under the premise of ensuring passing capacity; otherwise, the delay of the chamber velocity will increase.
Within the radius of the impeller, is approximately equal to the circumferential velocity of the impeller with the same radius. The is approximately 14.3 m/s, which is equal to the circumferential velocity ( = 14.326 m/s) near the outer edge of the impeller. It can be shown that numerical simulation achieves sufficient accuracy. In the range of 20 mm < r < 48 mm, the center of and of the impeller rise almost linearly, and the center of the vortex chamber rises to approximately a parabola. In the range of r > 48 mm, the center of the vortex chamber is of approximately parabolic descent, and the partial kinetic energy is converted to pressure energy.
Because of the pressure gradient, the r < 20 mm range of axial velocity distribution is not uniform, and the maximum velocity is in the position of 90° section y = 15 mm. At a maximum of about ≈ 5.0 m/s for the impeller, the volute chamber value is ≈ 3.0 m/s. In the r < 20 mm range, the radial velocity varies irregularly, while in the r > 20 mm range, the radial velocity fluctuates little.
The calculated pump head is slightly higher than the measured value. This is due to the neglect of mechanical loss and volumetric loss. The maximum hydraulic efficiency of the vortex pump is about 60%. Low hydraulic efficiency (i.e., flow loss) causes the main part of pump energy loss [
20].
8. Conclusions
In a small flow state, the qv-NPSHc curve of vortex pump is the opposite to the qv-NPSHc of the centrifugal pump, and NPSHc value is large, that is, the vortex pump has poor anti-cavitation. At present, in ISO9906 (rotodynamic pumps—hydraulic performance acceptance tests—grades 1, 2 and 3), the evaluation methods of the centrifugal pumps are used for assessment of the rotodynamic pumps. Hence, the rationality of the evaluation methods is worthy of further study.
The experimental performance test of the prototype of the vortex pump (32WB8-12) indicates the 32WB8-12 pump has good performance and its optimized hydraulic model is worth popularizing. The numerical calculation is verified by the local experiment, and it is concluded that the numerical calculation of the internal flow field of the vortex pump gives realistic results. The pressure distribution in the vortex chamber presents a uniform increasing trend, and the velocity distribution appears to be in a turbulent state, which is characterized by a combined vortex. It is suggested that the flow rate be controlled to improve the efficiency of the vortex pump [
21].
The stable pressure gradient is formed from the inlet of the vortex pump to the center of the impeller. The lowest pressure in the vortex chamber is −2.5 × 104 Pa and it is located at the center of the cavity of vortex chamber. The area of low pressure is greater than the pump entrance area. Therefore, the 32WB8-12 vortex pump has advanced suction and anti-cavitation performance.
Probes were used for measuring the pressure and velocity field synchronously, which had to be installed in the flow and caused slight interference to the original flow field and generated some error. In the following research, non-contact laser PIV technology will be applied to measure velocity, as it can obtain the velocity field accurately without disturbing the flow field. Non-contact laser PIV can be combined with the probe to measure the accurate internal flow field of the vortex pump. Furthermore, we will establish a flow model based on studies of precise numerical simulation technology and carry out studies on the two-phase flow medium and cavitation fluid field of the vortex pump. Thus, we can reveal the internal mechanism of the excellent performance of the vortex pump.