Eq.

(6) relates the work done by the unit of mass on the wheel with the wheel angular
frequency ω, the outlet and inlet fluid tangential velocities, and the outlet and inlet fluid
radial position. The power produced/absorbed by the wheel is

P = ṁ(−W12 )

If P > 0 the wheel is a turbine; if P < 0, the wheel is a compressor or pump.

2. Application of the Euler Equation

Consider air flowing inside the control volume without swirl (Cθ1 = 0) and leaving the
control volume at the radial position r2 = 15cm with a tangential velocity equal to 90%
of the rotor velocity. The air inlet temperature is 15 C and its pressure is 100kPa. The
rotor frequency is 20,000 revolutions/minute. Determine the compression ratio P out /Pin
assuming that a diffuser is present at the outlet and the velocity at the diffuser exit is
negligible. Also neglect the inlet kinetic energy with respect to the air enthalpy and the
heat transfer in the compressor. Assume that there are no sources of irreversibilities.

Using Euler equation, the work done on the unit of mass of air is W12 = Cθ2 r2 ω. Since
Cθ2 = 0.9ωr2 , then W12 = 0.9(ωr2 )2 . Using ω = 2πf = 2π × 20000/60s−1 , r2 = 0.15m,
we find W12 = 88.736kJ/kg. The energy balance between the inlet and the outlet of the
c.v. is
C12 C2
h1 + + W12 = h2 + 2 (7)
2 2
C22
The energy balance between the inlet and the outlet of the diffuser is h2 + 2
= hout .
Neglecting the kinetic energy at the inlet, Eq. (7) becomes

hin + W12 = hout (8a)

From the properties of perfect gases, dh = Cp dT and

hout − hin = Cp (Tout − Tin ) (8b)

2
Also dh = Cp dT = vdp and pv = RT . Thus,
  CRp
Pout Tout
= (9)
Pin Tin
Determining Tout /Tin from Eq. (8) and substituting into Eq. (9) yields
  Cp
Pout W12 R
= 1+ (10)
Pin Cp Tin
Using Cp = 1kJ/kg · K, T1 = 288K gives

Pout
= 2.55
Pin
The device studied is here is a compressor with compression ratio β = 2.55. It seems that
the only requirement for building such compressor is to design the wheel in such a way
to induce an outlet tangential velocity equal to 90% of the rotor speed and to design the
diffuser to convert kinetic energy into pressure.

How can we do that?????????

3.1 Blades in an Axial Compressors

We start by using the following shape of the control volume interior. In the control
volume there are radial blades as described in Fig. 2. The air flow enters the c.v. in the
axial direction (no swirl), the mass flow rate is ṁ = 0.9 kg/s, the blades radial length is
h = 1cm. The blades are installed on a cylindrical shaft of radius rs = 14.5 cm. The average
radius of the blades is rb = 15cm. The flow area is A1 = π[(rs + h)2 − rs2 ] = 94cm2 and
the air density is ρ1 = p1 /RT1 = 1.2kg/m3 . The average flow velocity of the air entering
the blades is C1 = ṁ/ρ1 A1 = 80m/s and the blade velocity is Cb = ωrb = 314m/s.

~1 − C
In the frame of reference of the rotating blades, the air velocity is C ~b = C
~ r1 .

Defining with α the deflection angle of Cr1 , the blade must be oriented with the same
angle in order to be aligned with the flow. Thus,

C1
α = tan−1 (11)
Cb

3
Since Cb = ωrb = 314m/s, Eq. (11) yields α = 14o .

~ r2 oriented
In the blades frame of reference, the air leaves the blades with a velocity C
~2 = C
like the blades. The air velocity in the rest frame of reference is C ~ r2 + C
~ b . Defining

with φ, the outlet deflection angle of the blade, it follows that

Cθ2 = Cb − Cr2 cos φ (12)

where Cθ2 is the tangential component of the velocity at the exit of the rotor. In order to
design the compressor described in Sec. 2, it is required that Cθ2 = 0.9Cb thus leading to
Cr2 cos φ = 0.1Cb . There are two unknowns in this equation: Cr2 and φ. Now remember
that according to the mass conservation, the axial flow of air is constant. That is

ρ2
ρ1 C1 = ρ2 Cr2 sin φ → C1 = Cr2 sin φ (13)
ρ1

Furthermore, energy conservation in the blade frame of reference requires

2
Cr1 C2
h1 + = h2 + r2 (14)
2 2

Using the perfect gas model (dh = Cp dT , Cp = constant), Eq. (14) can be rewritten as

T2 C 2 − Cr22
= 1 + r1 (15)
T1 2Cp T1

Neglecting heat exchanges and generation of irreversibilities, the pressure and temperature
are related by the following equation
 Cp /R
P1 T1
= (16a)
P2 T2

and the densities (ρ = P/RT ) obey

  CRp −1   CRv
ρ1 T1 T1
= = (16b)
ρ2 T2 T2

Substituting Eq. (16) and (15) into Eq. (13) yields

Cv /R
C 2 − Cr22

C1 = 1 + r1 Cr2 sin φ (17)
2Cp T1

4
 p
Since Cr1 = C12 + Cb2 = 324m/s and Cr2 = 0.1Cb / cos φ, Eq. (17) can be rewritten in
the following form
 Cv /R
(0.1Cb )2
C12 + Cb2 − cos2 φ
C 1 = 1 +  (0.1Cb ) tan φ (18)
2Cp T1

Eq. (18) is the equation for the blade deflection angle φ. Substituting C1 = 80m/s,
Cb = 314m/s, Cp = 1000J/kgK and T1 = 288 K, Cv /R = 2.48, Eq. (18) is numerically
solved and yields φ = 59o . Then using Cr2 = 0.1Cb / cos φ one finds Cr2 = 61m/s.

3.2 The Diffuser

Observe that at the exit of the rotor blades, the air velocity is large,
q
2
C2 = Cr2 − 2Cr2 Cb cos φ + Cb2 = 287m/s (19)

and so is the kinetic energy

C22
Ek = = 41kJ/kg (20)
2
The compression ratio between the inlet and outlet of the rotor is low. From Eq. (15) and
(16a), P2 /P1 = 1.76. In order to transform the kinetic energy into pressure, one can use a
set of diffusers in parallel by appropriately installing a set of static blades (stator blades).

The deflection angle of the air velocity at the exit of the rotor blades can be obtained by
the sine rule C2 / sin φ = Cr2 / sin δ yielding
 
Cr2
δ = sin −1
sin φ = 11o (20)
C2

In order to follow the flow, the stator blades are also deflected by 11o . It is usually desirable
~ 3 = C3 ez ). Defining with A2 and A3
to have axial flow at the exit of the stator blades (C
the flow areas at the inlet and outlet of the stator blades, the conservation of mass gives

C2 A2 ρ3 1
ρ2 A2 C2 sin δ = ρ3 C3 A3 → = (21)
C3 A3 ρ2 sin δ

5
The energy conservation in the stator reads as

C22 C2
h2 + = h3 + 3 (22)
2 2

and using the perfect gas model

Cp /R
C 2 − C32

T3 P3 T3
=1+ 2 = (23)
T2 2Cp T2 P2 T2

Combining Eq. (22) and Eq. (23) yields the following equation
C /R
C22 − C32 v

A2 C3 1
= 1+ (24)
A3 C2 sin δ 2Cp T2

Using Cr1 = 324m/s into Eq. (15) gives T2 = 338 K and Eq. (24) becomes
"  2 #2.5
A2 C3 C3
= 5.5 1.12 − 0.12 (25)
A3 C2 C2

For a given value of C3 /C2 , Eq. (25) yields the stator inlet/outlet ratio. For example:

C3 = 20m/s C3 /C2 = 0.07 A2 /A3 = 0.5

C3 = 10m/s C3 /C2 = 0.035 A2 /A3 = 0.2

For C3 = 20m/s, Eq. (23) yields T3 /T2 = 1.12 and P3 /P2 = 1.48 leading to a compression
ratio
P3
= 1.48 × 1.76 = 2.6
P1

as expected from Sec. 2.

4.1. Axial Flow Turbines

Remember the expression of the work done by the fluid on the Euler wheel (Eq. (6)).
If you can design the wheel in such a way that Cθ2 = 0 and Cθ1 > 0, then the work
−W12 = ωCθ1 r1 > 0 can be extracted from the fluid. Here Cθ is the velocity component
in the tangential direction or direction of rotation of the wheel. Thus, it is important

6
to induce a tangential component to the fluid entering the wheel. This can be achieved

through a set of stator blades. Fig. 5 shows a stator blade deflecting the fluid by 70 o

(α1 = 20o ). Fig. 4 shows the cross section of a turbine. Observe that the stator blades

are connected to the casing and do not move. Their length changes in the axial direction

(See Fig. 6. where R0 < R < R1 ). The Euler wheel of Fig. 1 is represented by the rotor

blades that are connected to the shaft and rotate with its angular frequency. At the rotor

inlet, the fluid has a tangential component Cθ1 . At the outlet of the rotor blades of Fig.

5. the fluid velocity is deflected by the angle α2 ' 90o and Cθ2 = 0.

The stator and rotor represent one stage of the turbine. There are several stages in

each turbine (three stages in Fig. 4). The rotor blades have constant length in the radial

direction (R2 = R1 ).

4.2 Example

Consider a flow of saturated steam at 70 bars and hin = h0 = 2772.1 kJ/kg entering
p
the turbine with velocity C0 = 0.25Cs where Cs = 5p/3ρ is the speed of sound

C0 = 140m/s

The shaft radius is Rs = 18cm and the mass flow rate is ṁ = 1200kg/s. The flow area is

ṁv0
A0 = = 0.23m2 (26)
C0

corresponding to a inlet radius R0 = 32cm. The steam pressure drops from 70 to 10 bars

at constant entropy. The quality at the exit of the turbine is

5.8133 − 2.1387
xout = = 0.826 (27)
6.5863 − 2.1387

The enthalpy at the exit of the turbine is

hout = 0.826 × 2778.1 + (1 − 0.826) × 762.81 = 2427kJ/kg (28)

7
In our design, the enthalpy drop is the same through each stage and the turbine has three

stages (Fig. 4)
hin − hout
h0 − h 2 = = 115kg/kJ (29)
3

The stage reaction is defined as the ratio of the enthalpy drop through the rotor and the

total enthalpy drop through the stage

h1 − h 2
Reaction = (30)
h0 − h 2

For a 50% reaction stage, h0 − h1 = 57.5kJ/kg and h1 − h2 = 57.5kJ/kg.

4.3 Stator Blades

The energy conservation through the stator yields

C02 C2
h0 + = h1 + 1 (31)
2 2

The mass conservation yields

ρ0 C0 A0 = ρ1 C1 A1 sin α1 (32)

From the energy conservation C1 = 367m/s. If we knew the pressure p1 , we could calculate

the density ρ1 to use in the mass conservation equation. The equation for p1 is obtained

by
s1 − sl (p1 ) h1 − hl (p1 )
x1 = = (33)
sv (p1 ) − sl (p1 ) hv (p1 ) − hl (p1 )

where s1 = s0 is the entropy of the mixture, sl , sv , hl , hv are the entropies and enthalpies

of the saturated liquid and vapor at p1 Obviously Eq. (33) has one unknown: the pressure

p1 . I try a few values of p1 and I found that p1 = 52 bars is close to the exact solution.

This value of p1 leads to x1 = 0.95 and

v1 = x1 vv (p1 ) + (1 − x1 )vl (p1 ) = 0.036m3 /kg (34)

8
In order to maximize Cθ1 , we choose a small deflection angle α1 = 20o thus leading to

v1 C0
A1 = A 0 = 1.48 × 0.23 = 0.34m2 (35)
v0 C1 sin α1

corresponding to R1 = 38cm.

4.4 Rotor Blades

The rotor blades are connected to the shaft that rotates at 6000rpm. The average
radius of the rotor blades (Fig. 6) is

Rs + R 1
Rb = = 28cm (36)
2

The blade velocity at the midpoint is Cb = ωRb = 176m/s. In the frame of reference of
the rotor (Fig. 5), the fluid velocity is
q
Cr1 = C12 + Cb2 − 2Cb C1 cos α1 = 210m/s (37)

The angle β1 is
 
C1
β1 = sin −1
sin α1 = 37o (38)
Cr1
In the rotor frame of reference, the conservation of energy reads as
2
Cr1 C2
h1 + = h2 + r2 (39)
2 2

This equation can be solved to find Cr2

q
Cr2 2 = 399m/s
= 2(h1 − h2 ) + Cr1

The angle β2 can be determined by the conservation of the axial flow

ρ2 Cr2 sin β2 A2 = ρ1 Cr1 sin β1 A1 (40)

where A2 = A1 (Fig. 6). To determine ρ2 we need the quality and the pressure at 2

ρ−1
2 = v2 = x2 vv (p2 ) + (1 − x2 )vl (p2 ) (41)

9
For isentropic flows, the pressure p2 can be determined from the quality equation

s2 − sl (p2 ) h2 − hl (p2 )
x2 = = (42)
sv (p2 ) − sl (p2 ) hv (p2 ) − hl (p2 )

where s2 = s1 = s0 = 5.8133 and h2 = h0 − 115 = 2657kJ/kg. I tried a few values of

p2 and p2 = 40 bars seems to be close enough to the exact solution. The corresponding

quality is x2 = 0.92 and the specific volume is v2 = 0.046. Eq. (40) can be solved to find
β2
 
v2 Cr1
β2 = sin −1
sin β1 = 24o (43)
v1 Cr2

The velocity of the fluid in the lab. frame of reference is

q
2 + C 2 − 2C C cos β = 248m/s
C2 = Cr2 (44)
b r2 b 2

and α2 = 40o . Observe that the angle α2 is not 90o as we wished. We could change the
degree of reaction until we find α2 = 90o . However, observe that values of α2 < 90 make
Cθ2 < 0 and therefore they increase the power transfer from the fluid to the rotor, thus

improving the performance of the turbine. The energy variation of the fluid is

C02 C2
∆E = h0 + − h2 − 2 = 94kJ/kg
2 2

This is the work transferred from the fluid to the turbine.

10