Chapter 9: External Incompressible Viscous Flows: (Unbounded)
Chapter 9: External Incompressible Viscous Flows: (Unbounded)
VISCOUS FLOWS
*(unbounded)
CHAPTER 9: EXTERNAL INCOMPRESSIBLE
VISCOUS FLOWS
= separated
bdy layer
adverse pressure
gradient
thicker leads to separation
difficult to use theory
Re = Ux/; Re = Uc/;
laminar and turbulent boundary layers
displaced inviscid outer flow
adverse pressure gradient and separation
Boundary Layer Provides Missing Link
Between Theory and Practice
- 1904 Prandtl
Fluid Motion with Very Small Friction
2-D boundary layer equations
- 1908 Blasius
The Boundary Layers in Fluids with Little Friction
Solution for laminar, 0-pressure gradient flow
Immersed Bodies
~ wall shear stress/drag?
~ lift?
~ minimize wake
breath
FLAT PLATE ZERO PRESSURE GRADIENT
Laminar Flow
/x ~ 5.0/Rex1/2
THEORY
Turbulent Flow
Rextransition > 500,000
u(y)/U = (y/)1/7
/x ~ 0.382/Rex1/5
EXPERIMENTAL
No simple theory
for Re < 1000;
(cant assume
is thin)
ReL = Ux/
*
* *
*
*
Rex 1000
breath
FLAT PLATE ZERO PRESSURE GRADIENT: (x)
Definition:
u(x,) = 0.99 of U=U=Ue
(within 1 % of U)
+Disturbance
is at y location where u(x,y) = 0.99 U
DISPLACEMENT THICKNESS
(x)
*
DISPLACEMENT THICKNESS
Displacement Thickness
(x)
*
Definition:
= 0 (1 u/U)dy
*
is displacement of outer
*
streamlines due to boundary layer
Displacement thickness *
x=L
* 0 (1 u/U)dy 0 (1 u/U)dy
function of x!
breath
Displacement Thickness *
Definition: * = 0 (1 u/U)dy
U*w = 0 (U u)dyw 0 (U u)dyw
Show that * =
=
Displacement Thickness *
Suppose u = 0 from y = 0 to and u = Ue for y > .
Show that * =
* = 0 (1 u/Ue)dy
= 0 (1 0/Ue)dy + (1 Ue/Ue)dy
* =0 dy
=
=
Laminar flow on flat plate in uniform free stream
Blasius = 5x/(Rex1/2)
Blasius* = 1.721x/(Rex)1/2
*(x)
(x)
(x)
*
PROBLEM: Find Ue as a function
of U, D and *
2-D DUCT
Ue(x)
Ue(x) = ?
U f(y) u(y) Ue(x)
UD =
-D2 D/2 U dy{
e (-D/2) (-D/2+) (U -u)dy+
e (D/2- ) D/2 (U -u)dy}
e
*
Continuity equation (per W):
UD =
-D2 D/2 U dy {
e (-D/2) (-D/2+) (U -u)dy +
e (D/2- ) D/2 (U -u)dy}
e
* = 0 (1 u/Ue)dy 0 (1 u/Ue)dy
INLET OF DUCT
*
UD =
-D2 D/2 U dyU {
e e (-D/2) (-D/2+) (1-u/U )dy+
e (D/2- ) D/2 (1-u/U )dy}
e
-{ * + *}
UD = Ue [D *
2 ]
We see that the free stream velocity in the duct
is given by the effective decrease in the cross
sectional area due to the growth of the
boundary layers, and this decrease in area is
measured by the displacement thickness!!!
FLAT PLATE ZERO PRESSURE GRADIENT: (x)
MOMENTUM THICKNESS
Momentum Thickness
* (x)
MOMENTUM THICKNESS
Momentum Thickness
(x)
Definition:
= 0 u/U(1 u/U)dy
UU w
Summary:
= 0.99U (within 1 %)
* = 0(1 u/Ue)dy
0(1 u/Ue)dy
=0 u/Ue(1 u/Ue)dy
0 u/Ue(1 u/Ue)dy
& Easier to calculate form data
*
and more physical significance
but can be dependent on
breath
Relate to drag
Relate to Drag
u(Vn)dy(w)
= 0h Uo(-Uo)wdy + 0 u(L,y) u(L,y)wdy
(Vn)dA= 0
h
- 0Uowdy + 0u(x,y)wdy = 0
Uoh = 0u(x,y)dy
-D = - Uow 0u(x,y)dy + 0 u(x,y)2wdy
-D = - Uow 0(x)u(y)dy + 0(x) u(x,y)2wdy
D = w [0(x)Uou(x,y)dy - u(x,y)2dy]
D(x) = Uo 2w (x)
(First derived by Von Karman in 1921)
breath
me
breath
me
0.08
0.06
0.04
0.8
0.02
0 y/(x)
600 650 700 750 800 850 900 950 1000 1050 0.6
u (cm/sec)
1 0.2
0.8
0
y (cm)
0.6
0.4
0.16
0.14
0.12
0.1
0.2
turbule nt bdy laye r v e locity profile
U = 10 m/s; x1=10cm
0.08
0.06
0.04
0.02
0
6000 650 700 750 800 850 900 950 1000 1050
u(x)/U
u (cm /se c)
600 650 700 750 800 850 900 950 1000 1050
u (cm/sec)
* du/dy|wall =
Blasius developed an exact solution (but numerical integration
was necessary) for laminar flow with no pressure variation.
Blasius could theoretically predict boundary layer thickness (x),
velocity profile u(x,y)/U vs y/, and wall shear stress w(x).
wall
Turbulent
forced Laminar
natural
Turbulent
Laminar
What is wrong with this figure?
What is wrong with this figure?
Are Antarctic Icebergs Towable Arctic News Record Summer 1984; 36
U14/5 = 0.03054
U = 0.288 m/s
Rex ~ 3x108
80% of fresh water found in world
Antarctic ice
Derivation
Blasius developed an exact solution (but numerical integration
was necessary) for laminar flow with no pressure variation.
Blasius could theoretically predict boundary layer thickness (x),
velocity profile u(x,y)/U vs y/, and wall shear stress w(x).
u
mf represents x-component of momentum flux
Fsx will be composed of shear force on boundary
and pressure forces on other sides of c.v.
Surface Momentum Flux Through Side ab
X-momentum
u
Flux = w
cvuVdA
Surface Momentum Flux Through Side cd
X-momentum
Flux = u
cvuVdA
Surface Momentum Flux Through Side bc
U=Ue=U
X-momentum
Flux = u
cvuVdA
X-Momentum Flux Through Control Surface
a-b c-d
c-d b-c
IN SUMMARY
X-Momentum Equation
X-Force on Control Surface
w is unit width
p(x+dx) into page
w
Surface x-Force On Side bc
p + (dp/dx)dx
Why w and
not w (bc)?
Surface x-Force On Side bc
p + dp/dx along bc
In x-direction: [p + (dp/dx)] wd
Surface x-Force On Side ad
-(w + dw/dx]xdx)wdx
p(x)
Fx = Fab + Fcd + Fbc + Fad
d << dw << w
=
U
ad = 0
ab -cd bc
Divide by wdx
Newtons Law: F = ma
u(u/x) + v(u/y) = -(1/) p/x + (2u/x2 + 2u/y2)
u(v/x) + v(v/y) = -(1/) p/y + (2v/x2 + 2v/y2)
Cons. of Mass
u/x + v/y = 0
Boundary Conditions:
u = v = 0 at y = 0
Ue outside boundary layer
Blasius Solution for Laminar, dP/dx=dUe/dx=0 Flow
Pressure can vary only along boundary layer not through it!
Blasius Solution for Laminar, dP/dx=dUe/dx=0 Flow
Prandtl made some simplifying assumptions, based on
the premise that the boundary layer is very thin (v and
small compared to U and L) so 2-D equations over flat
plate more tractable:
= 5 (x/U)1/2 or /x = 5(Rex)-1/2
u(x,y)/Ue =1 at y(Ue/[x])1/2 =
u(x,y)/Ue = 0.99 at y(Ue/[x])1/2 5.0
u(x,y)