Interleaved — Phase 3

Physics interleaved practice

printable — key stays hidden on paper

Instructions: Work each problem fully. Use γ=1.4\gamma = 1.4 and R=287 J/(kgK)R = 287\ \mathrm{J/(kg\cdot K)} for air unless stated otherwise. Show method selection reasoning. Marks in brackets. No calculators of shock tables provided — use the closed-form relations where possible.


1. [8] Air flows at V=250 m/sV = 250\ \mathrm{m/s} with static temperature T=280 KT = 280\ \mathrm{K}. Compute (a) the local speed of sound, (b) the Mach number, (c) the stagnation temperature T0T_0, and (d) classify the regime (subsonic/transonic/supersonic/hypersonic).

2. [6] Starting from the steady-flow energy equation for an adiabatic, no-work open system, derive the stagnation temperature relation T0T=1+γ12M2\dfrac{T_0}{T} = 1 + \dfrac{\gamma-1}{2}M^2. State every assumption.

3. [7] A normal shock stands in a duct with upstream Mach number M1=2.0M_1 = 2.0. Find M2M_2, the static pressure ratio P2/P1P_2/P_1, and the stagnation pressure ratio P02/P01P_{02}/P_{01}.

4. [8] A converging–diverging nozzle has throat area A=10 cm2A^* = 10\ \mathrm{cm}^2. At a section where A=15 cm2A = 15\ \mathrm{cm}^2 the flow is supersonic. Using the area–Mach relation, find the Mach number there (the supersonic root), and hence P/P0P/P_0 at that section.

5. [7] A supersonic stream at M1=3.0M_1 = 3.0 expands isentropically around a convex corner, turning through a deflection of Δθ=15\Delta\theta = 15^\circ. Using the Prandtl–Meyer function, find the downstream Mach number M2M_2.

6. [6] Derive the area–velocity relation dAA=(M21)dVV\dfrac{dA}{A} = (M^2 - 1)\dfrac{dV}{V} from continuity, the isentropic momentum (Euler) relation, and the definition of sound speed. Use it to explain physically why a de Laval nozzle needs a diverging section to accelerate flow beyond M=1M=1.

7. [7] A symmetric thin airfoil (zero camber) operates at angle of attack α=4\alpha = 4^\circ in incompressible flow. (a) Using thin-airfoil theory, compute the lift coefficient cc_\ell. (b) A finite wing of aspect ratio AR=6AR = 6 and elliptic loading has the same cc_\ell; compute its induced drag coefficient cD,ic_{D,i}.

8. [7] A converging (only) nozzle is fed from a reservoir at P0=500 kPaP_0 = 500\ \mathrm{kPa}, T0=300 KT_0 = 300\ \mathrm{K}, and discharges to a back pressure PbP_b. (a) Determine the back pressure at which the nozzle first becomes choked. (b) If PbP_b is lowered further, what happens to the exit Mach number and mass flow rate, and why?

9. [6] An oblique shock forms on a wedge in a M1=2.5M_1 = 2.5 flow with wave angle β=40\beta = 40^\circ. (a) Find the normal component Mn1M_{n1}. (b) Find the deflection angle θ\theta from the θ\thetaβ\betaMM relation. (c) State what happens to β\beta if the wedge half-angle exceeds θmax\theta_{max}.

10. [6] A de Laval nozzle is designed for perfectly expanded exit Mach Me=2.4M_e = 2.4. It is operated at a back pressure higher than the design exit pressure but low enough that the flow remains supersonic just inside the exit. (a) Name the operating condition (over- or under-expanded). (b) Describe the wave structure at the exit plane and the sign of PePbP_e - P_b.

Answer keyMark scheme & solutions

1. Tests 3.1.3 (sound speed) + 3.1.4 (Mach) + 3.1.2 (stagnation T). Why: Given VV and TT directly — you must build aa first, then MM, then T0T_0. Recognizing the chain (no area/nozzle info) is the method choice.

(a) a=γRT=1.4287280=112476=335.4 m/sa = \sqrt{\gamma R T} = \sqrt{1.4 \cdot 287 \cdot 280} = \sqrt{112476} = 335.4\ \mathrm{m/s}. (b) M=V/a=250/335.4=0.745M = V/a = 250/335.4 = 0.745. (c) T0=T(1+γ12M2)=280(1+0.20.555)=280(1.111)=311.1 KT_0 = T\left(1 + \tfrac{\gamma-1}{2}M^2\right) = 280(1 + 0.2\cdot0.555) = 280(1.111) = 311.1\ \mathrm{K}. (d) M<1M < 1subsonic.


2. Tests 3.1.1 (first law open system) + 3.1.2 (stagnation derivation). Why: Pure derivation — must start from energy conservation, not memorized formula.

Steady-flow energy equation for adiabatic (q=0q=0), no shaft work (w=0w=0), negligible gravity: h+12V2=h0=const.h + \tfrac{1}{2}V^2 = h_0 = \text{const}. Stagnation state defined by V=0V=0: h0=h+12V2h_0 = h + \tfrac12 V^2. For a calorically perfect gas h=cpTh = c_p T, cp=γRγ1c_p = \dfrac{\gamma R}{\gamma-1}: cpT0=cpT+12V2T0T=1+V22cpT.c_p T_0 = c_p T + \tfrac12 V^2 \Rightarrow \frac{T_0}{T} = 1 + \frac{V^2}{2 c_p T}. Since cpT=γRTγ1=a2γ1c_p T = \dfrac{\gamma R T}{\gamma-1} = \dfrac{a^2}{\gamma-1} and M2=V2/a2M^2 = V^2/a^2: T0T=1+V2(γ1)2a2=1+γ12M2.  \frac{T_0}{T} = 1 + \frac{V^2(\gamma-1)}{2 a^2} = 1 + \frac{\gamma-1}{2}M^2. \;\blacksquare Assumptions: adiabatic, no work, steady, calorically perfect gas, negligible potential energy.


3. Tests 3.1.12 (normal shock properties) + 3.1.11 (R-H relations). Why: "Normal shock, given M1M_1" triggers closed-form R–H relations, not isentropic tables (stagnation P is lost).

M22=1+γ12M12γM12γ12=1+0.241.440.2=1.85.4=0.3333M2=0.577M_2^2 = \dfrac{1 + \frac{\gamma-1}{2}M_1^2}{\gamma M_1^2 - \frac{\gamma-1}{2}} = \dfrac{1 + 0.2\cdot4}{1.4\cdot4 - 0.2} = \dfrac{1.8}{5.4} = 0.3333 \Rightarrow M_2 = 0.577.

P2P1=1+2γγ+1(M121)=1+2.82.4(3)=1+3.5=4.5\dfrac{P_2}{P_1} = 1 + \dfrac{2\gamma}{\gamma+1}(M_1^2 - 1) = 1 + \dfrac{2.8}{2.4}(3) = 1 + 3.5 = 4.5.

Stagnation pressure ratio: P02P01=[(γ+1)M12/21+γ12M12]γγ1[2γM12(γ1)γ+1]1γ1.\frac{P_{02}}{P_{01}} = \left[\frac{(\gamma+1)M_1^2/2}{1+\frac{\gamma-1}{2}M_1^2}\right]^{\frac{\gamma}{\gamma-1}}\left[\frac{2\gamma M_1^2 - (\gamma-1)}{\gamma+1}\right]^{-\frac{1}{\gamma-1}}. Term A =(2.44/21.8)3.5=(4.81.8)3.5=(2.6667)3.5=30.9=\left(\dfrac{2.4\cdot4/2}{1.8}\right)^{3.5}=\left(\dfrac{4.8}{1.8}\right)^{3.5}=(2.6667)^{3.5}=30.9. Term B =(11.20.42.4)2.5=(4.5)2.5=0.0233=\left(\dfrac{11.2-0.4}{2.4}\right)^{-2.5}=(4.5)^{-2.5}=0.0233. P02/P01=30.90.0233=0.720P_{02}/P_{01} = 30.9\cdot0.0233 = 0.720. (Standard table value ≈ 0.7209.)


4. Tests 3.1.6 (area–Mach) + 3.1.7 (isentropic P/P₀). Why: A/AA/A^* given → invert area–Mach; "supersonic" tells you which of two roots. Then use isentropic (no shock) for P/P0P/P_0.

AA=1.5=1M[2γ+1(1+γ12M2)]γ+12(γ1)\dfrac{A}{A^*} = 1.5 = \dfrac{1}{M}\left[\dfrac{2}{\gamma+1}\left(1+\tfrac{\gamma-1}{2}M^2\right)\right]^{\frac{\gamma+1}{2(\gamma-1)}}. Solving the supersonic branch: M1.854M \approx 1.854.

PP0=(1+γ12M2)γγ1=(1+0.23.437)3.5=(1.6875)3.5=0.1602\dfrac{P}{P_0} = \left(1+\tfrac{\gamma-1}{2}M^2\right)^{-\frac{\gamma}{\gamma-1}} = (1+0.2\cdot3.437)^{-3.5} = (1.6875)^{-3.5} = 0.1602.


5. Tests 3.1.16/3.1.17 (Prandtl–Meyer expansion). Why: "Expands around convex corner" = isentropic turning → use ν(M)\nu(M), not shock relations.

ν(M)=γ+1γ1tan1 ⁣γ1γ+1(M21)tan1 ⁣M21\nu(M) = \sqrt{\dfrac{\gamma+1}{\gamma-1}}\tan^{-1}\!\sqrt{\dfrac{\gamma-1}{\gamma+1}(M^2-1)} - \tan^{-1}\!\sqrt{M^2-1}. At M1=3M_1=3: ν(3)=49.76\nu(3) = 49.76^\circ. Downstream: ν(M2)=ν(M1)+Δθ=49.76+15=64.76\nu(M_2) = \nu(M_1) + \Delta\theta = 49.76 + 15 = 64.76^\circ. Inverting: M23.94M_2 \approx 3.94. (Flow accelerates, as expected for expansion.)


6. Tests 3.1.5 (area–velocity derivation). Why: Pure derivation; connects continuity + Euler + sound speed.

Continuity (log-differentiate ρAV=\rho A V = const): dρρ+dAA+dVV=0\dfrac{d\rho}{\rho} + \dfrac{dA}{A} + \dfrac{dV}{V} = 0. Euler (inviscid, isentropic): dP=ρVdVdP = -\rho V\,dV. Sound speed: a2=dP/dρdρ=dP/a2=ρVdV/a2a^2 = dP/d\rho \Rightarrow d\rho = dP/a^2 = -\rho V\,dV/a^2, so dρρ=M2dVV\dfrac{d\rho}{\rho} = -M^2\dfrac{dV}{V}. Substitute: M2dVV+dAA+dVV=0dAA=(M21)dVV.  -M^2\frac{dV}{V} + \frac{dA}{A} + \frac{dV}{V} = 0 \Rightarrow \frac{dA}{A} = (M^2-1)\frac{dV}{V}. \;\blacksquare Physics: For M<1M<1, (M21)<0(M^2-1)<0: to accelerate (dV>0dV>0) area must decrease (converging). For M>1M>1, (M21)>0(M^2-1)>0: to accelerate area must increase (diverging). Hence to pass through M=1M=1 (throat) and continue accelerating supersonically, the nozzle must first converge then diverge — the de Laval geometry.


7. Tests 3.1.21 (thin airfoil) + 3.1.22/3.1.23 (induced drag / AR). Why: Symmetric = zero camber so lift is purely α\alpha-driven; then finite-wing correction uses elliptic result.

(a) Thin airfoil: c=2παc_\ell = 2\pi\alpha (radians). α=4=0.0698 rad\alpha = 4^\circ = 0.0698\ \mathrm{rad}. c=2π(0.0698)=0.4386c_\ell = 2\pi(0.0698) = 0.4386. (b) Elliptic loading: cD,i=c2πAR=0.43862π6=0.192418.85=0.01021c_{D,i} = \dfrac{c_\ell^2}{\pi\,AR} = \dfrac{0.4386^2}{\pi\cdot6} = \dfrac{0.1924}{18.85} = 0.01021.


8. Tests 3.1.9 (converging nozzle) + 3.1.8 (choked flow). Why: Converging-only nozzle: exit Mach caps at 1; choke condition uses critical pressure ratio.

(a) Choking begins when Me=1M_e=1, i.e. Pe=PP_e = P^*: PP0=(2γ+1)γγ1=(0.8333)3.5=0.5283\dfrac{P^*}{P_0} = \left(\dfrac{2}{\gamma+1}\right)^{\frac{\gamma}{\gamma-1}} = (0.8333)^{3.5} = 0.5283. Pb=0.5283500=264.1 kPaP_b = 0.5283 \cdot 500 = 264.1\ \mathrm{kPa}. (b) Lowering PbP_b below 264264 kPa: exit stays at Me=1M_e=1 (cannot exceed 1 in converging duct), so mass flow rate stays constant (maximum, choked) and exit pressure remains PP^*; the further expansion occurs outside the nozzle. Exit Mach unchanged.


9. Tests 3.1.13/3.1.14 (oblique shock, θ–β–M). Why: Given β\beta and M1M_1 → forward θ–β–M gives deflection; normal component drives shock strength.

(a) Mn1=M1sinβ=2.5sin40=2.5(0.6428)=1.607M_{n1} = M_1\sin\beta = 2.5\sin40^\circ = 2.5(0.6428) = 1.607. (b) tanθ=2cotβM12sin2β1M12(γ+cos2β)+2\tan\theta = 2\cot\beta\dfrac{M_1^2\sin^2\beta - 1}{M_1^2(\gamma+\cos2\beta)+2}. M12sin2β=6.250.4132=2.583M_1^2\sin^2\beta = 6.25\cdot0.4132 = 2.583; numerator =2.5831=1.583= 2.583-1 = 1.583. cos2β=cos80=0.1736\cos2\beta = \cos80^\circ = 0.1736; denominator =6.25(1.4+0.1736)+2=6.25(1.5736)+2=11.835= 6.25(1.4+0.1736)+2 = 6.25(1.5736)+2 = 11.835. cot40=1.1918\cot40^\circ = 1.1918. $\tan\the