Compressible Flow & Aerodynamics
Level 2 — Recall & Standard Problems
Time Limit: 30 minutes Total Marks: 40
Use and for air unless otherwise stated.
Q1. Define the following (2 marks each): (a) stagnation temperature , (b) Mach number, (c) critical Mach number of an airfoil. (6 marks)
Q2. Air flows at a static temperature of . Calculate the speed of sound. If the flow velocity is , find the Mach number and classify the flow regime. (4 marks)
Q3. Starting from the energy equation for adiabatic flow, derive the stagnation-to-static temperature ratio (4 marks)
Q4. A reservoir has and . For a point in isentropic flow where , calculate the static temperature and static pressure . (5 marks)
Q5. State the area–velocity relation and use it to explain physically why a converging–diverging (de Laval) nozzle is required to accelerate a flow from subsonic to supersonic speed. (5 marks)
Q6. Define choked flow. State the condition at the throat of a nozzle at which choking occurs, and explain what "maximum mass flow" means in this context. (4 marks)
Q7. A normal shock occurs in air at upstream Mach number . Using the normal-shock relations, calculate: (a) the downstream Mach number , and (b) the static pressure ratio . (4 marks)
Q8. Sketch and label a cambered airfoil, marking the chord, camber line, leading edge, trailing edge, and angle of attack. Define the lift coefficient . (4 marks)
Q9. State Prandtl's result for the induced drag of a finite wing and explain, in one or two sentences, how increasing the aspect ratio affects induced drag. (4 marks)
End of paper
Answer keyMark scheme & solutions
Q1. (6 marks)
- (a) = temperature the flowing gas would attain if brought to rest adiabatically (isentropically); it accounts for both static temperature and kinetic energy: . (2)
- (b) Mach number , ratio of local flow speed to local speed of sound. (2)
- (c) Critical Mach number = the free-stream Mach number at which the flow first reaches locally somewhere on the airfoil surface. (2)
Q2. (4 marks)
- . (2)
- . (1)
- transonic/supersonic (close to 1, effectively transonic). (1)
Q3. (4 marks)
- Adiabatic energy eq: (const). (1)
- Divide by : . (1)
- Use and , so . (1)
- Hence . (1)
Q4. (5 marks)
- . (2)
- . (3)
Q5. (5 marks)
- State relation correctly. (1)
- Subsonic (): , so to increase () requires → area must converge. (2)
- Supersonic (): , so to increase requires → area must diverge. Sonic condition () forces (throat). Therefore accelerate subsonic→sonic in converging section, sonic→supersonic in diverging section → C-D nozzle. (2)
Q6. (4 marks)
- Choked flow: condition where mass flow rate through the nozzle reaches its maximum and becomes independent of downstream (back) pressure. (2)
- Occurs when at the throat (minimum area). (1)
- Once sonic at throat, further lowering back pressure cannot increase ; maximum mass flux is fixed by and throat area. (1)
Q7. (4 marks)
- . (2)
- . (2)
Q8. (4 marks)
- Correct sketch with chord (LE–TE straight line), camber line (mid-thickness curve), leading edge, trailing edge labelled. (2)
- Angle of attack = angle between chord and free-stream velocity. (1)
- (lift / dynamic pressure × reference area). (1)
Q9. (4 marks)
- Induced drag coefficient (elliptic loading): (general: ). (2)
- Increasing aspect ratio decreases induced drag (inversely proportional); higher AR wings have weaker tip vortices/smaller downwash. (2)
[
{"claim":"Speed of sound at 250 K is ~316.9 m/s","code":"a=sqrt(1.4*287*250); result = abs(float(a)-316.86) < 0.5"},
{"claim":"M=1.073 for V=340 at T=250K","code":"a=sqrt(1.4*287*250); M=340/a; result = abs(float(M)-1.073) < 0.01"},
{"claim":"Static T at M=2 with T0=500 is 277.8 K","code":"T0=500; T=T0/(1+0.2*4); result = abs(float(T)-277.78) < 0.1"},
{"claim":"Static P at M=2 with P0=600 kPa is ~76.7 kPa","code":"P0=600; P=P0/((1.8)**3.5); result = abs(float(P)-76.68) < 0.3"},
{"claim":"M2=0.577 for M1=2 normal shock","code":"M1=2; M2sq=(1+0.2*M1**2)/(1.4*M1**2-0.2); result = abs(float(sqrt(M2sq))-0.5774) < 0.001"},
{"claim":"P2/P1=4.5 for M1=2 normal shock","code":"M1=2; r=1+(2.8/2.4)*(M1**2-1); result = abs(float(r)-4.5) < 0.001"}
]