2.5 Lifting-line theory

Lifting-line theory fundamentals

Lifting-line theory gives you a mathematical framework for predicting how finite (real, 3D) wings generate lift and drag. While 2D airfoil theory tells you what happens at a single cross-section, it can't account for the spanwise effects that come with a wing that actually ends at two tips. Lifting-line theory fills that gap, and it remains one of the most important analytical tools in wing design.

Ludwig Prandtl developed this theory in 1918, and it still forms the backbone of how engineers think about finite wing performance.

Prandtl's lifting-line concept

Prandtl's key idea is to replace the entire wing with a single bound vortex filament running along the quarter-chord line. This filament is assumed straight and perpendicular to the freestream velocity.

The vortex strength isn't constant across the span. It varies from root to tip to capture how lift changes along the wing. Where the circulation changes, vorticity is shed downstream, forming trailing vortices. Together, the bound vortex and its trailing vortices create what's called a horseshoe vortex system.

Vortex filament representation

The vortex filament is an idealization of the circulation around the wing. Circulation ($\Gamma$) measures the vortex strength and directly relates to the lift the wing produces.

This filament induces a velocity field in the surrounding flow. The induced velocity acts perpendicular to the filament, and its magnitude scales with the local circulation strength. That induced velocity is what modifies the flow each wing section "sees."

Bound vortex vs trailing vortices

These two types of vortices play different roles:

• The bound vortex represents circulation around the wing itself and is responsible for generating lift.
• Trailing vortices shed from the wing tips and extend downstream, forming the vortex wake.

Trailing vortices exist because the wing has finite span. High pressure on the lower surface and low pressure on the upper surface drive air around the tips, rolling it into vortices. The strength of the trailing vorticity at any spanwise station equals the local rate of change of circulation, $d\Gamma/dy$.

Helmholtz's vortex theorems

Helmholtz's theorems govern vortex behavior in inviscid, incompressible flow and underpin the entire horseshoe vortex model:

1. First theorem: The strength of a vortex filament is constant along its entire length.
2. Second theorem: A vortex filament cannot end in the fluid. It must either form a closed loop or extend to infinity (or to a boundary).

These theorems explain why the bound vortex must trail vortices downstream. The bound vortex can't just stop at the wing tips; the circulation has to continue into the wake.

Kutta-Joukowski theorem

The Kutta-Joukowski theorem connects circulation to lift. The lift per unit span is:

$L' = \rho_{\infty} V_{\infty} \Gamma$

where $\rho_{\infty}$ is the freestream density, $V_{\infty}$ is the freestream velocity, and $\Gamma$ is the local circulation.

The theorem assumes the flow leaves the trailing edge smoothly (the Kutta condition). This relationship is what lets you go from a circulation distribution $\Gamma(y)$ to a complete lift distribution across the span.

Circulation distribution

The shape of the circulation distribution along the span determines both the lift distribution and the induced drag. Different distributions produce very different aerodynamic behavior.

Elliptical vs non-elliptical distributions

An elliptical circulation distribution is the theoretical optimum: it minimizes induced drag for a given total lift and wingspan. The reason is that elliptical loading produces a uniform downwash velocity across the entire span.

Non-elliptical distributions (triangular, rectangular, or anything else) produce non-uniform downwash and therefore higher induced drag. However, non-elliptical distributions aren't always bad. For example, a wing designed to stall at the root first (for better roll control) will intentionally deviate from elliptical loading.

Relation to lift distribution

Lift distribution is directly proportional to circulation distribution through the Kutta-Joukowski theorem. Elliptical $\Gamma(y)$ gives elliptical $L'(y)$, and so on.

The lift distribution determines the spanwise loading, which in turn drives structural requirements (bending moments) and overall aerodynamic efficiency. Getting the lift distribution right is one of the central goals of wing design.

Effect on induced drag

Induced drag arises because trailing vortices create downwash, which tilts the local lift vector backward, producing a drag component. For the ideal elliptical case, the induced drag coefficient is:

$C_{D,i} = \frac{C_L^2}{\pi AR}$

This is the minimum possible induced drag for a given $C_L$ and aspect ratio $AR$. Any non-elliptical distribution will have higher induced drag. To account for this, a span efficiency factor $e$ (where $e \leq 1$) is introduced:

$C_{D,i} = \frac{C_L^2}{\pi e AR}$

For a perfect elliptical distribution, $e = 1$. Real wings typically have $e$ in the range of 0.85 to 0.95.

Induced angle of attack

Trailing vortices don't just create drag directly. They also change the angle of attack that each wing section experiences.

Downwash velocity

The downwash velocity ($w$) is the downward component of velocity induced by the trailing vortex sheet. It reduces the angle at which air actually meets each wing section.

For an elliptical distribution, the downwash is constant across the span. For non-elliptical distributions, it varies, with the distribution shape determining where downwash is strongest. The downwash is related to the circulation distribution and the geometry of the trailing vortex system.

Effective angle of attack

The wing is set at some geometric angle of attack $\alpha$ relative to the freestream. But because of downwash, each section actually "sees" a reduced angle called the effective angle of attack:

$\alpha_{eff} = \alpha - \alpha_i$

where the induced angle of attack is:

$\alpha_i = \frac{w}{V_{\infty}}$

This is a critical concept. The local lift coefficient depends on $\alpha_{eff}$, not $\alpha$. So downwash effectively reduces the lift a finite wing produces compared to what 2D theory would predict at the same geometric angle.

Induced drag coefficient

You can also understand induced drag geometrically. The downwash tilts the local aerodynamic force vector backward by the angle $\alpha_i$. The rearward component of the tilted lift vector is the induced drag.

The induced drag coefficient for the general case is:

$C_{D,i} = \frac{C_L^2}{\pi e AR}$

Two key takeaways:

• Induced drag grows with the square of $C_L$, so it's most significant at high lift (low speed, heavy aircraft, steep climbs).
• Higher aspect ratio wings have lower induced drag at any given $C_L$, which is why gliders and long-range aircraft use high-AR wings.

Finite wing theory

Finite wing theory takes the lifting-line concepts and turns them into equations you can actually solve for the circulation distribution and aerodynamic forces.

Lifting-line equation derivation

The derivation applies the Biot-Savart law to the entire horseshoe vortex system. The Biot-Savart law gives the velocity induced at any point by a vortex filament of known strength and geometry.

By requiring that the flow at each spanwise station be consistent with the local airfoil properties (i.e., the local lift coefficient matches what the effective angle of attack demands), you get the fundamental equation of lifting-line theory. This is an integro-differential equation relating $\Gamma(y)$ to the wing geometry, twist distribution, and airfoil characteristics.

Fourier series solution

A standard approach is to express the circulation distribution as a Fourier sine series. The unknown coefficients $A_n$ are found by enforcing the lifting-line equation at a set of spanwise stations, which converts the integral equation into a system of linear algebraic equations.

Once you have the $A_n$ coefficients, you can compute lift, induced drag, and the span efficiency factor directly from the series.

Glauert's trigonometric series

Glauert formalized the Fourier approach using the coordinate transformation $y = -\frac{b}{2}\cos\theta$, where $\theta$ runs from $0$ to $\pi$ across the span. The circulation becomes:

$\Gamma(\theta) = 4bV_{\infty}\sum_{n=1}^{N}A_n\sin(n\theta)$

The first coefficient $A_1$ dominates and corresponds to the elliptical loading component. Higher-order coefficients ($A_2, A_3, \ldots$) represent deviations from elliptical loading and always increase induced drag.

This is the method you'll most commonly see in textbooks and homework problems. The steps are:

1. Apply the coordinate transformation to convert spanwise position $y$ to $\theta$.
2. Write the lifting-line equation at $N$ spanwise stations (values of $\theta$).
3. Solve the resulting $N \times N$ linear system for the coefficients $A_1, A_2, \ldots, A_N$.
4. Compute lift from $A_1$ and induced drag from all $A_n$.

Multhopp's quadrature method

Multhopp's method is a numerical technique that discretizes the span into control points and uses a quadrature rule to evaluate the downwash integrals. It reduces the lifting-line equation to a system of linear algebraic equations, similar to the Fourier approach but without assuming a specific series form.

The advantage is versatility: Multhopp's method handles arbitrary planforms, twist distributions, and even mixed airfoil sections along the span without requiring analytical simplifications.

Wing planform effects

The wing's planform shape has a major influence on its aerodynamic characteristics. Lifting-line theory lets you study how aspect ratio, taper, and sweep each affect lift distribution and induced drag.

Aspect ratio influence

Aspect ratio is defined as:

$AR = \frac{b^2}{S}$

where $b$ is the wingspan and $S$ is the wing area. Higher AR means a longer, narrower wing.

High aspect ratio reduces induced drag because the trailing vortices are farther apart, weakening the downwash over most of the span. This is why sailplanes have aspect ratios of 20+ while fighter jets may have AR below 4. The tradeoff is structural: long, slender wings experience higher bending loads and are heavier.

Taper ratio considerations

Taper ratio is:

$\lambda = \frac{c_{tip}}{c_{root}}$

A rectangular wing has $\lambda = 1$; a highly tapered wing might have $\lambda = 0.2$.

Taper shifts the lift distribution inboard, which reduces the wing root bending moment (a structural benefit). A taper ratio around $\lambda \approx 0.4$ produces a lift distribution close to elliptical for an untwisted wing, which is why moderate taper ratios (0.3 to 0.5) are common in practice. Very low taper ratios push the stall toward the tips, which is undesirable because it reduces aileron effectiveness.

Sweep angle impact

Sweep angle ($\Lambda$) is measured between the leading edge (or quarter-chord line) and a line perpendicular to the fuselage.

Sweep is primarily used to delay compressibility drag rise at high subsonic Mach numbers. The effective velocity component perpendicular to the leading edge is $V_{\infty}\cos\Lambda$, so sweep reduces the effective Mach number the airfoil sections experience.

Sweep also introduces a spanwise flow component, which affects the lift distribution by shifting load outboard on a swept-back wing. This can promote tip stall, so swept wings often use geometric or aerodynamic twist (washout) to compensate.

Winglets and end plates

Winglets are vertical or near-vertical surfaces at the wing tips designed to reduce the strength of trailing vortices. They work by disrupting the tip flow that rolls from the lower surface to the upper surface.

The net effect is similar to increasing the wing's effective aspect ratio without physically extending the span. Modern winglets can reduce induced drag by 3-5%, translating directly into fuel savings on transport aircraft. End plates serve a similar function but are less aerodynamically refined and produce more parasitic drag than well-designed winglets.

Limitations and extensions

Lifting-line theory is remarkably useful, but it rests on assumptions that break down in certain flight regimes and wing configurations.

Assumptions and simplifications

The key assumptions are:

• The wing has a high aspect ratio (typically AR > 4 for reasonable accuracy).
• Angles of attack are small so that linear airfoil theory applies.
• The flow is inviscid and incompressible.
• Wing thickness effects are neglected.
• The wake is planar and extends straight downstream.

These assumptions mean lifting-line theory loses accuracy for low aspect ratio wings, high angles of attack (near stall), transonic or supersonic flow, and wings with significant thickness effects.

Prandtl's classical lifting-line theory

Prandtl's original formulation uses all the assumptions above. It works well for high-AR, unswept wings at low angles of attack, which covers a wide range of practical cases (general aviation, gliders, transport aircraft wings at cruise).

It does not account for sweep, dihedral, or wake deformation. For unswept, high-AR wings, though, it remains one of the best tools for quick, accurate estimates.

Weissinger's extended lifting-line model

Weissinger extended the theory to handle swept and non-planar wings. His model accounts for the non-planar shape of the vortex sheet and places the boundary condition at the three-quarter chord rather than the quarter chord.

This makes it applicable to wings with moderate to high sweep angles where Prandtl's classical theory would give poor results. Weissinger's method is still relatively simple computationally while offering significantly improved accuracy for swept configurations.

Vortex lattice method comparison

The vortex lattice method (VLM) is a natural numerical extension of lifting-line ideas. Instead of a single vortex line, VLM distributes horseshoe vortices over a lattice of panels covering the entire wing surface.

VLM handles complex geometries (sweep, taper, twist, dihedral, multiple surfaces) and captures some effects that lifting-line theory misses. It's more computationally expensive but still far cheaper than CFD.