BSL 2B3 - Sarec/Fiq

BSL 2B.3 Laminar flow in a narrow slit
(a) A Newtonian fluid is in laminar flow in a narrow slit formed by two parallel walls a
distance 2B apart. It is understood that B

€
€
€
€
€
€
€
€
p(z). The only nonvanishing components of the stress tensor are τxz = τzx, which depend only
on x.
Consider now a thin rectangular slab (shell) perpendicular to the x−direction extending a
distance W in the y−direction and a distance L in the z−direction. A 'rate of z−momentum'
balance over this thin shell of thickness Δx in the fluid is of the form:
(LW) ( φxz − φ
x
xz
x +Δx ) + (WΔx) ( φzz − φ
z=0
zz
z= L)
+ (ΔxLW)ρg = 0
Dividing the above equation by LWΔx, and the limit taken as Δx approaches zero, we get
∂φxz ∂x − φzz − φ
z=0
zz
z= L − ρg = 0
L
At this point we have to write out explicitly what components φxz and φzz are, making use of
the definition of φ in Eqs. 1.7-1 to 3 and the expressions for τxz and τzz in Table B.1. This
ensures that we do not miss out any of the forms of momentum transport. Hence, we get
φ xz = τ xz + ρv xv z = τ xz
φ zz = p + τ zz + ρv zv z = p
In accordance with the postulates that vz = vz(x), vx = 0, vy = 0, and p = p(z), we see that (i)
since vx = 0, the term ρvxvz is zero; (ii) since vz = vz(x), the term τzz is zero; (iii) since vz =
vz(x), the term ρvzvz is the same at both ends of the slit, so the convective terms are cancel
out.
dτxz dx − p0 − pL − ρg = 0
L
dτxz dx = p0 − pL + ρg
L
The right – hand side may be compactly and conveniently written by introducing the modified
pressure P, which is the sum of the pressure and gravitational terms. The general definition of
the modified pressure is P = p + ρgh, where h is the distance upward (in the direction opposed
to gravity) from a reference plane of choice. Since the z−axis points downward in this
problem, h = − z and therefore P = p − ρgz. Thus, P0 = p0 at z = 0 and PL = pL − ρgL at z = L
giving p0 − pL + ρgL = P0 − PL.
dx = P ⎛⎛
0 − P ⎞⎞
L
⎜⎜ ⎟⎟
⎝⎝ L ⎠⎠
dτ xz
Integration leads to the following expression for the shear stress distribution
τxz = P0 − P ⎛⎛ ⎞⎞
L
⎜⎜ ⎟⎟ x + C1 ⎝⎝ L ⎠⎠
€
€
Copyright © 2005 by Apolinar Picado

€
€
€
€
€
€
It is worth noting that the two last equations apply to Newtonian and non – Newtonian fluid.
The constant of integration C1 is determined later using boundary conditions.
Substituting Newton's law of viscosity for τxz in above equation gives
−µ dv z
dx = P0 − P ⎛⎛ ⎞⎞
L
⎜⎜ ⎟⎟ x + C1 ⎝⎝ L ⎠⎠
The above first – order differential equation is simply integrated to obtain the following
velocity profile
v z = − P0 − P ⎛⎛ ⎞⎞
L
⎜⎜ ⎟⎟ x
⎝⎝ 2µL ⎠⎠
2 − C1 µ x + C2 Boundary conditions
BC1: at x = B; vz = 0
BC2: at x = − B; vz = 0
Using these, the integration constants may be evaluated as C1 = 0 and C2 = (P0 – PL)B 2 /(2µL).
Substituting C1 = 0 in the shear stress equation, the final expression is found to be linear as
given by
τxz = P0 − P ⎛⎛ ⎞⎞
L
⎜⎜ ⎟⎟ x
⎝⎝ L ⎠⎠
Further, substitution of the integration constants gives the final expression for the velocity
profile as
v z = (P 0 − P L )B 2
2µL
1− x
⎡⎡ 2
⎛⎛ ⎞⎞ ⎤⎤
⎢⎢ ⎜⎜ ⎟⎟ ⎥⎥
⎣⎣ ⎢⎢ ⎝⎝ B⎠⎠
⎦⎦ ⎥⎥
It is observed that the velocity distribution for laminar, incompressible flow of a Newtonian
fluid in a plane narrow slit is parabolic.
b)
The maximum velocity occurs at x = 0 (where dvz/dx = 0). Therefore,
v z,max = (P 0 − P L )B2
2µL
v z =v z,max 1− x ⎛⎛
⎜⎜
⎝⎝ B
⎞⎞ ⎡⎡ 2⎤⎤
⎢⎢ ⎟⎟ ⎥⎥
⎣⎣ ⎢⎢ ⎠⎠ ⎦⎦ ⎥⎥
Copyright © 2005 by Apolinar Picado

€
€
€
€
€
€
€
The average velocity is obtained by dividing the volumetric flow rate by the cross – sectional
area as shown below.
v z =
∫
B
∫
∫
W
−B 0
B
−B
v z = v z,max
B
v z = v z,max
B
v z
v z,max
= 2
3
∫
v z(x)dydx
∫
0
W
⎛⎛
dydx
1− x2
B 2
B
0 ⎜⎜
⎝⎝
2
3 B
⎛⎛ ⎞⎞
⎜⎜ ⎟⎟
⎝⎝ ⎠⎠
⎞⎞
⎟⎟ dx
⎠⎠
Thus, the ratio of the average velocity to the maximum velocity for Newtonian fluid flow in a
narrow slit is 2/3.
This seems reasonable since
€
v z /vz,max = ½ for flow in a circular pipe. For a slit a larger
cross − sectional area carries fluid flowing at the larger velocity than for flow in a circular
pipe. So v z /vz,max for slit > v z /vz,max for a circular pipe.
c)
€ The mass rate of flow € is the product of the density ρ, the cross − sectional area (2BW) and the
average velocity v z .
w = ρ (2BW) v z
€
w = ρ (2BW) 2 (P0 − PL )B
3
2
2µL
w = 2 (P0 − PL )B
3
3 Wρ
µL
The flow rate vs. pressure drop (w vs. ΔP) expression above is the slit analog of the Hagen –
Poiseuille equation (originally for circular tubes). It is a result worth noting because it
provides the starting point for creeping flow in many systems (e.g. radial flow between two
parallel circular disks; and flow between two stationary concentric spheres).
d)
The above analysis is not applicable if B = W, because of the presence of a wall at y = 0 and y
= B would cause vz to vary significantly in y in addition to x, then vz = vz (x, y).
Copyright © 2005 by Apolinar Picado

€
€
€
€
If W = 2B, then a solution can be obtained for flow in a square duct.
e)
In Eq. 2.5-20, set both viscosities equal to µ, p0 − pL = P0 – PL, and set b equal to B.
v z = (P 0 − P L )B2
12µL
⎛⎛ 8⎞⎞
⎜⎜ ⎟⎟ =
⎝⎝ 2⎠⎠
(P0 − PL )B2
3µL
v z = 2
2 ⋅ (P0 − PL )B2
=
3µL
2 (P0 − PL )B
3
2
2µL
v z = 2
3 v z,max
v z
v z,max
= 2
3
Copyright © 2005 by Apolinar Picado