Two-Dimensional Performance Analysis
Two-Dimensional Performance Analysis
3.1 Introduction
In this and the following chapter, we briefly survey the more detailed analyses of
the flow in axial and centrifugal pumps, and provide a survey of some of the mod-
els used to synthesize the noncavitating performance of these turbomachines. The
survey begins in this chapter with a summary of some of the results that emerge
from a more detailed analysis of the two-dimensional flow in the meridional plane
of the turbomachine, while neglecting most of the three-dimensional effects. In this
regard, sections 3.2 through 3.4 address the analyses of linear cascades for axial flow
machines, and section 3.5 summarizes the analyses of radial cascades for centrifugal
machines. Three-dimensional effects are addressed in the next chapter.
The fluid mechanics of a linear cascade will now be examined in more detail, so that
the role played by the geometry of the blades and information on the resulting forces
on individual blades may be used to supplement the analysis of section 2.7. Referring
to the periodic control volume indicated in figure 3.1, and applying the momentum
theorem to this control volume, the forces, Fx and Fy , imposed by the fluid on each
blade (per unit depth normal to the sketch), are given by
Fx = −(p2 − p1 )h (3.1)
Fy = ρhvm (w1 cos β1 − w2 cos β2 ) (3.2)
22
3.2 Linear Cascade Analyses 23
2 w2
h vm2
y, Fy
h P2
2
x, Fx
CONTROL
b2
VOLUME
P1
h
1 C D
w1
Fy
L
1
M
b1 – Fx
vm1
Figure 3.1. Schematic of a linear cascade showing the blade geometry, the periodic control volume and
the definition of the lift, L, and drag, D, forces on a blade.
It is conventional and appropriate (as discussed below) to define the lift, L, and the
1
drag, D, components of the total force on a blade, (Fx2 + Fy2 ) 2 , as the components
normal and tangential to the vector mean velocity, wM . More specifically, as shown
in figure 3.1,
where L and D are forces per unit depth normal to the sketch. Nondimensional lift
and drag coefficients are defined as
1 2 1 2
CL = L ρwM c; CD = D ρwM c (3.7)
2 2
The list of fundamental relations is complete if we write the expression for the pressure
difference across the cascade as
ρ 2
p1 − p2 = pLT + w1 − w22 (3.8)
2
24 Two-Dimensional Performance Analysis
where pLT denotes the total pressure loss across the cascade caused by viscous
effects. In frictionless flow, pLT = 0, and the relation 3.8 becomes the Bernoulli
equation in rotating coordinates (equation 2.30 with r1 = r2 as is appropriate here).
A nondimensional loss coefficient, f , is defined as:
1
f = pLT 2
ρwM (3.9)
2
Equations 3.1 through 3.9 can be manipulated to obtain expressions for the lift and
drag coefficients as follows:
CD = 2f sin βM s (3.10)
2 ψ f (φ − cos βM sin βM )
CL = sin βM + (3.11)
s φ sin βM
2 2
ψ is the head coefficient, (p2 − p1 ) ρ R , and
where s = c/h is the solidity, T T
In frictionless flow, when the discharge is parallel with the blades (β2 = βb2 ), this, of
course, reduces to the characteristic equation 2.33. Note that the use of the relation
3.13 allows us to write the expression 3.11 for the lift coefficient as
2
CL = [2 sin βM (cot β1 − cot βM ) − f cos βM ] (3.14)
s
Figure 3.2 presents examples of typical head/flow characteristics resulting from
equation 3.13 for some chosen values of β2 and the friction coefficient, f . It should be
noted that, in any real turbomachine, f will not be constant but will vary substantially
with the flow coefficient, φ, which determines the angle of incidence and other flow
characteristics. More realistic cases are presented a little later in figure 3.3.
3.2 Linear Cascade Analyses 25
1.0
2 = 30°
0
0.05
0.8 0.10
HEAD COEFFICIENT,
0.6
2 = 20°
2 = 10°
0
0.4 f=0
0.05
f = 0.05 0.10
f = 0.10
0.2
1.0
mD = 0
0.8
2 = 20º
HEAD COEFFICIENT,
0.6 mD = 2
0.4
2 = 10º
0.2
1.0
bM = 90º
60º
SHUTOFF HEAD COEFF., O
30º
0.5
10º
The observant reader will have noted that all of the preceding equations of this
section involve only the inclinations of the flow and not of the blades, which have
existed only as ill-defined objects that achieve the turning of the flow. In order to
progress further, it is necessary to obtain a detailed solution of the flow, one result of
which will be the connection between the flow angles (βM , β2 ) and the geometry of the
blades, including the blade angles (βb , βb1 , βb2 ). A large literature exists describing
methods for the solutions of these flows, but such detail is beyond the scope of
this text. As in most high Reynolds number flows, one begins with potential flow
solutions, for which the reader should consult a modern text, such as that by Horlock
(1973), or the valuable review by Roudebush (1965). König (1922) produced one
of the earliest potential flow solutions, namely that for a simple flat plate cascade
of infinitely thin blades. This was used to generate figure 3.4. Such potential flow
methods must be supplemented by viscous analyses of the boundary layers on the
blades and the associated wakes in the discharge flow. Leiblein (1965) provided an
excellent review of these viscous flow methods, and some of his basic methodology
will be introduced later.
To begin with, however, one can obtain some useful insights by employing our
basic knowledge and understanding of lift and drag coefficients obtained from tests,
both those on single blades (airfoils, hydrofoils) and those on cascades of blades. One
such observation is that the lift coefficient, CL , is proportional to the sine of the angle
of attack, where the angle of attack is defined as the angle between the mean flow
direction, βM , and a mean blade angle, βbM . Thus
CL = mL sin(βbM − βM ) (3.15)
3.3 Deviation Angle 27
where, for convenience, the first factor on the right-hand side is denoted by
2mL s sin βbM cot β2 − cot βb2 −1
ψ0 = = 1+ (3.17)
4 + mL s sin βbM cot β1 − cot β2
While the simple, empirical approach of the last section has practical and educational
value, it is also valuable to consider the structure of the flow in more detail, and to
examine how higher level solutions to the flow might be used to predict the perfor-
mance of a cascade of a particular geometry. In doing so, it is important to distinguish
between performance characteristics that are the result of idealized inviscid flow and
those that are caused by viscous effects. Consider, first, the inviscid flow effects.
König (1922) was the first to solve the potential flow through a linear cascade, in
particular for a simple cascade of infinitely thin, straight blades. The solution leads to
values of the deviation, δ, that, in turn, allow evaluation of the shut-off head coeffi-
cient, ψ0 , through equation 3.17. This is shown as a function of solidity in figure 3.4.
Note that for solidities greater than about unity, the idealized, potential flow exits the
blade passages parallel to the blades, and hence ψ0 → 1.
28 Two-Dimensional Performance Analysis
Another approach to the same issue of relating the flow angle, β2 , to the blade
angles, is to employ an empirical rule for the deviation angle, δ = βb2 − β2 (equation
2.2), in terms of other geometric properties of the cascade. One early empirical relation
suggested by Constant (1939) (see Horlock 1973) relates the deviation to the camber
angle, θc , and the solidity, s, through
1
δN = C θc s 2 (3.19)
SUCTION
SURFACE
w1
PRESSURE w2
SURFACE
BLADE
WAKE
Figure 3.5. Sketch of the boundary layers on the surfaces of a cascade and the resulting blade wakes.
3.4 Viscous Effects in Linear Cascades 29
0.06
THICKNESS\CHORD, *\c
WAKE MOMENTUM
0.04
0.02
resulting viscous effects on cascade performance. In this respect, the diffusion factor,
introduced by Lieblein et al. (1953), is a useful concept that is based on the following
approximations. First, we note that under normal operating conditions, the boundary
layer on the suction surface will be much thicker than that on the pressure surface of the
foil, so that, to a first approximation, we may neglect the latter. Then, the thickness
of the wake (and therefore the total pressure loss) will be primarily determined by
that fraction of the suction surface over which the velocity gradient is adverse, since
that is where the majority of the boundary layer growth occurs. Therefore, Lieblein
et al. argued, the momentum thickness of the wake, θ ∗ , should correlate with a
parameter they termed the diffusion factor, given by (wmax − w2 )/wmax , where wmax
is the maximum velocity on the suction surface. One should visualize deceleration or
diffusion of the flow from wmax to w2 , and that this diffusion is the primary factor
in determining the wake thickness. However, since wmax is not easily determined,
Lieblein et al. suggest an approximation to the diffusion factor that is denoted Df ,
and given by
w2 vθ 2 − vθ1
Df = 1 − +
w1 2sw1
sin β1 sin β1 (cot β1 − cot β2 )
= 1− + (3.20)
sin β2 2s
Figure 3.6 shows the correlation of the momentum thickness of the wake (normalized
by the chord) with this diffusion factor, Df , for three foil profiles. Such correlations
are now commonly used to determine the viscous loss due to blade boundary layers
and wakes. Note that, once θ ∗ /c has been determined from such a correlation, the
drag coefficient, CD , and the friction or loss coefficient follow from equations 3.7,
30 Two-Dimensional Performance Analysis
2 sin2 βM θ ∗ s sin βM θ ∗
CD = ; f= (3.21)
sin2 β2 c sin2 β2 c
The data shown in figure 3.6 were for a specific Reynolds number, Re, and the
correlations must, therefore, be supplemented by a statement on the variation of the
loss coefficient with Re. A number of correlations of this type exist (Roudebush
and Lieblein 1965), and exhibit the expected decrease in the loss coefficient with
increasing Re. For more detail on viscous losses in a cascade, the reader should
consult the aforementioned papers by Lieblein.
In an actual turbomachine, there are several additional viscous loss mechanisms
that were not included in the cascade analyses discussed above. Most obviously, there
are additional viscous layers on the inner and outer surfaces that bound the flow, the
hub and the shroud (or casing). These often give rise to complex, three-dimensional
secondary flows that lead to additional viscous losses (Horlock and Lakshminarayana
1973). Moreover, the rotation of other, “non-active” surfaces of the impeller will lead
to viscous shear stresses, and thence to losses known as “disk friction losses” in the
terminology of turbomachines. Also, leakage flows from the discharge back to the
suction, or from one stage back to a preceding stage in a multistage pump, constitute
effective losses that must be included in any realistic evaluation of the losses in an
actual turbomachine (Balje 1981).
The terms SfB and ψ0 result from quite separate and distinct fluid mechanical effects.
The term involving ψ0 is a consequence of the frictionless, potential flow head rise
through any simple, nonrotating cascade whether of axial, radial, or mixed flow
geometry. Therefore, ψ0 is identical to the quantity, ψ0 , defined by equation 3.17 in
3.5 Radial Cascade Analyses 31
2
STREAMLINE OF
FLOW RELATIVE
TO IMPELLER
DISCHARGE
R2
R1 1
INLET b
r
BLADE
Figure 3.7. Schematic of the radial cascade corresponding to the linear cascade of figure 3.1.
the context of a linear cascade. The values for ψ0 for a simple cascade of infinitely
thin blades, whether linear, radial or mixed flow, are given in figure 3.4. The ψ0 term
can be thought of as the “through flow” effect, and, as demonstrated by figure 3.4,
the value of ψ0 rapidly approaches unity when the solidity increases to a value a little
greater than one.
However, it is important to recognize that the ψ0 term is the result of a frictionless,
potential flow solution in which the vorticity is zero. This solution would be directly
applicable to a static or nonrotating radial cascade in which the flow entering the
cacade has no component of the vorticity vector in the axial direction. This would be
the case for a nonswirling axial flow that is deflected to enter a nonrotating, radial
cascade in which the axial velocity is zero. But, relative to a rotating radial cascade
(or centrifugal pump impeller), such an inlet flow does have vorticity, specifically a
vorticity with magnitude 2 and a direction of rotation opposite to the direction of
rotation of the impeller. Consequently, the frictionless flow through the impeller is
not irrotational, but has a constant and uniform vorticity of −2.
In inviscid fluid mechanics, one frequently obtains solutions for these kinds of
rotational flows in the following way. First, one obtains the solution for the irrotational
flow, which is represented by ψ0 in the current problem. Mathematically, this is the
complementary solution. Then one adds to this a particular solution that satisfies all the
same boundary conditions, but has a uniform vorticity, −2. In the present context,
this particular, or rotational, solution leads to the term, SfB , which, therefore, has a
32 Two-Dimensional Performance Analysis
BLADE
DISCHAGE
CIRCLE
d/2
d/2
DISCHARGE INLET
SURFACE SURFACE BLADES
Figure 3.8. A sketch of the displacement component of the inviscid flow through a rotating radial
cascade.
quite different origin from the irrotational term, ψ0 . The division into the rotational
solution and the irrotational solution is such that all the net volumetric flow through the
impeller is included in the irrotational (or ψ0 ) component. The rotational solution has
no through flow, but simply consists of a rotation of the fluid within each blade passage,
as sketched in figure 3.8. Busemann (1928) called this the displacement flow; other
authors refer to its rotating cells as relative eddies (Balje 1980, Dixon 1978). In his
pioneering work on the fluid mechanics of turbomachines, Stodola (1927) was among
the first to recognize the importance of this rotational component of the solution.
Busemann (1928) first calculated its effect upon the head/flow characteristic for the
case of infinitely thin, logarithmic spiral blades, in other words the simple cascade in
the radial configuration. For reasons which will become clear shortly, the function,
SfB , is known as the Busemann slip factor, and Busemann’s solutions lead to the
values presented in figure 3.9 when the solidity, s > 1.1. Note that the values of SfB
are invariably less than or equal to unity, and, therefore, the effect of the displacement
flow is to cause a decrease in the head. This deficiency can, however, be minimized
by using a large number of blades. As the number of blades gets larger, SfB tends to
unity as the rotational flow within an individual blade passage increasingly weakens.
In practice, however, the frictional losses will increase with the number of blades.
Consequently, there is an important compromise that must be made in choosing the
number of blades. As figure 3.9 shows, this compromise will depend on the blade angle.
Furthermore, the compromise must also take into account the structural requirements
for the blades. Thus, radial machines for use with liquids usually have a smaller number
of blades than those used for gases. The reason for this is that a liquid turbomachine
3.5 Radial Cascade Analyses 33
1.0
32
16
8
SLIP FACTOR, SfB
6
4
0.5
2
ZR = 1
Figure 3.9. The Busemann slip factor, SfB , plotted against the blade angle, βb , for various numbers of
blades, ZR . The results shown are for radial cascades of infinitely thin logarithmic spiral blades with
solidities, s > 1.1. Adapted by Sabersky, Acosta and Hauptmann (1989) and Wislicenus (1947) from
Busemann’s (1928) theory.
requires much thicker blades, and, therefore, each blade creates much more flow
blockage than in the case of a gas turbomachine. Consequently, liquid machines tend
to have a smaller number of blades, typically eight for the range of specific speeds for
which radial machines are designed (ND < 1.5) (Stepanoff 1948, Anderson). Another
popular engineering criterion (Stepanoff 1948) is that ZR should be one third of the
discharge blade angle, βb (in degrees).
The decrease in the head induced by the displacement flow is due to the nonuni-
formity in the discharge flow; this nonuniformity results in a mean angle of discharge
(denoted by β2 ) that is different from the discharge blade angle, βb2 , and, there-
fore, implies an effective deviation angle or slip, Sf (see section 2.1). In fact, it is
clear that the relations 2.16, 2.32, 3.22, and 2.4 imply that Sf = SfB , and, hence,
the terminology used above. Stodola (1927) recognized that slip would be a conse-
quence of the displacement flow, and estimated the magnitude of the slip velocity,
vθs , in the following approximate way. He argued that the slip velocity could be
roughly estimated as d/2, where d/2 is the radius of the blade discharge circle
shown in figure 3.8. He visualized this as representative of the rotating cell of fluid
in a blade passage, and that the rotation of this cell at would lead to the afore-
mentioned vθs . Then, provide ZR is not too small, d ≈ 2πR2 sin βb2 , and it follows
that
vθ s = πR2 sin βb2 /ZR (3.23)
and, from equation 2.4, that the estimated slip factor, SfS , is
π sin βb2
SfS = 1 − (3.24)
ZR
34 Two-Dimensional Performance Analysis
T
JE
E
AK
W
PRESSURE
SURFACE
SUCTION
SURFACE
At design flow, the wake or boundary layer on the suction surface may be quite
thin, but as the flow coefficient, φ, is decreased, the increased incidence leads to larger
wakes (Fischer and Thoma 1932, Johnston and Dean 1966). Clearly, the nonuniformity
of the discharge flow implies an “effective” slip due to these viscous effects. This slip
will not only depend on the geometry of the blades but will also be a function of
the flow coefficient and the Reynolds number. The change with flow coefficient is
particularly interesting. As φ is decreased below the design value and the wake grows
in width, an increasing fraction of the flow is concentrated in the jet. Johnston and
Dean (1966) showed that this results in a flow that more closely follows the geometry
of the pressure surface, and, therefore, to a decrease in the slip. This can be a major
effect in radial compressors. Johnston and Dean made measurements in an 18-bladed
radial compressor impeller with a 90◦ discharge blade angle (for which SfS = 0.825),
and found that the effective slip factor increased monotonically from a value of about
0.8 at φ2 = 0.5 to a value of 1.0 at φ2 = 0.15. However, this increase in the slip factor
did not produce an increase in the head rise, because the increase in the viscous losses
was greater than the potential gain from the decrease in the slip.
Finally, it is important to recognize that secondary flows can also have a substantial
effect on the development of the blade wakes, and, therefore, on the jet-wake structure.
Moreover, the geometric differences between the typical radial compressor and the typ-
ical centrifugal pump can lead to significant differences in the secondary flows, the loss
mechanisms, and the jet-wake structure. The typical centrifugal pump geometry was
illustrated in figure 2.7, to which we should append the typical number of blades,
ZR = 8. A typical example is the geometry at ND = 0.6, namely RT 1 /RT 2 ≈ 0.5
and B2 ≈ 0.2RT 2 . Assuming ZR = 8 and a typical blade angle at discharge of 25◦ , it
36 Two-Dimensional Performance Analysis
follows that the blade passage flow at discharge has cross-sectional dimensions nor-
mal to the relative velocity vector of 0.2RT 2 × 0.3RT 2 , while the length of the blade
passage is approximately 1.2RT 2 . Thus the blade passage is fairly wide relative to its
length. In contrast, the typical radial compressor has a much smaller value of B2 /RT 2 ,
and a much larger number of blades. As a result, not only is the blade passage much
narrower relative to its length, but also the typical cross-section of the discharge flow
is far from square, being significantly narrower in the axial direction. The viscous
boundary layers on the suction and pressure surfaces of the blades, and on the hub and
shroud (or casing), will have a greater effect the smaller the cross-sectional dimen-
sions of the blade passage are relative to its length. Moreover, the secondary flows
that occur in the corners of this passage amplify these viscous effects. Consequently,
the flow that discharges from a blade passage of a typical radial compressor is more
radically altered by these viscous effects than the flow discharging from a typical
centrifugal pump.