Introduction

When a freestream flow approaches a finite wing, characterized by a definite span from tip to tip, the downstream flow develops a trailing-vortex wake system, with the strongest vortices concentrated near the wing tips, as shown in Figure 1. These vortices induce a vertical downwash velocity over the wing’s surface, particularly near the wing tips.^[1] This downwash reduces the effective angle of attack of the wing relative to its geometric angle of attack and produces a rearward component of the aerodynamic force, which is called induced drag.

A wing’s aerodynamic performance, including its lift and drag characteristics, is significantly influenced by its geometric design. Factors such as planform shape (chord distribution), span and aspect ratio, and spanwise twist are all essential in determining the wing’s aerodynamic efficiency. The wing design process must balance these elements to optimize lift generation and minimize induced drag, among other factors, to achieve the wing’s best aerodynamic performance for the intended application.

Lifting line theory is another cornerstone of classical aerodynamics, explaining and predicting the aerodynamic behavior of finite wings at low speeds, i.e., without considering compressibility effects. Unlike two-dimensional airfoil theory, which assumes infinite span and aspect ratio, lifting-line theory accounts for the three-dimensional effects of finite wings, including trailing-vortex formation and induced drag. By modeling the wing as a bound vortex line with spanwise-varying circulation, along with the associated effects of the trailing-vortex wake, this theory provides a sound mathematical framework for determining the spanwise lift distribution, induced drag, and overall aerodynamic efficiency of a finite wing.

History

Lifting-line theory, initially developed by Ludwig Prandtl and his students in the early 1900s, revolutionized the understanding of finite-wing aerodynamics and marked a significant milestone in aeronautical research. Before Prandtl’s work, aerodynamic theories relied primarily on two-dimensional airfoil analysis, which assumed an infinite-span wing and ignored wingtip vortices and the wake’s three-dimensional effects. Frederick Lanchester’s work on the roll-up of wingtip vortices, illustrated in a figure from his book (Figure 2), helped establish the importance of adequately modeling finite-wing aerodynamics. Furthermore, as aviation advanced, particularly during WWI and the transition from biplanes to higher-performance monoplanes in the 1920s, the predictive limitations of existing aerodynamic theories became evident. Designers needed a way to better predict and optimize the performance of not just airfoils but also finite-span monoplane wings, particularly lift and drag, as well as other factors such as aeroelastic characteristics and stalling behavior.

Prandtl’s lifting-line theory was developed in its mature form during World War I and published in his 1918–1919 papers on Tragflügeltheorie, or wing theory.^[2] Prandtl’s lifting-line theory addressed this need by modeling a finite wing as a bound vortex line with a spanwise circulation distribution. Prandtl’s advance over Lanchester’s work was in the theoretical modeling of the wing, in which its wake was represented as a sheet of vortices trailing behind it. This wake produced a downwash, decreased the effective angle of attack along the wing, and generated induced drag, a component directly associated with lift production. Prandtl also demonstrated that an elliptical lift distribution along the wingspan minimizes induced drag, setting a practical goal for efficient wing design. Prandtl’s original published works also show an “ideal” elliptical wing planform, which may have influenced the design of the Supermarine Spitfire in the early 1930s.

In the 1920s, Hermann Glauert^[3] further refined and extended Prandtl’s work, making it more accessible and applicable to practical engineering problem-solving.^[4] The representation of the wing and its wake by a system of horseshoe vortices, as shown in Figure 3 from his book, allowed the induced velocity over the wing to be calculated using the Biot-Savart law, which formed a cornerstone of lifting line theory. He then clarified the fundamental importance of the elliptical lift distribution and introduced systematic methods for analyzing wings of arbitrary planform and twist using a Fourier series, much like the approach employed in the theory of thin airfoils.

Glauert’s detailed methods for calculating induced drag helped bridge the gap between theoretical predictions and practical wing design, making lifting-line theory one of the first valuable approaches for analyzing three-dimensional wings. Additionally, the theory was formulated in a way that later enabled its adaptation to swept wings and non-planar configurations, and it was extended to lifting-surface methods that model the chordwise distribution of circulation. Although more advanced techniques, such as computational fluid dynamics (CFD), are now widely used, lifting-line theory remains a cornerstone of aerodynamic modeling. Its mathematical simplicity and elegance make it valuable for preliminary wing design and essential to aerodynamic education.

Helmholtz’s Vortex Theorems

Helmholtz’s vortex theorems describe the fundamental properties of vortex motion in an ideal fluid. First formulated by Hermann von Helmholtz in 1858, they remain central to understanding vorticity dynamics in fluid mechanics and aerodynamics. These theorems are derived from the Euler equations and describe the behavior of vortex lines and vortex tubes in inviscid flow.

1. First Helmholtz Theorem (Strength of a Vortex Tube): The circulation around any closed curve enclosing a vortex tube is constant along the tube’s length. Therefore, a vortex tube has the same strength at every cross-section. This result follows from the divergence-free nature of the vorticity field, i.e., , which implies that vortex tubes cannot begin or end in the fluid interior.
2. Second Helmholtz Theorem (Motion of Vortex Lines): Vortex lines move with the fluid. If a vortex line passes through a fluid element at some initial time, it will continue to pass through that fluid element at all later times. This implies that vortex lines are “frozen” into the fluid; i.e., the topology of the vortex structure is conserved as it is convected by the flow.
3. Third Helmholtz Theorem (Persistence of Vorticity): In an inviscid fluid, a vortex tube cannot spontaneously appear or vanish within the fluid. A vortex tube must either form a closed loop, extend to infinity, or terminate at a boundary. This result is consistent with Kelvin’s Circulation Theorem, which states that the circulation around a material loop moving with the fluid remains constant over time, i.e., , for an inviscid fluid with conservative body forces.

Helmholtz’s theorems are foundational in the study of vortex dynamics, including the formation and evolution of wingtip vortices and the inviscid modeling of aerodynamic flows. They provide the physical basis for representing finite wings using bound vortices and trailing vortex sheets, which naturally leads to the lifting-line theory and other potential-flow models that underpin much of classical aerodynamics. While real fluids are viscous and may deviate from these idealized results under certain conditions, Helmholtz’s laws provide an essential framework for understanding vortex flow behavior.

Theory of Vortex Lines

The line vortex is a fundamental singularity in fluid dynamics used to model the velocity field induced by a circulation . The concepts of circulation and the effects associated with vortices have been introduced previously. A vortex line is a curve of any shape that is everywhere tangent to the local vorticity vector, as shown in Figure 4. A line vortex in three dimensions is analogous to a point vortex in two dimensions, but instead of a singular point, it extends along a curve. The circulation around any closed curve equals the sum of the strengths of the line vortices that intersect the surface bounded by the curve. Consequently, a line vortex cannot terminate within a fluid. It must form a closed loop, end at a solid boundary, or extend to infinity; Helmholtz’s theorems formally embody these principles.

The Biot-Savart law describes the velocity field induced by the circulation distribution . Regarding the situation shown in Figure 4, the increment in the downwash, , from the vortex element can be mathematically expressed as

(1)  

where is the circulation or vortex strength, is a differential element of the line vortex, is the distance of the point P from the element, and is the angle between the direction of the element and the straight line joining the element to the point P. Notice that the element of a line vortex cannot exist by itself and only forms the basis for integrating the effects along a line vortex of finite length.

Straight Line Vortices

Straight-line vortices are used to develop finite-wing theory. Consider the evaluation of the induced velocity of a straight-line vortex of finite length AB, as shown in Figure 5. If PN has a perpendicular length , then the induced velocity at P from the element at point Q is

(2)  

To show this result, the element may be geometrically expressed as

(3)  

where the angle is as shown in Figure 5, so that

(4)  

The total induced velocity is then obtained by integration along the length of the vortex line using

(5)  

which gives

(6)  

Doubly-Infinite Straight Line Vortex

If, as a special case, the line is of doubly infinite length, as shown in Figure 6, then it will be apparent that , so the general result reduces to

(7)  

which is consistent with the result of a two-dimensional point vortex.

Singly-Infinite Straight Line Vortex

For a singly infinite vortex, which starts at point N and then extends to infinity only in one direction, as shown in Figure 7, then and , and the downwash is given by

(8)  

This latter result for the downwash from a semi-infinite line vortex is widely used in developing lifting-line theory, in which several line vortices are appropriately combined to represent the flow produced by a finite wing.

Lifting Line Theory

The following exposition of lifting-line theory and the development of the monoplane wing equation follows Glauert’s work. In addressing the problem of a wing of finite span in three-dimensional flow, the following assumptions were made:

1. The wing’s chord is relatively small compared to its span, resulting in a high aspect ratio.
2. The wing’s spanwise axis can be taken as a straight line perpendicular to the freestream, i.e., the wing is unswept.
3. The wing planform is symmetric about its centerline; i.e., the left and right wing panels are mirror images of each other.

Apart from these restrictions, the planform (chord), pitch angle (including any twist), and zero-lift angle of the airfoil section may vary arbitrarily along the wing’s span.

Bound & Trailed Vortices

If a wing generates lift, there must be a flow circulation around each airfoil section comprising the wing, effectively creating a line vortex or a set of line vortices along the wing’s span. These line vortices, which are attached to the wing, as shown in Figure 8, are referred to as the bound vortices and create lift in accordance with the Kutta-Joukowski theorem. Physically, these bound vortices arise from the integrated effect of vorticity in the boundary layer, or from an equivalent vortex sheet surrounding the airfoil surface, as in the development of thin airfoil theory. According to the general theory of vortex flows and Helmholtz’s theorems, these bound vortices cannot terminate within the fluid and must continue into the wake as free line vortices. These are called trailing vortices.

This entire vortex system, which extends to infinity far behind the wing, is completed by a transverse vortex parallel to the span. However, this vortex has no consequence in the steady state as its induced effects are asymptotically zero. For practical purposes, the trailing vortices can be considered to align with the freestream flow and extend downstream to infinity as singly infinite line vortices.

Horseshoe Vortex System

Glauert describes the simplest type of vortex system as one in which the circulation around each airfoil section is assumed to be constant, , across the wing’s span. In this special case, the bound vortex system may be represented as a single “lumped” line vortex of strength , which, to be consistent with thin airfoil theory, can be positioned at the aerodynamic center at 1/4 chord of the wing, as shown in Figure 9. The trailing vortices will consist of two line vortices, each with the same strength magnitude but opposite sense, originating at the wing tips and extending downstream in the direction of the freestream flow. This representation is called a horseshoe vortex system.

In practice, the vortices trailing from the wing tips are generally not straight lines because downwash varies with distance behind the wing. Experiments show that as they develop downstream, the trailing vortices tend to descend below the wing and contract inward. However, as Glauert explains, they can be approximated as straight lines parallel to the direction of motion for most practical purposes. This assumption yields a simplified representation of the wing and its wake as a horseshoe vortex system, which underpins lifting-line theory.

Nested Horseshoe System

In reality, an airfoil’s vortex system is more complex because the circulation is not uniform across the span. Instead, it typically peaks at the center and tapers to zero at the tips. This distribution can be represented by superimposing multiple simple “horseshoe vortex systems,” as shown in Figure 10. This setup more accurately represents the vortex system over the wing and in its wake.

This model then comprises interlaced bound vortices, all located at the aerodynamic center on the 1/4-chord axis, and a sheet of trailing vortices extending from the wing’s trailing edge. Glauert used this vortex-system representation as the foundation for developing the mathematics of the lifting-line model.

Wake Rollup

As Glauert explains, the origin of the trailing vortex wake system can also be understood through the pressure difference between the wing’s upper and lower surfaces. Suppose the lift distribution across a wing’s span peaks at its centerline. In that case, pressure above the wing decreases significantly and pressure below it increases, as shown in Figure 11. These pressure differences decrease toward the wing tips as the lift distribution tapers to zero.

As the streamlines above the wing flow inward toward the center and those below the wing flow outward, they leave the trailing edge and form a discontinuity surface, as shown in Figure 11. The trailing vortices then represent the vorticity of this discontinuity surface. The sheet of trailing vortices rolls up into a pair of concentrated vortices downstream, forming the characteristic pair behind a wing, as observed in practice with flow visualization methods.

Closer to the wing, the trailing vortex system can be modeled by assuming the individual trailing vortices extend downstream as straight lines. For the wake farther from the wing, which tends to contract somewhat, it is perhaps more accurate to assume a horseshoe vortex system with a span a little shorter than the wingspan. Both representations are sufficiently accurate for modeling the near-wing induced velocities used in lifting-line theory.

Calculating the Induced Velocity

In general, the circulation varies along the wing span, is symmetric about the centerline, and decreases to zero at the tips. Consistent with the nested horseshoe vortex model, between the points and of the span of the wing, the circulation can be assumed to decrease by an amount

(9)  

Therefore, a trailed wake filament of strength originates from the span element , as shown in Figure 12.

Therefore, in the aggregate, a sheet of trailing vortices extends along the entire wingspan because varies continuously along the wing, i.e., . The induced downwash at any point along the span will then be the sum of the effects of all the trailing vortices of this sheet. For example, consider some reference point on the wing, say . Using Eq. 8, the downwash at this point from all of the trailing line vortices comprising the sheet will be

(10)  

where the numerator is the strength per unit span of the trailing vortex sheet. This result applies to every section along the wing span, and different values of result depending on the distribution and proximity to the wing tips. As noted above, the circulation gradient is usually much higher near the wing tips, so the trailed wake has greater local strength there. This behavior is consistent with Lanchester’s interpretation of the wake generated behind a wing. The question now is how to calculate this downwash distribution, as it affects the lift and drag at each wing section, , and the overall wing lift and drag.

Sectional Effects of the Induced Velocity

The induced velocity from the trailing wake will affect the aerodynamic angle of attack at each section across the wing. Adding the downwash flow vector to the freestream flow vector produces a resultant flow velocity (or local relative wind) that turns through an angle , called the induced angle of attack, as shown in Figure 13. It will be apparent that the resultant flow now approaches the wing at a different angle, thereby reducing the “effective” angle of attack, .

From the geometry, as shown above, the induced angle, , is given by

(11)  

which, in general, will differ at each station on the wing, i.e., and . For small angles, which are typical for a wing, it is sufficient to write that

(12)  

Therefore, the effective (and lower) angle of attack of the wing section is now

(13)  

This means the lift per unit span is reduced relative to its two-dimensional value, i.e., the lift obtained in the absence of downwash. Notice that because will typically be much smaller than , then the resultant velocity, , can be assumed to be equal to , i.e.,

(14)  

It is particularly significant in this context that, on a finite wing, the aerodynamic force at each section is slightly tilted rearward relative to the freestream because of the induced downwash, thereby reducing the vertical component of lift at a given angle of attack. As a result, the lift vector is locally tilted rearward by the induced angle of attack, , so a component of the lift acts downstream. This component is the induced drag, , i.e.,

(15)  

This result shows that induced drag directly results from the wing’s three-dimensional vortex system, which produces the downwash velocity and hence the induced angle of attack .

Development of the Monoplane Wing Equation

These results provide the elements needed to formulate the lifting-line problem for a finite wing. The circulation determines the local lift per unit span through the Kutta-Joukowski theorem, while the spanwise variation of generates a trailing vortex sheet that induces a downwash velocity over the wing. This downwash reduces the effective angle of attack at each spanwise station and, therefore, couples the local lift on the wing to the circulation distribution over the entire span. The objective is to combine these relationships into a single governing equation for , the fundamental equation of Prandtl’s lifting-line theory, commonly called the monoplane wing equation.

Connection Between Circulation and Downwash

As previously shown, if is the circulation produced around any section of the wing, the downwash from the trailing wake system at a point of the span is determined using

(16)  

The notation P.V. indicates that the singular integral is evaluated in the Cauchy principal-value sense. A representative wing section, therefore, experiences the lift force corresponding to the two-dimensional conditions at the effective angle of attack, i.e.,

(17)  

If the geometric angle of attack is measured relative to the wing reference line, and is the zero-lift angle of attack of the airfoil section, then the effective angle of attack relative to zero lift is

(18)  

Therefore, the lift coefficient of the wing section will be

(19)  

where is the slope of the lift coefficient curve with respect to angle of attack for the wing section in two-dimensional flow.

Furthermore, the circulation around the wing section can be connected to the lift coefficient, , using

(20)  

This means that

(21)  

This equation relates the local circulation to the local chord, the local geometric angle of attack, the section’s zero-lift angle, and the downwash induced by the entire trailing vortex sheet. Therefore, the circulation at any spanwise station depends not only on the local wing geometry but also on the circulation distribution over the complete span.

When the circulation and the downwash of any monoplane wing have been determined, the lift and induced drag are obtained by evaluating the integrals

(22)  

Method of Solution

Glauert used a transformation^[5] to replace the spanwise coordinate , measured from the wing centerline, by the angle , as defined by

(23)  

where is the semi-span. As varies from to , varies from to across the span from the left tip to the right tip. The circulation, , which is a function of , can then be expressed as the Fourier sine series

(24)  

where the values of the coefficients must be determined by connecting the distribution of “bound” to the wing with the downwash produced by the trailed wake consistently.

The chosen series for the circulation satisfies the condition that the circulation vanishes at the wing tips. As shown in Figure 14, if the wing geometry and loading are symmetric about the centerline, only odd integer values of will appear in the Fourier series.

The primary mode is the or mode, corresponding to an elliptical circulation distribution over the wing. Because , then

(25)  

Simplifying gives

(26)  

and so

(27)  

which is elliptical. The higher symmetric modes, i.e., , with and shown in Figure 14, represent a symmetric modification to the primary loading, i.e., a deviation from the elliptical distribution.

The downwash at the point or of the wing now becomes

(28)  

because

(29)  

Therefore, at the general point of the wing, the downwash is

(30)  

Monoplane Wing Equation

The equation connecting the circulation and the downwash velocity now becomes

(31)  

which gives

(32)  

and finally

(33)  

Equation 33 is called the fundamental equation of lifting line theory or the monoplane wing equation and can be used to determine the coefficients for any monoplane wing.^[6] It is generally solved numerically. The equation must be satisfied at all points of the wing, but because the wing is symmetrical about its centerline in level flight, it is usually sufficient to consider values of between and . In this case, only the symmetric modes, i.e., the odd values , must be evaluated. The dimensionless Glauert parameter is proportional to the local chord and may vary across the span. The geometric angle of attack and zero-lift angle may also vary across the span, for example, if the wing has twist or washout.

Lift and Induced Drag

The lift on the wing is given by

(34)  

Switching to the spanwise angle coordinate, , where and , the lift becomes

(35)  

Substituting the general series expression for circulation gives

(36)  

and so

(37)  

Using the orthogonality property of sine functions, i.e.,

(38)  

then the only term that contributes to the integral is , so

(39)  

The coefficient , which is the elliptical loading mode, is related to the lift coefficient using

(40)  

The induced drag is determined using

(41)  

Switching to the coordinate gives

(42)  

Recall that the downwash velocity is given by

(43)  

Substituting and gives

(44)  

Expanding the summation gives

(45)  

Again, the orthogonality of functions ensures that only terms where are retained, so that

(46)  

It is convenient to write

(47)  

so the induced drag, , may be written in force form as

(48)  

where is the wing span. Equivalently, because , this result may be written as

(49)  

In terms of the drag coefficient, then

(50)  

The total drag coefficient includes both profile drag and induced drag. If the profile drag coefficient varies across the span, the profile drag coefficient for the wing is

(51)  

Adding the induced drag gives

(52)  

Asymmetric Loading Modes

The even modes in lifting-line theory represent an asymmetric lift distribution about the wing centerline. These modes can arise from aileron deflections, rolling motion, yawing motion, sideslip, or other spanwise asymmetries. This lift asymmetry influences rolling moments, possible yawing moments, and induced drag. For example, ailerons are control surfaces near the wingtips, deflected asymmetrically to induce a rolling moment about the aircraft’s longitudinal axis. Figure 15 shows that their operation produces asymmetric modes in lifting-line theory, particularly and other higher even-numbered modes.

The total circulation distribution can be expressed as

(53)  

where represents the primary symmetric lift distribution associated with the freestream, and captures the first antisymmetric lift contribution associated with aileron deflection, as shown in Figure 16. The primary contribution from the ailerons is then

(54)  

The coefficient depends on the aileron deflection angle , which modifies the local angle of attack near the wingtips. For small deflection angles, is approximately linear in , but the proportionality depends on the wing’s aspect ratio, planform, and the spanwise extent and effectiveness of the ailerons.

The resulting rolling moment, , about the longitudinal axis is directly related to this antisymmetric lift distribution, i.e.,

(55)  

where is the differential change in lift per unit span caused by the aileron deflection at the spanwise location . This rolling moment scales with , i.e.,

(56)  

From a design perspective, proper aileron sizing and placement are crucial to ensure sufficient rolling authority without excessive induced drag or adverse yaw effects. To this end, lifting-line theory can provide valuable insights into the aileron design that may be required.

Elliptic Spanwise Loading

The lift and induced drag of a wing are given by

(57)  

For a wing of given span, planform area, and lift coefficient, the coefficient has a fixed value because

(58)  

Therefore, the induced drag is minimized when all higher-order coefficients () in the circulation series are zero, i.e., when . In this case, the distribution of circulation across the span of the wing simplifies to

(59)  

so the circulation is elliptically distributed. Putting into Eq. 33 gives

(60)  

This equation can be satisfied with a constant value of over the span only if the local value of varies consistently with the chord distribution. For an untwisted elliptical wing with a constant airfoil section, then

(61)  

and therefore

(62)  

Substitution then gives

(63)  

which is independent of spanwise position. The corresponding induced drag coefficient will be

(64)  

Notice that the downwash becomes

(65)  

and the equivalent angle of attack is

(66)  

along the span. This result implies that the local lift coefficient, , is constant across the span, as shown in Figure 17.

To produce a constant and an elliptical circulation distribution, the chord length must vary elliptically along the span. This result can be obtained by using

(67)  

where it will be apparent that

(68)  

so that a solution is

(69)  

where is the chord length at the wing root. This special case of elliptic loading is significant for several reasons:

1. It results in the minimum possible induced drag for a given lift and span, i.e., .
2. It produces a constant downwash over the span.
3. For an untwisted wing with a constant airfoil section, the sectional lift coefficient is constant along the span.
4. Elliptic loading, or a close approximation, can also be achieved with non-elliptic planforms by suitably varying the twist angle and local chord along the span.

Effect of Aspect Ratio

The results developed for the lift and drag coefficients of an elliptically loaded wing can be used to calculate the effect of a change in aspect ratio. If the aspect ratio is changed from to , the change in the induced drag coefficient at a given value of the lift coefficient is

(70)  

It is apparent that wings with a higher aspect ratio have aerodynamic advantages, all other factors being equal, such as wing area and planform shape.

While this formula for the induced drag coefficient applies exactly only to wings with elliptic loading, the lift distribution curves for rectangular and most airplane wings do not differ significantly from the elliptic form. Therefore, the formula may be used more generally to estimate the effect of a modest change in aspect ratio.

Effect on Lift Curve Slope

The assumptions are: 1. The two-dimensional lift curve slope for the airfoil is , which assumes no three-dimensional effects. 2. The lift curve slope for the finite wing is , which includes the effects of finite aspect ratio, . Lifting-line theory introduces an induced angle of attack caused by the downwash from the wing’s trailing vortices. If the angle of attack is measured relative to the zero-lift angle of the airfoil section, then the effective angle of attack is related to the geometric angle of attack and the induced angle by

(71)  

The total lift coefficient is then

(72)  

where denotes the two-dimensional lift curve slope of the airfoil section. If the geometric angle of attack is measured relative to a wing reference line, then in this derivation must be replaced by . For an elliptically loaded wing, the induced angle of attack is related to the lift coefficient and aspect ratio by

(73)  

Substituting this result back into the lift equation gives

(74)  

and rearranging gives

(75)  

Dividing through by gives

(76)  

Therefore, for wings with a high aspect ratio or a small , the lift curve slope approaches the two-dimensional value, corrected by an aspect ratio term to account for three-dimensional effects.

Numerical Solution to the Monoplane Wing Equation

For an elliptically loaded wing, the lifting-line equation has a simple analytical solution because only the first circulation mode, , is present and the downwash is constant across the span. However, most practical wings do not have exactly elliptic loading. The chord, twist, section lift-curve slope, and zero-lift angle may vary along the span, so the resulting circulation distribution generally contains several Fourier modes. In such cases, the monoplane wing equation cannot usually be solved in closed form. Instead, the circulation series is truncated after a finite number of terms, and the equation is enforced at selected spanwise stations to obtain a system of simultaneous linear equations for the unknown coefficients .

The governing equation for the monoplane wing is

(77)  

where:

1. are the coefficients of the circulation expansion.
2. is referred to as Glauert’s nondimensional parameter.
3. is the local geometric angle of attack, or more correctly, the local pitch angle of the wing section.
4. is the zero-lift angle of attack of the local airfoil section.
5. is the local chord length, which may vary with and so with .
6. is the wing semi-span.

To solve the equation numerically, the circulation expansion is truncated after terms, and the equation is enforced at selected spanwise stations. A convenient choice is

(78)  

which avoids the singular endpoints at and . At each discrete point , the governing equation becomes

(79)  

where , , and are evaluated at the spanwise station corresponding to . This equation can be rewritten in matrix form as

(80)  

where is the vector of unknown coefficients, is the right-hand side vector, and is the influence coefficient matrix.

The components of the matrix equation are

(81)  

and

(82)  

The matrix equation is expanded as

(83)  

where and .

The unknown coefficients are found by solving the linear system

(84)  

Efficient numerical methods, such as LU decomposition, are typically used to solve this system when is large.

Once the coefficients are determined, the following aerodynamic quantities can be computed:

1. The spanwise circulation distribution, , i.e.,

   (85)  

2. The sectional lift coefficient distribution, obtained from the Kutta-Joukowski theorem, i.e.,

   (86)  

3. The total lift coefficient, which is directly related to , i.e.,

   (87)  

4. The induced drag coefficient, which is computed using

   (88)  

Wing Shape & Loading Distributions

The wing’s planform, twist, and interference effects strongly influence the spanwise lift distribution and overall wing performance. The planform shape, including the wing’s aspect ratio, taper ratio, and sweep angle, determines the baseline spanwise lift distribution and affects induced drag and aerodynamic efficiency. Whether geometric or aerodynamic, wing twist lets designers tailor the angle-of-attack distribution, optimizing lift and reducing induced drag under specific flight conditions. Interference effects, such as those arising from wing-body junctions, nacelles, or wingtip devices, further modify the flow field around the wing, impacting both the spanwise lift distribution and the overall drag characteristics. Wing design must carefully consider these factors to balance aerodynamic efficiency, structural feasibility, and operational requirements.

Planform Effects

The theoretically optimal aerodynamic efficiency () and the lowest induced drag are achieved with an elliptically loaded wing. For an untwisted wing with a constant airfoil section, this condition is achieved with an elliptical planform. This wing shape produces elliptical spanwise aerodynamic loading (i.e., lift per unit span) and uniform downwash over the wing, which, as noted above, is theoretically the condition for minimum induced drag. However, practical wing design rarely achieves exact elliptic loading because of structural, manufacturing, control, and integration requirements. Values of for a plain wing can range from about 0.01 to 0.1, as shown in Figure 18; values greater than 0.1 are usually considered poorly designed. However, even a rectangular planform wing will have a reasonably good value of (probably between 0.05 and 0.1) if it has a decent aspect ratio, i.e., , and also uses some wing twist or washout to control the spanwise lift distribution.

The wing shape, i.e., its planform, affects the lift distribution, the formation of the trailing vortex system, and so the magnitude of the induced drag on the wing. The idea is to use variations in wing chord along the span, and perhaps wing twist and airfoil sections, to closely approximate the ideal elliptical spanwise loading and minimize . The effective aspect ratio can also be improved somewhat by paying attention to the shape of the wing tips, which can lower if they are suitably shaped (e.g., with more rounded contours) or by adding a winglet.

Incorporating Wing Twist

In lifting-line theory, wing twist affects the circulation distribution by altering the local geometric angle of attack along the span, which modifies the local lift coefficient. The effective angle of attack at any spanwise location is now given by

(89)  

where is the local geometric angle of attack, including the effect of twist, and the induced angle of attack is

(90)  

Wing twist, denoted as , modifies the geometric angle of attack by

(91)  

where is the root angle of attack and positive denotes washout, i.e., a reduction in local incidence toward the tip. Substituting this result, the effective angle of attack becomes

(92)  

Notice that in the case of no twist, i.e., (), the geometric angle of attack is constant, and the circulation distribution depends only on the wing shape and induced effects. In the case of washout, i.e., when the local incidence decreases toward the tip, lift at the wingtips decreases, lift is redistributed inboard, as shown in Figure 19, and the spanwise loading can be brought closer to the desired distribution. If the wing has wash-in, i.e., when the local incidence increases toward the tip, which is unusual, it increases lift at the wingtips, leading to a higher induced drag and less benign stall behavior. Twist lets designers optimize lift distribution and induced drag for specific flight conditions.

Interference Effects

The photograph in Figure 20 below illustrates the spanwise loading over a wing due to natural condensation in the low-pressure zones, which is nominally elliptically distributed, even on this relatively low-aspect-ratio wing. The highest lift, or highest pressure difference, is at midspan, and the lowest lift, or lowest pressure difference, is at the wing tips. The problem, however, is that even minor deviations from the ideal elliptical form can increase induced drag. Such variations can arise from flow around the fuselage, engine, undercarriage, external stores, and other components that interfere with the flow around the wing.

Examples of the effects of spanwise interference are shown in Figure 21, all at approximately the same total wing lift. Notice that the presence of the fuselage can cause the loading to deviate significantly from the ideal elliptical form, potentially increasing to over 1.15. However, local effects along the wing, such as aerodynamic interference from an engine, tend to have a minor impact on .

Wing Sweep

The sweep of a wing significantly affects its aerodynamic characteristics, especially at higher subsonic and transonic speeds. The sweep angle, , alters the effective flow direction, thereby affecting the airfoil sections and reducing the component of the freestream velocity normal to the leading edge, which modifies the aerodynamic response.

Lifting-line theory can incorporate sweep by modifying the freestream velocity components and the associated expressions for lift and circulation. For a wing with sweep angle , the local airfoil sections respond primarily to the component of the freestream velocity normal to the local sweep line, which in lifting-line theory is usually taken to be the 1/4-chord line. As shown in Figure 22, the component of the freestream velocity normal to the local quarter-chord is

(93)  

This normal-flow approximation modifies the sectional aerodynamic response, but the geometric angle of attack should not be interpreted as simply multiplied by . A more precise statement is that the section lift is governed primarily by the normal component of the flow, whereas the induced angle of attack is still associated with the downwash velocity.

Therefore, in an approximate swept-wing lifting-line treatment, the sectional lift per unit span may be written using the velocity component normal to the local quarter-chord line, i.e.,

(94)  

The Biot-Savart law still governs the induced velocity, or downwash, but must be paired with the modified form of lift. Therefore, the lifting-line equation for a swept wing is usually treated as an approximate extension of the unswept theory rather than as an exact result. In symbolic form, one approximate normal-flow correction may be represented by writing

(95)  

where the induced angle of attack is still related to the downwash by

(96)  

This formulation preserves the structure of the original theory while incorporating the effect of sweep into the normal-flow component. In practice, swept-wing loading is better treated using lifting-surface theory or vortex-lattice methods, especially when the sweep angle is large, or the Mach number is high.

It can be shown that the effect of a constant wing sweep angle on the lift curve slope gives the approximate result

(97)  

where and is the wing aspect ratio. This approximation assumes a two-dimensional incompressible lift-curve slope of approximately per radian and a span-efficiency correction of unity. At low Mach numbers, where , it reduces to the corresponding incompressible swept-wing approximation.

Sweep also changes the spanwise loading and the corresponding induced drag. However, induced drag is still governed primarily by the lift, span, and spanwise loading distribution of the complete wing. Therefore, a simple replacement of by an “effective” value such as should be regarded only as a crude heuristic, not as a general lifting-line result. For swept wings, the induced drag is better retained in the form

(98)  

where accounts for the departure of the actual spanwise loading from the elliptic-loading result.

Winglets

Winglets reduce the induced drag of finite wings by modifying the location and structure of the trailing vortex system generated by lift. Although winglets increase the wing’s wetted surface area and add a small amount of skin friction and profile drag, they typically reduce total airplane drag with little or no increase in wingspan. This trade-off explains the widespread adoption of winglets on modern commercial airliners as an effective strategy to reduce drag and save fuel.

In the following analysis, the wing and its wake are modeled using a classical lifting-line horseshoe-vortex system, modified to include winglets by geometrically displacing the trailing-vortex legs, as illustrated in Figure 23. This is a deliberately simplified model. A real winglet is a lifting surface with its own circulation, side force, viscous drag, interference effects, and wake roll-up. These effects are not included here. Instead, the model isolates the first-order influence of moving the trailing-vortex system farther from the lifting wing.

The induced velocity is evaluated at the wing’s mid-span. Although this effective-span winglet theory does not capture spanwise variation in downwash, it provides a useful approximation of total induced drag based on the Biot-Savart law. It highlights the aerodynamic benefits of winglets by showing how they influence the geometry of the trailing-vortex system.

In this simplified model, the winglet itself is assumed to generate no lift or drag; its effect is modeled purely as a geometric shift in the tip vortex’s starting position, effectively moving it farther from the main lifting wing. As such, this model isolates the influence of vortex positioning on the downwash field without accounting for the winglet’s detailed aerodynamic behavior or the tip vortex’s roll-up. The wing tip vortices are still modeled as semi-infinite straight filaments extending downstream, so the induced vertical velocity field can be calculated using the Biot-Savart law.

When the wing has regular tips, i.e., no winglets, the trailing vortices originate in the plane of the wing at , neglecting the effects of any dihedral angle. The induced vertical velocity at mid-span from one semi-infinite trailing vortex is . Therefore, including both tip vortices gives

(99)  

where is the semi-span and is the wingspan.

Now consider the case when the winglets relocate the vortex tips along a winglet length at a cant angle measured from the vertical. The trailed vortex filaments, i.e., the wing tip vortices, now start from

(100)  

With this convention, corresponds to a vertical winglet, and corresponds to a flat span extension. The spanwise projection of the winglet is , and the vertical projection is .

This approach assumes that the spatial separation of the trailing-vortex system dominates the change in induced drag. This formulation is consistent in spirit with Prandtl-Jones minimum induced drag theory,^[7] which shows that induced drag depends on the distribution and effective spacing of the trailing-vortex system. The present model is much simpler: it assumes that the main effect of the winglet is to displace the tip vortex away from the lifting wing.

Using the Biot-Savart law, the vertical velocity induced at mid-span by one vortex leg is

(101)  

where

(102)  

and

(103)  

Including both tip vortices, the induced velocity becomes

(104)  

which gives the induced angle of attack as

(105)  

The lift generated by the wing is

(106)  

As a consequence, the induced drag is

(107)  

Compared to the baseline wing without the winglet, the induced-drag ratio is

(108)  

or, after defining

(109)  

the result becomes

(110)  

The simplified horseshoe-vortex model is used here to predict the ratio of the induced drag with and without the winglet. This ratio may then be applied to the conventional induced-drag relation by defining an effective aspect ratio based on the inverse of the predicted drag ratio, i.e.,

(111)  

or equivalently,

(112)  

where is the baseline aspect ratio of the wing without the winglet. The model does not independently derive the absolute induced-drag coefficient because it represents the wing using a single constant-circulation horseshoe vortex. Instead, it applies the predicted drag ratio to the conventional finite-wing relation. Therefore,

(113)  

where is the span efficiency factor of the baseline wing. If elliptic loading is assumed, then .

Special cases include:

1. No winglet:

   (114)  

2. Vertical winglet:

   (115)  

3. Flat winglet, equivalent to a span extension in this simplified vortex-displacement model:

   (116)  

Minimizing induced drag from winglets corresponds to maximizing , which depends on the effective location of the trailed vortex system. Define the geometric function

(117)  

Differentiating with respect to gives

(118)  

Therefore, for practical winglet-length ratios satisfying , this function increases with . With the present convention, corresponds to a vertical winglet and corresponds to a flat span extension.

Therefore, in this simplified vortex-displacement model, a vertical winglet can reduce induced drag with essentially no increase in geometric span because the trailed vortex system is displaced vertically. As the winglet is canted outward, its spanwise projection increases, and the model predicts a larger effective aspect ratio. The limiting case corresponds to a flat span extension, which gives the largest value of in this idealized model. In practice, this increase must be balanced against the associated increases in wingspan, wetted area, profile drag, structural weight, and wing-root bending moment. The results summarized in Figure 24 show a simplified trend in that increasing the cant angle and winglet length increases the predicted induced-drag reduction, but the practical benefit diminishes once the additional geometric and structural penalties are included.

A commonly used approach in airplane performance analysis is to estimate the effect of a winglet by using a linear correction to the aspect ratio of the form

(119)  

where is the corrected or effective aspect ratio, is the winglet length used in the present geometric model, and is a coefficient that depends on the winglet geometry. This expression is widely used in preliminary wing design because of its simplicity.

This model, given by Eq. 119, can be reconciled with the theoretical effective-span winglet model by performing a Taylor expansion of the more general result. Recall that the effective aspect ratio based on the trailed wingtip vortex geometry is

(120)  

where . Expanding this expression for small gives

(121)  

Comparing this expansion with Eq. 119, the linear winglet correction coefficient in this particular geometric model is identified as

(122)  

This result shows that the linear model is a first-order approximation to the present vortex-displacement model only when is interpreted as the length of the canted winglet. With the convention used here, a typical cant angle of gives . If the correction is written in terms of the winglet’s vertical height, then a different coefficient is needed because the vertical height is . In practice, empirical values of can be larger or smaller depending on the winglet geometry, toe angle, airfoil section, loading, viscous drag penalty, and the winglet’s interaction with the rolled-up tip vortex. Therefore, while Eq. 119 provides a convenient approximation, the full geometric expression is more consistent with the simplified vortex-displacement model, especially for larger winglet lengths or varying cant angles. In both cases, increasing the spatial separation between the trailing vortices, whether through increased span and/or vertical offset, reduces induced drag on the wing, consistent with the well-established aerodynamic effects of winglets.

Formation Flight

Formation flight has long been recognized as a way to reduce aerodynamic drag, most notably in birds that migrate in stable V-formations. The same principle has been explored for long-range ferry flights, where small reductions in induced drag can add up to substantial fuel savings over many hours, thereby extending flight range. The underlying mechanism is geometric: a following flyer can position itself in the upwash region generated by the leader’s tip vortices, thereby reducing its own downwash, induced angle of attack, and induced drag. As with the simplified winglet analysis, this effect can be modeled using a lifting-line representation in which each flyer is modeled as a classical horseshoe vortex.

In practice, formation flight is fully three-dimensional: the follower is displaced laterally, vertically, and by a finite streamwise distance behind the leader. For a semi-infinite trailing vortex aligned with the -axis, the induced velocity at the follower depends on this streamwise separation. Therefore, the cross-plane formulation developed below corresponds to the far-wake approximation, in which the follower is sufficiently far downstream that each trailing vortex may be treated locally as an effectively infinite straight vortex. Consider a leader of semi-span whose trailing vortices originate at , as shown in the schematic in Figure 25.

The follower is located at a lateral and vertical offset relative to the leader’s mid-span and reference plane. In the far-wake approximation, the vertical velocity induced at the follower’s mid-span by one effectively infinite trailing vortex passing through is

(123)  

with the sign depending on the sense of the vortex filament. Because the two trailing vortices have opposite senses, their vertical velocity contributions must be combined with opposite signs. By introducing nondimensional offsets given by

(124)  

then the induced velocity from the leader can be written in terms of the geometric function

(125)  

With this sign convention, positive values correspond to additional downwash, and negative values correspond to upwash. Beneficial formation positions correspond to , which physically represent the upwash regions just outboard and slightly below the leader’s tip vortices.

Let denote the induced drag of the follower in isolation. The velocity induced by the leader is

(126)  

whereas the induced angle of attack of the isolated follower is

(127)  

Therefore, adding the downwash or upwash from the leader modifies the induced drag according to

(128)  

where

(129)  

accounts for differences in aspect ratio and lift coefficient. For similar aircraft operating at comparable , the value of is close to unity, so the change in induced drag is governed primarily by the geometric function . This expression naturally leads to an “effective aspect ratio” for the follower, analogous to the effective-span interpretation used for winglets, i.e.,

(130)  

Therefore, the drag benefit of formation flight emerges directly from the displacement of the leader’s tip-vortex system relative to the follower’s lifting line. Positions where increase the follower’s effective aspect ratio, thereby reducing its induced drag through the same geometric mechanism that winglets use on a single wing.

The preceding analysis describes the change in induced drag for a follower behind a single leader. In a formation containing several aircraft, however, each flyer is influenced by the vortex systems of all aircraft ahead or beside it. Because the induced velocity of a vortex system is linear in circulation, the net vertical velocity at any follower is obtained by superposition. If the airplanes are indexed by , with the follower denoted by , then each leader contributes a geometric function of the same form as , shifted to that aircraft’s lateral and vertical position. For identical aircraft operating at comparable lift coefficients, the total geometric interference function becomes

(131)  

where is the nondimensional position of airplane relative to the follower and where the original definition of is retained. Each term in Eq. 131 corresponds to a single horseshoe vortex system, so a follower experiences upwash or downwash resulting from the combined action of all aircraft ahead. A typical variation of is shown in Figure 26. Negative values correspond to upwash regions that reduce the follower’s induced drag.

The induced-drag ratio for the follower generalizes accordingly, i.e.,

(132)  

with the same coefficient

(133)  

As before, for similar aircraft at similar lift coefficients, is close to unity, so the total induced-drag modification is determined primarily by the combined geometric function .

Equation 132 naturally leads to a generalized effective aspect ratio for the follower in a multi-airplane formation, i.e.,

(134)  

This expression shows that each additional aircraft ahead contributes to the overall upwash or downwash experienced by the follower. When the follower occupies a region where , its effective aspect ratio is increased, and its induced drag is reduced. This is the same geometric mechanism by which a winglet increases the effective span of a single wing, but here the effect arises from the displacement of multiple tip-vortex systems in three-dimensional formation flight.

Consequently, optimal formation patterns, such as the familiar V-formation, emerge as aircraft, or birds, position themselves within the negative range of the lobes generated by the group ahead, thereby distributing themselves so that each follower benefits from the upwash produced by several leaders. The photograph below (Figure 27) shows five F-15Cs in formation during a ferry flight, a method used to conserve fuel over long distances. Because the lead airplane has the highest drag, the pilots periodically rotate positions.

Introduction to Lifting Surface Theory

The progression from lifting-line theory to lifting-surface theory and, ultimately, to thickness modeling represents a systematic refinement in aerodynamic analysis using vortex and potential-flow methods. Each step increases the model’s fidelity by incorporating additional aspects of the three-dimensional flow around finite wings. Lifting-line theory is the simplest of these models, representing the wing as a single bound vortex filament near the quarter-chord that extends along the span. While analytically elegant and computationally efficient, lifting-line theory does not account for chordwise loading variations or wing thickness. This limitation matters for swept, tapered, low-aspect-ratio, or non-planar wings, where a single spanwise circulation distribution cannot adequately represent the aerodynamic loading.

Lifting-surface theory is the natural two-dimensional extension of Prandtl’s lifting-line theory. In lifting-line theory, the circulation varies only along the span, i.e., , and the governing equation relates the spanwise circulation distribution to the induced downwash. As shown in Figure 28, lifting-surface theory allows circulation to vary over both span and chord, i.e., . In this way, the wing is treated as a continuous lifting surface rather than a single lifting line, allowing a more realistic representation of the aerodynamic loading and its chordwise variation.

The fidelity improves by resolving the variation in aerodynamic loading across the planform. In lifting-line theory, all bound circulation is collapsed onto a single spanwise line, so the theory can predict the spanwise lift distribution and induced drag but cannot represent chordwise pressure loading, local pitching moment, aerodynamic-center movement, or the detailed effects of sweep and low aspect ratio. In lifting-surface theory, the loading is distributed over both and , so the method can account for the fact that the leading edge, mid-chord region, and trailing edge do not generally carry the same aerodynamic load. This additional chordwise resolution is especially important for swept, delta, tapered, and low-aspect-ratio wings, where the flow over one part of the chord can strongly influence the loading elsewhere on the surface.

Classical Lifting-Surface Theory

Lifting-surface theory is a continuum theory, not just a numerical method. In its classical form, the wing is represented by a continuous distribution of vorticity, doublets, or related singularities over the planform or mean camber surface. The unknown loading is then determined by satisfying the flow-tangency condition over the lifting surface. For idealized planforms and loading distributions, this problem can sometimes be treated analytically or semi-analytically.

Early examples include Blenk’s treatment of the monoplane as a lifting vortex surface^[8] and Multhopp’s later subsonic lifting-surface method.^[9] Blenk’s method replaced the single spanwise lifting line with a distributed vortex surface over the wing planform. In modern notation, the induced velocity at a point on the lifting surface may be written as

(135)  

where is the vortex strength distributed over the lifting surface, and is the kernel giving the induced velocity at from a unit vortex element at . This form shows the essential difference from lifting-line theory: the loading is distributed over an area rather than along a spanwise line.

Multhopp’s method was a more systematic subsonic lifting-surface method for calculating the wing lift distribution. The continuous lifting surface was represented by a finite set of unknown load values, and the downwash condition was imposed at selected points on the wing planform. In matrix form, the method has the general structure

(136)  

where is an aerodynamic influence coefficient, is an unknown loading or circulation coefficient, and the imposed angle of attack, twist, camber, and local surface orientation determine . Multhopp’s method therefore serves as a bridge between analytical lifting-surface theory and subsequent numerical vortex-lattice methods.

These classical theories are important because they show that lifting-surface theory is neither inherently a vortex-lattice method nor inherently numerical. However, exact or closed-form solutions exist only for special geometries or loading assumptions. For practical aircraft wings with arbitrary sweep, taper, twist, camber, and planform shape, the lifting-surface equations are usually solved approximately.

Numerical Lifting-Surface Methods

Vortex-lattice methods are practical numerical implementations of lifting-surface theory. The bound circulation is distributed over both span and chord, allowing a spatial representation of aerodynamic loading over the planform. In practice, the surface is discretized into a lattice of quadrilateral panels, each containing a vortex element or vortex ring, as shown in Figure 29. The unknowns are the strengths of these vortex elements, and their values are determined by enforcing the flow-tangency condition at selected control points on the lifting surface.

The induced velocity at a control point at is the cumulative effect of all vortex elements. For the normal or vertical component of the induced velocity, this influence may be written schematically as

(137)  

where is the influence coefficient for a unit-strength vortex element on panel , and is the unknown circulation strength of that element. Each coefficient represents the velocity induced at the control point by that unit vortex element, thereby relating the distributed circulation over the lifting surface to the induced velocity field.

As in lifting-line theory, the unknown circulation strengths are determined by enforcing the flow-tangency (no-penetration) boundary condition at control points on the lifting surface. This condition requires that the normal component of velocity at the surface be zero, i.e.,

(138)  

where is the total velocity at the surface and is the local unit normal vector. The total velocity is the sum of the freestream and induced components, i.e.,

(139)  

so the boundary condition becomes

(140)  

or equivalently,

(141)  

For most classical vortex-lattice implementations, the wing is assumed thin and represented by a mean camber surface. The boundary condition is then applied at control points on this surface rather than on a finite-thickness body. For a locally flat surface at a small angle of attack, the normal component of the freestream velocity is approximately , where represents the local geometric angle of attack, including any effects of twist and camber slope under the small-disturbance approximation. The boundary condition then reduces to a scalar equation involving the induced vertical velocity component, i.e.,

(142)  

with the sign depending on the convention used for positive downwash and positive angle of attack.

Equivalently, in lifting-line notation, the induced angle of attack is

(143)  

so the local effective angle of attack is

(144)  

This expression has the same physical interpretation as in lifting-line theory: the induced velocity field changes the local angle at which the flow meets the lifting surface.

Substituting the influence-coefficient expression for the induced velocity into the boundary condition gives a linear system of simultaneous equations for the unknown circulation values over the lifting surface, i.e.,

(145)  

for the simple flat-surface convention used above. More generally, the right-hand side is obtained from , with the local normal vector accounting for angle of attack, camber, twist, dihedral, and local surface orientation. This system is the numerical lifting-surface analog of the lifting-line equation and yields the circulation distribution over the entire lifting surface.

Once the circulation distribution is obtained, the corresponding aerodynamic loading can be determined. The local pressure difference across the lifting surface is related to the local vortex strength, and the lift, pitching moment, rolling moment, and induced drag can be obtained by integrating or summing the contributions of all vortex elements. Compared with lifting-line theory, numerical lifting-surface methods better represent chordwise loading, aerodynamic-center movement, pitching moment, and the effects of planform geometry. For example, a swept or low-aspect-ratio wing may have strongly nonuniform chordwise loading, so a single spanwise lifting line cannot recover the local pressure forces and pitching moments. A lifting-surface method can represent these effects because each chordwise row of panels carries its own circulation or vortex strength.

Thickness and Surface-Panel Methods

Lifting-surface theory provides significantly greater fidelity than lifting-line theory, particularly for wings with moderate or low aspect ratios and for configurations with sweep, taper, or other planform complexity. However, in its classical thin-surface form, it still represents the wing by a mean lifting surface. It represents the lifting action of the wing and can include camber through the shape of the mean surface or through the boundary condition, but it does not directly represent the displacement effect of finite thickness.

To incorporate thickness effects, additional singularity distributions must be introduced, and the no-penetration boundary condition must be applied on the actual body surface rather than only on a mean lifting surface. This step leads to surface-panel methods. In modern panel methods, the aerodynamic surface is discretized into panels carrying source, doublet, or vortex singularities. Source distributions are commonly used to represent thickness, while doublet or vortex distributions are used to represent lift.

The total perturbation potential may be written schematically as

(146)  

where represents the lifting contribution from an equivalent vortex or doublet distribution, and represents the source contribution associated with thickness. The source contribution has the form

(147)  

where is the source strength distribution over the surface, and the distance between the field point and the source point is

(148)  

The distinction between theory and method is important. Lifting-surface theory is the continuum formulation in which vorticity, doublets, or related singularities are distributed over a thin wing surface. Classic examples, including Blenk’s lifting-surface calculation for a monoplane and Multhopp’s subsonic lifting-surface method, show that lifting-surface theory is broader than the vortex-lattice method. A vortex-lattice method is one numerical realization of lifting-surface theory for thin wings. A surface-panel method is a finite-thickness potential-flow method for bodies with actual three-dimensional shape. Both are based on singularity distributions and no-penetration boundary conditions, but they represent different levels of geometric fidelity. Classical vortex-lattice methods primarily account for the lifting effects of thin surfaces, whereas contemporary panel methods combine vortex, doublet, and source distributions to model camber, thickness, and complex aircraft geometry.

Such lifting-surface and surface-panel formulations form the theoretical foundation of many modern aerodynamic analysis methods for wings and complete aircraft configurations. Lifting-surface theory improves on lifting-line theory by resolving the loading over the planform. Thickness panel methods extend the approach further by representing the actual surface shape of wings, fuselages, nacelles, pylons, and other aircraft components.

Summary & Closure

Lifting-line theory provides a foundational framework for analyzing the aerodynamic behavior of finite wings. By modeling the wing as a spanwise distribution of bound circulation and accounting for the induced downwash generated by the trailing vortex system, the theory allows estimation of essential aerodynamic quantities, such as the lift distribution and induced drag. Although the theory relies on simplifying assumptions, including high aspect ratio, small angles of attack, and inviscid, incompressible flow, it represents the essential physics governing three-dimensional wing aerodynamics with reasonable accuracy. It remains a core analytical tool in aircraft design because of its simplicity, analytical tractability, and predictive capability.

Lifting-surface theory, therefore, represents the natural extension of lifting-line theory to a fully two-dimensional lifting surface. By allowing the circulation to vary over both span and chord, these methods provide a more complete representation of the three-dimensional aerodynamic loading on finite wings. They can accommodate sweep, taper, and other geometric complexities. When augmented with source and doublet or vortex distributions to account for thickness and lifting effects, lifting-surface formulations form the basis of modern aerodynamic panel methods used in preliminary aircraft design and analysis. They occupy an intermediate position between classical analytical methods and full computational fluid dynamics, providing substantially improved fidelity while retaining the physical transparency of potential-flow theory.