Category Fundamentals of Aerodynamics

Volume Integrals

Consider a volume V in space. Let p be a scalar field in this space. The volume integral over the volume V of the quantity p is written as

image99p dV = volume integral of a scalar p over the
volume V (the result is a scalar)

image100

Figure 2.9 Sketch for surface integrals. The

three-dimensional surface area S is bounded by the closed curve C.

 

image101

Figure 2.10 Volume V enclosed by the closed surface S.

 

Let A be a vector field in space. The volume integral over the volume V of the quantity A is written as

 

A dV = volume integral of a vector A over the volume V (the result is a vector)

 

2.2.11 Relations Between Line, Surface, and Volume Integrals

 

Consider again the open area S bounded by the closed curve C, as shown in Figure 2.9. Let A be a vector field. The line integral of A over C is related to the surface integral of A over S by Stokes ’ theorem:

 

[2.25]

 

Consider again the volume V enclosed by the closed surface S, as shown in Figure 2.10. The surface and volume integrals of the vector field A are related through the divergence theorem:

 

image104

s

 

(V. A) dV

 

[2.26]

 

v

 

If p represents a scalar field, a vector relationship analogous to Equation (2.26) is given by the gradient theorem:

 

[2.27]

 

VpdV

 

image102image103image105

2.2.12 Summary

This section has provided a concise review of those elements of vector analysis which we will use as tools in our subsequent discussions. Make certain to review these tools until you feel comfortable with them, especially the relations in boxes.

Infinity Boundary Conditions

Far away from the body (toward infinity), in all directions, the flow approaches the uniform freestream conditions. Let be aligned with the x direction as shown in Figure 3.18. Hence, at infinity,

Подпись: [3.47a]З Ф дф

и = — = — = Lx дх 3 у

Infinity Boundary Conditions

Подпись:Boundary conditions at infinity and on a body; inviscid flow.

Подпись: [3.47b]дф дф ду дх

Equations (3.47a and b) are the boundary conditions on velocity at infinity. They apply at an infinite distance from the body in all directions, above and below, and to the left and right of the body, as indicated in Figure 3.18.

Prandtl’s Classical Lifting-Line Theory

The first practical theory for predicting the aerodynamic properties of a finite wing was developed by Ludwig Prandtl and his colleagues at Gottingen, Germany, during the period 1911-1918, spanning World War I. The utility of Prandtl’s theory is so great

that it is still in use today for preliminary calculations of finite-wing characteristics. The purpose of this section is to describe Prandtl’s theory and to lay the groundwork for the modern numerical methods described in subsequent sections.

Prandtl reasoned as follows. A vortex filament of strength Г that is somehow bound to a fixed location in a flow—a so-called bound vortex—will experience a force L = Poo ‘oc Г from the Kutta-Joukowski theorem. This bound vortex is in contrast to a free vortex, which moves with the same fluid elements throughout a flow. Therefore, let us replace a finite wing of span b with a bound vortex, extending from у = —b/2 to у = b/2, as sketched in Figure 5.10. However, due to Helmholtz’s theorem, a vortex filament cannot end in the fluid. Therefore, assume the vortex filament continues as two free vortices trailing downstream from the wing tips to infinity, as also shown in Figure 5.10. This vortex (the bound plus the two free) is in the shape of a horseshoe, and therefore is called a horseshoe vortex.

Подпись: w(y) — -
Prandtl’s Classical Lifting-Line Theory Подпись: [5.13]

A single horseshoe vortex is shown in Figure 5.11. Consider the down wash w induced along the bound vortex from — b/2 to b/2 by the horseshoe vortex. Examining Figure 5.11, we see that the bound vortex induces no velocity along itself; however, the two trailing vortices both contribute to the induced velocity along the bound vortex, and both contributions are in the downward direction. Consistent with the xyz coordinate system in Figure 5.11, such a downward velocity is negative; that is, w (which is in the г direction) is a negative value when directed downward and a positive value when directed upward. If the origin is taken at the center of the bound vortex, then the velocity at any point у along the bound vortex induced by the trailing semi-infinite vortices is, from Equation (5.11),

In Equation (5.12), the first term on the right-hand side is the contribution from the left trailing vortex (trailing from —b/2), and the second term is the contribution from the right trailing vortex (trailing from b/2). Equation (5.12) reduces to

image415[5.13]

Подпись: LПодпись: tПодпись: Free-trailing vortexПодпись: Finite wingПодпись: Horseshoe vortexПодпись:image4164 Free-trailing vortex

image417

Figure 5.1 1 Downwash distribution along the у axis for a single horseshoe vertex.

This variation of w(y) is sketched in Figure 5.11. Note that w approaches —oo as у approaches —b/2 or b/2.

The downwash distribution due to the single horseshoe vortex shown in Figure 5.11 does not realistically simulate that of a finite wing; the downwash approaching an infinite value at the tips is especially disconcerting. During the early evolution of finite-wing theory, this problem perplexed Prandtl and his colleagues. After several years of effort, a resolution of this problem was obtained which, in hindsight, was simple and straightforward. Instead of representing the wing by a single horseshoe vortex, let us superimpose a large number of horseshoe vortices, each with a different length of the bound vortex, but with all the bound vortices coincident along a single line, called the lifting line. This concept is illustrated in Figure 5.12, where only three horseshoe vortices are shown for the sake of clarity. In Figure 5.12, a horseshoe vortex of strength <іГі is shown, where the bound vortex spans the entire wing from —b/2 to b/2 (from point A to point F). Superimposed on this is a second horseshoe vortex of strength (ІГ2, where its bound vortex spans only part of the wing, from point В to point E. Finally, superimposed on this is a third horseshoe vortex of strength tf Г3, where its bound vortex spans only the part of the wing from point C to point D. Asa result, the circulation varies along the line of bound vortices—the lifting line defined above. Along A В and EF, where only one vortex is present, the circulation is dV . However, along BC and DE, where two vortices are superimposed, the circulation is the sum of their strengths dV]+dr2- Along CD, three vortices are superimposed, and hence the circulation is dVi + dT2 + dT3. This variation of Г along the lifting line is denoted by the vertical bars in Figure 5.12. Also, note from Figure 5.12 that we now have a series of trailing vortices distributed over the span, rather than just two vortices trailing downstream of the tips as shown in Figure 5.11. The series of

oo

 

oo

 

Prandtl’s Classical Lifting-Line Theory

Figure 5.12 Superposition of a finite number of horseshoe vortices along the lifting line.

 

image418

trailing vortices in Figure 5.12 represents pairs of vortices, each pair associated with a given horseshoe vortex. Note that the strength of each trailing vortex is equal to the change in circulation along the lifting line.

Let us extrapolate Figure 5.12 to the case where an infinite number of horseshoe vortices are superimposed along the lifting line, each with a vanishingly small strength dV. This case is illustrated in Figure 5.13. Note that the vertical bars in Figure 5.12 have now become a continuous distribution of Г(у) along the lifting line in Figure 5.13. The value of the circulation at the origin is Г0. Also, note that the finite number of trailing vortices in Figure 5.12 have become a continuous vortex sheet trailing downstream of the lifting line in Figure 5.13. This vortex sheet is parallel to the direction of Too. The total strength of the sheet integrated across the span of the wing is zero, because it consists of pairs of trailing vortices of equal strength but in opposite directions.

Let us single out an infinitesimally small segment of the lifting line dy located at the coordinate у as shown in Figure 5.13. The circulation at у is Г(у), and the change in circulation over the segment dy is dV = (dГ/dy) dy. In turn, the strength of the trailing vortex at у must equal the change in circulation dГ along the lifting line; this is simply an extrapolation of our result obtained for the strength of the finite trailing vortices in Figure 5.12. Consider more closely the trailing vortex of strength dГ which intersects the lifting line at coordinate y, as shown in Figure 5.13. Also consider the arbitrary location yo along the lifting line. Any segment of the trailing vortex dx will induce a velocity at yo with a magnitude and direction given by the Biot-Savart law, Equation (5.5). In turn, the velocity dw at yo induced by the entire semi-infinite trailing vortex located at у is given by Equation (5.11), which in terms of the picture given in Figure 5.13 yields

Prandtl’s Classical Lifting-Line Theory

[5.14]

 

image419

line.

 

The minus sign in Equation (5.14) is needed for consistency with the picture shown in Figure 5.13; for the trailing vortex shown, the direction of dw at yo is upward and hence is a positive value, whereas Г is decreasing in the у direction, making dV/dy a negative quantity. The minus sign in Equation (5.14) makes the positive dw consistent with the negative dV/dy.

The total velocity w induced at y(l by the entire trailing vortex sheet is the summa­tion of Equation (5.14) over all the vortex filaments, that is, the integral of Equation

(5.14)

Подпись: w(yo) = Подпись: [5.15]

from – b/2tob/2:

Equation (5.15) is important in that it gives the value of the downwash at y(l due to all the trailing vortices. (Keep in mind that although we label w as downwash, w is treated as positive in the upward direction in order to be consistent with the normal convention in an xyz rectangular coordinate system.)

Pause for a moment and assess the status of our discussion so far. We have replaced the finite wing with the model of a lifting line along which the circulation Г (у) varies continuously, as shown in Figure 5.13. In turn, we have obtained an ex­pression for the downwash along the lifting line, given by Equation (5.15). However, our central problem still remains to be solved; that is, we want to calculate T(y) for a given finite wing, along with its corresponding total lift and induced drag. Therefore, we must press on.

Return to Figure 5.4, which shows the local airfoil section of a finite wing. Assume this section is located at the arbitrary spanwise station yo. From Figure 5.4,

Подпись:Prandtl’s Classical Lifting-Line Theory
the induced angle of attack a, is given by

-w(yp)

Prandtl’s Classical Lifting-Line Theory Prandtl’s Classical Lifting-Line Theory

[Note in Figure 5.4 that w is downward, and hence is a negative quantity. Since a, in Figure 5.4 is positive, the negative sign in Equation (5.16) is necessary for consistency.] Generally, w is much smaller than VJ, and hence a, is a small angle, on the order of a few degrees at most. For small angles, Equation (5.16) yields

that is, an expression for the induced angle of attack in terms of the circulation distribution Г(у) along the wing.

Consider again the effective angle of attack aeff, as shown in Figure 5.4. As explained in Section 5.1, oy-n is the angle of attack actually seen by the local airfoil section. Since the downwash varies across the span, then оуц is also variable; ауп = ateff (jo)• The lift coefficient for the airfoil section located at у = y0 is

Ci = ao[aeff(>’o) — «/.= ol = 27r[aeff(yo) — аТ=о] [5.19]

Подпись: і ■i In Equation (5.19), the local section lift slope ao has been replaced by the thin airfoil theoretical value of 2n-(rad-1). Also, for a wing with aerodynamic twist, the angle of zero lift ai=o in Equation (5.19) varies with vy. If there is no aerodynamic twist, aL=0 is constant across the span. In any event, aL=0 is a known property of the local airfoil sections. From the definition of lift coefficient and from the Kutta-Joukowski theorem, we have, for the local airfoil section located at yo.

L’ = {pccV£,c(yo)ci = РооЕооГ(уо) [5.20]

From Equation (5.20), we obtain

Подпись:_ 2Г (yo)

C> V’coc(y0)

Подпись: Г(уо) л-Еосс(у0) Подпись: <*L= 0 Подпись: [5.22]

Substituting Equation (5.21) into (5.19) and solving for оуц, we have

The above results come into focus if we refer to Equation (5.1):

Подпись: [5.1]«eff = a – a.

Prandtl’s Classical Lifting-Line Theory

Substituting Equations (5.18) and (5.22) into (5.1), we obtain

 

Prandtl’s Classical Lifting-Line Theory

[5.33]

 

the fundamental equation of Prandtl’s lifting-line theory, it simply states that the geometric angle of attack is equal to the sum of the effective angle plus the induced angle of attack. In Equation (5.23), aeff is expressed in terms of Г, and a, is expressed in terms of an integral containing dT/dy. Hence, Equation (5.23) is an integro – differential equation, in which the only unknown is Г; all the other quantities, a, c, Voo, and a£=o, are known for a finite wing of given design at a given geometric angle of attack in a freestream with given velocity. Thus, a solution of Equation (5.23) yields Г = Г(уо), where y0 ranges along the span from —h/2 to h/2.

The solution Г = Г(Ч’о) obtained from Equation (5.23) gives us the three main aerodynamic characteristics of a finite wing, as follows:

1. The lift distribution is obtained from the Kutta-Joukowski theorem:

 

Prandtl’s Classical Lifting-Line Theory

[5.34]

 

fb/ 2

L= L'(y) dy J-b/2

fb/2

L = PooVoo r(y)dy

J-b/2

 

[5.35]

 

or

 

(Note that we have dropped the subscript on y, for simplicity.) The lift coefficient follows immediately from Equation (5.25):

 

Prandtl’s Classical Lifting-Line Theory

CL

 

[5.36]

 

image424

An important geometric property of a finite wing is the aspect ratio, denoted by AR and defined as

Подпись: a; Подпись: CL л-AR Подпись: [5.42]

Hence, Equation (5.41) becomes

Equation (5.42) is a useful expression for the induced angle of attack, as shown below.

Prandtl’s Classical Lifting-Line Theory Подпись: 2a І Fab VXS 2 Подпись: sin2 9 d9 Подпись: TtajF{)b 2VXS Подпись: [5.42a]
image425

The induced drag coefficient is obtained from Equation (5.30), noting that a, is constant:

Substituting Equations (5.40) and (5.42) into (5.42a), we obtain

Подпись: or Подпись: CDJ Подпись: C2 ж AR Подпись: [5.43]

nb / CL 2VqcSCl 2 Voo S tr AR ) bn

Equation (5.43) is an important result. It states that the induced drag coefficient is directly proportional to the square of the lift coefficient. The dependence of induced drag on the lift is not surprising, for the following reason. In Section 5.1, we saw that induced drag is a consequence of the presence of the wing-tip vortices, which in turn are produced by the difference in pressure between the lower and upper wing surfaces. The lift is produced by this same pressure difference. Hence, induced drag is intimately related to the production of lift on a finite wing; indeed, induced drag is frequently called the drag due to lift. Equation (5.43) dramatically illustrates this point. Clearly, an airplane cannot generate lift for free; the induced drag is the price for the generation of lift. The power required from an aircraft engine to overcome the induced drag is simply the power required to generate the lift of the aircraft. Also, note that because Co. і cx C, the induced drag coefficient increases rapidly as CL increases and becomes a substantial part of the total drag coefficient when Cl is high (e. g., when the airplane is flying slowly such as on takeoff or landing). Even at relatively high cruising speeds, induced drag is typically 25 percent of the total drag.

Another important aspect of induced drag is evident in Equation (5.43); that is, Co. і is inversely proportional to aspect ratio. Hence, to reduce the induced drag, we want a finite wing with the highest possible aspect ratio. Wings with high and low aspect ratios are sketched in Figure 5.14. Unfortunately, the design of very high aspect ratio wings with sufficient structural strength is difficult. Therefore, the aspect ratio of a conventional aircraft is a compromise between conflicting aerodynamic and structural requirements. It is interesting to note that the aspect ratio of the 1903 Wright

AR = b2/S

 

c

 

High AR (low induced drag)

 

b

 

Low AR (high induced drag)

 

Figure 5.14 Schematic of high-and low-aspect-ratio wings.

 

image426

Prandtl’s Classical Lifting-Line Theory

Flyer was 6 and that today the aspect ratios of conventional subsonic aircraft range typically from 6 to 8. (Exceptions are the Lockheed U-2 high-altitude reconnaissance aircraft with AR = 14.3 and sailplanes with aspect ratios in the 10 to 22 range.)

Another property of the elliptical lift distribution is as follows. Consider a wing with no geometric twist (i. e., a is constant along the span) and no aerodynamic twist (i. e., a/ =o is constant along the span). From Equation (5.42), we have seen that o’, is constant along the span. Hence, — a — cti is also constant along the span. Since the local section lift coefficient q is given by

ci = a0(aeS – aL=o)

then assuming that ao is the same for each section (ag = 2л from thin airfoil theory), ci must be constant along the span. The lift per unit span is given by

Подпись: [5.44]L'(y) = qcaCCi

Prandtl’s Classical Lifting-Line Theory Подпись: [5.45]

Solving Equation (5.44) for the chord, we have

In Equation (5.45), q00 and c; are constant along the span. However, L'(y) varies elliptically along the span. Thus, Equation (5.45) dictates that for such an elliptic lift distribution, the chord must vary elliptically along the spam, that is, for the conditions given above, the wing planform is elliptical.

The related characteristics—the elliptic lift distribution, the elliptic planform, and the constant downwash—are sketched in Figure 5.15. Although an elliptical lift distribution may appear to be a restricted, isolated case, in reality it gives a reasonable approximation for the induced drag coefficient for an arbitrary finite wing. The form

Подпись: figure 5.1 5 Illustration of the related quantities: an elliptic lift distribution, elliptic planform, and constant downwash.

of Сил given by Equation (5.43) is only slightly modified for the general case. Let us now consider the case of a finite wing with a general lift distribution.

Thoughts

The term “aerodynamics” is generally used for problems arising from flight and other topics involving the flow of air.

Ludwig Prandtl, 1949

Aerodynamics: The dynamics of gases, especially atmospheric interactions with moving objects.

The American Heritage Dictionary of the English Language, 1969

1.1 Importance of Aerodynamics: Historical Examples

On August 8, 1588, the waters of the English Channel churned with the gyrations of hundreds of warships. The great Spanish Armada had arrived to carry out an invasion of Elizabethan England and was met head-on by the English fleet under the command of Sir Francis Drake. The Spanish ships were large and heavy; they were packed with soldiers and carried formidable cannons that fired 50 lb round shot that could devastate any ship of that era. In contrast, the English ships were smaller and lighter; they carried no soldiers and were armed with lighter, shorter-range cannons.

The balance of power in Europe hinged on the outcome of this naval encounter. King Philip II of Catholic Spain was attempting to squash Protestant England’s rising influence in the political and religious affairs of Europe; in turn, Queen Elizabeth I was attempting to defend the very existence of England as a sovereign state. In fact, on that crucial day in 1588, when the English floated six fire ships into the Spanish formation and then drove headlong into the ensuing confusion, the future history of Europe was in the balance. In the final outcome, the heavier, sluggish, Spanish ships were no match for the faster, more maneuverable, English craft, and by that evening the Spanish Armada lay in disarray, no longer a threat to England. This naval battle is of particular importance because it was the first in history to be fought by ships on both sides powered completely by sail (in contrast to earlier combinations of oars and sail), and it taught the world that political power was going to be synonymous with naval power. In turn, naval power was going to depend greatly on the speed and maneuverability of ships. To increase the speed of a ship, it is important to reduce the resistance created by the water flow around the ship’s hull. Suddenly, the drag on ship hulls became an engineering problem of great interest, thus giving impetus to the study of fluid mechanics.

This impetus hit its stride almost a century later, when, in 1687, Isaac Newton (1642-1727) published his famous Principia, in which the entire second book was devoted to fluid mechanics. Newton encountered the same difficulty as others before him, namely, that the analysis of fluid flow is conceptually more difficult than the dynamics of solid bodies. A solid body is usually geometrically well defined, and its motion is therefore relatively easy to describe. On the other hand, a fluid is a “squishy” substance, and in Newton’s time it was difficult to decide even how to qualitatively model its motion, let alone obtain quantitative relationships. Newton considered a fluid flow as a uniform, rectilinear stream of particles, much like a cloud of pellets from a shotgun blast. As sketched in Figure 1.1, Newton assumed that upon striking a surface inclined at an angle в to the stream, the particles would

Rectilinear stream of discrete particles

 

Upon impacting the body, the particles give up their momentum normal to the surface, and travel downstream along the surface.

 

Figure 1«1 Isaac Newton’s model of fluid flow in the year 1687.

This model was widely adopted in the seventeenth and eighteenth centuries but was later found to be conceptually inaccurate for most fluid flows.

 

image6

transfer their normal momentum to the surface but their tangential momentum would be preserved. Hence, after collision with the surface, the particles would then move along the surface. This led to an expression for the hydrodynamic force on the surface which varies as sin2 9. This is Newton’s famous sine-squared law (described in detail in Chapter 14). Although its accuracy left much to be desired, its simplicity led to wide application in naval architecture. Later, in 1777, a series of experiments was carried out by Jean LeRond d’Alembert (1717-1783), under the support of the French government, in order to measure the resistance of ships in canals. The results showed that “the rule that for oblique planes resistance varies with the sine square of the angle of incidence holds good only for angles between 50 and 90° and must be abandoned for lesser angles.” Also, in 1781, Leonhard Euler(1707-1783)pointed out the physical inconsistency of Newton’s model (Figure 1.1) consisting of a rectilinear stream of particles impacting without warning on a surface. In contrast to this model, Euler noted that the fluid moving toward a body “before reaching the latter, bends its direction and its velocity so that when it reaches the body it flows past it along the surface, and exercises no other force on the body except the pressure corresponding to the single points of contact.” Euler went on to present a formula for resistance which attempted to take into account the shear stress distribution along the surface, as well as the pressure distribution. This expression became proportional to sin2 в for large incidence angles, whereas it was proportional to sin 9 at small incidence angles. Euler noted that such a variation was in reasonable agreement with the ship-hull experiments carried out by d’Alembert.

This early work in fluid dynamics has now been superseded by modern concepts and techniques. (However, amazingly enough, Newton’s sine-squared law has found new application in very high-speed aerodynamics, to be discussed in Chapter 14.) The major point here is that the rapid rise in the importance of naval architecture after the sixteenth century made fluid dynamics an important science, occupying the minds of Newton, d’Alembert, and Euler, among many others. Today, the modem ideas of fluid dynamics, presented in this book, are still driven in part by the importance of reducing hull drag on ships.

Consider a second historical example. The scene shifts to Kill Devil Hills, 4 mi south of Kitty Hawk, North Carolina. It is summer of 1901, and Wilbur and Orville Wright are struggling with their second major glider design, the first being a stunning failure the previous year. The airfoil shape and wing design of their glider are based on aerodynamic data published in the 1890s by the great German aviation pioneer Otto Lilienthal (1848-1896) and by Samuel Pierpont Langley (1934-1906), secretary of the Smithsonian Institution—the most prestigious scientific position in the United States at that time. Because their first glider in 1900 produced no meaningful lift, the Wright brothers have increased the wing area from 165 to 290 ft2 and have increased the wing camber (a measure of the airfoil curvature—the larger the camber, the more “arched” is the thin airfoil shape) by almost a factor of 2. But something is still wrong. In Wilbur’s words, the glider’s “lifting capacity seemed scarcely one-third of the calculated amount.” Frustration sets in. The glider is not performing even close to their expectations, although it is designed on the basis of the best available aerodynamic data. On August 20, the Wright brothers despairingly pack themselves

aboard a train going back to Dayton, Ohio. On the ride back, Wilbur mutters that “nobody will fly for a thousand years.” However, one of the hallmarks of the Wrights is perseverance, and within weeks of returning to Dayton, they decide on a complete departure from their previous approach. Wilbur later wrote that “having set out with absolute faith in the existing scientific data, we were driven to doubt one thing after another, until finally after two years of experiment, we cast it all aside, and decided to rely entirely upon our own investigations.” Since their 1901 glider was of poor aerodynamic design, the Wrights set about determining what constitutes good aerodynamic design. In the fall of 1901, they design and build a 6 ft long, 16 in square wind tunnel powered by a two-bladed fan connected to a gasoline engine. An original photograph of the Wrights’ tunnel in their Dayton bicycle shop is shown in Figure 1.2a. In this wind tunnel they test over 200 different wing and airfoil shapes, including flat plates, curved plates, rounded leading edges, rectangular and curved planforms, and various monoplane and multiplane configurations. A sample of their test models is shown in Figure 1.2b. The aerodynamic data is taken logically and carefully. It shows a major departure from the existing “state-of-the-art” data. Armed with their new aerodynamic information, the Wrights design a new glider in the spring of 1902. The airfoil is much more efficient; the camber is reduced considerably, and the location of the maximum rise of the airfoil is moved closer to the front of the wing. The most obvious change, however, is that the ratio of the length of the wing (wingspan) to the distance from the front to the rear of the airfoil (chord length) is increased from 3 to 6. The success of this glider during the summer and fall of 1902 is astounding; Orville and Wilbur accumulate over a thousand flights during this period. In contrast to the previous year, the Wrights return to Dayton flushed with success and devote all their subsequent efforts to powered flight. The rest is history.

The major point here is that good aerodynamics was vital to the ultimate success of the Wright brothers and, of course, to all subsequent successful airplane designs up to the present day. The importance of aerodynamics to successful manned flight goes without saying, and a major thrust of this book is to present the aerodynamic fundamentals that govern such flight.

Consider a third historical example of the importance of aerodynamics, this time as it relates to rockets and space flight. High-speed, supersonic flight had become a dominant feature of aerodynamics by the end of World War II. By this time, aerody – namicists appreciated the advantages of using slender, pointed body shapes to reduce the drag of supersonic vehicles. The more pointed and slender the body, the weaker the shock wave attached to the nose, and hence the smaller the wave drag. Conse­quently, the German V-2 rocket used during the last stages of World War II had a pointed nose, and all short-range rocket vehicles flown during the next decade fol­lowed suit. Then, in 1953, the first hydrogen bomb was exploded by the United States. This immediately spurred the development of long-range intercontinental ballistic missiles (ICBMs) to deliver such bombs. These vehicles were designed to fly outside the region of the earth’s atmosphere for distances of 5000 mi or more and to reenter the atmosphere at suborbital speeds of from 20,000 to 22,000 ft/s. At such high velocities, the aerodynamic heating of the reentry vehicle becomes severe, and this heating problem dominated the minds of high-speed aerodynamicists. Their first

image7

(a)

image8

(b)

Figure 1 Л (a) Wind tunnel designed, built, and used by the Wright brothers in Dayton, Ohio, during 1901-1902. (b) Wing models tested by the Wright brothers in their wind tunnel during 1901-1902.

thinking was conventional—a sharp-pointed, slender reentry body. Efforts to min­imize aerodynamic heating centered on the maintenance of laminar boundary layer flow on the vehicle’s surface; such laminar flow produces far less heating than turbu­lent flow (discussed in Chapters 15 and 19). However, nature much prefers turbulent flow, and reentry vehicles are no exception. Therefore, the pointed-nose reentry body was doomed to failure because it would burn up in the atmosphere before reaching the earth’s surface.

However, in 1951, one of those major breakthroughs that come very infrequently in engineering was created by H. Julian Allen at the NACA (National Advisory Com­mittee for Aeronautics) Ames Aeronautical Laboratory—he introduced the concept of the blunt reentry body. His thinking was paced by the following concepts. At the beginning of reentry, near the outer edge of the atmosphere, the vehicle has a large amount of kinetic energy due to its high velocity and a large amount of potential en­ergy due to its high altitude. However, by the time the vehicle reaches the surface of the earth, its velocity is relatively small and its altitude is zero; hence, it has virtually no kinetic or potential energy. Where has all the energy gone? The answer is that it has gone into (1) heating the body and (2) heating the airflow around the body. This is illustrated in Figure 1.3. Here, the shock wave from the nose of the vehicle heats the airflow around the vehicle; at the same time, the vehicle is heated by the intense frictional dissipation within the boundary layer on the surface. Allen reasoned that if more of the total reentry energy could be dumped into the airflow, then less would be available to be transferred to the vehicle itself in the form of heating. In turn, the way to increase the heating of the airflow is to create a stronger shock wave at the nose, i. e., to use a blunt-nosed body. The contrast between slender and blunt reentry bodies is illustrated in Figure 1.4. This was a stunning conclusion—to minimize aerodynamic heating, you actually want a blunt rather than a slender body. The result was so important that it was bottled up in a secret government document. Moreover, because it was so foreign to contemporary intuition, the blunt-reentry-body concept was accepted only gradually by the technical community. Over the next few years, additional aerodynamic analyses and experiments confirmed the validity of blunt reentry bodies. By 1955, Allen was publicly recognized for his work, receiving the Sylvanus Albert Reed Award of the Institute of the Aeronautical Sciences (now the American Institute of Aeronautics and Astronautics). Finally, in 1958, his work was made available to the public in the pioneering document NACA Report 1381 entitled “A Study of the Motion and Aerodynamic Heating of Ballistic Missiles Entering the

Подпись: Hot boundary layerПодпись: Energy of reentry goes into heating both the body and the air around the body.image9"Very high­speed flow

Figure 1.3

Very high­speed flow

 

Thoughts

image10

Thoughts

Earth’s Atmosphere at High Supersonic Speeds.” Since Harvey Allen’s early work, all successful reentry bodies, from the first Atlas ICBM to the manned Apollo lunar capsule, have been blunt. Incidentally, Allen went on to distinguish himself in many other areas, becoming the director of the NASA Ames Research Center in 1965, and retiring in 1970. His work on the blunt reentry body is an excellent example of the importance of aerodynamics to space vehicle design.

In summary, the purpose of this section has been to underscore the importance of aerodynamics in historical context. The goal of this book is to introduce the fundamentals of aerodynamics and to give the reader a much deeper insight to many technical applications in addition to the few described above. Aerodynamics is also a subject of intellectual beauty, composed and drawn by many great minds over the centuries. If you are challenged and interested by these thoughts, or even the least bit curious, then read on.

Physical principle Energy can be neither created nor destroyed; it can only change in form

This physical principle is embodied in the first law of thermodynamics. A brief review of thermodynamics is given in Chapter 7. Thermodynamics is essential to the study of compressible flow; however, at this stage, we will only introduce the first law, and we defer any substantial discussion of thermodynamics until Chapter 7, where we begin to concentrate on compressible flow.

Consider a fixed amount of matter contained within a closed boundary. This matter defines the system. Because the molecules and atoms within the system are constantly in motion, the system contains a certain amount of energy. For simplicity, let the system contain a unit mass; in turn, denote the internal energy per unit mass by e.

The region outside the system defines the surroundings. Let an incremental amount of heat 8q be added to the system from the surroundings. Also, let 8w be the work done on the system by the surroundings. (The quantities 8q and 8w are discussed in more detail in Chapter 7.) Both heat and work are forms of energy, and when added to the system, they change the amount of internal energy in the system. Denote this change of internal energy by de. From our physical principle that energy is conserved, we have for the system

8q + 8w = de [2.85]

Equation (2.85) is a statement of the first law of thermodynamics.

Let us apply the first law to the fluid flowing through the fixed control volume shown in Figure 2.17. Let

B = rate of heat added to fluid inside control volume from surroundings

Z?2 = rate of work done on fluid inside control volume

Вт, = rate of change of energy of fluid as it flows through control volume

From the first law,

Note that each term in Equation (2.86) involves the time rate of energy change; hence, Equation (2.86) is, strictly speaking, a power equation. However, because it is a statement of the fundamental principle of conservation of energy, the equation is conventionally termed the “energy equation.” We continue this convention here.

Подпись: Rate of volumetric
Подпись: heating = Подпись: qpdV Подпись: [2.87]

First, consider the rate of heat transferred to or from the fluid. This can be visual­ized as volumetric heating of the fluid inside the control volume due to absorption of radiation originating outside the system or the local emission of radiation by the fluid itself, if the temperature inside the control volume is high enough. In addition, there may be chemical combustion processes taking place inside the control volume, such as fuel-air combustion in a jet engine. Let this volumetric rate of heat addition per unit mass be denoted by q. Typical units for q are J/s • kg or ft • lb/s • slug. Examining Figure 2.17, the mass contained within an elemental volume is p dV hence, the rate of heat addition to this mass is q{p dV). Summing over the complete control volume, we obtain

In addition, if the flow is viscous, heat can be transferred into the control volume by means of thermal conduction and mass diffusion across the control surface. At this stage, a detailed development of these viscous heat-addition terms is not warranted; they are considered in detail in Chapter 15. Rather, let us denote the rate of heat addition to the control volume due to viscous effects simply by Qviscous – Therefore, in Equation (2.86), the total rate of heat addition is given by Equation (2.87) plus

Q viscous*

Physical principle Energy can be neither created nor destroyed; it can only change in form

[2.88]

 

в 1

 

V

 

Before considering the rate of work done on the fluid inside the control volume, consider a simpler case of a solid object in motion, with a force F being exerted on the object, as sketched in Figure 2.23. The position of the object is measured from a fixed origin by the radius vector r. In moving from position Г] to Г2 over an interval

 

Later time, t + dt

 

Figure 2.23 Schematic for the rate of doing work by a force F exerted on a moving body.

 

image137image138

of time dt, the object is displaced through dr. By definition, the work done on the object in time dt is F • dr. Hence, the time rate of doing work is simply F • йт/dt. However, йт/dt = V, the velocity of the moving object. Hence, we can state that

Rate of doing work on moving body = F • V

In words, the rate of work done on a moving body is equal to the product of its velocity and the component of force in the direction of the velocity.

Physical principle Energy can be neither created nor destroyed; it can only change in form

This result leads to an expression for B2, as follows. Consider the elemental area dS of the control surface in Figure 2.17. The pressure force on this elemental area is —p dS. From the above result, the rate of work done on the fluid passing through dS with velocity V is (—p dS) • V. Hence, summing over the complete control surface, we have

If the flow is viscous, the shear stress on the control surface will also perform work on the fluid as it passes across the surface. Once again, a detailed development of this term is deferred until Chapter 15. Let us denote this contribution simply by VTViSCOUS. Then the total rate of work done on the fluid inside the control volume is the sum of Equations (2.89) and (2.90) and VTviscou

Подпись: ШПодпись: [2.91]Подпись: p(f • Y)dV + Wvi +

S V

To visualize the energy inside the control volume, recall that in the first law of thermodynamics as stated in Equation (2.85), the internal energy e is due to the random motion of the atoms and molecules inside the system. Equation (2.85) is written for a stationary system. However, the fluid inside the control volume in Figure 2.17 is not stationary; it is moving at the local velocity V with a consequent kinetic energy per unit mass of Vі/2. Hence, the energy per unit mass of the moving fluid is the sum of both internal and kinetic energies e + Vі/2. This sum is called the total energy per unit mass.

We are now ready to obtain an expression for ZL, the rate of change of total energy of the fluid as it flows through the control volume. Keep in mind that mass flows into the control volume of Figure 2.17 bringing with it a certain total energy; at the same time mass flows out of the control volume taking with it a generally different amount of total energy. The elemental mass flow across dS is pV • dS, and therefore

Подпись: Net rate of flow of total energy across control surface Подпись: [2.92]
image139 image140

the elemental flow of total energy across dS is (pV • dS)(e + V2/2). Summing over the complete control surface, we obtain

In addition, if the flow is unsteady, there is a time rate of change of total energy inside the control volume due to the transient fluctuations of the flow-field variables. The total energy contained in the elemental volume dV is p(<? + V2/2) dV, and hence the total energy inside the complete control volume at any instant in time is

Therefore,

Подпись: dVimage141"Подпись: v and (2.93): Подпись: [2.93]Time rate of change of total energy
inside V due to transient variations
of flow-field variables

Подпись: S3 Physical principle Energy can be neither created nor destroyed; it can only change in form Physical principle Energy can be neither created nor destroyed; it can only change in form Подпись: [2.94]
image142

In turn, S3 is the sum of Equations (2.92)

image143

Repeating the physical principle stated at the beginning of this section, the rate of heat added to the fluid plus the rate of work done on the fluid is equal to the rate of change of total energy of the fluid as it flows through the control volume; i. e., energy is conserved. In turn, these words can be directly translated into an equation by combining Equations (2.86), (2.88), (2.91), and (2.94):

Equation (2.95) is the energy equation in integral form; it is essentially the first law of thermodynamics applied to a fluid flow.

For the sake of completeness, note that if a shaft penetrates the control surface in Figure 2.17, driving some power machinery located inside the control volume (say, a compressor of a jet engine), then the rate of work delivered by the shaft, lVShafb must be added to the left side of Equation (2.95). Also note that the potential energy does not appear explicitly in Equation (2.95). Changes in potential energy are contained in the body force term when the force of gravity is included in f. For the aerodynamic

problems considered in this book, shaft work is not treated, and changes in potential energy are always negligible.

Following the approach established in Sections 2.4 and 2.5, we can obtain a partial differential equation for total energy from the integral form given in Equation

(2.95) . Applying the divergence theorem to the surface integrals in Equation (2.95), collecting all terms inside the same volume integral, and setting the integrand equal to zero, we obtain

3

Г ( v2X

г / v2 1

dt

/Т + т)_

+ V.

4c+t)v

= pq-V • (p) + p(f ■ V)

+ &. + w!-

1 ^ VISCOUS 1 VISCOUS

[2.96]

where <2viscous and ^viscous represent the proper forms of the viscous terms, to be obtained in Chapter 15. Equation (2.96) is a partial differential equation which relates the flow-field variables at a given point in space.

If the flow is steady (3/31 = 0), inviscid ((/viscous = 0 and Wviscous = 0), adiabatic (no heat addition, q = 0), without body forces (f = 0), then Equations

(2.95) aEquations (2.97) and (2.98) are discussed and applied at length beginning with Chapter 7.

With the energy equation, we have introduced another unknown flow-field vari­able e. We now have three equations, continuity, momentum, and energy, which involve four dependent variables, p, p, V, and e. A fourth equation can be obtained from a thermodynamic state relation for e (see Chapter 7). If the gas is calorically perfect, then

e = cvT [2.99]

where cv is the specific heat at constant volume. Equation (2.99) introduces temper­ature as yet another dependent variable. However, the system can be completed by using the perfect gas equation of state

p — pRT [2.100]

where R is the specific gas constant. Therefore, the continuity, momentum, and energy equations, along with Equations (2.99) and (2.100) are five independent equations

for the five unknowns, p p, V, e, and T. The matter of a perfect gas and related equations of state are reviewed in detail in Chapter 7; Equations (2.99) and (2.100) are presented here only to round out our development of the fundamental equations of fluid flow.

Applied Aerodynamics: the Flow over a Circular Cylinder—the Real Case

The inviscid, incompressible flow over a circular cylinder was treated in Section 3.13. The resulting theoretical streamlines are sketched in Figure 3.26, characterized by a

symmetrical pattern where the flow “closes in” behind the cylinder. As a result, the pressure distribution over the front of the cylinder is the same as that over the rear (see Figure 3.29). This leads to the theoretical result that the pressure drag is zero—d’Alembert’s paradox.

The real flow over a circular cylinder is quite different from that studied in Section 3.13, the difference due to the influence of friction. Moreover, the drag coefficient for the real flow over a cylinder is certainly not zero. For a viscous incompressible flow, the results of dimensional analysis (Section 1.7) clearly demonstrate that the drag coefficient is a function of the Reynolds number. The variation of Co = /(Re) for a circular cylinder is shown in Figure 3.44, which is based on a wealth of experimental data. Here, Re = (PooFoo^O/P-co. where d is the diameter of the cylinder. Note that Co is very large for the extremely small values of Re < 1, but decreases monotonically until Re ^ 300, 000. At this Reynolds number, there is a precipitous drop of CD from a value near 1 to about 0.3, then a slight recovery to about 0.6 for Re = 107. (Note: These results are consistent with the comparison shown in Figure 139d and e, contrasting С о for a circular cylinder at low and high Re.) What cause this precipitous drop in Co when the Reynolds number reaches about 300,000? A detailed answer must await our discussion of viscous flow in Part 4. However, we state now that the phenomenon is caused by a sudden transition of laminar flow within the boundary layer at the lower values of Re to a turbulent boundary layer at the higher values of Re. Why does a turbulent boundary layer result in a smaller Co for this case? Stay tuned; the answer is given in Part 4.

image287

Figure 3.44 Variation of cylinder-drag coefficient with Reynolds number.

(Source: Experimental data as compiled in Panton, Ronald, Incompressible Flow, Wiley-lnterscience, New York, 1 984.)

The variation of Сд shown in Figure 3.44 across a range of Re from 10-1 to 107 is accompanied by tremendous variations in the qualitative aspects of the flow field, as itemized, and as sketched in Figure 3.45.

image288

(el

image289

image290

(0

image291

image292

Figure 3.45 Various types of flow over a circular cylinder.

[Source: Panton, Ronald, Incompressible Flow, Wiley-lnterscience, New York, 1984.)

1. For very low values of Re, say, 0 < Re < 4, the streamlines are almost (but not exactly) symmetrical, and the flow is attached, as sketched in Figure 3.45a. This regime of viscous flow is called Stokes flow, it is characterized by a near balance of pressure forces with friction forces acting on any given fluid element; the flow velocity is so low that inertia effects are very small. A photograph of this type of flow is shown in Figure 3.46, which shows the flow of water around a circular cylinder where Re = 1.54. The streamlines are made visible by aluminum powder on the surface, along with a time exposure of the film.

2. For 4 < Re < 40, the flow becomes separated on the back of the cylinder, forming two distinct, stable vortices that remain in the position shown in Figure 3.45b. A photograph of this type of flow is given in Figure 3.47, where Re = 26.

3. As Re is increased above 40, the flow behind the cylinder becomes unstable; the vortices which were in a fixed position in Figure 3.45b now are alternately shed from the body in a regular fashion and flow downstream. This flow is sketched in Figure 3.45c. A photograph of this type of flow is shown in Figure 3.48, where Re = 140. This is a water flow where the streaklines are made visible by the electrolytic precipitation method. (In this method, metal plating on the cylinder

image293

Figure 3.46 Flow over a circular cylinder. Re = 1.54. (Source: Photograph by

Sadatoshi Taneda, from Van Dyke, Milton, An Album of Fluid Motion, The Parabolic Press, Stanford, Calif., 1982.)

Подпись: Figure 3.47 Flow over a circular cylinder. Re = 26. (Source: Photograph by Sadatoshi Taneda, from Van Dyke, Milton, An Album of Fluid Motion, The Parabolic Press, Stanford, Calif., 1982.)
surface acts as an anode, white particles are precipitated by electrolysis near the anode, and these particles subsequently flow downstream, forming a streakline. The definition of a streakline is given in Section 2.11) The alternately shed vortex pattern shown in Figures 3.45c and 3.48 is called a Karman vortex street, named after Theodore von Karman, who began to study and analyze this pattern in 1911 while at Gottingen University in Germany, (von Karman subsequently had a long and very distinguished career in aerodynamics, moving to the California Institute of Technology in 1930, and becoming America’s best-known aerodynamicist in the mid-twentieth century. An autobiography of von Karman was published in 1967; see Reference 49. This reference is “must” reading for anyone interested in a riveting perspective on the history of aerodynamics in the twentieth century.)

4. As the Reynolds number is increased to large numbers, the Karman vortex street becomes turbulent and begins to metamorphose into a distinct wake. The laminar boundary layer on the cylinder separates from the surface on the forward face, at a point about 80° from the stagnation point. This is sketched in Figure 3.45J. The value of the Reynolds number for this flow is on the order of 105. Note, from Figure 3.44, that Co is a relatively constant value near unity for 103 < Re < 3 x 105.

5. For 3 x 105 < Re < З X 106, the separation of the laminar boundary layer still takes place on’the forward face of the cylinder. However, in the free shear layer over the top of the separated region, transition to turbulent flow takes place. The

Подпись: Figure 3.48 Flow over a circular cylinder. Re = 140. A Karman vortex street exists behind the cylinder at this Reynolds number. (Source: Photograph by Sadatoshi Taneda, from Van Dyke, Milton, An Album of Fluid Motion, The Parabolic Press, Stanford, Calif., 1982.)

flow then reattaches on the back face of the cylinder, but separates again at about 120° around the body measured from the stagnation point. This flow is sketched in Figure 3.45e. This transition to turbulent flow, and the corresponding thinner wake (comparing Figure 3.45e with Figure 3.45d), reduces the pressure drag on the cylinder and is responsible for the precipitous drop in Co at Re = 3 x 105 shown in Figure 3.44. (More details on this phenomenon are covered in Part 4.)

6. For Re < 3 x 106, the boundary layer transits directly to turbulent flow at some point on the forward face, and the boundary layer remains totally attached over the surface until it separates at an angular location slightly less than 120° on the back surface. For this regime of flow, CD actually increases slightly with increasing Re because the separation points on the back surface begin to move closer to the top and bottom of the cylinder, producing a fatter wake, and hence larger pressure drag.

In summary, from the photographs and sketches in this section, we see that the real flow over a circular cylinder is dominated by friction effects, namely, the separation of the flow over the rearward face of the cylinder. In turn, a finite pressure drag is created on the cylinder, and d’Alembert’s paradox is resolved.

Let us examine the production of drag more closely. The theoretical pressure distribution over the surface of a cylinder in an inviscid, incompressible flow was given

image296

Figure 3.49 Pressure distribution over a circular cylinder in low-speed flow. Comparison of the theoretical pressure distribution with two experimental pressure distributions—one for a subcritical Re and the other for a supercritical Re.

in Figure 3.29. In contrast, several real pressure distributions based on experimental measurements for different Reynolds numbers are shown in Figure 3.49, and are compared with the theoretical inviscid flow results obtained in Section 3.13. Note that theory and experiment agree well on the forward face of the cylinder, but that dramatic differences occur over the rearward face. The theoretical results show the pressure decreasing around the forward face from the initial total pressure at the stagnation point, reaching a minimum pressure at the top and bottom of the cylinder (i9 — 90° and 270°), and then increasing again over the rearward face, recovering to the total pressure at the rear stagnation point. In contrast, in the real case where flow separation occurs, the pressures are relatively constant in the separated region over the rearward face and have values slightly less than freestream pressure. (In regions of separated flow, the pressure frequently exhibits a nearly constant value.) In the separated region over the rearward face, the pressure clearly does not recover to the higher values that exist on the front face. There is a net imbalance of the pressure distribution between the front and back faces, with the pressures on the front being higher than on the back, and this imbalance produces the drag on the cylinder.

Return to Figure 3.44, and examine again the variation of CD as a function of Re. The regimes associated with the very low Reynolds numbers, such as Stokes flow for Re «5 1, are usually of no interest to aeronautical applications. For example, consider a circular cylinder in an airflow of 30 m/s (about 100 ft/s, or 68 mi/h) at standard sea level conditions, where Poo = 1-23 kg/m3 and = 1.79 x 10“5 kg/(m • s). The smaller the diameter of the cylinder, the smaller will be the Reynolds number.

speed of the SPAD was 130 mi/h, or 57.8 m/s. For this velocity at standard sea level, we have

With this value of Re, we are beginning to enter the world of practical aerodynamics

Applied Aerodynamics: the Flow over a Circular Cylinder—the Real Case

image297

Figure 3.50 The French SPAD XIII, an example of a strut-and-wire biplane from World War I. Captain Eddie Rickenbacker is shown at the front of the airplane. (Courtesy of U. S. Air Force.)

 

for the flow over cylinders. It is interesting to note that, from Figure 3.44, Сд = 1 for the wires on the SPAD. In terms of airplane aerodynamics, this is a high drag coefficient for any component of an aircraft. Indeed, the bracing wires used on biplanes of the World War I era were a source of high drag for the aircraft, so much so that early in the war, bracing wire with a symmetric airfoil-like cross section was utilized to help reduce this drag. Such wire was developed by the British at the Royal Aircraft Factory at Farnborough, and was first tested experimentally as early as 1914 on an SE-4 biplane. Interestingly enough, the SPAD used ordinary round wire, and in spite of this was the fastest of all World War I aircraft.

This author was struck by another example of the effect of cylinder drag while traveling in Charleston, South Carolina, shortly after hurricane Hugo devastated the area on September 28, 1989. Traveling north out of Charleston on U. S. Route 17, near the small fishing town of McClellanville, one passes through the Francis Marion National Forest. This forest was virtually destroyed by the hurricane; 60-ft pine trees were snapped off near their base, and approximately 8 out of every 10 trees were down. The sight bore an eerie resemblance to scenes from the battlefields in France during World War I. What type of force can destroy an entire forest in this fashion? To answer this question, we note that the Weather Bureau measured wind gusts as high as 175 mi/h during the hurricane. Let us approximate the wind force on a typical 60-ft pine tree by the aerodynamic drag on a cylinder of a length of 60 ft and a diameter of 5 ft. Since V = 175 mi/h = 256.7 ft/s, Poo = 0.002377 slug/ft3, and

(0.002377) (256.7)(5)
3.7373 x 10-[11]

= 8.16 x 10[12]

Подпись: Poo vxd До©

= 3.7373 x 10-7 slug/(ft • s), then the Reynolds number is

Examining Figure 3.44, we see that Сд = 0.7. Since Сд is based on the drag per unit length of the cylinder as well as the projected frontal area, we have for the total drag exerted on an entire tree that is 60 ft tall

D = qocSCo = PooVl(d)(60)CD

= і (0.002377)(256.7)2(5)(60)(0.7) = 16,446 lb

a 16,000 lb force on the tree—it is no wonder a whole forest was destroyed. (In the above analysis, we neglected the end effects of the flow over the end of the vertical cylinder. Moreover, we did not correct the standard sea level density for the local reduction in barometric pressure experienced inside a hurricane. However, these are relatively small effects in comparison to the overall force on the cylinder.) The aerodynamics of a tree, and especially that of a forest, are more sophisticated than discussed here. Indeed, the aerodynamics of trees have been studied experimentally with trees actually mounted in a wind tunnel.6

Three-Dimensional Doublet

Подпись: or image471 Подпись: [6.8]

Consider a sink and source of equal but opposite strength located at points О and A, as sketched in Figure 6.2. The distance between the source and sink is l. Consider an arbitrary point P located a distance r from the sink and a distance rx from the source. From Equation (6.7), the velocity potential at P is

Let the source approach the sink as their strengths become infinite; that is, let

In the limit, as / —> 0, r — rx —> OB = l cosd, and rrx —»■ r2. Thus, in the limit, Equation (6.8) becomes

Подпись: X lcos 9 An r2, X r-rx

Ф = — hm————–

An rrx

Three-Dimensional Doublet

Three-Dimensional Doublet

or

 

[6.9]

 

image472

Figure 6.2 Source-sink pair. In the limit as / -»• 0, a three-dimensional doublet is obtained.

 

where д = Л/. The flow field produced by Equation (6.9) is a three-dimensional dou­blet’ д is defined as the strength of the doublet. Compare Equation (6.9) with its two­dimensional counterpart given in Equation (3.88). Note that the three-dimensional effects lead to an inverse r-squared variation and introduce a factor 47Г, versus 2тг for the two-dimensional case.

Three-Dimensional Doublet Подпись: [6.10]

From Equations (2.18) and (6.9), we find

The streamlines of this velocity field are sketched in Figure 6.3. Shown are the streamlines in the zr plane; they are the same in all the zr planes (i. e., for all values of Ф). Hence, the flow induced by the three-dimensional doublet is a series of stream surfaces generated by revolving the streamlines in Figure 6.3 about the г axis. Compare these streamlines with the two-dimensional case illustrated in Figure 3.18; they are qualitatively similar but quantitatively different.

Note that the flow in Figure 6.3 is independent of Ф; indeed, Equation (6.10) clearly shows that the velocity field depends only on r and 9. Such a flow is defined as axisymmetric flow. Once again, we have a flow with two independent variables. For this reason, axisymmetric flow is sometimes labeled “two-dimensional” flow. However, it is quite different from the two-dimensional planar flows discussed ear­lier. In reality, axisymmetric flow is a degenerate three-dimensional flow, and it is somewhat misleading to refer to it as “two-dimensional.” Mathematically, it has only two independent variables, but it exhibits some of the same physical characteristics as general three-dimensional flows, such as the three-dimensional relieving effect to be discussed later.

Three-Dimensional Doublet

image473

Figure 6.3 Sketch of the streamlines in the zr plane

(Ф = constant plane) for a three-dimensional doublet.

 

Mach Number Regimes

Of all the ways of subdividing and describing different aerodynamic flows, the dis­tinction based on the Mach number is probably the most prevalent. If M is the local Mach number at an arbitrary point in a flow field, then by definition the flow is locally:

Subsonic if M < I Sonic if M = 1 Supersonic if M > 1

Looking at the whole field simultaneously, four different speed regimes can be iden­tified using Mach number as the criterion:

1. Subsonic flow (M < 1 everywhere). A flow field is defined as subsonic if the Mach number is less than 1 at every point. Subsonic flows are characterized by smooth streamlines (no discontinuity in slope), as sketched in Figure 1.37a. Moreover, since the flow velocity is everywhere less than the speed of sound, disturbances in the flow (say, the sudden deflection of the trailing edge of the airfoil in Figure 1.37a) propagate both upstream and downstream, and are felt throughout the entire flow field. Note that a freestream Mach number M^ less than 1 does not guarantee a totally subsonic flow over the body. In expanding over an aerodynamic shape, the flow velocity increases above the freestream value, and if Moo is close enough to 1, the local Mach number may become supersonic in certain regions of the flow. This gives rise to a rule of thumb that Мж <0.8 for subsonic flow over slender bodies. For blunt bodies, Мж must be even lower to ensure totally subsonic flow. (Again, emphasis is made that the above is just a loose rule of thumb and should not be taken as a precise quantitative definition.) Also, we will show later that incompressible flow is a special limiting case of subsonic flow where M -> 0.

2. Transonic flow (mixed regions where M < landM > 1). As stated above, if

is subsonic but is near unity, the flow can become locally supersonic (M > 1). This is sketched in Figure 1.37b, which shows pockets of supersonic flow over both the top and bottom surfaces of the airfoil, terminated by weak shock waves behind which the flow becomes subsonic again. Moreover, if Mx is increased slightly above unity, a bow shock wave is formed in front of the body; behind this shock wave the flow is locally subsonic, as shown in Figure 1.37c. This subsonic flow subsequently expands to a low supersonic value over the airfoil. Weak shock waves are usually generated at the trailing edge, sometimes in a “fishtail” pattern as shown in Figure 1.37c. The flow fields shown in Figure 1.37b and c are characterized by mixed subsonic-supersonic flows and are dominated by the physics of both types of flow. Hence, such flow fields are called transonic flows. Again, as a rule of thumb for slender bodies, transonic flows occur for freestream Mach numbers in the range 0.8 < Мж < 1.2.

3. Supersonic flow (M > 1 everywhere). A flow field is defined as supersonic if the Mach number is greater than 1 at every point. Supersonic flows are frequently characterized by the presence of shock waves across which the flow properties

(a) Subsonic flow

 

Mach Number Regimes

Thin, hot shock layer with viscous interaction and chemical reactions

 

Moo> 5

 

(e) Hypersonic flow

 

image65image66image68

Mach Number Regimes

Figure 1.37 Different regimes of flow.

and streamlines change discontinuously (in contrast to the smooth, continuous variations in subsonic flows). This is illustrated in Figure 1.37d for supersonic flow over a sharp-nosed wedge; the flow remains supersonic behind the oblique shock wave from the tip. Also shown are distinct expansion waves, which are common in supersonic flow. (Again, the listing of Мж > 1.2 is strictly a rule of thumb. For example, in Figure 131 d, if 9 is made large enough, the oblique shock

wave will detach from the tip of the wedge, and will form a strong, curved bow shock ahead of the wedge with a substantial region of subsonic flow behind the wave. Hence, the totally supersonic flow sketched in Figure .31d is destroyed if 9 is too large for a given M^. This shock detachment phenomenon can occur at any value of M^ > 1, but the value of в at which it occurs increases as Moo increases. In turn, if 9 is made infinitesimally small, the flow field in Figure l.31d holds for Moo > 1.0. These matters will be considered in detail in Chapter 9. However, the above discussion clearly shows that the listing of Moo > 1-2 in Figure 1.31 d is a very tenuous rule of thumb, and should not be taken literally.) In a supersonic flow, because the local flow velocity is greater than the speed of sound, disturbances created at some point in the flow cannot work their way upstream (in contrast to subsonic flow). This property is one of the most significant physical differences between subsonic and supersonic flows. Is is the basic reason why shock waves occur in supersonic flows, but do not occur in steady subsonic flow. We will come to appreciate this difference more fully in Chapters 7 to 14.

4. Hypersonic flow (very high supersonic speeds). Refer again to the wedge in Fig­ure .31d. Assume 9 is a given, fixed value. As Мж increases above 1, the shock wave moves closer to the body surface. Also, the strength of the shock wave in­creases, leading to higher temperatures in the region between the shock and the body (the shock layer). If Мх is sufficiently large, the shock layer becomes very thin, and interactions between the shock wave and the viscous boundary layer on the surface occur. Also, the shock layer temperature becomes high enough that chemical reactions occur in the air. The O2 and N2 molecules are torn apart; i. e., the gas molecules dissociate. When Мх becomes large enough such that viscous interaction and/or chemically reacting effects begin to dominate the flow (Figure 1.37e), the flow field is called hypersonic. (Again, a somewhat arbitrary but frequently used rule of thumb for hypersonic flow is Мж > 5.) Hypersonic aerodynamics received a great deal of attention during the period 1955-1970 because atmospheric entry vehicles encounter the atmosphere at Mach numbers between 25 (ICBMs) and 36 (the Apollo lunar return vehicle). Again during the period 1985-1995, hypersonic flight received a great deal of attention with the concept of air-breathing supersonic – combustion ramjet-powered transatmo­spheric vehicles to provide single-stage-to-orbit capability. Today, hypersonic aerodynamics is just part of the whole spectrum of realistic flight speeds. Some basic elements of hypersonic flow are treated in Chapter 14.

In summary, we attempt to organize our study of aerodynamic flows according to one or more of the various categories discussed in this section. The block diagram in Figure 1.38 is presented to help emphasize these categories and to show how they are related. Indeed, Figure 1.38 serves as a road map for this entire book. All the material to be covered in subsequent chapters fits into these blocks, which are lettered for easy reference. For example, Chapter 2 contains discussions of some fundamental aerodynamic principles and equations which fit into both blocks C and D. Chapters 3 to 6 fit into blocks D and E, Chapter 7 fits into blocks D and F, etc. As we proceed

image69

Figure 1.38 Block diagram categorizing the types of aerodynamic flows.

with our development of aerodynamics, we will frequently refer to Figure 1.38 in order to help put specific, detailed material in proper perspective relative to the whole of aerodynamics.

Theoretical (Analytical) Solutions

Students learning any field of physical science or engineering at the beginning are usually introduced to nice, neat analytical solutions to physical problems that are simplified to the extent that such solutions are possible. For example, when Newton’s second law is applied to the motion of a simple, frictionless pendulum, students in elementary physics classes are shown a closed form analytical solution for the time period of the pendulum’s oscillation, namely,

T = 2njq~g

where T is the period, і is the length of the pendulum, and g is the acceleration of gravity. However, a vital assumption behind this equation is that of small amplitude oscillations. Similarly, in studying the motion of a freely falling body in a gravitational field, the distance у through which the body falls in time t after release is given by

image182

However, this result neglects any aerodynamic drag on the body as it falls through the air. The above examples are given because they are familiar results from elementary physics. They are examples of theoretical, closed-form solutions—straightforward algebraic relations.

The governing equations of aerodynamics, such as the continuity, momentum, and energy equations derived in Sections 2.4, 2.5, and 2.7, respectively, are highly non-linear, partial differential, or integral equations; to date, no general analytical solution to these equations has been obtained. In lieu of this, two different philoso­phies have been followed in obtaining useful solutions to these equations. One of these is the theoretical, or analytical, approach, wherein the physical nature of certain aerodynamic applications allows the governing equations to be simplified to a suffi­cient extent that analytical solutions of the simplified equations can be obtained. One such example is the analysis of the flow across a normal shock wave, as discussed in Chapter 8. This flow is one-dimensional, i. e., the changes in flow properties across the shock take place only in the flow direction. For this case, the у and z derivatives in the governing continuity, momentum, and energy equations from Sections 2.4, 2.5, and 2.7 are identically zero, and the resulting one-dimensional equations, which are still exact for the one-dimensional case being considered, lend themselves to a direct analytical solution, which is indeed an exact solution for the shock wave properties. Another example is the compressible flow over an airfoil considered in Chapters 11 and 12. If the airfoil is thin and at a small angle of attack, and if the freestream Mach number is not near one (not transonic) nor above five (not hypersonic), then many of the terms in the governing equations are small compared to others and can be neglected. The resulting simplified equations are linear and can be solved ana­lytically. This is an example of an approximate solution, where certain simplifying assumptions have been made in order to obtain a solution.

The history of the development of aerodynamic theory is in this category—the simplification of the full governing equations apropos a given application so that analytical solutions can be obtained. Of course this philosophy works for only a limited number of aerodynamic problems. However, classical aerodynamic theory is built on this approach and, as such, is discussed at some length in this book. You can expect to see a lot of closed-form analytical solutions in the subsequent chapters, along with detailed discussions of their limitations due to the approximations necessary to obtain such solutions. In the modem world of aerodynamics, such classical analytical solutions have three advantages:

1. The act of developing these solutions puts you in intimate contact with all the physics involved in the problem.

2. The results, usually in closed-form, give you direct information on what are the important variables, and how the answers vary with increases or decreases in these variables. For example, in Chapter 11 we will obtain a simple equation for the compressibility effects on lift coefficient for an airfoil in high-speed subsonic flow. The equation, Equation (11.52), tells us that the high-speed effect on lift coefficient is governed by just Мж alone, and that as increases, then the lift coefficient increases. Moreover, the equation tells us in what way the lift coefficient increases, namely, inversely with (1 — M^)1/2. This is powerful information, albeit approximate.

3. Finally, the results in closed-form provide simple tools for rapid calculations, making possible the proverbial “back of the envelope calculations” so important in the preliminary design process and in other practical applications.

The Cambered Airfoil

image332

Thin airfoil theory for a cambered airfoil is a generalization of the method for a symmetric airfoil discussed in Section 4.7. To treat the cambered airfoil, return to Equation (4.18):

image333,image334,image335,image336,image337
Recall that Equation (4.46) was obtained directly from Equation (4.42), which is the transformed version of the fundamental equation of thin airfoil theory, Equa­tion (4.18). Furthermore, recall that Equation (4.18) is evaluated at a given point, r along the chord line, as sketched in Figure 4.19. Hence, Equation (4.46) is also eval­uated at the given point x; here, dz/dx and во correspond to the same point л: on the chord line. Also, recall that dz/dx is a function of во, where л: = (c/2)(l — cos0o) from Equation (4.21).

Examine Equation (4.46) closely. It is in the form of a Fourier cosine series ex­pansion for the function of dz/dx. In general, the Fourier cosine series representation of a function f (в) over an interval 0 < в < ж is given by

OO

/((9) = B0 + Bn cos пв [4.47]

n=I

where, from Fourier analysis, the coefficients Bo and Bn are given by

Bo = – f f (в) dO [4.48]

л Jo

2 Г

and Bn — — / f {6)cosn6 d6 [4.49]

л – Jo

(See, e. g., page 217 of Reference 6.) In Equation (4.46), the function dz/dx is analo­gous to f(6) in the general form given in Equation (4.47). Thus, from Equations (4.48) and (4.49), the coefficients in Equation (4.46) are given by

1 Г dz

a – A0 = — —dOо Jr Jo dx

or

1 [л dz

Ao = a—— / — dOo

Jr Jo dx

[4.50]

and

2 Cn dz

An = — / — cos пво dOо n Jo dx

[4.51]

Keep in mind that in the above, dz/dx is a function of в0. Note from Equation (4.50) that Ao depends on both a and the shape of the camber line (through dz/dx), whereas from Equation (4.51) the values of An depend only on the shape of the camber line.

Pause for a moment and think about what we have done. We are considering the flow over a cambered airfoil of given shape dz/dx at a given angle of attack a. In order to make the camber line a streamline of the flow, the strength of the vortex sheet along the chord line must have the distribution у (в) given by Equation (4.43), where the coefficients A0 and An are given by Equations (4.50) and (4.51), respectively. Also, note that Equation (4.43) satisfies the Kutta condition y(n) = 0. Actual numbers for Ao and An can be obtained for a given shape airfoil at a given angle of attack simply by carrying out the integrations indicated in Equations (4.50) and (4.51). For an example of such calculations applied to an NACA 2412 airfoil, see pages 120­125 of Reference 13. Also, note that when dz/dx = 0, Equation (4.43) reduces to Equation (4.24) for a symmetric airfoil. Hence, the symmetric airfoil is a special case of Equation (4.43).

Let us now obtain expressions for the aerodynamic coefficients for a cambered airfoil. The total circulation due to the entire vortex sheet from the leading edge to the trailing edge is

Г= f y(^)d^=C – [ y(0)sm0de JO ^ Jo


The Cambered Airfoil

[4.53]

 

The Cambered Airfoil

and

 

The Cambered Airfoil

Hence, Equation (4.53) becomes

Подпись: [4.54]Г — cVqo (ttAq +

From Equation (4.54), the lift per unit span is

L — Poo 1^00 Г — Рсо^оо^ (жА0 3“

In turn, Equation (4.55) leads to the lift coefficient in the form

L’

 

[4.55]

 

[4.56]

 

The Cambered Airfoil

ci

 

The Cambered Airfoil

Recall that the coefficients A0 and A i in Equation (4.56) are given by Equations (4.50) and (4.51), respectively. Hence, Equation (4.56) becomes

Подпись: dci da [4.57]

and

 

Lift slope =

 

[4.58]

 

— 2tz

 

Equations (4.57) and (4.58) are important results. Note that, as in the case of the symmetric airfoil, the theoretical lift slope for a cambered airfoil is 2jt. It is a general result from thin airfoil theory that dci/da = 2jt for any shape airfoil. However, the expression for c; itself differs between a symmetric and a cambered airfoil, the difference being the integral term in Equation (4.57). This integral term has physical significance, as follows. Return to Figure 4.4, which illustrates the lift curve for an airfoil. The angle of zero lift is denoted by <xl=q and is a negative value. From the geometry shown in Figure 4.4, clearly

dci, ,

ci = — (a – aL=Q) da

Substituting Equation (4.58) into (4.59), we have

ci = 2n{a — aL=o) [4.60]

Comparing Equations (4.60) and (4.57), we see that the integral term in Equa­tion (4.57) is simply the negative of the zero-lift angle; that is

Подпись: 1 UL= о = n Подпись:image340[4.61]

Hence, from Equation (4.61), thin airfoil theory provides a means to predict the angle of zero lift. Note that Equation (4.61) yields aL=o = 0 for a symmetric airfoil, which is consistent with the results shown in Figure 4.20. Also, note that the more highly cambered the airfoil, the larger will be the absolute magnitude of c/l=o-

Returning to Figure 4.21, the moment about the leading edge can be obtained by substituting у (в) from Equation (4.43) into the transformed version of Equa­tion (4.35). The details are left for Problem 4.9. The result for the moment coefficient is

Аг

~2

 

The Cambered Airfoil

[4.62]

 

Substituting Equation (4.56) into (4.62), we have

Г сі n і

Cm, le = — 1^— + ~(Al ~ A2)J

 

[4.63]

 

Note that, for dz/dx = 0, A = Аг = 0 and Equation (4.63) reduces to Equa­tion (4.39) for a symmetric airfoil.

The moment coefficient about the quarter chord can be obtained by substituting Equation (4.63) into (4.40), yielding

[4.64]

 

Cm, с/4 — ^ (Аг A])

 

Unlike the symmetric airfoil, where ст с/4 = 0, Equation (4.64) demonstrates that cm, c/4 is finite for a cambered airfoil. Therefore, the quarter chord is not the center of pressure for a cambered airfoil. However, note that A1 and Аг depend only on the shape of the camber line and do not involve the angle of attack. Hence, from Equation (4.64), сш, с/4 is independent of a. Thus, the quarter-chord point is the theoretical location of the aerodynamic center for a cambered airfoil.

The location of the center of pressure can be obtained from Equation (1.21):

Cm, leC

Подпись: Clu

[4.66]

Подпись:image341Equation (4.66) demonstrates that the center of pressure for a cambered airfoil varies with the lift coefficient. Hence, as the angle of attack changes, the center of pressure also changes. Indeed, as the lift approaches zero, xcp moves toward infinity; that is, it leaves the airfoil. For this reason, the center of pressure is not always a convenient point at which to draw the force system on an airfoil. Rather, the force-and-moment system on an airfoil is more conveniently considered at the aerodynamic center. (Re­turn to Figure 1.19 and the discussion at the end of Section 1.6 for the referencing of the force-and-moment

image342

Section angle of attack a, degree

Figure 4.33 Lift – and moment-coefficient data for an NACA 23012 airfoil, for comparison with the theoretical results obtained in Example 4.2.

Note that the results from thin airfoil theory for a cambered airfoil agree very well with the experimental data. Recall that excellent agreement between thin airfoil theory for a symmetric airfoil and experimental data has already been shown in Figure 4.20. Hence, all of the work we have done in this section to develop thin airfoil theory is certainly worth the effort. Moreover, this illustrates that the development of thin airfoil theory in the early 1900s was a crowning achievement in theoretical aerodynamics and validates the mathematical approach of replacing the chord line of the airfoil with a vortex sheet, with the flow tangency condition evaluated along the mean camber line.

This brings to an end our introduction to classical thin airfoil theory. Returning to our road map in Figure 4.2, we have now completed the right-hand branch.