From Bernoulli's principle to hypersonic shock physics: a complete engineering reference to aerodynamics covering subsonic, transonic, supersonic, and hypersonic regimes.
Every aircraft in steady level flight is in equilibrium under lift, weight, thrust, and drag. Understanding how these forces arise and interact is the foundation of all aerodynamics.
Lift is the aerodynamic force perpendicular to the freestream velocity. It arises from a net pressure difference between the upper and lower wing surfaces. The Kutta-Joukowski theorem states L′ = ρV∞Γ, where Γ is the circulation around the aerofoil. The Kutta condition requires flow to leave the trailing edge smoothly, uniquely determining Γ and therefore lift.
In practice, the upper surface of a cambered wing accelerates flow (lower pressure) while the lower surface decelerates it (higher pressure). This pressure asymmetry integrated over the chord produces lift. Lift coefficient CL = L / (½ρV²S) increases linearly with angle of attack until stall.
Total drag is the sum of parasite drag (profile drag, skin friction, interference) and induced drag. The drag polar is:
CD = CD₀ + CL² / (π · AR · e)
Parasite drag CD₀ dominates at high speed. Induced drag dominates at low speed and high lift, it arises because a finite wing sheds trailing vorticity all along its span wherever the bound circulation changes, dΓ/dy ≠ 0, not only at the tips. That sheet rolls up into the two concentrated tip vortices you can see, and the downwash it induces across the span tilts the local lift vector rearward. Minimum induced drag for a given lift and span corresponds to elliptical loading; e measures how far a real wing falls short of it. Winglets reduce this by limiting vortex strength. High aspect ratio wings (large span relative to chord) minimise induced drag for a given lift coefficient.
The angle of attack (AoA) is the angle between the chord line and the incoming freestream. Increasing AoA raises CL linearly at a rate of 2π per radian (0.110 per degree) for a thin aerofoil in two-dimensional flow. A finite wing is always shallower, because trailing vorticity induces downwash that reduces the effective incidence:
a = a₀ / (1 + a₀/(π·AR·e))
An AR 9.5 airliner wing lands near 0.088 per degree, about 20% under the 2D figure. Use 2π/rad on a real wing and your lift comes out a quarter too high.
At the critical stall angle the boundary layer separates and CL drops abruptly. There is no single stall angle: a thick, round-nosed section typically reaches 15 to 18° through progressive trailing-edge separation, while a thin sharp-nosed section stalls earlier and far more abruptly at 10 to 12° via leading-edge or thin-aerofoil stall (McCullough & Gault, NACA TN 2502). Which mechanism you get depends on leading-edge radius, thickness and Reynolds number, and it determines whether the stall is gentle or sudden. Stall speed Vs = √(2W / ρSCLmax). High-lift devices (flaps, slats) increase CLmax and reduce Vs.
An aerofoil is characterised by its chord length, camber (maximum distance of the mean camber line from the chord, as a % of chord), and thickness (maximum thickness as a % of chord). The NACA 4-digit series encodes these: NACA 2412 has 2% maximum camber at 40% chord and 12% maximum thickness.
Camber makes the aerofoil generate lift at zero AoA. Thickness affects the pressure distribution and maximum lift. The leading edge radius determines stall behaviour, sharper leading edges produce more abrupt stalls. The NACA 6-series (e.g. NACA 63-215) was designed for delayed transition to turbulent flow, reducing skin friction drag.
The boundary layer is where all skin friction drag, stall, and separation originate. Understanding it is essential for any aerodynamic design.
In laminar flow the fluid moves in smooth parallel layers with no lateral mixing. Skin friction coefficient Cf is low. Laminar boundary layers are fragile. Transition is not set by Reynolds number alone: pressure gradient (favourable stabilises, adverse destabilises), freestream turbulence intensity, surface roughness and wall temperature all move it, and on a swept wing crossflow instability usually dominates. The Re ≈ 5×10⁵ figure is a smooth flat-plate rule of thumb, not a regime boundary, and airliner wings routinely hold laminar flow well past it. Laminar aerofoils like the NACA 6-series delay this transition to achieve lower cruise drag.
Turbulent boundary layers have chaotic velocity fluctuations and strong lateral mixing. Skin friction is higher but the layer is more resistant to separation. The additional momentum from the turbulent outer layer energises the near-wall flow, allowing it to survive adverse pressure gradients that would separate a laminar layer.
When the pressure gradient becomes sufficiently adverse (pressure rising in the flow direction), boundary layer fluid loses momentum and the flow separates from the surface. This creates a large low-pressure wake behind the body, massively increasing pressure drag. On a wing, trailing-edge separation growing toward the leading edge is the mechanism behind the stall.
The Reynolds number Re = ρVL/μ is the ratio of inertial to viscous forces. It determines whether the boundary layer is laminar or turbulent, governs scale effects in wind tunnel testing, and is the primary reason aircraft aerodynamics cannot simply be scaled down to model size without correction. Commercial aircraft wings operate at Re ≈ 10⁷ to 10⁸. Model aircraft operate at Re ≈ 10⁴ to 10⁵. The flow physics are fundamentally different at these scales.
As flight speed increases toward and beyond the speed of sound, compressibility effects dominate aerodynamics. New phenomena appear: wave drag, shock waves, and thermal loads that have no subsonic equivalent.
The Mach number M = V/a (ratio of flow speed to local speed of sound) defines the flow regime. Below M = 0.3 compressibility is negligible. From M = 0.8 to 1.2 the flow is transonic, local shock waves appear on the upper wing surface even though the aircraft flies below Mach 1, causing wave drag and potentially buffet. Above M = 1.2 the flow is supersonic, and above M = 5 it is hypersonic, where dissociation and ionisation of air molecules become significant.
When flow accelerates past Mach 1, pressure disturbances cannot propagate upstream. They coalesce into shock waves: thin discontinuities across which pressure, temperature, and density jump abruptly. A normal shock is perpendicular to the flow and always decelerates the flow to subsonic; total pressure is lost irreversibly. An oblique shock forms at an angle and can leave the flow supersonic. Shocks are the dominant source of drag on supersonic aircraft.
The critical Mach number Mcrit is the freestream Mach at which local flow first reaches M = 1 on the upper wing surface. For a typical commercial wing this is around M = 0.72 to 0.78. Beyond Mcrit, wave drag increases steeply and buffet may onset. Swept wings delay Mcrit by reducing the effective chord-wise Mach component. The drag divergence Mach number (where CD increases by 0.002 above the low-speed value) is typically 0.04 to 0.08 above Mcrit.
For isentropic (adiabatic, reversible) flow the total conditions (total pressure p₀, total temperature T₀) are constant along a streamline. The Prandtl-Glauert compressibility correction adjusts subsonic pressure coefficients: Cp = Cp_incomp / √(1 − M²). This correction breaks down near M = 1. The Karman-Tsien and Laitone rules provide more accurate corrections at higher subsonic Mach numbers.
Shock waves are not just "loud bangs". They are physical boundaries across which everything changes instantly. Pressure, temperature, density, and velocity all jump discontinuously. Understanding the types, causes, and consequences of shock waves is fundamental to supersonic and hypersonic design.
At subsonic speeds, air molecules ahead of an aircraft receive a pressure warning and move out of the way. At supersonic speeds, the aircraft is moving faster than those pressure signals can travel, so the air gets no warning at all.
The result: air molecules are suddenly compressed into a thin layer, the shock wave. Across this boundary, kinetic energy is irreversibly converted to heat. Pressure and density jump up sharply. Velocity drops sharply. This energy conversion is permanent, unlike isentropic compression, it cannot be recovered. This is why wave drag exists and why supersonic flight is inefficient.
Normal shock wave, stands perpendicular to the flow. The flow goes from supersonic to subsonic across it. Always. No exceptions. Maximum total pressure loss, maximum entropy increase. Found at the front of blunt bodies and inside engine inlets.
Oblique shock wave, stands at an angle to the flow. The flow can remain supersonic after passing through it (weak shock) or become subsonic (strong shock). Less total pressure loss than a normal shock. Found on sharp-nosed aircraft, wedges, and compression ramps.
The key difference: a normal shock is the most violent, highest temperature rise, highest drag penalty. Oblique shocks are "gentler" and are what aerodynamicists try to use instead whenever possible.
A strong oblique shock occurs at steep shock angles. The flow crosses from supersonic to subsonic, just like a normal shock, but at an angle. It has a large deflection angle and a large total pressure loss. These are two separate phenomena and they are commonly confused. For a given M₁ and a deflection θ below the maximum θ_max, the θ-β-M relation admits two solutions, weak and strong; the weak one is what you observe in most external flows, because the strong solution requires a downstream back pressure that an external flow does not impose. Ducts and inlets, which do impose back pressure, can and do run the strong solution. Detachment is different again: it happens when θ exceeds θ_max, and the shock stands off the body as a curved bow shock with a subsonic pocket behind its normal portion. Designers avoid the strong solution in external flow because the total pressure loss and drag penalty are severe.
A weak oblique shock occurs at shallow shock angles. The flow remains supersonic after passing through it, but at a lower Mach number and a slightly different direction. It has a small total pressure loss. This is what engineers design for, sharp-nosed aircraft and pointed compression surfaces deliberately create weak oblique shocks to minimise drag while still compressing the flow. The Concorde intake used a series of weak oblique shocks for exactly this reason.
Aircraft motion is described about three axes. Each axis has a primary control surface and a stability characteristic that determines how the aircraft responds to disturbances.
Rotation about the lateral axis (wing tip to wing tip). Positive pitch is nose-up. Controlled by the elevator or stabilator. Longitudinal stability is provided by the horizontal tail. The tail volume coefficient determines the restoring moment back to trimmed AoA after a disturbance.
Rotation about the longitudinal axis (nose to tail). Controlled by ailerons acting differentially, one down (more lift), one up (less lift). Dihedral (wings angled upward) provides natural roll stability. Roll generates adverse yaw on the rising wing, which must be corrected by the rudder for coordinated flight.
Rotation about the normal axis (top to bottom). Controlled by the rudder on the vertical tail fin. Directional stability comes from the fin acting as a weathervane. The rudder must counteract adverse yaw from ailerons and correct for engine-out asymmetric thrust on multi-engine aircraft.
Drag the sliders to change angle of attack and Mach number: the lift, drag, and pressure distribution update instantly.
Thin-aerofoil theory lives in a 2D slice and knows nothing about span. Everything that makes a real wing different happens in the third dimension: circulation varies along the span, vorticity trails off behind it, and the downwash that creates tilts the lift vector backwards into drag. Drag the model to rotate it.
This is not a canned curve. It solves Prandtl's monoplane equation for the circulation distribution on every frame, expanding Γ as a Fourier sine series in θ where y = (b/2)cos θ, then solving the resulting system for the odd coefficients An:
Σ An sin nθ (sin θ + nμ) = μ α sin θ
with μ = a₀c(θ)/8s. Everything else falls out: CL = π·AR·A₁, CDi = π·AR·Σ n An², and e = 1/(1+δ) where δ = Σn>1 n(An/A₁)². Change the taper and you change c(θ), which changes every coefficient.
A 2D aerofoil has infinite span, so no vorticity trails and there is no downwash. Give the wing an end and circulation must fall to zero there. Wherever dΓ/dy is non-zero the wing sheds a trailing vortex filament, and the sheet of filaments induces a downward velocity at the wing itself.
That downwash reduces the angle the aerofoil actually sees from α to α − εi, which is why the measured lift slope is always below 2π/rad. It also tilts the local lift vector backwards by εi, and the rearward component of that vector is induced drag. Induced drag is not friction: it is the price of making lift with a wing of finite span, and it would exist in a perfectly inviscid fluid.
How thrust is generated: turbofans, rockets, and the thermodynamic cycles that power every aircraft.
Simulate the aerodynamic forces you've just learned about. Pressure fields, velocity contours, and turbulence models.
How lift and drag translate into stability, control, and aircraft performance: pitch, roll, yaw, and trim.
Aerodynamics doesn't stop at the atmosphere: fins, drag, and stability govern rockets all the way to orbit.
New to aerospace engineering? Start with the big picture: disciplines, career paths, and how everything connects.
Aerofoil optimisation, turbulence transition studies, and supercritical wing analysis: aerodynamics project ideas with full methodologies.
SheCodes Lab teaches Python and C++ from scratch, side by side, free, no experience needed. Includes an engineering module covering NumPy, pandas, ISA models, cost index, and flight data analysis. The same tools used to build the calculators on this site.
shecodeslab.com →