ISRN Mechanical EngineeringVolume 2011 (2011), Article ID 794341, 16 pagesdoi:10.5402/2011/794341

Research Article

Numerical Identification of Key Design Parameters Enhancing the Centrifugal Pump Performance: Impeller, Impeller-Volute, and Impeller-Diffuser

Massinissa Djerroud, Guyh Dituba Ngoma, and Walid Ghie

Department of Applied Sciences, University of Quebec in Abitibi-Témiscamingue, 445 Boulevard de l'Université, Rouyn-Noranda, QC, Canada, J9X 5E4

Received 8 July 2011; Accepted 11 August 2011

Academic Editors: S. W. Chang, A. M. Sadegh, and S. C. M. Yu

Copyright © 2011 Massinissa Djerroud 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.

This paper presents the numerical investigation of the effects that the pertinent design parameters, including the blade height, the blade number, the outlet blade angle, the blade width, and the impeller diameter, have on the steady state liquid flow in a three-dimensional centrifugal pump. Three cases were considered for this study: impeller, combined impeller and volute, and combined impeller and diffuser. The continuity and Navier-Stokes equations with the k-ε turbulence model and the standard wall functions were used by means of ANSYS-CFX code. The results achieved reveal that the selected key design parameters have an impact on the centrifugal pump performance describing the pump head, the brake horsepower, and the overall efficiency. To valid the developed approach, the results of numerical simulation were compared with the experimental results considering the case of combined impeller and diffuser.

1. Introduction

At the present time, single and multistage centrifugal pumps are widely used in industrial and mining enterprises. One of the most important components of a centrifugal pump [1] is the impeller. The performance characteristics related to the pump comprising the head, the brake horsepower, and the overall efficiency rely a great deal on the impeller. To achieve better performance for a centrifugal pump, design parameters such as the number of blades for the impeller and the diffuser, the impeller blade angle, the blade height for the impeller and the diffuser, the impeller blade width, the impeller diameter, and the volute radius must be accurately determined, due to the complex liquid flow through a centrifugal pump. This liquid flow is three-dimensional and turbulent. It is therefore important to be aware of the liquid flow’s behavior when traveling through an impeller. This can be done by accounting for the volute and/or the diffuser in the planning, design, and optimization phases at conditions of design and off-design. Many experimental and numerical studies have been carried out on the liquid flow through a centrifugal pump [2–21], where the effects of the number of impeller blades on the pump’s performance were examined experimentally in [11, 12]. The effects of the impeller outlet blade angle on the pump’s performance were also investigated numerically [13, 14], using a CFD code and experimentally in [15]. In [16] the dynamic effects due to the impeller-volute interaction within a centrifugal pump were numerically investigated, whereas the effects of the volute on velocity and pressure fields were examined in [17, 18]. Additional experimental investigation carried out [19] consisted of measuring unsteady velocity, the pressure and flow angle at the centrifugal pump’s impeller outlet, with and without volute casing. The liquid flow and head distribution within a centrifugal pump’s volute were compared with the impeller’s characteristics, without the volute casing. Moreover, two centrifugal pump impellers with different outlet diameters for the same volute were examined both experimentally and numerically [20], to evaluate the influence the radial gap between the impeller exit and the volute tongue had on the unsteady radial forces acting upon the impeller of a centrifugal pump with volute casing. Additionally, the effects of flow behavior in a centrifugal pump, whose diffuser was subjected to different radial gaps, were investigated numerically in [21] using a CFD code. The analysis of previous works clearly demonstrated that research results obtained are specific to the centrifugal pump design parameter values and thus cannot be generalized. In this work therefore a numerical study was performed using a finite volume method according to the CFX code [22] to gain further insight into the characteristics of the three-dimensional turbulent liquid flow through an impeller, a combined impeller and volute, and a combined impeller and diffuser, while also considering various flow conditions and pump design parameters: blade heights of 12 mm, 18 mm, and 24 mm; blade numbers of 5, 7, and 9 for the impeller and 5, 8, and 12 for the diffuser; outlet blade angles of 9°, 28°, and 60°, blade widths of 4 mm, 10 mm, and 15 mm; impeller outer diameters of 285 mm and 320 mm. The reference dimensions selected for the impeller and diffuser were based on the existing impeller and diffuser [23]. Upon applying the continuity and Navier-Stokes equations, the liquid flow velocity and the liquid pressure distributions in an impeller, a combined impeller and volute, and a combined impeller and diffuser were determined, while accounting for boundary conditions. Since the rotating speed of the centrifugal pump under consideration was constant a valve installed on the pump’s discharge side was used to regulate the volume flow rate. We accounted for suction pressure variation as a function of the valve volume flow rate in the numerical simulations being run. The pump head, brake horsepower, and efficiency were represented as a function of the volume flow rate, where the objective was to identify the values of selected key design parameter that might improve pump performance with respect to their value ranges.

2. Mathematical Formulation

Figure 1 shows a centrifugal pump consisting of three components, including an impeller, a diffuser, and a volute [24]. The models selected for the liquid flow in an impeller, a combined impeller and volute, and a combined impeller and diffuser are depicted in Figure 2, placing greater emphasis on the fluid domain.

Figure 1: Centrifugal pump.

Figure 2: Models of centrifugal pump components.

In the governing equations for liquid flow in the centrifugal pump components, the following assumptions were made: (i) a steady state, three-dimensional, and turbulence flow using the k-ε model; (ii) it was an incompressible liquid; (iii) it was a Newtonian liquid; (iv) the liquid’s thermophysical properties were constant with temperature.

To account for these assumptions, the theoretical analysis of the liquid flow in an impeller, a combined impeller and volute, and a combined impeller and diffuser was based on the continuity and Navier-Stokes [22] equations. For the three-dimensional liquid flow through the components of a centrifugal pump as shown in Figure 2, the continuity equations are expressed by and the Navier-Stokes equations are given by where is the liquid flow velocity vector, is the pressure, ρ is the density, is the effective viscosity accounting for turbulence, ⊗ is a tensor product, and is the source term. More particularly, for flows in an impeller rotating at a constant speed ω, the source term can be written as follows: where is the location vector.

In addition, is defined as where is the dynamic viscosity and is the turbulence viscosity.

Since the k-ε turbulence model is used in this work because convergence is better than with other turbulence models, is linked to turbulence kinetic energy k, equation (6), and dissipation , equation (7), via the relationship where is a constant.

The values for k and ε come directly from the differential transport equations for turbulence kinetic energy and turbulence dissipation rates: where , , andare constants. is the turbulence production due to viscous and buoyancy forces, which is modeled using: where can be neglected for the k-ε turbulence model.

Moreover, for the modeling of flow near the wall, the logarithmic wall function is used to model the viscous sublayer [22].

2.1. Impeller

Three velocity types are involved when considering the flow through a centrifugal pump impeller: the tangential velocity ω, the relative velocity , and the absolute velocity . The last is expressed in vector format as follows: Figure 3 shows the velocity triangles at the impeller inlet and outlet at the design conditions where the liquid enters and leaves the impeller at the blade angles and , respectively. The components of and in the direction of are (swirl velocity) and , respectively, while those normal to are and .

Figure 3: Velocity triangles.

Moreover, according to the Euler equation [1], the energy transfer per unit mass of flow for a centrifugal pump can be formulated as where is the ideal pump total head.

Neglecting the swirl velocity at the impeller inlet, (10) can be expressed as follows: When accounting for the hydraulic efficiency, , the actual pump head rise is given by Also, the hydraulic efficiency can be calculated using the following empirical formula [1]: where is the volume flow rate in m³/s. It is given by with A as the flow passage area normal to the meridional direction.

Since in reality the flow through a centrifugal pump is turbulent and three-dimensional, the actual relative flow direction at the impeller exit is different from that of the blade angle. As depicted in Figure 4, the flow angle is always less than the blade angle . This can lead to secondary flows in the flow passage, from the pump inlet to discharge [1].

Figure 4: Flow angle and blade angle.

As such, the slip factor is used to take into account the difference between and , which is formulated as where is the actual swirl flow velocity at the impeller exit and is the ideal swirl flow velocity at the impeller exit.

In addition, the slip velocity is given by Taking into account the slip factor, (12) can be expressed as Moreover, to account for the leakage flow from the impeller, the volumetric efficiency is defined by where is the leakage flow from the impeller exit back to the inlet through the clearance.

In addition, the pump’s mechanical efficiency is formulated as follows: where is the brake horsepower and the power delivered by the impeller to the fluid.

is globally expressed by where is the pump shaft torque, is the centrifugal pump horsepower. It is expressed as is the loss power due to the friction, which is given by is the loss power due to leakage, which is defined as: is the friction loss power in bearings and seals, and is the disk friction power due to impeller shrouds.

in (18) can be formulated as follows: Furthermore, (23) can be rewritten as Accounting for (23), (18) can be expressed as Thus, the overall efficiency of a centrifugal pump can be formulated as Finally, the overall efficiency can also be formulated in terms of the other efficiencies as

2.2. Volute Parameters

Figure 5 shows the parameters of a volute without diffuser defined by the radius of volute basic circle , the radius of volute cut water circle , the volute angle , the volute cross-sectional area , which depends on the angle , and the volute outlet cross-sectional area [1].

Figure 5: Impeller-volute without diffuser.

The average flow velocity at the volute outlet is given by where the volute velocity constant is an empirical parameter correlated with the specific speed, as shown in Figure 6 along with other volute parameters such as the volute angle and the volute basic circle diameter [1].

Figure 6: Volute velocity constant, volute angle, and diameter of volute basic circle versus specific speed.

In addition, the volute cross-sectional area can be formulated as where is the centroid radius of the volute cross-sectional area, is the angular momentum of flow at the impeller outlet, it is expressed as , and to account for friction loss.

2.3. Diffuser Velocity and Pressure Difference

Figure 7 shows the velocities at the inlet and the outlet of a vaned diffuser immediately downstream from the impeller.

Figure 7: Velocities at the diffuser inlet and outlet.

The outlet velocity can be determined using [24] where where is the flow passage area normal to the meridional direction at the diffuser outlet.

Finally, the pressure difference between the diffuser outlet and inlet is given by To solve (1) and (2) numerically while accounting for the boundary conditions and the turbulence model k-ε, the computational fluid dynamics ANSYS-CFX code, based on the finite volume method, was used to obtain the liquid flow velocity and the pressure distributions. In the cases examined involving the impeller, combined impeller and volute, and combined impeller and diffuser, the boundary conditions were formulated as follows: the static pressure provided was given at the inlet, while the flow rate provided was specified at the outlet. The frozen rotor condition was used for both the impeller-volute and the impeller-diffuser interfaces. A no-slip condition was set for the flow at the wall boundaries. Figure 8 shows the inlet, outlet, and interface domains for the selected centrifugal pump components.

Figure 8: Domains of inlet, outlet, and interface.

Accounting for the fact that the pump rotating speed was constant, the volume flow rate was controlled by a regulator valve, which had an influence on the pressure at the pump inlet as indicated in Figure 9 [23]. This was accounted for in the numerical simulations performed.

Figure 9: Pressure at the pump inlet versus valve volume flow rate regulation.

Furthermore, the ANSYS-CFX code comprises by geometry (DesignModeler), CFX-pre, CFX-solver, and CFX-post modules. According to the applied ANSYS-CFX code, Figure 10 depicts the steps specifically used to obtain the numerical simulation results from the geometry models to the numerical models for the impeller, the combined impeller and volute, and the combined impeller and diffuser.

Figure 10: Steps from 3D geometry model to 3D numerical model and to numerical simulation results.

3. Results and Discussion

Water was used as the working liquid for all simulations run and for use in this study considered to have the following reference values: temperature of 25°C for water, density of ρ = 997 kg/m³, and dynamic viscosity of μ = 8.899 × 10⁻⁴ Pa s. The main data for the reference impeller, volute, and diffuser are given in Tables 1, 2, and 3.

Table 1: Main data of the reference impeller [23].

Inlet diameter (mm)
145

Outlet diameter (mm)
320

Inlet blade angle (degree)
11.69

Outlet blade angle (degree)
28

Inlet blade width (mm)
12

Blade thickness (mm)
4

Number of blades
7

Rotating speed (rpm)
1800

Table 2: Main data of the reference volute [1].

Volute angle (degree)
Volute radius (mm)
Volute angle (degree)
Volute radius (mm)

0
165
225
278.96

45
183.79
270
302.76

90
207.58
315
326.55

135
231.38
360
350.35

180
255.17

Table 3: Main data of the diffuser [23].

Inlet diameter (mm)
321.536

Outlet diameter (mm)
441.77

Blade width (mm)
12

Blade thickness (mm)
3.401

Inlet blade angle (degree)
11.07

Outlet blade angle (degree)
39.42

Number of blades
8

3.1. Case Studies

Six key design parameters of a centrifugal pump were selected for an examination of their effects mainly on the pump performance: impeller blade height without volute, impeller blade width without volute, impeller blade angle without volute, impeller blade number with volute, impeller diameter with volute, and diffuser blade number with impeller. For the highest accuracy of numerical simulation results, mesh-independent solution tests were conducted in each case study by varying the number of mesh elements. Table 4 indicates the required number of mesh elements to achieve mesh-independent results.

Table 4: Number of mesh elements.

Case study
Mesh element number

Impeller outlet blade height ( in mm)

    12
43380

    18
48776

    24
74210

Impeller blade width (e in mm)

    4
43380

    10
75841

    15
119945

Impeller outlet blade angle ( in degree)

    9
38508

    28
43380

    60
42288

Impeller blade number when accounting for volute

    5
28287

    7
47979

    9
132852

Impeller diameter when accounting for volute ( in mm)

    285
47377

    320
47979

Diffuser blade number

    5
39126

    8
43380

    12
45117

3.1.1. Effect of Impeller Outlet Blade Height

To analyze the outlet blade height’s effect on the pump head, the pump brake horsepower, and the overall pump efficiency, the values 0.012 m, 0.018 m, and 0.024 m were selected for outlet blade height, while keeping the other parameters constant. Figure 11 shows the pump head as a function of the volume flow rate with the outlet blade height as a parameter. There, it can be clearly observed that the pump head decreases with increasing volume flow rate due to decreasing liquid pressure. In addition, the pump head increases with increasing outlet blade height. This can be explained by the fact that, when the volume flow rate is kept constant, the increased outlet blade height leads to the decreasing meridional velocity, which increases the pump head since the outlet tangential velocity and the outlet blade angle remain constant. In other words, the liquid pressure drop in the impeller decreases as a function of the increase in outlet blade height.

Figure 11: Pump head versus volume flow rate (parameter: blade height).

The curves expressing the pump brake horsepower as a function of the volume flow rate are shown in Figure 12, illustrating that the brake horsepower increases with increasing volume flow rate. This can be explained by the additional decrease in liquid pressure relative to the volume flow rate. Also, the brake horsepower increases relative to the impeller blade height due to the requested increase in pump shaft torque relative to the increased blade height.

Figure 12: Brake horsepower versus volume flow rate (parameter: blade height).

As depicted in Figure 13, the curves representing overall pump efficiency as a function of volume flow rate illustrate that the overall efficiency for = 12 mm decreases rapidly to the right of the best efficiency point (BEP). The overall efficiency for = 18 is highest when the volume flow rate reaches 0.08 m³/s.

Figure 13: Overall efficiency versus volume flow rate (parameter: blade height).

Figures 14, 15, and 16 show the corresponding contours for static pressure, liquid flow velocity vectors, and streamlined liquid flow velocities for Q = 0.065 m³/s. From these figures it can be observed that the static pressure is higher at the impeller outlet than at the impeller inlet. This is due to the decrease in liquid flow velocity at the impeller outlet. As such these figures clearly illustrate the correlation between variations in liquid flow velocity and static pressure. Moreover, Figures 14–16 illustrate the impact of variations in blade height on static pressure, liquid flow velocity, and streamlined liquid velocity, respectively, where average liquid flow velocities at the impeller outlet were 15.92 m/s, 12.64 m/s, and 10.56 m/s for = 12 mm, = 18 mm, and = 24 mm, respectively.

(a)  mm  Pa

(b)  mm Pa

(c)  mm Pa

Figure 14: Static pressure contour.

(a)  mm

(b)  mm

(c)  mm

Figure 15: Liquid flow velocity vector.

(a)  mm

(b)  mm

(c)  mm

Figure 16: Streamlines of liquid flow velocity.

3.1.2. Effect of Impeller Blade Width

To investigate the effect that the impeller blade width has on the pump head, the pump brake horsepower, and the pump overall efficiency, the blade widths of 4 mm, 10 mm, and 15 mm were selected, while the other parameters were keep constant. Figure 17 shows the pump head as a function of the volume flow rate, illustrating that the pump head decreases with increased blade width. This is due to augmenting the liquid pressure drop with increasing blade width. Also, the required pump brake horsepower decreases when the blade width rises, as indicated in Figure 18. The corresponding overall efficiency curves are shown in Figure 19, illustrating that the blade width’s impact on the overall efficiency is more pronounced in at high volume flow rates. In other words, the overall efficiencies for the three blade widths decrease rapidly to the right side of the BEP and the lowest overall efficiency is obtained when e = 15 mm.

Figure 17: Pump head versus volume flow rate (parameter: blade width).

Figure 18: Pump brake horsepower versus volume flow rate (parameter: blade width).

Figure 19: Overall efficiency versus volume flow rate (parameter: blade width).

3.1.3. Effect of Impeller Outlet Blade Angle

Three impeller outlet blade angle values of 9°, 28°, and 60° were selected to investigate their influence on the pump head, the pump brake horsepower, and the pump’s overall efficiency. Figure 20 depicts the distribution of the pump head as a function of volume flow rate and with outlet blade angle as a parameter. This figure thus shows that the pump head increases with increasing outlet blade angle, which can be explained by the increased outlet cross-section size relative to the increased outlet blade angle, thus leading to diminution of liquid pressure drop in flow passage between blades.

Figure 20: Pump head versus volume flow rate (parameter: outlet blade angle).

In addition, Figure 21 depicts the corresponding brake horsepower curves as a function of the volume flow rate, illustrating that the pump brake horsepower increases relative to the augmenting outlet blade angle. This is due to the increase in the requested shaft torque, along with the augmented outlet blade angle.

Figure 21: Brake horsepower versus volume flow rate (parameter: outlet blade angle).

Moreover, the efficiency curves shown in Figure 22 illustrate that the overall efficiency for decreases rapidly to the right of the BEP.

Figure 22: Overall efficiency versus volume flow rate (parameter: outlet blade angle).

Additionally, Figures 23 and 24 show the static pressure contour and the liquid flow velocity vector for Q = 0.065 m³/s. From these figures, it can be observed that the static pressure difference between the impeller outlet and inlet increases with the augmented blade angle, due to a decrease in the liquid flow velocity at the impeller outlet, as indicated in Figure 26. The average liquid flow velocities at the impeller outlet are 21.06 m/s, 15.92 m/s, and 10.09 m/s for ,, and , respectively.

(a)  Pa

(b)  Pa

(c)  Pa

Figure 23: Static pressure contour.

(a)

(b)

(c)

Figure 24: Liquid flow velocity vector.

3.1.4. Effect of Impeller Blade Number When Accounting for Volute

To investigate the effect of the impeller blade number on the pump head, the pump brake horsepower and the overall pump efficiency, three impellers whose blade number were 5, 7, and 9 were selected, while the other parameters were kept constant. Figure 25 shows the pump head as a function of the volume flow rate, illustrating that the pump head increases with a greater blade number. This is explained by the decrease in the liquid pressure drop in the flow passage with an augmented impeller blade number, keeping the same total volume flow rate. Also, as shown in Figure 26, the pump brake horsepower increases relatively with the augmented blade number. This is due to the increase in the request pump shaft torque, as the pump blade number also increases.

Figure 25: Pump head versus volume flow rate (parameter: impeller blade number).

Figure 26: Brake horsepower versus volume flow rate (parameter: impeller blade number).

In addition, Figure 27 shows the overall efficiency curves, showing that the impeller having 5 blades has the lowest overall efficiency.

Figure 27: Overall efficient versus blade number (parameter: blade number).

Moreover, Figures 28 and 29 depict the corresponding static pressure contour and liquid flow velocity vector for Q = 0.065 m³/s, respectively. These figures thus clearly show the increased static pressure difference between the volute outlet and the impeller inlet relative to the increasing blade number. This confirms the reduction in the liquid flow velocity at the impeller outlet relative to the greater blade number, as represented in Figure 29, where the average liquid flow velocities at the impeller outlet were 16.06 m/s, 15.40 m/s, and 12.53 m/s for 5 blades, 7 blades, and 9 blades, respectively.

(a) 5 blades  Pa

(b) 7 blades  Pa

(c) 9 blades  Pa

Figure 28: Static pressure contour.

(a) 5 blades

(b) 7 blades

(c) 9 blades

Figure 29: Vectors of liquid flow velocity contour.

3.1.5. Effect of Impeller Diameter When Accounting for Volute

The impeller outlet diameter values of 285 mm and 320 mm were selected to investigate their effects on pump performance when keeping the other parameters constant. Figure 30 shows that the pump head increases with increasing impeller diameter, which can be explained by the fact that the liquid static pressure drop in impeller decreases with increasing impeller diameter. In other words, for a given volume flow rate, the pressure difference between the volute outlet and the impeller inlet is higher for an impeller with a greater diameter. In addition, Figure 31 shows that the brake horsepower increases relative to the increasing impeller diameter, due to the requested augmented impeller shaft torque relative to the size of the impeller diameter.

Figure 30: Pump head versus volume flow rate (parameter: impeller diameter).

Figure 31: Brake horsepower versus volume flow rate (parameter: impeller diameter).

Moreover, the corresponding overall efficiency curves shown in Figure 32 indicate that the impeller having a great diameter has better overall efficiency with volume flow rates greater than 0.02 m³/s.

Figure 32: Overall efficiency versus volume flow rate (parameter: impeller diameter).

3.1.6. Effect of Diffuser Blade Number

To analyze the effect the diffuser blade number has on the pump head, the pump brake horsepower, and the overall pump efficiency, three diffuser models with blade numbers of 5, 8, and 12 were selected, while the other parameters were kept constant. Figure 33 shows the pump head as a function of the volume flow rate, where it is observed that the impact of the diffuser number on the pump head is small, even if the pump head for the diffuser blade number of 8 is highest for the Q between 0.012 m³/s and 0.055 m³/s. As depicted in Figure 34, the variation in brake horsepower due to diffuser blade number is also small, even if the diffuser blade number of 12 corresponds to a lowest brake horsepower.

Figure 33: Pump head versus volume flow rate (parameter: diffuser blade number).

Figure 34: Brake horsepower versus volume flow rate (parameter: diffuser blade number).

Furthermore, Figure 35 shows that, for the low and the high volume flow rates, the overall efficiency for the diffuser blade number 12 is highest whereas the overall efficiencies for diffuser blade numbers of 5 and 8 are nearly the same. This figure also indicates that the overall efficiency is lowest for diffuser blade number 5.

Figure 35: Overall efficiency versus volume flow rate (parameter: diffuser blade number).

Additionally, Figures 36 and 37 depict the corresponding static pressure contour and the liquid flow velocity vector for Q = 0.065, respectively, clearly showing for these figures the correlation between the increased static pressure difference and decreased liquid flow velocity at the diffuser outlet, with augmented diffuser blade number. Thus, the average values obtained for the static pressure difference between the diffuser outlet and the impeller inlet are 3.428 × 10⁵ Pa, 3.49 × 10⁵ Pa, and 3.65 × 10⁵ Pa for blade numbers of 5, 8, and 12, respectively, as represented in Figure 36. Also, the average liquid flow velocity values at the diffuser outlet of 15.13 m/s, 12.22 m/s, and 9.06 m/s were found for the blade numbers of 5, 8, and 12, respectively, as shown in Figure 37.

(a) 5 blades Pa

(b) 8 blades Pa

(c) 12 blades Pa

Figure 36: Static pressure contour.

(a) 5 blades

(b) 8 blades

(c) 12 blades

Figure 37: Vectors of liquid flow velocity.

3.2. Model Comparison

The impeller and diffuser combination was selected to validate the numerical approach developed, since the experimental results of this case were available from Technosub Inc. When accounting for experimental boundary conditions for the numerical simulations run, Figures 38–40 show the comparison between the experimental and the numerical results for the pump head, the brake horsepower, and the overall efficiency. The discrepancies observed could be explained by the fact that lost mechanical power, power lost due to leakage, and the pump casing were not taken into account in the numerical simulations carried out. The horsepower for experimental pump brake was therefore higher than the numerical brake horsepower obtained, as illustrated in Figure 39.

Figure 38: Pump head versus volume flow rate (parameter: numerical and experimental results).

Figure 39: Brake horsepower versus volume flow rate (parameter: numerical and experimental results).

Figure 40: Overall efficiency versus volume flow rate (parameter: numerical and experimental results).

4. Conclusion

In this study, a steady-state liquid flow in a three-dimensional centrifugal pump was numerically investigated. Models of impeller, combined impeller and volute, and combined impeller and diffuser were developed to analyze the effects the key design parameters, including the blade height, the outlet blade angle, the blade width, the blade number, and the impeller outer diameter, had on the pump head, the brake horse power, and the overall efficiency. The obtained results demonstrate, among others, that the pump head and the brake horsepower increase with increasing impeller blade number and impeller blade height, while they decrease with increasing impeller blade width. Also, the interaction between the impeller and the volute reveals that the decrease of the impeller outer diameter keeping the volute dimensions constant leads to the reduction of the pump head and the brake horsepower. The pump overall efficiency is also influenced by the selected key design parameter. A relatively good agreement was observed comparing the developed numerical approach with the experimental results for the case of the combined impeller and diffuser obtained from a pump manufacturer.

Nomenclature

B: | Source term (Nm⁻³)
b: | Height (m)
d: | Diameter (m)
e: | Width (m)
g: | Acceleration of gravity (m s⁻²)
H: | Head (m)
P: | Power (W)
p: | Pressure (Nm⁻²)
: | Turbulence production due to viscous and buoyancy forces
Q: | Volume flow rate (m³ s⁻¹)
r: | Radial coordinate (m)
U: | Velocity or tangential velocity (m s⁻¹)
u: | Flow velocity in direction (m s⁻¹)
V: | Absolute velocity (m s⁻¹)
v: | Flow velocity in direction (m s⁻¹)
W: | Relative velocity (m s⁻¹)
w: | Flow velocity in direction (m s⁻¹)
x: | x-coordinate (m)
y: | y-coordinate (m)
z: | z-coordinate (m).

Greek Symbols

α: | Angle between and (degree)
β: | Blade angle between and (degree)
Δ: | Difference
ε: | Turbulence dissipation (m² s⁻³)
η: | Efficiency
k: | Turbulence kinetic energy (kg m⁻² s⁻²)
θ: | Angle (°)
ρ: | Fluid density (kg m⁻³)
μ: | Dynamic viscosity (Pa s)
: | Effective viscosity (Pa s)
: | Slip factor
: | Turbulence viscosity (Pa s)
ω: | Angular velocity (rad s⁻¹).

Subscripts

1: | Inlet
2: | Outlet or diffuser inlet
3: | Volute outlet or diffuser outlet
b: | Blade
df: | Disk friction
f: | Flow
h: | Hydraulic
i: | Inlet or ideal
imp: | Impeller to fluid
L: | Leakage
m: | Mechanical
o: | Outlet
r: | Radial or perpendicular to the vector
s: | Shaft or slip
u: | Direction of vector
v: | Volumetric or volute
w: | Wall.

Acknowledgments

The authors are grateful to the Foundation of University of Quebec in Abitibi-Temiscamingue (FUQAT) and the company Technosub Inc.

References

W. W. Peng, Fundamentals of Turbomachinery, John Wiley and Sons, Hoboken, NJ, USA, 2008.

W. Zhou, Z. Zhao, T. S. Lee, and S. H. Winoto, “Investigation of flow through centrifugal pump impellers using computational fluid dynamics,” International Journal of Rotating Machinery, vol. 9, no. 1, pp. 49–61, 2003.

J. S. Anagnostopoulos, “A fast numerical method for flow analysis and blade design in centrifugal pump impellers,” Computers & Fluids, vol. 38, no. 2, pp. 284–289, 2009.

B. Cui, Z. Zhu, J. Zhang, and Y. Chen, “Flow simulation and experimental study of low-specific-speed high-speed complex centrifugal impellers,” Chinese Journal of Chemical Engineering, vol. 14, no. 4, pp. 435–441, 2006.

S. Derakhshan and A. Nourbakhsh, “Theoretical, numerical and experimental investigation of centrifugal pumps in reverse operation,” Experimental Thermal and Fluid Science, vol. 32, no. 8, pp. 1620–1627, 2008.

R. Spence and J. Amaral-Teixeira, “Investigation into pressure pulsations in a centrifugal pump using numerical methods supported by industrial tests,” Computers & Fluids, vol. 37, no. 6, pp. 690–704, 2008.

K. W. Cheah, T. S. Lee, S. H. Winoto, and Z. M. Zhao, “Numerical flow simulation in a centrifugal pump at design and off-design conditions,” International Journal of Rotating Machinery, vol. 2007, Article ID 83641, 8 pages, 2007.

M. H. Shojaee Fard, F. A. Boyaghchi, and M. B. Ehghaghi, “Experimental study and three-dimensional numerical flow simulation in a centrifugal pump when handling viscous fluids,” International Journal of Engineering Science, vol. 17, no. 3-4, pp. 53–60, 2006.

R. K. Byskov, C. B. Jacobsen, and N. Pedersen, “Flow in a centrifugal pump impeller at design and off-design conditions—part II: large eddy simulations,” Journal of Fluids Engineering, vol. 125, no. 1, pp. 73–83, 2003.

K. Majidi, “Numerical study of unsteady flow in a centrifugal pump,” Journal of Turbomachinery, vol. 127, no. 2, pp. 363–371, 2005.

W.-G. Li, F.-Z. Su, and C. Xiao, “Influence of the number of impeller blades on the performance of centrifugal oil pumps,” World Pumps, vol. 2002, no. 427, pp. 32–35, 2002.

H. Liu, Y. Wang, S. Yuan, M. Tan, and K. Wang, “Effects of blade number on characteristics of centrifugal pumps,” Chinese Journal of Mechanical Engineering, vol. 23, no. 6, pp. 742–747, 2010.

M. H. Shojaee Fard and F. A. Boyaghchi, “Studies on the influence of various blade outlet angles in a centrifugal pump when handling viscous fluids,” American Journal of Applied Sciences, vol. 4, no. 9, pp. 718–724, 2007.

E. C. Bacharoudis, A. E. Filios, M. D. Mentzos, and D. P. Margaris, “Parametric study of a centrifugal pump impeller by varying the outlet blade angle,” The Open Mechanical Engineering Journal, vol. 2, pp. 75–83, 2008.

M. Gölcü and Y. Pancar, “Investigation of performance characteristics in a pump impeller with low blade discharge angle,” World Pumps, vol. 2005, no. 468, pp. 32–40, 2005.

J. González, J. Fernández-Francos, E. Blanco, and C. Santolaria-Morros, “Numerical simulation of the dynamic effects due to impeller-volute interaction in a centrifugal pump,” Journal of Fluids Engineering, vol. 124, no. 2, pp. 348–355, 2002.

M. Asuaje, F. Bakir, S. Kouidri, F. Kenyery, and R. Rey, “Numerical modelization of the flow in centrifugal pump: volute influence in velocity and pressure fields,” International Journal of Rotating Machinery, vol. 2005, no. 3, pp. 244–255, 2005.

K. A. Kaupert and T. Staubli, “The unsteady pressure field in a high specific speed centrifugal pump impeller—part I: influence of the volute,” Journal of Fluids Engineering, vol. 121, no. 3, pp. 621–626, 1999.

P. Hergt, S. Meschkat, and B. Stoffel, “The flow and head distribution within the volute of a centrifugal pump in comparison with the characteristics of the impeller without casing,” Journal of Computational and Applied Mechanics, vol. 5, no. 2, pp. 275–285, 2004.

J. González Pérez, J. Parrondo, C. Santolaria, and E. Blanco, “Steady and unsteady radial forces for a centrifugal pump with impeller to tongue gap variation,” Journal of Fluids Engineering, vol. 128, no. 3, pp. 454–462, 2006.

A. Ozturk, K. Aydin, B. Sahin, and A. Pinarbasi, “Effect of impeller-diffuser radial gap ratio in a centrifugal pump,” Journal of Scientific and Industrial Research, vol. 68, no. 3, pp. 203–213, 2009.

Ansys inc., ANSYS-CFX, User Manual, USA, 2008

Technosub Inc., http://technosub.net/

G. Lemasson, “Les machines transformatrices d'énergie,” TOME II, Turbomachine-Machine alternatives, 1967