22-Mec-B6 Advanced Fluid Mechanics · December 2017
Question 5 of 8: Compressor First Stage
Nivaar worked solution (AI-drafted; not reviewed by a licensed engineer)
Notes on this paper
Paper format. 16-Mec-B6 Fluid Machinery,
National Examinations December 2017 — three hours, closed book. Section A is
calculative (Questions 1–5) and Section B descriptive (Questions 6–8); the rubric asks for
four of Section A and two of Section B, six questions of ten marks each for a sixty-mark paper.
Reference data for individual questions are supplied as Attachments (pages 10–17) and a general
nomenclature/constants/equations sheet occupies pages 18–22. All eight questions are
solved here, because the set is a study resource rather than a timed attempt.
Reference texts.
S. L. Dixon & C. A. Hall, Fluid Mechanics and Thermodynamics of Turbomachinery, 7th ed.;
R. K. Turton, Principles of Turbomachinery, 2nd ed.;
H. Cohen, G. F. C. Rogers & H. I. H. Saravanamuttoo, Gas Turbine Theory, 6th ed.;
F. M. White, Fluid Mechanics, 8th ed.;
R. W. Fox, A. T. McDonald & P. J. Pritchard, Introduction to Fluid Mechanics, 9th ed.;
Y. A. Çengel & M. A. Boles, Thermodynamics: An Engineering Approach, 9th ed.
Constants are those printed on page 19 of the paper
(g = 9.81 m/s², ρwater = 1000 kg/m³,
ρair = 1.21 kg/m³ at 15 °C and 1.19 kg/m³
at 20 °C, cp = 1.005 kJ/kg°C, cv = 0.718 kJ/kg°C,
patm = 100 kPa, pvapour = 1.71 kPa at 15 °C).
Every reference equation quoted below is one of those printed on pages 20–22, and is identified
as such where it is first used.
Subject note. Page 1 of the examination reads
16-MEC-B6 FLUID MACHINERY, and every question is a turbomachine question. The
solutions below answer the paper as printed. The three Canadian hydro stations named in Questions 1
and 2 (Mactaquac, Churchill Falls) are used by the examiner as settings for stated
“hypothetical measurements”; the numbers solved here are the paper's, not plant
records.
Given. First (N1) rotor stage of an industrial gas-turbine compressor:
Quantity
Symbol
Value
Hub diameter at inlet
Dh
480 mm
Tip diameter at inlet
Dt
1120 mm
Guide-vane exit (absolute) angle
α1
30° from axial
Rotor blade exit (relative) angle
β2
40° from axial
Rotational speed
N
6800 rev/min
Air mass flow
M
136 kg/s
Inlet temperature
T1
15°C = 288 K
Inlet air density (page 19)
ρ
1.21 kg/m³
Specific heats (page 19)
cp, cv
1.005, 0.718 kJ/kg°C
The page-17 attachment defines both angles from the axial direction, and both velocity
triangles share the blade speed at the mean diameter. The axial velocity is taken constant through
the rotor, which is the standard assumption for a first stage of constant annulus area and the only
one that makes part (b) solvable.
Find. Mean blade speed and axial velocity, the four velocities on the two
triangles, the specific work and stage power, the enthalpy and temperature rise, and the isentropic
pressure ratio of the stage.
Approach. Mean diameter and continuity give U and
Cx; the two given angles then close the inlet and outlet triangles; page 21's
compressor work w = U(CY2 −
CY1) gives the work, which for an adiabatic rotor is the stagnation
enthalpy rise, hence the temperature rise; the isentropic relation converts that into a pressure
ratio.
Rotor inlet and outlet velocity triangles at the mean diameter, to a common scale of 10 m/s = 4 mm. The rotor turns the relative flow from 55.6° to 40°; the swirl it adds to the absolute flow is the work.
Part (a) — mean blade speed. The mean diameter of the annulus is the
arithmetic mean of hub and tip:
$$D_m=\frac{0.480+1.120}{2}=0.800\ \text{m}$$
and the blade speed there is
$$U=\frac{\pi D_mN}{60}=\frac{\pi(0.800)(6800)}{60}=\boxed{284.8\ \text{m/s}}$$
Part (a) — axial velocity from continuity. The annulus flow area is the
difference of two circles:
$$A=\frac{\pi}{4}\left(D_t^{2}-D_h^{2}\right)=\frac{\pi}{4}\left(1.120^{2}-0.480^{2}\right)
=\frac{\pi}{4}(1.0240)=0.8042\ \text{m}^2$$
and continuity M = ρACx1 (page 20) gives
$$C_{x1}=\frac{M}{\rho A}=\frac{136}{1.21\times 0.8042}=\boxed{139.8\ \text{m/s}}$$
The page-19 density of 1.21 kg/m³ is consistent with
p/RT = 100 000/(287 × 288) = 1.210 kg/m³,
so the two data sources agree to three figures.
Part (b) — close the rotor inlet triangle. The guide vanes leave the air
at 30° to the axial, so the absolute velocity and its swirl component are
$$C_1=\frac{C_{x1}}{\cos\alpha_1}=\frac{139.8}{\cos 30^{\circ}}=\boxed{161.4\ \text{m/s}},\qquad
C_{Y1}=C_{x1}\tan\alpha_1=139.8\tan 30^{\circ}=80.7\ \text{m/s}$$
Subtracting that swirl from the blade speed gives the relative velocity seen by the rotor:
$$W_{Y1}=U-C_{Y1}=284.8-80.7=204.2\ \text{m/s},\qquad
W_1=\sqrt{139.8^{2}+204.2^{2}}=\boxed{247.4\ \text{m/s}}$$
at a relative angle of
β1 = tan−1(204.2/139.8) = 55.6°.
Part (b) — close the rotor outlet triangle. Holding the axial velocity
constant and applying the given relative exit angle of 40°,
$$W_2=\frac{C_{x1}}{\cos\beta_2}=\frac{139.8}{\cos 40^{\circ}}=\boxed{182.4\ \text{m/s}},\qquad
W_{Y2}=C_{x1}\tan\beta_2=139.8\tan 40^{\circ}=117.3\ \text{m/s}$$
Adding the blade speed back gives the absolute state leaving the rotor:
$$C_{Y2}=U-W_{Y2}=284.8-117.3=167.6\ \text{m/s},\qquad
C_2=\sqrt{139.8^{2}+167.6^{2}}=\boxed{218.2\ \text{m/s}}$$
leaving at α2 = 50.2° to the axial. Note the two things
that matter physically: the relative velocity has fallen from 247.4 to 182.4 m/s, which is
the diffusion that raises the pressure, while the absolute velocity has risen from 161.4 to
218.2 m/s, which is the swirl the stator must later remove.
Part (c) — specific work and stage power. Page 21 gives the compressor
work as
$$w=U(C_{Y2}-C_{Y1})=284.8(167.6-80.7)=284.8\times 86.9=24\,750\ \text{J/kg}
=\boxed{24.75\ \text{kJ/kg}}$$
so the power the shaft must supply to this one stage is
$$P=wM=24.75\times 136=\boxed{3366\ \text{kW}}$$
The energy form of the same equation on page 21,
w = [(C2² − C1²) +
(W1² − W2²)]/2, returns the same
24.75 kJ/kg, which checks all four velocities at once.
Part (d) — enthalpy and temperature rise. The rotor is adiabatic and does
work on the air, so the whole of that work appears as a rise in stagnation enthalpy:
$$\Delta h_0=w=\boxed{24.75\ \text{kJ/kg}}$$
and with cp = 1.005 kJ/kg°C from page 19,
$$\Delta T=\frac{\Delta h_0}{c_p}=\frac{24.75}{1.005}=\boxed{24.6\ \text{K}}$$
so the stage leaves at T2 = 288 + 24.6 = 312.6 K. About
25 K per stage is exactly the loading a subsonic axial stage is designed for — push much past
30 K and the diffusion becomes too severe for the boundary layer.
Part (e) — isentropic stage pressure ratio. The ratio of specific heats
follows from the paper's own constants,
k = cp/cv = 1.005/0.718 =
1.400, and page 21's isentropic relation inverted gives
$$\frac{p_2}{p_1}=\left(\frac{T_2}{T_1}\right)^{k/(k-1)}
=\left(\frac{312.6}{288.0}\right)^{3.50}=(1.0855)^{3.50}=\boxed{1.333}$$
A first-stage ratio of 1.33 against the machine's overall gas-generator pressure ratio of 14.1
(page 16) implies roughly nine equivalent stages spread across the two spools, which is consistent
with the cross-section drawn in the attachment.
Quantity
Symbol
Result
(a) Mean diameter
Dm
0.800 m
(a) Mean blade speed
U
284.8 m/s
(a) Annulus area
A
0.8042 m²
(a) Axial velocity
Cx1
139.8 m/s
(b) Absolute inlet velocity
C1
161.4 m/s at 30°
(b) Relative inlet velocity
W1
247.4 m/s at β1 = 55.6°
(b) Absolute outlet velocity
C2
218.2 m/s at α2 = 50.2°
(b) Relative outlet velocity
W2
182.4 m/s at 40°
(c) Specific work
w
24.75 kJ/kg
(c) Stage power
P
3366 kW
(d) Enthalpy rise
Δh0
24.75 kJ/kg
(d) Temperature rise
ΔT
24.6 K (to 312.6 K)
(e) Stage pressure ratio
p2/p1
1.333
Check: constant axial velocity, and angles from the axial.
The paper gives no outlet annulus dimensions, so the outlet triangle can only be closed by taking
Cx2 = Cx1 — the usual
first-stage design assumption, and the one the question's own wording (“ideal conditions”)
invites. Second, both angles are measured from the axial direction, as the page-17 attachment
defines them. Reading them from the tangential instead swaps every sine for a cosine and yields
CY1 = 242 m/s and a negative stage work, which is
physically impossible for a compressor — that sign is the quickest way to catch the error.