# Effect of Blade Outlet Angle on Unsteady Hydrodynamic Force of Closed-Type Centrifugal Pump with Single Blade

Authors: Yasuyuki Nishi; Junichiro Fukutomi
Journal: International Journal of Rotating Machinery
 (2014)

---

## Abstract

Geometrically, the single-blade centrifugal impeller, commonly used today as a sewage pump, is not axially symmetric. For this reason, the static pressure around the impeller fluctuates greatly when the impeller is rotating, and not only the radial thrust but also the axial thrust shows large fluctuations. Therefore, it is extremely important for the improvement of pump reliability to quantitatively grasp these fluctuating hydrodynamic forces. In this study, we investigated the unsteady hydrodynamic forces in a closed-type centrifugal pump with a single blade for different blade outlet angles using a numerical analysis that takes into account both experiment and the leakage flow. The results clearly showed the effect of the blade outlet angle on that act on the impeller. The root-mean-square value of the fluctuating component of the total radial thrust was roughly the same for whichever impeller at low flow rate, but at high flow rates, the value increased for impellers with larger blade outlet angles. Moreover, when the leakage flow rate increased with increasing static pressure around the impeller, such that the rear and front shroud parts were subject to high pressure, the absolute value of the axial thrust on both these parts increased.


---

## Body

## 1. Introduction

Because a sewage pump is used for the transfer of liquids containing foreign bodies, such as solid and fibrous material, it must not only be highly efficient but also allow the easy passage of foreign bodies. The passed particle size relative to the pump bore (the minimum particle size in the flow path inside the pump) exists as an index of the ease with which these foreign bodies can pass through [1]. It therefore follows that the single-blade centrifugal pump, consisting of just a single blade of broad blade width, has been widely used for this type of pump so that a large passed particle size is formed. Nevertheless, even for centrifugal pumps with many blades, geometrical manufacturing errors cause disproportional fluid force [2]. In contrast, with the single-blade centrifugal pump, by its very nature, the impeller itself is not axially symmetric from a geometrical viewpoint. Additionally, because a volute casing is installed on the outside of the impeller that is not axially symmetric, the static pressure distribution surrounding the impeller is extremely nonuniform, and an excessively fluctuating radial thrust acts on the impeller [3, 4]. Because such excessively fluctuating radial thrust can lead to a drop in the axial strength, a shortening of the lifetime of the bearings and shaft sealing equipment, and the generation of axial vibrations, elucidating the generation mechanism behind this has become an important challenge.Regarding the radial thrust in single-blade centrifugal pumps, such subjects have been investigated as the effect of the number of blades [3, 5], winding angle of the blade [3], distribution of the blade angle [6] and its relationship with the static pressure between the blades [4], effect of the casing shape [3, 6], effectiveness of adding extra mass [4, 7], and relationship to vibration of the casing [7]. In addition, the relationship between the hydrodynamic forces acting on the impeller and the whirling of the impeller has been investigated [5, 8–10]. We have also proposed a design method for single-blade centrifugal impellers and have discovered, for a closed-type centrifugal impeller with a single blade that was designed using this design method as a basis, how the performance characteristics relate to the radial thrust [11]. While modeling the components of the radial thrust using an unsteady conservation-of-momentum law to illustrate their behavior [12], we also explained how they were affected by the blade outlet angle [13]. However, hardly any research has been made on the axial thrust. In our previous studies, the hydrodynamic forces acting on the rear and front shroud parts were neglected in the numerical analyses. Accordingly, by performing a numerical analysis that takes into account the leakage flow on the rear surface of the rear and front shroud parts, we were able to demonstrate that not only did the radial thrust fluctuate with fluctuations in the static pressure distribution around the impeller but also the axial thrust fluctuated greatly [14]. Furthermore, modeling the components of the axial thrust, we showed how they behaved [15].In this study, we implemented a numerical analysis that takes into account the leakage flow on the rear surface of the rear and front shroud parts of a closed-type single-blade centrifugal pump for different blade outlet angles. And with a comparative investigation, we could clarify how the blade outlet angle affects both the radial and axial thrust in the closed-type single-blade centrifugal pump.


## 2. Experimental Apparatus and Method

In this study, we used three types of single-blade centrifugal impellers with respective blade outlet anglesβb2 of 8°, 16°, and 24°. However, in the experiments, we used only impellers with βb2 = 8 and βb2 = 16. An overview of the test impeller (βb2 = 8) is shown in Figure 1, and its specifications are given in Table 1. An overview of the test pump is also given in Figure 2. The impeller with βb2 = 8 had a passed particle size of 56 mm and had been designed on the basis of the design method reported previously [11]. This design method stipulates the passed particle size beforehand, and the impeller design is implemented using, as design parameters, first, the relative velocity ratio of the blade inlet and blade outlet and, second, the absolute velocity ratio of the impeller suction inlet and the blade inlet. The impellers with βb2 = 16 and βb2 = 24 have meridian planes that are equivalent to that for the impeller with βb2 = 8, and, as shown in Table 1, only the blade outlet angle is varied. This pump has three narrow gaps between the impeller and the casing, as seen in Figure 2. These are the liner-ring gap (the gap with the casing wearing ring), the gap with the rear shroud part, and the gap with the front shroud part. The gap widths are, respectively, 0.75, 0.5, and 0.5 mm and the corresponding gap lengths are 17, 10.1, and 10.5 mm. It should be noted that, because a mechanical tool has been adopted for the shaft sealing equipment in this pump, almost no leakage to the exterior occurs on the rear shroud part side.Table 1
Specifications of test impellers.
Impeller suction diameterD0
67 mm
Impeller inner diameterD1
58 mm
Impeller outer diameterD2
194 mm
Blade inlet widthb1
56 mm
Blade outlet widthb2
56 mm
Blade inlet angleβb1
13°
Blade outlet angleβb2
8°
16°
24°Figure 1
Test impeller.Figure 2
Test pump.Figure3 gives an overview of the experimental apparatus. The experiment was conducted at a rotational speed of n=1740 min−1, and the flow rate Q was adjusted with the gate valve mounted on the discharge side and measured using an electromagnetic flow rate meter. The total head H was found from the static pressure measured before and after the pump with a strain-gauge-type pressure transducer and flow rate. The shaft power L was found by measuring the rotational speed n and torque T with a magnetoelectric-type rotation detector and torque detector, respectively.Figure 3
Experimental apparatus.The coordinate system is defined by thex-, y-, and z-axes, as shown in Figure 2. The blade phase angle θ0 is taken in the anticlockwise direction from the x-axis, and the point where the end of the blade outlet passes through the x-axis is assumed as θ0 = 0° (the position shown in Figure 2). Eight monitor points were provided on the wall surface on the side of the front shroud part at positions 10 mm radially outwards from the impeller outlet in an anticlockwise direction from the x-axis as θ = 0°. The pressure measurement was synchronized with the rotation of the impeller using the torque signal for the θ0 = 0° position of the end of the blade outlet as a reference.For the radial thrust measurement, anX-Y load cell was incorporated in the periphery of the bearings, as shown in Figures 4 and 5, and it was synchronized with the rotation of the impeller in the same way as that for the pressure measurement [11–15]. The “four-active-gauge” method was used for the X-Y load cell by attaching the connected strain gauges together. The dynamic calibration method [16] was adopted for calibration. The relationship between the output voltage and the force was tested in air for a known centrifugal force acting on a calibrating circular plate, which has a mass found by subtracting from the impeller mass; that mass corresponds to the buoyancy force that normally acts on it. The calibration curves for the X- and Y-directions when the force acts in the Y-direction are shown in Figure 6. The output voltage and the force are proportionally related. Moreover, when the force acts in the Y-direction, its influence on the X-direction is small. It should be noted that the dynamic balance was taken between the impeller under test and calibrating circular plate. The effect of the heat from the bearings was eliminated using a method where the pump was stopped over a short interval of around 15 s under thermally stable conditions and was operated immediately afterwards, adjusting the zero point to take measurements under the same conditions.Figure 4
X-Y load cell.Figure 5
X-Y load cell incorporated in periphery of bearings.Figure 6
Calibration curves forX- and Y-directions when force acts in the Y-direction.


## 3. Numerical Analysis Method

For the numerical analysis, we performed a three-dimensional unsteady-flow analysis using the general-purpose thermofluid analysis code ANSYS CFX13.0. Reynolds-averaged Navier-Stokes equations were used as a basis, and the SST (shear stress transport) model was used to model the turbulent flow.The computational domain and computational grid for the pump part are shown, respectively, in Figures7(a) and 7(b). This study differs from past research [11–13] in that the gaps between the outer surfaces of the rear and front shroud parts of the impeller and the casing are actually modeled. The computation domain consists of the suction pipe, impeller, volute casing, regions in the gaps with the rear and front shroud parts, and discharge pipe. The computational grids for the impellers with βb2 = 8, βb2 = 16, and βb2 = 24 have 690,000, 705,000, and 707,000 elements, respectively, where the total number of elements in the grids is 2,328,000, 2,343,000, and 2,345,000, respectively. Here, the computational grids for the suction pipe, volute casing, regions in the gaps with the rear and front shroud parts, and the discharge pipe are identical in all cases, with, respectively, 149,000, 462,000, 340,000, 572,000, and 115,000 elements. The number of elements in the gaps counted radially is five for the liner-ring gap and four for both the rear shroud part and front shroud part gaps. A nonstructured tetrahedral grid is used as the computational grid for the impeller and volute casing, whereas a nonstructured hexahedral grid is used for the others.Computational model.
(a)
Computational domain
(b)
Computational grids of pump partTo study the dependence of the computational model on the grid used, we performed analyses for theβb2 = 8 case, where there were, respectively, 8, 6, and 6 elements in the radial direction in the liner-ring gap, rear shroud part gap, and front shroud part gap, and computational grids for the rear shroud part gap and front shroud part gap regions had 561,000 and 1,017,000 elements, respectively, with a total grid number of approximately 2,994,000 elements. The results showed that the errors in the flow rate, at the maximum efficiency point in the analysis results for the βb2 = 8 case, were 0.2% for the head, 0.3% for the shaft power, 0.7% for the root-mean-square value of the total radial thrust, and 4.1% for the time-averaged value of the total axial thrust. For an impeller with βb2 = 8, we performed analyses with volute casing computational grid numbers of 1,035,000 and 693,000 elements, with a total grid number of approximately 2,904,000 elements. The results showed that the errors in the flow rate, at the maximum efficiency point in the analysis results for the βb2 = 8 case, were 1.8% for the head, 1.6% for the shaft power, 0.6% for the root-mean-square value of the total radial thrust, and 3.7% for the time-averaged value of the total axial thrust. As mentioned above, we verified that the grid dependence for this computational model was comparatively small [15].The boundary conditions for the suction pipe inlet were set as those on the inlet boundary and for the discharge pipe as those on the outlet boundary and were given as the mass flow rate on the inlet boundary, a static pressure (gauge pressure) on the outlet boundary of 0 Pa, and the rotational speed in the impeller domain. The wall surface was assumed to be a nonslip condition, and the boundary between the rotating and static domains was conjugated using the transient rotor-stator technique [17]. The time step was set such that 360 steps were a single rotation of the impeller, and the calculations were continued until the fluctuations in the flow became practically constant.


## 4. Experimental and Analysis Results and Considerations

### 4.1. Comparison of Pump Performances

The pump performance for each different impeller based on both experiment and numerical analysis is shown in Figures8(a)–8(c). Qualitatively, the experimental and calculated values for the cases βb2 = 8 and βb2 = 16 showed good agreement for the head coefficient ψ, shaft power coefficient λ, and pump efficiency η. The calculated values of ψ are seen to increase in impellers for which the βb2 increases with increasing flow rate. With the flow rate at the point of maximum efficiency, φ = 0.030, ψ for βb2 = 16 and βb2 = 24 increased by 13.7% and 33.8%, respectively, relative to βb2 = 8. This could be attributed to the increase in the circumferential component of the absolute velocity of the blade outlet as the blade outlet angle increases, whereas, for the calculated values of η, the flow rate at the point of maximum efficiency with βb2 = 8 is φ = 0.030 and that for βb2 = 16 and βb2 = 24 migrated to the larger side is φ = 0.039. The η for βb2 = 16 and βb2 = 24 at φ = 0.039 was 0.735 and 0.767, respectively, and we find that η increases for impellers with larger βb2. The flow rate at the point of maximum efficiency is thought to be moving toward a larger flow rate because the theoretical pump head for a finite number of blades increases as the blade outlet angle increases, and the matching point with Worster’s characteristic formula for the volute [18] then moves toward a larger flow rate.Performance curves.
(a)
Head coefficient
(b)
Shaft power coefficient
(c)
Pump efficiency


### 4.2. Comparison of Radial Thrust

In this study, the impeller surface is divided into three regions, as shown in Figure9. The blade surface is assumed to be the hydraulic part, the outer surface of the rear shroud part including its outer periphery as the rear shroud part, and the outer surface of the front shroud part including its outer periphery as the front shroud part. The total radial thrust Fr acting on the impeller is expressed by the following formula [14, 15]:
(1)Fr=FrH+FrR+FrF.
Here, FrH is the radial thrust on the hydraulic part [N], FrR is the radial thrust on the rear shroud part [N], and FrF is the radial thrust on the front shroud part [N].Figure 9
Impeller surfaces.The radial thrustFr can be defined by the following formula as the vector sum of the time-averaged value (an average component) and a fluctuating component [3, 4, 16]:
(2)Fr=F-r+ΔFr.
Here, F-r is the time-averaged value and ΔFr is the fluctuating component.The time-averaged valuef-r and the root-mean-square value of fluctuating component Δfr of the total radial thrust coefficient on the impeller from both experiment and numerical analysis are shown in Figure 10. The calculated value for radial thrust was found by integrating, over each of the impeller parts, the pressure acting on the wall surface and the shear stress on the wall surface. From Figure 10, the experimental and calculated values for f-r and the root-mean-square value of Δfr are seen to agree qualitatively. Although all the impellers show roughly the same values for the root-mean-square value of Δfr with a small flow rate of φ = 0.008, at large flow rates, it increases for impellers with larger βb2. This is thought to be because, although it also depends on the load distribution on the blade, for large flow rates, there is a relative increase in the load on the blade for impellers with larger βb2. On the other hand, with f-r, it increases for either a low flow rate or a high flow rate for whichever of the impellers, and there are comparatively small differences due to changes in βb2. This is believed to be due to the circumferential distribution of the time-averaged pressure coefficient CP at the monitor points in the impeller outlet shown in Figure 11, which are the results of numerical analysis; CP gets higher for impellers with larger βb2 for any of the flow rates, and that value increases uniformly. Thus, we have confirmed that there is good agreement between the experimental and calculated values of CP at the impeller outlet [15].Figure 10
Time-averaged values of total radial thrust.Figure 11
Time-averaged static pressure distributions (cal.).Subsequently, we investigate the behavior of the radial thrust acting on each part of the different impellers. The Lissajous figures for the radial thrust coefficient for the hydraulic part, rear shroud part, and front shroud part of each impeller are shown in Figures12(a)–12(c) for φ = 0.030. For each radial thrust coefficient, the values at θ0 = 0, 90, 180, and 270 are shown, respectively, by ●, ▲, ■, and ◆, and the time-averaged values are given using the same symbols as those for the radial thrust coefficient. These are represented similarly in Figures 13(a) and 13(b). From Figures 12(a)–12(c), the time-averaged value of the radial thrust acting on each part is located near the origin for any of the impellers. The fluctuating component of the radial thrust for both the hydraulic and rear shroud parts increases for impellers with larger βb2. However, in the case of radial thrust for the front shroud part, it is smaller for the impeller with βb2 = 24 than for the one with βb2 = 16. Because a leakage flow exists on the side of the front shroud part, the effect of an increase in the leakage flow rate for the impeller with βb2 = 24 (discussed later) can reasonably be considered as the reason for this finding. Moreover, we find that the radial thrust for the rear shroud and front shroud parts is small in relation to the radial thrust for the hydraulic part for all the impellers. It therefore follows that the radial thrust for the hydraulic part dominates the total radial thrust.Lissajous figures of radial thrust which acts on each part (φ = 0.030, cal.).
(a)
Hydraulic part
(b)
Rear shroud part
(c)
Front shroud partLissajous figures of hydraulic part radial thrust (cal.).
(a)
φ = 0.008
(b)
φ = 0.047Thus, concentrating on the radial thrust coefficient for the hydraulic part, its Lissajous figure is shown in Figures13(a) and 13(b) for the low flow rate φ = 0.008 and high flow rate φ = 0.047. With the low flow rate φ = 0.008, the radial thrust for the hydraulic part is large in the vicinity of θ0 = 180°, and it is largest for the impeller with βb2 = 24. On the contrary, the radial thrust for the hydraulic part is small in the vicinity of θ0 = 0°, and it is smallest for the impeller with βb2 = 24.On the other hand, with the high flow rateφ = 0.047, the radial thrust for the hydraulic part is large in the vicinity of θ0 = 0°, and it is largest for the impeller with βb2 = 24. Conversely, the radial thrust for the hydraulic part is small in the vicinity of θ0 = 180° and is smallest for the impeller with βb2 = 8. As shown with the numerical analyses in Figures 14(a)–14(d), 15(a)–15(d), 16(a)–16(d), and 17(a)–17(d), from the static pressure distribution on the impellers with βb2 = 8 and βb2 = 24, at both φ = 0.008 and φ = 0.047, the following observations can be made. The cross-section is in the center of the blade width. At the low flow rate φ = 0.008, the time-averaged CP increases as the end of the volute winding is approached. Moreover, if θ0 = 180°, because CP for a pressure surface with 180° < θ < 360° increases more, the radial thrust for the hydraulic part acts in the +y-axis direction. For a large blade outlet angle of βb2 = 24, both the CP for the pressure surface and time-averaged CP at the end of the volute winding are the largest, and the radial thrust for the hydraulic part in that case is the largest. However, if θ0 = 0°, the CP for a pressure surface with 0° < θ < 180° is increased, and the CP for βb2 = 24 is the largest. However, owing to the effect of the time-averaged CP, which increases as the end of the volute winding is approached, the time-averaged value of the radial thrust for the hydraulic part acts largely in the +y-axis direction. Therefore, the radial thrust for the hydraulic part can be considered to act in the –y-axis direction in terms of its time-averaged value. For βb2 = 24, the time-averaged CP is the largest, and the radial thrust for the hydraulic part is the smallest. On the other hand, with a flow rate φ = 0.047, the time-averaged CP increases as the beginning of the volute winding is approached. Moreover, if θ0 = 0°, because the CP for a pressure surface with 0° < θ < 180° increases more, the radial thrust for the hydraulic part acts in the -Y direction. For βb2 = 24, both the CP for the pressure surface and time-averaged CP at the beginning of the volute winding are the largest, and the radial thrust for the hydraulic part is the largest. On the contrary, if θ0 = 180°, CP for a pressure surface with 180° < θ < 360° is increased. However, due to the effect of the time-averaged CP, which decreases as the end of the volute winding is approached, the time-averaged value of the radial thrust for the hydraulic part acts largely in the −y-axis direction. Therefore, the radial thrust for the hydraulic part can be considered to act in the +y-axis direction in terms of its time-averaged value. For βb2 = 8, both the CP for the pressure surface and time-averaged CP at the end of the volute winding are the smallest, and the radial thrust for the hydraulic part can then be considered to have its smallest value.Static pressure distributions at blade width central (φ = 0.008, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at blade width central (φ = 0.008, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at blade width central (φ = 0.047, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at blade width central (φ = 0.047, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°


### 4.3. Comparison of Axial Thrust

Similar to the radial thrust in the previous section, the total axial thrustFa acting on the impeller is expressed by the following formula [14, 15]:
(3)Fa=FaH+FaR+FaF.
Here, FaH is the axial thrust for the hydraulic part [N], FaR is the axial thrust for the rear shroud part [N], and FaF is the axial thrust for the front shroud part [N].The time-averaged values of the axial thrust coefficient for the hydraulic part, rear shroud part, and front shroud part are shown for each of the impellers in Figures18(a) and 18(b). The calculated value for the axial thrust was found in a similar manner to that for the radial thrust, where the pressure acting on the wall surfaces and the shear stress on the wall surfaces is integrated overall into each part of the impeller. Comparing the absolute values of the axial thrust coefficient for each part and for all the impellers, the axial thrust coefficient faH for the hydraulic part is the smallest. The axial thrust coefficients for the rear shroud part faR and front shroud part faF are dominant, unlike the cases for the radial thrust. faR and faF are in opposing directions relative to each other, and their absolute values increase for smaller flow rates. Moreover, faR is larger than faF, from which it follows that the total axial thrust coefficient fa acts in the −z-axis direction and its absolute value increases for smaller flow rates.Time-averaged values of axial thrust which acts on each part (cal.).
(a)
Total and hydraulic part
(b)
Rear and front shroud partsComparing the individual impellers, the absolute value of the total axial thrust coefficientfa is the largest for βb2 = 24. Little difference exists between βb2 = 8 and βb2 = 16; however, for small flow rates, the difference is slightly larger for βb2 = 8. The larger the βb2 for the impeller is, the more the absolute values of faR and faF increase at high flow rates. However, when the blade outlet angle βb2 is varied from 8° to 16°, the increase in the absolute value of faR is relatively small compared with the increase in the absolute value of faF, and for φ = 0.008 in particular, the absolute values of faR for βb2 = 8 and βb2 = 16 approximately agree.Concentrating on the total axial thrust and the axial thrust on the rear and front shroud parts forφ = 0.008 and φ = 0.047, the fluctuation of fa, faR, and faF during a single rotation of the impeller is shown in Figures 19(a)–19(c). From Figure 19(a), the axial thrust coefficient fa for all the impellers shows a similar behavior and fluctuates greatly during a single rotation of the impeller. The absolute value of fa is largest in the vicinity of θ0 = 90°–180° for a low flow rate φ = 0.008, but conversely, for a high flow rate φ = 0.047, the largest value is at θ = 0°(360°). From Figures 19(b) and 19(c), the absolute values of faR and faF also behave in a manner similar to the absolute value of fa. Therefore, for all the impellers, the total axial thrust is partially offset by the axial thrust of the front shroud part, but we see a behavior roughly similar to the axial thrust for the rear shroud part.Axial thrust fluctuation which acts on each part (cal.).
(a)
Total
(b)
Rear shroud part
(c)
Front shroud partFor impellers withβb2 = 8 and βb2 = 24, we investigated the axial thrust for the rear and front shroud parts in relation to the pressure field. For a low flow rate φ = 0.008, the static pressure distribution in the clearance domain for the rear and front shrouds for the βb2 = 8 case is shown in Figures 20(a)–20(d) and 21(a)–21(d), respectively, whereas the same pressure distributions for the case βb2 = 24 are shown in Figures 22(a)–22(d) and 23(a)–23(d), respectively. Furthermore, the same pressure distributions are shown for a high flow rate φ = 0.047 in Figures 24(a)–24(d) and 25(a)–25(d) for βb2 = 8 and in Figures 26(a)–26(d) and 27(a)–27(d) for βb2 = 24. Here, the static pressure distribution in the clearance domain for the front and rear shrouds is for the cross-sections C-C and D-D, respectively, of Figure 3. For any of the impellers, the CP in the clearance domain of either the front or rear shroud fluctuates greatly during rotation of the impeller in the same way as the static pressure distribution in the interior of the pump, as shown in Figures 14–17. For a low flow rate φ = 0.008, the CP for the clearance domain of the front shroud takes a high value for all the impellers in the region θ0 = 90°–180°. Owing to the effect of leakage flow, the CP for the clearance domain of the front shroud drops overall below that for the clearance domain of the rear shroud but similarly takes a high value in the region θ0 = 90°–180°, which corresponds to the behavior of the axial thrust for the rear and front shroud parts, shown, respectively, in Figures 19(b) and 19(c). At this point, if we examine the static pressure distribution in the interior of the pump, as shown in Figures 14 and 15, the static pressure around the impeller is greater overall in the region θ0 = 90°–180° than in the region θ0 = 180°–270°. That is to say, we can reason that when the static pressure around the impeller increases overall, the static pressure in the clearance domain for the front and rear shrouds increases in accordance with the increase in the leakage flow rate. This is evident from the fact that if we consider the fluctuation in the leakage flow rate coefficient Cq from the liner ring as shown in Figure 28(a) with φ = 0.008, Cq for any of the impellers attains its maximum value in the vicinity of θ0 = 90°–180°, which is larger; the larger value is the βb2 for the impeller.Static pressure distributions at rear shroud (φ = 0.008, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at front shroud (φ = 0.008, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at rear shroud (φ = 0.008, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at front shroud (φ = 0.008, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at rear shroud (φ = 0.047, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at front shroud (φ = 0.047, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at rear shroud (φ = 0.047, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at front shroud (φ = 0.047, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Leakage flow rate fluctuation (cal.).
(a)
φ = 0.008
(b)
φ = 0.047On the other hand, for a large flow rateφ = 0.047, the CP for the clearance domain of the front and rear shrouds takes a high value at θ0 = 0°(360°), which corresponds to the behavior of the axial thrust for the rear and front shroud parts, as shown in Figures 19(b) and 19(c). At this point, looking at the static pressure distribution in the interior of the pump as shown in Figures 16 and 17, the static pressure around the impeller at θ0 = 0° is higher overall than in the region θ0 = 180°–270°. It therefore follows that, as shown in Figure 28(b), for φ = 0.047, Cq attains its maximum value in the vicinity of θ0 = 0°, which is larger; the larger value is βb2 for the impeller.From the above considerations, because the time-averagedCP at the end of the volute winding for low flow rates is large, in the region θ0 = 90°–180°, and because the time-averaged CP at the beginning of the volute winding for large flow rates is large, at θ0 = 0°, the high CP of the pressure side propagates in the circumferential direction and the CP around the impeller increases. We can assume that accompanying the increase in this CP, the leakage flow rate also increases, and owing to the high pressure acting on the outer surface of the rear and front shrouds, the axial thrust for the rear and front shroud parts therefore increases. Because the larger the blade outlet angle for the impeller, the higher the CP of the pressure side and time-averaged CP at either the beginning or end of the volute winding. We can argue that this tendency is very pronounced.


## 4.1. Comparison of Pump Performances

The pump performance for each different impeller based on both experiment and numerical analysis is shown in Figures8(a)–8(c). Qualitatively, the experimental and calculated values for the cases βb2 = 8 and βb2 = 16 showed good agreement for the head coefficient ψ, shaft power coefficient λ, and pump efficiency η. The calculated values of ψ are seen to increase in impellers for which the βb2 increases with increasing flow rate. With the flow rate at the point of maximum efficiency, φ = 0.030, ψ for βb2 = 16 and βb2 = 24 increased by 13.7% and 33.8%, respectively, relative to βb2 = 8. This could be attributed to the increase in the circumferential component of the absolute velocity of the blade outlet as the blade outlet angle increases, whereas, for the calculated values of η, the flow rate at the point of maximum efficiency with βb2 = 8 is φ = 0.030 and that for βb2 = 16 and βb2 = 24 migrated to the larger side is φ = 0.039. The η for βb2 = 16 and βb2 = 24 at φ = 0.039 was 0.735 and 0.767, respectively, and we find that η increases for impellers with larger βb2. The flow rate at the point of maximum efficiency is thought to be moving toward a larger flow rate because the theoretical pump head for a finite number of blades increases as the blade outlet angle increases, and the matching point with Worster’s characteristic formula for the volute [18] then moves toward a larger flow rate.Performance curves.
(a)
Head coefficient
(b)
Shaft power coefficient
(c)
Pump efficiency


## 4.2. Comparison of Radial Thrust

In this study, the impeller surface is divided into three regions, as shown in Figure9. The blade surface is assumed to be the hydraulic part, the outer surface of the rear shroud part including its outer periphery as the rear shroud part, and the outer surface of the front shroud part including its outer periphery as the front shroud part. The total radial thrust Fr acting on the impeller is expressed by the following formula [14, 15]:
(1)Fr=FrH+FrR+FrF.
Here, FrH is the radial thrust on the hydraulic part [N], FrR is the radial thrust on the rear shroud part [N], and FrF is the radial thrust on the front shroud part [N].Figure 9
Impeller surfaces.The radial thrustFr can be defined by the following formula as the vector sum of the time-averaged value (an average component) and a fluctuating component [3, 4, 16]:
(2)Fr=F-r+ΔFr.
Here, F-r is the time-averaged value and ΔFr is the fluctuating component.The time-averaged valuef-r and the root-mean-square value of fluctuating component Δfr of the total radial thrust coefficient on the impeller from both experiment and numerical analysis are shown in Figure 10. The calculated value for radial thrust was found by integrating, over each of the impeller parts, the pressure acting on the wall surface and the shear stress on the wall surface. From Figure 10, the experimental and calculated values for f-r and the root-mean-square value of Δfr are seen to agree qualitatively. Although all the impellers show roughly the same values for the root-mean-square value of Δfr with a small flow rate of φ = 0.008, at large flow rates, it increases for impellers with larger βb2. This is thought to be because, although it also depends on the load distribution on the blade, for large flow rates, there is a relative increase in the load on the blade for impellers with larger βb2. On the other hand, with f-r, it increases for either a low flow rate or a high flow rate for whichever of the impellers, and there are comparatively small differences due to changes in βb2. This is believed to be due to the circumferential distribution of the time-averaged pressure coefficient CP at the monitor points in the impeller outlet shown in Figure 11, which are the results of numerical analysis; CP gets higher for impellers with larger βb2 for any of the flow rates, and that value increases uniformly. Thus, we have confirmed that there is good agreement between the experimental and calculated values of CP at the impeller outlet [15].Figure 10
Time-averaged values of total radial thrust.Figure 11
Time-averaged static pressure distributions (cal.).Subsequently, we investigate the behavior of the radial thrust acting on each part of the different impellers. The Lissajous figures for the radial thrust coefficient for the hydraulic part, rear shroud part, and front shroud part of each impeller are shown in Figures12(a)–12(c) for φ = 0.030. For each radial thrust coefficient, the values at θ0 = 0, 90, 180, and 270 are shown, respectively, by ●, ▲, ■, and ◆, and the time-averaged values are given using the same symbols as those for the radial thrust coefficient. These are represented similarly in Figures 13(a) and 13(b). From Figures 12(a)–12(c), the time-averaged value of the radial thrust acting on each part is located near the origin for any of the impellers. The fluctuating component of the radial thrust for both the hydraulic and rear shroud parts increases for impellers with larger βb2. However, in the case of radial thrust for the front shroud part, it is smaller for the impeller with βb2 = 24 than for the one with βb2 = 16. Because a leakage flow exists on the side of the front shroud part, the effect of an increase in the leakage flow rate for the impeller with βb2 = 24 (discussed later) can reasonably be considered as the reason for this finding. Moreover, we find that the radial thrust for the rear shroud and front shroud parts is small in relation to the radial thrust for the hydraulic part for all the impellers. It therefore follows that the radial thrust for the hydraulic part dominates the total radial thrust.Lissajous figures of radial thrust which acts on each part (φ = 0.030, cal.).
(a)
Hydraulic part
(b)
Rear shroud part
(c)
Front shroud partLissajous figures of hydraulic part radial thrust (cal.).
(a)
φ = 0.008
(b)
φ = 0.047Thus, concentrating on the radial thrust coefficient for the hydraulic part, its Lissajous figure is shown in Figures13(a) and 13(b) for the low flow rate φ = 0.008 and high flow rate φ = 0.047. With the low flow rate φ = 0.008, the radial thrust for the hydraulic part is large in the vicinity of θ0 = 180°, and it is largest for the impeller with βb2 = 24. On the contrary, the radial thrust for the hydraulic part is small in the vicinity of θ0 = 0°, and it is smallest for the impeller with βb2 = 24.On the other hand, with the high flow rateφ = 0.047, the radial thrust for the hydraulic part is large in the vicinity of θ0 = 0°, and it is largest for the impeller with βb2 = 24. Conversely, the radial thrust for the hydraulic part is small in the vicinity of θ0 = 180° and is smallest for the impeller with βb2 = 8. As shown with the numerical analyses in Figures 14(a)–14(d), 15(a)–15(d), 16(a)–16(d), and 17(a)–17(d), from the static pressure distribution on the impellers with βb2 = 8 and βb2 = 24, at both φ = 0.008 and φ = 0.047, the following observations can be made. The cross-section is in the center of the blade width. At the low flow rate φ = 0.008, the time-averaged CP increases as the end of the volute winding is approached. Moreover, if θ0 = 180°, because CP for a pressure surface with 180° < θ < 360° increases more, the radial thrust for the hydraulic part acts in the +y-axis direction. For a large blade outlet angle of βb2 = 24, both the CP for the pressure surface and time-averaged CP at the end of the volute winding are the largest, and the radial thrust for the hydraulic part in that case is the largest. However, if θ0 = 0°, the CP for a pressure surface with 0° < θ < 180° is increased, and the CP for βb2 = 24 is the largest. However, owing to the effect of the time-averaged CP, which increases as the end of the volute winding is approached, the time-averaged value of the radial thrust for the hydraulic part acts largely in the +y-axis direction. Therefore, the radial thrust for the hydraulic part can be considered to act in the –y-axis direction in terms of its time-averaged value. For βb2 = 24, the time-averaged CP is the largest, and the radial thrust for the hydraulic part is the smallest. On the other hand, with a flow rate φ = 0.047, the time-averaged CP increases as the beginning of the volute winding is approached. Moreover, if θ0 = 0°, because the CP for a pressure surface with 0° < θ < 180° increases more, the radial thrust for the hydraulic part acts in the -Y direction. For βb2 = 24, both the CP for the pressure surface and time-averaged CP at the beginning of the volute winding are the largest, and the radial thrust for the hydraulic part is the largest. On the contrary, if θ0 = 180°, CP for a pressure surface with 180° < θ < 360° is increased. However, due to the effect of the time-averaged CP, which decreases as the end of the volute winding is approached, the time-averaged value of the radial thrust for the hydraulic part acts largely in the −y-axis direction. Therefore, the radial thrust for the hydraulic part can be considered to act in the +y-axis direction in terms of its time-averaged value. For βb2 = 8, both the CP for the pressure surface and time-averaged CP at the end of the volute winding are the smallest, and the radial thrust for the hydraulic part can then be considered to have its smallest value.Static pressure distributions at blade width central (φ = 0.008, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at blade width central (φ = 0.008, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at blade width central (φ = 0.047, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at blade width central (φ = 0.047, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°


## 4.3. Comparison of Axial Thrust

Similar to the radial thrust in the previous section, the total axial thrustFa acting on the impeller is expressed by the following formula [14, 15]:
(3)Fa=FaH+FaR+FaF.
Here, FaH is the axial thrust for the hydraulic part [N], FaR is the axial thrust for the rear shroud part [N], and FaF is the axial thrust for the front shroud part [N].The time-averaged values of the axial thrust coefficient for the hydraulic part, rear shroud part, and front shroud part are shown for each of the impellers in Figures18(a) and 18(b). The calculated value for the axial thrust was found in a similar manner to that for the radial thrust, where the pressure acting on the wall surfaces and the shear stress on the wall surfaces is integrated overall into each part of the impeller. Comparing the absolute values of the axial thrust coefficient for each part and for all the impellers, the axial thrust coefficient faH for the hydraulic part is the smallest. The axial thrust coefficients for the rear shroud part faR and front shroud part faF are dominant, unlike the cases for the radial thrust. faR and faF are in opposing directions relative to each other, and their absolute values increase for smaller flow rates. Moreover, faR is larger than faF, from which it follows that the total axial thrust coefficient fa acts in the −z-axis direction and its absolute value increases for smaller flow rates.Time-averaged values of axial thrust which acts on each part (cal.).
(a)
Total and hydraulic part
(b)
Rear and front shroud partsComparing the individual impellers, the absolute value of the total axial thrust coefficientfa is the largest for βb2 = 24. Little difference exists between βb2 = 8 and βb2 = 16; however, for small flow rates, the difference is slightly larger for βb2 = 8. The larger the βb2 for the impeller is, the more the absolute values of faR and faF increase at high flow rates. However, when the blade outlet angle βb2 is varied from 8° to 16°, the increase in the absolute value of faR is relatively small compared with the increase in the absolute value of faF, and for φ = 0.008 in particular, the absolute values of faR for βb2 = 8 and βb2 = 16 approximately agree.Concentrating on the total axial thrust and the axial thrust on the rear and front shroud parts forφ = 0.008 and φ = 0.047, the fluctuation of fa, faR, and faF during a single rotation of the impeller is shown in Figures 19(a)–19(c). From Figure 19(a), the axial thrust coefficient fa for all the impellers shows a similar behavior and fluctuates greatly during a single rotation of the impeller. The absolute value of fa is largest in the vicinity of θ0 = 90°–180° for a low flow rate φ = 0.008, but conversely, for a high flow rate φ = 0.047, the largest value is at θ = 0°(360°). From Figures 19(b) and 19(c), the absolute values of faR and faF also behave in a manner similar to the absolute value of fa. Therefore, for all the impellers, the total axial thrust is partially offset by the axial thrust of the front shroud part, but we see a behavior roughly similar to the axial thrust for the rear shroud part.Axial thrust fluctuation which acts on each part (cal.).
(a)
Total
(b)
Rear shroud part
(c)
Front shroud partFor impellers withβb2 = 8 and βb2 = 24, we investigated the axial thrust for the rear and front shroud parts in relation to the pressure field. For a low flow rate φ = 0.008, the static pressure distribution in the clearance domain for the rear and front shrouds for the βb2 = 8 case is shown in Figures 20(a)–20(d) and 21(a)–21(d), respectively, whereas the same pressure distributions for the case βb2 = 24 are shown in Figures 22(a)–22(d) and 23(a)–23(d), respectively. Furthermore, the same pressure distributions are shown for a high flow rate φ = 0.047 in Figures 24(a)–24(d) and 25(a)–25(d) for βb2 = 8 and in Figures 26(a)–26(d) and 27(a)–27(d) for βb2 = 24. Here, the static pressure distribution in the clearance domain for the front and rear shrouds is for the cross-sections C-C and D-D, respectively, of Figure 3. For any of the impellers, the CP in the clearance domain of either the front or rear shroud fluctuates greatly during rotation of the impeller in the same way as the static pressure distribution in the interior of the pump, as shown in Figures 14–17. For a low flow rate φ = 0.008, the CP for the clearance domain of the front shroud takes a high value for all the impellers in the region θ0 = 90°–180°. Owing to the effect of leakage flow, the CP for the clearance domain of the front shroud drops overall below that for the clearance domain of the rear shroud but similarly takes a high value in the region θ0 = 90°–180°, which corresponds to the behavior of the axial thrust for the rear and front shroud parts, shown, respectively, in Figures 19(b) and 19(c). At this point, if we examine the static pressure distribution in the interior of the pump, as shown in Figures 14 and 15, the static pressure around the impeller is greater overall in the region θ0 = 90°–180° than in the region θ0 = 180°–270°. That is to say, we can reason that when the static pressure around the impeller increases overall, the static pressure in the clearance domain for the front and rear shrouds increases in accordance with the increase in the leakage flow rate. This is evident from the fact that if we consider the fluctuation in the leakage flow rate coefficient Cq from the liner ring as shown in Figure 28(a) with φ = 0.008, Cq for any of the impellers attains its maximum value in the vicinity of θ0 = 90°–180°, which is larger; the larger value is the βb2 for the impeller.Static pressure distributions at rear shroud (φ = 0.008, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at front shroud (φ = 0.008, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at rear shroud (φ = 0.008, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at front shroud (φ = 0.008, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at rear shroud (φ = 0.047, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at front shroud (φ = 0.047, βb2 = 8).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at rear shroud (φ = 0.047, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Static pressure distributions at front shroud (φ = 0.047, βb2 = 24).
(a)
θ
0 = 0°
(b)
θ
0 = 90°
(c)
θ
0 = 180°
(d)
θ
0 = 270°Leakage flow rate fluctuation (cal.).
(a)
φ = 0.008
(b)
φ = 0.047On the other hand, for a large flow rateφ = 0.047, the CP for the clearance domain of the front and rear shrouds takes a high value at θ0 = 0°(360°), which corresponds to the behavior of the axial thrust for the rear and front shroud parts, as shown in Figures 19(b) and 19(c). At this point, looking at the static pressure distribution in the interior of the pump as shown in Figures 16 and 17, the static pressure around the impeller at θ0 = 0° is higher overall than in the region θ0 = 180°–270°. It therefore follows that, as shown in Figure 28(b), for φ = 0.047, Cq attains its maximum value in the vicinity of θ0 = 0°, which is larger; the larger value is βb2 for the impeller.From the above considerations, because the time-averagedCP at the end of the volute winding for low flow rates is large, in the region θ0 = 90°–180°, and because the time-averaged CP at the beginning of the volute winding for large flow rates is large, at θ0 = 0°, the high CP of the pressure side propagates in the circumferential direction and the CP around the impeller increases. We can assume that accompanying the increase in this CP, the leakage flow rate also increases, and owing to the high pressure acting on the outer surface of the rear and front shrouds, the axial thrust for the rear and front shroud parts therefore increases. Because the larger the blade outlet angle for the impeller, the higher the CP of the pressure side and time-averaged CP at either the beginning or end of the volute winding. We can argue that this tendency is very pronounced.


## 5. Conclusions

From the results of investigating the effect of the blade outlet angle on the fluctuating hydrodynamic forces in a closed-type centrifugal pump with a single blade, using both experiment and numerical analysis, we have found the following.(1)
At low flow rates, all the impellers give a root-mean-square value for the fluctuating component of the total radial thrust that is roughly the same, but for high flow rates, it is larger for impellers with larger blade outlet angles. On the other hand, although the time-averaged value of the total radial thrust also increases with the flow rate for all the impellers, for either low or high flow rates, the difference due to variations in the blade outlet angle is comparatively small.(2)
At low flow rates, the radial thrust for the hydraulic part increases for a blade phase angle in the vicinity ofθ0 = 180° and decreases in the vicinity of θ0 = 0°. These tendencies are more significant in impellers with larger blade outlet angles. On the contrary, for high flow rates, the radial thrust for the hydraulic part increases in the vicinity of θ0 = 0° and is largest in impellers with large blade outlet angles. Conversely, the radial thrust for the hydraulic part decreases in the vicinity of θ0 = 180° and is smallest in impellers with small blade outlet angles.(3)
For the time-averaged value of the axial thrust for the rear and front shroud parts, the absolute value gets larger for impellers with large blade outlet angles at high flow rates. Therefore, for the time-averaged value of the total axial thrust, the absolute value is largest for the impeller with a blade outlet angle of 24°, but hardly any difference is seen between the impellers of 8° and 16°.(4)
The absolute value of the axial thrust for the rear and front shroud parts is largest, for low flow rates, in the vicinity ofθ0 = 90°–180° and for high flow rates, in the vicinity of θ0 = 0°. This is because the leakage flow rate increases with an increase in the static pressure around the impeller, so that a high pressure acts on the outer surface of the rear and front shrouds and there is a greater tendency for this to occur in impellers with large blade outlet angles.



---

*Source: 838627-2014-11-18.xml*