Research Article | Open Access
A One-Dimensional Flow Analysis for the Prediction of Centrifugal Pump Performance Characteristics
A one-dimensional flow procedure for analytical study of centrifugal pump performance is done applying the principle theories of turbomachines. Euler equation and energy equation are manipulated to find pump performance parameters at different discharge coefficients. Fluid slippage loss at impeller exit and volute loss are estimated. The fluid slippage is modeled by the slip factor approach using Wiesner empirical expression. The volute loss model counts friction loss associated with the volute throw flow velocity, diffusion friction loss due to circulation associated with volute flow, loss due to vanishing of radial flow at volute outlet, and loss inside pump volute throat. Models for impeller hydraulic friction power loss, disk friction power loss, internal flow leakage power loss, and inlet shock circulation power loss are considered by suitable models. Pump internal volumetric flow leakage and volumetric efficiency are related to pump geometry and flow properties. The procedure adopted in this paper is capable of obtaining performance characteristic curves of centrifugal pump in a dimensionless form. Pump head coefficient, manometric efficiency, power coefficient, and required NPSH are characterized. The predicted coefficients and obtained performance curves are consistent with experimental characteristics of centrifugal pump.
Centrifugal pumps are used in various applications and are integral to many industries. Yet, in spite of their prevalence and relatively simple configurations compared to other turbomachines, designing an efficient and durable pump remains a challenge.
The design of centrifugal pumps is still determined empirically because it relies on the use of a number of experimental and statistical rules. However, during the last few years, the design and performance analysis of turbomachinery have experienced great progress due to the joint evolution of computer power and the accuracy of numerical methods.
The one-dimensional performance analysis has proved to be an effective and important approach on pump design . Analytical calculations of pump characteristics depend on geometrical dimensions of pump and losses models in different parts of pump. A series of formulae for calculating losses exist [2–5], but they lack accuracy when applied to centrifugal pumps.
In this work, suggested models for calculating several losses in pump are introduced to examine its validity in evaluating pump performance. This paper is an effort towards theoretically obtaining accurate centrifugal pump performance characteristics. Pump characteristics and parameters are presented in dimensionless forms. It presents a one-dimensional flow analysis procedure towards obtaining optimum centrifugal pump design parameters.
2. Theoretical Analysis
The pump flow coefficient and pump speed coefficient are defined as where is the pump manometric head, is the flow radial velocity at impeller outlet, and is the tangential velocity at outlet of impeller.
2.1. Pump Leakage
Due to the difference between the outlet pressure and the inlet pressure of the impeller, a portion of impeller outlet flow rate, , returns to the impeller inlet from the existing clearances between the impeller and the casing, Figure 1. This internal discharge leakage, , causes some losses as the flow rate through the impeller is greater than the pump useful outlet discharge .
The volumetric efficiency of pump is defined as the ratio of pump outlet discharge to the impeller discharge: in which where is the blade width at inlet, is the impeller diameter at inlet, is the blade thickness coefficient at impeller inlet, , is the blade width at outlet, is the impeller diameter at outlet, is the blade thickness coefficient at impeller outlet, , is the blade thickness, is the number of blades, is the blade angle at impeller outlet, and is the blade angle at impeller inlet. The value of is about 0.95 . The relationship between and with constant blade thickness is given as
The impeller inlet flow velocity coefficient is calculated from (3) dividing both sides by and assuming that the flow enters the impeller without swirl ():
2.2. Outlet Velocity Diagram
A fluid slippage occurs at impeller exit due to the relative rotation of fluid in a direction opposite to that of impeller. A slip factor defined as could be estimated later in Section 2.3.
From the velocity diagram at the impeller outlet, Figure 2, the tangential component of the outlet flow absolute velocity is given as
The ratio of outlet swirl velocity to outlet tangential velocity is and thus,
The outlet tangential velocity (from (7)) and the pump speed coefficient are given as
2.3. Slip Factor
2.3.1. The Relative Eddy (Eddy Circulation)
A simple explanation for the slip effect in an impeller is obtained from the idea of a relative eddy. Suppose that an irrotational and frictionless fluid flow is possible which passes through an impeller. If the absolute flow enters the impeller without spin, then at outlet the spin of the absolute flow must still be zero. The impeller itself has an angular velocity so that, relative to the impeller, the fluid has an angular velocity of ; this is termed the relative eddy. At outlet of impeller, the relative flow, , can be regarded as a through flow on which a relative eddy is superimposed. The net effect of these two motions is that the average relative flow emerging from the impeller passages is at an angle to the vanes and in a direction opposite to the blade motion.
One of the earliest and simplest expressions for the slip factor was obtained by Stodola, cited in . Referring to Figure 3, the slip velocity due to relative eddy, , is considered to be the product of the relative eddy and the radius of a circle which can be inscribed within the channel.
Thus, An approximate expression for can be written if the number of blades is not small: Since the relative eddy angular velocity , then The slip factor is given by Substituting for from (16) and for from (8), the formulae proposed by Stodola, cited in  for the calculation of slip factor, , are obtained which are At () Therefore, where
Wiesner  introduced an empirical expression which extremely well fits the experimental results of slip factor for wide range of practical blade angles and number of blades. It is used in this work and given as where is the limiting diameter ratio,
It is assumed that the water entering the pump impeller is purely in the radial direction. Relative to the impeller, the fluid has an angular velocity of , relative eddy, and thus the relative flow at blade inlet acquires an additional component opposite to rotational direction, as seen in Figure 4, which is where Hence,
When the pump that operates at a discharge differs from that at designed condition , the relative flow velocity at blade inlet tends to acquire an additional component in counter of the rotational direction, . So, the flow enters the blade passage tangent to the blade surface, and a shock eddy or a shock circulation exists prior to the blade leading edge inside pump eye. As could be noticed from Figure 5, the relative velocity additional rotational speed at blade inlet equals
2.4. Euler Equation of Turbomachines
In the case that the fluid entering the pump impeller is purely in the radial direction without swirl, the pump Euler head is given as [6, 7] Due to the fluid slippage at impeller exit, the actual head given to fluid by the impeller, , is calculated from [6, 7]: And so, the slippage head loss is
2.5. Pump Volute
The flow that discharges from the impeller requires careful handling in order to preserve the gains in energy imparted to the fluid. This requires the conversion of velocity head to pressure head by means of a diffuser, and this inevitably implies hydraulic losses. The application of mass conservation to a volute element, , reveals that the discharge flow from impeller is matched to the flow in the volute if . This requires a circumferentially uniform rate of increase of the volute area over the entire development of the spiral (). Consequently, for a given impeller, there exist a specific volute angle and a specific volute throat inlet area for the volute geometry.
The volute angle , Figure 6, is chosen to match the angle of flow entering the volute at a certain pump operating point. With a uniform rate of increase of the volute area, the volute angle is The volute outlet flow velocity and the volute outlet flow velocity coefficient are The throat is assumed to have an expanding angle , and hence the throat diameter equals With the assumption that the throat height , thus The throat outlet flow velocity , and using (2) and (3), then The throat outlet flow velocity coefficient equals It is assumed that the throat diameter equals the eye diameter ; that is, , which equals the suction pipe diameter . Thus, the throat outlet flow velocity, , equals the flow velocity at suction pipe, . Therefore, And the eye velocity coefficient is The eye diameter relative to the impeller inlet diameter is therefore The ratio must not exceed 1, and thus an upper limit is imposed on the volute angle: The throat has a cone angle , where Substituting from (35),
2.6. Volute Loss Model
Prediction models that account for the main features of the swirling flow in volutes, reviewed in , do not account for the circulatory flow initiated in volute at off-design discharge operation of pump. In this work, a model is proposed to account the volute head loss at off-design pump operation.
The flow enters the volute with a through-velocity at an angle (which may differ from that of the volute angle ) on which an eddy of a tangential velocity is superimposed opposite to impeller motion. This inlet volute velocity is decomposed into a velocity parallel to the direction of volute, , and another one in the tangential direction of impeller, (Figure 6). The velocity component is the motive of a second circulatory motion given to the volute flow in direction of impeller motion. Thus, the net circulation velocity of flow in volute is in direction of impeller motion. The component of the velocity in the tangential direction denoted as equals the volute outlet velocity .
The volute loss model counts friction loss associated with the volute throw flow velocity, diffusion friction loss due to circulation associated with volute flow, loss due to vanishing of radial flow at volute outlet, and loss inside pump volute throat. Consequently, the volute head loss and the volute head loss relative to the pump head are written, respectively, as where is the volute friction loss coefficient: where is the volute friction coefficient, is the volute hydraulic diameter, and is the average volute length: The average volute hydraulic diameter relative to impeller diameter is The volute friction coefficient which corresponds to a pipe flow is function of the volute Reynolds number and the roughness , : where The volute Reynolds number is calculated as where In the first term of (44) and (45), is the volute throw flow velocity (Figure 6) and given as In the second term of (44) and (45), is the volute diffusion loss coefficient which could be assumed to have the value , and is the volute circulatory velocity component: The third term of (45) is In the fourth term of (44) and (45), the volute throat head loss could be calculated as where is the volute throat friction loss coefficient assumed to be, , where is the throat cone angle.
2.7. Pump Eye Head Loss
2.8. Pump Manometric Head
The pump manometric head is the difference in static pressure heads between pump outlet and pump eye: With the assumption that the throat diameter equals the eye diameter , the flow velocity is the same at pump suction pipe and pump delivery pipe and neglecting the difference in elevation head across the pump, then or Define a parameter as Using (45) as well as (57) and (61) then The sum of volute and eye head losses, from (59), is written as Using (9), (29), and (59), the pump manometric head is
Coefficients. There are two additional groups of coefficients, namely, the pump head coefficients and head loss coefficients.
The pump head coefficients are the pump manometric head coefficient, , Euler head coefficient, , and the head coefficient at impeller outlet, . They are defined and given, respectively, as According to (21) and (65b), .
The head loss coefficients include the slippage head loss coefficient, the volute-eye head loss coefficient, , the eye head loss coefficient, , and the volute head loss coefficient, . They are given, respectively, as The relationship between the pump speed coefficient and the pump flow coefficient is derived as follows.
Using (11), (64) becomes Dividing both sides by , noting that , and using (11), (64) becomes a quadric equation for : which has a solution (since should be greater than , whence only the +ve sign is considered) Multiplying (68) by and then subtracting from both sides and using (11) yield the following relation for :
2.9. Pump Shaft Power , and Pump Shaft Head
The total shaft power required to drive the impeller is where is the impeller power given to water, is the power lost in friction inside impeller, is the power needed for given circulation to flow at impeller inlet, and is the power lost in friction on outside surface of impeller disks.
The total shaft power and pump shaft head can be simplified, respectively, to in which where is the impeller skin friction head loss, is the inlet shock circulation head loss, is the volumetric (leakage) head loss, and is the disk friction head loss.
The shaft head and the shaft head coefficient are given, respectively, as where is the impeller skin friction head loss coefficient, is the inlet shock circulation head loss coefficient, is the volumetric head loss coefficient, and is the disk friction head loss coefficient.
2.10. Pump Efficiency (Manometric Efficiency)
The pump manometric efficiency is the ratio of gained water power () to the pump shaft power () supplied to pump impeller. According to its definition, it takes the following forms:
2.11. Pump Shaft Power Coefficient
The pump shaft power coefficient and the pump water power coefficient are given, respectively, as
2.12. Impeller Skin Friction Power Head Loss
The impeller hydraulic friction head loss is estimated by the theory of flow through pipes and is given by : where is the dissipation coefficient, is the blade length, is the hydraulic diameter, and is the average relative velocity. Therefore, the impeller skin friction head loss coefficient could be calculated from The hydraulic diameter and the average relative velocity are given, respectively, as, Gülich  Substituting for , (26), , (15), and , (3), yield: The blade length is , and thus where The impeller dissipation coefficient is given as, Gülich  The impeller friction coefficient is function of the average impeller Reynolds number Re and the roughness : where where is defined by (53).
2.13. Disk Friction Head Loss
The disk friction power loss is the power loss in the fluid between external surfaces of the impeller disks and internal walls of the pump casing, Figure 7. The can be estimated as with where is the shear stress in circumferential direction and is the disk skin friction coefficient, and it is assumed constant along the disk surface.
Thus, Therefore, the final expression after the integration for the disk friction power loss and, consequently, the impeller disk friction head loss are
Correlations for were obtained by Kruyt, cited in . Four different regimes were identified, Figure 8: Regime I (laminar flow, boundary layers have merged), Regime II (laminar flow with two separate boundary layers), Regime III (turbulent flow, boundary layers have merged), and Regime IV (turbulent flow with two separate boundary layers). These regimes are characterized by the Reynolds number, , and a nondimensional gap parameter, , where is the axial gap between impeller disk and casing (Figure 7). The equations of curves 1 to 5 in Figure 8 are , , , , and , respectively. The disk friction coefficient for each regime is given as
2.14. Volumetric Head Loss
The volumetric (leakage) head loss and the volumetric head loss coefficient are given as follows after using (2): The leakage flow rate can be estimated using orifice formula : where is the clearance area of wearing ring , is the clearance of wearing ring, Figure 7, is the leakage discharge coefficient (≈0.6), and is the pressure head drop across the clearance: where is the pressure before clearance and is the pressure after clearance at suction side of the impeller (Figure 7).
The pressure head difference between volute and pump suction side equals The pressure distribution along the radial direction of impeller shroud is parabolic , and thus the pressure head difference between impeller outlet and before clearance is given by
2.15. Inlet Shock Circulation Head Loss
At impeller inlet, a shock circulation exists inside pump eye when the flow discharge differs from that of designed one. As discussed before in Section 2.3, the tangential velocity responsible for this shock eddy is , which could be estimated from (28).
Therefore, the inlet shock circulation head loss equals Thus, the inlet shock circulation head loss coefficient is
2.16. Required Net Positive Suction Head
The pump net positive suction head is defined as the difference between the fluid inlet stagnation pressure head and vapour pressure head : where , are the absolute pressure and absolute velocity of fluid at impeller inlet and is the absolute vapour pressure of fluid at the corresponding fluid temperature.
In the vicinity of the leading edge of the impeller blades, the fluid has to accelerate in order to follow the rotating movement of the blades. This acceleration leads to a drop of the static pressure, which results in a local minimum pressure at blade inlet: Therefore, the pump required net positive suction head is The pump coefficient is defined as Referring to Figure 5, , and thus The term is a design parameter for the pump and could be assumed to be 0.02.
From (3) , and therefore the coefficient is
2.17. Pump Specific Speed
The pump specific speed is defined with suitable units as Substitute for Therefore,
3. Calculation Procedure and Results
An iterative procedure using Microsoft Office EXCEL program is performed to calculate the pump design parameters described by the preceding derived equations. The input parameters are the blades number, , impeller inlet diameter to outlet diameter ratio, , volute diameter to impeller outlet diameter ratio, , ratio of blade width at impeller outlet to impeller outlet diameter, , blade width at impeller inlet to blade width at impeller outlet ratio, , volute width to blade width at impeller outlet ratio, , blade angle at impeller outlet, , blade angle at impeller inlet, , volute angle, , throat angle, , ratio of wearing ring clearance/impeller eye diameter, , leakage discharge coefficient, , volute diffusion loss coefficient, , eye loss coefficient, , blade thickness coefficient at impeller outlet, , with constant-thickness blades, ratio of axial gap between impeller disk and casing to impeller outlet radius, , Reynolds number at impeller outlet, , and relative roughness, approximated values are used by putting , , , : , and .
For a certain flow coefficient, , all pump performance parameters and coefficients are calculated after iterations.
Case Study. For reason of comparison between results of present analytical study and experimental pump performance obtained by Baun and Flack , calculations of centrifugal pump performance are performed with the following input parameters given in Table 1.
Other pump constant-parameters are calculated according to the present procedure as shown in Table 2.