Abstract

Slip factor is an important parameter in the hydraulic design of centrifugal pump impeller for handling viscous oils. How to extract the factor from CFD computational results and how flow rate and liquid viscosity to affect it remain unclear. In the present paper, the flip factor was estimated by means of two approaches: one is from the velocity triangles at the impeller outlet and the other is due to the impeller theoretical head of 3D turbulent viscous fluid. The velocity of water and viscous oils in the impeller and volute computed by CFD was validated with LDV measurements at the best efficiency point. The effect of exit blade angle on slip factor was clarified. It was shown that the two approaches result into two different slip factors. The factors are significantly dependent of flow rate; however, the liquid viscosity seems to take less effect on them. Volute is responsible for reduction in tangential velocity of liquid at the outlet of impeller at low flow rates. The slip factor of impeller with large exit blade angle is not sensitive to flow rate.

1. Introduction

The performance of centrifugal pump handling water and viscous oils was investigated numerically by using a CFD code FLUENT based on a steady, 3D, and incompressible turbulent flow. The turbulence effect was involved with the standard turbulence model and wall roughness was taken into account with the nonequilibrium wall function in [1]. The effect of liquid viscosity on pump performance was clarified through observing the pump head and hydraulic efficiency as well as hydraulic loss coefficient in terms of flow rate. A comparison of computed and experimental overall performance of the pump was made. It was confirmed that the “sudden-rising head effect” exists and is caused from the high viscosity and certain large surface roughness. The volute results in an increasing influence on the flow around the impeller exit at a low flow rate.

Slip factor is one important design parameter for deciding a correct impeller diameter of centrifugal pump. The factor can be obtained theoretically and experimentally. The typical investigations include those conducted by Kasai [2, 3], Sakai and Watanabe [4], Noorbakhsh [5], Whitfield [6], Murata et al. [7], Harada and Senoo [8], Visser et al. [9], von Backstron [10], Hassenpflug [11], Ji et al. [12], Qiu et al. [13], and so on. The slip factor depends on impeller geometry [5, 7, 9, 10, 14], flow rate [24, 6, 13, 15], and viscosity of the liquid pumped [16, 17].

Recently, Slip factor is increasingly estimated by using CFD approach [1822]. Thus it is necessary to explore effects of viscosity and flow rate on slip factor when a centrifugal pump handles viscous oils. The results can be useful and helpful to primary hydraulic design and performance prediction of a viscous oil transporting centrifugal pump by means of a simple and quick method such as the mean line method proposed by Veres [23].

In this paper, it is intended to establish the slip factor dependence upon liquid viscosity and flow rate and show the difference in theoretical head of impeller between CFD and traditional one estimated by the Euler turbomachinery equation and a slip factor correction. Besides, a comparison of LDV measured velocity and CFD computed one in the impeller and volute was presented for water and highly viscous oil at the best efficiency point (BEP) to validate CFD results in the flow field.

2. Computational Models and Slip Factor Extraction

2.1. Computational Models

The computational physical model is the end-suction, single stage centrifugal pump in [1]. The pump flow rate, head, and rotational speed are = 25 m3/h, = 8 m and = 1450 r/min, = (r/min, m3/s, m) = 93 at design point. The impeller eye and outlet diameters are 62 mm and 180 mm, the number of blades is 4, the exit blade angle is 20°, and its warp angle is 140°. The diameter of base circle and width of the volute are 190 mm and 40 mm, and the cross-sectional area of volute throat is 1440 mm2.

The computational fluid domain consists of suction pipe, rotating impeller, and volute, as shown in Figure 1(a). The rotor-stator interfaces are established at the entrance to the impeller and the impeller outlet to couple the flow variables in the impeller and volute. The impeller rotational effect is involved by means of MRF (multiple reference frame) method proposed in FLUENT. For such a flow, the CFD computations were carried out at = 45°, which is the angle between the volute tongue tip and the outer edge of pressure surface of impeller blade.

3D, steady, and incompressible turbulent flow in the pump was solved by FLUENT. The turbulence effect was taken into account with the standard turbulence model, and the wall roughness was involved by the nonequilibrium wall function. The finite volume method, SIMPLE algorithm, and the second-order upwind scheme for the convection terms in the governing equations were adopted. The viscosity and density of liquids are presented in Table 1.

The fluid domain of suction pipe was discretized with hexahedral cells, but those of impeller and volute were meshed with tetrahedral cells. It was found that when the total number of cells was increased up to 910 k (10 k in the suction pipe, 380 k in the impeller, and 520 k in the volute), the pump performance became independent of mesh size.

The relaxation factors for the pressure correction equation, momentum equations, and and transport equations were 0.3, 0.5, 0.8, and 0.8, respectively. The residual tolerances for the pressure correction equation, momentum equations, and transport equations were 1 × 10−4.

There was an inlet velocity boundary condition at the entrance to the suction pipe. The reference pressure is zero. There was an outlet pressure boundary condition at the volute nozzle discharge. The rest boundaries were no-slip solid walls with the roughness Ra = 50 μm.

2.2. Slip Factor Extraction

Slip factor is used to specify the flow slip effect in the exit of a centrifugal pump impeller and is a key parameter to establish a theoretical head developed by the impeller. There are several versions of definition for the slip factor. However, the following definition in [24] is accepted in this paper: where the slip velocity of fluid (Figure 2) depends on the impeller geometry and flow rate as well as fluid viscosity, and so forth, and , is the relative tangential velocity with slip, whilst is the relative tangential velocity without slip. Traditionally, the slip factor of a centrifugal pump impeller is frequently estimated by means of the Wiesner or Stodola formulas. The Wiesner expression is read as in [25] and the Stodola’s relation is written as [25] where is the impeller blade discharge angle. Note that those equations do not involve the effects of both fluid viscosity and flow rate on the slip factor.

A slip factor can be calculated by using a velocity triangle at the impeller discharge in [4]. In addition, it was proposed that the slip factor can be estimated by using the Euler equation of turbomachinery with an actual relative flow angle at the impeller exit given by CFD computations in [19]. In that case, the slip factor is defined by the ratio of the impeller theoretical head over the head of the impeller with infinite number of blades. In that proposal, it was assumed that the 3D viscous and 1D ideal flow share the same radial flow velocity component. In experiments in [3], however, 3D viscous fluid flow showed an increasing radial velocity at the impeller exit compared to 1D uniform inviscid flow sketched in Figure 3. The same effect is also observed in CFD computational results in the paper. Thus, we have to apply the velocity triangle in Figure 3 to extract a slip factor.

For a uniform 1D ideal flow, the tangential component of the relative velocity is given by the following expression: where is the impeller tip speed, , the impeller diameter, the volumetric efficiency, = 0.869 based on the experiments in [26, 27], the impeller exit area, , the blade exit width, the blade blockage coefficient at exit, , the tangential thickness of blade at exit, and  mm.

Based on the CFD computational results, the velocity components , , and are available; thus, the slip velocity can be calculated with the following expression:

With the slip factor definition (1), the slip factor extracted from the velocity triangles is expressed by

It is important to point out that for 3D viscous fluids the slip factor estimated based on the velocity triangle at the impeller outlet constructed with CFD outcomes does not seem to be linked directly to the impeller theoretical head because the variable at the impeller outlet in the Euler equation no longer represents the total energy of the fluids; the details can be found in Section 4.1.

Additionally, we are able to evaluate a slip factor by using the impeller theoretical head of viscous fluids obtained from CFD computations. The impeller theoretical head of 1D fluid flow is written by

Substituting the theoretical head of viscous fluids given by CFD computations into (7) for , then a slip factor corresponding to is calculated by the following equation:

3. Results

3.1. Velocity Validation

For simplicity, the fluid relative velocity vectors measured by using LDV in [27] and predicted by CFD here are shown in Figure 4 for water (1cSt) and oil (48cSt) at BEP and Ra = 50 μm. The vectors lay in the plane that is through the mid-span of blade at trailing edge and perpendicular to the pump shaft. For the experimental and computational velocities, the velocity profile shows a low velocity near the pressure side of blade and a high velocity near the suction side at a small radius. However, as the radius increases, the velocity profile exhibits a high velocity near the pressure side and a low velocity near the suction side.

It is seen that the predicted flow fields of oil appear to show slight difference from water at the same working condition. Importantly, the computed flow fields are insensitive to the change in fluid viscosity. However, this is not true for the measured flow fields of oil.

The 2D, time-averaged, circumferential components of liquid velocity in the volute throat, section VIII, see Figure 1(b), given by LDV measurement in [28] and CFD computation, are shown in Figure 5 at BEP and Ra = 50 μm for water (1cSt) and oil (48cSt). In Figures 1(b) and 1(c), is the width of volute, the distance between the impeller outlet and the observation point, the distance from the one side wall of volute to the observation point, and the circumferential component of fluid velocity in the volute throat. Compared with the velocity profiles measured, the computed velocity distributions show considerable change along the spanwise direction () of volute near the impeller outlet and fall off substantially in the radial () direction. The velocity on the line of = 2 mm agrees well with the LDV observation. Note that the predicted velocity profile of viscous oil slightly differs from that of water once again.

Because the flows in the side chambers between the impeller outside surface and the volute casing are omitted in CFD computations, the effect of impeller outside surfaces driving fluid particles to rotate in the chambers is excluded. As a result, the computational velocity shows a quick decline in the radial direction and along the volute span.

3.2. Slip Factor

The slip factors extracted by using the foregoing methods are shown as a function of flow rate in Figure 6 at four viscosities and the surface roughness of Ra = 50 μm. The slip factor and largely depend on flow rate rather than viscosity. and show very different relationships with flow rate. The slip factors is reduced with increasing flow rate. This characteristic is the same as those revealed in the experiments conducted in [2, 3]. The slip factor approaches asymptotically the Wiesner slip factor when flow rate is more than 6.0 L/s.

The slip factor estimated by using the impeller theoretical head of 3D viscous flow is close to the Stodola slip factor at BEP. Moreover, the slip factor exhibits a large value at high and low flow rates.

The slip factors and are plotted against the Reynolds number (= , is the kinematic viscosity of liquid) in Figure 7 to show their dependence on liquid viscosity. The dependence of on viscosity is more evident than . Even though is less viscosity dependent, it still declines with increasing viscosity (reducing Re).

In Figure 7, the Stodola, Wiesner, Pfleiderer, and Ohta slip factors have been presented as well. The Pfleiderer slip factor is not a slip factor but a head correction, which is related to the slip factor defined in (1) by the following expression: where and expressions are taken from [29], and the first coefficient in is chosen to be 0.6. is the radius of blade leading edge at the mid-span, and is the radius of blade trailing edge at the mid-span too, .

The Ohta slip factor is totally empirical based on the flow fields behind 20 semiopen impellers measured by means of three-hole probe when the pump handling water and viscous oils and is given by the following expressions [17]:

In Figure 7, the Stodola slip factor is nearly the same as the Pfleiderer factor but is slightly larger than the Wiesner factor. The Ohta slip factor is valid just as [17]. When is higher than , the Ohta slip factor seems too large compared with the rest slip factors; further it depends largely on Reynolds number.

3.3. Impeller Theoretical Head and Hydraulic Efficiency Curves

The impeller theoretical head hydraulic efficiency are demonstrated in Figure 8 in terms of flow rate at the kinematical viscosities of 1cSt, 24cSt, 48ct, and 60cSt and with the wetted surface roughness of Ra = 50 μm. For comparison, the theoretical head of one-dimensional inviscid fluid flow derived from the Euler’s equation for turbomachinery and the Wiesner and Stodola slip factors are presented too.

The impeller theoretical head of 3D viscous flow significantly differs from 1D inviscid flow. The impeller theoretical head of viscous flow no longer increases linearly with decreasing flow rate; as a result of this, a maximum head occurs at the flow rate of 4.0 L/s (around 57% the flow rate at BEP) instead, causing a stall operational condition emerges there. Such a theoretical head variation characteristic agrees well with that in the experiments carried out by Acosta and Hollander [30] and Peck [31].

The impeller theoretical head of 1D inviscid flow based on the Stodola slip factor is in good agreement with those of 3D viscous flow given by CFD, especially nearly at = 6.0 L/s. Note that the head based on the Wiesner factor is overpredicted substantially, suggesting the Stodola’s slip factor is more suitable for the impellers with less number of blades and small exit blade angle.

The impeller theoretical head is viscosity dependent. Except at a viscosity of = 24cSt, the theoretical head of 3D viscous flow is degraded with increasing viscosity. At = 24cSt, however, the situation is quite different. In that case, the theoretical head is higher than water ( = 1cSt), causing the sudden-rising head effect, as flow rate is higher than 4.0 L/s.

The impeller hydraulic efficiency achieves maximum values ranged in 0.75–0.9 at = 6.0−7.0 L/s. Further the impeller hydraulic efficiency is better at = 24cSt when the flow rate is more than 5.0 L/s.

4. Discussions

4.1. Euler Head against Total Energy Head

At the outlet of impeller, we have the following expression for fluid energy head by using the cosine rule in the velocity triangle ΔABC shown in Figure 2 without considering elevation potential energy of fluid [32]

Considering the relations and , then (11) takes the following form: Rearranging the equation The term is the relative total head or rothalpy of fluid and is the Euler head at the outlet; consequently (13) can be rewritten as

Likewise, at the inlet of impeller, the fluid energy head can be expressed as the following: where is the rothalpy at the inlet of impeller. Based on (14) and (15), the relation between the Euler head rise and the total head rise can be expressed as where the fluid energy head rise, that is, impeller theoretical head . For an ideal (inviscid and incompressible) fluid flow, there are no hydraulic losses in the impeller at all, suggesting and ; for a viscous fluid flow, however, and .

, , and are shown in Figure 9 in terms of flow rate at the viscosities of 1cSt and 24cSt. Clearly, differs from in the figure. Moreover is somehow less than zero as the flow rate is lower than 5.0 L/s. Just beyond that flow rate becomes positive. In fact, represents the hydraulic losses in the impeller and should be positive at all the flow rates. Unfortunately, it is really hard to understand that FLUENT somehow gives a negative , which was observed in [20] as well. At the rest viscosity, the same situation remains.

Since the Euler head is different from the energy head for a viscous fluid, a slip factor from velocity triangle must be different from that based on the theoretical head of impeller even at design or BEP condition. The experimental data in [4, 17] confirmed that there is a difference in the two slip factors.

4.2. Effect of Exit Blade Angle on Slip Factor

The flow field of water and viscous oils in the pump with 44° exit blade angle was computed to identify the effect of the blade angle on slip factor. Except the exit blade angle and wrap angle, the rest of geometrical parameters of the impeller and volute were kept the same as the pump with 20° exit blade angle. Two slip factors extracted from the CFD results are illustrated in Figure 10 as a function of flow rate and in terms of Reynolds number at = 6.0 L/s (BEP), respectively.

Compared with the slip factors shown in Figure 6, the two slip factors of the impeller with larger exit blade angle are less dependent on flow rate. The Wienser rather than Stodola slip factor seems more close to CFD prediction. In Figure 10(c), the slip factors estimated with the Stodola, Wiesner, Pfleiderer, and Ohta correlations are presented as well. The Stodola and Pfleiderer slip factors seem overpredicted.

Compared with Figure 7, a larger blade angle leads to an increased slip factor. Moreover a higher viscosity can reduce slip factor considerably for an impeller with large exit blade angle.

The relative velocity vectors in the mid-span plane of the impellers with exit blade angle 20° and 44° at low flow rate ( = 3.4 L/s) and BEP for the viscosities of 1cSt and 48cSt are presented in Figure 11. In the impeller with large blade angle ( = 44°) there are clear vortices in the flow passages even at the BEP. In the impeller with small blade angle ( = 20°) the vortices emerge just at the low flow rate. A vortex increases the curvature of relative flow stream lines, causing a severe slip effect exhibits at the impeller discharge and an increased slip factor.

Since the vortices always stay in the passages of the impeller with large blade angle ( = 44°), the slip factor is less sensitive to flow rate and it is higher than that in the impeller with small blade angle ( = 20°) at BEP.

4.3. Reason for Degraded Euler Head at Low Flow Rate

In Figure 9, the Euler head is lower than the total energy head when < 5.0 L/s. To clarify the reasons behind this, the impeller static pressure head rise (= ) and kinetic head rise (= ) versus flow rate are illustrated in Figure 12 at various liquid viscosities and = 50 μm. It can be seen that the impeller static pressure head rise increases continuously with decreasing flow rate; contrarily, the kinetic head rise achieves a maximum at around = 4.0 L/s, then exhibits a dramatic decrease. Obviously, this too low kinetic head at the low flow rate is responsible for the lower Euler head at underloading condition.

To disclose the root-cause for the small kinetic head at a low flow rate, the computed tangential velocity of liquid at the impeller outlet is illustrated in Figure 12(c) as well. The velocity component has been made dimensionless with . It is clear that the tangential velocity is dramatically decreased at low flow rate. Note that the LDV measurements also revealed this effect in [28].

The reduction in tangential velocity implies that the volute has applied a strong influence on the flow around the impeller. To confirm this, the flow field of water in the impeller without a volute was computed and the tangential velocity at the impeller outlet is shown in Figure 12(d). The tangential velocity increases steadily with declining flow rate. This effect was observed in the experiment conducted by Kurokawa and Hode [33], where there was a vaneless diffuser rather than a volute behind a centrifugal impeller. Hence, it is the volute that is responsible for the reduction of tangential velocity at impeller discharge at a low flow rate.

5. Conclusions

Slip factor was extracted from the CFD computational outcomes by using two ways; one is from the velocity triangles at the impeller discharge and one is from the impeller theoretical head of 3D viscous turbulent flow. The impeller theoretical head against flow rate was compared between traditional 1D inviscid flow and 3D viscous one. The CFD computed velocity profiles are basically in agreement with LDV observations. Clearly, the slip factors given by the two methods are not consistent in numerical value. Flow rate rather than liquid viscosity show a remarked effect on the two slip factors. Volute is responsible for reduction of Euler head at low flow rates. The slip factor of impeller with large exit blade angle is less sensitive to change in flow rate because of existing vortices.