Document 7154188

Download Report

Transcript Document 7154188

External Flows CEE 331 April 29, 2020

Overview  Non-Uniform Flow  Boundary Layer Concepts  Viscous Drag  Pressure Gradients: Separation and Wakes  Pressure Drag  Shear and Pressure Forces  Vortex Shedding

Non-Uniform Flow  In pipes and channels the velocity distribution was uniform (beyond a few pipe diameters or hydraulic radii from the entrance or any flow disturbance)  In external flows the boundary layer is always growing and the flow is non uniform

Boundary Layer Concepts  Two flow regimes  Laminar boundary layer  Turbulent boundary layer  with laminar sub-layer  Calculations of  boundary layer thickness  Shear (as a function of location on the surface)  Drag (by integrating the shear over the entire surface)

Flat Plate: Parallel to Flow U y t o U U U boundary layer thickness d x shear Why is shear maximum at the leading edge of the plate?

du

is maximum

dy

Laminar Boundary Layer: Shear and Drag Force d

x

= 5 Re

x

Re

x

=

Ux

n d = 5 n

x U

t 0  0 .

332 

U

3

x F d

w

0 

l

t 0

dx

w

0 

l

0 .

332 

U

3

dx x

Based on momentum and mass conservation and assumed velocity distribution Integrate along length of plate

F d

 0 .

664

w



U

3

l

On one side of the plate!

Laminar Boundary Layer: Coefficient of Drag

F d

 0 .

664

w



U

3

l

C

d

= r 2F

D

2

U A

=

f

(Re) C

d

 2 ( 0 .

664 )

w

U

2

lw



U

3

l

Dimensional analysis C

d

 1 .

328

w

U

 2

lw U

3

l

C

d

 1 .

328 

Ul

 C

d

= 1.328

Re Re

l

=

Ul

n

Transition to Turbulence  The boundary layer becomes turbulent when the Reynolds number is approximately 500,000 (based on length of the plate)  The length scale that really controls the transition to turbulence is the Re

x

=

Ux

n Re d =

U

d n Re d Re

x

= d

x

= d

x

= 5 Re

x

Re d = 5 Re

x Re

d = 3500

U y x t o Transition to Turbulence U U d U turbulent This slope (du/dy) controls t 0 .

Viscous sublayer Transition (analogy to pipe flow)

Turbulent Boundary Layer: (Smooth Plates) d

x

= 0.37

Re

x

1/ 5 Re

x

=

Ux

n d  0 .

37

x

4 / 5   

U

  1 / 5 t 0

F d

  0 .

029 

U

x 5/4 2   

Ux

  1 / 5

w

l

0 t 0

dx

 Derived from momentum conservation and assumed velocity distribution 0 .

036 

U

2

wl

  

Ul

  1 / 5 Integrate shear over plate C

d

= 0.072 Re

l

1/ 5 5 x 10 5 <

Re

l

< 10 7 C

d

= r 2F

D

2

U A

=

f

ж Re, и e

l

ц ш

Boundary Layer Thickness  Water flows over a flat plate at 1 m/s. Plot the thickness of the boundary layer. How long is the laminar region?

0.1

10,000,000 laminar 0.09

9,000,000 turbulent Re

x

0.08

Reynolds Number 8,000,000 =

Ux

n 0.07

7,000,000 0.06

0.05

d  0 .

37

x

4 / 5    

U

   1 / 5 6,000,000 5,000,000

x

= n Re

x U

0.04

0.03

4,000,000 3,000,000

x

 1

x

10  6

m

2 /

s

( 500 , 000 ) 1

m

/

s

0.02

d = 5 n

U x

2,000,000 x = 0.5 m 0.01

1,000,000 0 0 1 2 3 length along plate (m) 4 5 6 -

Grand Coulee

Flat Plate Drag Coefficients

C Df

0.01

C Df

= [ ( e /

l

) ] 2.5

0.001

C Df

= 1.328

( Re

l

) 0.5

C Df

= [ 0.455

(

l

) ] 2.58

1700 Re

l

1e +04 1e +05 1e +06 1e +07 Re

l

=

Ul

n 1e +08 1e +09 1e +10 1 x 10 -3 5 x 10 -4 2 x 10 -4 1 x 10 -4 5 x 10 -5 2 x 10 -5 1 x 10 -5 5 x 10 -6 2 x 10 -6 1 x 10 -6

C Df

e

l

= 0.072 Re

l

0.2

C Df

= [ 0.455

l

) ] 2.58

Example: Solar Car  Solar cars need to be as efficient as possible. They also need a large surface area for the (smooth) solar array. Estimate the power required to counteract the viscous drag at 40 mph  Dimensions: L: 5.9 m W: 2 m H: 1 m  Max. speed: 40 mph on solar power alone  Solar Array: 1200 W peak  air = 14.6 x10 -6 m 2 /s  air = 1.22 kg/m 3

Viscous Drag on Ships  The viscous drag on ships can be calculated by assuming a flat plate with the wetted area and length of the ship

F D

F viscous

F wave

C

d

 2 F

d

U

2

A F wave

F

viscous

= C

d

r 2

U A

2 L r 3

C Df

[ ( ( e e / /

l l

) ) ] ] ( ( ) ) 0.5

[ [ ( ( ) ) ] ] 2.58

1e +04 1e +05 1e +06 1e +07 Re

l

= = n n 1e +08 1e +09 1e +10 1 x 10 5 x 10 -3 -4 2 x 10 -4 1 x 10 5 x 10 2 x 10 1 x 10 5 x 10 2 x 10 1 x 10 -4 -5 -5 -5 -6 -6 -6 e e e e [ ( ( ) ) ] ] 2.58

Re

l

=

Ul

n

Separation and Wakes  Separation often occurs at sharp corners  fluid can’t accelerate to go around a sharp corner  to the free stream velocity)  (determined by the pressure in the adjacent flow)

0 U Flat Plate: Streamlines 2 1 3 4 Point 1 2 3 4

v

0 U

C p

 1 

v

2

U

2

C p

1 >0 <0

p

>p 0 >p 0

p

U p

0 2    Points outside boundary layer!

Application of Bernoulli Equation

p

1  

h

1 

v

1 2 2

g

p

2  

h

2 

v

2 2 2

g p

0  

U

2 2

g

 

p

v

2 2

g U

2 2

g

v

2 2

g

 

p

p

0  In air pressure change due to elevation is small U = velocity of body relative to fluid 1 

v U

2 2  2

p

U p

0 2 

C p

1 0.8

Flat Plate:

C p

 1 

v

2

U

>U 2

p

U p

0 2   

C p

 2   

p

U p

0 2    

U

2

C p

p

p

0 2

F d F d

  

p F d front

 

F d rear p rear

A front

0 1 0

C p

-1 -1.2

F d

 

C p front

C p rear

 

U

2 2

A F d

  0 .

8  1 .

2  

U

2 2

A

C d = 2

Drag of Blunt Bodies and Streamlined Bodies  Drag dominated by viscous  Drag dominated by pressure  Whether the flow is viscous drag dominated or pressure drag dominated depends entirely on the shape of the body. Bicycle page at Princeton

Flow and Drag: Spheres C

d

=

f

ж e и

D

, Re,

M

,

shape orientation

ц ш General relationship for submerged objects Spheres only have one shape and orientation!

C

d

 2 F

D

U

2

A

r 2F

D

2

U A

= C

d

=

f

F

D

 C

d

U

2

A

2 Where

C d

is a function of Re

How fast do particles fall in dilute suspensions?

 What are the important parameters?

 Initial conditions Acceleration due to gravity  After falling for some time...

drag  What principle or law could help us?

Newton’s Second Law...

Sedimentation: Particle Terminal Fall Velocity 

F

ma F d

W F b

 

W

p

  0

p g F b

 

p

w g F d

C D A P

w V t

2 2 

p

 4 

r

3 3

A p

 

r

2

F d W F b

p

 particle volume

A p

 particle cross sectional area

ρ p

 particle density

ρ w g

  water density accelerati on due to gravity

C D V t

  drag coefficien t particle terminal velocity C

d

 2 F

D

U

2

A

Particle Terminal Fall Velocity (continued)

F d

W

F b C D V t A P

2  

w V t

2 2  2

p

( 

C D

p

A P

p

(  

w

w p

 )

g

w

)

g

General equation for falling objects 

p

 2

d

Relationship valid for spheres

A p V t

2  4

gd

3

C D

 3 

p

 

w

w

V t

 4

gd

3

C D

 

p

 

w

w

Drag Coefficient on a Sphere 1000 100 10 Stokes Law 1 0.1

0.1

1 24

C D

= Re 10 10 2 10 3 10 4 Reynolds Number 10 5 10 6 10 7 Re=500000 Turbulent Boundary Layer

Drag Coefficient: Equations Laminar flow R < 1

C D

= Transitional flow 1 < R < 10 4

C D

= Fully turbulent flow R > 10 4

C D

= 24 Re

V t

24 3 + Re Re 1/ 2 + 0.34

d

2

g

 

p

18   

w

 0.4

V t

»

gd

( r 0.3

p

r

w

r

w

) Re =

V d t

r m

Example Calculation of Terminal Velocity Determine the terminal settling velocity of a

cryptosporidium

oocyst having a diameter of 4  m and a density of 1.04 g/cm 3

ρ p ρ w g d

 1040 kg/m 3  9 99 kg/m 3  9.

8 1 m/s 2  4x10  6 m

V t

in water at 15°C.

V t

d

2

g

  18

p

  

w

   4x10  6 m 2   9.

8 1 m/s 2  1040 kg/m 18   1.14x10

 3 s kg  m   3  9 99 kg/m 3    1.14x10

 3 kg s  m

V t

 3 .

14

x

10  7 m/s

V t

 2 .

7 cm/day

Reynolds

Pressure Gradients: Separation and Wakes Diverging streamlines Van Dyke, M. 1982. An Album of Fluid Motion. Stanford: Parabolic Press.

Adverse Pressure Gradients   Streamlines diverge behind object Increasing pressure in direction of flow

p

Fluid is being decelerated g  than the main flow and may be completely

V

2 2

g

=

C

stopped.

 If boundary layer stops flowing then separation occurs

Point of Separation  Predicting the point of separation on smooth bodies is beyond the scope of this course.

 Expect separation to occur where streamlines are diverging (flow is slowing down)  Separation can be expected to occur around any sharp corners (where streamlines diverge rapidly)

Drag on Immersed Bodies (more shapes)  Figures 9.19-21 bodies with drag coefficients on p 392-394 in text.

 hemispherical shell  hemispherical shell  cube  parachute 0.38

1.42

1.1

1.4

Why?

Velocity at separation point determines pressure in wake Vs ?

Shear and Pressure Forces: Horizontal and Vertical Components F

D

F

L

U   

p

sin    

p

cos   t 0 cos   t 0 sin  

dA

dA

lift drag Parallel to the approach velocity lift p < p 0 Normal to the approach velocity negative pressure

F d

C d A

U

2

A

2 defined as projected 

F L

C L A

U

2 2 drag p > p 0 positive pressure

Shear and Pressure Forces  Shear forces  viscous drag, frictional drag, or skin friction  caused by shear between the fluid and the solid surface   Pressure forces  pressure drag or form drag   function of area normal to the flow

Example: Beetle Power

C d

= 0.38

Height = 1.511 m Width = 1.724 m Length = 4.089 m Ground clearance = 15 cm?

85 kW at 5200 rpm Where does separation occur?

Calculate the power required to overcome drag at 60 mph and 120 mph.

Is the beetle streamlined?

Electric Vehicles  Electric vehicles are designed to minimize drag.  Typical cars have a coefficient drag of 0.30-0.40.  The EV1 has a drag coefficient of 0.19. Smooth connection to windshield

Drag on a Golf Ball DRAG ON A GOLF BALL comes mainly from pressure drag. The only practical way of reducing pressure drag is to design the ball so that the point of separation moves back further on the ball. The golf ball's dimples increase the turbulence in the boundary layer, and delay the onset of separation. The effect is plotted in the chart, which shows that for Reynolds numbers achievable by hitting the ball with a club, the coefficient of drag is much lower for the dimpled ball.

Why not use this for aircraft or cars?

Effect of Turbulence Levels on Drag  Flow over a sphere: (a) Reynolds number = 15,000; (b) Reynolds number = 30,000, with trip wire. Causes boundary layer to become turbulent Point of separation

Effect of Boundary Layer Transition Ideal (non viscous) fluid Real (viscous) fluid: laminar boundary layer Real (viscous) fluid: turbulent boundary layer No shear!

Spinning Spheres  What happens to the separation points if we start spinning the sphere?

LIFT!

Vortex Shedding  Vortices are shed alternately from each side of a cylinder  The separation point and thus the resultant drag force oscillate  Dimensionless frequency of shedding given by Strouhal number S  S is approximately 0.2 over a wide range of Reynolds numbers (100 - 1,000,000)

S

nd U

Summary: External Flows  Spatially varying flows  boundary layer growth  Example: Spillways  Two sources of drag  shear (surface area of object)  pressure (projected area of object)  Separation and Wakes  Interaction of viscous drag and adverse pressure gradient

Solution: Solar Car F

d

= 0.455

C Df

= [

l

) ] 2.58

2 ( C F

d d

  C

d

2 F

d

U

2 

U A

2

A

2 3 )( 1.22

/ 3 ) ( 2

C Df

1e +04

F d

=14 N

( ( ) ) 0.5

[ [ [ [ ( ( ) ) ] ] 2.58

1e +05 1e +06 1e +07 Re

l l

= = 1e +08

R x

( ( e e / /

l l

) ) ] ] 1e +09 1e +10

Ul

1 x 10 5 x 10 -3 -4 2 x 10 -4 1 x 10 5 x 10 2 x 10 1 x 10 5 x 10 2 x 10 1 x 10 -4 -5 -5 -5 -6 -6 -6 e e e e

[ [ ( ( ) ) ] ] 2.58

U = 17.88 m/s

l

= 5.9 m  air  air = 14.6 x 10 -6 m 2 /s = 1.22 kg/m 3 ) 2 ( 11.88

m

2 )

Re l

= 7.2 x 10 6

A C d

= 3 x 10 -3 = 5.9 m x 2 m = 11.8 m 2

P

=F*U=250 W

Reynolds Number Check

R

Vd

 

R

  3 .

14

x

10  7 m/s  4

x

10  6 m  999 kg/m 3  1.14x10

 3 kg s  m R = 1.1 x 10 -6 R<<1 and therefore in Stokes Law range

Solution: Power a New Beetle at 60 mph C

d

 2 F

D

U

2

A

 F

D

 C

d

U

2

A

2

f

(

R

)

P

 C

d

U

3

A

2

P

 ( 0 .

38 )( 1 .

2

kg

/

m

3 )( 26 .

82

m

/

s

) 3 ( 2 .

0

m

2 ) 2 P = 8.8 kW at 60 mph P = 70 kW at 120 mph

Solution: Power a Plymouth Voyager at 60 mph C

d

 2 F

D

U

2

A

 F

D

 C

d

U

2

A

2

f

(

R

) Frontal Area 30.15 sq ft (2.80 sq m) Drag Coefficient 0.35

P

 C

d

U

3

A

2 How can

C d

larger object?

be smaller for

P

 ( 0 .

35 )( 1 .

2

kg

/

m

3 )( 26 .

82

m

/

s

) 3 ( 2 .

80

m

2 ) 2 P = 11.3 kW at 60 mph

Grand Coulee Dam Turbulent boundary layer reaches surface!

Drexel SunDragon IV  http://cbis.ece.drexel.edu/SunDragon/Cars.html

Vehicle ID: SunDragon IV (# 76) Dimensions: L: 19.2 ft. (5.9 m) W: 6.6 ft. (2 m) H: 3.3 ft. (1 m) Weight: 550 lbs. (249 kg) Solar Array: 1200 W peak; 8 square meters terrestrial grade solar cells; manf: ASE Americas Batteries: 6.2 kW capacity lead-acid batteries; manf: US Battery Motor: 10 hp (7.5 kW) brushless DC; manf: Unique Mobility Range: Approximately 200 miles (at 35 mph on batteries alone) Max. speed: 40 mph on solar power alone, 80 mph on solar and battery power.

Chassis: Graphite monocoque (Carbon fiber, Kevlar, structural glass, Nomex) Wheels: Three 26 in (66 cm) mountain bike, custom hubs Brakes: Hydraulic disc brakes, regenerative braking (motor)

Pressure Coefficients on a Wing

C p

 1 

v

2

U

2  2   

p

U p

0 2   

NACA 63-1-412 Flowfield Cp's

Head Loss in Bends High pressure  Head loss is a function of the ratio of the bend radius to the pipe diameter (R/D) R n Possible separation from wall  Velocity distribution returns to normal several pipe diameters downstream D

K b

varies from 0.6 - 0.9

p

+ r у ф х

V

2 R

dn

+ g

z

= Low pressure

C h b

K b V

2 2

g