Ludwing Prandtl BL JH Arakeri
Ludwing Prandtl BL JH Arakeri
--------~--------
48 RESONANCE I December 2000
GENERAL I ARTICLE
flows, such as the flow past a bluff body like a sphere, d'Alembert's
viscous effects need not be confined to such thin layers;
paradox was
viscosit.y t.hen has a more dramatic effect than what its
resolved by
low value might suggest. The key concept of boundary
Prandtl who
layers has now spread to many other fields; boundary
showed that
layers often arise in what are known as singular pertur-
viscous effects, no
bation problems.
matter how small
In this article we illustrate the boundary-layer concept the viscosity, can
by considering flows around three representative bodies, never be
namely a thin plate aligned with the flow, an airfoil and neglected.
a circular cylinder.
Viscous Stress and the No-slip Condition
It is useful when studying fluid motion to consider the
motion of a fluid particle or a small element of fluid.
Although forces such as that due to gravity are at times
important we ignore them here. For the purposes of this
article the only forces that we will consider are those due
to pressure and viscosity. These forces can accelerate or
decelerate a fluid particle.
(a) (b)
line of dye at Figure 1. Illustration of the
initial time line of dye at
flow flow short time later no-slip condition. Fluid is
flowing past a stationary
solid surface. The figure
shows how an initially ver-
tical/ine ofdye is displaced
at a later time by the flow.
The fluid and the dye next
wall wall to the wall do not move.
-RE-S-O-N-A-N-C-E--I-D-e-ce-m-b-e-r-2-0-0-0---------~~----------------------------4-9
GENERAL I ARTICLE
-so--------------------------------------------~-------------------------------------------
RESONANCE I December 2000
GENERAL I ARTICLE
.. laminar .!
'
transi tion l'
zone
. turbulent ..
-RE-S-O-N-A-N-C-E--I-D-e-c-em---be-r--2-00-0----------~------------------------------------~
GENERAL I ARTICLE
--------~--------
52 RESONANCE I December 2000
GENERAL I ARTICLE
Box 1. The Equations Governing Fluid Flow and their Approximate Forms.
Here we summarize for the mathematically minded reader the relevant two-
dimensional field equations in cartesian coordinates. In order to simplify
matters we only consider the flow over a thin flat plate aligned with the
flow. Let p and f1 be the density and viscosity of an incompressible fluid
and let ('11" v) be the (x, y) velocity components and p be the pressure. The
non-linear Navier-Stokes equations, which govern viscous fluid motion then
take the form
'/J. x + Vy = 0 (la)
+ V11. y ) =
p('/J.11. x -Px + f1(V,xx + 11. yy ) (lb)
p(v,vx + vV y ) = -Py + f1(V xx + V yy ) (Ie)
where the subscripts indicate partial differentiation. The first of the above
equations is the continuity or the conservation of mass equation. The next
two are the x and y momentum equations, respectively. For a viscous fluid
1J, and v have to vanish on stationary solid boundaries.
-RE-S-O-N-A-N-C-E--I-D-e-ce-m-b-e-r-2-0-0-0---------~-----------------------------~-
GENERAL I ARTICLE
d 2 71. du
f -2 + (1 + f)- + u = 0, 71(0) = 0, u(l) = 1. (1)
dy dy
If we wish to draw an analogy with the fluid flow equations, we may consider
f, U and y to be analogous to the viscosity, strean1wise velocity and the direction
nonnal to the plate, respectively. Note that the E in the coefficient of the middle
term is there only to silnplify the solution. The reader must be warned however
that the analogy is only a very crude, qualitative one. In any case, since the
differential equation is a linear one with constant coefficients it can be solved
exactly by, for example, assun1ing exponential solutions. The exact solution
which satisfies the above boundary conditions is given by
e-Y - e- Y/ E
u(y)= e -1 - e -1/' (2)
E
Now note that when E -r 0 (i.e. the viscosity vanishes)-l/E -r -00 and for
y > 0 -yiE -r -00; as a consequence for every positive y, u(y) -r e 1- y . This
is the outer solution. As E -r 0 this solution is a very good approximation to
the exact solution over most of the field, 0 (e) < y ~ 1. But it does not satisfy
the boundary condition at y=O since it gives the value e for u instead of O. This
pathology immediately suggests that although the outer solution is valid almost
everywhere a boundary layer is required near y = o.
We now outline the boundary layer analysis of (1) as E -r O. First, to directly get
the outer solution we let f -r 0 in (1) and get the equation
duo
- +710 = O. (3)
dy
for the outer solution uo(Y). Note that the second derivative term has been
dropped as has the small term multiplying the first derivative. Equation (3) can
be easily integrated and made to satisfy the boundary condition at x = 1. We
then get the outer solution: 71 o (y) = e 1 - y , just as we had found from the exact
solution.
Box 2. continued...
-54-----------------------------~------------------------------
RESONANCE I December 2000
GENERAL I ARTICLE
The trickier part, the Prandtl part, is to get the boundary layer or inner solution.
We know from our discussion above that although the outer solution goes toe as
y -+ 0, the correct solution must actually go to zero very rapidly in a thin layer
of thickness Eo This suggests scaling or magnifying the region close to y = O. Let
'" = YIE be the inner variable; note that for fixed y , '" -+ 00 as € -+ 0 and so '"
tends to magnify y. Now we assume that the inner solution Ui(",) isa fUnction of
the inner variable", alone. Note that now -Jy =€-l!:ry and ~ ::::: E-2 £,.
Now if
Ui("') is assumed to satisfy (1) with", defined as above, and if the substitutions
are made and the leading terms in € alone are collected we find
d 2 u'
+ _duo
1 . = O.
(4)
-;;;;r
1,
d",
1.00 l
Figure A. Comparison of the exact y ~
solution u(y) (--) with the inner and
outer solutions u, (y) (---) and u/y)
(- - - -) for two valu(!s of & • Note that
as & becomes smaller the perturba-
j
tion solution becomes better. Com-
pare with Figure 5. 0.00 1
0.0
i "
Box 2. continued...
-R-ES-O-N-A-N-C-E--I-o-e-c-em--be-r--2-00-0---------~~-------------------g-
GENERAL I ARTICLE
(6a)
-56-----------------------------~------------------------------
RESONANCE I December 2000
GENERAL I ARTICLE
Tw ~ 3/2jE!lP
0.332Uoo - ex /-LU oo / 8.
x
-RE-S-O-N-A-N-C-E-\--D-eC-e-m-b-er--2-00-0---------~------------------------------
GENERAL I ARTICLE
b)
--------~--------
58 RESONANCE I December 2000
GENERAL I ARTICLE
t~\,·--···.L.,,~~.~.~,:~~>:~~l
112 P Uoc2
'>"''"1
/'. . +
....."...,/".r
/'
,...,.;:i4- upper surface
streanlline picture (Fig'uTe 6) and t.he surface pressure Figure 7. The pressure dis-
tributions on the upper and
(Figure 7) are very dose to those obtained in the ideal
lower surfaces of an airfoil.
flow case. That is why the lift. force, which is lllainly
There is negligible differ-
caused by pressure difference between thc top and bot- ence between the pressure
tOln surfaces. predicted using ideal flow theory is close distribution calculated as-
to the rneasured force. suming non-viscous flow
and the measured pressure
On the ot.her hand, the drag is caused by viscous stresses distribution in flow ofa low-
at the wall and canllot be predicted just based on an viscosity fluid. The higher
ideal flow calculation. But a calculation based on bound- pressure on the lower sur-
ary layer t.heory is effective. face and the lower pres-
sure on the upper surface
Flow around Bluff Bodies: Boundary Layer Sep- result in lift.
aration.
--------~--------
RESONANCE I December 2000 59
GENERAL I ARTICLE
a) (b)
-J-··--r·-.. .,---.--"I·-...........--·.,E
ABC D
Figure 9. Streamline picture of the flow of a low-viscosity fluid past a circular cylinder. Except for
the boundary layers the flow in the front portion of the cylinder is nearly identical to the ideal flow
shown in Figure 8. Just after the maximum thickness point the increasing pressure with distance
causes the boundary-layer fluid to reverse direction and the boundary layer separates. In the rear
portion of the cylinder the flow, called the wake, is completely different from ideal flow; it is
unsteady and turbulent.
--------~--------
60 RESONANCE I December 2000
GENERAL I ARTICLE
-RE-S-O-N-A-N-C-E--I-D-e-c-e-m-b-e-r-2-0-0-0----------~-------------------------------6-1
GENERAL I ARTICLE
-62------------------------------~-------------------------------
RESONANCE I December 2000
GENERAL \ ARTICLE
are approximately 1000 kg/m 3 and 1 x 10- 3 kg/ms, re- Suggested Reading
spectively; the corresponding values for air are approx-
imately 1 kg/m 3 and 0.15 x 10- 4 kg/ms. This implies [1] S Goldstein (ed), Mod-
that if U is in km/hr and l is in m, Reair ~ 1.85 x 104Ul ern Development in
and Rewater ~ 2.8 x 10 5 Ul. We immediately see that un- Fluid Dynamics, Vol.l,
Dover Publ., 1965.
der most normal circumstances encountered in everyday
[2] D J Tritton, Physical
life the Reynolds numbers are likely to be very large; Fluid Dynamics, Van
thus boundary layers are likely to be present and the Nostrand, 1977
flows are likely to be turbulent.
The boundary layer on a thin lIn x 1m flat plate trav-
elling at 60 knlph (~17m/s) in air will remain larni-
nar nearly till the end of the plate, at which point the
boundary layer thickness is about 4mm. The drag on the
plate ~ 0.002 x 2 x (1/2pU 2 ) X 1 = 0.55N ~ 0.057kgf.
The same plate placed head-on will have a drag force
~ 1/2(pU 2 A) ~ 140N~14kgf. A person with a frontal
area ~ 0.317},2 travelling at 60kmph on say a motorcycle
will feel a drag force of about 4kgf.
It is close to a hundred years since Ludwig Prandtl in-
troduced the idea of a boundary layer. By doing so he
resolved a paradox that had been puzzling scientists for
almost two centuries. It took another 50 years to un-
derstand mathematically what Prandtl had intuitively,
through his own genius, seen to be true. Now his ideas
are not only routinely utilized in many flow situations
but are being applied to many other branches of the
physical and engineering sciences.
Acknow ledgements
Address for Correspondence
We thank M urali Cholemari and Shyama Prasad Das Jaywant H Arakeri
of Mechanical Engineering Departnlent, IISc, and V S Department of Mechanical
Engineering
Narasimha of NAL, for their help in the manuscript Indian Institute of Science
preparation. Bangalore 560 012, India.
PN Shankar
CTFD Division, NAL
Bangalore 560017, India.
-RE-S-O-N-A-N-C-E--I-D-ec-e-m-b-e-r-2-0-0-0---------~-----------------------------6-3