TRANSPORT PHENOMENA

Prof. Sunando DasGupta (3 classes)

and
Prof. Manish Kaushal (1 class)
 TRANSPORT PHENOMENA
Basic Overview
Over the last few decades, the subject has revolutionized the way
Chemical Engineering science is taught. This course deals with the
unified treatment of the different transport processes, ubiquitous in
industry as well as in nature. Momentum, heat and mass transfer are
taught together due to the underlying similarities of the mathematics and
molecular mechanisms describing such processes.

Books

1. Transport Phenomena by Bird, Stewart and Lightfoot (Wiley)

2. Transport Phenomena Fundamentals by J. L. Plawsky (CRC Press)
3. Heat and Mass Transfer – A Transport Phenomena Approach by K. S. Gandhi (New Age)
4. Fluid Mechanics by R. W. Fox, P. J. Pritchard and A. T. McDonald (Wiley)
5. Introduction to Heat and Mass Transfer by F. P. Incropera & D. P. DeWitt (Wiley)

Grading: 3 tests of 1 hour duration, on Moodle

 Course Plan
Formulation and solution of momentum transfer in laminar flow.
Navier-Stokes equation and its applications
Boundary Layer concepts, boundary layer thicknesses (disturbance, displacement
and momentum), Blasius solution for flow over a flat plate
Use of momentum integral equation, turbulent boundary layers, fluid flow about
immersed bodies, drag
Formulation and solution of heat transfer in laminar flow
Development and use of energy equation
Transient conduction - lumped capacitance, analytical solutions and other methods.
Formulation and solution of mass transfer in laminar flow. Development and use of
species balance equation
Introduction to convective flow, natural convection, relevant examples from heat
and mass transfer
Mathematical treatment of the similarities between heat, mass and momentum
transfer, similarity parameters, and relevant analogies.
Solution of coupled heat, mass and momentum transfer problems based on analogy.
 TRANSPORT PHENOMENA

Application of Navier Stokes Equations

Continuity Equation Equation of Motion –
Navier Stokes equation
𝜕𝜕𝑣𝑣𝑥𝑥 𝜕𝜕𝑣𝑣𝑦𝑦 𝜕𝜕𝑣𝑣𝑧𝑧
+ + =0
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 ρ Dv / Dt = −∇p − [∇.τ ] + ρ g
Cartesian Coordinate – z component

𝜕𝜕 𝑣𝑣𝑧𝑧 𝜕𝜕 𝑣𝑣𝑧𝑧 𝜕𝜕 𝑣𝑣𝑧𝑧 𝜕𝜕 𝑣𝑣𝑧𝑧

𝜌𝜌 + 𝑣𝑣𝑥𝑥 + 𝑣𝑣𝑦𝑦 + 𝑣𝑣𝑧𝑧
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝜕𝜕 𝑃𝑃 𝜕𝜕 2 𝑣𝑣𝑧𝑧 𝜕𝜕 2 𝑣𝑣𝑧𝑧 𝜕𝜕 2 𝑣𝑣𝑧𝑧
= − + 𝜇𝜇 2
+ 2
+ 2
+ 𝜌𝜌 𝑔𝑔𝑧𝑧
𝜕𝜕𝜕𝜕 𝜕𝜕𝑥𝑥 𝜕𝜕𝑦𝑦 𝜕𝜕𝑧𝑧

Relevant Boundary Conditions

No slip at the liquid-solid interface

No shear at the liquid-vapor interface

6
 Transport Phenomena - Momentum Transfer

The record-write head for a computer disk storage system floats

above the spinning disk on a very thin film of air (the film
thickness is 0.5micron). The head location is 150mm from the
disk centre line, the disk spins at 3600rpm. The record-write
head is 10mm x 10mm square. Determine for standard air
(density=1.23kg/m3, viscosity=1.78x10-5 kg/(m.s)) in the gap
between the head and the disk -
(a) the Reynold's number of the flow, (b) the viscous shear stress
and (c) the power required to overcome viscous shear stress.

For such a small gap the flow can be considered as flow

between parallel plates.

7
 Flow between parallel plates

U
V dP/dx < 0 V

Couette Flow Pressure gradient Driven Flow

U
dP/dx < 0 V

Couette and Pressure gradient Driven Flow

8
Solution

Since the gap is too small, the situation can be approximated by the
case of flow between two parallel plates with one plate moving and
the other stationary with no applied pressure gradient.

Here,

V=Rω = 0.15 m x 3600 rev/min x 2 π rad/rev x 1/60 min/s=56.5 m/s

Re = (ρ V a)/µ = (V a)/(μ/ρ) = 56.5 m/s x 0.5 x 10 -6 m/(1.45 x 10 -5

m2 /s) = 1.95 (Laminar Flow)

For Couette flow between two parallel plates, with no applied

pressure gradient, the velocity profile will be linear with no-slip
boundary conditions at the bottom (stationary) and top (moving)
plates.
9
Therefore, the shear stress, τ, can be expressed as

τ = μdu/dy = μ V/a (linear velocity profile)

τ=1.78 x 10 – 5 kg /(m.s) x 56.5 m/s x 1/(0.5x10 – 6 m) = 2.01 KN/m2

Therefore,

Force, F = τ A = τ x (W x L), Torque, T = F x R = τ x (WL) x R

Power dissipation rate. P = T x ω = τ WL R ω

P = 2.01x10- 3 x 0.01x 0.01 x 0.15 x 3600 x 2π x 1/60 = 11.4 W

Thus the power required to overcome viscous shear stress is 11.4 W

10
 Transport Phenomena - Momentum Transfer
Problem 2

An oil skimmer uses a 5 m

wide x 6 m long moving
belt above a fixed platform
(θ = 30º) to skim oil off of
rivers (T = 10ºC). The belt
travels at 3 m/s. The
distance between the belt
and the fixed platform is 2
mm.
The belt discharges into an open tank on the ship. The fluid is actually a
mixture of oil and water. To simplify the analysis, assume crude oil
dominates. Find the discharge of oil into the tank on the ship, the force
acting on the belt and the power required (kW) to move the belt. For oil:
ρ = 860 kg/m3, viscosity, µ = 1x10 – 2N.s/m2 11
The x-component of the equation of motion is given as
𝜕𝜕 𝑢𝑢 𝜕𝜕 𝑢𝑢 𝜕𝜕 𝑣𝑣 𝜕𝜕 𝑤𝑤
𝜌𝜌 + 𝑢𝑢 + 𝑣𝑣 + 𝑤𝑤
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝜕𝜕 𝑃𝑃 𝜕𝜕 2 𝑢𝑢 𝜕𝜕 2 𝑢𝑢 𝜕𝜕 2 𝑢𝑢
= − + 𝜇𝜇 2
+ 2
+ 2
+ 𝜌𝜌 𝑔𝑔𝑥𝑥
𝜕𝜕𝜕𝜕 𝜕𝜕𝑥𝑥 𝜕𝜕𝑦𝑦 𝜕𝜕𝑧𝑧

12
Considering one-dimensional flow (v = 0 = w), no
applied pressure gradient, u is a function of only y, and a
steady state process, the above equation reduces to the
following governing equation

𝑑𝑑 2 𝑢𝑢 𝜌𝜌𝜌𝜌
= sinθ
𝑑𝑑𝑦𝑦 2 𝜇𝜇

The relevant boundary conditions (no slip) are

BC 1 At y = 0 u = 0
BC 2 At y = h, u = U
Therefore,

𝑑𝑑𝑑𝑑 𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌 𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌 𝑦𝑦 2

= 𝑦𝑦 + 𝐴𝐴 𝑎𝑎𝑎𝑎𝑎𝑎 𝑢𝑢 = + 𝐴𝐴𝐴𝐴 + 𝐵𝐵
𝑑𝑑𝑑𝑑 𝜇𝜇 𝜇𝜇 2
𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌 ℎ 𝑢𝑢
Applying the BC s  B = 0 and 𝐴𝐴 = +
𝜇𝜇 2 ℎ

Thus

𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌 ℎ𝑦𝑦 𝑦𝑦 2 𝑈𝑈𝑈𝑈

𝑢𝑢 = − − +
𝜇𝜇 2 2 ℎ

The volumetric flow rate per unit width of the film is given by

ℎ ℎ
𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌 ℎ𝑦𝑦 𝑦𝑦 2 𝑈𝑈𝑈𝑈 𝜌𝜌𝜌𝜌ℎ3 𝑈𝑈𝑈
𝑄𝑄 = � 𝑢𝑢𝑢𝑢𝑢𝑢 = − � − + 𝑑𝑑𝑑𝑑 = − 𝑠𝑠𝑠𝑠𝑠𝑠𝜃𝜃 +
𝜇𝜇 2 2 ℎ 12𝜇𝜇 2
0 0

14
 860 𝑥𝑥 9,81 𝑥𝑥 0.002 3 𝑥𝑥 0.002 𝑚𝑚2
𝑄𝑄 = − sin 30 + = 0.0027 (𝑝𝑝𝑝𝑝𝑝𝑝 𝑢𝑢𝑢𝑢𝑢𝑢𝑢𝑢 𝑤𝑤𝑤𝑤𝑤𝑤𝑤𝑤𝑤)
12 𝑥𝑥 10−2 2 𝑠𝑠

For the width of 5 m, the volumetric flow rate will be equal to 0.0135 m3/s.

Evaluate τ = μ du/dy at the moving belt from the expression for velocity as
𝑑𝑑𝑑𝑑 𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌 ℎ 𝑈𝑈
=− − 𝑦𝑦 +
𝑑𝑑𝑑𝑑 𝜇𝜇 2 ℎ
And at the moving belt,
𝑑𝑑𝑑𝑑 𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌𝜌 𝜇𝜇𝜇𝜇
𝜏𝜏 = 𝜇𝜇 𝑦𝑦 = ℎ = ℎ+
𝑑𝑑𝑑𝑑 2 ℎ

Using the values h = 0.002m, μ = 10-2N.s/m2 and U = 3 m/s

τ│at the belt = 19.21 N/m2

Power = (τ x L x W) U = 19.21 x 6 x 5 x 3 = 1.73 KW

15
 Transport Phenomena - Momentum Transfer
Problem 3
Consider two concentric
cylinders with a Newtonian
liquid of constant density,
ρ, and viscosity, μ,
contained between them.
The outer pipe, with radius,
Ro, is fixed while the inner
pipe, with radius, Ri, and
mass per unit length, m,
falls under the action of

gravity at a constant speed. There is no pressure gradient within the flow and no
swirl velocity component. Determine the vertical speed, V, of the inner cylinder as a
function of the following parameters: g, Ro, Ri, m, ρ, and μ. The space between the
two cylinders is not ‘too small’ compared to the radii of the cylinders
16
Solution

For the inner cylinder moving at constant velocity, the downward force
is exactly balanced by the viscous force as

𝜏𝜏𝑤𝑤 𝐴𝐴𝑤𝑤 │𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼𝐼 𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶𝐶 = 𝑚𝑚 𝐿𝐿 𝑔𝑔

The z component of the equation of motion in cylindrical coordinate is

𝜕𝜕𝑣𝑣𝑧𝑧 𝜕𝜕𝑣𝑣𝑧𝑧 𝑣𝑣𝜃𝜃 𝜕𝜕𝑣𝑣𝑧𝑧 𝜕𝜕𝑣𝑣𝑧𝑧

𝜌𝜌 + 𝑣𝑣𝑟𝑟 + + 𝑣𝑣𝑧𝑧
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝑟𝑟 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝜕𝜕𝜕𝜕 1 𝜕𝜕 𝜕𝜕𝑣𝑣𝑧𝑧 1 𝜕𝜕 2 𝑣𝑣𝑧𝑧 𝜕𝜕 2 𝑣𝑣𝑧𝑧
=− + µ 𝑟𝑟 + 2 + + 𝜌𝜌𝑔𝑔𝑧𝑧
𝜕𝜕𝜕𝜕 𝑟𝑟 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝑟𝑟 𝜕𝜕𝜃𝜃 2 𝜕𝜕𝑧𝑧 2

Cancelling the terms with the observations i) steady state, ii) vz is a

function of r only, not of z or θ, iii) no applied pressure gradient

The following simplified form of the NS equation can be obtained.

17
 Governing equation

1 𝑑𝑑 𝑑𝑑𝑣𝑣𝑧𝑧 𝜌𝜌𝜌𝜌
𝑟𝑟 =−
𝑟𝑟 𝑑𝑑𝑑𝑑 𝑑𝑑𝑑𝑑 𝜇𝜇

Upon integration (with C1 and C2 being the constants of integration)

𝜌𝜌𝜌𝜌𝑟𝑟 2
𝑣𝑣𝑧𝑧 = − + 𝐶𝐶1 ln 𝑟𝑟 + 𝐶𝐶2
4𝜇𝜇
Boundary Conditions
𝜌𝜌𝜌𝜌
vz = V at r = Ri => 𝑉𝑉 = − 𝑅𝑅𝑅𝑅 2 + 𝐶𝐶1 ln 𝑅𝑅𝑖𝑖 + 𝐶𝐶2
4𝜇𝜇

𝜌𝜌𝜌𝜌
vz = 0 at r = Ro => 0=− 𝑅𝑅𝑜𝑜 2 + 𝐶𝐶1 ln 𝑅𝑅𝑜𝑜 + 𝐶𝐶2
4𝜇𝜇

18
 𝜌𝜌𝜌𝜌 2 2 𝑅𝑅𝑖𝑖
Therefore 𝑉𝑉 = 𝑅𝑅𝑅𝑅 − 𝑅𝑅𝑅𝑅 + 𝐶𝐶1 ln
4𝜇𝜇 𝑅𝑅𝑜𝑜

1 𝜌𝜌𝜌𝜌
𝐶𝐶1 = 𝑉𝑉 − 𝑅𝑅𝑅𝑅 2 − 𝑅𝑅𝑅𝑅 2
𝑅𝑅𝑖𝑖 4𝜇𝜇
ln
𝑅𝑅𝑜𝑜
𝑑𝑑𝑣𝑣𝑧𝑧 𝜌𝜌𝜌𝜌𝜌𝜌 𝐶𝐶1 𝑑𝑑𝑣𝑣𝑧𝑧 𝜌𝜌𝜌𝜌𝜌𝜌 𝐶𝐶1 𝜇𝜇
=− + 𝜏𝜏𝑟𝑟𝑟𝑟 = 𝜇𝜇 =− +
𝑑𝑑𝑑𝑑 2𝜇𝜇 𝑟𝑟 𝑑𝑑𝑑𝑑 2 𝑟𝑟
Since the force on the inner cylinder = force due to gravity

𝜏𝜏│𝑟𝑟=𝑅𝑅𝑅𝑅 2𝜋𝜋𝑅𝑅𝑖𝑖 𝐿𝐿 = 𝑚𝑚 𝐿𝐿 𝑔𝑔

Upon substitution for the expression of the shear stress

𝑅𝑅𝑖𝑖 𝜌𝜌𝜌𝜌𝑅𝑅𝑖𝑖 𝑚𝑚𝑚𝑚 𝜌𝜌𝜌𝜌

𝑉𝑉 = 𝑅𝑅𝑖𝑖 ln − − 𝑅𝑅𝑖𝑖 2 − 𝑅𝑅0 2
𝑅𝑅𝑜𝑜 2𝜇𝜇 2𝜋𝜋𝑅𝑅𝑖𝑖 𝜇𝜇 4𝜇𝜇
19
 Order of magnitude analysis of NS Equation

Is it possible to identify the relative magnitudes of the different terms

(even approximately)?

It may then be possible to neglect the term(s) that may not play a
crucial role in the transport process thereby simplifying NS equations.
 Flow between two parallel disks with liquid entry through
a small hole at the centre of the top plate

Flow In
Disk radius = r2
Hole radius = r1
Z=+b

Z=-b

Flow Out
The pressure is p1 at r = r1 and p2 at r = r2
 ∂ρ 1 ∂ 1 ∂ ∂
+ ( ρ r vr ) + ( ρ vθ ) + ( ρ vz ) =
0
∂t r ∂r r ∂θ ∂z

 ∂vr ∂vr vθ ∂vr vθ 2 ∂vr 

ρ + vr + − + vz 
 ∂ t ∂r r ∂θ r ∂z 

∂P  ∂ 1 ∂  1 ∂ vr 2 ∂vθ ∂ vr 
2 2
=
− +µ  ( rvr )  + 2 2 − 2 + 2  + ρ gr
∂r  ∂r  r ∂r  r ∂θ r ∂θ ∂z 
 The governing equation

∂ vr ∂P ∂ 2 vr φ ( z)
ρ vr =
− +µ vr =
∂r ∂r ∂ z2 r

Major assumption on the nature of the flow

Assume a constant applied pressure difference
Boundary Conditions

∆P b 2  z 
2

=
vr ( r , z ) 1 −   
r2   b  
2 µ r ln  
r1
 +b +b
4 π ∆ P b3
= π ∫ r vr dz 2 π=
Q 2= ∫ φ ( z ) dz
r2
−b −b 3 µ ln
r1
An incompressible fluid flows between two porous, parallel flat plates
as shown in the figure. An identical fluid is injected at a constant
speed V through the bottom plate and simultaneously extracted from
the upper plate at the same velocity. Assume the flow to be steady,
fully-developed, the pressure gradient in the x-direction is a constant,
and neglect body forces.

Determine expressions for the y component of velocity.

Show that the x component of velocity can be expressed as
𝜌𝜌 𝑉𝑉 𝑦𝑦
ℎ 𝜕𝜕𝜕𝜕 1 − exp 𝑦𝑦
𝜇𝜇
𝑢𝑢𝑥𝑥 = −
𝜌𝜌𝜌𝜌 𝜕𝜕𝜕𝜕 𝜌𝜌 𝑉𝑉 ℎ ℎ
1 − exp 𝜇𝜇
The equation of continuity for fully developed, steady flow in x-direction

𝜕𝜕𝑣𝑣𝑥𝑥 𝜕𝜕𝑣𝑣𝑦𝑦 𝜕𝜕𝑣𝑣𝑧𝑧

+ + =0
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕

The x component of the NS equation

𝜕𝜕 𝑣𝑣𝑥𝑥 𝜕𝜕 𝑣𝑣𝑥𝑥 𝜕𝜕 𝑣𝑣𝑥𝑥 𝜕𝜕 𝑣𝑣𝑥𝑥

𝜌𝜌 + 𝑣𝑣𝑥𝑥 + 𝑣𝑣𝑦𝑦 + 𝑣𝑣𝑧𝑧
𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕 𝜕𝜕𝜕𝜕
𝜕𝜕 𝑃𝑃 𝜕𝜕 2 𝑣𝑣𝑥𝑥 𝜕𝜕 2 𝑣𝑣𝑥𝑥 𝜕𝜕 2 𝑥𝑥
= − + 𝜇𝜇 2
+ 2
+ 2
+ 𝜌𝜌 𝑔𝑔𝑥𝑥
𝜕𝜕𝑥𝑥 𝜕𝜕𝑥𝑥 𝜕𝜕𝑦𝑦 𝜕𝜕𝑧𝑧
 The governing equation

∂ vx ∂P ∂ 2 vx
ρ vy =
− +µ
∂y ∂x ∂ y2

Assumptions:

Constant pressure gradient in the x direction

Boundary Conditions
A thrust bearing as shown in the figure is lubricated by pumping oil at a high
pressure of p0. The angular velocity is equal to ω. Note that under laminar
conditions, both Vr and Vθ in the thin gap will be non-zero and p is a function of
r only. Neglect convective and body force terms in equations of motion. You
may also assume that the pressure at r = R is equal to patm, whereas the pressure
from r = 0 till r = R0 is equal to p0.

(i) Start with the equation of continuity to obtain the functional form of Vr.

ω r  z  satisfies the θ component of the NS equation.

(ii) Show that= Vθ 1 + 
2  h
(iii) Write the r component of NS equation to show that r ∂ p is a constant.
∂r
(iv) Evaluate the pressure distribution using the boundary conditions

(v) Find the vertical load the bearing can support and the flow rate of oil
required.
(i) Start with the equation of continuity to obtain the functional form of Vr.

The Equation of Continuity

 ωr z  satisfies the θ component of the NS equation.
(ii) =
Show that Vθ  1 + 
2  h
 ∂p
(iii) Write the r component of NS equation to show that r is a constant
∂r
(iv) Evaluate the pressure distribution using the boundary conditions
(v) Find the vertical load the bearing can support and the flow rate of
oil required.
Consider the spin-coating process used to coat silicon wafers with photoresist. The
process is designed to produce a very thin, uniform coating by spinning a viscous,
Newtonian, liquid onto a substrate (wafer). The process has angular symmetry, the
rotation rate is constant, and since the film is thin, there are no real pressure gradients or
fluid accelerations and body forces to speak of. The thin film also moves with the
substrate as if it were a rigid body, vθ ≠ f(z)

a) Show that the continuity and momentum equations reduce to:

vθ2
−ρ =µ
d 2 vr
;
1∂
( r vr ) + ∂ v z = 0
r dz 2
r ∂r ∂z

b) What are the boundary conditions for this problem?

c) Solve the equations for vr and vz.
d) The velocity, vz, at the film/air interface is just the change in film thickness
with time. Use this to obtain a differential equation for h, and integrate this
equation to obtain the solution as:

1 1 1   2ρω 2 
 2 + 2  =   t
2 h ho   3µ 
Where ω is the rate of rotation (vθ = 2πrω) and h0 is the initial height of the film.
 In this case, the following conditions will hold:

vθ= Cons tan t , g=

r g=
θ g=z 0, ρ , µ constant
∂ ∂ ∂P ∂P ∂P
= 0, = 0, = = = 0
∂θ ∂t ∂r ∂θ ∂z

The continuity and momentum equations in cylindrical coordinates become

1 ∂ ∂ vz
( r)
r v + =
0
r ∂r ∂z
Since the fluid moves essentially as a rigid body, vr = f (r, z). This film is
very thin with respect to the radial or angular dimensions so that
vθ , vr  vz

Thus since the angular velocity is constant and vz is small, we only need to deal
with the momentum equation for vr
Neglecting all the obvious terms,

 ∂ vr vθ2 ∂v   ∂ 1 ∂  ∂ 2 vr 
ρ  vr −= + vz r  µ  ( r vr )  + 2 
 ∂r r ∂z   ∂r  r ∂r  ∂z 
  ∂ 1 ∂  ∂ 2 vr   ∂ vr vθ2 ∂ vr 
µ  ( r vr )  + = 2 
ρ  vr − + vz 
∂
 r r ∂ r  ∂ z   ∂ r r ∂ z 

Since the film is very thin, vr changes rapidly with the film thickness, much
more rapidly than it changes with the radial position. vz is quite small as
compared to the other components

Thus 1st term on LHS is neglected w.r.t. the second term.

For the same reason, 1st and 3rd terms on the rhs can be neglected.

In one case, vz is very small and in the second

∂ vr is also very small.
∂r

vθ is quite large
Thus
d 2 vr ρ vθ2
µ 2 + =
0
dz r
The boundary conditions are

=
z 0; v=
r vθ= 0
∂ vr
=z h=
; 0
∂z

Noting that vθ = 2π rω , the above equation is solved

 4 π 2ω 2  z2 
=vr r  hz − 
γ  2 

Using the continuity equation

1∂ ∂ vz
( r)
r v + =
0
r ∂r ∂z

4 π 2ω 2  z 3 2 
vz  −hz 
γ 3 

BC: vz = 0 at z = 0
 1
h
π 2ω 2 h3 d h
vz =
h0∫ vz dz =
−
γ
=
dt

=t 0=
At h h0
1 1 2π 2 ω 2
− 2 = t
h 2
h0 γ

This can also be solved by assuming vz at z = h to be equal to dh/dt

Problem: One disk approaching another displacing a liquid in between

Two parallel, plane circular disks (of radius R) lie one above the other a
small distance apart. The space between them is filled with a liquid. The
upper disk approaches the lower at a constant (small) velocity U,
displacing the liquid. The imposed pressure at r = R is po and pr is not a
function of z. Simplify the basic equations (continuity and motion, r
component) using an analysis based on your understanding and estimate
of the magnitudes of the different terms in the original equations.

Solve the these equations with proper boundary conditions to evaluate a

expression of the resistance to motion on the moving disk in terms of the
separation between the two plates. Take the origin to be at the centre of the
lower (stationary) plate.
Assumptions and special considerations

1. R is very large compared to h

2. Leakage rate is small (Small vr)

3. ∂p/∂z is negligible

4. vz is very small as compared to vr

5. ∂vz/∂z may not be small (as z is very small)

𝜕𝜕 𝜕𝜕𝑣𝑣𝑧𝑧
6. ( )
𝜕𝜕𝑧𝑧 𝜕𝜕𝜕𝜕
can be appreciable

7. vr is not small as compared to vz

𝜕𝜕𝑣𝑣𝑟𝑟 𝜕𝜕𝑣𝑣𝑟𝑟
8. 𝜕𝜕𝑟𝑟
is small compared to 𝜕𝜕𝑧𝑧

9. Vθ = 0, no dependence of vr , vz
The equation of continuity

∂ρ 1 ∂ 1 ∂ ∂v
+ ( r vr ) + ( ρ vθ ) + z =
0 ( A)
∂t r ∂ r r ∂θ ∂z

v r ∂ vr ∂ vz
+ + =
0
r ∂r ∂z

The terms in the equation

Relative magnitudes of 1,2 w.r.t. 3

The equation of motion (r component)

1 ∂p  1 ∂ vr vr ∂ 2 vr ∂ 2 vr 
0=
− +γ  − 2 + 2 + 2
ρ ∂r  r ∂r r ∂r ∂z 
 1 ∂p  1 ∂ vr vr ∂ 2 vr ∂ 2 vr 
0=
− +γ  − 2 + 2 + 2
ρ ∂r  r ∂r r ∂r ∂z 

Compare terms to arrive at the following equation

1 ∂p ∂ 2 vr
0=
− +γ ( B)
ρ ∂r ∂ z2

1 ∂P 2
=
vr z + C1 z + C2
2µ ∂ r

BC : z=0; vr =0 ⇒ C2 =0
1 dp
z=
h; vr =
0 ⇒ C1 =
− h
2µ dr
 1 ∂P
=vr ( z − h ) z (C )
2µ ∂ r

1 ∂ ∂ vz
( r)
r v + =
0 ( A)
r ∂r ∂z

Putting the expression of vr from (C ) to (A) and integrating w.r.t. z and as

BC : z= 0; vr= 0= vz
z = h; vr = 0; vz = − U

1 d  dp 1 h3 1 d  d p 
h h
1 d
U = ∫ r vr dz =  r  ∫ z ( z − h ) dz =
− r 
r dr 0 r dr  dr  2µ 0 12 µ r d r  d r 
 h3 1 d  d p 
U= − r 
12 µ r d r  d r 
Integrating w.r.t. r
12 µ r 2 dp
− 3 =
r + C1
h 2 dr

As dp/dr is finite at r = 0 C1= 0 and p = po at r = r

3µ U 2
p=
− po
h3
( R − r 2
)
3µ U 2
R R

∫ 2π rdr ( p −= 3 ( )
=
Force,F po ) 2π ∫ R r − r 3
dr
0 0
h

3 π µ U R4
Force, F =
2 h3
A fluid (of constant density ρ) is in
incompressible, laminar flow through a
tube of length L. The radius of the tube
of circular cross section changes linearly
from R0 at the tube entrance (z = 0) to a
slightly smaller value RL at the tube exit
(z = L).

Using the lubrication approximation,

determine the mass flow rate vs. pressure
drop (ΔP) relationship for a Newtonian
fluid (of constant viscosity μ).

The approximation where a flow

between non-parallel surfaces is treated
locally as a flow between parallel
surfaces is commonly called the
lubrication approximation because it is
often employed in the theory of
lubrication.
The mass flow rate versus the pressure drop (W versus ΔP) relationship for a Newtonian
fluid in a circular tube is given by (where ΔP = P0 – PL)

π ∆ P R4 ρ ∆P 8 µ W
=W = or
8µ L L π ρ R4
 π ∆ P R o4 ρ  3 ( λ − 1)  RL
=W =  −3 
; λ
8µ L  1 − λ 
Ro

Taper correction to Hagen-Poiseuille equation

Lubrication flows are characterized by incompressible fluids to thin gaps e.g.,
the layer of water between the ice skate and the ice, the oil that lubricates the
moving parts of an internal combustion engine etc. The fluid flows at
relatively small velocities and so the inertial terms in the Navier-Stokes
equations are insignificant compared to the diffusive terms. Since the gaps are
so thin, lubrication flows can be treated as two-dimensional. Finally, if the
gap is very thin, there is little effect of gravity.
Consider the lubrication flow
occurring between the fixed
and moving components of a
bearing as shown in the figure
in presence of a pressure
gradient with one of the plate
moving with a constant
velocity Uo. Assume that
there is a constant volumetric
flow rate V of the lubricant in
the system as a result of this
motion of the plate and the
applied pressure gradient.
Solve the governing equations with appropriate boundary conditions to
obtain the x-component of the velocity, u, as a function of the pressure
gradient and H(x).

Use relevant conditions to obtain an expression for dP/dx and local

pressure P(x) in terms of H(x), Uo and V.

Finally obtain expressions for the (i) volumetric flow rate V and (ii) the
load the bearing can support in terms of the system parameters (Hi, Ho, θ,
Uo) and µ.
A cell separation (fractionation) system is based on cell density. Cells are
injected at the centre of a tube of radius R, and are carried by fluid flowing at
a flow rate of Q. Dense cells fall quickly under the action of gravity, adhere
to the tube wall and hence do not pass out of the tube. Assume that the
concentration of the cell is low enough that the laminar flow in the tube is
not perturbed by the presence of the cells. Let the cells be spherical, with
radius ‘a’ and let them have a density ρ + Δρ, where ρ is the density of the
flowing fluid.
Find the axial distance (L) that a cell travels
before it hits the bottom wall. You need to use
the fact that, in fully developed flow in the
presence of gravity, the pressure distribution
in the vertical distance is hydrostatic and the
axial velocity profile is the familiar parabolic
shape. You may assume that the cell is
spherical, that it reaches its terminal (falling)
velocity nearly immediately after injection
into the tube, and that it gets carried axially at
the local fluid velocity in the tube. State any
other assumptions you make.
We need to find the axial distance L that the cell travels before it hits the bottom wall.

The expressions of the forces acting on the spherical cell are:

4
Gravitational force =
F g m=
cell g π a3 ( ρ + ∆ ρ ) g
3

4
Buoyancy force =
F b m=
fluid g π a3 ρ g
3

The cell is moving downward with a

Drag force F d = 6 π µ a vr velocity vr in the r direction

4 4
π a ( ρ=
3
+ ∆ρ ) g π a 3 ρ g + 6 π µ a vr
3 3
2 a2 g ρ
vr =
9µ
Find the radial location as a function of time as

2 a2 g ∆ ρ t
t
r ( t ) ∫=
= vr dr
0
9µ

The time T, when the cell will hit the wall

R 9µ R
= =
T
vr 2 a 2 g ∆ ρ

∆ P R2   r 2 
=vx 1 −   
4µ L   R  
∆P R 2 ∆P π R 4
=vx =; Q π=
R 2 vx
8µ L 8µ L
2 Q ∆P R 2
=
πR 2
4µ L
 2 Q   r (t )  
2

=vx 1 −   
πR   R  
2
 

2 Q  [ r (t ) ] 
T T 2

=L ∫=vx dt 2 ∫
1−  dt
0
π R 0 
R 2


2Q   2 a2 g ∆ ρ  T 3 
2

=L T −   
π R 
2
 9 µ R  3 
2

9µ R
As T = 2
2a g ρ
6µQ
L=
π R a2 g ∆ ρ
Epithelium is one basic type of animal tissue, which lines the cavities and surfaces
of structures throughout the body. Epithelial layers contain no blood vessels, so they
must receive nourishment via diffusion of substances from the underlying
connective tissue, through the basement membrane.

An apparatus has been built for testing the effect of various drugs on the rate at
which an epithelium can pump fluid from its luminal side ( the side facing the fluid )
to its basal side (which lies on the channel wall). The cells line the top and bottom
surface of a flow channel that has a separation of h (from top plate to bottom plate;
ignore the thickness of cells), a length L, and a depth into the page of W. Each of
these walls is porous so that any fluid pumped by the cells can leave the channel. Let
each cell layer (top and bottom) pump fluid at a rate of q per unit area of the channel
walls (q has thus units of length/time). The height of the channel is much less than its
length (h≪L).
Fluid enters the channel at the left at a flow rate Q0 and a gauge pressure of P0. Because
of the pumping action of the cells, the flow rate through the channel decreases as a
function of x, the distance from the beginning of the channel. To determine the rate at
which the cells are pumping fluid out of channel, the channel is instrumented with
pressure transducers that can measure P(x). We would like to use this information to find
the rate at which the cells pump fluid. The fluid in the channel has a density of ρ and a
viscosity of μ. The flow is dominated by viscous effects and is steady.

(a) Find the pressure distribution, P(x) , in the flow channel if q=0.

(b) Find the pressure distribution P(x) in the channel for q≠0.

(c) Given that P(x=L)=Pe , find q .

(d) Find the criterion necessary for the assumption that viscous flow dominates to be valid.

All answers must be given in terms of the known quantities e.g., x,L,W,h,Q0 , P0, Pe, ρ and
μ (not all of these parameters need necessarily be used).
Find the pressure distribution, P(x) , in the flow channel if q=0
Find the pressure distribution P(x) in the channel for q≠0.
Given that P(x=L)=Pe , find q .
Find the criterion necessary for the assumption that viscous flow dominates to be valid.
The Equation of Continuity
The Equation of Motion for a Newtonian Fluid with constant μ, ρ

The Navier Stokes Equation