Worksheet Flow
Worksheet Flow
Page 1
3/10/2016
2. System Diagram
Velocity Profile
yx
y
x
xy
xy
yx
3. Assumptions
Steady state condition
Constant fluid properties (density and viscosity)
2-D problem (xy plane)
Laminar, fully developed flow
4. Governing Equations
u v
0
x y
Ghosh - 550
Page 2
3/10/2016
DV
g p 2V
Navier-Stokes Equation:
Dt
v u
Shear Stress (for xy plane): xy yx
x y
5. Detailed Solution
Shear stress as a function of y
u = u(y) and v 0 because the fluid flows between two parallel plates.
Considering this fact, the continuity equation becomes:
0
u v
u
0
0
x y
x
This tells that the flow is fully developed, which is given in the problem and assumed to be,
because the velocity doesn't change in the x direction. Therefore, the velocity u is function of
only y, u = u(y).
Navier-Stokes equation in y direction is useless on this problem because v 0.
Navier-Stokes equation in x direction:
Continuity equation
gx = 0
2u 2u
u
u
u
p
u v g x 2 2
x
y
x
y
t
x
Steady state
v0
Continuity equation
Therefore,
p
2u
0 2
x
y
2 u 1 p
d 2 u 1 p
because u = u(y).
y 2 x
dy 2 x
du
1 p
d dy x dy
du 1 p
y Const
dy x
Ghosh - 550
Page 3
3/10/2016
This problem deals with the flow between two parallel plates, so the maximum velocity
occurs at y = 0 (centerline). This implies
du
dy
du 1 p
du p
y
dy x
dy x
yx
u
du
yx
because u = u(y).
y
dy
We already have this expression above; therefore, the shear stress is:
yx
du p
y
dy x
Ghosh - 550
Page 4
-3
2.5
3/10/2016
x 10
Upper Wall
2
1.5
1
y (ft)
0.5
yx = (p/x)y
0
-0.5
-1
-1.5
-2
Lower Wall
-2.5
-0.01 -0.008 -0.006 -0.004 -0.002
0.01
yx (lbf/ft )
2
yx,max
p h
x 2
4.167 10 3 ft
4.0lbf / ft 2 / ft
2
2
= 0.008333 lbf/ft
6. Critical Assessment
The shear stress is negative for y > 0 (see the plot). This is correct because the way shear
stress is defined. See the "System Diagram." The upper shear stress is positive toward
right, but on the upper wall, the shear stress works toward left in reality (since the shear
Ghosh - 550
Page 5
3/10/2016
force must be positive and the upper wall is a negative surface). This is why it is
negative.
[Note: Instead of making the fully developed flow assumption, if only the parallel flow
assumption is made, the analysis will follow the same pattern since parallel flows are
fully developed. Simplifications will be easier since u0, but v=w=0.]
2. System Diagram
u(y)
y
x, u
umax
3. Assumptions
Steady state condition
Incompressible fluid flow
4. Governing Equations
Ghosh - 550
Page 6
3/10/2016
Q
A
5. Detailed Solution
It is obvious to tell the units of constants (A, B, and C) by observing the given velocity profile
function: u = umax(Ay2 + By + C). Note that each term inside the parenthesis has to have no
units. Therefore the units on A, B and C are as follows.
Units
m-2
A
m-1
B
1
C
Note: This will be verified after the constants (A, B, and C) are evaluated in terms of h.
Evaluation of Constants
Realizing the fact that the maximum velocity occurs at y = 0 in this problem condition and
geometry, the constant C can be evaluated as follows:
u = umax(Ay2 + By + C)
umax = umax[A(0)2 + B(0) + C]
C=1
Because the maximum velocity occurs at y = 0, this relation must be satisfied:
du
dy
du
dy
y 0
d
u max Ay 2 By C
dy
y 0
u max 2 Ay B y 0 B u max 0
u m ax 0, B 0
Ghosh - 550
Page 7
3/10/2016
h 2
h
0 u max A B C
2
2
1 2 1
Ah Bh C
4
2
1 2
4
Ah 1 A 2 .
4
h
u u m ax 2 y 2 1
h
or
u
u max
y
1 4
h
Q 2h u b dy
2h u m ax 2 y 2 1 b dy
2
h
4 1
2
u m ax b 2 y 3 y
h 3
h
2
u m ax b h
3
Q 2
hu max
b 3
Ratio V / u m ax
Q Q Q1 2
1 2
humax u max
A bh b h 3
h 3
V
2
u m ax 3
6. Critical Assessment
The problem illustrates how to determine velocity profiles from boundary conditions
and calculate volumetric flow rate and average velocity. An important point on this
Ghosh - 550
Page 8
3/10/2016
problem is to realize where the maximum velocity occurs and how the mathematical
expression can be written for it.
2. System Diagram
a
U
r1
r2
a
u
3. Assumptions
Steady state condition
Incompressible fluid flow
Laminar flow
Ghosh - 550
Page 9
3/10/2016
4. Governing Equations
u v
0
x y
DV
g p 2V
Navier-Stokes Equation:
Dt
v u
Shear Stress (for xy plane): xy yx
x y
2-D Incompressible Continuity Equation:
5. Detailed Solution
Viscosity of the Oil
The problem gives the velocity profile to be linear, so we can simply say that the velocity
profile is (See "Critical Assessment" for a detailed discussion of this velocity profile):
u u( y)
U
y
a
yx
u
U
y
y a
r
U
1
a
a
Torque is given by T F R yx A R yx 2 r1 L r1
2 r13 L
Equating and gives: T
a
Therefore,
y .
Ghosh - 550
Page 10
Ta
2 r13 L
3/10/2016
0.2 N m 1 10 3 m
3
rev 1 min 2rad
3
2 25 10 3 m 2800
100 10 m
min 60s 1rev
= 0.0695 Ns/m2
6. Critical Assessment
The problem gives the velocity profile as linear, so we simply said that the velocity
profile is
u u( y)
U
y
a
But here is the detailed demonstration of how this linear velocity profile was obtained.
u v
u
0
0
x y
x
This shows that the flow is fully developed, which is assumed to be, because the velocity
doesn't change in the x direction. Therefore, the velocity u is function of y only, u = u(y).
Navier-Stokes equation in y direction is useless on this problem because v 0.
Navier-Stokes equation in x direction:
Continuity equation
gx = 0
2u 2u
u
u
u
p
u v g x 2 2
x
y
x
y
t
x
Steady state
v0
Continuity equation
Ghosh - 550
Page 11
3/10/2016
Therefore,
p
2u
2
x
y
2 u 1 p
d 2 u 1 p
because u = u(y).
y 2 x
dy 2 x
Integrating this equation twice gives: u
Boundary Conditions
u = 0 at y = 0 C2 = 0
u = U at y = a C1
1 p 2
y C1 y C 2
2 x
U
1 p
a
a 2 x
Finally,
2
U
a 2 p y y
u u( y) y
a
2 x a a
u u( y)
U
p
y because
0 (See assumptions listed above).
a
x
dy
y
Bath
Belt
Ghosh - 550
Page 12
3/10/2016
2. System Diagram
U0
xy
h
p = patm
yx
yx
x
y
xy
,
3. Assumptions
Steady state condition
Constant property fluid
v0
2-D problem (xy plane). That is, the z-depth is infinite u
0
z
Ghosh - 550
Page 13
4. Governing Equations
Navier-Stokes Equation
3/10/2016
DV
g p 2V
Dt
xy yx
x y
5. Detailed Solution
The problem gives the flow is fully developed. Thus,
u
u(y) only.
0 . Also, . u
0
x
z
Steady state
Fully developed
2u 2u
u
u
u
p
u v g x 2 2
x
y
x
y
t
x
Fully developed
Given
Therefore,
0 g
x
Shear Stress
We are interested in yx and v 0 in this problem. Thus,
xy yx yx
y
x y
Boundary Conditions
@ y = 0, u = U0
@ y = h, yx
u
u
0
0 because the atmosphere produces no shear stress at
y
y
Ghosh - 550
Page 14
3/10/2016
0 g
x
du
du
0 g
dy
dy
2
d u g
dy
y C
dy
1
g
h.
y
h
dy
g 2 gh
y
y C2 .
2
u( y)
1 g 2 gh
y
y U0
2
or
g 2 1 y y
h U 0
2 h h
2
u( y)
6. Critical Assessment
Here, the velocity profile has been obtained by a direct application of Navier-Stokes
equation. The same result can be obtained by applying the momentum equation for
integral control volume:
FS FB VdV VV dA
CS
t CV
Ghosh - 550
Page 15
3/10/2016
128L
D 4
which is independent of flow rate. Find the maximum pressure drop for which this relation
is valid for a tube 50 mm long with 0.25 mm inside diameter for both kerosine and caster oil
at 40 C.
Ohm's Law
Fluid Flow
Driving Force
V (voltage)
p (pressure drop)
Current
I (current)
Q (volume flow rate)
Resistance
R (resistance)
Rf (flow resistance)
L = 50 mm = 0.05 m
D = 0.25 mm = 2.5 10-4 m
Kerosine (at 40 C)
= 1.1 10-3 Ns/m2
= 1.35 10-6 m2/s
Caster Oil (at 40 C)
S.G. = 969
= 2.4 10-1 Ns/m2
= 2.48 10-4 m2/s
b) Find
Expression of flow resistance for laminar flow using the analogy to Ohm's law
Maximum pressure drop in the pipe for both kerosine and castor oil as working fluids
2. System Diagram
Ohm's Law
Fluid Flow
I
Q
R
V
Ghosh - 550
Page 16
3/10/2016
3. Assumptions
Steady state condition
Constant fluid properties
Laminar, fully developed flow
4. Governing Equations
1 p 2
2
r R
4 x
Q V dA
A
5. Detailed Solution
Flow Resistance
The problem gives the resistance to flow is R f
Volume flow rate Q can be calculated as follows:
p
.
Q
Ghosh - 550
Page 17
3/10/2016
Q V dA
A
u 2r dr
0
1 p 2
r R 2 2r dr
0 4 x
p R 3
r R 2 r dr
2 x 0
p 1 4 R 2 2
r
r
2 x 4
2
0
p R 4
2 x 4
p D
8 x 2
D 4 p
Q
128 x
Because
p
p p
p
p p p
lim
x
x
L
L
L
2
drop
x 0
be written as:
D 4 p D 4 p
Q
128 L 128 L
Finally, the flow resistance is
Rf
p
p
128L
4
Q D p
D 4
128L
This expression, flow resistance, is independent of the flow rate because there is no Q or V
included. It is consisted of only fluid property and the geometry of a pipe.
Rf
p
128 L
128 L 2 32 L
V A
V D
V
p R f Q
4
4
4
Q
D2
D
D
Ghosh - 550
Page 18
3/10/2016
Examining this pressure drop expression, one can say that the faster the average velocity is,
the bigger the pressure drop is; therefore, to obtain the maximum pressure drop, we want the
fastest flow velocity in laminar flow constraint.
VD
Re
VD
Re cr V
Re cr
D
Substituting this fastest average flow velocity into the pressure drop expression, we get:
32L Re cr 32L Re cr
D
D2
D3
Continue