2.3 Thin airfoil theory

Fundamentals of thin airfoil theory

Thin airfoil theory simplifies airfoil analysis by replacing the actual airfoil shape with a vortex sheet placed along the camber line. This lets you calculate lift, moment, and pressure distribution using analytical methods rather than numerical simulation. The trade-off is a set of restrictive assumptions, but for preliminary design and building intuition about how airfoils behave, it's one of the most useful tools in aerodynamics.

The theory rests on potential flow theory and assumes inviscid, incompressible, irrotational flow around a thin airfoil at small angles of attack.

Assumptions and limitations

• The airfoil is thin: both camber and thickness are small relative to the chord length.
• The flow is inviscid (no viscosity), incompressible (constant density), and irrotational (no vorticity in the freestream).
• Boundary layer separation and stall are completely neglected.
• Valid only for small angles of attack (typically below about 10°) and subsonic flow.

These assumptions mean the theory won't tell you anything about drag (since drag in this framework is zero by d'Alembert's paradox) or about what happens near or beyond stall. Keep that in mind whenever you apply it.

Airfoil geometry

The shape of an airfoil directly controls its aerodynamic characteristics. Three geometric features matter most: the camber line, the thickness distribution, and the shapes of the leading and trailing edges.

Camber line and chord line

The camber line (also called the mean line) is the curve that runs midway between the upper and lower surfaces of the airfoil. The chord line is the straight line connecting the leading edge to the trailing edge, and the chord length $c$ is the distance between them.

Camber is the maximum perpendicular distance between the camber line and the chord line, usually expressed as a percentage of $c$. A NACA 2412 airfoil, for example, has 2% maximum camber located at 40% chord. Positive camber is what allows a cambered airfoil to produce positive lift even at zero angle of attack.

Thickness distribution

Thickness distribution describes how the airfoil's thickness varies from leading edge to trailing edge. Symmetric airfoils (like the NACA 0012) have zero camber and equal thickness above and below the chord line. Thickness affects structural strength, stall characteristics, and drag, but in thin airfoil theory, thickness is neglected entirely. The airfoil is collapsed down to just its camber line.

Leading and trailing edges

The leading edge is where the oncoming flow first meets the airfoil. Its radius of curvature influences stall behavior and maximum lift coefficient: a blunter leading edge generally delays leading-edge stall.

The trailing edge is where the upper and lower surface flows rejoin. Its geometry is critical because it's where the Kutta condition is enforced, which directly determines the circulation and lift.

Vortex sheet representation

In thin airfoil theory, you replace the physical airfoil with a vortex sheet distributed along the camber line. This sheet consists of infinitesimal point vortices whose strengths vary with position. The total circulation produced by this sheet determines the lift.

Continuous vortex distribution

The vortex strength per unit length is described by a function $\gamma(x)$, where $x$ is the chordwise position. This function varies continuously along the camber line, and its distribution depends on both the camber line shape and the angle of attack.

To make the math cleaner, a change of variables is used. The chordwise coordinate $x$ is related to an angular variable $\theta$ by:

$x = \frac{c}{2}(1 - \cos\theta)$

where $\theta$ ranges from $0$ (leading edge) to $\pi$ (trailing edge). This transformation clusters points near the leading and trailing edges, where the flow changes most rapidly.

Kutta condition at trailing edge

The Kutta condition requires that the flow leaves the trailing edge smoothly with finite velocity. Mathematically, this means:

$\gamma(\theta = \pi) = 0$

Without this condition, the vortex strength distribution would not be unique. Physically, the Kutta condition represents the role that viscosity plays in real flows: it prevents the physically impossible infinite velocity that would otherwise occur at a sharp trailing edge. Even though the theory is inviscid, this condition smuggles in the most important viscous effect.

Derivation of thin airfoil equations

The thin airfoil equations come from enforcing the flow-tangency (no-penetration) boundary condition on the camber line, then solving for the unknown vortex distribution $\gamma(\theta)$.

Governing equations and boundary conditions

Three conditions govern the problem:

1. No-penetration condition: The flow velocity normal to the camber line must be zero at every point on the camber line. This is the fundamental boundary condition.
2. Kutta condition: $\gamma = 0$ at the trailing edge.
3. Far-field condition: Disturbances caused by the airfoil must vanish at infinity (the flow returns to the undisturbed freestream far from the airfoil).

Linearization of boundary conditions

For a thin airfoil at a small angle of attack, the boundary conditions can be linearized. Instead of applying the no-penetration condition on the actual curved camber line, you apply it on the chord line ($y = 0$). The camber line's effect enters as a known function in the boundary condition.

The linearized boundary condition takes the form:

$\frac{1}{2\pi}\int_0^c \frac{\gamma(\xi)}{x - \xi} d\xi = V_\infty\left(\alpha - \frac{dy_c}{dx}\right)$

The left side is the downwash induced by the vortex sheet (computed via the Biot-Savart law), and the right side combines the angle of attack $\alpha$ with the local slope of the camber line $dy_c/dx$. This is a singular integral equation for $\gamma$.

Solution using Fourier series

The standard solution method uses the $\theta$ substitution and expands $\gamma(\theta)$ as a Fourier-type series:

$\gamma(\theta) = 2V_\infty\left[A_0\frac{1 + \cos\theta}{\sin\theta} + \sum_{n=1}^{\infty} A_n \sin(n\theta)\right]$

The first term (with $A_0$) automatically satisfies the Kutta condition and represents the flat-plate contribution. The remaining terms capture the effect of camber. The coefficients are found by substituting into the linearized boundary condition and using orthogonality of trigonometric functions:

$A_0 = \alpha - \frac{1}{\pi}\int_0^{\pi}\frac{dy_c}{dx}d\theta$

$A_n = \frac{2}{\pi}\int_0^{\pi}\frac{dy_c}{dx}\cos(n\theta)\,d\theta$

Lift and moment coefficients

Thin airfoil theory gives closed-form expressions for the lift and moment coefficients in terms of the Fourier coefficients $A_0$, $A_1$, $A_2$, etc.

Expressions for lift and moment

The lift coefficient is:

$C_L = 2\pi\left(\alpha - \alpha_{L=0}\right)$

where the zero-lift angle of attack is:

$\alpha_{L=0} = -\frac{1}{\pi}\int_0^{\pi}\frac{dy_c}{dx}(\cos\theta - 1)\,d\theta$

This can also be written as:

$C_L = \pi(2A_0 + A_1)$

The moment coefficient about the leading edge is:

$C_{M_{LE}} = -\frac{C_L}{4} + \frac{\pi}{4}(A_1 - A_2)$

The moment coefficient about the quarter-chord ($x = c/4$) is:

$C_{M_{c/4}} = \frac{\pi}{4}(A_2 - A_1)$

This quarter-chord moment depends only on camber (through $A_1$ and $A_2$), not on angle of attack. That's why the quarter-chord is the aerodynamic center for a thin airfoil in incompressible flow.

Effect of angle of attack

The angle of attack $\alpha$ is the angle between the chord line and the freestream direction. Key results:

• $C_L$ increases linearly with $\alpha$ at a rate of $2\pi$ per radian (about 0.11 per degree). This is the lift-curve slope, often written $a_0 = 2\pi$.
• A symmetric airfoil (zero camber) produces zero lift at $\alpha = 0$.
• A positively cambered airfoil has a negative $\alpha_{L=0}$, meaning it still produces lift at zero geometric angle of attack.
• The linear relationship breaks down at high angles of attack when the real flow separates (stall), which thin airfoil theory cannot predict.

Thin airfoil approximations

For a symmetric airfoil ($dy_c/dx = 0$ everywhere), all the camber-related Fourier coefficients vanish, and the results simplify to:

• $C_L = 2\pi\alpha$
• $C_{M_{LE}} = -\frac{\pi\alpha}{2} = -\frac{C_L}{4}$
• $C_{M_{c/4}} = 0$

These are the simplest and most commonly tested results from thin airfoil theory.

Pressure distribution on airfoil

The pressure distribution across the airfoil surface reveals where lift is generated and how the aerodynamic loading is distributed.

Calculation using thin airfoil theory

The pressure coefficient at any point is:

$C_p = 1 - \left(\frac{V}{V_\infty}\right)^2$

In thin airfoil theory, the local velocity difference between upper and lower surfaces is directly related to the local vortex strength:

$\Delta V = \gamma(\theta)$

So the pressure difference between lower and upper surfaces (the loading) is proportional to $\gamma(\theta)$. Integrating this loading along the chord gives the total lift and moment, which should match the coefficient expressions derived earlier.

For a symmetric airfoil at angle of attack, the loading is highest near the leading edge and drops to zero at the trailing edge (enforced by the Kutta condition). Camber shifts and reshapes this loading distribution.

Comparison with experimental data

Wind tunnel experiments measure surface pressures using pressure taps or pressure-sensitive paint. Thin airfoil theory agrees well with experimental data for:

• Thin airfoils (thickness ratio below about 12%)
• Small angles of attack (well below stall)
• Regions away from the leading edge

Discrepancies appear near the leading edge (where the theory predicts an inverse-square-root singularity in $\gamma$ that real viscous flows smooth out) and at angles of attack approaching stall, where boundary layer separation becomes significant.

Limitations of thin airfoil theory

Validity for small angles of attack

The linear relationship $C_L = 2\pi\alpha$ holds only for small $\alpha$, typically below about 8-10°. Beyond this range, viscous effects cause boundary layer separation, and the actual lift curve departs from the linear prediction. At stall, lift drops sharply. Thin airfoil theory has no mechanism to predict stall or post-stall behavior.

Inaccuracies near leading edge

The theory predicts that $\gamma \to \infty$ at the leading edge (an integrable singularity). In reality, the leading edge has a finite radius, and viscous effects prevent infinite velocities. This means the predicted pressure distribution near the leading edge is unreliable. For airfoils with very small leading-edge radii, this can lead to poor predictions of suction peaks and the onset of leading-edge separation.

Neglect of viscous effects

Because the theory is inviscid, it predicts zero drag (d'Alembert's paradox). Real airfoils experience both skin friction drag and pressure drag due to boundary layer effects. The theory also cannot account for:

• Boundary layer transition (laminar to turbulent)
• Flow separation and stall
• Reynolds number effects on airfoil performance

These limitations become more pronounced at low Reynolds numbers and near the maximum lift coefficient.

Extensions and modifications

Inclusion of thickness effects

Thin airfoil theory treats the airfoil as having zero thickness. To account for thickness, you can decompose the problem: the camber line generates lift (handled by thin airfoil theory), and the thickness distribution is modeled separately as a source/sink distribution that doesn't produce lift but modifies the velocity field. First-order thickness corrections adjust the predicted pressures and velocities. Higher-order methods, such as van Dyke's second-order theory, improve accuracy for thicker airfoils.

Correction for boundary layer

Viscous-inviscid interaction (VII) methods couple the inviscid thin airfoil solution with a boundary layer calculation. The boundary layer's displacement thickness effectively changes the airfoil shape seen by the outer flow. Steps in a typical VII approach:

1. Solve the inviscid problem (thin airfoil theory) to get the pressure distribution.
2. Use that pressure distribution to compute the boundary layer development.
3. Calculate the displacement thickness and add it to the airfoil shape.
4. Re-solve the inviscid problem with the modified shape.
5. Iterate until convergence.

This approach improves predictions of lift, drag, and stall onset, especially at moderate angles of attack.

Application to cascades and propellers

Thin airfoil theory extends naturally to cascades, which are arrays of airfoils found in compressors and turbines. Cascade theory accounts for the mutual interaction between adjacent blades and introduces the parameter solidity (chord/spacing ratio), which affects the lift and deflection of the flow.

For propellers and rotors, thin airfoil theory combines with blade element momentum theory to predict thrust and torque. The rotating reference frame introduces additional effects (centrifugal and Coriolis forces), and three-dimensional flow corrections become necessary. These extended theories are widely used in the design of axial compressors, turbines, helicopter rotors, and marine propellers.