M M P P: Chapter 5 Solutions Problem 5.1
M M P P: Chapter 5 Solutions Problem 5.1
Problem 5.1
pt 0 pt 2 0.1q 0
2
q0 p0 M 0
2
pt 0
p0
(1 ( 1) M 0 / 2) /( 1)
2
Therefore,
2
p M0
d t 2 1 0.1 ≈0.9685
2 (1 ( 1) M 0 / 2) /( 1)
2
pt 0
( 1) /
( d ) r 1
d 0.9278
r 1
Problem 5.2
From the flight Mach number, we calculate τr and πr. The flight total pressure and
temperature are calculated using isentropic formulas. With pt2 specified, we calculate πd.
Then we calculate ηd. The area ratio is determined on continuity equation, i.e.,
( 1) / 2 ( 1)
A2 pt 0 M 0 1 ( 1) M 2 / 2
2
A0 pt 2 M 2 1 ( 1) M 0 / 2
2
From Mach number and total pressure at 2, we calculate p2. Also, we calculate the non-
dimensional entropy rise following the Gibbs equation:
s p
ln t 2
R pt 0
2 pt 0
From C p ,stag 2
1 we calculate the Cp,stag=1.170
M 0 p 0
From isentropic table for M0=0.80, we get A0/A*=1.0382 and since we know the external
diffusion area ratio, A0/A1=0.92, we calculate:
We use continuity to get the diffuser area ratio, A2/Ath, similar to problem 5.2. We get:
A2/Ath=1.317
From M1≈0.66 and pt1=pt0=15.24 kPa, we calculate p1≈11.38 kPa. Therefore the non-
dimensional additive drag becomes:
Dadd
0.0061
p 0 A1
Problem 5.4
Also since πd=1.348/1.498=0.8998, the non-dimensional entropy rise is minus the natural
log of πd, therefore,
Δs/R≈0.1056
Problem 5.6
2
pt 0.1q1 0.1 p1 M 1 p t1 pt 2
2
2
pt 2 M1
1 0.1
Therefore, /( 1) ≈0.974
p t1 2 1 ( 1) M 1 2 / 2
From continuity equation, we have all the RHS terms of the following equation:
( 1) / 2 ( 1)
M2 A1 pt1 1
M 1 0.4228
1 ( 1) M 2 / 2
2
( 1) / 2 ( 1)
A2 p t2 1 ( 1) M 1
2
/ 2
≈1.17157
p1 pt1 1 ( 1) M 2 2 / 2
Therefore CpR≈0.473
Problem 5.7
From M0=0.85 and isentropic table, we get A0/A*=1.0207 and pt0/p0=1.6038, which gives
pt0=40.1 kPa and since A1=3 m2, A0=2.7 m2
Therefore, A1/A*=1.134 which from isentropic table gives M1≈0.65
From isentropic table, we get pt1/p1=1.328
Since pt1=pt0, we get p1=30.18 kPa
1
1 2 2 ( 1) 1 2 1 2 1
1 M0 1 M0 1 M0
Dadd M
M 1 2 2 M0 2 1 =0.012
p0 A1 1 1 2 1
1 2 1 2
M1 1 M1 1 M1
2 2 2
46
Therefore Dadd=0.01246(25 kPa)(3 m2)≈935 N
pAM
m AV
a
The speed of sound, a0=315.59 m/s and thus mass flow rate is
m 0 254.5 kg / s
Ram drag is Dr m 0V0 68.28 kN
From continuity, we can calculate A2 since in the following equation; the only unknown
is the area, A2.
m
pt 2
R Tt 2
A2 M 2 1 ( 1) M 2 / 2
2
( 1) / 2 ( 1)
Problem 5.8
Problem 5.9
For M1=0.3, we get p1=25.01 kPa from isentropic tables, using the same pt0 that still
prevails at the inlet lip.
1
1 2 2 ( 1) 1 2 1 2 1
1 M0 1 M0 1 M0
Dadd M
M 1 2 2 M0 2 1 ≈0.08
p0 A1 1 1 2 1
1 2 1 2
M1 1 M1 1 M1
2 2 2
1
Therefore, Dadd≈5376 N
Problem 5.10
From the flight Mach number and isentropic table, we get A1/A*=1.483 which from
isentropic table gives a corresponding Mach number at the lip of M1=0.436.
Since A1/Ath=1.15, then Ath/A*=1.2895, which gives Mth≈0.53. Below critical throat
Mach number of 0.75, therefore, there is no adverse shock boundary layer interaction due
to the convex throat geometry.
Problem 5.11
Problem 5.12
rt = 2 r h
rsp = rh+0.25(rt-rh) = 0.75 rh+0.25 rt = 1.25 rh this is the radius of the spoiled sector.
pt p spoiled
V 2
/2 spoiled Vspoiled
2
0.85 2 0.7225
pt p uniform V 2
/ 2 uniform
V
uniform
2 2
(V A V2 A2 )
pt ,area avg p 1 1
2 A1 A2
The area A1 is the outer flow area where the speed of the gas is V1 and A2 is the spoiled
area of the compressor annulus.
In terms of percentages of the entire compressor face area, A1=81.25% Aface and
A2=18.75% Aface, Also A2=0.2308 A1
2 2
(V A V2 A2 ) 2 2
pt ,area avg p 1 1 p (0.8125 V1 0.1875 V2 )
2 A1 A2 2
2
pt ,area avg p (0.948 V1 )
2
2 2
(V m V2 m 2 )
pt ,mass avg p 1 1
2 m 1 m 2
We replace the mass flow rate by the product of density area and velocity to get:
2 1 (V2 / V1 ) ( A2 / A1 )
3 3
(V1 A1 V2 A2 )
3
pt ,mass avg p p V1
2 V1 A1 V2 A2 2 1 (V2 / V1 )( A2 / A1 )
Again, substitute for A2=0.2308 A1 and V2=0.85 V1 to simplify the above expression to:
2
pt , mass avg p (0.9545V1 )
2
Problem 5.13
Problem 5.14
We apply continuity to the control volume to get the side mass flow rate (equation 5.12),
0 0V0 ( A A0 ) m
m s 0V0 ( A AM )
1 m
s 0V0 ( AM A0 )
m
Applying the momentum equation to the control volume in the x-direction gives:
0 (V1 V0 ) Fx ,cowl ( p1 p 0 ) A1
m
Fx ,cowl m
0 (V1 V0 ) ( p1 p 0 ) A1 = Dcowl
The cowl drag coefficient is based on AM, therefore we construct the following:
Dcowl AM A V 2 p1
2 0 1 1 1
2
q 0 A1 AM A1 V0 M 0 p 0
AM A V 2 p1
C D ,cowl 2 0 1 1 2
1
A1 A1 V0 M 0 p0
A A1 Aspil V1 2 p1
C D ,cowl 2 1 1 1
2
A
M A1 V
0 M 0 0
p
A Aspil V1 2 p1
C D ,cowl 2 1 1 1 1
2
A
M A 1 0 V M 0 0
p
Note that the first bracket on the RHS is the nacelle geometry, the second bracket is the
contribution of supersonic spillage, the third bracket produces a negative quantity related
to the flow deceleration, V1/V0<1. The pressure drag contribution is given by the last
term on the RHS.
Problem 5.15
Problem 5.16
From M0=2.0 and β=45o using the O.S. chart, we get θ=14.74o and
Mn0=1.414 which from a NS table gives Mn1=0.733 and pt1/pt0=0.954. Also
M1=Mn1/sin(β-θ)=1.456
From NS tables pt2/pt1=0.9430 and My=0.7171
The shock system total pressure recovery is:
pt2/pt0=(0.954)(0.943)=0.900
Also, we may use isentropic table between the compressor face and station “y”, which is
downstream of the terminal shock. Based on two Mach numbers, My and M2 we
calculate
A2/Ath=1.2378
Problem 5.17
We may develop a simple spreadsheet (for a NS inlet in super-critical mode) that
produces the following output:
A0/A1
M0 (spec) M-y pt-y/pt-x (NS) (spec) A-y/A* A1/A*
1.6 0.668437365 0.89520026 0.9 1.119211 1.243567
We write a simple Excel spreadsheet program that calculates the following variables:
M-Des A1/A-th (geo) M0 sub. (choked) M0>M0 sub A0/Ath sub % spillage (sub)
1.6 1.250235 0.56 0.7 1.0943727 12.46664199
Problem 5.19
M-Des A1/A-th (geo) M0 sub. (choked) M0>M0 sub A0/Ath sub % spillage (sub)
1.75 1.386492048 0.48 0.7 1.0943727 21.06895383
Problem 5.20
M-
Design A-1/A* A-1/A-th M-y (NS) A-y/A* A-1/A-th' % opening
4 10.71875 10.71875 0.43495884 1.4872932 1.487293243 620.688408
3.5 6.789620536 6.789620536 0.45115388 1.4458332 1.445833224 369.59915
3 4.234567901 4.234567901 0.47519096 1.3903945 1.390394489 204.558737
Problem 5.21
M- %
Design A-1/A* A-1/A-th M-y (NS) A-y/A* A-1/A-th' opening
3.5 6.789620536 6.789620536 0.45115388 1.4458332 1.445833224 369.59915
A-th'/A* M-th'
4.695991502 3.1
Problem 5.22
pt-
M-design M-y (NS) A1/Ath (design) M-th (started) M-y (NS-th) recovery(best)
2 0.577350269 1.216474641 1.55 0.6841011 0.913188272
Problem 5.23
M-design M-y (NS) A1/Ath (design) M-th (started) M-y (NS-th) pt-recovery(best)
1.7 0.640543841 1.144617989 1.45 0.7195619 0.944836685
Problem 5.24
This is a K-D inlet designed for M0=3.2. We again use the spreadsheet that we designed
for the K-D inlets to get:
M-design M-y (NS) A1/Ath (design) M-th (started) M-y (NS-th) pt-recovery(best)
3.2 0.464348678 1.414554828 1.75 0.6280887 0.834565392
Problem 5.25
Problem 5.26
M-
Design A-1/A* A-1/A-th M-y (NS) A-y/A* A-1/A-th'
1.6 1.250235 1.250235 0.668437365 1.1192107 1.119210698
Problem 5.27
M-Des A1/A-th (geo) M0 sub. (choked) M0>M0 sub A0/Ath sub % spillage (sub)
1.5 1.176167052 0.62 0.7 1.0943727 6.954316032
A/A*-
M-y (NS)des Ay/A* % spill.(sup) My-overspeed oversp M-th-oversp
0.701088742 1.093584262 7.02134877 1.784524596 1.422344 1.75
Problem 5.28
pt-
M-design M-y (NS) A1/Ath (design) M-th (started) M-y (NS-th) recovery(best)
2 0.577350269 1.216474641 1.55 0.6841011 0.913188272
Problem 5.29
Problem 5.30
M-
Design A-1/A* A-1/A-th M-y (NS) A-y/A* A-1/A-th'
2.6 2.895975385 2.895975385 0.503871026 1.3325049 1.332504859
Problem 5.31
Problem 5.32
This is a K-D inlet with a design Mach number of 3.0. We use our spreadsheet to
calculate the following variables:
A-shock/Ath
(spec) Ath/A* A-shock/A* M-x shock M-y-shock pt-y/pt-x
1.25 1.386492048 1.73311506 2 0.5773503 0.72087386
ptx=pt0=10 kPa(1+0.2(9))3.5=367.3 kPa
pty=0.7209(367.3 kPa)=264.8 kPa
py=pty/(1+0.2My2)3.5=211.26 kPa
A-shock/A1
M-0 A1/A* (spec) A-shock/A* M-x shock M-y shock pt-recov. Shock
1.86 1.5068943 1.2 1.808273161 2.05 0.569062515 0.697508439
Problem 5.34
A spreadsheet program is written that marches through two oblique shocks (double-ramp
case). The results are tabulated below:
Beta-1
M0 (deg) Theta 1 (deg) Mn-0 pt-1/pt-0 Mn-1 M1
2.5 30.01 8.005285065 1.250377855 0.987006 0.812424258 2.168298255
Theta 2 (over
M1 Beta-2 theta-1) Mn-1 pt2/pt1
2.168298 38.3995 12.00828019 1.346818818 0.9704077
Problem 5.35
We use the spreadsheet program for the analysis of double-ramp inlet in a reverse
iteration mode to do a quick ramp design angle study to achieve the conditions specified
in the problem. Namely, we wish to have the Mach number downstream of the terminal
NS to be 0.7 and the first two oblique shocks be of the same strength. This requires the
normal component of the Mach number into the two OS be the same. The results of the
spreadsheet iteration are summarized below:
Beta-1
M0 (deg) Theta 1 (deg) Mn-0 pt-1/pt-0 Mn-1 M1
2.5 34 12.19723956 1.397982257 0.9586961 0.74056051 1.993902668
Theta 2 (over
M1 Beta-2 theta-1) Mn-1 pt2/pt1 Mn-2 M2
1.993903 44 13.86815412 1.385081177 0.9618333 0.746077084 1.486234426
Problem 5.36
A9/A-th M9 Neu-10
(spec) (isen) p9/pt-9 Neu-9 (deg) Plume angle (deg)-given (deg)
2 2.15 0.1011317 30.42527507 15 45.4252751
In the over-expanded mode where shocks appear in the exhaust flow, we develop a
simple spreadsheet program that calculates the NPR and the plume contraction angle, etc.
In this problem, a NS appears at the exit. The results are:
A9/A-th (spec) M9 (isen) p9/pt-9 NPR (perfect) NPR (NS at the lip)
2 2.15 0.1011317 9.888096679 1.892006062
Problem 5.37
This is an over-expanded nozzle. We use our spreadsheet to get the basics of this nozzle:
A9/A-th (spec) M9 (isen) p9/pt-9 NPR (perfect) NPR (NS at the lip)
2.4 2.35 0.0739569 13.52138738 2.154373612
We have throat Mach number, pt,Tt and area, therefore we substitute in continuity
36.96 kg/s
equation to get m
Problem 5.38
b) This is the under-expanded case. With the specified 10o plume angle, we use the
spreadsheet for under-expanded nozzles to get:
c) We may use the O.S. chart (or look-up array in a spreadsheet) to get the shock angle at
the lip that turns the flow 4o inward when the Mach number is 2.75. We get β~23.8o.
Therefore, M9sinβ≈1.111, which gives a static pressure ratio py/px=1.2728.
NPR=(pt/p9)/(py/px)=19.75
The ambient static pressure is 1900 kPa/19.75=96.2 kPa which gives a delta,
δ=96.2/101=0.952 From altitude table (Appendix A) we get ~ 1400 ft.
Problem 5.39
From isentropic relation between pt and p, we calculate the Mach number, then A/A*
follows from isentropic relations. Finally static temperatures, speed of sound and exit
velocity are calculated. The spreadsheet results are:
pt (kPa) Tt (K) p-0 (kPa) NPR M9
101 2000 1 101 3.7001393
Problem 5.40
We use the over-expanded spreadsheet program and obtain the following results:
Problem 5.41
δ=50.5/101=0.50 which from altitude tables gives 18 kft as the altitude where the nozzle
is perfectly expanded.
For p0=1/2(50.5 kPa)=25.25 kPa, the Mach number downstream of the expansion waves
is M10 which from isentropic pt/p=38.13(101 kPa)/25.25 kPa=152.52. the corresponding
Mach number is
M10=4.0
The Prandtl-Meyer angle for M10=4 is ν10=65.79o
The Prandtl-Meyer angle for M9=3.5 is ν9=58.53o
Problem 5.42
The nozzle area ratio is A9/A*=3.5. From isentropic table, we get M9≈2.8
NPRcrit=(1/πn)[(γ+1)/2]γ/(γ-1) ≈1.87
% Thrust
Eta-n V9-s (m/s) M9 (perfect) V9 perfect F-g-conv/m-dot gain
0.9699239 564.3301 1.6044184 812.2086 792.9710877 2.426007
CD8=πcn=0.98
Convergent
Nozzle
Problem 5.44
CA(conical)=(1+cosα)/2=0.953
Problem 5.45
T-hot/T-
cold T-cold/T-cold Alpha T-mixed/T-cold % Fg-improvement
3 1 3.5 1.444444444 1.033691589
Problem 5.46
d) From My=0.447 we get Tty/Ty=1.0387 and thus Ty=2500 K/1.0387=2406 K and ay=983
m/s and thus Vy=V9=440 m/s
e) For As/A8=4.23 we use isentropic table to get Ms=3.0. From NS tables, we get
py/px=10.33, pty/ptx=0.328 and My=0.4752. From My=0.4752, we get Ay/A*=1.38
Tt9/T9=1.011, therefore T9=2500 K/1.011=2472 K and a9=996.5 m/s and V9=240 m/s
f) Based on My=0.4752, Tty/Ty=1.046 and thus Ty=2500 K/1.046=2390 K and ay=980 m/s
Therefore, Vy=466 m/s and since A9e=As V9=Vy=466 m/s
g) Gross thrust nearly doubled with flow separation in the divergent portion of the nozzle.
The flow downstream of the NS in the nozzle is a diffusing flow, therefore velocity drops
and static pressure is gained. With separation velocity remains nearly constant and thus
nozzle exit velocity is higher than attached flow.
Problem 5.47
From the jet total pressure ptj=3(100 kPa)=300 kPa and Ttj=4 (288 K)= 1152 K, the jet
Mach number Mj=1.0 and mass flow rate, we can use continuity equation to calculate the
jet area, Aj.
1
1 2 2 ( 1)
p t AM 1 M
2
m
( 1)c p Tt
Therefore
1
1 2 2 ( 1)
Am
1 M ( 1)c p Tt /( p t M ) Aj=27.98 cm2
2
Therefore A=4 Aj A=111.92 cm2
p1 A1 (1 M 1 ) p j A j (1 M j ) p 0 A(1 M 2 )
2 2 2
s c p1Tt1 m
m j c pj Ttj (m
s m
j )c p 2Tt 2
From isentropic expansion between stations 0 and 1 we can relate static pressure and
temperature at 1 to M1 and p0 and T0, i.e.,
p0
p1 /( 1)
1 2
1 M1
2
T0
T1
1 2
1 M1
2
Now, we set up a spreadsheet calculation that assumes an M1 and proceeds to calculate all
the remaining unknowns, namely M2, T2 etc. then we check to see if the mass flow arte
out of the box is equal to the mass flow arte intro the box (continuity check). Then
iterate, until the mass is balanced.
T0 (K) p0 (kPa) Ttj (K) ptj (kPa) Aj (m2) M-j p-j (kPa)
288 100 1152 300 0.002798 1 158.48
a-j
T-j (K) rho-j (kg/m3) (m/s) V-j (m/s) m-dot-j (kg/s) I-j (kN) F-gross (N)
960 0.57551 620.9 620.9 1 1.064426 784.5815
Here only the M1 choice that balanced the mass flow rate is listed.
T-2 (K) rho-2 (kg/m3) a-2 (m/s) V-2 (m/s) m-dot (kg/s)-new
572.14 0.60930 479.3 362.6 2.4734
Percent increase in mass flow rate = 100(2.471-1)/1= 147.1% mass flow rate increase
Percent drop in exit temperature = 100(960-572.14)/969= 40.4% drop in jet exit temp.
Gross thrust with ejector = 896.2 N
Gross thrust of the jet = 784.6 N
Caution: the static pressure loss in the ejector due to wall friction and total pressure loss
due to turbulent jet mixing were not included in the analysis.
Problem 5.48
pt9/p9=(pt9/pt7)(pt7/p0)(p0/p9)=(0.96)(6)(1/1.5)=3.84
T9s=Tt7(1/4)0.3/1.3≈0.7262 Tt7
T9=Tt9/(1+0.15M92) ≈0.733 Tt9 note that Tt7=Tt9
Tt 7 T9 1 0.733
n 0.975
Tt 7 T9 s 1 0.7262
Problem 5.49
a) From area ratio A9/Ath and isentropic table (in supersonic section), we get M9=3.35,
which also gives pt9/p9=61.5
pt/p0=61.5
From the shock jump condition in static pressure, we calculate the NPR with an OS
Problem 5.50
pt0/p0=[1+0.15(36)]1.3/0.3 ≈ 3115
pt4/p0= (pt4/pt2)( pt2/pt0)( pt0/p0)= 0.42(0.5)(3115) ≈654
Therefore a Mach-6 scramjet (with then specified losses) experiences an NPR of ~ 654
Problem 5.51