# Component Analysis of Unsteady Hydrodynamic Force of Closed-Type Centrifugal Pump with Single Blades of Different Blade Outlet Angles

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

---

## Abstract

Single-blade centrifugal impellers for sewage systems undergo both unsteady radial and axial thrusts. Therefore, it is extremely important for the improvement of pump reliability to quantitatively grasp these fluctuating hydrodynamic forces and determine the generation mechanism behind them. In this study, we conducted component analyses of radial and axial thrusts of closed, single-blade centrifugal pumps with different blade outlet angles by numerical analysis while considering leakage flow. The results revealed the effect of the blade outlet angle on the components of radial and axial thrusts. For increased flow rates, the time-averaged values of the pressure component were similar for all impellers, although its fluctuating components were higher for impellers with larger blade outlet angles. Moreover, the fluctuating inertia component of the impeller with a blade outlet angle of 8° decreased as the flow rate increased, whereas those with 16° and 24° angles increased. Therefore, the radial thrust on the hydraulic part was significantly higher for impellers with high blade outlet angles.


---

## Body

## 1. Introduction

Sewage pumps for transporting fluids with solid and fibrous foreign bodies are required to have a high foreign body passage diameter; thus, a centrifugal pump with a single wide blade is often used. Currently, foreign body passage is evaluated according to the passed particle size compared to the pump diameter [1]. Nevertheless, even for centrifugal pumps with many blades, geometrical manufacturing errors cause disproportional fluid force [2]. The single-blade centrifugal pumps in current usage are equipped with an axially nonsymmetric volute casing around an axially nonsymmetric impeller. For this reason, the static pressure distribution surrounding the impeller is extremely nonuniform, and an excessively fluctuating radial thrust acts on the impeller [3, 4]. Excessive radial thrust shortens the lifetime of the bearing and shaft seals and lessens the pump’s reliability by causing shaft vibration; thus, it is essential to determine the development mechanism of such vibration.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. For this purpose, we carried out numerical analyses with consideration for leakage flows behind the rear and front shrouds and revealed the behaviors of radial and axial thrusts [14]. We also examined the effect of the blade outlet angel on the radial and axial thrusts [15].As a continuation of a previous report [15], we conducted numerical analyses herein with consideration of leakage flows behind the rear and front shrouds for closed, single-blade centrifugal pumps with different blade outlet angles. We conducted component analyses of radial and axial thrusts and revealed the effect of the blade outlet angle on the thrust components.


## 2. Test Pumps

The three closed, single-blade centrifugal impellers used as test pumps for this study had blade outlet anglesβb2 of 8°, 16°, and 24°, respectively [15]. 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’s rotational speed was 1740 min−1. 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. Moreover, no leakage flow from the pumps at the rear shroud side was noted.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.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 monitoring points of the static pressure 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°.


## 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 [15]. The governing equations include a continuous equation [16] and Reynolds-Averaged Navier-Stokes equations [16], and a SST (Shear Stress Transport) [16] was used for the turbulence model.The computational domain and computational grid for the pump part are shown, respectively, in Figures3(a) and 3(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.Figure 3
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 differences between the two meshes 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 differences between the two meshes 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 [14].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 (n=1740 min−1) 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 flow became nearly periodic.


## 4. Modeling of the Components of Radial and Axial Thrusts

To understand the behavior of unsteady hydrodynamic forces on each part of the impellers, we divided the impeller surface into three parts, as shown in Figure4. We classified the blade surface and the rear and front shroud inner surfaces as the hydraulic part, the rear shroud outer surface, including the outer circumference surface, as the rear shroud part, and the front shroud outer surface, including the front shroud circumference surface, as the front shroud part. The total radial thrust Fr and the total axial thrust Fa applied on the impeller are expressed as follows [14]:(1)Fr=FrH+FrR+FrFFa=FaH+FaR+FaF.Here, FrH is the radial thrust on the hydraulic part [N], FrR is the radial thrust on the rear shroud part [N], FrF is the radial thrust on the front shroud part [N], FaH is the axial thrust on the hydraulic part [N], FaR is the axial thrust on the rear shroud part [N], and FaF is the axial thrust on the front shroud part [N]. For the radial and axial thrusts on the hydraulic part, we modeled the components to examine their behaviors in detail. FrH and FaH are expressed as follows:(2)FrH=FHx2+FHy2FaH=FHz.Figure 4
Impeller surfaces.Let the fluid domain inside the impeller be the test volume, let the impeller suction inlet be the test inlet, and let the impeller outlet be the test outlet. Applying the unsteady momentum conservation law, thex, y, and z components of the hydrodynamic force on the hydraulic part inside the impeller FH can be expressed as follows [14]:(3)FHx=FIx+FMx+FPx+FτxFHy=FIy+FMy+FPy+FτyFHz=FIz+FMz+FPz+Fτz.Here, FI is the inertia component [N], FM is the momentum component [N], FP is the pressure component [N], and Fτ is the shear stress component [N].The inertia componentFI is caused by the time-varying momentum of the fluid inside the impeller. Its x, y, and z components can be obtained as follows:(4)FIx=-ρ∫VdvxdtdVFIy=-ρ∫VdvydtdVFIz=-ρ∫VdvzdtdV.Here, V is the test volume, which is the fluid domain inside the impeller.The momentum componentFM is caused by the difference between the fluid momentum at the impeller outlet and inlet. Its x, y, and z components can be obtained as follows:(5)FMx=-ρ∫A2vx2vr2dA-∫A0vx0vz0dAFMy=-ρ∫A2vy2vr2dA-∫A0vy0vz0dAFMz=-ρ∫A2vz2vr2dA-∫A0vz0vz0dA.Here, A2 is the impeller outlet area [m2] and A0 is the area of the impeller suction inlet [m2].The pressure componentFP is caused by the difference between the fluid momentum at the impeller outlet and inlet. Its x, y, and z components can be obtained as follows:(6)FPx=-∫A2P2-Pscos⁡θdAFPy=-∫A2P2-Pssin⁡θdAFPz=-∫A0P0-PsdA.The shear stress componentFτ is applied on the fluid surface at the impeller outlet and at the suction inlet. Its x, y, and z components can be obtained as follows:(7)Fτx=-∫A2τ2sin⁡γsin⁡θdA+∫A0τ0cos⁡θdAFτy=-∫A2τ2sin⁡γcos⁡θdA+∫A0τ0sin⁡θdAFτz=-∫A2τ2cos⁡γdA.Here, γ is the angle between the z-axis and the direction of the shear stress on the impeller outlet surface.


## 5. Analysis Results and Considerations

### 5.1. Comparisons of Pump Characteristics

Figure5 shows the comparisons of the pump characteristics for βb2=8, βb2=16, and βb2=24. Moreover, it has been confirmed that the calculated values of the head coefficient ψ and the pump efficiency η qualitatively match the experimental results [15]. As the flow rate increased, the head coefficient ψ was higher for impellers with larger βb2. For the flow rate at the maximum efficiency point φ=0.030 for βb2=8, values of ψ for βb2=16 and 24 were larger than that of βb2=8 by 13.7% and 33.8%, respectively. This result is presumably due to the increase in the circumferential component of the absolute velocity at the blade outlet with greater blade outlet angles. For the pump efficiency η, the flow rate with maximum efficiency was φ=0.030 for βb2=8, whereas it moved to a higher rate of φ=0.039 for βb2=16 and 24. At φ=0.039, η was observed to be 0.735 and 0.767 for βb2=16 and 24, respectively, indicating that higher η values occurred for impellers with higher βb2 values. 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.Figure 5
Performance curves.


### 5.2. Radial Thrust and Comparison of Its Components

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, 19]:(8)Fr=F¯r+ΔFr.Here, F¯r is the time-averaged value and ΔFr is the fluctuating component.Figures6(a)–6(d) and 7(a)–7(d) show the Lissajous figures of the total, hydraulic part, rear shroud part, and front shroud part with the radial thrust coefficient. Here, each radial thrust was obtained by integrating the pressure and the shear stress on the wall of each impeller part. 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 9–11. The time-averaged values f¯r of the radial thrust coefficient and the root-mean-square values of the fluctuating component Δfr for the total, hydraulic part, rear shroud part, and front shroud part are shown in Figures 8(a)–8(d). The difference between the calculated and experimental values for the time-averaged values f¯r of the radial thrust coefficient and the root-mean-square values of the fluctuating component Δfr (where βb2=8) was found to be 15.4% and 1.8% at φ=0.008, respectively. Moreover, at φ=0.047, it was −24.4% and −9.7%, respectively. However, we verified that the experimental results qualitatively correlated with the calculated results for the time-averaged value f¯r of radial thrust coefficient and the root-mean-square values of the fluctuating component Δfr [15].Figure 6
Lissajous figures of radial thrust coefficient on each part (βb2=8).
(a)
Total
(b)
Hydraulic part
(c)
Rear shroud part
(d)
Front shroud partFigure 7
Lissajous figures of radial thrust coefficient on each part (βb2=24).
(a)
Total
(b)
Hydraulic part
(c)
Rear shroud part
(d)
Front shroud partFigure 8
Radial thrust coefficient (fr) on each part as functions of flow rate coefficient (φ).
(a)
Total
(b)
Hydraulic part
(c)
Rear shroud part
(d)
Front shroud partFigure 9
Lissajous figures of hydraulic part radial thrust coefficient (φ=0.030).
(a)
β
b
2
=
8
(b)
β
b
2
=
24Figure 10
Lissajous figures of each component of radial thrust coefficient (βb2=8).
(a)
Inertia component
(b)
Momentum component
(c)
Pressure componentFigure 11
Lissajous figures of each component of radial thrust coefficient (βb2=24).
(a)
Inertia component
(b)
Momentum component
(c)
Pressure componentForβb2=8 and 24, the Lissajous figures of the total radial thrust coefficient moved to the +y and -y directions for decreased and increased flow rates, respectively. The Lissajous figure of the radial thrust coefficient for each part showed similar behavior. In particular, the hydraulic part behaved almost identically.For the time-averaged values,f¯r of βb2=8 was located near the origin at each part for φ=0.019–0.030, whereas f¯r of βb2=24 behaved similarly for φ=0.030. At these flow rates, f¯r was minimized at each part and increased as the flow rate deviated from these points. For f¯r of each part, in every impeller, f¯r was highest for the hydraulic part, followed by the rear shroud and front shroud parts. The value of f¯r at each part of the three impellers was approximately the same. In particular, f¯r was almost identical for βb2=8 and 16 in the hydraulic part for the total flow rate domain. This result is due to the pressure component, which will be subsequently discussed.For the root-mean-square values of the fluctuating componentΔfr for every impeller, a minimum of φ=0.019 was set for the total, hydraulic part and rear shroud part. The value increased as the flow rate deviated from the minimum point. For φ=0.008, the root-mean-square values of Δfr were approximately the same for the total, hydraulic part and rear shroud part, whereas the values were larger for higher values of βb2 at increased flow rates. This result is presumably due to the increased load on the blade for increased flow rates for impellers with higher βb2. In the front shroud part, the root-mean-square value of Δfr was highest for βb2=24 at φ=0.047. However, as the flow rate decreased from this point, it fell to below those of βb2=16 and 8 at φ=0.019. This result is believed to have occurred because lower flow rates lead to increased leakage flow rates for impellers with larger βb2 [15]. Moreover, for every impeller, the root-mean-square values of Δfr for the rear and front shroud parts were smaller than that of the hydraulic part.From the aforementioned findings, it can be concluded that the hydraulic part is dominant for each impeller regarding time-averaged values and fluctuating component of the total radial thrust. Next, we examine the components of the radial thrust on the hydraulic part as follows.Figures9(a) and 9(b) show the component analysis values at their maximum efficiency point of φ=0.030 for βb2=8 and 24 and the Lissajous figures of the above-calculated radial thrust coefficient of the hydraulic part. For both impellers, the component analysis values and the Lissajous figures nearly matched. The same approximate match was noted for other flow rates, verifying the validity of the component analysis.Figures10(a)–10(c) and 11(a)–11(c) show the Lissajous figures of the inertia, momentum, and pressure components of the radial thrust of the hydraulic part for βb2=8 and 24. Figures 12(a)–12(c) show the time-averaged f¯r and the fluctuating and root-mean-square values of the inertia, momentum, and pressure components for βb2=8, 16, and 24. We omitted the time-averaged and fluctuating values of the shear stress component because they were considerably small.Figure 12
Each component of radial thrust coefficient (fr) as functions of flow rate coefficient (φ).
(a)
Inertia component
(b)
Momentum component
(c)
Pressure componentFor the time-averaged valuesf¯r, the inertia component of f¯r was near the origin for different flow rates and for each impeller. That is, f¯r of the inertia component was small and constant for different flow rates and for each impeller. Moreover, f¯r of the momentum component in each impeller reached its minimum at φ=0.030 and increased as the flow rate deviated from that point. For points other than φ=0.030, impellers with higher βb2 had higher f¯r. In each impeller, f¯r was located in the direction of the beginning of the volute winding for a low flow rate of φ=0.008 and moved in the -x direction or the direction opposite the outlet, as the flow rate increased. To investigate the behavior of the momentum component of f¯r, we determined the absolute velocity vectors at the impeller outlet for βb2=8 and 24 in Figures 13(a), 13(b), 14(a), and 14(b), respectively. The blade phase angle was θ0=0°. For a low flow rate of φ=0.008, the stagnation point was located toward the outlet rather than the volute tongue for βb2=8 and 24, indicating a strong recycle flow. Therefore, flows from the impeller to the beginning of the volute winding are prevented and flows toward parts other than the beginning of the volute winding become relatively large; moreover, f¯r, in turn, presumably moves toward the beginning of the volute winding. For a high flow rate of φ=0.047, however, the stagnation point for βb2=8 was near the beginning of the volute winding due to the high flow rate. The stagnation point for βb2=24 was located near the volute tongue but on the side of the winding start of the volute. Therefore, flows from the impeller largely proceeded toward the outlet, including the beginning of the volute winding. The movement of f¯r in the -x direction, or opposite the outlet, may have been due to the counteraction to this phenomenon. Similar to the radial thrust in the hydraulic part shown in Figures 6(b) and 7(b), f¯r of the pressure component in every impeller was located near the origin for φ=0.019–0.030 and moved in the +y direction for decreased flow rates and in the -y direction for increased flow rates. That is, f¯r of the pressure component in each impeller achieved its minimum at φ=0.019–0.030 and increased as the flow rate deviated from that point. Therefore, the pressure component of the time-averaged value of the radial thrust in the hydraulic part was dominant for each impeller. Figure 15 shows the circumferential direction distribution of the time-averaged pressure coefficient CP at the monitoring point in the impeller outlet shown in Figure 2. Here, the pressure coefficient CP is the time-averaged value for a single rotation of the runner. At φ=0.030, CP was nearly constant for each circumferential direction in each impeller. The value of CP was high near the winding end of the volute at φ=0.008 and near the beginning of the volute winding at φ=0.047, which corresponds to the behavior of f¯r of the pressure component. The value of f¯r of the pressure component was low for each βb2 and nearly the same likely because CP was higher for impellers with higher βb2 at any flow rate but increased uniformly in each circumferential direction (Figure 15). It has been confirmed that the calculated values of CP at the impeller outlet matched well the experimental values [14].Figure 13
Absolute velocity vectors (βb2=8, θ0=0°).
(a)
φ
=
0.0
08
(b)
φ
=
0.047Figure 14
Absolute velocity vectors (βb2=24, θ0=0°).
(a)
φ
=
0.0
08
(b)
φ
=
0.047Figure 15
Circumferential direction distributions of time-averaged pressure coefficient (CP).For the root-mean-square values of the fluctuating componentΔfr, the Δfr of the inertia component in each impeller with the same blade phase angle θ0 was in a similar direction for different flow rates. The root-mean-square value of Δfr of the inertia component for βb2=8 decreased with increasing flow rates. In contrast, it achieved its minimum at φ=0.019 for βb2=16 and 24 and increased as the flow rate deviated from that point. The value of Δfr was nearly the same for all impellers at values of φ=0.019 and below; for φ=0.019 and above, however, impellers with higher βb2 had higher values. The root-mean-square values of the momentum component were small compared to those of other components; however, they increased for βb2=8 and 16 when the flow rate decreased. Impellers with higher βb2 had higher values, which can likely be attributed to the increased circumferential component vu2 of the absolute velocity of the blade outlet for higher βb2. For every impeller, the root-mean-square value of Δfr of the pressure component for φ=0.019 increased with an increase in flow rate. At φ=0.008, the root-mean-square values Δfr of the pressure and inertia components were similar, although the pressure component was higher than the other components for any other flow rate. Therefore, the pressure component of the fluctuating component of the radial thrust in the hydraulic part was dominant for each impeller. The value of Δfr was nearly the same for any impeller for φ=0.008 and below. For φ=0.019 and above, however, impellers with higher βb2 had higher values, similar to the radial thrust on the hydraulic part. This result is likely due to the larger load on the blade for higher βb2 for the increased flow rates and increased pressure distribution around the impeller.On the basis of these findings, it can be concluded that, for low flow rates, each component has a similar magnitude for every impeller. Moreover, the radial thrust on the hydraulic part is nearly the same magnitude. For increased flow rates, the time-averaged values of the pressure component were similar for all impellers, although the fluctuating component was larger for impellers with larger blade outlet angles. Moreover, the fluctuating inertia component ofβb2=8 decreased as the flow rate increased, whereas those of βb2=16 and 24 increased. Therefore, the radial thrust on the hydraulic part is significantly higher for impellers with high blade outlet angles.


### 5.3. Axial Thrust and Comparisons of Its Components

The time-averaged values of the axial thrust coefficients for the total, hydraulic part, rear shroud part, and front shroud part for each impeller are shown in Figures16(a)–16(d). Here, each axial thrust was obtained by integrating the pressure and the shear stress on the wall of each impeller part.Figure 16
Time-averaged values of axial thrust coefficient (fa) on each part as functions of flow rate coefficient (φ).
(a)
Total
(b)
Hydraulic part
(c)
Rear shroud part
(d)
Front shroud partThe absolute value of the total axial thrust coefficient increased with decreasing flow rates, and the impeller withβb2=24 had the highest values. No significant differences were noted between βb2=8 and 16; βb2=8 had a higher value for low flow rates. For every impeller, the axial thrust coefficient was highest in the -z direction. The axial thrust coefficient on the front shroud part was applied in the +z direction and was slightly smaller than that on rear shroud part. The axial thrust coefficient on the hydraulic part was also applied in the +z direction but was smaller in magnitude. The axial thrust coefficients for the rear and front shroud parts were higher for larger βb2 values; moreover, that for the hydraulic part did not change significantly for different βb2 values. However, unlike the total axial thrust coefficient, that of the hydraulic part increased with an increase in flow rate, comprising 25.5% of the total axial thrust coefficient at φ=0.047 for βb2=8. Subsequently, we examine the components of the axial thrust on the hydraulic part.Figure17 shows the component analysis values of the time-averaged values of the axial thrust for the hydraulic part for βb2=8 and 24. For both impellers, the component analysis values and the calculated values nearly matched, which validates the analysis results.Figure 17
Time-averaged values of hydraulic part axial thrust coefficient (fa) as functions of flow rate coefficient (φ).Figures18(a)–18(c) show the time-averaged values of the inertia, momentum, and pressure components for the hydraulic part for βb2=8, 16, and 24. We omitted the shear stress component because it was quite small.Figure 18
Time-averaged values of each component of axial thrust coefficient (fa) as functions of flow rate coefficient (φ).
(a)
Inertia component
(b)
Momentum component
(c)
Pressure componentA comparison of each component revealed that the inertia component was relatively small forβb2 values other than 16°. For every impeller, the momentum component was in the +z direction. It was the largest for flow rates other than φ=0.008 and increased with increasing flow rates. Flow into the impeller shifted from the axial direction to the radial direction and received a counteracting force in the +z direction. The momentum component was the force on the impeller; the momentum increased with an increase in flow rate, thus increasing the momentum component. For every impeller, the pressure component was in the -z direction, and its absolute value increased with an increase in flow rate. Therefore, for each impeller, the momentum component of the time-averaged values of the axial thrust on the hydraulic part was dominant and was applied in the +z direction. For higher flow rates, however, more than half of the momentum components were counteracted by the pressure components.


## 5.1. Comparisons of Pump Characteristics

Figure5 shows the comparisons of the pump characteristics for βb2=8, βb2=16, and βb2=24. Moreover, it has been confirmed that the calculated values of the head coefficient ψ and the pump efficiency η qualitatively match the experimental results [15]. As the flow rate increased, the head coefficient ψ was higher for impellers with larger βb2. For the flow rate at the maximum efficiency point φ=0.030 for βb2=8, values of ψ for βb2=16 and 24 were larger than that of βb2=8 by 13.7% and 33.8%, respectively. This result is presumably due to the increase in the circumferential component of the absolute velocity at the blade outlet with greater blade outlet angles. For the pump efficiency η, the flow rate with maximum efficiency was φ=0.030 for βb2=8, whereas it moved to a higher rate of φ=0.039 for βb2=16 and 24. At φ=0.039, η was observed to be 0.735 and 0.767 for βb2=16 and 24, respectively, indicating that higher η values occurred for impellers with higher βb2 values. 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.Figure 5
Performance curves.


## 5.2. Radial Thrust and Comparison of Its Components

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, 19]:(8)Fr=F¯r+ΔFr.Here, F¯r is the time-averaged value and ΔFr is the fluctuating component.Figures6(a)–6(d) and 7(a)–7(d) show the Lissajous figures of the total, hydraulic part, rear shroud part, and front shroud part with the radial thrust coefficient. Here, each radial thrust was obtained by integrating the pressure and the shear stress on the wall of each impeller part. 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 9–11. The time-averaged values f¯r of the radial thrust coefficient and the root-mean-square values of the fluctuating component Δfr for the total, hydraulic part, rear shroud part, and front shroud part are shown in Figures 8(a)–8(d). The difference between the calculated and experimental values for the time-averaged values f¯r of the radial thrust coefficient and the root-mean-square values of the fluctuating component Δfr (where βb2=8) was found to be 15.4% and 1.8% at φ=0.008, respectively. Moreover, at φ=0.047, it was −24.4% and −9.7%, respectively. However, we verified that the experimental results qualitatively correlated with the calculated results for the time-averaged value f¯r of radial thrust coefficient and the root-mean-square values of the fluctuating component Δfr [15].Figure 6
Lissajous figures of radial thrust coefficient on each part (βb2=8).
(a)
Total
(b)
Hydraulic part
(c)
Rear shroud part
(d)
Front shroud partFigure 7
Lissajous figures of radial thrust coefficient on each part (βb2=24).
(a)
Total
(b)
Hydraulic part
(c)
Rear shroud part
(d)
Front shroud partFigure 8
Radial thrust coefficient (fr) on each part as functions of flow rate coefficient (φ).
(a)
Total
(b)
Hydraulic part
(c)
Rear shroud part
(d)
Front shroud partFigure 9
Lissajous figures of hydraulic part radial thrust coefficient (φ=0.030).
(a)
β
b
2
=
8
(b)
β
b
2
=
24Figure 10
Lissajous figures of each component of radial thrust coefficient (βb2=8).
(a)
Inertia component
(b)
Momentum component
(c)
Pressure componentFigure 11
Lissajous figures of each component of radial thrust coefficient (βb2=24).
(a)
Inertia component
(b)
Momentum component
(c)
Pressure componentForβb2=8 and 24, the Lissajous figures of the total radial thrust coefficient moved to the +y and -y directions for decreased and increased flow rates, respectively. The Lissajous figure of the radial thrust coefficient for each part showed similar behavior. In particular, the hydraulic part behaved almost identically.For the time-averaged values,f¯r of βb2=8 was located near the origin at each part for φ=0.019–0.030, whereas f¯r of βb2=24 behaved similarly for φ=0.030. At these flow rates, f¯r was minimized at each part and increased as the flow rate deviated from these points. For f¯r of each part, in every impeller, f¯r was highest for the hydraulic part, followed by the rear shroud and front shroud parts. The value of f¯r at each part of the three impellers was approximately the same. In particular, f¯r was almost identical for βb2=8 and 16 in the hydraulic part for the total flow rate domain. This result is due to the pressure component, which will be subsequently discussed.For the root-mean-square values of the fluctuating componentΔfr for every impeller, a minimum of φ=0.019 was set for the total, hydraulic part and rear shroud part. The value increased as the flow rate deviated from the minimum point. For φ=0.008, the root-mean-square values of Δfr were approximately the same for the total, hydraulic part and rear shroud part, whereas the values were larger for higher values of βb2 at increased flow rates. This result is presumably due to the increased load on the blade for increased flow rates for impellers with higher βb2. In the front shroud part, the root-mean-square value of Δfr was highest for βb2=24 at φ=0.047. However, as the flow rate decreased from this point, it fell to below those of βb2=16 and 8 at φ=0.019. This result is believed to have occurred because lower flow rates lead to increased leakage flow rates for impellers with larger βb2 [15]. Moreover, for every impeller, the root-mean-square values of Δfr for the rear and front shroud parts were smaller than that of the hydraulic part.From the aforementioned findings, it can be concluded that the hydraulic part is dominant for each impeller regarding time-averaged values and fluctuating component of the total radial thrust. Next, we examine the components of the radial thrust on the hydraulic part as follows.Figures9(a) and 9(b) show the component analysis values at their maximum efficiency point of φ=0.030 for βb2=8 and 24 and the Lissajous figures of the above-calculated radial thrust coefficient of the hydraulic part. For both impellers, the component analysis values and the Lissajous figures nearly matched. The same approximate match was noted for other flow rates, verifying the validity of the component analysis.Figures10(a)–10(c) and 11(a)–11(c) show the Lissajous figures of the inertia, momentum, and pressure components of the radial thrust of the hydraulic part for βb2=8 and 24. Figures 12(a)–12(c) show the time-averaged f¯r and the fluctuating and root-mean-square values of the inertia, momentum, and pressure components for βb2=8, 16, and 24. We omitted the time-averaged and fluctuating values of the shear stress component because they were considerably small.Figure 12
Each component of radial thrust coefficient (fr) as functions of flow rate coefficient (φ).
(a)
Inertia component
(b)
Momentum component
(c)
Pressure componentFor the time-averaged valuesf¯r, the inertia component of f¯r was near the origin for different flow rates and for each impeller. That is, f¯r of the inertia component was small and constant for different flow rates and for each impeller. Moreover, f¯r of the momentum component in each impeller reached its minimum at φ=0.030 and increased as the flow rate deviated from that point. For points other than φ=0.030, impellers with higher βb2 had higher f¯r. In each impeller, f¯r was located in the direction of the beginning of the volute winding for a low flow rate of φ=0.008 and moved in the -x direction or the direction opposite the outlet, as the flow rate increased. To investigate the behavior of the momentum component of f¯r, we determined the absolute velocity vectors at the impeller outlet for βb2=8 and 24 in Figures 13(a), 13(b), 14(a), and 14(b), respectively. The blade phase angle was θ0=0°. For a low flow rate of φ=0.008, the stagnation point was located toward the outlet rather than the volute tongue for βb2=8 and 24, indicating a strong recycle flow. Therefore, flows from the impeller to the beginning of the volute winding are prevented and flows toward parts other than the beginning of the volute winding become relatively large; moreover, f¯r, in turn, presumably moves toward the beginning of the volute winding. For a high flow rate of φ=0.047, however, the stagnation point for βb2=8 was near the beginning of the volute winding due to the high flow rate. The stagnation point for βb2=24 was located near the volute tongue but on the side of the winding start of the volute. Therefore, flows from the impeller largely proceeded toward the outlet, including the beginning of the volute winding. The movement of f¯r in the -x direction, or opposite the outlet, may have been due to the counteraction to this phenomenon. Similar to the radial thrust in the hydraulic part shown in Figures 6(b) and 7(b), f¯r of the pressure component in every impeller was located near the origin for φ=0.019–0.030 and moved in the +y direction for decreased flow rates and in the -y direction for increased flow rates. That is, f¯r of the pressure component in each impeller achieved its minimum at φ=0.019–0.030 and increased as the flow rate deviated from that point. Therefore, the pressure component of the time-averaged value of the radial thrust in the hydraulic part was dominant for each impeller. Figure 15 shows the circumferential direction distribution of the time-averaged pressure coefficient CP at the monitoring point in the impeller outlet shown in Figure 2. Here, the pressure coefficient CP is the time-averaged value for a single rotation of the runner. At φ=0.030, CP was nearly constant for each circumferential direction in each impeller. The value of CP was high near the winding end of the volute at φ=0.008 and near the beginning of the volute winding at φ=0.047, which corresponds to the behavior of f¯r of the pressure component. The value of f¯r of the pressure component was low for each βb2 and nearly the same likely because CP was higher for impellers with higher βb2 at any flow rate but increased uniformly in each circumferential direction (Figure 15). It has been confirmed that the calculated values of CP at the impeller outlet matched well the experimental values [14].Figure 13
Absolute velocity vectors (βb2=8, θ0=0°).
(a)
φ
=
0.0
08
(b)
φ
=
0.047Figure 14
Absolute velocity vectors (βb2=24, θ0=0°).
(a)
φ
=
0.0
08
(b)
φ
=
0.047Figure 15
Circumferential direction distributions of time-averaged pressure coefficient (CP).For the root-mean-square values of the fluctuating componentΔfr, the Δfr of the inertia component in each impeller with the same blade phase angle θ0 was in a similar direction for different flow rates. The root-mean-square value of Δfr of the inertia component for βb2=8 decreased with increasing flow rates. In contrast, it achieved its minimum at φ=0.019 for βb2=16 and 24 and increased as the flow rate deviated from that point. The value of Δfr was nearly the same for all impellers at values of φ=0.019 and below; for φ=0.019 and above, however, impellers with higher βb2 had higher values. The root-mean-square values of the momentum component were small compared to those of other components; however, they increased for βb2=8 and 16 when the flow rate decreased. Impellers with higher βb2 had higher values, which can likely be attributed to the increased circumferential component vu2 of the absolute velocity of the blade outlet for higher βb2. For every impeller, the root-mean-square value of Δfr of the pressure component for φ=0.019 increased with an increase in flow rate. At φ=0.008, the root-mean-square values Δfr of the pressure and inertia components were similar, although the pressure component was higher than the other components for any other flow rate. Therefore, the pressure component of the fluctuating component of the radial thrust in the hydraulic part was dominant for each impeller. The value of Δfr was nearly the same for any impeller for φ=0.008 and below. For φ=0.019 and above, however, impellers with higher βb2 had higher values, similar to the radial thrust on the hydraulic part. This result is likely due to the larger load on the blade for higher βb2 for the increased flow rates and increased pressure distribution around the impeller.On the basis of these findings, it can be concluded that, for low flow rates, each component has a similar magnitude for every impeller. Moreover, the radial thrust on the hydraulic part is nearly the same magnitude. For increased flow rates, the time-averaged values of the pressure component were similar for all impellers, although the fluctuating component was larger for impellers with larger blade outlet angles. Moreover, the fluctuating inertia component ofβb2=8 decreased as the flow rate increased, whereas those of βb2=16 and 24 increased. Therefore, the radial thrust on the hydraulic part is significantly higher for impellers with high blade outlet angles.


## 5.3. Axial Thrust and Comparisons of Its Components

The time-averaged values of the axial thrust coefficients for the total, hydraulic part, rear shroud part, and front shroud part for each impeller are shown in Figures16(a)–16(d). Here, each axial thrust was obtained by integrating the pressure and the shear stress on the wall of each impeller part.Figure 16
Time-averaged values of axial thrust coefficient (fa) on each part as functions of flow rate coefficient (φ).
(a)
Total
(b)
Hydraulic part
(c)
Rear shroud part
(d)
Front shroud partThe absolute value of the total axial thrust coefficient increased with decreasing flow rates, and the impeller withβb2=24 had the highest values. No significant differences were noted between βb2=8 and 16; βb2=8 had a higher value for low flow rates. For every impeller, the axial thrust coefficient was highest in the -z direction. The axial thrust coefficient on the front shroud part was applied in the +z direction and was slightly smaller than that on rear shroud part. The axial thrust coefficient on the hydraulic part was also applied in the +z direction but was smaller in magnitude. The axial thrust coefficients for the rear and front shroud parts were higher for larger βb2 values; moreover, that for the hydraulic part did not change significantly for different βb2 values. However, unlike the total axial thrust coefficient, that of the hydraulic part increased with an increase in flow rate, comprising 25.5% of the total axial thrust coefficient at φ=0.047 for βb2=8. Subsequently, we examine the components of the axial thrust on the hydraulic part.Figure17 shows the component analysis values of the time-averaged values of the axial thrust for the hydraulic part for βb2=8 and 24. For both impellers, the component analysis values and the calculated values nearly matched, which validates the analysis results.Figure 17
Time-averaged values of hydraulic part axial thrust coefficient (fa) as functions of flow rate coefficient (φ).Figures18(a)–18(c) show the time-averaged values of the inertia, momentum, and pressure components for the hydraulic part for βb2=8, 16, and 24. We omitted the shear stress component because it was quite small.Figure 18
Time-averaged values of each component of axial thrust coefficient (fa) as functions of flow rate coefficient (φ).
(a)
Inertia component
(b)
Momentum component
(c)
Pressure componentA comparison of each component revealed that the inertia component was relatively small forβb2 values other than 16°. For every impeller, the momentum component was in the +z direction. It was the largest for flow rates other than φ=0.008 and increased with increasing flow rates. Flow into the impeller shifted from the axial direction to the radial direction and received a counteracting force in the +z direction. The momentum component was the force on the impeller; the momentum increased with an increase in flow rate, thus increasing the momentum component. For every impeller, the pressure component was in the -z direction, and its absolute value increased with an increase in flow rate. Therefore, for each impeller, the momentum component of the time-averaged values of the axial thrust on the hydraulic part was dominant and was applied in the +z direction. For higher flow rates, however, more than half of the momentum components were counteracted by the pressure components.


## 6. Conclusions

We studied the effects of blade outlet angles of closed, single-blade centrifugal pumps on radial and axial thrusts and their components. The results are presented in the following points:(1)
The pressure component of the time-averaged value and fluctuating component of the radial thrust in the hydraulic part were dominant for each impeller.(2)
For small flow rates each component had similar magnitudes for every impeller, and the radial thrust on the hydraulic part was approximately the same magnitude. For increased flow rates, the time-averaged values of the pressure component were similar for all impellers, although its fluctuating component was higher for impellers with larger blade outlet angles. Moreover, the fluctuating inertia component of the impeller with a blade outlet angle of 8° decreased as the flow rate increased, whereas those with 16° and 24° angles increased. Therefore, the radial thrust on the hydraulic part was significantly higher for impellers with high blade outlet angles.(3)
For each impeller, the momentum component of the time-averaged values of the axial thrust on the hydraulic part was dominant and was applied in the+z direction. For higher flow rates, however, more than half of the momentum components were counteracted by the pressure components.



---

*Source: 419736-2015-08-09.xml*