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 carefully balance these elements to optimize lift generation and minimize induced drag, among other factors, thereby achieving 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 the two-dimensional airfoil theory, which assumes infinite span and aspect ratio, the lifting line theory accounts for the three-dimensional effects of finite wings, including the formation of trailing vortices and the creation of 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, marking a significant milestone in aeronautical research. Before Prandtl’s work, aerodynamic theories were primarily based on two-dimensional airfoil analysis, which assumed a wing of infinite span and ignored the effects of the wingtip vortices and the three-dimensional impact caused by the wake. Frederick Lanchester’s work on the roll-up of wingtip vortices, as illustrated in a figure from his book (Figure 2), was influential in establishing the importance of adequately modeling finite-wing aerodynamics. Furthermore, as aviation advanced, particularly during WWI, including the transition from biplanes to higher-performance monoplanes by the 1920s, the predictive limitations of existing aerodynamic theories quickly became evident. Designers needed a way to better predict and optimize the performance not just of airfoils but of 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^[2] 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 over the span of the wing minimized 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. Though more advanced techniques, such as computational fluid dynamics (CFD), are now widely used, the 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 the dynamics of vorticity in fluid mechanics and aerodynamics. These theorems are derived from the Euler equations and express how vortex lines and vortex tubes behave in an 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 length of the tube. 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 convects with 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 velocity fields 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. Notice that a line vortex in three dimensions is analogous to a point vortex in two dimensions, but instead of being a singular point, it extends along a curve. The circulation around any closed curve is equal to the sum of the strengths of the line vortices intersecting 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; the principles are formally embodied in Helmholtz’s theorems.

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 in the development of the 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 the development of 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 the lifting line theory and the development of the monoplane wing equation follows the work of Glauert. 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 mirror each other.

Apart from these restrictions, the planform (chord), pitch angle (including any twist angle), and zero-lift angle of the airfoil section may vary arbitrarily across 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 are formed by the integrated effect of the vorticity in the boundary layer or equivalent vortex sheet surrounding the airfoil’s surface, as in the manner used to develop the thin airfoil theory. According to the general theory of vortex flows and Helmholtz’s theorems, these bound vortices cannot terminate within the fluid but must be continued into the wake as free line vortices. These are referred to as 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 occurring when the circulation around each airfoil section is assumed to have a constant value across the wing’s span. In this special case, the bound vortex system may be represented as a single “lumped” line vortex with strength , which, to be consistent with the thin airfoil theory, can be positioned at the aerodynamic center along the 1/4-chord of the wing, as shown in Figure 9. The trailing vortices will consist of two line vortices, each with the same magnitude of strength but opposite sense, originating from the tips of the wing and extending downstream in the direction of the freestream flow. This representation is referred to as a “horseshoe vortex system.”

In practice, the vortices trailing from the wing tips will generally not be straight lines because of variations in the downwash at different distances 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, the vortex system of an airfoil is more complex because the circulation is not uniform across the span. Instead, it typically reaches its maximum value 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 yields a more accurate representation of 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 representation of the vortex system as the fundamental premise 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 from the perspective of the pressure difference between the upper and lower surfaces of the wing. Suppose the lift distribution across the span of a wing reaches a maximum at its centerline. In that case, there will be a significant decrease in pressure above the wing and an increase in pressure below it, as shown in Figure 11. These pressure differences decrease toward the tips of the wing as the lift distribution tapers toward zero.

As the streamlines pass above the wing and flow inward toward the center, while those below the wing flow outward, they leave the trailing edge and form a discontinuity surface, as shown in Figure 11. The trailing vortices from the wing 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 of vortices behind a wing, as observed in practice using flow visualization methods.

Nearer to the wing, the influence of the trailing vortex system can be modeled by assuming that 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 will vary across the span of a wing, being symmetrical about the centerline and decreasing 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 element of the span, as shown in Figure 12.

Therefore, in the aggregate, a sheet of trailing vortices will extend behind the entire wingspan because of the continuous change in along the wing, i.e., . The induced downwash at any point of the span will then be obtained as 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 latter result will then apply at every section, , along the span of the wing, and different values of will be obtained depending on the distribution, , and the proximity to the wing tips. It will be apparent from the previous discussion that the magnitude of the circulation gradient, i.e., , is usually much higher near the wing tips, so the trailed wake will have greater local strength there. This behavior is consistent with Lanchester’s interpretation of the wake generated behind a wing. The question now is how this downwash distribution can be calculated, 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 will, in general, be different 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 outcome means that the lift per unit span will be 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, in two-dimensional flow, the lift vector at each section changes orientation and is slightly tilted rearward relative to its original direction, thereby reducing the vertical component of lift at a given angle of attack. The consequence is that the lift vector is locally tilted rearward by the induced angle of attack, , so that a component of the lift now acts in the downstream direction. This component is the induced drag, , i.e.,

(15)  

This result shows that the induced drag is a direct consequence of the three-dimensional vortex system of the wing, which produces the downwash velocity and hence the induced angle of attack .

Development of the Monoplane Wing Equation

The preceding results now 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 connects the local circulation to the local chord, the local geometric angle of attack, the zero-lift angle of the section, 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. It will be apparent that as varies from to , then 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.

Notice that the series chosen for the circulation satisfies the condition that the circulation decreases to zero at the wing tips. As shown in Figure 14, if the wing geometry and loading are symmetrical about the centerline, then only odd integral values of will occur in the Fourier series.

It will be apparent that 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., , for which and are shown in Figure 14, represent a symmetric modification to the primary loading, i.e., a deviation from the elliptically distributed form.

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^[6] wing. 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 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 distribution of lift about the wing centerline. These modes can arise from aileron deflections, rolling motion, yawing motion, sideslip, or other asymmetries along the wing’s span. This situation leads to lift asymmetry, which influences rolling moments, possible yawing moments, and induced drag. For example, ailerons are control surfaces located near the wingtips, which are 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 sizing and placement of ailerons 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 and lift coefficient, the coefficient has a fixed value because

(58)  

Therefore, the induced drag will be minimized when all higher-order coefficients () in the series for circulation 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 elliptic 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, it gives a constant sectional lift coefficient along the span.
4. Elliptic loading, or a close approximation to it, can also be achieved with wings of non-elliptic planform 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 instead measured from a wing reference line, then in this derivation must be replaced by . 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, and the resulting circulation distribution must generally contain 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 non-dimensional 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 semi-span of the wing.

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 enables designers to tailor the angle-of-attack distribution, thereby 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. These factors must be carefully considered in wing design to balance aerodynamic efficiency, structural feasibility, and operational requirements.

Planform Effects

The theoretically best aerodynamic efficiency () and lowest induced drag are obtained with an elliptically loaded wing. For an untwisted wing with a constant airfoil section, this condition is obtained with a wing that is elliptical in planform shape. This wing shape gives an elliptical spanwise aerodynamic loading, i.e., lift per unit span, and uniform downwash over the wing, which, as previously mentioned, is theoretically the minimum induced drag condition. However, exact elliptic loading is rarely achieved in practical wing design 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; anything more than 0.1 would usually be classified as 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 shape of the wing, i.e., its planform, will affect the distribution of lift and the formation of the trailing vortex system, hence 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 approximate the ideal elliptical spanwise loading closely and minimize the value of . The actual effective aspect ratio of the wing can also be improved somewhat by paying attention to the shape of the wing tips, which can lower the values of if they are suitably shaped, e.g., with more rounded contours or perhaps with the addition of 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 allows designers to optimize the lift distribution and induced drag for specific flight conditions.

Interference Effects

The photograph in Figure 20 below illustrates the nature of the spanwise loading over a wing from the process of 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 is, however, that even minor deviations from the ideal elliptical form can result in higher induced drag from the wing. 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, which are for the same approximate value of total wing lift. Notice that the fuselage can produce significant deviations from the ideal elliptical form, perhaps increasing the value of to over 1.15. However, local effects along the wing, such as those caused by the aerodynamic interference effects produced by an engine, tend to have a minor impact on the value of .

Wing Sweep

The sweep of a wing has a significant effect on the aerodynamic characteristics of a wing, 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.

Sweep can be included approximately within the framework of lifting-line theory 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 freestream velocity normal to the leading edge, i.e., the 1/4-chord line in the lifting line theory. 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 itself should not be interpreted simply as being multiplied by . A more careful statement is that the section lift is governed mainly by the normal component of the flow, while 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 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 more accurately 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’s aspect ratio. At low Mach numbers, where , this 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 trailed vortex system generated as a consequence of lift generation. Although winglets increase the wing’s wetted surface area and give a small increment to overall skin friction and profile drag, they typically produce a net reduction in total airplane drag for 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 mid-span of the wing. While this effective-span winglet theory does not represent the spanwise variation in downwash, it provides a useful approximation of the total induced drag using the principles of the Biot-Savart law. It highlights the aerodynamic benefits of winglets through their influence on 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 starting position of the tip vortex, effectively moving it farther from more of the 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 , the effects of any dihedral angle being neglected. 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 is the dominant effect governing 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, its predicted drag ratio is applied 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. <p>No winglet:</p>

   (114)  

2. <p>Vertical winglet:</p>

   (115)  

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

   (116)  

Minimizing induced drag from the use of a winglet 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, the spanwise projection of the winglet 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 this simplified trend: increasing the cant angle and winglet length increases the predicted induced-drag reduction, but with diminishing practical benefit 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 instead written in terms of the vertical height of the winglet, then a different coefficient must be used because the vertical height is . In practice, empirical values of can be larger or smaller depending on the actual winglet geometry, toe angle, airfoil section, loading, viscous drag penalty, and the interaction of the winglet 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 means of reducing 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 tip vortices of a leader, 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 of 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 of 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 induced-drag change 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 by which winglets act 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 every leader contributes a geometric function of the same form as but 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 one horseshoe vortex system, so a follower experiences upwash or downwash resulting from the combined action of all aircraft in front. A typical variation of is shown in Figure 26. Negative values are upwash regions that reduce induced drag for the follower.

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

(132)  

with the same coefficient

(133)  

As before, for similar aircraft flying at similar lift coefficients, the value of 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, in Figure 27, shows five F-15Cs in formation during a ferry flight, a method used to conserve fuel over long distances. Because the leader 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, in which the wing is represented by a single bound vortex filament located approximately at the quarter-chord and extending along the span. While analytically elegant and computationally efficient, lifting-line theory does not account for chordwise loading variations or the presence of wing thickness. This limitation becomes important for swept, tapered, low-aspect-ratio, or non-planar wings, where the aerodynamic loading cannot be represented adequately by a single spanwise circulation distribution.

Lifting-surface theory may be viewed as 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, in lifting-surface theory the circulation is allowed to vary over both the span and the chord, i.e., in the form . 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 improvement in fidelity comes from resolving how aerodynamic loading varies 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 wings, delta wings, tapered wings, 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, not just along a spanwise line.

Multhopp’s method was a more systematic subsonic lifting-surface method for calculating the lift distribution of wings. The continuous lifting surface was represented by a finite set of unknown loading values, and the downwash condition was imposed at selected pivotal points over 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 not inherently a vortex-lattice method and is not inherently numerical. However, exact or closed-form solutions are available only for special geometries or special 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 located at is obtained from 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 associated with a unit-strength vortex element located 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 and thereby relates 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, or 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 to be thin and is 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 distribution of circulation 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 from all of the vortex elements. Compared with lifting-line theory, numerical lifting-surface methods give a better representation of chordwise loading, aerodynamic-center movement, pitching moment, and the effects of planform geometry. For example, a swept or low-aspect-ratio wing may have a strongly nonuniform chordwise loading, so the local pressure forces and pitching moments cannot be recovered from a single spanwise lifting line. 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 account primarily for the lifting effects of thin surfaces, while 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 and 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.