Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V...

23
1 Applied Aerodynamics Def: Mach Number (M), a V M = Ratio of flow velocity to the speed of sound Compressibility Effects Def: Reynolds Number (Re), = µ ρ c V R e Ratio of inertial forces to viscous forces Viscous Effects If R is a reaction force acting on a body and C R is its dimensionless coefficient: S V R C R 2 2 1 = ρ then dimensional analysis (Buckingham’s PI theorem) tells us ) , ( 1 = M R f C e R which implies ) , ( 2 = M R f C e L ) , ( 3 = M R f C e D ) , ( 4 = M R f C e M

Transcript of Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V...

Page 1: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

1

Applied Aerodynamics Def: Mach Number (M),

aVM =∞

Ratio of flow velocity to the speed of sound Compressibility Effects Def: Reynolds Number (Re),

∞∞=µ

ρ cVRe

Ratio of inertial forces to viscous forces Viscous Effects If R is a reaction force acting on a body and CR is its dimensionless coefficient:

SV

RCR2

21

∞∞

then dimensional analysis (Buckingham’s PI theorem) tells us

),(1 ∞= MRfC eR which implies

),(2 ∞= MRfC eL ),(3 ∞= MRfC eD ),(4 ∞= MRfC eM

Page 2: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

2

Your Fluid Mechanics Courses: Basic Fluid Mechanics

• Application of F=ma to fluids • Imcompressible Flows • Mostly Inviscid (some pipe flow if you were lucky)

Fundamental Aerodynamics

• Incompressible Flow Theory (no compressibility, no viscousity)

• Sources, sinks and vortices • Panel Methods (2D only)

Gas Dynamics

• Mach Number effects (mostly shock/expansion theory, supersonic flows)

• Possible corrections for high subsonic flow Applied Aerodynamics

• 3D Bodies • Viscous Flows • Compressibility Effects • Aerodynamic Design Issues

Inviscid 3D Wing Theory Viscous Airfoil Codes – Xfoil from MIT

Page 3: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

3

Class Plan: • Review 2D Airfoil Theory (Chapters 3&4 of Anderson) • Introduce 3D Wing Theory (Chapters 5&6) • Navier-Stokes and Boundary Layer Equations • XFOIL • Applied Aerodynamics Design Concepts (along the way)

Page 4: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

4

Review of Airfoil Theory and 2D Aerodynamics Elemental Solutions Flows about 2D airfoils are built from four elemental solutions: uniform flow, source/sink, doublet and vortex. These elemental solutions are solutions to the governing equations of incompressible flow, Laplace’s equation.

02 =∇ φ 02 =∇ ϕ (3.1) which relies on the flow being irrotational

0=×∇ Vr

(3.2)

Equations (3.1) are solved for N - the velocity potential

R - the stream function The velocity potential is such that

φ∇=Vr

, i.e., y

vx ∂

u ∂=

∂∂

=φφ , (3.3)

whereas the stream function is always orthogonal (perpendicular) to the velocity potential,

xv

yu

∂∂

−=∂∂

=ϕϕ , (3.4)

Page 5: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

5

Convenient Coordinate Systems Cartesian: ),,( zyxφφ =

02

2

2

2

2

22 =

∂∂

+∂∂

+∂∂

=∇zyxφφφφ (3.5)

Cylindrical: ),,( zr θφφ =

0112

2

2

2

2

2 =∂∂

+∂∂

+

∂∂

∂∂

=∇zrr

rrr

φθφφφ (3.6)

Spherical: ),,( Φ= θφφ r

0sin

1sinsinsin1 2

2

2 =

Φ∂∂

Φ∂∂

+

∂∂

∂∂

+

∂∂

∂∂

=∇φ

θθφθ

θφθ

θφ

rr

rr (3.7) Big Advantage: Laplace’s Equation is LINEAR

∴ Superposition is possible

Page 6: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

6

Boundary Conditions Solution about a specific body requires that the body be a streamline, i.e., flow tangency condition.

0ˆ =⋅nVr

(3.8) or

0

0

=∂∂

=∂∂

s

φ

where n- surface normal, s – surface tangent

Another way to say it is that the surface has to be a streamline:

uv

dxdyb = streamline equation

Page 7: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

7

Uniform Flow

21 CxC +=φ 31 CyC +=ϕ (3.9)

0

1

=∂∂

−=∂∂

=

∂∂

===∂∂

= ∞

xyv

yVC

xu

ϕφ

ϕφ

In cylindrical coordinates θφ cosrV∞= θϕ sinrV∞= (3.10)

Page 8: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

8

Source/Sink

0, == θVrcVr (3.11)

Source Strength rrVlv π2=≡Λ&

(3.12)

Where v& is the volume flow rate (sometimes called Q) ρm&

& =v

So r

Vr π2Λ

= (3.13)

012

ln2

=∂∂

=

Λ=

∂∂

=

Λ=

θφ

πφ

πφ

θ rV

rrV

r

r

0

21

2

=∂∂

−=

Λ=

∂∂

=

Λ=

rV

rrVr

ϕπθ

ϕ

θπ

ϕ

θ

(3.14)

Def: Circulatio ': rn

( ) dsVdsVC∫∫∫ ⋅−=⋅×∇−=Γr

(3.15)

if S is an area that includes the source, 0sin/ =Γ ksource

Page 9: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

9

Superposition of Source and Uniform Flow

θϕπ

θθϕ

θπ

θϕ

θ sin

2cos12

sin

−=∂∂

−=

Λ+=

∂∂

=

Λ+=

Vr

V

rV

rV

rV

r (3.16)

Stagnation Point – Set 0== θVrV

Λ

ππ

,2 V

(3.17)

The stagnation point can now be used to define the body streamline, i.e., insert the stagnation point coordinates into the stream function equation to get its constant value. That new equation defines the body streamline.

Page 10: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

10

Superposition of Source, Sink and Uniform Flow

21 22sin θ

πθ

πθϕ Λ

−Λ

+= ∞rV (3.18)

Doublet Flow The source and sink of the last section were separated by a distance l and they had identical strengths. If one imagines that their stream function is variable and their distance can be altered a new condition would result if they are moved together while holding constant the product l7. In the limit, as l 0 a doublet is formed.

rK θπ

φ cos2

= r

K θπ

ϕ sin2

−= (3.19)

where K= l7. Streamlines are found by choosing R=constant so that we get

θθπ

sinsin2

dcKr =−=

(3.20)

Page 11: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

11

Doublet + Uniform Flow When a doublet is superposed together with a uniform flow we get the flow over a circular cylinder.

rKrV θπ

θϕ sin2

sin −= ∞ (3.21)

or if we let ∞

= VKR π2

2

−= ∞ 2

2

1sinrRrV θϕ

(3.22)

Velocities become:

θ

θθϕ

θ

θθϕ

θ

θ

sin1

sin12sin

cos1

1cos11

2

2

2

2

3

2

2

2

2

2

∞∞

+−=

−+−=∂∂

−=

−=

−=∂∂

=

VrRV

VrR

rRrV

rV

VrRV

rRrV

rrV

r

r

(3.23-3.24)

Page 12: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

12

Stagnation Points are again found by setting the velocities to zero.

0cos1 2

2

=

− ∞ θVrR

0sin1 2

2

=

+ ∞ θV

rR

(3.25)

Which gives (R,0) and (R,B). The stagnation streamline then becomes

01sin 2

2

=

−= ∞ rRrV θϕ

, i.e., r=R, a circle (3.26)

Page 13: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

13

Pressure Coefficient From the Bernoulli Equation we get:

2

1

−=

∞VVCp (3.27)

Which becomes on the cylinder

θ2sin41−=pC (3.28)

From this we see that the axial and normal force coefficients become:

( )

( )dyCCc

C

dxCCc

C

TE

LEbpfpa

c

uplpn

−=

−=

,,

0,,

1

1

(3.29)

Both turn out to be zero since the pressure coefficient is symmetric. , i.e., D’Alembert’s paradox!

Page 14: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

14

Vortex Flow

θ

πφ

−=

rln2π

ϕ Γ=

(3.30)

1 Γ−=

∂=V φ

θ

0

2

=∂∂

=

rV

rrφ

πθ

θ

Circulation rVdsV

C

πθ 2−=⋅−= ∫Γr

(3.31)

But we also know that ( ) dsV ⋅×∇−=Γ ∫∫

r (3.32)

How does this square if the flow is everywhere irrotational? Answer, it doesn’t, the flow is irrotational everywhere except at the center of the vortex. The angular velocity becomes infinite at that singularity.

Page 15: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

15

Doublet + Uniform Flow + Vortex

( )Rr

rRrV ln

21sin 2

2

πθϕ Γ

+

−= ∞ (3.33)

rV

rR

rV

VrR

rVr

πθϕ

θθϕ

θ 2sin1

cos11

2

2

2

2

Γ−

+−=∂∂

−=

−=∂∂

=

(3.34)

This time the stagnation points depend upon '.

RV

RV

RV

π

π

π

4

4

4

( )

Γ±

Γ

Γ−

∞∞

23,

42

23,

4sin,

22

1

πππ

π

π

RVV

R

RVR

Page 16: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

16

Γ==′=

Γ=

Γ+

Γ+−=

∞∞∞

∞∞

VSCqLC

RVC

RVRVC

l

d

l

p

ρ

ππθθ

0

2sin2sin41

22

(3.35)

Lift is directly proportional to circulation

Page 17: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

17

Panel Methods The primary technique for low speed analysis (incompressible flow) is the panel method. We will discuss two, the source panel method and the vortex panel method. Source Panel Method To begin consider what would happen if you stretched a point source/sink out over some curve in space

Defined by the equation

rdsd ln2πλφ = (3.36)

where )(sλλ = , s being the direction along the “source sheet.” To determine the potential at some point P(x,y) induced by the source sheet, one integrates Eq. 3.36 along the sheet.

∫=b

a

rdsyx ln2

),(π

λφ (3.37)

Page 18: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

18

When combined with a uniform flow this yields the flow over an arbitrary body if )(sλλ = is prescribed in such a way that the desired body shape becomes a streamline.

In general, it is difficult to find an analytical expression for

)(sλλ = , so it is approximated either by a series of constant source “panels” or by linearly varying source “panels.”

The equation for the potential induced by the jth panel is then

jj ijj

j dsr∫=∆ ln2πλ

φ (3.38)

Which gives the contribution to N at point P from panel j. The potential for the entire flow field then comes by summing the contributions from all the panels. At the center of each panel is a control point at which the basic equation for the surface normal velocity can be found from

Page 19: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

19

( )[ ] ( )∑ ∫≠= ∂

∂+=

∂∂

=n

ijj

jj

iji

jiii

in dsr

nyx

nV

i1

ln22

λλφ (3.39)

This is the normal velocity contribution at the control point from only the source sheet. However, what we want is the body to be a streamline, so we must have:

0=+∞ inn VV (3.40) This leads to the equation that can be solved for the jλ unknowns

0cos22 1

=++ ∑≠=

∞ i

n

ijj

ijji VI β

πλλ

(3.41)

where I is the integral in Eq. 3.39. The surface velocity on the body at the control point can then be found from

( )[ ] ( )∑ ∫= ∂

∂=

∂∂

=n

jj

jij

jiis dsr

syx

sV

i1

ln2

λφ (3.42)

once the jλ are found. Lastly, the integral can be evaluated as on page 256 of the Anderson text and C found from the Bernoulli equation.

p

Unfortunately, since we are dealing only with sources and sinks the lift and drag are both zero!!!

Page 20: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

20

Lifting Airfoils The problem of lift in general is bound up in the idea of circulation and hence requires vortices in the panel method, which we showed earlier to be the only elemental solution to exhibit circulation. To accommodate this we use an idea similar to the source sheet, the vortex sheet

Whose potential is given by the equation

( ) ∫−=b

a

dszx θγπ

φ21, (4.1)

where )(sγγ = . Circulation comes about from the vortices and is related directly to the distribution of γ , via:

∫=Γb

a

dsγ (4.2)

From this point the entire thin airfoil theory can be developed (which should be reviewed) or one can follow the same approach

Page 21: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

21

used to develop the source panel method using the vortex sheet equation to produce a vortex panel method for arbitrary lifting bodies. However, one thing missing is the specification of stagnation points.

Which requires the Kutta Condition. Kutta Condition The basic idea is related to the fact that an infinite number of stagnation point locations can exist on a lifting airfoil if only the earlier theory is applied, and because of that, virtually any circulation (lift) is possible. However, Wilhelm Kutta observed that the flow leaves a pointed trailing edge airfoil smoothly from the trailing edge, placing the stagnation point at the trailing edge. This is the so-called Kutta Condition and it provides the additional information we need to determine the circulation and hence lift of an arbitrary airfoil.

Page 22: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

22

Start-Up Vortex Another interesting, yet apparently contradictory, idea associated with airfoils comes about because of Kelvin’s Theorem, which says circulation does not change along a fixed fluid element. If it starts out zero around an airfoil at rest, it must remain zero along that fluid element.

For circulation to exist about the airfoil (and hence lift) it must coexist with an equal strength negative circulation shed at airfoil start-up and convecting downstream bounded by the original fluid element.

Page 23: Applied Aerodynamics - University of Cincinnati · 1 Applied Aerodynamics Def: Mach Number (M), a V M∞ = Ratio of flow velocity to the speed of sound Compressibility Effects Def:

23

Vortex Panel Method So just as in the source panel method, the vortex sheet can be approximated by a series of discrete vortex panels of constant or varying strength. Again the body must be a streamline and hence the flow tangency condition can be applied at the control points. The resulting equation for vortex strength is

02

cos1

=∂∂

− ∑ ∫=

n

jj

j i

ijji ds

nV

θπ

γβ (4.3)

This produces a system of n unknowns and n equations. We think “good” but not so, because it still doesn’t incorporate the Kutta Condition. A way to do this is to make the trailing edge panels very small and enforce

1−−= ii γγ (4.4) as illustrated in the Figure below.

Unfortunately, now the system is over determined. To obtain a determinate system one must eliminate one of the panels from the linear equation set. Then, once the system is solved, the lift can be determined directly from the known coefficients, i.e.,

∑=

∞∞=′n

jjj sVL

1γρ (4.5)

Next up…Wing Theory!!!!