**1. Introduction**

In Prandtl’s classical lifting line theory, a finite wing is replaced by a single spanwise lifting line along which the circulation (*y*) varies; coordinate *y *denotes the spanwise distance from the center of the wing. At some arbitrary distance *y _{0}* along the lifting line, the trailing vortex system that accompanies the lifting line will induce an infinitesimal velocity

*dw*given by

The total velocity *w *induced by the entire trailing vortex sheet is obtained by integrating (1) over the entire wingspan, that is, from –*b* to *b:*

This induced velocity is generally in the downward direction and is termed *downwash*. Most importantly, downwash has the effect of tilting the resultant wind through an angle such that (Figure 1)

where we have used the small-angle approximation because the downwash *w *is generally much less than the freestream airspeed *U*. Per equation (3), the downwash reduces the effective incidence so that for the same lift as an equivalent airfoil at incidence , an incidence of = + is required at that section of the airfoil. This modification of angle of attack is a *kinematic *consequence of downwash, but, as we’ll see in a moment, downwash also has crucial *dynamic *(i.e., force-related) implications for wing aerodynamics.

*Figure 1. Angles and velocities on an airfoil section. ** = Geometric angle of attack; * *= Induced angle of attack; ** *= *Effective angle of attack. U *= *Freestream velocity; w = Downwash; U _{R} = Resultant velocity.*

**2. Elliptic circulation distribution and downwash**

In order to ascertain an induced angle of attack or a lift distribution on a wing, we must know the form of the circulation distribution (*y*). In theoretical terms, the most important such distribution is given by

where is the circulation at midspan. Since the circulation described by (4) varies elliptically along the span, it is termed an *elliptic circulation distribution* (Fig. 2).

*Figure 2. Elliptically loaded wing.*

Obtaining the downwash associated with an elliptic distribution is straightforward. The first step is to differentiate (4) with respect to *y*, giving

Then, we substitute into (2) and set up the integral

Adding and subtracting *y _{0}* to the numerator, we can restate (6) as

The first integral in brackets is a standard result that equals , while the second may be replaced with some symbol *K, *so that

Now, as this is a symmetric flight case, the vorticity shed is the same from each side of the wing; therefore, the downwash induced at some spanwise coordinate *y _{0}* in one half of the wing must be the same as that in coordinate –

*y*in the opposite half of the wing. In mathematical terms,

_{0}The only way to satisfy this equality is to have *K = *0, so that

This important result implies that the downwash on an elliptically loaded wing is constant along the wing span.

**3. Lift and drag in an elliptic circulation distribution**

According to the Kutta-Joukowski theorem, the contribution to lift *L’ *at a local airfoil section located at *y _{0}* is given by

where is the freestream air density and is the freestream velocity. Replacing (*y _{0}*) with (4) indicates that the lift distribution of an elliptically loaded wing is also elliptic in nature:

To find the *total *lift afforded by the circulation distribution, we integrate (12) along the wingspan:

Equation (13) gives the total lift as a function of midspan circulation and wingspan. However, we know that the lift is also given by the classic relationship

where *S *= (*span ** chord*) is the projected area of the wing in the direction normal to the velocity *U*. Combining (13) and (14) and solving for , we get

Now, the induced drag caused by the downwash is given by

so that, replacing *w *with (10) and with (4), we write

However, we can replace midspan circulation with equation (15), giving

But drag can also be expressed as

where *C _{D,I}* is the induced drag coefficient. Equating (17) to (19) and manipulating, we find that

where *A _{R}* is the aspect ratio of the wing, namely

As can be seen from (20), the variation of *C _{L}* with

*C*for an elliptically loaded airfoil is a parabola, provided the only source of drag is the induced velocity – that is, the downwash. The plot of

_{D,I}*C*versus

_{L}*C*is called the

_{D,I}*polar curve*of the airfoil (Fig. 3).

Thus, the induced drag coefficient is directly proportional to the square of the lift coefficient. The dependence of induced drag on lift is not surprising. In phenomenological terms, induced drag is a consequence of the pressure of the wing-tip vortices, which in turn are produced by the difference in pressure between the lower and upper wing surfaces. But the lift is produced by this same pressure difference, hence induced drag is closely related to the production of lift on a finite wing – so much so, in fact, that induced drag is also known as *drag due to lift*. As Anderson (1984) puts it, an airplane cannot generate lift for free; the induced drag is the price for generation of lift.

In practice, in addition to induced drag there is also profile drag due to skin friction and wake. These phenomena are accounted for by adding a constant *C _{D,0}* to the right-hand side of the drag polar, giving

Since profile drag coefficient *C _{D,0}* is independent of

*C*, it is also known as the

_{L}*zero-lift drag coefficient*. A more rigorous derivation, which is beyond the scope of this post, leads to the modified drag polar

where *e *is a parameter known as Oswald’s wing efficiency and is a nonnegative parameter that equals zero only in the case of elliptic loading.

*Figure 3. Drag polar. *

For a wing under elliptic circulation distribution, in straight level flight, what is the condition to minimize drag? We can answer this question using the drag polar (22). Note first that the drag *D *can be recast as

But at level flight, lift *L *equals the weight *W*, so that

Thus, for a wing of given weight *W*, the minimum drag occurs when the drag-to-lift ratio is a minimum (or, equivalently, the lift-to-drag ratio is a maximum). We proceed to restate *D/L *as

where *k *is a coefficient that condenses all factors that multiply *C _{L}* in the drag polar (22). For minimum drag, we differentiate (26) with respect to

*C*and set the ensuing expression to zero:

_{L}This derivative is straightforward to obtain. We can speed things up using Mathematica:

Clearly,

Setting this to zero yields

which means that, for minimum drag in an elliptical load distribution, the zero-lift drag coefficient must be equal to the lift-dependent contribution to drag. (The reader is invited to differentiate (28) a second time, substitute from (29), and check that the second derivative at the point in question is positive, which indicates that the point we’ve discussed is indeed a local minimum, not a local maximum.)

**Solved example 1**

A wing with elliptic loading has span equal to 18 m, planform area 56 m^{2}, and moves in level flight at 792 km/h at an altitude where air density equals 0.76 kg/m^{3}. Assuming the induced drag on the wing equals 3400 N, determine:

*(a)* The lift coefficient.

*(b)* The downwash velocity.

*(c)* The wing loading.

*Part a. *Before anything else, we compute the aspect ratio *A _{R}*:

With a flight velocity *V *= 792 km/h = 220 m/s, we proceed to calculate the induced drag coefficient:

Then, the corresponding lift coefficient is determined as

*Part b. *The downwash for an elliptically loaded wing is constant and given by (10):

where

so that

*Part c. *In level flight, lift equals weight, or

Lastly, the wing loading is

**Solved example 2**

Obtain an expression for the downward induced velocity behind a wing of span 2*b *at a point *y *from the center of the span, the circulation around the wing at any point *y *being . If the circulation is parabolic, that is,

calculate the value of the induced velocity *w *at midspan, and compare this value with that obtained when the same lift is distributed elliptically.

*Solution. *The derivative of the given span circulation we were given is

so that, substituting in (2),

At midspan this reduces to

It is pertinent to ask how this compares with an elliptically loaded wing with the same lift. Note first that the lift for the parabolic circulation distribution is such that

In contrast, the lift for an elliptically loaded wing is given by eq. (13), namely

For the same lift in both cases, we must have

Further, noting that

and substituting from (30) and (31), we get

which implies that the downwash associated with an elliptically loaded wing is 2/3 of the downwash of a parabolically loaded wing.

**References**

- ANDERSON, J.D. Jr. (1984).
*Fundamentals of Aerodynamics*. New York: McGraw-Hill. - HOUGHTON, E.L., CARPENTER, P.W., COLLICOTT, S.H. and VALENTINE, D.T. (2017).
*Aerodynamics for Engineering Students*. 7th edition. Oxford: Butterworth-Heinemann. - RATHAKRISHNAN, E. (2013).
*Theoretical Aerodynamics*. Hoboken: John Wiley and Sons.