3.1.21Compressible Flow & Aerodynamics

Thin airfoil theory — lift per unit span = πρV²(α + 2β - π c)

1,782 words8 min readdifficulty · medium4 backlinks

WHAT are we computing?

The prime (LL') means per unit span — units of N/m, not N. We will derive this, never dump it.


WHY does a wing lift at all? (first principles)


HOW we derive it — the vortex sheet

Step 1 — Replace the airfoil with a vortex sheet. Place a sheet of vorticity γ(x)\gamma(x) along the chord (0xc0 \le x \le c). Each element acts like a tiny vortex of strength γdx\gamma\,dx.

Why this step? A thin airfoil barely disturbs the flow, so the camber line itself can carry the vorticity. We don't need the real curved surface — a line of vortices reproduces the same outer flow.

Step 2 — Flow-tangency (boundary) condition. The flow must be tangent to the camber line. The velocity normal to the mean line from (a) freestream and (b) the vortex sheet must cancel: V ⁣(αdzdx)=w(x)V\!\left(\alpha - \frac{dz}{dx}\right) = w(x) where z(x)z(x) is the camber shape and w(x)w(x) is the downwash induced by the sheet.

Why this step? Air cannot flow through the surface. This is the condition that pins down γ(x)\gamma(x).

Step 3 — Downwash from Biot–Savart. A vortex element γ(ξ)dξ\gamma(\xi)d\xi at position ξ\xi induces at xx: w(x)=12π0cγ(ξ)dξxξw(x) = -\frac{1}{2\pi}\int_0^c \frac{\gamma(\xi)\,d\xi}{x-\xi}

Why this step? It is just the 2-D vortex velocity v=Γ/2πrv=\Gamma/2\pi r summed over the sheet.

Step 4 — Substitute the cosine transform. Let x=c2(1cosθ)x = \tfrac{c}{2}(1-\cos\theta). The standard Glauert solution is γ(θ)=2V ⁣(A01+cosθsinθ+n=1Ansinnθ)\gamma(\theta) = 2V\!\left(A_0\frac{1+\cos\theta}{\sin\theta} + \sum_{n=1}^\infty A_n\sin n\theta\right) satisfying the Kutta condition (finite γ\gamma at trailing edge). Plugging into the boundary condition and using Glauert's integral 0πcosnθcosθcosθdθ=πsinnθsinθ\int_0^\pi \frac{\cos n\theta'}{\cos\theta'-\cos\theta}d\theta' = \frac{\pi\sin n\theta}{\sin\theta} gives the Fourier coefficients A0=α1π0πdzdxdθ,An=2π0πdzdxcosnθdθA_0 = \alpha - \frac{1}{\pi}\int_0^\pi \frac{dz}{dx}\,d\theta, \qquad A_n = \frac{2}{\pi}\int_0^\pi \frac{dz}{dx}\cos n\theta\,d\theta

Step 5 — Total circulation and lift. Γ=0cγdx=cVπ ⁣(A0+A12)\Gamma = \int_0^c \gamma\,dx = c\,V\,\pi\!\left(A_0 + \frac{A_1}{2}\right) L=ρVΓ=πρcV2(A0+A12)\boxed{L' = \rho V \Gamma = \pi \rho c V^2\left(A_0 + \frac{A_1}{2}\right)}

Why this step? Only A0A_0 and A1A_1 survive the integration for Γ\Gamma (the orthogonality of sinnθ\sin n\theta kills n2n\ge2). This is the general result.

Step 6 — The simple cambered case. For a symmetric flat plate, dz/dx=0A0=α, A1=0dz/dx=0 \Rightarrow A_0=\alpha,\ A_1=0, giving the famous L=πρcV2αc=2πα.L' = \pi\rho c V^2 \alpha \quad\Rightarrow\quad c_\ell = 2\pi\alpha. For the parabolic camber line whose combined camber contribution equals A0+A1/2=α+2βπA_0+A_1/2 = \alpha + \dfrac{2\beta}{\pi}, we recover the target: L=πρV2c(α+2βπ)\boxed{L' = \pi\rho V^2 c\left(\alpha + \frac{2\beta}{\pi}\right)}

Figure — Thin airfoil theory — lift per unit span = πρV²(α + 2β - π c)

Worked examples



Flashcards

Kutta–Joukowski theorem for 2-D lift
L=ρVΓL' = \rho V \Gamma
What does the prime in LL' mean?
Lift per unit span (N/m), not total lift.
Lift-curve slope from thin airfoil theory
dc/dα=2πdc_\ell/d\alpha = 2\pi per radian.
Section lift coefficient with camber
c=2πα+4βc_\ell = 2\pi\alpha + 4\beta.
What sets the circulation magnitude on a sharp airfoil?
The Kutta condition (smooth flow off the trailing edge).
Which Fourier coefficients survive in the lift integral?
Only A0A_0 and A1A_1: Γ=πcV(A0+A1/2)\Gamma = \pi c V(A_0 + A_1/2).
Zero-lift angle of attack
αL=0=2β/π\alpha_{L=0} = -2\beta/\pi.
Why must α\alpha be in radians?
The small-angle linearisation sinαα\sin\alpha\approx\alpha assumed radians.
What replaces the airfoil in the model?
A vortex sheet γ(x)\gamma(x) along the camber line.
Does camber change the lift-curve slope?
No — it only shifts the curve (constant offset).

Recall Feynman: explain to a 12-year-old

Imagine throwing a Frisbee. To make the air push it up, the wing has to shove air downwards. A flat tilted card pushes air down a little; a card that's gently curved (cambered) scoops even more air down even when it isn't tilted much. We pretend the wing is a row of tiny spinning straws that grab the air and twist it down. Count up all the spinning, multiply by how fast the air rushes and how heavy the air is, and that's your lift. Tilt (α\alpha) and curve (β\beta) are the two ways to make it scoop harder.


Connections

  • Kutta–Joukowski theorem — supplies L=ρVΓL'=\rho V\Gamma.
  • Circulation and bound vortices — origin of Γ\Gamma.
  • Kutta condition — selects the physical γ(x)\gamma(x).
  • Biot–Savart law — gives the induced downwash w(x)w(x).
  • Glauert's integral and Fourier coefficients — solves the sheet.
  • Lift coefficient and angle of attack — the cc_\ellα\alpha curve.
  • Compressibility corrections (Prandtl–Glauert) — next step beyond incompressible.

Concept Map

reduces to finding

forces net

replaces

integrates to

carries vorticity on

pins down

involves

involves

feeds into

solves for

via Kutta-Joukowski gives

adds slice to

adds slice to

Kutta-Joukowski L' = rho V Gamma

Circulation Gamma

Kutta condition

Vortex sheet gamma of x

Thin cambered airfoil

Camber line z of x

Flow-tangency condition

Angle of attack alpha

Camber slope beta

Biot-Savart downwash w of x

Glauert cosine transform

Lift per span L' = pi rho V^2 c times alpha + 2 beta over pi

Hinglish (regional understanding)

Intuition Hinglish mein samjho

Dekho, thin airfoil theory ka core idea bahut simple hai: ek patli wing ko hum maan lete hain ki wo ek line hai jisme bahut saare chhote-chhote "vortices" (ghoomne wali air ke tukde) lage hain. Inka total ghoomna circulation Γ\Gamma kehlata hai. Aur Kutta–Joukowski theorem kehta hai ki lift per unit span L=ρVΓL' = \rho V \Gamma. Matlab poora khel sirf Γ\Gamma nikalne ka hai — baaki sab usi se aata hai.

Ab Γ\Gamma do cheezon se badhta hai: ek to angle of attack α\alpha (wing kitni tilt hai) aur doosra camber β\beta (wing kitni curved hai). Inn dono ko jodke final formula banta hai: L=πρV2c(α+2βπ)L' = \pi \rho V^2 c (\alpha + \tfrac{2\beta}{\pi}). Isse c=2πα+4βc_\ell = 2\pi\alpha + 4\beta milta hai. Yaad rakho — slope hamesha 2π2\pi per radian rehta hai; camber sirf line ko upar shift karta hai, slope ko nahi badalta.

Practical baat: camber ki wajah se wing zero angle pe bhi lift de deti hai (thoda negative angle pe lift zero hota hai, jise αL=0=2β/π\alpha_{L=0} = -2\beta/\pi kehte hain). Isiliye real aeroplane wings flat nahi, gently curved hoti hain — kam tilt pe zyada lift. Ek common galti: α\alpha ko degree me daal dena. Hamesha radian me convert karo, kyunki derivation me sinαα\sin\alpha \approx \alpha radian me hi sahi hota hai.

Go deeper — visual, from zero

Test yourself — Compressible Flow & Aerodynamics

Connections