Ghosh - 550 Page 1 3/19/2013

Solutions of Boundary Layer Equations

Now we develop the solution strategies for the boundary layer equations given by
Prandtl. Remember in his 2 equations of continuity and x-momentum, the pressure
gradient term is assumed to be known.
Continuity: 0
y
v
x
u

(A)
x-momentum:
2
2
1
y
u
v
dx
dp
y
u
v
x
u
u

(B)
Thus this set becomes mathematically solvable. There are two approaches to solve
boundary layer equations. We shall present both here. However the emphasis will be in
the second approach since it is easier to work with and gives an insight to the behavior of
fluid particles in the boundary layer. The standard approaches are:
(i) Exact solution method (Blasius Solution)
(ii) Approximate Solution Method (Karman-Pohlhausen Method)
The second approach is also called the momentum integral method. We begin with the
exact solution method given by Blasius.
Exact Solution Method
Blasius performed a transformation technique to change the set of two partial
differential equations (A and B) into a single ordinary differential equation. He solved
the boundary layer over a flat plate in external flows. If we assume the plate is oriented
along the x-axis, we may neglect the pressure gradient term, i.e., 0
x
p

. The traditional
approach before Blasius was to drop out the continuity equation from the set by the
introduction of the stream function (x,y). With this definition:
(B)
3
3
2
2 2
y y x y x y

,
_

,
_

However, Blasius used this equation in non-dimensional variables. Let us define (x,y)
as a single variable by:
U
x
y

or
x
U
y

. It is easy to verify that will be non-

dimensional by substituting the units of , U, and y. He also introduced the non-
dimensional stream function given by: (Correction: Please replace by x in the section
below)
Ghosh - 550 Page 2 3/19/2013
U
) ( f

Using mathematical manipulation from calculus, we may write:
( ) ( )

f U
U
f U
y y
u
( )
U
u
f = non-dimensional velocity function
( )
1
]
1

x
f
U f
U
2
1
x
v (by chain rule)
( )
1
]
1

x
f
U f
U
2
1
But since
x 2
x
2
1 U
y
x
2
3

,
_

, we may simplify v into [ ] f - f

U
2
1
v .
Similarly we can show:
f
U
U
y
u
, f
U
2
1
x
u

and f
U
y
u
2
2
2

,
_

Therefore, the original x-momentum can be written as

0 f f f + 2
upon
simplifications. Note that this is an ordinary differential equation with as the
independent variable and f is the dependent variable. To solve this third order equation
we need three boundary conditions. Let us check the figure below.
(x)
At y = , u = U
At x = 0, (x) = 0 At y = 0 (on wall),
u = v = 0 (no slip condition)
U
y
x
Ghosh - 550 Page 3 3/19/2013
We may write at y = 0, u = v = 0 into a different form:
0 f f , 0 At
. Similarly
as y , u = U may be written as
1 f , At
Using these three boundary conditions the solution of the governing equation may be
obtained by the use of power series solution and shown in the table below. The important
things to note are the points corresponding to the edge of the boundary layer. Since u
U, f 1, (we choose the value of .9915, since is defined at u = .99U). Thus
0 . 5
U
y

from the table below:

Using the alternate definition,

y
, we get
x
Re
0 . 5
, or
U
0 . 5

Now, ( ) u v
x
v
where
x
v
y
u
yx
<<

,
_

0
Therefore, the wall shear stress,
w
, may be written as
x
2
0 y
w
Re
U 332 .
y
u

We define Skin Friction Coefficient as the non-dimensional wall shear stress, given by:
Ghosh - 550 Page 4 3/19/2013
x
2
w
f
Re
664 . 0
U
2
1
C

In this case both

x
Re
0 . 5
) x (


and
x
Re
664 . 0
C
f are claimed to be exact solution of
steady, laminar boundary layer over a flat plate oriented along the x-axis.
We notice from the above expressions that both (x) and C
f
change along the plate.
While (x) increases (boundary layer grows) with , C
f
0 as x . Both quantities
depend on the variable Reynolds number, Re
x

1
]
1

U
. If the plate length is not infinite,
how do we obtain the shear force on it? We may do this by integrating directly or, the
use of the concept of overall Skin Friction Coefficient. For example, for a finite length,
L, of the plate, the shear force

A
0 y
yx yx
dA F
where, wdx dA .
but:

,
_

2
f w
0 y
yx
U
2
1
C
Therefore the
x
Re
664 . 0
) x ( C
f may be substituted above and F
yx
obtained by integration.
Alternately, define
f
C
= Overall Skin Effect Coefficient

L
0
dx ) x ( C
L
1
f
. Thus the
f
C

is nothing but length-averaged friction coefficient. Unlike C
f
(x),
f
C
is a constant
value for the whole plate. Similarly the average shear stress for the plate may be defined
as

L
0
w w
dx ) x (
L
1
. Finally, the shear force on the plate may be written as the
product of
w

and the plate area.

L
w
dA = w dx
x
y
Ghosh - 550 Page 5 3/19/2013
Approximate Solution Method
Unlike the Blasius solution, which is exact, approximate solution method assumes an
approximate shape of the velocity profile. This velocity profile is then utilized to
evaluate quantities related to the governing differential equation, given below by Karman
and Pohlhausen. This method, which is called the momentum integral method, changes
the two equations given by Prandtl into a single differential equation. This equation over
a flat plate may be written as:
dx
d
U
2
w


where,
x
=

Wall shear stress
= Density of fluid
U = Free stream velocity
and, = Momentum thickness of the boundary layer
= dy
U
u
1
U
u
0

,
_

The above equation is applicable only when the pressure gradient term is zero. For
the case of non-zero pressure gradients you should use
( )
dx
dU
U U
dx
d
* 2 w
+

Velocity Profiles
Since the Karman-Pohlhausen method requires an assumed velocity profile, let us
explore some velocity profiles and their characteristics (see example problem 1). For
example, suppose we assume the velocity profile to be a second order polynomial
2
Cy By A ) y ( u + +
where A, B, and C are constants.
To evaluate velocity profile constants A, B, and C, we must use boundary conditions. The
following three conditions may be used:
1) No-slip: y = 0, u = 0
2) B.L. Edge Velocity: y = , u = U
3) B.L. Edge Shear: y = ,
dy
du
= 0
Ghosh - 550 Page 6 3/19/2013
Note that at the edge of the defined edge of the boundary layer u = .99U and
0
dy
du

.
However we approximate them with the rounded values. This is the reason the solution
method by Momentum Integral Method is considered an approximate one.
With the above profile,
1) 0 = 0 A ) 0 ( C ) 0 ( B A
2
+ +
2) U =
+

+ + C B
U
) ( C ) ( B A
2
[ 0 A ]
3)
0 C 2 B
dy
du
y
+

Subtracting the second condition from the third,
2
U
C
U
C

Using this in the second condition,

U 2
B
U
B
U
2
2
y
U
y
U 2
) y ( u

or,
2
y y
2 ) y (
U
u

,
_

,
_

Parabolic Profile (see plot in the

example)
Remember the use of the boundary layer velocity profile is only meaningful when
y 0
. The use of this velocity profile may now be made to obtain and
w

,
_

,
_

,
_

y
d
U
u
1
U
u
dy
U
u
1
U
u
1
0 0
or,

1
1
]
1

,
_

1
1
]
1

,
_

y
d
y y
2 1
y y
2
2
1
0
2
Note that defining a new variable

y
makes the evaluation much easier.
( )( ) +

d 2 1 2
2
1
0
2
Ghosh - 550 Page 7 3/19/2013
( )

+ + d 2 2 4 2
1
0
4 3 3 2 2
( )

+ d 4 5 2
1
0
4 3 2
15
2
15
3 15 25 15
5
1
1
3
5
1
+
+
Similiarly,
0 y
w
dy
du

[ 0 v in boundary layer]
0 0 y
) U u ( U
) y (
) U u ( U

[ ]

U 2
2 2
U
0
for the parabolic profile
Using the above results for and
w
in the momentum integral equation for a flat plate
gives
dx
d
15
2
U
U 2
2

,
_

Separating the variables and x, and integrating

dx
U
15 d
x
0 x
) x (
0


,
_

,
_

U
15
2
2
U
x 30
2

To express the boxed equation in a non-dimensional form divide both sides by x
2
,
x
2
Re
30
Ux
30
x

,
_

, where Re
x
=

Ux
is the Reynolds number based upon the
variable x.
x x
Re
48 . 5
Re
30
x

Compare this result with the earlier exact solution obtained under the Blasius method.
A C
Ghosh - 550 Page 8 3/19/2013
x
Re
0 . 5
x

We therefore see the popularity of the parabolic velocity profile. Although the solution by
Karman-Pohlhausen method is approximate it gives less than 10% error when compared
with the exact solution is laminar flows over a flat plate.
Now that we have obtained (x), the shear stress,
w,
and skin friction coefficient, C
f
, may
be obtained for the parabolic profile.
x
x
w
Re
x
U
365 .
x 48 . 5
Re U 2
U 2

x
x
2
2
w
f
Re
73 .
Re
x
U
U
73 .
U
2
1
) x ( C

This is comparable with
x
f
Re
664 .
) x ( C
obtained earlier in the exact solution method.
To summarize, we have obtained the growth of the boundary layer (x) and the skin
friction characteristic C
f
(x) as a solution of the boundary layer equations by the exact and
approximate methods. Once C
f
(x) is known, the shear stress and skin friction force may
be evaluated (see examples).
As stated before, the frictional forces are not the dominant forces in high-speed flows.
The component of drag due to skin friction is called the friction drag. Thus friction drag
is significantly lower than pressure drag in boundary layers of high Reynolds number
flows. However, Prandtl found a very important influence of these small frictional
forces in controlling the pressure drag. To understand this we must investigate the
phenomenon of flow separation.
Flow Separation and Boundary Layer Control
Earlier we noted that as the boundary layer over a flat plate grows, the value of the skin
friction coefficient goes down. This may be explained from the fact that as more fluid
layers are decelerated due to shear at the plate shear values near the plate need to be as
large compared to the entrance region of the plate.
Point of
Stagnation S
A C
Ghosh - 550 Page 9 3/19/2013
Compare the station (2) ( =
2
) with station (1) ( =
1
). The shear on the plate at (2) is
smaller since the shear angle
) 1 ( ) 2 (
y
u
y
u

<

. Mathematically, we know C
f
(x) 0 as Re
x
. But can the shear go to zero on the flat plate, and if so, what are the physical
implications? The answer depends on the physical configurations. For a flat plate, shear
may never go to zero as C
f
0 only when Re
x
or x . However if we get some
assistance from the pressure gradient, C
f
can be zero much earlier. Consider, for this
purpose, flow over a circular cylinder.
(1) (2)
U
x
y
A
B
C
D
Point of
Stagnation S
) 2 (
y
u

) 1 (
y
u

U
y
x

2
A C
(x)
B
Ghosh - 550 Page 10 3/19/2013
The figure above shows a circular cylinder in steady, ideal flow, U. The stagnation points
are A and C, while the maximum velocity points are B and D. Since the regions A to B
and A to D accelerate the flow, 0
dx
dp
< (note x is in the tangential direction along the
cylinder). Similiarly the regions B to C and D to C are the adverse pressure gradient
regions ( 0
dx
dp
< ). Now imagine if this cylinder was placed in a real flow, viscous
boundary layer will start to grow from the front stagnation point A, slowing the fluid
particles.
However, fluid pressure field still naturally pushes the particles near the surface to
proceed toward B. This is not the case between B and C though, where the natural
tendency of the fluid is to flow C to B due to the adverse pressure gradient. Thus the
boundary layer slow down that started in the region A to B due to viscous effects
bringing C
f
toward 0, gets compounded by the reverse push due to the adverse
pressure gradient in the region B to C. This brings the flow to separation. Flow
separation point is defined as the point on the surface where C
f
= 0, or,
w
= 0, or
0
y
u
0 y

.
At flow separation, fluid particles rest on the solid surface but there is no hold on them
due to shear from the surface. There is however shearing action from the high-speed
flow a little away from the surface, which drags these stagnant particles away into
the main flow stream due to viscosity. This creates a partial void inside the boundary
layer, which is promptly filled by particles traveling upstream creating a reverse flow
near the surface.
The figure shows real flow separation over a circular cylinder with the separation point
and reverse flow after separation. Due to symmetry, the exact same processes are
repeated on the lower surface ADC. The reverse flow near the surface is the cause of
x
y
y
x
A
B
D
C
Point of
Stagnation S
Vortex formation
C
p
= 1- 4 sin
2
(Ideal flows)
Ghosh - 550 Page 11 3/19/2013
vortex formation. Two symmetric vortices appear first in the downstream of the cylinder
following flow separation.
Real Flow Over the Cylinder
These vortices occupy the wake region since they are shed behind the cylinder due to the
forward fluid motion. As that process happens the shed vortices grow in size and start
interacting with each other creating an alternating vortex pattern known as the
Karman Vortex Street. These create oscillatory flows behind the cylinder.
Eventually all the vortices break down due to viscous interactions creating a region of
chaos, which is characteristic of a turbulent mixing.
In the initial phase a laminar separated flow is not necessarily turbulent. It creates a
large region of low pressure behind the body called the wake region. Due to the
Laminar
Wake
U
Symmetric
Vortices
Karman Vortex Street
C
p
= 1- 4 sin
2
(Ideal flows)
Ghosh - 550 Page 12 3/19/2013
separation process, the pressure never recovers its stagnation value in laminar separated
flows. If instead of a laminar follow, we had placed the cylinder in a turbulent flow,
separation will occur with a much narrower wake behind the body. This is due to the
fact that turbulent flows have flatter velocity profiles with rapid mixing and a lot more
momentum in the boundary layer. This gives turbulent flows much better chance to
resist separation in the region behind the body (B to C or, D to C). The late separation
gives a much smaller wake size with a much better pressure recovery as shown in the
figure below:
Laminar
Wake Turbulent
Wake
A
B
C
D

1
-3
C
p
= 1- 4 sin
2
(Ideal flows)
Laminar Flows
Turbulent Flows
Ghosh - 550 Page 13 3/19/2013
Thus the drag calculated in the turbulent flows will be much smaller compared to laminar
flows (Recall that ideal flow drag is zero due to 100% pressure recovery). This is the
reason a turbulent flow separation is preferred over a laminar flow separation (see
example of flow momentum calculation). The drag coefficient versus Reynolds number
for the flow over a sphere is shown below.
The figure shows that drag coefficient drops as the Reynolds number increases in the low
speed range. In this range, drag on the sphere is directly proportional to the
diameter of the sphere (F
D
= 3VD) as was shown by Stokes. On the other hand, for
high-speed flows,
D
C
b
A
2
U
2
1
D
F . Thus, if C
D
is constant,
2
U
D
F . In the low speed
range, drag on the sphere is mostly due to friction, whereas in the high-speed range drag
is mostly (due to flow separation) from the pressure drag. The sharp drop in the C
D
curve around Re
D
= 2x10
5
is due to the transition from laminar to turbulent flows.
As we saw earlier, transition into turbulence brings smaller wake size and a lower overall
drag. This feature is often incorporated into design. For example, golf balls are dimpled
Ghosh - 550 Page 14 3/19/2013
to take advantage of this fact. The dimples cause early tripping of the flow into
turbulence. This would reduce the drag and will produce longer flights of the ball.
Drag reduction is an active design topic for aerodynamicists and fluid mechanists. A
major controlling feature of laminar flow separation is by removal of stagnant fluid
particles near the walls by suction. Similarly by blowing into boundary layer, we
may be able to energize the stagnant particles and prevent separation. Control of
separation and drag reduction in various applied problems is an active area of research.
Turbulent Boundary Layers
We know that turbulent flow occurs if the flow velocity is large enough (or, viscosity is
small enough) to create a Reynolds number greater than the critical Reynolds number
over an object. For spheres or circular cylinders this critical Reynolds number is between
2 to 4x10
5
. For flat plate flows this is around 500,000. We also discussed the implications
of turbulent flows in drag reduction. What characterizes such flow is a flatter, fuller
velocity profile. It is important to recognize that turbulent flows have two components:
(i) a mean,
u
, and (ii) a random one, u . u u u + . Similarly, v v v + and
w w w + . The random u cannot be determined without statistical means. Therefore
for turbulent fluid flows, we usually work with a time averaged mean flow u . Remember
that when we speak of turbulent velocity profiles it is this u that we are considering. To
avoid confusion with this rotation (we earlier indicated v as area averaged velocity, not
time average velocity), we shall write turbulent flow velocities without the bars.
You understand that whenever we speak about turbulent flows here, we are representing
the mean flow. Turbulent flows in boundary layers over flat plates may be represented by
the power law velocity profile:
n
1
n
1
y
) y (
U
u

,
_

[where,

y
]
This profile covers a fairly broad range of turbulent Reynolds numbers for 6 < n <10. The
most popular one is n = 7. Although this velocity profile is an excellent representation of
the real turbulent flow, this may not be used to calculate skin friction coefficient in the
approximate solution method seen earlier (since
0 y
y
u

will be negligible for this

profile). For the purpose of calculating shear stress we use an experimental result:
4
1
2
w
U
U 0233 . 0
,
_

for the 1/7 power law profile. To obtain the skin friction
coefficient, we must first evaluate (x) from the solution of Karman-Pohlhausen:
Ghosh - 550 Page 15 3/19/2013
x
U
2
w

(A)
dy
U
u
1
U
u
dy
U
u
1
U
u
1
0 0

,
_

,
_

Q
Using
7
1
U
u
,

,
_

d 1
7
1
1
0
7
1
72
7
9
7
8
7
d
1
0
7
2
7
1

,
_

) A (
dx
d
72
7
U
U
U 0233 . 0
2
4
1
2
w

,
_

,
_

or,
dx
U
24 . d
x
0 x
4
1
) x (
0
4
1

,
_

5 1
x
Re
382 .
x

(Skipping the integral evaluation)

Terms for the 1/7 power law velocity profile gives:
5
1
x
f
Re
0594 .
C
Stokes Flows
We have so far discussed very high-speed flows in which the boundary layers are very
thin regions near the body. However for very low speed flows boundary layers dont
exist. Viscous effects are felt everywhere (recall the heat transfer analogy). External
flow applications at very slow speeds (or, highly viscous flows) may be solved by
neglecting the inertia force term in the Newtons second law. For example, if you drop a
steel ball into glycerin, how can you calculate drag on it? In this context, let us introduce
the concept of terminal velocity. When any object starts its motion in any fluid medium,
there may be a period of acceleration of motion. However, if we are interested in steady
flows, if one exists in such configurations, there must be a time when the fluid forces
around the body are balanced providing it a constant velocity. We call this velocity
terminal velocity of the body. For the ball dropped in glycerin, the free body diagram
shows
Ghosh - 550 Page 16 3/19/2013
Note that the V
t
(Terminal velocity) is not a force, and shown on the sketch (just for
reference) using dashed lines. Since the body is traveling at constant speed, the inertia
force term is zero. Thus, all external forces are balanced and in the Newtons second law
in the y-direction:
0 ma F F
y
y
S
y
B

(A)
0 F F W
D B

where, mg W
V g F
B

D V 3 F
t D

6
D
V
3

D = Diameter, = Density of glycerin

g = Acceleration due to gravity
We may only use the above equation to calculate the terminal velocity V
t
.
In the above drag representation of Stokes flow, F
D
V
t
. This behavior is in contrast
with high-speed flows, where drag is usually proportional to the square of velocity.
Terminal velocity concept is similar to fully developed flows in internal flow
configuration. Notice that before the terminal velocity is developed (in the internal flow
case, in the entrance length region), the inertia force term in the above equation (A) is not
negligible. In that case, the only way to solve the equation will be by integration or, using
differential equations approach (see example 8).
y
w
g
v
V
t
F
B
F
D
F
D
= Drag on the body
F
B
= Buoyancy force on the body
w = weight of the ball
Ghosh - 550 Page 17 3/19/2013
For engineering design purposes, handbooks list a large variety of objects in different
orientations and their drag coefficients. Rather than solving each problem from first
principles, you may be able to utilize these tables and charts. Just make sure that you note
the range of applicability of these. They need to be verified during problem solving.
continue