Module 1: Flight Aerodynamics
Bird flight spans an enormous range of modes — from the hovering figure-of-8 kinematics of a hummingbird to the months-long dynamic soaring of an albatross over the Southern Ocean. This module derives the fundamental aerodynamic forces from first principles and builds the quantitative framework for understanding each mode.
1. Lift: From Bernoulli to the Kutta Condition
Bernoulli's Equation
Along a streamline in steady, inviscid, incompressible flow, energy is conserved:
\[ p + \tfrac{1}{2}\rho u^2 + \rho g z = \text{const} \]
On the upper wing surface, flow accelerates (speed \(u_{upper} > U_\infty\)), so pressure drops. On the lower surface, flow decelerates, pressure rises. The net pressure difference integrated over the wing area gives lift.
Kutta–Joukowski Theorem and Circulation
The Kutta condition states that the flow must leave the wing smoothly at the trailing edge — the sharp edge cannot sustain infinite velocity. This forces the generation of bound circulation \(\Gamma\) around the wing section:
\[ \Gamma = \oint_C \mathbf{u} \cdot d\mathbf{l} \]
The Kutta–Joukowski theorem then gives lift per unit span:
\[ L' = \rho U_\infty \Gamma \]
This is exact for 2D inviscid flow around any shape. For a thin airfoil at angle of attack\(\alpha\): \(\Gamma = \pi c U_\infty \sin\alpha\) (thin airfoil theory, \(c\) = chord), giving \(dC_L/d\alpha = 2\pi\,\text{rad}^{-1}\).
Finite Wing: The Lift Equation
For a finite wing of span \(b\) and area \(S\), integrating the spanwise circulation and accounting for downwash from trailing vortices:
\[ L = \frac{1}{2} \rho U^2 S \, C_L \]
where \(C_L\) (lift coefficient) depends on wing geometry, angle of attack, and Re. Typical values: \(C_L \approx 0.4\) (cruise), \(C_L \approx 1.5\) (near stall).
Wing Cross-Section: Streamlines and Pressure Distribution
2. Drag: The Drag Polar
Total aerodynamic drag on a bird has three components:
Parasitic Drag
\( D_{par} = \tfrac{1}{2}\rho U^2 C_{D,par} S_{ref} \)
Skin friction + pressure drag on body and wing surfaces. Scales as U². Dominant at high speed.
Induced Drag
\( D_{ind} = \tfrac{L^2}{\tfrac{1}{2}\rho U^2 \pi e b^2} \)
Cost of generating lift — trailing vortices extract kinetic energy. Scales as U⁻². Dominant at low speed.
Profile Drag
\( D_{pro} = \tfrac{1}{2}\rho U^2 C_{D0} S \)
Pressure + friction drag of the wing itself moving at the local air speed (important in flapping flight).
The Drag Polar
The drag polar relates total drag coefficient to lift coefficient. For a parabolic polar (Prandtl lifting-line theory):
\[ C_D = C_{D0} + \frac{C_L^2}{\pi e \, AR} \]
\(C_{D0}\) = zero-lift (parasite) drag coefficient; \(e\) = Oswald span efficiency (0.8–0.95);\(AR = b^2/S\) = aspect ratio. The second term is induced drag, minimized at high AR.
Aspect Ratio: Why It Matters
| Bird | AR = b²/S | Best L/D | Flight strategy |
|---|---|---|---|
| Wandering Albatross | ~18 | ~60 | Ocean dynamic soaring, minimal flapping |
| Common Crane | ~9 | ~20 | Long-distance thermal soaring & flapping |
| Red-tailed Hawk | ~6 | ~12 | Thermal soaring, slow maneuvering |
| Rock Pigeon | ~6 | ~10 | Fast sustained flapping |
| House Sparrow | ~5 | ~8 | Short bursts, dense habitat maneuver |
| Hummingbird | ~3 | ~5 | Hovering, slow precise flight |
3. Gliding Flight: The Glide Polar
In steady gliding, thrust is zero and gravity provides the propulsive force. The glide angle\(\gamma\) satisfies:
\[ \tan\gamma = \frac{D}{L} = \frac{1}{L/D} \]
The glide ratio \(L/D\) equals the horizontal distance per unit altitude loss. An albatross with \(L/D = 60\) descends only 1 m for every 60 m forward.
Minimum Sink Rate
Sink rate \(w_s = U \sin\gamma \approx U \cdot D/L\). Substituting the drag polar and level-flight constraint \(L = W\):
\[ w_s = \frac{P_{mech}}{W} = \frac{\frac{1}{2}\rho U^3 C_D S}{W} \]
Minimum sink speed: \( U_{ms} = \left(\frac{2W}{\rho S}\right)^{1/2} \left(\frac{C_{D0}}{3 C_{Di}^*}\right)^{1/4} \)occurs at \(C_L^* = \sqrt{3\pi e\, AR\, C_{D0}}\)
Maximum Range Speed
Maximum range (horizontal distance per unit altitude) occurs at maximum \(L/D\), where:
\[ (L/D)_{max} = \frac{1}{2}\sqrt{\frac{\pi e\, AR}{C_{D0}}} \]
This occurs at \(C_L = \sqrt{\pi e\, AR\, C_{D0}}\) — the induced and parasite drags are equal.
4. Flapping Flight: Vortex Wake Models
Bound Circulation and Wake Shedding
As the wing flaps, the bound circulation \(\Gamma(y,t)\) varies along the span and in time. By the Kelvin circulation theorem, every change in bound circulation must be compensated by vorticity shed into the wake:
\[ \frac{d\Gamma_{bound}}{dt} + \frac{d\Gamma_{wake}}{dt} = 0 \]
Two competing wake models: (1) Continuous vortex sheet (Prandtl–Philips) — shed vorticity forms a continuous ribbon behind the wing. (2) Ring vortex model(Rayner 1979) — each downstroke sheds a discrete closed vortex ring of momentum \(J = \rho \Gamma A_{ring}\).
Actuator Disk Model (Hover / Slow Flight)
In the actuator disk approximation, the wing sweep area \(A_{disk} = \pi b^2/4\) acts as a propeller disk accelerating air downward. By momentum theory:
\[ P_{ind} = \frac{(mg)^{3/2}}{\sqrt{2\rho A_{disk}}} \]
Induced power at hover — the theoretical minimum. Real birds have additional profile power from wing drag, typically increasing total hover cost by 30–50%.
Hummingbird Hovering: Figure-of-8 Kinematics
Unlike insects, hummingbirds generate lift on both the downstroke AND the upstroke. The wing traces a horizontal figure-of-8 when viewed from the side:
- Downstroke: wing sweeps forward with leading edge up, \(\alpha \approx +25°\), generating upward lift
- Upstroke: wing rotates 180° (supinated), sweeps back with leading edge still leading airflow, \(\alpha \approx +15°\), generating upward lift (60–70% of downstroke value)
- Stroke reversal: rapid wing rotation (<10 ms) involving wake capture — the wing re-encounters the vortex from its previous stroke and extracts additional lift
Why Hummingbirds Are Unique Among Birds
Most birds generate negligible lift on the upstroke — the wing is partially folded to reduce drag. Hummingbirds evolved a ball-and-socket shoulder joint with rotational freedom 50× greater than other birds, enabling the full upstroke inversion. This allows 100% utilisation of the stroke cycle for lift generation — a necessity for hovering, where there is no forward speed to provide translational lift.
5. The U-Shaped Power Curve
Total mechanical power for forward flight is the sum of three components, each with a different speed dependence:
\[ P_{total} = P_{ind} + P_{pro} + P_{par} \]
\(P_{ind} \propto U^{-1}\)
Induced: decreases with speed (less induced downwash needed at high speed)
\(P_{pro} \approx \text{const}\)
Profile: approximately speed-independent in basic models
\(P_{par} \propto U^3\)
Parasite: increases steeply with speed (dominant at high speed)
The superposition of \(U^{-1}\) and \(U^3\) terms creates the characteristic U-shaped power curve. This has two important optima:
Minimum Power Speed \(U_{mp}\)
Speed at which \(P_{total}\) is minimized — the bird uses minimum power per unit time. Relevant for maximum hover endurance, minimum metabolic cost (e.g., a sitting bird is just above \(U_{mp}\)). Typically \(\approx 0.76 U_{mr}\).
Maximum Range Speed \(U_{mr}\)
Speed that minimizes \(P/U\) (power per unit distance) — the bird travels furthest per joule expended. Critical for migration. Found graphically as the tangent from the origin to the \(P(U)\) curve.
Simulation: Power Curve and Glide Polars
Flight Power Curve (30g Bird) and Glide Polars (Albatross, Hawk, Sparrow)
PythonClick Run to execute the Python code
Code will be executed with Python 3 on the server
Power Curve Anatomy
Module Summary
Kutta–Joukowski
L = ρU∞Γ per unit span; bound circulation enforced by Kutta condition at trailing edge
Lift Equation
L = ½ρU²SC_L; C_L from angle of attack, wing shape, and Re
Drag Polar
C_D = C_D0 + C_L²/(πeAR); induced drag dominates at low speed, parasite at high speed
Aspect Ratio
High AR (albatross AR=18) → low induced drag → efficient soaring; low AR → maneuverability
Glide Ratio
L/D = C_L/C_D; max L/D ≈ ½√(πeAR/C_D0); albatross achieves L/D ≈ 60
U-Shaped Power
P_total = P_ind + P_pro + P_par; minimum power speed U_mp; max range speed U_mr
Hummingbird Hover
Figure-of-8 kinematics; lift on upstroke AND downstroke; wake capture at stroke reversal
Vortex Wake
Continuous sheet (gliding) vs ring vortex model (flapping); wake momentum = 2× lift impulse
References
- Pennycuick, C. J. (2008). Modelling the Flying Bird. Academic Press.
- Norberg, U. M. (1990). Vertebrate Flight. Springer.
- Tobalske, B. W. (2007). Biomechanics of bird flight. Journal of Experimental Biology, 210, 3135–3146.
- Videler, J. J. (2005). Avian Flight. Oxford University Press.
- Anderson, J. D. (2011). Fundamentals of Aerodynamics, 5th ed. McGraw-Hill.
- Hedenström, A. & Alerstam, T. (1995). Optimal flight speed of birds. Philosophical Transactions of the Royal Society B, 348, 471–487.
- Tucker, V. A. (1973). Bird metabolism during flight: evaluation of a theory. Journal of Experimental Biology, 58, 689–709.
- Gill, F. B. (2007). Ornithology, 3rd ed. W. H. Freeman.