Advances in Materials Science and Engineering
Volume 2020 (2020), Article ID 6665520, 17 pages
https://doi.org/10.1155/2020/6665520

Research Article

Influence of Blade Profiles on Plastic Centrifugal Pump Performance

Hui Zhang, Lingfeng Tang, and Yongtao Zhao

School of Mechanical Engineering, Anhui Polytechnic University, Wuhu, Anhui 241000, China

Correspondence should be addressed to Hui Zhang; 2190110111@stu.ahpu.edu.cn

Received 9 October 2020; Revised 26 November 2020; Accepted 2 December 2020; Published 22 December 2020

Academic Editor: Jinyang Xu

Copyright © 2020 Hui Zhang et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Abstract.

To study the influence of blade profiles of the plastic centrifugal pump on pump performance, the impeller blade profiles were designed and drawn by the single arc method, double arc method, logarithmic spiral method, and B-spline curve method, respectively, with the known structural parameters.The shape and size of four profiles were drawn, and two-dimensional models and three-dimensional models of four impellers and volute were completed, respectively. The impeller models were printed by 3D printing technology, and the performance experiments of the plastic centrifugal pump were carried out. The numerical simulation of the internal flow field was performed. From the contours of the velocity and pressure, the vapor volume fraction distribution, and fluid-structure interaction analysis of impellers, the impeller drawn by the logarithmic spiral method was better than others. The optimization of the logarithmic spiral method was completed. The impeller inlet and outlet diameters (D₁ and D₂) and impeller inlet and outlet installation angles (β₁ and β₂) were taken as control variables, and the total power loss and the minimum NPSHr of the pump were taken as the objective functions. The optimization results were that D₁ = 58 mm and D₂ = 162 mm and β₁ = 17° and β₂ = 31°. The hydraulic efficiency was increased by 1.68%.

1. Introduction

As a type of general machinery, centrifugal pumps are widely used such as in mechanical engineering, aerospace, and petrochemical industries [1–3]. The flow components of plastic centrifugal pumps are made of engineering plastic, with less manufacturing costs and better corrosion resistance. Plastic centrifugal pumps are suitable for most corrosive environments such as acid, alkali, and salt and can be applied to the transportation of corrosive materials in chemical and petroleum industries. However, it is well known that inevitable cavitation in centrifugal pumps will affect the pump performance and cause noise and vibrations [4–10].

The impeller is the core part of the centrifugal pump, and the impeller blade profile plays a vital role in the fluid flow which will directly impact the performance of centrifugal pumps. At present, there have been many studies on the influence of the impeller parameters on the flow characteristics of centrifugal pumps [11–20]. Tao et al. [21] used the numerical method to study how the blade thickness influenced the flow characteristics in the impeller passages and volute casing. They concluded that the pressure pulsations increased at the leading edges but decreased at the trailing edges according to the increase in the blade thickness. However, their work was only helpful for the design of ceramic pumps. Tao et al. [22] studied the influence of blade leading-edge shape on the cavitation and concluded it delayed the inception cavitation and critical cavitation. Peng et al. [23] analyzed various blade outlet angles without changing other impeller parameters and concluded that if the blade outlet angle was increased, the flow separation and the axial vortex along the blade working surface would be more serious. However, their work was aimed at low specific speed submersible pumps. Yu-Qin and Ze-Wen [24] studied how different blade numbers influenced the flow-induced noise of the centrifugal pump based on the idea of acoustic-vibration coupling. They found that the 6-blade model pump was the optimal value. Han et al. [25] presented an impeller and volute shape optimization design method that improved the hydraulic performance of the pump. Zhang et al. and Wang et al. [26, 27] explored the effects of slotted blades on the performance of centrifugal pumps. Li et al. [28] studied the influence of the blade angle of a low specific speed plastic centrifugal pump on the pump performance. They applied the one-way fluid-structure interaction but did not take the impact of structural deformation on performance into consideration.

The blade profile directly affects the bending of the blade and has an important impact on the fluid flow in the impeller channel. Studying the impact of the blade profile on the performance of the centrifugal pumps is beneficial to improve the performance and efficiency of the centrifugal pumps. Hu et al. [29] designed an iterative method of double arc blade profile for drawing cylindrical blade profile which overcame the shortcomings of blade profile drawing methods. El-Gazzar and Hawash [30] used the splitting technique to change the shape of the impeller and concluded that the impeller blade split caused the uneven mass distribution around the impeller, which led to the imbalance of the impeller. Ming et al. [31] proposed a new variable-angle spiral equation that was complementary to the traditional variable-angle spiral equation by changing the blade placement angle function.

In existing research studies, the studies on blade profile are mainly focused on metal pumps and a single profile, while only a few pieces of research are on the performance of plastic centrifugal pumps. Moreover, there is no regularity between the various blade profiles and the performance parameters of centrifugal pumps. In this paper, four design methods of blade profiles, which are the single arc method, double arc method, logarithmic spiral method, and B-spline curve method, are proposed and the influence of four kinds of blade profile on plastic centrifugal pump performance is studied. The flow field simulation of pumps with various impellers is analyzed in detail.

2. Design Methods of Blade Profiles

2.1. Single Arc Method

Figure 1 shows the design principle of the single arc method. In Figure 1, D₁ and D₂ are impeller inlet and outlet diameters, respectively, β₁ and β₂ are the impeller inlet and outlet installation angles, respectively, and arc AB is the blade profile.

Figure 1: Blade profile of single arc method.

Take a random point (point A) on the entrance circle D₁. Take point A as the fixed point. Draw a line with a β₁ angle and a length of ρ in the negative direction of the x-axis. Take point E as the endpoint. Then, take point E as the center and ρ as the radius to make the arc intersect the inlet and outlet circles at A and B. Therefore, the arc AB is the blade profile of the single arc method.

The radius of curvature of the single arc:

where ρ = EA—radius of curvature of the single arc (mm) and R₁, R₂—impeller inlet and outlet radii (mm).

2.2. Double Arc Method

The principle of the double arc method is to find the radius of rotation according to the drawing method to determine the first arc and then calculate the value of the radius of rotation to determine the second arc according to the formula.

In Figure 2, D₁ and D₂ are impeller inlet and outlet diameters and β₁ and β₂ are the impeller inlet and outlet installation angles. Take D₁ as the diameter to make a circle and divide the circle into equal parts according to the number of blades Z. The equal points are A₁, A₂, A₃, A₄, A₅, and A₆. Take the diameter δ = D₁sin β₁ as the concentric circle with the diameter D₁. Make tangent lines A₁E₁, A₂E₂, A₃E₃, A₄E₄, A₅E₅, and A₆E₆ from each equal point to the circle of diameter δ. Take E₁ as the center of the circle (∠E₁A₁O = β₁) and E₁A₁ as the radius to make an arc, which intersects the auxiliary circle (point F); then, the arc FA₁ is the first profile of the blade. The second arc needs to meet three conditions which determine the position of point P. Three conditions are that (1) F, E1, and M points are collinear; (2) ∠OPM = β₂; and (3) PM = MF = ρ.

Figure 2: Blade profile of double arc method.

The value of ρ is

where ρ—radius of rotation of the second arc (mm); Rf—radius of rotation of the point F (mm); and βf—blade angle of point F (°).

After the position of point P is determined, the center point M can be determined and the second arc can be determined according to ρ.

2.3. Logarithmic Spiral Method

After determining the impeller inlet and outlet installation angles and the impeller inlet and out diameters, the blade wrap angle has a larger value range. Therefore, the advantage of the logarithmic spiral method in the shape of the cylindrical impeller’s blade is also more prominent. Design and draw the impeller blade profile according to the equivariant angle spiral, and its design principle of the profile is shown in Figure 3.

Figure 3: Blade profile of logarithmic spiral method.

In Figure 3, β is impeller installation angle (°); β₁ and β₂ are impeller inlet and outlet installation angles (°); D₁ and D₂ are the impeller inlet and outlet diameters (mm); and θ is the variation of wrap angle (°).

When β₁ ≠ β₂, the change law of β is

Profile equation:

Wrap angle:

Substitute β₁ = 20, β₂ = 30, D₁ = 58 mm, and D₂ = 162 mm into the program written by MATLAB software, and the resulting blade profile is shown in Figure 4(a).

(a)

(b)

Figure 4: Blade profiles drawn by MATLAB software. (a) Logarithmic spiral method. (b) B-spline curve method.

2.4. B-Spline Curve Method

Use 5-point (or 4-point) Bezier curve to fit the studied impeller profile, which eliminates the inconvenience of curve adjustment caused by traditional cubic polynomial fitting and makes the design of blade profile simpler, smoother, and easier to control. Determine the coordinates of the control points first. According to the Bezier equation, n = 4,

The coordinate of any point P (xp, yp) on the Bezier curve is

Figure 5 shows the design principle of the B-spline curve method. In Figure 5, the impeller inlet and outlet installation angles β₁ and β₂ and the impeller inlet and outlet diameters D₁ and D₂ are the known parameters.

Figure 5: Blade profile of B-spline curve method.

Substitute β₁ = 20, β₂ = 30, D₁ = 58 mm, D₂ = 162 mm, and ψ = 88.3 into the calculation and obtain P₀ (0, 29), P₁ (−26.32, 38.58), P₂ (60.87, 12.77), P₃ (−66.79, 25.49), and P₄(−81, 2.83). From the obtained coordinates of P₀, P₁, P₂, P₃, and P₄, the designed blade profile can be drawn by MATLAB software, as shown in Figure 4(b).

3. Numerical Simulation Method

3.1. Computation Model

Take the plastic centrifugal pump designed by a certain unit as an example, and its main design parameters are shown in Table 1. According to the known structural parameters of a plastic centrifugal pump, the impeller blade profiles were designed and drawn by the single arc method, double arc method, logarithmic spiral method, and B-spline curve method, respectively.

Table 1: Design parameters of plastic centrifugal pump.

Performance parameters
Geometrical dimension

Q (m³/h)
20
Impeller inlet diameter (mm)
58

H (m)
30
Impeller outlet diameter (mm)
162

n (r/min)
2900
Impeller inlet width (mm)
23

ns
60.871
Impeller outlet width (mm)
10

Pt (kW)
5.5
Impeller inlet installation angle (°)
20

NPSHa (m)
4
Impeller outlet installation angle (°)
30

Volute width (mm)
33
Base circle diameter of volute (mm)
175

Figure 6 shows a two-dimensional diagram of the impeller. According to the parameters of the impeller in Figure 6, the UG software is used to design the three-dimensional model of the impeller. The three-dimensional models of the impeller with four profile design methods are shown in Figure 7.

Figure 6: Two-dimensional diagram of the impeller.

(a)

(b)

(c)

(d)

Figure 7: Models of impellers drawn by four profile design methods. (a) Single arc. (b) Double arc. (c) Logarithmic spiral. (d) B-spline curve.

Figure 8 shows a two-dimensional diagram of the volute. According to the parameters of the volute in Figure 8, the UG software is used to design the three-dimensional model of the volute. The three-dimensional model of the volute is shown in Figure 9.

Figure 8: Two-dimensional diagram of the volute.

Figure 9: Model of volute.

3.2. Meshing

According to the above completed three-dimensional models, ANSYS-ICEM software is employed to utilized to mesh the computational domains. When we mesh the models, the actual grid cannot reach the ideal shape. If the grid deforms or deformation exceeds a certain limit, the accuracy of the calculation results will change accordingly. Therefore, in the initial division of the grid, it is necessary to use appropriate measures to control or measure the quality of the grid and try to achieve the best grid. The final determination of the number of grid cells in each part of the computational domains is shown in Table 2.

Table 2: The number of grid cells.

Computational domains
Inlet extension
Impeller
Volute
Outlet extension

The number of grid cells
581102
397526
393266
1502055

The final meshing is shown in Figure 10.

Figure 10: Meshing.

3.3. Computational Setup

Numerical calculations were performed in FLUENT software. The turbulence model was a k-epsilon (equation (2)) turbulence model. The computational impeller domain was set to a rotational domain. All other computational domains were set to static.

The boundary conditions were set to velocity inlet and outflow. The relevant setting parameters were set by default. The reference pressure was set to standard atmospheric pressure. The wall surface was placed under a non-slip boundary condition, and a standard wall surface function was applied. The calculation method was SIMPLEC. The convergence accuracy was set to 10⁻⁴.

4. Test Verification

4.1. 3D Printing of Plastic Centrifugal Pump Impellers

UTR9000 material was used for 3D printing of impellers, which was an ABS-like stereo light modeling resin with accurate and durable characteristics. The durability of components made of UTR9000 resin was more than 6.5 months. According to the three-dimensional models of the impellers drawn by four profile design methods (Figure 6), the impeller models printed by 3D printing technology are shown in Figure 11.

(a)

(b)

(c)

(d)

Figure 11: 3D printing of impellers. (a) Single arc. (b) Double arc. (c) Logarithmic spiral. (d) B-spline curve.

4.2. Performance Tests of Plastic Centrifugal Pump

Figure 12 shows the pump performance test bench, which mainly included a motor and a flow control valve.

Figure 12: Pump performance test bench.

The pump test system (electrical measurement method) V8.97 software on the data acquisition platform was used to collect data at each operating point by adjusting the flow control valve. The pump performance experiments were performed on the impellers of the four types of profiles. The parameter settings before the experiment are shown in Table 3, and the collected data are shown in Tables 4–7.

Table 3: Parameter settings of pump performance experiment.

Product number
Fluid density (kg/m³)
Fluid viscosity (CST)
Meter factor (f/1)
Inlet diameter (m)
Outlet diameter (m)

50UHB15-32
Water
Water
79.89
0.05
0.032

Table 4: Performance experiment of single arc method.

Measured value
n = 2900 (r/min)

Inlet pressure
Outlet pressure
Q
n
Pin
Q
H
Pin
η

(kPa)
(kPa)
(m³/h)
(r/min)
(kW)
(m³/h)
(m)
(kW)
(%)

1
−70.00
170.00
20.95
2998.60
5.85
20.26
25.72
5.29
26.83

2
−80.00
220.00
19.93
3001.44
5.83
19.26
31.19
5.26
31.11

3
−70.00
240.00
18.92
3002.21
5.81
18.28
31.94
5.24
30.37

4
−60.00
260.00
16.35
3001.81
5.56
15.80
32.47
5.01
27.87

5
−50.00
280.00
14.86
3001.81
5.37
14.36
33.21
4.84
26.82

6
−40.00
300.00
12.93
3002.06
5.12
12.49
33.90
4.62
24.99

7
−30.00
310.00
10.91
3002.36
4.96
10.54
33.66
4.47
21.62

8
−20.00
330.00
8.28
3001.63
4.68
8.00
34.39
4.22
17.76

9
−20.00
340.00
6.33
3001.54
4.53
6.12
35.21
4.09
14.36

10
−10.00
350.00
2.44
3001.34
4.26
2.36
35.06
3.84
5.86

11
−10.00
360.00
1.49
3001.16
4.22
1.44
36.00
3.81
3.71

12
−10.00
360.00
0.21
3000.81
4.15
0.20
35.99
3.75
0.53

13
0.00
380.00
0.00
3000.90
4.03
0.00
36.94
3.64
0.00

Table 5: Performance experiment of double arc method.

Measured value
n = 2900 (r/min)

Inlet pressure
Outlet pressure
Q
n
Pin
Q
H
Pin
η

(kPa)
(kPa)
(m³/h)
(r/min)
(kW)
(m³/h)
(m)
(kW)
(%)

1
−80.00
190.00
21.19
3000.59
5.47
20.48
28.60
4.94
32.31

2
−80.00
220.00
19.98
2998.58
5.33
19.32
31.26
4.82
34.13

3
−70.00
240.00
18.09
2997.92
5.21
17.50
31.89
4.72
32.23

4
−70.00
250.00
16.89
2998.15
5.18
16.34
32.64
4.69
30.98

5
−50.00
280.00
14.25
2998.20
5.01
13.78
33.20
4.53
27.50

6
−40.00
290.00
12.63
2998.18
4.82
12.22
33.00
4.36
25.17

7
−30.00
310.00
10.85
2997.94
4.60
10.50
33.76
4.16
23.18

8
−20.00
330.00
8.56
2998.14
4.44
8.28
34.50
4.02
19.36

9
−20.00
340.00
6.70
2998.48
4.32
6.48
35.31
3.91
15.95

10
−10.00
350.00
4.83
2998.58
4.20
4.67
35.20
3.80
11.79

11
−10.00
360.00
2.06
2998.72
3.99
1.99
36.06
3.61
5.42

12
−10.00
370.00
0.24
2999.02
3.82
0.23
36.99
3.45
0.68

13
−10.00
380.00
0.00
2999.52
3.71
0.00
37.93
3.35
0.00

Table 6: Performance experiment of logarithmic spiral method.

Measured value
n = 2900 (r/min)

Inlet pressure
Outlet pressure
Q
n
Pin
Q
H
Pin
η

(kPa)
(kPa)
(m³/h)
(r/min)
(kW)
(m³/h)
(m)
(kW)
(%)

1
−70.00
180.00
21.33
2998.60
5.58
20.63
26.75
5.05
29.79

2
−80.00
190.00
21.08
2998.07
5.36
20.39
28.62
4.85
32.77

3
−80.00
210.00
20.24
3001.45
5.21
19.56
30.30
4.70
34.35

4
−80.00
230.00
18.98
3002.29
5.20
18.33
31.95
4.69
34.05

5
−60.00
250.00
16.50
3001.64
5.01
15.94
31.55
4.52
33.32

6
−50.00
280.00
14.65
3001.62
4.80
14.15
33.18
4.33
29.55

7
−30.00
310.00
11.30
3001.65
4.48
10.92
33.72
4.04
24.82

8
−30.00
320.00
1.54
3001.62
4.42
10.18
34.59
3.99
24.08

9
−20.00
330.00
8.46
3001.77
4.20
8.17
34.40
3.79
20.23

10
−20.00
340.00
7.35
3001.72
4.11
7.10
35.27
3.71
18.41

11
−20.00
340.00
5.93
3001.50
3.98
5.73
35.19
3.59
15.30

12
−10.00
350.00
2.63
3001.25
3.74
2.54
35.06
3.37
7.19

13
−10.00
370.00
0.00
3001.34
3.43
0.00
36.93
3.09
0.00

Table 7: Performance experiment of B-spline curve method.

Measured value
n = 2900 (r/min)

Inlet pressure
Outlet pressure
Q
n
Pin
Q
H
Pin
η

(kPa)
(kPa)
(m³/h)
(r/min)
(kW)
(m³/h)
(m)
(kW)
(%)

1
−80.00
180.00
21.27
2998.15
5.40
20.57
27.70
4.89
31.77

2
−80.00
220.00
20.05
3001.80
5.31
19.37
31.21
4.79
34.39

3
−80.00
240.00
17.69
3001.60
5.13
17.09
32.69
4.63
32.90

4
−80.00
250.00
16.65
3001.51
5.03
16.09
33.48
4.54
32.34

5
−60.00
280.00
14.91
3001.36
4.81
14.41
34.17
4.34
30.91

6
−40.00
300.00
12.75
3001.13
4.60
12.32
33.90
4.15
27.41

7
−30.00
310.00
11.08
3000.93
4.44
10.71
33.71
4.01
24.54

8
−30.00
330.00
8.29
3000.63
4.26
8.01
35.37
3.85
20.07

9
−20.00
340.00
7.06
3000.61
4.17
6.82
35.28
3.76
17.42

10
−20.00
350.00
5.81
3000.56
4.05
5.62
36.16
3.66
15.13

11
−10.00
350.00
4.75
3000.45
3.92
4.59
35.16
3.54
12.42

12
−10.00
360.00
3.29
3000.55
3.82
3.18
36.05
3.45
9.05

13
−10.00
370.00
2.14
3000.56
3.72
2.07
36.97
3.36
6.20

4.3. Data Analysis

Figures 13 and 14 show the comparison of head and efficiency in performance experiments of four types of impellers. From the shut-off point to the maximum flow condition during the test, 13 flow condition points were selected to completely reflect the external characteristic variation curve. It can be seen that the efficiency of the logarithmic spiral method was better than that of the other three methods, and the efficiency of the single arc method was the worst.

Figure 13: Q-H curves of four profiles.

Figure 14: Q-η curves four profiles.

5. Analysis of Numerical Simulation Results

5.1. Velocity in the Flow Field

Figure 15 shows the contours of the velocity of the four impellers, respectively. As shown in Figure 15, when the other parameters were consistent, the impellers drawn by four profile design methods had insufficient fluid flow, and the internal flow fields in the impellers drawn by the single arc method and B-spline curve method had obvious vortex. The internal flow fields in the impellers drawn by the double arc method and logarithmic spiral method were significantly better than those drawn by the B-spline curve method and single arc method.

(a)

(b)

(c)

(d)

Figure 15: Contours of the velocity of the four impellers. (a) Single arc. (b) Double arc. (c) Logarithmic spiral. (d) B-spline curve.

5.2. Pressure in the Flow Field

Figure 16 shows the contours of the pressure of the four impellers, respectively. As shown in Figure 16, when the other parameters were consistent, the internal pressure trends of the impellers drawn by four design methods were correct, and the stratification was obvious. There was negative pressure near the impeller inlet in the flow fields of impellers drawn by the single arc method and B-spline method, while negative pressure did not occur in the other two flow fields of impellers. The maximum pressure values corresponding to the four methods were 437563 Pa, 441628 Pa, 438368 Pa, and 435518 Pa. The pressure fields in the impellers drawn by the double arc method and logarithmic spiral method were significantly better.

(a)

(b)

(c)

(d)

Figure 16: Contours of the pressure of the four impellers. (a) Single arc. (b) Double arc. (c) Logarithmic spiral. (d) B-spline curve.

5.3. Cavitation in the Flow Field

Figure 17 shows the vapor volume fraction distribution in four impellers when NPSHr was 1 m. As shown in Figure 17, when NPSHr was 1 m, the bubbles almost filled the entire impellers. The vapor volume fractions of the impellers drawn by the double arc method and logarithmic spiral method were 0.9455 and 0.9449, which were significantly lower than those in the other two impellers. The impeller drawn by the single arc method had the largest value.

(a)

(b)

(c)

(d)

Figure 17: Vapor volume fraction distribution (NPSHr = 1 m). (a) Single arc. (b) Double arc. (c) Logarithmic spiral. (d) B-spline curve.

Figure 18 shows the vapor volume fraction distribution in four impellers when NPSHr was 2 m. As shown in Figure 18, when NPSHr was 2 m, bubbles filled one third of the entire impellers. The vapor volume fractions of the impellers drawn by the double arc method and logarithmic spiral method were 0.861 and 0.8174, which were significantly lower than those in the other two impellers. The impeller drawn by the B-spline curve had the largest value.

(a)

(b)

(c)

(d)

Figure 18: Vapor volume fraction distribution (NPSHr = 2 m). (a) Single arc. (b) Double arc. (c) Logarithmic spiral. (d) B-spline curve.

In summary, when NPSHr = 1 m and 2 m, the impeller drawn by the logarithmic spiral method is significantly better than the other three impellers, which means that the logarithmic spiral method is the best design method.

5.4. Analysis of Efficiency

According to the above conclusion, CFD-Post software was used to calculate the performance parameters of the four profiles, as shown in Table 8. From Table 8, it can be concluded that the impeller of the logarithmic spiral method was the optimal impeller under the design condition.

Table 8: Performance parameters comparison of four methods.

Method
Single arc
Double arc
Logarithmic spiral
B-spline curve

H (m)
35.4435
35.4653
35.541
35.3835

P (kW)
3.50485
3.3662
3.30089
3.51129

η (%)
76.99
81.53
81.98
76.72

5.5. Analysis of Blade Profiles Based on Fluid-Structure Interaction

The purpose of fluid-structure interaction was to find how the flow field influenced the impeller, which was mainly reflected in the deformation of the impeller caused by the fluid, and the force on the impeller. ANSYS-CFX software was used to complete the fluid-structure interaction analysis.

Figure 19 shows the total deformation of the impellers drawn by four profile design methods under the same working condition. As a whole, the total deformation distribution of the four impellers was symmetrical and stratified obviously. The deformation of the impellers increased continuously with the increase of the radius, the maximum deformation occurred at the edge of the impellers, and the minimum deformation occurred at the inlet of the impellers.

(a)

(b)

(c)

(d)

Figure 19: Total deformation of the impellers. (a) Single arc. (b) Double arc. (c) Logarithmic spiral. (d) B-spline curve.

The maximum total deformation of the impellers drawn by four profile design methods was 0.61503 mm, 0.58614 mm, 0.57442 mm, and 0.61609 mm, respectively. The following can be concluded for total deformation: B-spline curve method > single arc method > double arc method > logarithmic spiral method.

Figure 20 shows the equivalent stress of the impellers drawn by four profile design methods under the same working conditions. As a whole, the impeller stress was mainly concentrated on the blade root, where the equivalent stress was significantly greater than that in other areas. Therefore, compared to other positions, the blade root near the impeller inlet was more prone to stress concentration.

(a)

(b)

(c)

(d)

Figure 20: Equivalent stress of the impellers. (a) Single arc. (b) Double arc. (c) Logarithmic spiral. (d) B-spline curve.

The maximum equivalent stress of the impellers drawn by four profiles was 41.357 MPa, 41.43 MPa, 40.909 MPa, and 44.891 MPa, respectively. The following can be concluded for equivalent stress: B-spline curve method > double arc method > single arc method > logarithmic spiral method.

6. Optimization of the Logarithmic Spiral Method

6.1. Objective Function

Based on the above analysis, it is shown that the logarithmic spiral method has the highest efficiency. To maximize the efficiency, the impeller inlet and outlet diameters (D₁ and D₂) and impeller inlet and outlet installation angles (β₁ and β₂) were taken as control variables and the total power loss and the minimum NPSHr of the pump were taken as the objective functions to optimize the impeller structure of the logarithmic spiral method.

6.1.1. Minimum NPSHr

The basic equation of NPSHr:

Assuming that there is no prerotation at impeller inlet,

where k₁ = 0.876; k₂ = 0.91; Dj—impeller inlet diameter (m); dh—impeller hub diameter (m); ; λ—blade inlet pressure drop coefficient, λ = 0.15–0.3; —the absolute velocity of the slightly forward part of the blade inlet (m/s); and ω₀—the relative velocity of the slightly forward part of the blade inlet (°/s).

Then,

Therefore, the minimum objective function of NPSHr is min [f₁(x) = NPSHr].

6.1.2. Minimum Power Loss

The centrifugal pump losses mainly include volume, mechanical, and hydraulic losses. The volume loss accounts for a small proportion of the total loss, so the volume loss can be ignored. The optimization model is established with the mechanical loss and hydraulic loss.

where P—total loss power (kW); Pm1—hydraulic loss power (kW); and Pm2—mechanical loss power (kW).

where t₁—velocity energy loss coefficient in the pressurized water chamber, t₁ = 0.15∼0.25; t₂—dimensionless disc friction coefficient, t₂ = 0.8∼1.0; v₂—the absolute velocity of impeller outlet (m/s); ω—the rotational angular velocity of the shaft (°/s); and R₂—impeller radius (m).

Assuming that there is no prerotation at impeller inlet, u₁ = 0. According to Stodola’s equation,

Then,

where u₁, u₂—the circumferential velocity of inlet and outlet (m/s); , —the circumferential absolute velocity of inlet and outlet (m/s); Ht—theoretical head of the pump (m); n—rotational speed (rad/s); ψ₂—crowding coefficient of impeller outlet, ψ₂ = 0.842; and σ—Stodola’s slip coefficient.

Impeller outlet speed triangle:

Then,

Therefore, the minimum objective function of P is min[f₂ (x) = P].

6.1.3. Unified Objective Function

The weighting coefficient transformation method is used to deal with the two objective functions. According to the actual situation, the power loss of the plastic centrifugal pump has a greater impact on the total efficiency of the pump than NPSHr. Therefore, it is assumed that NPSHr accounts for 30% and the power loss accounts for 70% for optimization and the unified objective function is min [f(x) = 0.3P + 0.7NPSHr].

6.2. Constraint Conditions

The constraint ranges of design variables were determined by combining the statistical formula of the velocity coefficient method with the actual excellent model pump. Then, we increased the constraint ranges appropriately.

Constraint formulas:

6.3. Multiobjective Function Optimization

According to the creation of a multiobjective optimization program, the four parameters of impeller inlet and outlet diameters and impeller inlet and outlet installation angles were optimized. The final optimization results were that the impeller inlet and outlet diameters were 57.778 mm and 161.57 mm and the impeller inlet and outlet installation angles were 16.836° and 30.68°. The optimization results were rounded and compared with those before optimization, and the design variables are shown in Table 9.

Table 9: Comparison of optimization design variables.

Parameters
D₁ (mm)
D₂ (mm)
β₁ (°)
β₂ (°)

Before optimization
58
162
20
30

After optimization
58
162
17
31

6.4. Simulation Comparison and Verification

The optimized parameters were used to remodel the impeller. Under the same working conditions, the internal flow field simulation was performed. The results after analysis are shown in Figure 21. It can be seen that after optimization, the flow in the impeller was more sufficient. The efficiency was 83.66% after calculation, which was 1.68% higher than before. Therefore, the performance of the optimized impeller had been improved.

(a)

(b)

Figure 21: Optimization results.

7. Conclusions

(1)According to the basic parameters of the centrifugal pump and the four different design methods of impeller profile, the shapes and sizes of several profiles were determined. Three-dimensional models were drawn by UG software. ANSYS-ICEM software was utilized to mesh the computational domains.(2)According to the three-dimensional models drawn by four profile design methods, the impeller models were printed by 3D printing technology. The performance experiments of the plastic centrifugal pumps were carried out. The experimental results showed that the logarithm spiral method was the best and the single arc method was the worst.(3)The analysis of numerical simulation was performed by FLUENT software. From the contours of velocity and pressure, the flow field in the impeller drawn by the logarithmic spiral method was significantly better. The maximum pressure values corresponding to the four methods were 437563 Pa, 441628 Pa, 438368 Pa, and 435518 Pa. From the vapor volume fraction distribution, when NPSHr = 1 m and 2 m, the impeller drawn by the logarithmic spiral method was significantly better than the other three impellers. According to the calculation of performance parameters, it can be concluded that the impeller drawn by the logarithmic spiral method was the optimal impeller under the design condition.(4)ANSYS-CFX software was used to complete the fluid-structure interaction analysis. The maximum total deformation of the impellers drawn by the four profile design methods was 0.61503 mm, 0.58614 mm, 0.57442 mm, and 0.61609 mm, respectively. Deformation: B-spline curve method > single arc method > double arc method > logarithmic spiral method. The maximum equivalent stress of the impellers drawn by the four profile design methods was 41.357 MPa, 41.43 MPa, 40.909 MPa, and 44.891 MPa, respectively. Equivalent stress: B-spline curve method > double arc method > single arc method > logarithmic spiral method.(5)Based on the analysis, the impeller drawn by the logarithmic spiral method had the highest efficiency. To maximize the efficiency, the impeller inlet and outlet diameters (D₁ and D₂) and impeller inlet and outlet installation angles (β₁ and β₂) were taken as control variables and the total power loss and the minimum NPSHr of the pump were taken as the objective functions. The final optimization results were that the impeller inlet and outlet diameters were 58 mm and 162 mm and the impeller inlet and outlet installation angles were 17°and 31°. The hydraulic efficiency was increased by 1.68% after calculation, which indicated that the impeller structure had been improved.

Abbreviations

D₁: | Impeller inlet diameter, mm
D₂: | Impeller outlet diameter, mm
β₁: | Impeller inlet installation angle, °
β₂: | Impeller outlet installation angle, °
ρ: | Radius of curvature, mm
R₁: | Impeller inlet radius, mm
R₂: | Impeller outlet radius, mm
Rf: | Radius of rotation of the point F, mm
βf: | Blade angle of point F, °
Z: | Number of blades
β: | Impeller installation angle, °
: | Wrap angle, °
θ: | Variation of wrap angle , °
Q: | Flow rate, m³/h
H: | Pump head, m
n: | Rotational speed, r/min
ns: | Specific speed
Pt: | Shaft power, kW
Dj: | Impeller inlet diameter, m
dh: | Impeller hub diameter, m
λ: | Blade inlet pressure drop coefficient
: | Absolute velocity of the slightly forward part of the blade inlet, m/s
ω₀: | Relative velocity of the slightly forward part of the blade inlet, °/s
P: | Total loss power, kW
Pm1: | Hydraulic loss power, kW
Pm2: | Mechanical loss power, kW
t₁: | Velocity energy loss coefficient in the pressurized water chamber
t₂: | Dimensionless disc friction coefficient
: | Absolute velocity of impeller outlet, m/s
u₁: | Circumferential velocity of inlet, m/s
u₂: | Circumferential velocity of outlet, m/s
: | Circumferential absolute velocity of inlet, m/s
: | Circumferential absolute velocity of outlet, m/s
Ht: | Theoretical head of the pump, m
ψ₂: | Crowding coefficient of impeller outlet
σ: | Stodola’s slip coefficient.

Data Availability

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

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This work was supported by the University Synergy Innovation Program of Anhui Province (GXXT-2019-004), the Teaching Research Project of Anhui Education Department (2019jyxm0229), and the Natural Science Foundation of Anhui Province Education Department (KJ2017ZD12).

References

M. Tan, Y. Lu, X. Wu, H. Liu, and X. Tian, “Investigation on performance of a centrifugal pump with multi-malfunction,” Journal of Low Frequency Noise, Vibration and Active Control, Article ID 146134842094234, 2020.

M. Liu, L. Tan, and S. Cao, “Theoretical model of energy performance prediction and BEP determination for centrifugal pump as turbine,” Energy, vol. 172, pp. 712–732, 2019.

X. Li, B. Chen, X. Luo, and Z. Zuchao, “Effects of flow pattern on hydraulic performance and energy conversion characterisation in a centrifugal pump,” Renewable Energy, vol. 151, pp. 475–487, 2020.

D. Wu, Y. Ren, J. Mou et al., “Unsteady flow and structural behaviors of centrifugal pump under cavitation conditions[J],” Chinese Journal of Mechanical Engineering, vol. 32, no. 01, pp. 127–141, 2019.

X. Li, Z. Jiang, and Z. Zhu, “Entropy generation analysis for the cavitating head-drop characteristic of a centrifugal pump,” Archive Proceedings of the Institution of Mechanical Engineers Part C Journal of Mechanical Engineering Science, vol. 203-210, Article ID 095440621775345, 2018, https://www.semanticscholar.org/author/Q.-Si/2164738https://www.semanticscholar.org/author/Y.-Li/47002268.

P. A. Kumar, R. J. Shruti, and T. Rajiv, “Prediction of flow blockages and impending cavitation in centrifugal pumps using Support Vector Machine (SVM) algorithms based on vibration measurements,” Measurement, vol. 130, pp. 44–56, 2018.

X. He, W. Jiao, C. Wang, and C. Weidong, “Influence of surface roughness on the pump performance based on computational fluid dynamics,” IEEE Access, vol. 99, p. 1, 2019.

H. S. Shim, K. Y. Kim, and Y. S. Choi, “Three-objective optimization of a centrifugal pump to reduce flow recirculation and cavitation,” Journal of Fluids Engineering, vol. 140, no. 9, Article ID 091202, 14 pages, 2018.

Y. Li, G. Feng, and X. LiZ. Zuchao, “An experimental study on the cavitation vibration characteristics of a centrifugal pump at normal flow rate,” Journal of Mechanical Science and Technology, vol. 32, no. 10, pp. 4711–4720, 2018, https://www.researchgate.net/profile/Quiorui_Si.

A. R. Al-Obaidi, “Investigation of effect of pump rotational speed on performance and detection of cavitation within a centrifugal pump using vibration analysis,” Heliyon, vol. 5, no. 6, 2019.

D. O. W. Dimar, V. D. S. Martijn, and N. P. Kruyt, “Prediction of the effect of impeller trimming on the hydraulic performance of low specific-speed centrifugal pumps,” Journal of Fluids Engineering, vol. 140, no. 8, 2018.

Y. Yuan, S. Yuan, and L. Tang, “Investigation on the effect of complex impeller on vibration characteristics for a high-speed centrifugal pump,” Proceedings of the Institution of Mechanical Engineers Part A Journal of Power and Energy, vol. 234, no. 5, Article ID 095765091987831, 2019.

B. Cui, C. Zhang, Y. Zhang, and Z. Zhu, “Influence of cutting angle of blade trailing edge on unsteady flow in a centrifugal pump under off-design conditions,” Applied Sciences, vol. 10, no. 2, p. 580, 2020.

Y. Song, Z. Yu, G. Shi et al., “Influence of impeller staggered arrangement on radial force and pressure fluctuation for a double-suction centrifugal pump,” Advances in Mechanical Engineering, vol. 10, p. 6, 2018.

J. M. Jin, K. J. Jung, Y. H. Kim, and Y.-J. Kim, “Effects of gap between impeller and volute tongue on a pressure fluctuation in double-suction pump,” Journal of Mechanical Ence and Technology, pp. 1–7, 2020.

D. C. Fu, F. J. Wang, P. J. Zhou et al., “Impact of impeller stagger angles on pressure fluctuation for a double-suction centrifugal pump,” Chinese Journal of Mechanical Engineering, vol. 31, no. 1, p. 10, 2018.

H. Yousefi, Y. Noorollahi, M. Tahani, R. Fahimi, and S. Saremian, “Numerical simulation for obtaining optimal impeller’s blade parameters of a centrifugal pump for high-viscosity fluid pumping,” Sustainable Energy Technologies and Assessments, vol. 34, pp. 16–26, 2019.

N. A. Mandhare, K. Karunamurthy, and S. Ismail, “Computational investigation of flow over nose cap of closed impeller of mixed flow centrifugal pump,” Journal of Pressure Vessel Technology, vol. 142, no. 3, 2019.

L. Dong, Y. Zhao, H. Liu et al., “The effect of front streamline wrapping angle variation in a super-low specific speed centrifugal pump,” Journal of Mechanical Engineering Science, vol. 232, no. 23, Article ID 095440621877260, 2018.

D. Zindani, A. K. Roy, and K. Kumar, “Design of mixed flow pump impeller blade using mean stream line theory and its analysis,” Scientia Iranica, 2018.

Y. Tao, S. Yuan, and J. LiuT. Jianping, “The influence of the blade thickness on the pressure pulsations in a ceramic centrifugal slurry pump with annular volute,” Journal of Power and Energy, vol. 231, no. a5, pp. 415–431, 2017, https://www.researchgate.net/profile/Fan_Zhang410.

R. Tao, R. Xiao, and Z. Wang, “Influence of blade leading-edge shape on cavitation in a centrifugal pump impeller,” Energies, vol. 11, no. 10, p. 11, 2018.

G. Peng, Q. Chen, L. Zhou, B. Pan, and Y. Zhu, “Effect of blade outlet angle on the flow field and preventing overload in a centrifugal pump,” Micromachines, vol. 11, no. 9, p. 811, 2020.

W. Yu-Qin and D. Ze-Wen, “Influence of blade number on flow-induced noise of centrifugal pump based on CFD/CA,” Vacuum, vol. 172, Article ID 109058, 2019.

X. Han, Y. Kang, J. Sheng et al., “Centrifugal pump impeller and volute shape optimization via combined NUMECA, genetic algorithm, and back propagation neural network,” Structural and Multidiplinary Optimization, vol. 61, no. 3, 2019.

R. H. Zhang, L. C. Yun, and J. Li, “Investigation on the effect of impeller slot jet on centrifugal pump performance,” Journal of Hydrodynamics, vol. 31, no. 3, 2018.

H. Wang, B. Long, C. Wang et al., “Effects of the impeller blade with a slot structure on the centrifugal pump performance,” Energies, vol. 13, no. 7, 2020, https://www.researchgate.net/deref/http%3A%2F%2Fdx.doi.org%2F10.3390%2Fen13071628.

J. Li, L. Tang, and Y. Zhang, “The influence of blade angle on the performance of plastic centrifugal pump,” Advances in Materials Science and Engineering, vol. 2020, Article ID 7205717, 20 pages, 2020.

J. Hu, Y. Chen, Y. Wang, J. Cheng, and J. Zhao, “New method for drawing two arc cylindrical blade of centrifugal pump,” Chemical Equipment Technology, vol. 37, no. 2, pp. 24–26, 2016.

D. M. S. El-Gazzar and S. A. A. F. Hawash, “Optimizing centrifugal pump performance by different blade configuration patterns,” American Journal of Mechanical and Industrial Engineering, vol. 3, no. 1, p. 1, 2018.

J. Ming, L. Dong, H. Liu et al., “Research of variable angle spiral method for cylindrical blade profile in centrifugal pump,” Water Resources and Power, vol. 35, no. 1, pp. 180–183, 2017.