Critical & Drag-Divergence Mach Number
Unified master study module covering Critical Mach Number ($M_{cr}$), Drag-Divergence Mach Number ($M_{dd}$), transonic shock wave formation, wave drag divergence, Whitcomb area rule, supercritical airfoils, and high-speed Mach tuck safety implications for the UPSC Air Safety Officer examination.
Standard Engineering Criteria for $M_{dd}$
- Boeing / Douglas Criterion: $M_{dd}$ is the Mach number at which the rate of change of drag coefficient with Mach number reaches $ rac{dC_D}{dM} = 0.10$.
- Alternative Standard Criterion: $M_{dd}$ is the Mach number at which the total drag coefficient increases by $Delta C_D = 0.0020$ (20 drag counts) above its low-speed subsonic baseline ($C_{D0}$).
Transonic Milestone Sequence
For any conventional or supercritical airfoil, the Mach number milestones always occur in strict hierarchical order:
Side-by-Side Comparison: $M_{cr}$ vs $M_{dd}$
| Parameter / Aspect | Critical Mach Number ($M_{cr}$) | Drag-Divergence Mach Number ($M_{dd}$) |
|---|---|---|
| Physical Trigger | Local $M_{max} = 1.0$ at point of peak velocity | Rapid formation of strong normal shock & boundary layer separation |
| Drag Behavior | Drag coefficient remains nearly constant ($C_D approx C_{D0}$) | $C_D$ experiences steep exponential rise ($ rac{dC_D}{dM} ge 0.10$) |
| Shock Waves | No shock wave yet; incipient sonic point | Distinct, oscillating normal shock waves on upper & lower surfaces |
| Typical Commercial Value | $0.70 le M_{cr} le 0.76$ (conventional) | $0.78 le M_{dd} le 0.86$ (supercritical) |
| Offset ($Delta M = M_{dd} - M_{cr}$) | Typically $Delta M approx 0.06 ext{ to } 0.10$ for high-performance transport wings | |
1. Transonic Flowfield Evolution
2. Aerodynamic Design Strategies to Increase $M_{cr}$ and $M_{dd}$
- Wing Sweepback ($Lambda$): Sweeping the wing by angle $Lambda$ reduces the effective normal Mach number to $M_n = M_infty cosLambda$. This delays both $M_{cr}$ and $M_{dd}$ by a factor of approximately $ rac{1}{cosLambda}$.
- Thin Airfoil Sections (Reduced $t/c$): Decreasing thickness-to-chord ratio from 15% to 9% significantly reduces supervelocities over the crest, increasing $M_{cr}$.
- Supercritical Airfoil Design: Features a flattened upper surface (limiting supersonic expansion to low Mach numbers over a broad chord), aft camber (cusped lower trailing edge) to restore lift, and a larger leading-edge radius to prevent leading-edge stall.
- Whitcomb Transonic Area Rule: For minimum wave drag near Mach 1.0, the total cross-sectional area distribution $A(x)$ of the aircraft along the longitudinal axis must be smooth and continuous (Coke-bottle fuselage waisting).
1. Subsonic Compressibility Corrections ($C_p$ vs $C_{p0}$)
In subsonic compressible flow ($0.3 < M_infty < M_{cr}$), surface pressure coefficients are magnified relative to incompressible values ($C_{p0}$):
2. Derivation of Critical Pressure Coefficient ($C_{p,cr}$)
Using the isentropic pressure-Mach relation for an ideal gas ($gamma = 1.4$):
$$ rac{p_0}{p_infty} = left(1 + rac{gamma - 1}{2} M_infty^2 ight)^{ rac{gamma}{gamma - 1}}$$ $$ rac{p_0}{p^*} = left(1 + rac{gamma - 1}{2} (1)^2 ight)^{ rac{gamma}{gamma - 1}} = left( rac{gamma + 1}{2} ight)^{ rac{gamma}{gamma - 1}}$$Dividing the two equations yields the critical pressure ratio at sonic speed:
$$ rac{p^*}{p_infty} = left[ rac{2 + (gamma - 1)M_infty^2}{gamma + 1} ight]^{ rac{gamma}{gamma - 1}}$$Substituting into the definition of pressure coefficient $C_p = rac{p - p_infty}{ rac{1}{2} gamma p_infty M_infty^2}$ gives the exact formula for $C_{p,cr}$:
For air ($gamma = 1.4$):
Graphical Solution for $M_{cr}$: Plot the Prandtl-Glauert curve $C_p(M_infty) = rac{C_{p0,min}}{sqrt{1 - M_infty^2}}$ and the critical curve $C_{p,cr}(M_infty)$. Their exact intersection gives $M_{cr}$.
3. Wing Sweep Delay Derivation
For a wing swept by angle $Lambda$:
$$M_n = M_infty cosLambda$$ $$M_{cr, ext{swept}} = rac{M_{cr, ext{unswept}}}{cosLambda} quad ext{and} quad M_{dd, ext{swept}} = rac{M_{dd, ext{unswept}}}{cosLambda}$$Wave drag scaling is reduced by the cube of the cosine: $C_{D,w, ext{swept}} = C_{D,w, ext{unswept}} cos^3Lambda$.
Near Mach 1.0, the linearized potential equation fails because transonic terms are non-linear. The transonic small-disturbance (TSD) equation introduces the Transonic Similarity Parameter ($K$ or $chi$):
Where $ au = t/c$ is the maximum thickness-to-chord ratio. Two affinity flows with identical $K$ have geometrically scaled pressure fields:
Crucial Examination Takeaway: Wave drag in the transonic regime scales as $ au^{5/3}$ (or $approx (t/c)^2$). Doubling airfoil thickness increases transonic wave drag by nearly 400%!
Key Aerodynamic Phenomenon in ASO Investigations
- Cause: When $M_infty > M_{cr}$, the normal shock forms and moves aft. The center of pressure shifts from $0.25c$ (subsonic) towards $0.50c$ (supersonic), creating a strong nose-down pitching moment ($Delta C_{m,cg} < 0$). Simultaneously, the downwash angle $epsilon$ behind the wing decreases, reducing tailplane downward download and exacerbating the dive.
- Safety Mitigation: Automatic Mach Trim systems that schedule incremental nose-up horizontal stabilizer trim as Mach number increases above $M_{cr}$.
- At high cruising altitudes (FL370 - FL450), the margin between $V_S$ (low-speed stall where $C_L = C_{L,max}$) and $M_{buffet}$ (high-speed shock-induced separation) contracts to as little as 10-15 knots.
- Any sudden turbulence, bank angle increase ($n = 1/cosphi$), or thermal inversion can pitch the aircraft simultaneously into a low-speed stall and high-speed shock stall.
Engineering Data Table: Transonic Aircraft Profiles
| Aircraft / Airfoil Type | Thickness Ratio ($t/c$) | Sweep Angle ($Lambda_{c/4}$) | Estimated $M_{cr}$ | Certified $M_{MO}$ | Estimated $M_{dd}$ |
|---|---|---|---|---|---|
| NACA 0012 (Unswept Baseline) | 12% | 0° | 0.70 | N/A | 0.74 |
| NACA 64A-010 (Laminar) | 10% | 0° | 0.76 | N/A | 0.80 |
| Boeing 737 Next Gen | 10.5% | 25° | 0.74 | 0.82 | 0.84 |
| Airbus A320neo | 11.0% (Supercritical) | 25° | 0.76 | 0.82 | 0.85 |
| Boeing 787 Dreamliner | 10.0% (Supercritical) | 32.2° | 0.81 | 0.89 | 0.91 |
Pillar 6: Aircraft Airworthiness, Systems & DGCA CAR / ICAO Compliance
Pillar 7: High-Frequency UPSC Examiner Traps & Misconceptions
❌ Examiner Trap: Assuming that a normal shock wave already exists when an aircraft reaches $M_{cr}$.
✅ UPSC Truth: At $M_{cr}$, the flow is purely sonic ($M_{local} = 1.0$) at only a single point on the suction surface. Shock waves only form at speeds above $M_{cr}$ ($M_infty > M_{cr}$).
❌ Examiner Trap: Treating $M_{cr}$ and $M_{dd}$ as interchangeable.
✅ UPSC Truth: $M_{cr}$ is the onset of local sonic speed with virtually no drag penalty ($Delta C_D approx 0$). $M_{dd}$ occurs at a higher Mach number ($M_{dd} > M_{cr}$) where drag rises precipitously ($ rac{dC_D}{dM} ge 0.10$ or $Delta C_D = 0.0020$).
❌ Examiner Trap: Assuming wave drag scales linearly with thickness ratio ($t/c$).
✅ UPSC Truth: Transonic wave drag scales non-linearly as $(t/c)^{5/3}$ or approximately $(t/c)^2$. Doubling airfoil thickness quadruples wave drag.
Pillar 8: Step-by-Step Worked Aircraft Numerical Problem
A jet transport aircraft is designed with a straight wing having a critical Mach number $M_{cr} = 0.72$ and an incompressible minimum pressure coefficient $C_{p0} = -0.55$. To increase the high-speed cruise performance, aerodynamicists apply a quarter-chord wing sweep of $Lambda = 30^circ$. Calculate: (a) the compressible pressure coefficient at $M_infty = 0.60$ using the Prandtl-Glauert rule, and (b) the new critical Mach number with the swept wing.
Baseline $C_{p0} = -0.55$
Sweep angle $Lambda = 30^circ$ ($cos 30^circ = 0.8660$)
Flight Mach $M_infty = 0.60$
(a) The compressible peak suction coefficient increases in magnitude to $C_p = -0.688$.
(b) Sweeping the wing by $30^circ$ increases the critical Mach number from $0.720 o 0.831$, allowing an additional 60+ knots of safe high-speed cruising speed without encountering shock buffet.
High-Yield Practice MCQs: Critical & Drag-Divergence Mach Number
Curated directly from UPSC Engineering Services, GATE Aerospace, and DGCA Air Safety Officer past question papers.
Valid up to $M_infty approx M_{cr}$. Singular at $M_infty = 1.0$.