Fundamentals of aerodynamics

Aeronautical / Aerospace Engineering (B.E./B.Tech) · Aerodynamics · free study notes for exams, placements and interviews

In short

Aerodynamics is the study of air (and other gases) in motion around bodies, and of the forces and moments the flow produces on those bodies. It rests on three conservation laws - mass, momentum and energy - applied to a fluid treated as a continuum, plus dimensionless groups (Reynolds number, Mach number) that decide which physical effects matter. Every later aerodynamics topic - airfoil theory, compressible flow, boundary layers, CFD - is built on these fundamentals, so they carry heavy weight in GATE Aerospace (Section: Fluid Mechanics and Aerodynamics), ISRO/DRDO/BARC scientist exams, and the Engineering Services route. Design and test roles at HAL, ISRO (VSSC/LPSC), DRDO (ADA, ADE, DRDL), NAL, Boeing/Airbus India engineering centres and drone companies assume this vocabulary from day one.

Explained

What aerodynamics studies, and the continuum assumption

Aerodynamics - the branch of fluid mechanics dealing with gas flow over bodies and the resulting forces, moments, heat transfer and pressure distribution. Divided into external (aircraft, missiles, cars) and internal (intakes, ducts, nozzles) aerodynamics.
Continuum assumption: air is treated as a continuous medium with properties (p, ρ, T, V) defined at every point, not as discrete molecules. Valid when the mean free path λ is far smaller than the body length L.
Knudsen number Kn = λ/L classifies the regime: Kn < 0.01 continuum (classical aerodynamics); 0.01-0.1 slip flow; 0.1-10 transitional; Kn > 10 free-molecular. Re-entry capsules (Gaganyaan crew module) pass through all four on descent; at sea level λ ≈ 68 nm, so continuum holds comfortably.
Fluid properties: density ρ (kg/m³), static pressure p (Pa), temperature T (K), dynamic viscosity μ (Pa·s or N·s/m²), kinematic viscosity ν = μ/ρ (m²/s), speed of sound a (m/s). ISA sea level: p = 101,325 Pa, T = 288.15 K, ρ = 1.225 kg/m³, μ = 1.789 × 10⁻⁵ Pa·s, a = 340.3 m/s.
Viscosity: Newton's law of viscosity τ = μ (du/dy). Air is Newtonian. Gas viscosity rises with temperature (Sutherland's law) - the opposite of liquids.
No-slip condition: at a solid wall the fluid velocity equals the wall velocity. This single condition creates the boundary layer, skin friction and (via separation) most of pressure drag.

Aerodynamic forces, moments and coefficients

All aerodynamic force on a body arises from only two sources distributed over its surface: pressure acting normal to the surface and shear stress acting tangential to it. Nothing else.
The resultant is resolved two ways: along the flight path into drag D (parallel to free stream) and lift L (perpendicular to free stream); or along the body into normal force N and axial force A. With angle of attack α: L = N cos α − A sin α, D = N sin α + A cos α.
Pitching moment M is taken about a stated reference point (leading edge, quarter chord, or centre of gravity) - the point must always be named, since the value changes with it.
Centre of pressure (cp): point where the resultant force acts and the moment is zero; x_cp = −M_LE/N. Aerodynamic centre (ac): point about which the moment is independent of α - near the quarter chord for subsonic airfoils, moving towards mid-chord supersonically.
Non-dimensional coefficients (reference area S, chord or length l, dynamic pressure q∞ = ½ρ∞V∞²): C_L = L/(q∞S), C_D = D/(q∞S), C_M = M/(q∞S l), C_p = (p − p∞)/q∞, skin friction coefficient c_f = τ_w/q∞.
Why coefficients: they collapse all geometrically similar cases onto one curve - a wind-tunnel model and a full-scale aircraft at the same Re and M share the same C_L and C_D. Reference area must be stated (wing planform area for aircraft, frontal area for cars and missiles).
At a stagnation point in incompressible flow C_p = 1 (its maximum); C_p = 0 where p = p∞; C_p is negative wherever the flow is faster than free stream. The theoretical minimum in incompressible flow is unbounded; in compressible flow C_p at stagnation exceeds 1.

Flow classification and the key dimensionless numbers

By viscosity: viscous (real) versus inviscid (ideal). Inviscid theory predicts lift well but gives zero drag for a closed body - d'Alembert's paradox, resolved by viscosity and separation.
By compressibility: incompressible (ρ constant, acceptable for M < 0.3 where density change < 5%) versus compressible.
By speed regime: subsonic M < 0.8; transonic 0.8-1.2 (mixed subsonic/supersonic pockets, shocks form); supersonic 1.2-5; hypersonic M > 5 (real-gas, high-temperature effects).
By rotationality: irrotational (vorticity ω = ∇ × V = 0) versus rotational. Flow outside boundary layers and shocks is usually irrotational.
By time dependence: steady (∂/∂t = 0) versus unsteady; and 1-D, 2-D, axisymmetric or 3-D.
Reynolds number Re = ρVl/μ = inertia/viscous forces - governs boundary layer behaviour, transition and separation. Flat plate transition near Re_x ≈ 5 × 10⁵; a transport aircraft cruises at Re ≈ 10⁷-10⁸, a small drone at 10⁴-10⁵ (hence the low-Reynolds-number airfoil problem).
Mach number M = V/a, a = √(γRT) - ratio of flow speed to sound speed, equivalently √(inertia/elastic forces); governs compressibility.
Others: Prandtl number Pr = μc_p/k ≈ 0.71 for air; Froude number (free-surface and water-channel work); Strouhal number St = fl/V for vortex shedding and flutter; Euler number and pressure coefficient for pressure scaling.

Flow description: Lagrangian, Eulerian and the flow lines

Lagrangian description follows an individual fluid particle; Eulerian description records properties at fixed points in space. Aerodynamics almost always uses the Eulerian form.
Substantial (total) derivative: D/Dt = ∂/∂t + (V·∇) - the rate of change experienced by a moving particle, equal to the local (unsteady) plus convective rates of change. Central to deriving the momentum equation.
Streamline: curve everywhere tangent to the velocity vector at a given instant (dy/dx = v/u in 2-D). Pathline: actual path traced by one particle. Streakline: locus of particles released earlier from one point (what smoke in a tunnel shows). All three coincide only in steady flow.
Stream tube: bundle of streamlines; no flow crosses its walls, so mass flow through it is constant.
Vorticity ω = ∇ × V = 2 × angular velocity of a fluid element. Circulation Γ = −∮ V·ds around a closed curve = −∬ ω·dS (Stokes' theorem). Lift and circulation are linked by the Kutta-Joukowski theorem L′ = ρ∞V∞Γ per unit span.
Stream function ψ (2-D incompressible): u = ∂ψ/∂y, v = −∂ψ/∂x; ψ = constant along a streamline, and ψ₂ − ψ₁ equals the volume flow between the streamlines. Velocity potential φ (irrotational flow): V = ∇φ. Both satisfy Laplace's equation ∇²ψ = 0, ∇²φ = 0, which is linear - so elementary solutions superpose.
Elementary potential flows and their superposition (examinable): uniform flow, source/sink (strength Λ), doublet (κ), point vortex (Γ). Uniform + doublet = non-lifting cylinder; add a vortex = lifting cylinder, giving L′ = ρV∞Γ and illustrating the Magnus effect. Non-lifting cylinder surface pressure: C_p = 1 − 4 sin²θ, with maximum surface speed 2V∞ at θ = 90°.

Governing equations: continuity, momentum, energy

Control volume approach: apply conservation laws to a fixed or moving finite region, then shrink to a point for the differential form. Integral form answers force and thrust problems; differential form is what CFD solves.
Continuity (mass): integral ∂/∂t ∭ρ dV + ∬ρV·dS = 0; differential ∂ρ/∂t + ∇·(ρV) = 0; incompressible steady form ∇·V = 0 and A₁V₁ = A₂V₂ for a stream tube.
Momentum (Newton's second law): the net surface plus body force on a control volume equals the net rate of momentum outflow. Differential form with viscous terms gives the Navier-Stokes equations; dropping viscous terms gives the Euler equations; adding steady, inviscid, incompressible flow along a streamline integrates to Bernoulli's equation.
Energy (first law): for adiabatic inviscid flow, total enthalpy is constant: h + V²/2 = h₀, i.e. c_pT + V²/2 = c_pT₀. Needed only for compressible problems; for incompressible flow continuity and momentum suffice.
Closure: the perfect gas equation of state p = ρRT, with R = 287 J/kg·K and γ = 1.4 for air, closes the system.
Boundary conditions: no-slip (viscous wall) or flow tangency V·n = 0 (inviscid wall); specified free-stream conditions far away; adiabatic or fixed-temperature wall for heat transfer.

Viscous effects: boundary layer, separation and drag breakdown

Boundary layer (Prandtl, 1904): a thin layer near the surface where viscosity matters and velocity rises from 0 to ~0.99 V_e; outside it the flow is effectively inviscid. This split made aerodynamics tractable.
Laminar versus turbulent: laminar layers are thinner with lower skin friction but separate early; turbulent layers have higher skin friction but resist separation because of better momentum mixing.
Thicknesses: boundary layer thickness δ, displacement thickness δ* (the distance the wall would be shifted to give the same mass deficit in inviscid flow), momentum thickness θ (momentum deficit, directly related to drag).
Separation: an adverse pressure gradient (dp/dx > 0) decelerates the near-wall fluid until the wall shear reverses; the flow leaves the surface, creating a wake, high pressure drag and loss of lift (stall).
Drag breakdown: skin friction drag + pressure (form) drag = profile drag; plus induced drag (from wing-tip vortices, 3-D only) and wave drag (supersonic/transonic). Streamlining shrinks pressure drag at small cost in friction drag.
A golf-ball dimple or turbulator trips the boundary layer turbulent, delaying separation and cutting total drag despite more friction - the classic exam counter-intuitive case.

Similarity, testing and Indian practice

Dimensional analysis / Buckingham π theorem: with n variables and k independent dimensions, there are n − k dimensionless groups. For aerodynamic force F(ρ, V, l, μ, a) → C_F = f(Re, M).
Dynamic similarity requires geometric similarity plus matched Re and M (and Pr, γ for heated/hypersonic flows). Matching both at model scale is usually impossible, so pressurised, cryogenic or variable-density tunnels are used, or one parameter is prioritised.
Indian facilities: NAL Bengaluru (1.2 m trisonic, 1.5 m low-speed), the National Wind Tunnel Facility at IIT Kanpur (3 m × 2.25 m low-speed), VSSC Thiruvananthapuram (hypersonic and shock tunnels), IISc Bengaluru shock tunnels (HST series), ADA/HAL for LCA Tejas development testing.
Measurement: pitot-static tubes and multi-tube manometers for pressure; pressure-sensitive scanners and taps for C_p distributions; internal/external strain-gauge balances for forces; hot-wire anemometry, PIV and oil-flow or smoke visualisation for flow fields.
Regulation and standards: civil airworthiness in India is under the DGCA (Aircraft Act, 1934 and the Aircraft Rules; the Bharatiya Vayuyan Adhiniyam, 2024 replaced the 1934 Act, with rules transitioning - check the current notification before citing a section), with CAR Section 2 and the Indian equivalents of FAR/CS-23 and 25 setting flight-performance requirements. Drones follow the Drone Rules, 2021 as amended, with airspace zones on the Digital Sky platform. ISA atmosphere properties follow ISO 2533; aerospace terminology and test practice are covered by BIS aerospace-sector standards and ARDB/CEMILAC documents for military programmes.

Key points

Aerodynamics - gas flow over bodies and the forces and moments produced; external and internal branches.
Continuum assumption - valid when Kn = λ/L < 0.01; beyond that, rarefied gas dynamics is needed.
Source of all aerodynamic force - surface pressure distribution plus surface shear stress, nothing more.
No-slip condition - the physical origin of the boundary layer, skin friction and separation.
Reynolds number Re = ρVl/μ - inertia to viscous forces; decides transition, separation and viscous drag.
Mach number M = V/a - decides compressibility; incompressible treatment valid up to M ≈ 0.3.
Dynamic pressure q = ½ρV² - the natural scale for all aerodynamic loads (not a real static pressure).
Pressure coefficient C_p = (p − p∞)/q∞; equals 1 at an incompressible stagnation point.
Dynamic similarity - geometric similarity + equal Re and M gives equal force coefficients.
Substantial derivative D/Dt = ∂/∂t + V·∇ - local plus convective change following a particle.
Streamline, pathline, streakline - coincide only in steady flow.
Circulation and lift - Kutta-Joukowski L′ = ρ∞V∞Γ links circulation to lift per unit span.
Stream function and velocity potential - both satisfy Laplace's equation; elementary flows superpose.
Euler versus Navier-Stokes - Euler drops viscous terms; Navier-Stokes retains them and needs no-slip.
d'Alembert's paradox - inviscid theory gives zero drag on a closed body; viscosity and separation resolve it.
Boundary layer types - laminar: less friction, early separation; turbulent: more friction, delayed separation.
Drag components - skin friction + form (= profile), plus induced (3-D) and wave (transonic/supersonic).
Aerodynamic centre - moment independent of α; near quarter chord subsonic, near mid-chord supersonic.

Formulas and examples

Dynamic pressure q∞ = ½ ρ∞ V∞² [ρ in kg/m³, V in m/s, q in Pa]
C_L = L/(q∞S), C_D = D/(q∞S), C_M = M/(q∞S l) [L, D in N; M in N·m; S in m²; l = chord or reference length, m]
C_p = (p − p∞)/q∞ ; incompressible C_p = 1 − (V/V∞)²
Re = ρVl/μ = Vl/ν [μ in Pa·s, ν in m²/s; dimensionless]
M = V/a, a = √(γRT) [γ = 1.4, R = 287 J/kg·K for air, T in K]
Kn = λ/L [λ = molecular mean free path, m]
Continuity: ∂ρ/∂t + ∇·(ρV) = 0 ; steady incompressible: A₁V₁ = A₂V₂ [A in m², V in m/s]
Substantial derivative: D/Dt = ∂/∂t + u ∂/∂x + v ∂/∂y + w ∂/∂z
Vorticity ω = ∇ × V [s⁻¹] ; Circulation Γ = −∮ V·ds [m²/s]
Kutta-Joukowski: L′ = ρ∞ V∞ Γ [N per metre of span]
Newton's viscosity law: τ = μ (du/dy) [τ in N/m²] ; c_f = τ_w/q∞
Energy (adiabatic): c_p T + V²/2 = c_p T₀ [c_p = 1005 J/kg·K for air]
Non-lifting cylinder: V_θ = −2V∞ sin θ, C_p = 1 − 4 sin²θ
Resolution: L = N cos α − A sin α ; D = N sin α + A cos α ; x_cp = −M_LE/N

Example 1 - Dynamic pressure, lift coefficient and wing loading. A trainer aircraft of mass 1,400 kg flies straight and level at 180 km/h at sea level (ρ = 1.225 kg/m³). Wing area S = 17 m².

V = 180/3.6 = 50 m/s.
q∞ = 0.5 × 1.225 × 50² = 0.5 × 1.225 × 2,500 = 1,531.25 Pa.
Lift required L = W = 1,400 × 9.81 = 13,734 N.
C_L = 13,734/(1,531.25 × 17) = 13,734/26,031.25 = 0.528.
If C_D = 0.045 at this C_L, D = 0.045 × 26,031.25 = 1,171 N, so L/D = 13,734/1,171 ≈ 11.7, and thrust needed = 1,171 N → power = 1,171 × 50 = 58.6 kW.

Example 2 - Reynolds number, Mach number and regime. Same aircraft, mean chord l = 1.6 m, μ = 1.789 × 10⁻⁵ Pa·s, a = 340.3 m/s.

Re = (1.225 × 50 × 1.6)/(1.789 × 10⁻⁵) = 98/(1.789 × 10⁻⁵) ≈ 5.48 × 10⁶ - fully turbulent over most of the chord.
M = 50/340.3 = 0.147 (< 0.3) → incompressible treatment is justified.
At 11 km (T = 216.65 K): a = √(1.4 × 287 × 216.65) = √(87,041) ≈ 295 m/s. A jet at 250 m/s there has M = 250/295 = 0.847 - transonic, so compressibility and wave drag cannot be ignored.

Example 3 - Pressure coefficient and local speed. A wind tunnel runs at V∞ = 60 m/s, ρ = 1.225 kg/m³, p∞ = 101,325 Pa. A tap on an airfoil reads p = 99,000 Pa.

q∞ = 0.5 × 1.225 × 3,600 = 2,205 Pa.
C_p = (99,000 − 101,325)/2,205 = −2,325/2,205 = −1.054 (suction).
Local speed from C_p = 1 − (V/V∞)²: (V/V∞)² = 1 + 1.054 = 2.054 → V/V∞ = 1.433 → V = 86.0 m/s.
At the stagnation point p = p∞ + q∞ = 103,530 Pa with C_p = 1 and V = 0.

Example 4 - Wind-tunnel similarity (Reynolds scaling). A 1:8 scale model of the above wing (chord 0.2 m) is tested in a sea-level tunnel. What speed matches full-scale Re = 5.48 × 10⁶?

Re equality: ρV_m l_m/μ = 5.48 × 10⁶ → V_m = (5.48 × 10⁶ × 1.789 × 10⁻⁵)/(1.225 × 0.2) = 98.04/0.245 ≈ 400 m/s.
But M = 400/340.3 = 1.18 - supersonic. Re and M cannot both be matched at atmospheric conditions.
Fix: raise tunnel density. At 5 atm (ρ = 6.125 kg/m³), V_m = 98.04/(6.125 × 0.2) = 80 m/s, M = 0.235 - both parameters acceptable. This is why pressurised and cryogenic tunnels exist.

Example 5 - Circulation and lift per unit span. A 2-D airfoil section in air (ρ = 1.225 kg/m³) at V∞ = 70 m/s has circulation Γ = 25 m²/s, chord 1.2 m.

L′ = ρ∞V∞Γ = 1.225 × 70 × 25 = 2,143.75 N/m.
q∞ = 0.5 × 1.225 × 4,900 = 3,001.25 Pa; c_l = 2,143.75/(3,001.25 × 1.2) = 2,143.75/3,601.5 = 0.595.
Check against thin-airfoil theory c_l = 2πα: α = 0.595/6.283 = 0.0947 rad ≈ 5.4°.

Quantity	Symbol	Physical ratio	Decides
Reynolds number	Re	Inertia / viscous	Transition, separation, friction drag
Mach number	M	Inertia / elastic (compressibility)	Shocks, wave drag, density change
Knudsen number	Kn	Mean free path / body size	Continuum validity
Strouhal number	St	Unsteady / convective	Vortex shedding, buffet, flutter
Prandtl number	Pr	Momentum / thermal diffusivity	Thermal boundary layer thickness

Speed regime	Mach range	Dominant physics	Indian example
Low subsonic	M < 0.3	Incompressible, viscous	HAL Dornier 228, agri drones
High subsonic	0.3-0.8	Compressibility corrections	Airbus A320neo on Indian routes
Transonic	0.8-1.2	Local shocks, wave drag rise	LCA Tejas in manoeuvre
Supersonic	1.2-5	Oblique/normal shocks, expansion fans	BrahMos cruise missile
Hypersonic	M > 5	Real-gas, high-temperature, ablation	RLV-TD, HSTDV, Gaganyaan re-entry

Where it is used

Aircraft and helicopter design: HAL (Tejas, HTT-40, LUH, Dhruv), ADA and NAL use these fundamentals for load estimation, performance, drag build-up and CFD setup; aerodynamics and loads engineer roles.
Space launch and re-entry: ISRO's VSSC computes launch-vehicle aerodynamic loads, fin sizing, base drag and aerodynamic heating for PSLV, GSLV/LVM3, SSLV and Gaganyaan; shock-tunnel and CFD validation teams.
Missiles and UAVs: DRDL, RCI and ADE size fins and control surfaces using force coefficients, stability derivatives and centre-of-pressure travel; drone startups (ideaForge, Garuda, NewSpace) work in the low-Re regime.
Wind engineering and civil structures: wind loads on buildings, bridges, towers and solar trackers per IS 875 (Part 3) - wind load code; boundary-layer wind-tunnel studies at CSIR-SERC Chennai and IITs.
Automotive and rail: drag coefficient and cooling-flow work at Maruti Suzuki, Tata Motors, Mahindra, Ashok Leyland and ICF/BEML (Vande Bharat aerodynamic nose), using the same C_D and separation concepts.
Energy and process industry: wind-turbine blade aerodynamics (Suzlon, NIWE), cooling towers, fans and blowers; CFD analyst roles across Infosys, TCS, L&T Technology Services, Cyient and Altair/Ansys India.

Interview questions

Q: What are the only two physical sources of aerodynamic force on a body? A: Pressure distribution acting normal to the surface and shear stress acting tangential to it; lift, drag and moment are just resolved integrals of these two.
Q: Why is air treated as incompressible below Mach 0.3? A: Below M ≈ 0.3 the density change from isentropic relations is under about 5%, so treating ρ as constant introduces negligible error and allows Bernoulli-based methods.
Q: What does the Reynolds number physically represent, and what does it decide? A: The ratio of inertial to viscous forces; it decides boundary-layer transition, separation behaviour and the magnitude of skin-friction drag.
Q: State d'Alembert's paradox and its resolution. A: Inviscid, irrotational flow over a closed body predicts zero drag; in reality viscosity creates skin friction and causes separation, which produces pressure drag.
Q: Difference between a streamline, a pathline and a streakline? A: A streamline is tangent to the velocity field at one instant, a pathline is one particle's actual trajectory, and a streakline joins particles released successively from one point; they coincide only in steady flow.
Q: Why can a wind-tunnel test rarely match both Reynolds and Mach number? A: At model scale, matching Re needs a much higher speed, which pushes Mach number too high; so tunnels are pressurised or cryogenic, or one parameter is matched and the other corrected.
Q: What is the difference between the centre of pressure and the aerodynamic centre? A: The centre of pressure is where the resultant force acts (zero moment) and it moves with angle of attack; the aerodynamic centre is the fixed point about which the moment coefficient does not change with angle of attack.
Q: What is the substantial derivative and why does it matter? A: D/Dt = ∂/∂t + V·∇ gives the rate of change following a fluid particle, combining local unsteady and convective change; it converts Newton's law for a particle into the Eulerian momentum equation.
Q: Why does a turbulent boundary layer reduce total drag on a bluff body? A: Turbulent mixing carries momentum to the wall, letting the layer resist the adverse pressure gradient longer; separation is delayed, the wake shrinks and the drop in pressure drag outweighs the rise in friction drag.
Q: How are circulation and lift related? A: By the Kutta-Joukowski theorem, lift per unit span L′ = ρ∞V∞Γ - lift exists only if there is net circulation around the section.
Q: Can the pressure coefficient exceed 1? A: In incompressible flow its maximum is exactly 1, at a stagnation point; in compressible flow the stagnation value exceeds 1 and grows with Mach number. Negative values (suction) have no lower bound in incompressible theory.
Q: What assumptions turn the Navier-Stokes equations into Bernoulli's equation? A: Inviscid, steady, incompressible flow with no body-force work, integrated along a single streamline (or everywhere if the flow is also irrotational).

Common mistakes

Treating dynamic pressure ½ρV² as a measurable static pressure, or using sea-level ρ for an altitude problem - always use the local ρ from ISA.
Mixing units: km/h with m/s, mm of mercury with Pa, cm² with m² - convert first, then compute coefficients.
Quoting C_L or C_D without stating the reference area, or comparing a car's C_D (frontal area) with an aircraft's (planform area).
Reporting a pitching moment or centre of pressure without naming the reference point.
Assuming inviscid theory means no drag in practice, or using potential-flow results in a separated region.
Confusing vorticity (a local, point property, = 2ω) with circulation (an integral around a closed curve), and forgetting the sign convention in Γ = −∮V·ds.
Applying Bernoulli across a boundary layer, a shock, a fan or a region with losses - it holds only along a streamline in steady, inviscid, incompressible flow.