Axisymmetric Flows and Slender Body Theories

Axisymmetric flows are needed to analyze flow over bodies of revolution. Extensions to bodies of revolution with fins (missiles) or slender wings (low aspect ratio) are also possible, using perturbations of axisymmetric flows. In this chapter, it is natural to use cylindrical coordinates. After the derivation of the governing equations at different Mach number regimes (subsonic, supersonic, transonic), lift and drag are calculated for standard shapes (Fig. 7.1).

7.1 Governing Equations in Cylindrical Coordinates

Assuming steady, inviscid, adiabatic, irrotational and isentropic flow with uniform upstream conditions, the governing equations, in vector notation, are

V.(pV) = 0, VxV = to = 0 (7.1)

In cylindrical coordinates, Fig. 7.2, the conservation of mass is

Подпись: d 1 d 1 d

(p Ux) + (prur) + (pue) = 0 (7.2)

dx r dr r 89

where (ux, ur, u9 ) are the velocity components in x, r and 9-directions, respectively. The vorticity components are

1 d 1 d

Mx = (rue) – (ur) (7.3)

r dr r d9

1 d d

Mr = (ux) – (u9) (7.4)

r d9 dx

© Springer Science+Business Media Dordrecht 2015

J. J. Chattot and M. M. Hafez, Theoretical and Applied Aerodynamics,

DOI 10.1007/978-94-017-9825-9_7

Notice that Vм = 0, since м = Vx V. The condition is a consequence of the vector identity: the divergence of a curl vanishes.

The density can be calculated from Bernoulli’s law, assuming isentropic conditions

Подпись: P_ Po Axisymmetric Flows and Slender Body Theories(7.6)

Axisymmetric Flows and Slender Body Theories

Axisymmetric Flows and Slender Body Theories

Fig. 7.2 Slender body of revolution and coordinate system

and

Y p V2

+ = H0 = const. (7.7)

Y – 1 P 2

In the latter equation of isoenergetic flow, the potential energy is assumed constant (and does not appear in the total energy).

Подпись: P_ P0 Подпись: (1 - ^M2(V2 Подпись: (7.8)
Axisymmetric Flows and Slender Body Theories

From the above two equations, one obtains

For the irrotational flow, the velocity can be described by the gradient of a potential function

Подпись:Подпись: (7.9)V = VФ

Axisymmetric Flows and Slender Body Theories Axisymmetric Flows and Slender Body Theories Подпись: (7.10)
Axisymmetric Flows and Slender Body Theories

The potential equation in cylindrical coordinates is given by

The tangency boundary condition on a solid surface is VФ. п = 0, where n is the unit vector normal to the surface. Together with the far field condition, Ф ^ 0, as x2 + r2 ^ to, the formulation is complete.

Axisymmetric Flows and Slender Body Theories Axisymmetric Flows and Slender Body Theories

For a body of revolution at zero angle of attack, there is no dependency on в, and the governing equation reduces to

or

Подпись: (7.13)гдФ _ дФ dR _ 1 дФ dS дг дx dx 2п дx dx

Axisymmetric Flows and Slender Body Theories Подпись: (7.14)

where S(x) = п R2(x), is the cross-sectional area of the body. The pressure coefficient is

Axisymmetric Flows and Slender Body Theories
and

respectively. Also, the pressure coefficient becomes

(7.21)

Подпись: 2 дф Cp = ф p U д x Axisymmetric Flows and Slender Body TheoriesNote that, for body of revolution, д/дв = 0, and the last term in the above equations vanishes.

For incompressible flows (M0 = 0), the fundamental solution of the Laplace equation is

Подпись: (7.22)1 1

4n Vx2 + r2

Подпись: ф.ч Подпись: 1 1 4п |lx2 + e2r2 Подпись: в2 = 1 - M02 Подпись: (7.23)

For subsonic flows (M0 < 1), using the Prandtl-Glauert transformation, the fun­damental solution becomes

Подпись: ф Подпись: 1 4n Axisymmetric Flows and Slender Body Theories Подпись: (7.24)

For a source distribution q(x) per unit length along the x-axis, the perturbation potential at a point (x, r) reads

Подпись: дф lim rur = lim r r^0 r^0 д r Axisymmetric Flows and Slender Body Theories Подпись: (7.25)

The strength of the sources is determined through imposing the boundary condition

where S(x) = n R2(x) represents the cross section of the body of revolution. Notice from mass balance, one obtains (see Fig. 7.3)

Подпись: Fig. 7.3 Mass balance involving the source intensity

q(x)Ax = 2nRAR U = UAS (7.26)

Подпись: and Подпись: Ф Axisymmetric Flows and Slender Body Theories Axisymmetric Flows and Slender Body Theories

hence

Also, integrating the tangency condition in the proximity of the x-axis (small r)

дф U, U ,

r = S'(x), ^ ф ~ S(x) ln r + g(x) (7.29)

д r 2n 2n

Подпись: g(x) Подпись: U S (x) ln Подпись: S'(§) |x Axisymmetric Flows and Slender Body Theories Подпись: (7.30)
Axisymmetric Flows and Slender Body Theories

Oswatitsch and Keune [2] provided an approximation of the singular integral leading to

Note that the term g(x) depends in general on M0. For a body with a pointed nose and ending at a point or as a cylinder (S'(l) = 0)

U. (в U fx,, U fl,,

g(x) ~ 2ПS'(x) ln 2 – __ s"(§) ln(x -§)d§ + — S (§) ln(§ – x)d§

0 x (7.31)

Подпись: Фз Подпись: 1 1 2n л/x2 - в2r2 Подпись: в2 = M2 - 1 Подпись: (7.32)

For supersonic flow (M0 > 1), the fundamental solution is different. The solution for a supersonic source is

Подпись: Ф Подпись: 1 г--в ssjm m 2n 0 y/(x - §)2 - в2r2 Подпись: (7.33)

The potential is a real number only within the Mach cone, Fig. 7.4. A source distribution will give the solution

where the integral is evaluated only for values of § < x – вr, Fig. 7.5.

Von Karman and Moore [3] analyzed this singular integral. Their results are given by

Axisymmetric Flows and Slender Body Theories

и

 

Fig. 7.4 Mach cone

 

Axisymmetric Flows and Slender Body Theories

Fig. 7.5 Traces of Mach cones in meridian plane

 

U fx S'(£) – S(x)d£ 2n J0 x – £ £

 

U, U, в

S (x) ln t + S (x)

2n 2n 2x

 

(7.34)

 

Ф

 

Thus, for supersonic flows

 

Axisymmetric Flows and Slender Body Theories
Axisymmetric Flows and Slender Body Theories

U

2n

 

(7.35)

 

g(x)

 

Axisymmetric Flows and Slender Body Theories

Leave a reply

You may use these HTML tags and attributes: <a href="" title=""> <abbr title=""> <acronym title=""> <b> <blockquote cite=""> <cite> <code> <del datetime=""> <em> <i> <q cite=""> <s> <strike> <strong>