Category Aerodynamics of Wings and Bodies

Integral Conservation Theorems for Inviscid Fluid

For later use in connection with the calculation of forces and moments on wings and bodies, we wish to be able to express these quantities in terms of the fluxes of linear and angular momentum through arbitrary control surfaces S. This approach will often be found to have a special convenience, because singularities which occur in velocities and pressures at the surface of a vehicle may not persist at great distances in the flow field around it, so that the integrations which must be carried out are facilitated. The basic tools for carrying out this task are irftegrated forms of Newton’s law of motion known as momentum theorems.

In connection with the presentation of the momentum theorems, we take the opportunity to discuss the question of conservation of other flow properties, as expressed in integral form.

System boundary at a slightly later time




Fixed control surface S and system boundar; at time t


Fig. 1-5. Control surface surroui___ 0 _

fluid volume and immersed body.


Immersed body surface a



Подпись: written
Consider any quantity E which is characteristic of the fluid particles contained within a fixed control volume V. Let V be bounded on the inside by one or more impermeable bodies, whose collective surfaces are denoted by the symbol a, and bounded on the outside by a larger fixed surface S. See Fig. 1-5. The closed system under examination is that fixed mass of fluid that happens to be contained within V at a certain instant of time. To find the rate of change of the total quantity for this system at the instant it coincides with V, we observe that this change is made up of the sum of all the local changes at points within V plus changes which occur as a result of the motion of the system boundary. Supposing that E is referred to unit volume of the fluid, the former rate of change can be




where the interchange of the operations of differentiation and integration is permissible in view of the constancy of the volume V. There is also an increment to the total amount of the quantity E as a result of the fact that the fluid is moving across the bounding surfaces S and a with a normal velocity component (Q • n). At points where this scalar product is positive, the fluid adjacent to an area element dS of the boundary takes up new positions outside this boundary, the volume per unit time passing outside the boundary being given by (Q • n) dS. Thus the rate of change of E for the system due to passage of fluid across the boundary is given by

<$E(Q ■ n) dS.


Combining these last two results, we find for the total rate of change of this generalized property for the system,

= fffd-~dV + <£f>E(Q-n)dS. (1-44)

v s+o-

This general result is now specialized to several cases of interest.

1. Conservation of Mass or Continuity. To derive an integral continuity equation we replace E by the mass per unit volume p and observe that, in the absence of sources and sinks, the total mass of the system must remain constant. Thus we are led to

/// Ы dV + ^P(Q ‘ n) dS = °- (b45)

V S+cr

For steady flow around an impermeable body of fixed position, of course, the first integral in (1-45) vanishes, and the contribution to the second integral from the inner surface must be zero because the quantity in parentheses vanishes.

We note incidentally how from Gauss’ theorem, (1-14),

jfip(Q-n)dS = IJjv-(pQ)dV. (1-46)

S+tr V

Подпись: dp dt Подпись: V • (pQ) Подпись: dV Подпись: 0. Подпись: (1-47)

Substituting into (1-45), we obtain

Since the volume V is arbitrary, the only way that this integral can vanish is for its integrand to be everywhere zero. Thus, the differential form of the continuity equation, (1-1), is confirmed.

2. Linear Momentum. Let F,- represent the vector sum of all forces applied by the surroundings to the system. According to Newton’s second law, this sum is equal to the time rate of change of linear momentum of the system, which corresponds to replacing E with the quantity pQ in (1—44). Consequently, we obtain the following generalized version of the law of conservation of momentum:

+ <ffpQ(Q ■ n) (IS. (1-48)

* V S+a

We now examine the various contributions that might appear to the force system in (1-48). If there is a conservative body force field, the left- hand side will include a quantity

jfjpFdV = JJJpVQ dV. (1-49)

v v

This will be omitted from what follows because of its relative unimportance in aeronautical applications.

The remaining external force will then be broken into two parts: the reaction (—Fbody) to the force exerted by the fluid on the body, and the force exerted across the outer boundary S by the surroundings. Recalling that n is the outward-directed normal, this latter might take the form

ffi—pn + r] dS,


where p is the pressure across S, and r is the sum of shear stress and devi – atoric normal stress exerted by the surroundings, if these are significant. We may write r as the dot product of a dyadic or tensor of deviatoric stress by the unit normal n. Since we are dealing generally with a non­viscous fluid, however, the question will not be elaborated here. Leaving out effects of shear stress, (1-48) can be modified to read

Again we remark that if the body is fixed in our coordinate system, the contribution to the second integral on the right from a will vanish. Also in steady flow the last term on the right is zero, leaving

Fbody = — <jij>lpn + pQ(Q ■ n)] dS. (1-51)


Equation (1-51) is actually the most useful form for practical applications. The specialized versions of (1-51) which occur when the flow involves
small perturbations will be discussed in a later chapter. It usually proves convenient to use an integrated form of the equations of motion (Ber­noulli’s equation) to replace the pressure in terms of the velocity field.

3. Angular Momentum. Let r be a vector of position measured from the origin about which moments are to be taken. Then it is an easy matter to derive the following counterpart of the first form of the linear momentum theorem:

Er*x F* = /// Jt(pr XQ)# + ^p(r X QKQ • n) dS- (1-52)

* V S+<7

The summation of moments on the left here can once more be broken up into a body-force term, a reaction to the moment exerted by the fluid through a and a pressure or shear moment exerted on the system over the outer boundary. Substitutions of this sort, neglecting the deviatoric stress, lead to the working form of the theorem of angular momentum:

Afbody = — ^ p(r X n) dS


S+a V

The steady-flow simplification involves dropping the integral over V and over the inner boundary a.

Подпись: E Подпись: h + у + gz Подпись: (1-54)

4. Thermodynamic Energy. Integral forms of the laws of thermo­dynamics will be found developed in detail in Chapter 2 of Shapiro (1953). Since these will have little direct usefulness in later applications and since many new definitions are involved, none of these results are reproduced here. Shapiro’s equation (2.20), for instance, provides an excellent working form of the first law. It is of interest that, when the pressure work exerted on the boundaries is included, the quantity E in the second or boundary term on the right of (1-44) is found to be

Here z is the distance vertically upward in a parallel gravity field, and h is the enthalpy per unit mass, which proves to be the effective thermo­dynamic energy in steady flow.

As a sidelight on the question of energy conservation we note that, in a constant-density fluid without body forces, the only way that energy can be stored is in kinetic form. Hence, a very convenient procedure for calculating drag, or fluid resistance, is to find the rate of addition of kinetic
energy to the fluid per unit time and to equate this to the work done by the drag. This represents a balance of mechanical rather than thermo­dynamic energy. When the fluid is compressible and there are still no dissipative mechanisms present, energy can be radiated away by compres­sion waves in an acoustic fashion. Therefore, the problem of computing drag from energy balance becomes a good deal more complicated.

Effects of Viscosity

3- 1 Introduction

Viscous flows at high Reynolds numbers constitute the most obvious example of singular perturbation problems. The viscosity multiplies the highest-order derivative terms in the Navier-Stokes equations, and these will therefore in the limit of zero viscosity degenerate to a lower order. The boundary condition of zero tangential velocity on the body (or of continuous velocity in the stream) is therefore lost. This necessitates the introduction of a thin boundary layer next to the body constituting the inner region where the inviscid equations are not uniformly valid. Unfor­tunately, only a very restricted class of viscous-flow problems can be analyzed by the direct use of the method of matched asymptotic expansions. First, only for a very limited class of bodies will the boundary layer remain attached to the body surface. When separation occurs, the location of the region where viscosity is important is no longer known a priori. The second difficulty is that for very high Reynolds numbers the flow in the boundary layer becomes unstable and transition to turbulence occurs. As yet, no complete theory for predicting turbulent flow exists.

Before considering some of the model problems that may be analyzed, we will give a short description of the qualitative effects of viscosity.

Three-Dimensional Wings in Steady, Subsonic Flow

7- 1 Compressibility Corrections for Wings

This chapter deals with the application to finite, almost-plane wings of the linearized, small-perturbation techniques introduced in Chapter 5. By way of introduction, we first review the similarity relations which govern variations in the parameter M, the flight Mach number.

In the light of the asymptotic expansion procedure, the principal un­known, from which all other needed information can be calculated, is the first-order term* Ф? in the outer expansion for the velocity potential. The term Ф)’ is connected to the more familiar perturbation velocity potential <p(x, y, z) by (5-28). The latter is governed by the differential equations and boundary conditions (5-29)-(5-30), which we reproduce here (see also Fig. 5-1):

(1 M )<Pxx “b <Pyy “f" <Pzz ——– 0,

for (x, y) on S. (5-30)

The pressure coefficient at any point in the field, including the upper and lower wing surfaces z = 0±, is found from

Cp — 2 <px.

Extending a procedure devised by Prandtl and Glauert for two-dimen­sional airfoils (see Fig. 7-1), Gothert (1940) introduced a transformation of independent and dependent variables which is equivalent to

* The zeroth-order term is, of course, the free stream Фо = x.


Where /3 = /l — M2, as in Chapter 5. Equation (7-1) converts (5-29) into the constant-density perturbation equation

(<Po)x0x0 + (‘Po)y0Vo + (<Po)z0zo = 0. (7-2)

Some care must be observed when interpreting the transformed boundary condition at the wing surface. Thus, for example, the first of (5-30) states that just above the wing’s projection on the x, у-plane the vertical velocity component produced by the sheet of singularities representing the wing’s disturbance must have certain values, say Fu(x, y). After transformation, we obtain

(<Po)z0 = Fu(0x о, у о) s FUo(x0, Vo)

at z0 = 0—f—, for (x0, y0) on »S’o, where <S0 is an area of the x0, уо-ріапе whose lateral dimensions are the same as the original planform projection S, but which is stretched chordwise by a factor l/13. (See Fig. 7-2.)

Equation (7-3) and the equiva­lent form for the lower surface state, however, that the “equivalent ” wing in zero-M, constant-density flow has (at corresponding stations) the same thickness ratio t, fractional camber

в, and angle of attack a as the origi – ——–

nal wing in the compressible stream. _____

Fig. 7-2. Equivalent wing planform in zero-Mach-number flow. If sweep is present, tanAo = (1 //3) tan A. The aspect ratio is А о = /ЗА.

The similarity law might be abbreviated

where the semicolon is used to separate the independent variables from the parameters.

By way of physical explanation,[6] Gothert’s extended Prandtl-Glauert law states that to every subsonic, compressible flow over a thin wing there exists an equivalent flow of constant density liquid (at the same flight speed and free-stream ambient conditions) over a second wing, obtained from the first by a chordwise stretching 1//3 without change of surface slope distribution. It is obvious from (5-31) that pressure coefficients at corresponding points in the two flows are related by


Since they are all calculated from similar dimensionless chordwise and spanwise integrations of the Cp-distribution, quantities like the sectional lift and moment coefficients Ci(y), Cm{y), the total lift and moment co­efficients Cl, Cm, and the lift-curve slope дСь/да are found from their constant-density counterparts by the same factor 1//3 as in (7-5). It is of interest in connection with spanwise load distribution, however, that the total lift forces and running lifts per unit //-distance are equal on the two wings, because of the increased chordwise dimensions at M = 0.

Unfortunately, when one is treating a given three-dimensional configura­tion, the foregoing transformation requires that a different planform be analyzed (or a different low-speed model be tested) for each flight Mach number at which loading data are needed. This is not true for two-dimen­sional airfoils, since then the chordwise distortion at fixed a, etc., is no more than a change of scale on an otherwise identical profile; we have already seen (Section 1-4) that such a change has no effect on the physical flow quantities at fixed M.

Measurements like those of Feldman (1948) correlate with the Gothert – Prandtl-Glauert law rather well up to the vicinity of critical Mach number, where sonic flow first appears at the wing surface. They also verify what we shall see later theoretically, that the coefficient of induced drag should be unaffected by M-changes below Afcrit. There exist, of course, more accurate compressibility corrections based on nonlinear considerations which are successful up to somewhat higher subsonic M.

Inasmuch as (5-29) applies also to small-perturbation supersonic flow, M > 1, one might suspect that the foregoing considerations could be extended directly into that range. This is an oversimplification, however, since the boundary conditions at infinity undergo an essential change— disturbances are not permitted to proceed upstream but may propagate only downstream and laterally in the manner of an outward-going sound wave. (The behavior is connected with a mathematical alteration in the nature of the partial differential equation, from elliptical to hyperbolic or “wavelike.”) What one does discover is the existence of a convenient reference Mach number, M = /2, which plays a role similar to M = 0 in the subsonic case. When M = /2, the quantity В = y/M2 – 1 becomes unity and all flow Mach lines are inclined at 45° to the flight direction. Repetition of the previous reasoning leads to a supersonic similarity law

<p(x, V, 2; M, А, г, 0, a) = <P (jj ‘ V, 2; M’ = л/2, BA, T, 0, a) . (7-6)

Pressure coefficients at corresponding points, lift coefficients, etc., are related by

CP = (Cp)*_* (7-7)

Once more the equivalent planform at M = /2 is obtained from the original by chordwise distortion, but now this involves a stretching if the original M < /2 and a shrinking if M > y/2. The process has been likened to taking hold of all Mach lines and rotating them to 45°, while chordwise dimensions vary in affine proportion.

Clearly, Eqs. (7-5), (7-7), and the associated transformation techniques fail in the transonic range where M ^ 1. It has been speculated, because the equivalent aspect ratio approaches zero as M —> 1 and slender-body – theory results for lift are independent of Mach number (Chapter 6), that linearized results for three-dimensional wings might be extended into this range. This is, unfortunately, an oversimplification. Starting from the proper, nonlinear formulation of transonic small-disturbance theory, Chapter 12 derives the actual circumstances under which linearization is permissible and gives various similarity rules. It is found, for instance, that loading may be estimated on a linearized basis whenever the param­eter At1/3 is small compared to unity.

Irrotational Flow

Enough has been said about the subject of vorticity, its conservation and generation, that it should be obvious that an initially irrotational, uniform, inviscid flow will remain irrotational in the absence of heat transfer and of strong curved shocks. One important consequence of per­manent irrotationality is the existence of a velocity potential. That is, the equation

f = VxQ = 0 (1-55)

is a necessary and sufficient condition for the existence of a potential Ф such that

Q = V<f>, (1-56)

where Ф(х, у, г, t) or Ф(г, t) is the potential for the velocity in the entire flow. Its existence permits the replacement of a three-component vector by a single scalar as the principle dependent variable or unknown in theo­retical investigations.

Given the existence of Ф, we proceed to derive two important conse­quences, which will be used repeatedly throughout the work which follows.

1. The Bernoulli Equation for Irrotational Flow (Kelvin’s Equation).

This integral of the equations of fluid motion is derived by combining (1-3) and (1-28), and assuming a distant acting force potential:

DQ dQ IQ2′

Dt dt 2 у dt 2 у



Under our present assumptions,





„ /эф

= VW’


so that (1-57) can be rearranged into


^ – ol = 0.

p J


The vanishing of the gradient implies that, at most, the quantity involved will be a function of time throughout the entire field. Hence the least restricted form of this Bernoulli equation is

f + T + ly-11-™- <1Hi0)

In all generality, the undetermined time function here can be eliminated by replacing Ф with

Ф’ = Ф — jF(t) dt. (1-61)

This artificiality is usually unnecessary, however, because conditions are commonly known for all time at some reference point in the flow. For instance, suppose there is a uniform stream I7„ at remote points. There Ф will be constant and the pressure may be set equal to p„ and the force potential to SI» at some reference level.

F(t) = + j ~ — Q, * = const. (1-62)

The simplified version of (1-60) reads


~ + m2 – Vl]+ ^ + [0. – Q] = 0. (1-63)

dt JPtx p

= V — P°° ip*>Ul

2 П~ 7-1 /дФ

~ УМ* IL1 al " dt

Irrotational Flow
image15 image16

In isentropic flow with constant specific heat ratio 7, (1-63) is easily reorganized into a formula for the local pressure coefficient

Подпись: where image18,image19 Подпись: (1-65) (1-66)

Here a is the speed of sound and M = Ux/ax is Mach number. For certain other purposes, it is convenient to recognize that

(Note that, here and below, the particular form chosen for the barotropic relation is isentropic. Under this restriction, dp/dp, (dp/Sp),, and a2 all have the same meaning.)

This substitution in (1-63) provides a convenient means of computing the local value of a or of the absolute temperature T,

a2 – al = -(7 – 1) + i(Q2 – ul) + (Q. – Q)j • (1-67)

Finally, we remind the reader that the term containing the body-force potential is usually negligible in aeronautics.

2. The Partial Differential Equation for Ф. By substituting for p and Q in the equation of continuity, the differential equation satisfied by the velocity potential can be developed:

lf + V. Q_0. (Ь68)

The second term here is written directly in terms of Ф, as follows:

V ■ Q = V • (УФ) = V4,

і V2* =

Irrotational Flow Подпись: Э2Ф dz2 : Подпись: in Cartesians J • Подпись: (1-69)

which we identify as the familiar Laplacian operator,

To modify the first term of (1-68), take the form of Bernoulli’s equation appropriate to uniform conditions at infinity, for which, of course, a special case would be that of fluid at rest, U«, = 0. The body-force term is dropped for convenience, leaving

[ ~ ~ – UQ2 – ul). (1-70)

j Poo P 01

By the Leibnitz rule for differentiation of a definite integral,

Подпись:d fP dp __ 1 dP h* P ~ P

We then apply the substantial derivative operator to (1-71) and make use of the first three members of (1-66),

Подпись: Dp dP Jvoo P J Dt ГР г – ГР ■

dp _ d I dp

Подпись: a2 Dp p Dt Подпись: (1-72) Подпись: 1 dp Dp p dp Dt

Jpc p ~ Up Jp« p.

In view of (1-67) and of the simple relationship between the velocity vector and Ф, this is essentially the desired differential equation. If it is multiplied through by a2, one sees that it is of third degree in the unknown dependent variable and its derivatives. It reduces to an ordinary wave equation in a situation where the speed of sound does not vary signifi­cantly from its ambient values, and where the squares of the velocity components can be neglected by comparison with a2.

It is of interest that Garrick (1957) has pointed out that (1-74) can be reorganized into

v2*=Ml+о -v) (f+«• – v#) – h Ш *■ (i-re>

where the subscript conQc and on the substantial derivative is intended to indicate that this velocity is treated as a constant during the second application of the operators d/dt and (Q • V). Equation (1-75) is just a wave equation (with the propagation speed equal to the local value of a) when the process is observed relative to a coordinate system moving at the local fluid velocity Q.

The question of boundary conditions and the specialization of (1-74) and (1-75) for small-disturbance flows will be deferred to the point where the subject of linearized theory is first taken up.

1-8 The Acceleration Potential

image22,image23 Irrotational Flow

It is of interest that when the equations of fluid motion can be simplified to the form

In a manner paralleling the treatment of irrotational flow, we can conclude that

Подпись: (1-78)a = V¥,

where ¥ (r, t) is a scalar function called the acceleration potential. Clearly,

Подпись:¥ = П – ^ + G(t),


G(t) being a function of time that is usually nonessential.

Подпись: ІЕ I Pm p Подпись: Pm Подпись: (1-80)

The acceleration potential becomes practically useful when disturbances are small, so that

In the absence of significant body forces, we then have

Подпись: (1-81)– V


This differs only by a constant from the local pressure, and doublets of ¥ prove a very useful tool for representing lifting surfaces. The authors have not been able to construct a suitable partial differential equation for the acceleration potential in the general case, but it satisfies the same equation as the disturbance velocity potential in linearized theory.

Qualitative Effects of Viscosity

It is a common feature of most flows of engineering interest that the viscosity of the fluid is extremely small. The Reynolds number, Re = TJJ,/v, which gives an overall measure of the ratio of inertia forces to viscous forces, is in typical aeronautical applications of the order 10® or more. For large ships, Reynolds numbers of the order 109 are common. Viscosity can then only produce significant forces in regions of extremely high shear, i. e., in extremely thin shear layers where there is a substantial variation of velocity across the streamlines. The thickness of the laminar boundary layer on a flat plate of length l is approximately ~ o///Re. For Re = 108, 8/1 ~ 0.005, which is so thin that it cannot even be illustrated in a figure without expanding the scale normal to the plate. A turbulent boundary layer has a considerably greater thickness, reflecting its higher drag and therefore larger momentum loss; an approximate formula given in Schlichting (1960), p. 38, is 8/1 ~ 0.37/Re1/s. For


Re = 10® this gives 6/1 = 0.023, which is still rather small. Viscosity is only important in a very small portion of the turbulent boundary layer next to the surface, in the “viscous sublayer. ”

The transition from a laminar to a turbulent boundary layer is a very complicated process that depends on so many factors that precise figures for transition Reynolds numbers cannot be given. For a flat plate in a very quiet free stream (i. e., one having a rms turbulent velocity fluctua­tion intensity of 0.1% or less) transition occurs at approximately a distance from the leading edge corresponding to Re ~ З X 10®. With a turbulent intensity of only 0.3% in the oncoming free stream the transition Reynolds number decreases to about 1.5 X 10®. These distances are far beyond that for which the boundary layer first goes unstable. Stability calculations show that on a flat plate this occurs at Re ~ 10s. The complicated series of events between the point where instability first sets in until transition occurs has only recently been clarified (see Klebanoff, Tidstrom, and Sargent, 1962). Transition is strongly influenced by the pressure gradient in the flow; a negative (“favorable”) pressure gradient tends to delay it and a positive (“adverse ”) one tends to make it occur sooner. As a practical rule of thumb one can state that the laminar boundary layer can only be maintained up to the point of minimum pressure on the airfoil. On a so-called laminar-flow airfoil one therefore places this point as far back on the airfoil as possible in order to try to achieve as large a laminar region as possible. Laminar-flow airfoils work successfully for moderately high Reynolds numbers (<107) and low lift coefficients but require ex­tremely smooth surfaces in order to avoid premature transition. At very high Reynolds numbers transition starts occurring in the region of favor­able pressure gradient.

Both laminar and turbulent boundary layers will separate if they have to go through extensive regions of adverse pressure gradients. Separation will always occur for a subsonic flow at sharp corners, because there the pressure gradient would become infinite in the absence of a boundary layer. Typical examples of unseparated and separated flows are shown in Fig. 4-1. In the separated flow there will always be a turbulent wake

behind the body. In principle one could have instead a region of fairly quiescent flow in the wake, separated from the outer flow by a thin laminar shear layer attached to the laminar boundary layer on a body. However, a free shear layer, lacking the restraining effect of a wall, will be highly unstable and will therefore turn turbulent almost immediately. Because of the momentum loss due to the turbulent mixing in the wake the drag of the body will be quite large. On a thin airfoil at a small angle of attack the boundary layer will separate at the sharp trailing edge but there will be a very small wake so that a good model for the flow is the attached flow with the Kutta condition for the inviscid outer flow determining the circulation.

Constant-Density Flow; the Thickness Problem

Having shown how steady, constant-density flow results are useful at all subcritical M, we now elaborate them for the finite wing pictured in Fig. 7-1. As discussed in Section 5-2 and elsewhere above, it is convenient to identify and separate portions of the field which are symmetrical and
antisymmetrical in z, later adding the disturbance velocities and pressures in accordance with the superposition principle. The separation process involves rewriting the boundary conditions (5-30) as

«* = Tfx + etx~a at г = 0+ ,

do. „dh –

<pz — — Гу + в——– a at z = 0—

dx dx

where h(x, y) is proportional to the ordinate of the mean camber surface, while 2i? (ж, у) is proportional to the thickness distribution. The differential equation is, of course, the three-dimensional Laplace equation

V2<p = 0. (7-9)

As a starting point for the construction of the desired solutions, we adapt (2-28) to express the perturbation velocity potential at an arbitrary field point (ж, у, z),

Here n is the normal directed into the field, and the integrals must be carried out over the upper and lower surfaces of S. Dummy variables (xi, 2/i, zi) will be employed for the integration process, so the scalar distance is properly written

r = V(x — xi)2 + (y — yi)2 + (z — zi)2. (7-11)

In wing problems Zj = 0 generally.

Considering the thickness alone, we have

ip{x, y, z) = <p(X, y, —z)

for all z and the boundary condition

Moreover, no discontinuities of <p or its derivatives are expected anywhere else on or off the ж, y-plane. In (7-10), dS = dzj dyb the values of <p(xi, //і, 0+) and <p(xi, yi,0—) appearing in the integrals over the upper and lower surfaces are equal, while the values of д/дп(1/4ят) are equal and opposite. Hence the contributions from the first term in brackets cancel,

where (7-13) and (7-11) have been – employed. Physically, (7-14) states that the flow due to thickness can be represented by a source sheet over the planform projection, with the source strength per unit area being pro­portional to twice the thickness slope dg/dx. [Compare the two-dimensional counterpart, (5-50).]

Examination of (7-14) leads to the conclusion that the thickness prob­lem is a relatively easy one. In the most common situation when the shape of the wing is known and the flow field constitutes the desired information, one is faced with a fairly straightforward double integration. For certain elementary functions g(x, y) this can be done in closed form; otherwise it is a matter of numerical quadrature, with careful attention to the pole singularity at aq = x, у = y, when one is analyzing points on the wing z = 0. The pressure can be found from (5-31) and (7-14) as

There is no net loading, since Cp has equal values above and below the wing. Also the thickness drag works out to be zero, in accordance with d’Alem­bert’s paradox (Section 2-5). Finally, it should be mentioned that, for any wing with closed leading and trailing edges,

avydaq = Рте — !7le — 0. (7-16)

J chord OX

This means that the total strength of the source sheet in (7-14) is zero. As a consequence, the disturbance at long distances from the wing ap­proaches that due to a doublet with its axis oriented in the flight direction, rather than that due to a point source.

Constant-Density Inviscid Flow

1- 1 Introduction

For the present chapter we adopt all the limitations listed in Section 1-1, plus the following:

(1) p = constant everywhere; and

(2) The fluid was initially irrotational.

The former assumption implies essentially an infinite speed of sound, while the latter guarantees the existence of a velocity potential. Turning to (1-74), we see that the flow field is now governed simply by Laplace’s equation

V4 = 0. (2-1)

Associated with this differential equation, the boundary conditions prescribe the values of the velocity potential or its normal derivative over the surfaces of a series of inner and outer boundaries. These conditions may be given in one of the following forms:

(1) The Neumann problem, in which vn = дФ/дп is given.

(2) The Dirichlet problem, in which the value of Ф itself is given.

(3) The mixed (Poincare) problem, in which Ф is given over certain portions of the boundary and дФ/дп is given over the remainder.

A great deal is known about the solution of classical boundary-value problems of this type and more particularly about fluid dynamic applica­tions. Innumerable examples can be found in books like Lamb (1945) and Milne-Thompson (1960). The subject is by no means closed, how­ever, as will become apparent in the light of some of the applications presented in later sections and chapters. A useful recent book, which combines many results of viscous flow theory with old and new develop­ments on the constant-density inviscid problem, is the one edited by Thwaites (1960).

Among the very many concepts and practical solutions that might be considered worthy of presentation, we single out here a few which are especially fundamental and which will prove useful for subsequent work.


2- 2 The Three-Dimensional Rigid Solid Moving Through a Liquid

Let us consider a single, finite solid S moving through a large mass of constant-density fluid with an outer boundary 2, which may for many purposes be regarded as displaced indefinitely toward infinity. For much of what follows, S may consist of several three-dimensional solids rather than a single one. The direction of the normal vector n will now be re­garded as from the boundary, either the inner boundary S or the outer boundary 2, toward the fluid volume.

In the early sections of Chapter 3, Lamb (1945) proves the following important but relatively straightforward results for a noncirculatory flow, which are stated here without complete demonstrations:

(1) The flow pattern is determined uniquely at any instant if the boundary values of Ф or дФ/дп are given at all points of S and 2. One important special case is that of the fluid at rest remote from <S; then Ф can be equated to zero at infinity.

(2) The value of Ф cannot have a maximum or minimum at any interior point but only on the boundaries. To be more specific, the mean value of Ф over any spherical surface containing only fluid is equal to the value at the center of this sphere. This result is connected with the interpreta­tion of the Laplacian operator itself, which may be regarded as a measure of the “lumpiness” of the scalar field; Laplace’s equation simply states that this “lumpiness ” has the smallest possible value in any region.

(3) The magnitude of the velocity vector Q = |Q| cannot have a max­imum in the interior of the flow field but only on the boundary. It can have a minimum value zero at an interior stagnation point.

(4) If Ф = 0 or дФ/дп = 0 over all of S and 2, the fluid will be at rest everywhere. That is, no boundary motion corresponds to no motion in the interior.

Next we proceed to derive some less straightforward results.

1. Green’s Theorem. For the moment, let the outer boundary of the flow field remain at a finite distance. We proceed from Gauss’ theorem for any vector field A, (1-14),

j^A-ndS = – jjjv ■ AdV. (2-2)


(The minus sign here results from the reversal of direction of the normal vector.) Let Ф and Ф’ be two continuous functions with finite, single­valued first and second derivatives throughout the volume V. We do not yet specify that these functions represent velocity potentials of a fluid flow. Let


Constant-Density Inviscid Flow

Equations (2-4) and (2-5) are now substituted into Gauss’ theorem. After writing the result, we interchange the functions Ф and Ф’, obtaining two alternative forms of the theorem:

2. Kinetic Energy. As a first illustration of the application of Green’s theorem, let Ф in (2-6) be the velocity potential of some flow at a certain instant of time and let Ф’ = Ф. Of course, it follows that

Подпись: (2-8)у2Ф = VV = 0.

We thus obtain a formula for the integral of the square of the fluid particle speed throughout the field

“ -///lvt>’dr – "///«*"• 0«>

S+2 V V

Moreover, if we multiply the last member of (2-9) by one-half the fluid density p and change its sign, we recognize the total kinetic energy T of the fluid within V. There results


Such an integral over the boundary is often much easier to evaluate than a triple integral throughout the interior of the volume. In particular, ЭФ/дп is usually known from the boundary conditions. If the solid S is moving through an unlimited mass of fluid, with Q = 0 at infinity, it is a simple matter to prove that the integral over S vanishes. It follows that (2-10) need be integrated over only the inner boundary at the surface of the solid itself,


3. A Reciprocal Theorem. Another interesting consequence of Green’s theorem is obtained by letting Ф and Ф’ be the velocity potentials of two different constant-density flows having the same inner and outer bounding surfaces. Then, of course, the two Laplacians in (2-6) and (2-7) vanish, and the right-hand sides of these two relations are found to be equal. Equating the left-hand sides, we deduce


S+S S+2

Constant-Density Inviscid Flow Подпись: Cv ~ P«) +^Q2 + P(O« - o) Подпись: (2-13)

4. The Physical Interpretation of Ф. To assist in understanding the significance of the last two results and to give a meaning to the velocity potential itself, we next demonstrate an artificial but nevertheless meaning­ful interpretation of Ф. We begin with Bernoulli’s equation in the form (1-63), assuming the fluid at infinity to be at rest and evaluating the pres­sure integral in consequence of the constancy of p,

Imagine a process in which a system of very large impulsive pressures P = Ґ p dr (2-14)

is applied to the fluid, starting from rest, to produce the actual motion existing at a certain time t. In (2-14), т is a dummy variable of integration. We can make the interval of application of this impulse arbitrarily short, and integrate Bernoulli’s equation over it. [Incidentally, the same result is obtainable from the basic equations of fluid motion, (1-3) and (1-4), by a similar integration over the interval dt.] In the limit, the integrals of the pressure p„, pQ2/2, and p(Sl„ — fi) become negligible relative to that of the very large p, and we derive

pf = “p – (2-15)

J u-bt) от

Hence, P = — рф constitutes precisely the system of impulsive pressures required to generate the actual motion swiftly from rest. This process might be carried out, for example, by applying impulsive force and torque to the solid body and simultaneously a suitable distribution of impulsive pressures over the outer boundary 2. The total impulse thus applied will equal the total momentum in the instantaneous flow described by Ф. Unfortunately, both this momentum and the impulse applied at the outer boundary become indeterminate as 2 spreads outward toward infinity, so

that there are certain problems of physical interpretation when dealing with an externally unbounded mass of liquid.

In the light of this interpretation of the velocity potential, we reexamine the kinetic energy in (2-11), rewriting this result

Подпись: S+2Подпись: S+2image25(2-16)

The work done by an impulse acting on a system which starts from rest is known to be the integral over the boundary of the product of the impulse by one-half the final normal velocity at the boundary. Thus, the starting impulses do a total amount of work given by exactly the last member of (2-16). Since the system is a conservative one, this integral would be expected to equal the change of kinetic energy which, of course, is the final kinetic energy T in the present case. The difficulty in connection with carrying 2 to infinity disappears here, since the work contribution at the outer boundary can be shown to approach zero uniformly. Hence, the kinetic energy of an unbounded mass of constant-density fluid without circulation can be determined entirely from conditions at the inner bound­ary, and it will always be finite if the fluid is at rest at infinity.

The reciprocal theorem, (2-12), can be manipulated, by multiplication with the density, into the form


As such, it becomes a special case of a fairly familiar theorem of dynamics which states that, for any two possible motions of the same system, the sum over all the degrees of freedom of the impulse required to generate one motion multiplied by the velocity in the second motion equals the same summed product taken with the impulses and velocities interchanged.

One final important result is stated without proof: it can be shown that for a given set of boundary conditions the kinetic energy Г of a liquid in a finite or infinite region is a minimum when the flow is acyclical and irrotational, relative to all other possible motions.