Slip Factor Prediction for Impellers
Slip Factor Prediction for Impellers
Abstract: Slip factor accounts for the deviation of the flow Φ flow coefficient = 𝑉𝑉𝑟𝑟 /𝑈𝑈𝑒𝑒
angle from the trailing edge blade angle at the exit of radial σ slip factor = 1 − ∆𝑊𝑊/𝑈𝑈𝑒𝑒
impellers. Accurate values are required to predict impeller Ω rotor angular velocity
torque and energy input. The slip factor prediction method
for back-swept radial bladed impellers is based on the so- Subscripts
called single relative eddy (SRE) method, which is an 1-5 version of F in SRE method
approximation of the classical, two-dimensional analytical c critical (beyond which slip factor depends on RR)
solution to the inviscid flow problem. The relatively cor corrected for radius ratio beyond the critical
accurate prediction of the slip factors of 19 impellers found e exit
in four data sets published since its formulation, shows the eq equivalent
reliability of the SRE method. The characteristics of i inlet
straight-bladed impeller layouts are explored and r radial
incorporated into the SRE method. It turns out that the ϑ circumferential
SRE method as developed for logarithmic spiral blades Superscript
predicts the inviscid flow slip factor of 42 different straight ˊ ideal (non-slip) value
bladed impeller geometries accurately, when a correction
for the critical radius ratio at low blade numbers is
1 Introduction
introduced. Predicted slip factor values are also compared
Designers of centrifugal impellers need methods to
to new experimental data for five different impellers with
accurately predict the impeller drive shaft torque and power
straight, back-swept blades, over a range of three blade
input into the flow. This is especially important in the design
angles and three blade numbers. Agreement is excellent
of simple centrifugal fans where the cost of the electric motor
near the flow coefficient corresponding to the volute design
may exceed the cost of the fan itself. No method developed
angle.
specifically for back-swept, straight-bladed centrifugal
impellers could be found in the literature. The question is
Additional keywords: Slip factor, Relative eddy, Rotor then whether classical, inviscid, two dimensional analytical
blades, Straight blades, Back-swept blades, SRE methods for logarithmic spiral blades apply to these blades,
and whether generally used, semi-empirical slip factor
Nomenclature prediction methods are adequate.
Lewis [1], Visser et al. [2] and Hassenpflug [3]
Roman
recalculated the classical, two-dimensional analytical
c blade chord length
solutions to the inviscid flow problem solved by Busemann
BH Busemann [4] data from Hassenpflug [3]
[4], Theodorsen [5, 6], and others. Visser et al. [2]
F solidity influence coefficient in SRE method
summarized the previous work, and applied conformal
F effective area coefficient in method of Qiu [16]
mapping to find exact solutions in terms of Fourier
f length (figure 1)
coefficients given by the Gauss hypergeometric series and the
RMS root mean squared (equation 8)
Beta function, for flow through two-dimensional radial
RR radius ratio; 𝑟𝑟𝑖𝑖 /𝑟𝑟𝑒𝑒
impellers with logarithmic spiral blades. These authors did
RRCF radius ratio correction factor
consider unswept, straight (radial) blades insofar as they
r radial position from axis
represent a limiting case of the logarithmic spiral blade.
SRE single relative eddy
Hassenpflug [3], however, was the only one, after
s blade spacing (pitch)
Theodorsen [5], to extend his approach from logarithmic
TH Theodorsen [5,6] data from Hassenpflug [3]
spiral blade shapes to straight, back-swept blades.
TSF theoretical slip factor, from Memardezfouli [21]
The rigorous, two-dimensional, inviscid, analytical
U rotor tangential velocity
solutions of Busemann [4], Visser et al. [2] and Hassenpflug
V absolute velocity
[3] depend on the number of blades, Z, the blade exit angle,
x length (figure 1)
βe and the impeller radius ratio, RR. RR is the ratio of the
Z number of blades
radial location of the rotor blade leading edge to that of the
Greek trailing edge. These methods apply the Kutta condition (local
β blade angle (from radial direction) flow parallel to the blade at the blade exit) and a similar, so
called shock-free condition, at the blade inlet. To ensure zero-
a. Department of Mechanical and Mechatronic incidence flow at the rotor blade leading edges at all through-
Engineering, Stellenbosch University, South Africa. flow conditions, the use of adjustable inlet guide vanes would
E-mail: [email protected] be implied in practice. To avoid this complication, many
R & D Journal of the South African Institution of Mechanical Engineering 2019, 35, 55-63 55
http://dx.doi.org/10.17159/2309-8988/2019/v35a7
http://www.saimeche.org.za (open access) © SAIMechE All rights reserved.
Slip Factor Prediction for Impellers with Straight, Back-swept Blades
practical radial turbomachines have impellers with inducers the single relative eddy (SRE) method of Von Backström [12,
that guide the flow into the impellers, but sometimes the high 13, 14], while some recent methods, e.g. Ji et al. [15] and Qiu
solidity near the rotor entrance fulfills this function et al. [16] add the effects of blade curvature.
automatically. Like all the semi-empirical slip factor
prediction methods mentioned in this paper, the SRE neglects 3 The SRE method
the flow conditions at the impeller inlet. Von Backström [12] showed that five of the common
approximate methods (Stodola [10], Stanitz [11], Wiesner
2 Slip factor [9], Eck [17] and Csanady [18] are actually simplified
Slip factor accounts for the flow deviation with respect to the versions of the SRE method. It emerged during its derivation
blade trailing edge angle at the impeller exit. In the that the eddy-induced slip factor, equation 4 above), depends
fundamental case of two-dimensional, inviscid flow, slip on the blade solidity, c/se at the impeller exit multiplied by a
factor is determined by the vorticity of the internal flow universal coefficient, F, which is a function of the impeller
relative to the rotating impeller, resulting in a so-called blade exit angle and the number of blades:
relative eddy or eddies. These present as a reduction in the 1
absolute circumferential velocity component at the rotor exit σ = 1 − ∆𝑊𝑊𝑒𝑒 ⁄Ω𝑟𝑟𝑒𝑒 = 1 − (5)
1 + 𝐹𝐹 × 𝑐𝑐/𝑠𝑠𝑒𝑒
compared to the ideal case where the exit relative flow
follows the local blade direction. Engineering and teaching The coefficient, F, is the ratio of the sum of the average
texts like Aungier [7], Lewis [1], Whitfield and Baines [8] eddy-induced velocities on the suction and pressure surfaces,
typically non-dimensionalise the slip velocity through to that along the impeller exit rim. This coefficient depends
division by the ideal absolute circumferential velocity on βe and Z, and has been developed further in successive
component at the blade trailing edges: papers by Von Backström [12, 13, 14] to improve the
′
agreement between the SRE method and the numerical values
σ = 𝑉𝑉θ𝑒𝑒 ⁄𝑉𝑉𝛳𝛳𝛳𝛳 (1) from the two-dimensional analytical solution of Hassenpflug
The ideal circumferential component is determined from [3], henceforth called the BH data. The best empirical
the ideal, no-slip rotor exit velocity triangles as: approximation was [13]:
′
𝑉𝑉𝛳𝛳𝛳𝛳 = 𝑈𝑈𝑒𝑒 − 𝑉𝑉𝑟𝑟𝑟𝑟′ 𝑡𝑡𝑡𝑡𝑡𝑡𝛽𝛽𝑒𝑒′ (2) 3
𝐹𝐹 = 𝐹𝐹4 = 2 + �2.7 + � � 𝑐𝑐𝑐𝑐𝑐𝑐𝛽𝛽𝑒𝑒 � 𝑐𝑐𝑐𝑐𝑐𝑐𝛽𝛽𝑒𝑒 (6)
Assuming no pre-rotation at the rotor inlet, the non- 𝑍𝑍
dimensionalised slip velocity then equals the fractional As explained by Von Backström [14], the digit 2, outside
reduction in rotor torque and in power demand due to slip. the brackets in equation 6 ensures that F4 = 2 when 𝛽𝛽𝑒𝑒 =
The slip factor equals the ratio of torque or power, with slip 90° , resulting in excellent agreement with the BH data at
included, to torque or power, without slip. high values of βe. A remark of Hassenpflug [3] that, for an
In this paper, however, the slip factor is written in terms impeller with 2 straight radial blades (𝑐𝑐𝑐𝑐𝑐𝑐𝛽𝛽𝑒𝑒 = 1.0), the
of the eddy-induced slip velocity at the rotor blade exit, given inviscid flow slip factor should be exactly 1/2, led to further
by: investigation. For RR = 0.5, and the conditions mentioned
′
∆𝑊𝑊𝑒𝑒 = ∆𝑉𝑉𝛳𝛳𝛳𝛳 = 𝑉𝑉𝛳𝛳𝛳𝛳 − 𝑉𝑉𝛳𝛳𝛳𝛳 (3) above, c/se = 1/ (2π), and so F must be equal to 2π to make F
× c/se equal to the integer, 1. This is achieved by replacing
The slip velocity is then normalised through division by
the (3/ Z) in equation 6 by the more meaningful (π /Z), the
the circumferential speed, Ue of the rotor blade trailing edges,
half angle between adjacent blades, and by replacing the 2.7
to give the slip factor as:
by (1.5 π-2). These substitutions allow us to rewrite F in
σ = 1 − ∆𝑊𝑊𝑒𝑒 / 𝑈𝑈𝑒𝑒 = 1 − ∆𝑊𝑊𝑒𝑒 ⁄Ω𝑟𝑟𝑒𝑒 (4) terms of π as:
This definition is chosen because, unlike equation 2, it is π
𝐹𝐹 = 𝐹𝐹5 = 2 + �(1.5π − 2) + � � 𝑐𝑐𝑐𝑐𝑐𝑐𝛽𝛽𝑒𝑒 � 𝑐𝑐𝑐𝑐𝑐𝑐𝛽𝛽𝑒𝑒 (7)
independent of the through-flow, as it does not contain the 𝑍𝑍
radial velocity component. This makes it useful for This empirical correction ensures perfect agreement with
comparing rotors that do not operate at exactly the same flow the BH data for 2 straight radial blades (βe = 0), and an error
coefficients, as does Wiesner [9]. Also, Whitfield and Baines of 0.3 % for Z > 2. When Z = 1, the error is < 2.5 %. The
[8] state that the definition in equation 4 is more commonly comparative curve fit of Paeng and Chung [19] requires five
adopted than equation 1, and support the notion that it is the equations for similar accuracy, and is not defined when Z =
slip velocity, ∆We that is important. Visser et al. [2] state 1.
that for inviscid methods the use of equation 4 results in lower The current paper applies the SRE method to impellers
sensitivity to through-flow. Such solutions disregard effects with straight, back-swept blades, by making the initial
such as boundary layer growth and separation as well as wake assumption that the slip factor depends primarily on the blade
formation that contribute to flow deviation at blade trailing solidity, as defined by the sum of the blade lengths divided
edges. Since these effects are dependent on flow rate, they by the impeller circumference. The initial objective is then
may contribute to some extent to slip in practice. to write the solidity for straight blades in terms of the number
Simple, approximate slip factor prediction methods like of blades, their exit angle and the radius ratio for straight,
those of Stodola [10] and Stanitz [11] consider the slip factor back-swept blades. As a first approximation the solidity
to be dependent on the number of blades and blade trailing influence coefficient, F5, as developed for logarithmic spiral
edge angle only, but to be independent of blade shape. Some blades is applied, and the resulting prediction compared to the
other methods include corrections for the effect of the experimental data.
impeller inner to outer radius ratio (RR), e.g. Wiesner [9] and
R & D Journal of the South African Institution of Mechanical Engineering 2019, 35, 55-63 56
http://dx.doi.org/10.17159/2309-8988/2019/v35a7
http://www.saimeche.org.za (open access) © SAIMechE All rights reserved.
Slip Factor Prediction for Impellers with Straight, Back-swept Blades
The validity of the solidity as a primary variable that Table 1 Comparison of slip factor prediction methods:
determines slip, has however been questioned by Huang et al. Huang [20]
[20]. This really relates to the effect of radius ratio, RR, as it Case 1 2 3 4 5 % RMS
is the only relevant variable that appears in the solidity but difference
not in its influence coefficient, F. This concern had been Z 13 15 15 15 17 -
addressed in the SRE method by assuming RR = 0.5 𝛽𝛽𝑒𝑒 (°) 30.0 30.0 41.8 47.0 41.8 -
Huang et 0.868 0.866 0.883 0.895 0.895 0.0
whenever RR < 0 [12]. It implies that slip factor does not al. data
depend on solidity when RR < 0.5. It also assumes that the Wiesner 0.845 0.860 0.870 0.876 0.881 1.6
drop-off in slip factor at larger radius ratios starts at RR = 0.5 SRE-F5 0.844 0.861 0.869 0.874 0.882 1.7
and that the drop-off rate is correctly predicted by the basic
SRE method. This assumption needs reconsideration, (RMS) value of the difference, to well below 2 %, as does
however, as discussed later. To check the validity of the SRE Wiesner [9] (table 1).
method, its performance will be compared to that of other Memardezfouli and Nourbakhsh [21] compare various
more recent data sets and prediction methods. slip factor prediction methods to their measured data for five
pump impellers. They compare the methods of Wiesner [9],
4 Objectives Stodola [10] and Stanitz [11], as well as an early version of
The general objective of this study is to investigate the slip the SRE [12] and their own theoretical slip factor (TSF),
factor of impellers with straight, back-swept blades, and to which is a superposition of an axi-symmetrical meridional
develop a method to predict their performance. Specific plane analysis and a blade-to-blade singularity method.
objectives of the work are to: Unfortunately they did not, in their application of the SRE
1. Investigate the geometric peculiarities of straight- method, correct the blade reference angle from radial to
bladed impellers tangential (i.e. change cosines in the SRE equations to sines).
2. Derive the SRE slip factor equations for the blade The SRE equation they list under their table also has 1 – RR2
solidity of straight bladed impellers. in it instead of 1 – RR. Table 2 shows the correct version of
3. Compare predicted and published slip factors. the SRE method, and compares, in the last column, the RMS
4. Select, design and manufacture a family of straight- differences. At their relatively large exit blade angles (60 to
bladed radial impellers. 72 °.) the SRE method agrees slightly better with their data
5. Measure the slip factor of the family of straight- than their own TSF method (2.1 % RMS difference compared
bladed back-swept impellers to 2.4), with both showing deviations of about half of that of
6. Compare predicted and measured slip factors Wiesner [9].
7. Make recommendations regarding the use of the Table 2 Comparison of slip factor prediction methods:
SRE and other methods for straight, back-swept Memardezfouli [21]
blades.
Case 1 2 3 4 5 % RMS
difference
5 Review of recent literature Z 6 8 6 6 4 -
Slip factors predicted by the SRE method will be compared 𝛽𝛽𝑒𝑒 (°) 66.0 66.0 60.0 72.0 66.0 -
Experiment 0.77 0.80 0.74 0.82 0.72 0.0
in terms of RMS differences, to those from four data sets
TSF 0.79 0.83 0.78 0.82 0.71 2.4
published since the SRE papers. The RMS value of the Wiesner 0.82 0.85 0.80 0.84 0.76 4.6
difference between xi and yi for i = 1 to N is defined as: SRE-F5 0.789 0.832 0.769 0.817 0.717 2.1
1 𝑁𝑁
𝑅𝑅𝑅𝑅𝑅𝑅 = � � (𝑥𝑥𝑖𝑖 − 𝑦𝑦𝑖𝑖 )2 (8) Ji et al. [15] developed a slip factor correlation for radial
𝑁𝑁 1 and mixed-flow impellers, based on the Euler turbo-
machinery equation. It incorporates the dependence of slip
The reference row in each of the tables below is the one
factor on the flow coefficient to better represent experimental
showing the zero RMS difference. The objective is to gauge
data. The relative eddy effect is introduced through a Stodola
how well the SRE predicts slip factor compared to these data
[10] type approach, where a circular eddy is fitted into each
sets. Due to space constraints the prediction methods from
blade passage near its exit. The Stodola model is inaccurate
literature are not discussed in detail.
Huang et al. [20] used commercial centrifugal compressor for low blade numbers (Z < 16 and high blade angle (βe > 50°)
software to design six impellers with 13, 15 and 17 blades (Von Backström, [14]). To this they then add the contribution
with exit angles of 30, 41.8 and 47 °, and modelled their the of blade curvature. Table 3 shows that the SRE with F5
performance at mass flows from 80 to 120 % of the design does about as well as their method (2.7 % RMS difference
flow, at zero tip clearance, using CFD. They investigated the compared to 2.8) with the Wiesner RMS difference even
effect of blade curvature, referring to a change in blade angle lower at 2.0 %.
relative to the local radial direction along the blade length. Qiu et al. [16] presented the development and validation
They found that it could affect the variation of slip factor with of a slip factor model for design and off-design, similar to that
flow rate, although, in their examples, the presented of Ji et al. [15]. They incorporate the flow coefficient, φe at
maximum variation in slip factor over their flow range is the rotor exit to cater for off-design variation in slip factor,
about 2 %. For the examples we could reproduce, the SRE and present four test cases that can be used as benchmarks.
predicted the slip factors, in terms of the root mean squared The Eckardt [22] O and A rotors, as well as a pump and
blower test case, all have relatively high sensitivity to flow
R & D Journal of the South African Institution of Mechanical Engineering 2019, 35, 55-63 57
http://dx.doi.org/10.17159/2309-8988/2019/v35a7
http://www.saimeche.org.za (open access) © SAIMechE All rights reserved.
Slip Factor Prediction for Impellers with Straight, Back-swept Blades
1
𝜎𝜎 = 1 − (17)
𝐹𝐹𝐹𝐹
1+ [𝑐𝑐𝑐𝑐𝑐𝑐𝛽𝛽𝑒𝑒 − 𝑅𝑅𝑅𝑅 𝑐𝑐𝑐𝑐𝑐𝑐𝛽𝛽𝑖𝑖 ]
2𝜋𝜋
The intention is to use the coefficient, F = F5, equation 7
listed above.
leaving the impeller. The velocity components were then coefficient of 0.14, the measured slip factors are expected to
corrected to account for the difference between the radial be the most reliable near that value. The SRE method for the
position of the measuring station and the exit radius of the
impeller, by considering the flow areas and conservation of
angular momentum. It was found that the average mass flow
balance between the bell mouth and the exit probe
measurement was within 0.1 %, but the standard deviation of
the difference was a relatively high 8.8 %. Since the
circumferential velocity component was about four times as
large as the radial component, a small error in measured angle
has a large effect on the radial component, and virtually no
effect on the circumferential component.
R & D Journal of the South African Institution of Mechanical Engineering 2019, 35, 55-63 61
http://dx.doi.org/10.17159/2309-8988/2019/v35a7
http://www.saimeche.org.za (open access) © SAIMechE All rights reserved.
Slip Factor Prediction for Impellers with Straight, Back-swept Blades
eddy-induced slip factor is not sensitive to flow coefficient, Table 8 Comparison of slip factor prediction methods:
but the Qiu et al. [16] method predicts only a small effect for Present
these cases, and the measured data show no clear trends for Case 1 2 3 4 5 % RMS
our data. The equations used to calculate the Qiu et al. [16] difference
slip factors are given in the appendix of their paper. For Z 4 8 8 8 16 -
application here, their equations were simplified to remove 𝛽𝛽𝑒𝑒 (°) 22.52 0.0 22.52 28.02 22.52 -
RR 0.50 0.50 0.50 0.50 0.50 -
the effects of disc coning, blade thickness and fluid density
Experimental 0.617 0.734 0.734 0.741 0.854 0.0
ratio. φ 0.169 0.186 0.186 0.178 0.190 -
𝛽𝛽𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒𝑒 (°) 73 76 71 73 78 -
F (Qiu et al.) 0.396 0.854 0.750 0.739 0.900 -
𝑑𝑑𝑑𝑑⁄𝑑𝑑𝑑𝑑 (m-1) 5.308 0.0 5.308 6.862 5.308 -
Qiu et al. 0.727 0.665 0.752 0.774 0.851 4.7
Wiesner 0.591 0.767 0.776 0.781 0.862 2.4
SRE-F5 0.641 0.765 0.769 0.772 0.866 2.1
RRc 0.281 - - - - -
SRE-F5Cor 0.597 0.765 0.769 0.772 0.866 2.0
10 Conclusions
The equations for the SRE method applied to impellers with
straight, back-swept blades have been derived. It was found
that in straight-bladed impellers the possible exit blade angles
are constrained by the impeller radius ratio. This means that
when radius ratios are moderate to low, say below 0.5, the
exit blade angle is then limited to a maximum of 30 °.
The previously published universal solidity influence
Figure 12 Measured and predicted slip factors versus exit coefficient in the SRE method has been slightly adjusted to
flow coefficient: 𝑍𝑍 = 8, 𝛽𝛽𝑒𝑒 = 28.02° achieve better numerical agreement with exact analytical data
for two-dimensional inviscid flow through logarithmic spiral
blades. Comparison with methods and measurements
published since the previous SRE papers show that, in general
the SRE method predicts slip factors with a typical RMS
difference of about 2 %. This establishes the use of an
accurate, inviscid relative eddy model as a basis for a
relatively accurate practical slip factor prediction method,
and questions in principle the use of the inaccurate Stodola
relative eddy model as basis, before adding addition terms to
account for flow variations and blade curvature.
A new radius ratio correction for cases with low numbers
of blades ensures excellent agreement to the analytical
prediction of slip factor for impellers with straight, back-
swept blades. The slip factor, as defined and used in these
cases is independent of flow coefficient. A set of impellers
with straight, back-swept blades was designed, manufactured
Figure 13 Measured and predicted slip factors versus exit and tested. It was found that their slip factors vary very little
flow coefficient: 𝑍𝑍 = 16, 𝛽𝛽𝑒𝑒 = 22.52° when the flow coefficient is close to the volute design flow
coefficient. The SRE method predicts these slip factors to
Slip factors and absolute flow angles at rotor exit for each
within 2 % RMS.
fan are presented at the same four throttle settings. The flow
In general it is recommended that standard slip factor
angles agree more or less with the volute angle of 72 °. For
prediction methods be used for impellers with straight, back-
each impeller the middle two flow coefficients resulted in
swept blades. For the range for radius ratio, blade number and
more or less the same slip factor, with the one at the higher
blade exit angle of straight-bladed impellers, it turns out that
flow corresponding better with the design flow coefficient.
accounting for the critical radius ratio is more important than
For these cases (table 8) the RMS difference between
to model the sensitivity to flow coefficient.
experiment and prediction was 2.0 % for the SRE compared
to 4.7 % for Qiu et al. [16] and 2.4 % for Wiesner [9]. Since
the method of Qiu et al. is blind to the effect of radius ratio, Acknowledgement
it over predicts the slip factor for the low blade number of 4 Stephen Ryall designed and installed the volute, and did the
in figure 9, where the critical radius ratio is less than 0.5. experimental work.
Also, since radial blades (βe =βi = 0) have no effective
curvature, unlike back-swept, ”straight” blades, the
curvature term is zero, and so it under predicts in figure
10.
R & D Journal of the South African Institution of Mechanical Engineering 2019, 35, 55-63 62
http://dx.doi.org/10.17159/2309-8988/2019/v35a7
http://www.saimeche.org.za (open access) © SAIMechE All rights reserved.
Slip Factor Prediction for Impellers with Straight, Back-swept Blades
R & D Journal of the South African Institution of Mechanical Engineering 2019, 35, 55-63 63
http://dx.doi.org/10.17159/2309-8988/2019/v35a7
http://www.saimeche.org.za (open access) © SAIMechE All rights reserved.