0% found this document useful (0 votes)
47 views5 pages

Answers1 (1) Fluid Mechanics

The document summarizes key concepts related to laminar boundary layers and the Navier-Stokes equations. 1. It provides the equations for 2D, incompressible flow and solves for the Couette-Poiseuille flow velocity profile between two parallel plates, showing it is the sum of shear and pressure-driven parabolic components. 2. For an oscillating plane, it solves the boundary value problem to find an exponentially damped transverse wave velocity profile propagating into the fluid. 3. It briefly discusses the conservation of energy equation and angular momentum of fluid elements.

Uploaded by

joe val
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
47 views5 pages

Answers1 (1) Fluid Mechanics

The document summarizes key concepts related to laminar boundary layers and the Navier-Stokes equations. 1. It provides the equations for 2D, incompressible flow and solves for the Couette-Poiseuille flow velocity profile between two parallel plates, showing it is the sum of shear and pressure-driven parabolic components. 2. For an oscillating plane, it solves the boundary value problem to find an exponentially damped transverse wave velocity profile propagating into the fluid. 3. It briefly discusses the conservation of energy equation and angular momentum of fluid elements.

Uploaded by

joe val
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 5

Laminar Boundary Layers

Answers to problem sheet 1: Navier-Stokes equations

The Navier–Stokes equations for 2d, incompressible flow are


∂u ∂v
+ = 0 (1)
∂x ∂y
∂2u ∂2u
! !
∂u ∂u ∂u ∂p
ρ +u +v = − +µ + (2)
∂t ∂x ∂y ∂x ∂x2 ∂y 2
∂2v ∂2v
! !
∂v ∂v ∂v ∂p
ρ +u +v = − − ρg + µ + (3)
∂t ∂x ∂y ∂y ∂x2 ∂y 2

1. Couette-Poiseuille flow
(a) For steady flow
∂u ∂v
= = 0,
∂t ∂t
i.e., we can ignore the time derivatives in equations 2 and 3.
The pressure gradient is applied along x and the upper plate is moving in the x
direction, so the vertical velocity v = 0. The continuity equation 1 therefore yields,
∂u
= 0,
∂x
which tells us that the horizontal velocity u = u(y) only. We can therefore ignore all
x−derivatives of u in Eqns. 2 and 3, which are now reduced to

∂p ∂2u
− +µ 2 = 0 (4)
∂x ∂y
∂p
− − ρg = 0. (5)
∂y
These must be solved subject to the following boundary conditions:

u = U, v = 0 for y = a

u = v = 0 for y = 0.
Integrating equation 5 gives,
p = −ρgy + A(x),
where A(x) is an integration constant, i.e. a function of x only. We identify −ρgy as
the hydrostatic pressure and note that it does not interact with the flow: it does not
appear in the equation for u since its x−derivative is zero.

To find the velocity profile we integrate equation 4, keeping in mind that the applied
pressure gradient ∂p/∂x is constant:

1 ∂p 2
u= y + By + C.
2µ ∂x
To determine the constants B and C we apply the boundary conditions:

u(y = 0) = 0 ⇒ C = 0,

U 1 ∂p
u(y = a) = U ⇒ B = − a.
a 2µ ∂x
Hence, the solution is
1 ∂p U
u(y) = y(y − a) + y.
2µ ∂x a

(b) The velocity profile is a superposition of two components, shown as dotted lines in
the figures. The straight line Uy/a is the shear flow resulting from the movement of
the upper plate y = a at speed U. The parabola is due to the pressure gradient. It
either points forward (favourable gradient ∂p/∂x < 0) or backwards (adverse gradient
∂p/∂x > 0). In the former case the maximum velocity occurs above the centre line a/2;
in the latter case below. This can be seen from
∂u a U µ
=0⇒y= −
∂y 2 a ∂p/∂x

If ∂p/∂x < 0

y
x

If ∂p/∂x > 0

y
x
(c) Total rate of throughout:
Z a
Q = u(y)dy
0
#a #a
1 ∂p y 3 ay 2 U y2
" "
= − +
2µ ∂x 3 2 0
a 2 0
3
1 ∂p a Ua
= − +
2µ ∂x 6 2
Mean velocity:
Q 1 ∂p a2 U
< u >= =− +
a 2µ ∂x 6 2

Per unit area, the drag force on the plane y = 0 is


! !
∂u ∂p a µU µU 1 ∂p
(Πxy )y=0 = µ = (y − ) + = − a.
∂y y=0
∂x 2 a y=0
a 2 ∂x

2. Oscillatory plane
(a) The oscillating plane is infinite, so there is a translational invariance in the x di-
rection: the velocity field does not depend on x. There is no applied pressure gradient
∂p/∂x. So we can set all x−derivatives equal to zero in Eqns. 1 to 3.
From the continuity equation, we then get
∂v
= 0,
∂y
which tells us that the vertical velocity v = 0 everywhere, since v = 0 on the boundary.
We therefore only need worry about the horizontal velocity field u = u(y, t). To
calculate this, we use the x-momentum equation

∂u ∂2u
= ν 2, (6)
∂t ∂y
where ν = µ/ρ is the kinematic viscosity. (You might identify this as a diffusion equa-
tion.) This must be solved subject to the boundary condition

u(t, y = 0) = U cos ωt

We seek a solution of the form

u(t, y) = f (y) cos ωt = Re(f (y)eiωt ),

where f (y) is a complex function. Substituting this into equation 6, we get

∂2f
iωf (y) = ν ,
∂y 2
for which the general solution is

f (y) = Ae−(1+i)ky + Be(1+i)ky ,


q
with k = ω/2ν.

Hence,
u(y, t) = Re{Ae−ky ei(ωt−ky) + Beky ei(ωt+ky) }.
Physically, the velocity must remain finite for all positive y no matter how far we are
from the oscillating plane, so B = 0. Also, on the solid boundary, u(t, 0) = U cos ωt, so
A = U, and we finally get

u(t, y) = Ue−ky cos(ωt − ky).

This indicates that the velocity oscillation imposed at the surface of the plane propagates
towards the interior of the viscous fluid as an exponentially damped, transverse wave
with velocity ω/k.

(b) The penetration depth δ of the velocity field is the distance from the plane for which
the amplitude of the disturbance is decreased by a factor e−1 . Thus
s
1 2ν
δ= = .
k ω
For water, ν = 10−6 m2 /s and with ω = 106 Hz, we get δ = 1.4 µm.

(c) The drag force on the plate per unit area is


!
∂u ω 3π
r
(Πxy )y=0 =µ = µUk(sin ωt − cos ωt) = µU cos(wt − ).
∂y y=0
2ν 4

3. ! !
∂ ∂ ∂ρ ∂ρui ∂E ∂E DE
(ρE) + (ρui E) = E + +ρ + ui =ρ ,
∂t ∂xi ∂t ∂xi ∂t ∂xi Dt
because
∂ρ ∂ρui
+ =0
∂t ∂xi
by the continuity equation.

4. Consider a two-dimensional section of a cubic element of fluid (i.e. a rectangle of fluid


of sides δx, δy). Taking moments about the centre 0 of the rectangle, gravity forces and
normal stresses (Πxx , Πyy ) contribute nothing since their line of action passes through
0. The net moment about 0 is due to shear stresses only, and in a clockwise sense this
moment is,
∂Πxy δy ∂Πyx δx
(Πxy + δy + Πxy )δx − (Πyx + δx + Πyx )δy
∂y 2 ∂x 2
( )!
1 ∂Πxy ∂Πyx
= Πxy − Πyx + δy − δx δxδy.
2 ∂y ∂x

The moment of inertia of the cube about 0 is, per unit depth in the paper,
1
ρδxδy[(δx)2 + (δy)2].
12
y 6 ∂Πyy 6
Πyy + δy
∂y Πxy
-
Πyx ∂Πxx
?
O Πxx + δx
 δy . - ∂x

Πxx δx 6 ∂Πyx
Πyx + δx
 ∂x
∂Πxy
Πxy + δy Πyy
∂y ?
-
x

Denoting the angular acceleration about 0 by a, Newton’s second law for a system in
rotation gives
( )!
1 1 ∂Πxy ∂Πyx
ρδxδy[(δx)2 + (δy)2 ]a = Πxy − Πyx + δy − δx δxδy.
12 2 ∂y ∂x

Taking δx and δy to both have the same order of magnitude δ, we notice that the LHS
scales as aδ 4 and the RHS as (Πxy − Πyx )δ 2 + O(δ 3 ). To avoid a becoming infinite in the
limit as δ → 0, the prefactor to the δ 2 term on the RHS must be zero, so we have

Πxy = Πyx .

This reasoning can be extended to a three-dimensional elementary cube of fluid to show


that the stress tensor [Π] is symmetric.

You might also like