FLUID DYNAMICS: Complete Syllabus Index

I. The Heavy Hitters

Part 1: Shell Momentum Balance (Falling Film & Pipe Flow)

• Detailed derivations of Velocity Profiles and Flow Rates.

Part 2: The General Equation of Motion

• Continuity, Stress Tensor, Cauchy Eq, Navier-Stokes.

II. The Mathematical Engine

Part 3: Tensor Analysis for Fluids

• Dyadics, Stress Tensors, Generalised Newton’s Law.

Part 4: Conservation Laws & Continuity

• Reynolds Transport, Convective Derivatives.

III. Microscopic Origins

Part 5: Kinetic Theory & Transport

• Viscosity origin, Temperature dependence.

IV. Exam Strategy Suite

Part 9: Strategy & Formula Sheet
Part 10: Essential Visuals
Part 11: Comprehensive Problems
Part 12: Conceptual Intuition
Part 13: HOTS Questions
Part 14: Deep Dive

1
Part 1: Shell Momentum Balance
Topic A: Falling Film (Flow on an Inclined Plane)
The Scenario: Imagine pouring honey down a cutting board that is tilted. The
liquid forms a smooth layer of constant thickness and slides down. We want to
find the equation for how fast it moves at every depth.

Step 1: The Physical Setup (The ”Definitions”)

Write this paragraph first to secure your ”definition” marks.
Consider a fluid flowing down an inclined plane of infinite width W and length
L.

• Geometry: We use a Cartesian Coordinate System (x, y, z).

– z-axis: The direction of flow (down the ramp).

– x-axis: Perpendicular to the ramp. We define x = 0 at the liquid-air
surface and x = δ at the solid wall.

• The Fluid: It is Newtonian (viscosity η is constant) and Incompressible

(density ρ is constant).

• Assumptions:

1. Laminar Flow: The fluid moves in smooth parallel layers, not mixing
(No turbulence).
2. Steady State: Nothing changes with time (∂/∂t = 0).
3. Fully Developed: The velocity profile vz does not change as we go
down the ramp. It only depends on the depth x. Thus, vz = vz (x).
4. No Edge Effects: The width W is very large, so we ignore what
happens at the side edges (vy = 0).

Step 2: Choosing the ”Shell”

To apply Newton’s Second Law (F = ma), we need to pick a specific ”chunk” of
fluid. We choose a Rectangular Shell (a thin slice) oriented parallel to the plate.

• Location: Between x and x + ∆x.

• Dimensions: Length L (along flow), Width W , Thickness ∆x.

• Volume: V = L · W · ∆x.

2
Step 3: Momentum Balance (Force Balance)
Since the flow is steady (acceleration is zero), the sum of forces acting on our shell
must be zero. X
Fz = 0
The forces acting in the z-direction are:
1. Viscous Shear Force (Molecular Transport): Fluid layers rub against
each other.

• In at x: Momentum enters the shell from the layer above.

Force In = (Area) × (Stress) = (LW ) · τxz |x

• Out at x + ∆x: Momentum leaves the shell to the layer below.

Force Out = (LW ) · τxz |x+∆x

2. Gravity Force: Gravity pulls the shell down the ramp. We need the
component of g parallel to the ramp (g cos β).

Fg = (Mass) × gz = (ρ · Volume) · (g cos β)

Fg = ρ(LW ∆x)g cos β

3. Pressure Force: Since the top of the liquid is open to air, the pressure is
atmospheric everywhere along the length L. The pressure at z = 0 pushes same
as z = L pushes back. Net Pressure Force = 0.
The Balance Equation:

(Shear In) − (Shear Out) + (Gravity) = 0

(LW )τxz |x − (LW )τxz |x+∆x + (LW ∆x)ρg cos β = 0

Step 4: Deriving the Differential Equation

Divide the entire equation by the volume (LW ∆x):

τxz |x − τxz |x+∆x

+ ρg cos β = 0
∆x
Rearrange to look like a derivative definition:
 
τxz |x+∆x − τxz |x
− + ρg cos β = 0
∆x
dτxz
Take the limit as ∆x → 0. The term in brackets becomes the derivative dx
.

dτxz dτxz
− + ρg cos β = 0 =⇒ = ρg cos β - - - #
dx dx

3
Step 5: Solving for Shear Stress
Integrate the equation with respect to x:
Z
τxz = (ρg cos β)dx = (ρg cos β)x + C1

To find the constant C1 , we use a Boundary Condition.

• Condition: At the free surface (x = 0), the air exerts negligible friction on
the fluid.

• Math: At x = 0, τxz = 0.

0 = (ρg cos β)(0) + C1 =⇒ C1 = 0

So, the shear stress distribution is linear:

τxz (x) = ρgx cos β - - - #

Step 6: Solving for Velocity

Now we convert stress to velocity using Newton’s Law of Viscosity.
dvz
τxz = −η
dx
Substitute our stress equation into this:
dvz
−η = ρgx cos β
dx
Rearrange to solve for the velocity gradient:
 
dvz ρg cos β
=− x
dx η

Integrate with respect to x:

x2
 
ρg cos β
vz (x) = − + C2
η 2

To find C2 , we use the second Boundary Condition.

• Condition: ”No Slip”. At the solid wall (x = δ), the fluid sticks to the wall.

• Math: At x = δ, vz = 0.
ρg cos β 2
0=− δ + C2
2η
ρg cos β 2
C2 = δ
2η

4
 Now plug C2 back into the velocity equation:
ρg cos β 2 ρg cos β 2
vz (x) = − x + δ
2η 2η
Factor out the common terms to make it look pretty:

ρgδ 2 cos β
  x 2 
vz (x) = 1− ---#
2η δ

Interpretation: This is a parabola. Velocity is max at x = 0 and zero at x = δ.

Step 7: Volumetric Flow Rate (Q)

Q is the total volume of fluid passing a cross-section per second. We sum up the
velocity over the area. Z
Q = vz dA

Since velocity only varies with x, the area element is a strip of width W : dA =
W dx. Z δ
Q= vz (x)W dx
0
2 Z δ
x2

ρgδ cos β
Q=W 1 − 2 dx
2η 0 δ
Do the integral:
Z δ δ
x2 x3 δ3
   
δ 2
1− 2 dx = x − 2 = δ− 2 =δ− = δ
0 δ 3δ 0 3δ 3 3

Combine everything:
ρgδ 2 cos β
  
2
Q=W δ
2η 3
W ρgδ 3 cos β
Q= ---#
3η

5
Topic B: Flow Through a Circular Tube (Hagen-Poiseuille)
The Scenario: Water flowing through a horizontal pipe. The walls slow it down,
the pressure pushes it forward.

Step 1: Physical Setup (Definitions)

Consider a fluid flowing through a circular tube of Radius R and Length L.

• Coordinate System: Cylindrical (r, θ, z).

– z: Direction of flow (along the pipe).

– r: Radial direction (0 at center, R at wall).

• Forces:

– No gravity effect in the flow direction (horizontal).

– Flow is driven by a Pressure Difference: ∆P = P0 − PL .

• Assumptions: Laminar, Steady, Incompressible, Axisymmetric (vz depends

only on r), No Swirl (vθ = 0).

Step 2: The Shell Balance

We define our control volume as a Cylindrical Shell (a hollow tube within the
flow).

• Location: Radius r, Thickness ∆r, Length L.

• Surface Area (of the cylinder wall): Acyl = 2πrL.

• Cross-Section Area (face of the ring): Af ace = 2πr∆r.

Force Balance: 1. Pressure Force: Pressure pushes on the entrance face

(P0 ) and pushes back on the exit face (PL ).

FP = (P0 − PL ) × (Area of Face) = (P0 − PL )(2πr∆r)

2. Viscous Shear Force: Shear acts on the cylindrical surface area (2πrL).

• In at r: (2πrLτrz )|r

• Out at r + ∆r: (2πrLτrz )|r+∆r

• Note: The area 2πrL changes as r changes! This is why r stays inside the
derivative.

6
 Balance Equation:

(Shear In − Shear Out) + (Pressure Force) = 0

[(2πrLτrz )|r − (2πrLτrz )|r+∆r ] + (P0 − PL )2πr∆r = 0

Step 3: Differential Equation

Divide by volume element 2πL∆r: (Note: We approximate r∆r in the pressure
term as just the volume component contribution).
 
−(rτrz )|r+∆r + (rτrz )|r P0 − P L
+ r=0
∆r L

Take limit ∆r → 0:  
d ∆P
− (rτrz ) + r=0
dr L
 
d ∆P
(rτrz ) = r ---#
dr L

Step 4: Solving for Stress

Integrate with respect to r:

r2
 
∆P
rτrz = + C1
L 2

Divide by r:
∆P C1
τrz = r+
2L r
Boundary Condition 1 (Finite Stress): At the center of the pipe (r = 0),
the shear stress cannot be infinity. If C1 ̸= 0, the term C1 /r blows up. Therefore,
C1 must be 0.
∆P
τrz = r ---#
2L

Step 5: Solving for Velocity

Use Newton’s Law: τrz = −η dv
dr
z
.

dvz ∆P
−η = r
dr 2L
dvz ∆P
=− r
dr 2ηL
Integrate:
∆P r2 ∆P 2
vz (r) = − + C2 = − r + C2
2ηL 2 4ηL

7
 Boundary Condition 2 (No-Slip): At the pipe wall (r = R), fluid velocity
is zero.
∆P 2 ∆P 2
0=− R + C2 =⇒ C2 = R
4ηL 4ηL
Substitute C2 back:

(P0 − PL )R2
  r 2 
vz (r) = 1− ---#
4ηL R

Step 6: Mass Flow Rate (w)

Mass flow rate is Density × Volumetric Rate.
Z
w = ρvz dA

For a pipe, the area element is a ring: dA = 2πrdr.

Z R
w= ρvz (r) · 2πrdr
0

R
2πρ∆P R2
Z   r 2 
w= 1− rdr
4ηL 0 R
Let u = r/R. The math simplifies to:

π(P0 − PL )R4 ρ
w= ---#
8ηL

Part 2: The General Equation of Motion

Exam Tip: Professor-Proofing
Instructor Requirement: When deriving these equations, you MUST
define the control volume first. Do not just write ∇ · ⃗v = 0. Start with
”Consider a fixed cubic element of fluid...”

1. The Equation of Continuity (Mass Conservation)

Physical Setup: Consider a small stationary volume element fixed in space (Eu-
lerian approach) with dimensions ∆x, ∆y, ∆z. The fluid has density ρ(x, y, z, t)
and velocity ⃗v = vx î + vy ĵ + vz k̂.
The Principle:

Rate of Mass Accumulation = Rate of Mass In − Rate of Mass Out

Derivation Steps:

8
 1. Accumulation: The rate of change of mass inside the volume.
∂ ∂ρ
Accumulation = (ρ∆x∆y∆z) = ∆x∆y∆z
∂t ∂t

2. Net Flow (x-direction): Mass flux is ρvx .

In − Out = [(ρvx )|x − (ρvx )|x+∆x ] ∆y∆z

Using Taylor expansion or derivative definition:

∂(ρvx )
≈− ∆x∆y∆z
∂x

3. Total Balance: Summing x, y, and z directions and dividing by volume

∆V :  
∂ρ ∂(ρvx ) ∂(ρvy ) ∂(ρvz )
=− + +
∂t ∂x ∂y ∂z

Vector Form (The Result):

∂ρ
+ ∇ · (ρ⃗v ) = 0 — #
∂t
Special Case: Incompressible Fluid If density ρ is constant in time and
space, then ∂ρ
∂t
= 0.
∇ · ⃗v = 0 — #

2. Momentum Balance (Equation of Motion)

This corresponds to Pages 54-56 of your notes. We apply Newton’s Second Law
to the moving fluid.
Physical Setup: Same control volume. We use the **Substantial Derivative**
(D/Dt) to follow a specific fluid particle.

A. Inertial Term (Ma)

The acceleration of a fluid particle is not just the change with time, but also the
change due to moving to a new location.
D⃗v ∂⃗v
⃗a = = + (⃗v · ∇)⃗v
Dt ∂t
Force = Mass × Acceleration (per unit volume):
 
∂⃗
v
F⃗inertial = ρ + ⃗v · ∇⃗v
∂t

9
B. Surface Forces (Stress Tensor)
Forces acting on the surfaces of the element are due to:

1. Pressure (P ): Normal, compressive force.

2. Viscous Stress (τ ): Friction/shear forces.

The net force per unit volume is the negative divergence of the Total Stress Tensor
Π:
F⃗surf ace = −[∇ · Π] = −∇P − [∇ · τ ]

C. Body Forces
Usually Gravity. Force per unit volume:

F⃗body = ρ⃗g

D. The Cauchy Momentum Equation

Combining Inertia = Surface Forces + Body Forces:
D⃗v
ρ = −∇P − [∇ · τ ] + ρ⃗g —#
Dt
Note: This equation is valid for ANY continuum fluid (Newtonian or Non-Newtonian).

3. The Navier-Stokes Equation

To get the famous Navier-Stokes equation, we must assume the fluid is **Newto-
nian** and **Incompressible**.
Step 1: Constitutive Law We use the Generalized Newton’s Law (from Part
3):
τ = −µ[∇⃗v + (∇⃗v )T ]
Step 2: Divergence of Stress We need to calculate ∇ · τ .

∇ · τ = −µ[∇ · (∇⃗v ) + ∇ · (∇⃗v )T ]

= −µ[∇2⃗v + ∇(∇ · ⃗v )]
Since the fluid is Incompressible, ∇ · ⃗v = 0. The second term vanishes.

∇ · τ = −µ∇2⃗v

Step 3: Substitution Substitute this back into the Cauchy Equation.

D⃗v
ρ = −∇P − (−µ∇2⃗v ) + ρ⃗g
Dt

10
 The Navier-Stokes Equation (Vector Form):
 
∂⃗v
ρ + ⃗v · ∇⃗v = −∇P + µ∇2⃗v + ρ⃗g —#
∂t

4. Special Case: Creeping Flow (Stokes Flow)

Condition: Very low velocity, high viscosity (Re ≪ 1). The inertial terms (LHS)
are negligible compared to viscous terms.
D⃗v
ρ ≈0
Dt
Resulting Equation:

0 = −∇P + µ∇2⃗v + ρ⃗g —#

This is a linear equation and is used for things like the ”falling raindrop” example
in your notes (Page 43).

Part 3: Tensor Analysis for Fluids

1. Why do we need Tensors?
The Concept: In basic physics, Force is a vector (Magnitude + Direction). In
fluids, we deal with Stress (Force per unit Area). Stress is more complex because
it depends on two directions:

1. The direction of the Force.

2. The orientation (Normal vector) of the Surface it acts on.

A quantity with two associated directions is a Second-Order Tensor.

A. Notation Definitions
• Scalar (0th Order): Temperature T , Pressure P . (1 component).

• Vector (1st Order): Velocity ⃗v . (3 components).

• Tensor (2nd Order): Stress τ . (9 components).

The Dyadic Product (⃗a⃗b): Your notes use ”Dyadic” notation often. If you
place two vectors side-by-side with no dot or cross, you create a tensor.

T⃗ = ⃗u⃗v

11
In component form (matrix form):

Tij = ui vj —#

For example, the ”Momentum Flux Tensor” is ρ⃗v⃗v .

2. Tensor Operations (The ”Grammar”)

Your instructor loves asking about identities (Pages 16-17 of notes). Here are the
essential ones.
A. The Unit Tensor (δ or I): ⃗ Also called the Kronecker Delta or Idemfactor.
(
1 if i = j
δ ij = —#
0 if i ̸= j

B. Tensor-Vector Dot Product (τ · ⃗v ): This results in a Vector. (Order

2 − 1 = 1). Physically, this often represents ”Force per area” on a surface.
X
[τ · ⃗n]i = τij nj — #
j

C. Double Dot Product (τ : ∇⃗v ): This results in a Scalar. (Order 2 −

2 = 0). This usually appears in the Energy Equation representing ”Viscous
Dissipation” (friction turning into heat).
XX
A:B= Aij Bji — #
i j

3. The Molecular Stress Tensor (Π)

This is the specific tensor used in the Momentum Balance derivation. It combines
two effects.
Definition: Forces acting on a surface element dS are due to:

1. Pressure (P ): Always acts Normal to the surface, compressive.

2. Viscous Stress (τ ): Caused by velocity gradients (friction).

The Combined Formula:

Π = Pδ + τ —#

In index notation:
Πij = P δij + τij

4. Generalized Newton’s Law of Viscosity

This is the most important derivation in this section (Pages 12-13 of your notes).
∂vi
We need a general formula relating Stress (τij ) to Velocity Gradients ( ∂xj
).
Logic Chain:

12
 1. Stress must be linear in velocity gradients (Newtonian Assumption).

2. Stress should vanish if velocity is constant (or pure rotation).

3. Stress must be symmetric (τij = τji ) to prevent infinite angular acceleration.

The Rate of Deformation Tensor (∆): We define a symmetric tensor

describing how fluid stretches/shears:

∂vi ∂vj
∆ij = +
∂xj ∂xi

The General Law: For an Isotropic Newtonian fluid, the stress is:
   
∂vi ∂vj 2
τij = −µ + + µ − κ (∇ · ⃗v )δij
∂xj ∂xi 3

• µ: Dynamic Viscosity (Shear viscosity).

• κ: Bulk Viscosity (usually ≈ 0 for dense gases/liquids).

Vector Form (Memorize This):

 
 T
 2
τ = −µ ∇⃗v + (∇⃗v ) + µ − κ (∇ · ⃗v )δ —#
3

Simplification for Incompressible Fluids: If the fluid is incompressible

(ρ = const), then ∇ · ⃗v = 0. The second term vanishes.

τ = −µ ∇⃗v + (∇⃗v )T
 
—#

This is the form you substitute into the Cauchy Equation to get Navier-Stokes.

Part 3 Supplement: Deep Tensor Algebra

Exam Tip: Professor-Proofing
The ”Index” Trick: If you get stuck on a vector identity, switch to Index
Notation (Einstein Summation). Rules:

• Repeated index (e.g., ai bi ) means Summation ( ).

P

• Free index (e.g., cj ) means it survives in the result.

• δij swaps indices: ui δij = uj .

13
1. The Dyadic Product (The ”Open” Product)
Definition: A dyadic product of two vectors ⃗u and ⃗v forms a Second-Order Tensor
T = ⃗u⃗v . Unlike the dot product (⃗u · ⃗v = scalar), the dyadic product places them
side-by-side.
Matrix Representation:
   
ux   ux vx ux vy ux vz
⃗u⃗v = uy  vx vy vz = uy vx uy vy uy vz  — #
uz uz vx uz vy uz vz

Index Notation:
[⃗u⃗v ]ij = ui vj —#
Physical Meaning in Fluids: The Momentum Flux Tensor is Φ = ρ⃗v⃗v .

• Component Φxy = ρvx vy .

• Meaning: Momentum directed in x (vx ) being transported across a face

normal to y (carried by velocity vy ).

2. The Double Dot Product (A : B)

This is often the most confusing operation. It is essentially a ”Scalar Product for
Tensors.”
Definition:
X3 X 3
A:B= Aij Bji — #
i=1 j=1

Note the order of indices in B! Aij matches with Bji . This is a trace operation.
Example 1: Viscous Dissipation Term (τ : ∇⃗v ) Let’s calculate the double
dot of the Stress Tensor with the Velocity Gradient Tensor.

• τ ij is the stress.
∂v
• [∇⃗v ]ij = ∂xji (Note: some texts define it as ∂vi /∂xj , check your notes care-
fully. Your notes imply standard definition).
X X  ∂vi 
τ : ∇⃗v = τij
i j
∂xj

This results in a single Scalar number representing Energy (Heat generated per
unit volume).
Example 2: Interaction with Unit Tensor (δ : A)
XX
δ:A= δij Aji
i j

14
Since δij = 0 unless i = j:
X X
= δii Aii = (1)Aii = A11 + A22 + A33 = Trace(A)
i i

3. Divergence of a Tensor (∇ · τ )
This appears in the Navier-Stokes derivation. The divergence of a Tensor reduces
it to a Vector.
Math Setup: Let τ be the tensor. ∇ acts as a vector operator.
3
X ∂
[∇ · τ ]j = τij —#
i=1
∂xi

Interpretation: We sum the derivatives along the rows.

Example: Divergence of the Unit Tensor (∇ · (P δ))
X ∂
[∇ · (P δ)]j = (P δij )
i
∂xi

Since δij is constant and only exists when i = j:

∂P
=
∂xj

In vector form:
∇ · (P δ) = ∇P —#
This is why the pressure term in Navier-Stokes is just ∇P .

4. The Convective Term Identity

This is a high-level proof often asked to separate the top students. Show that
∇ · (ρ⃗v⃗v ) = ⃗v (∇ · ρ⃗v ) + ρ(⃗v · ∇)⃗v .
Step 1: Index Notation We want the j-th component of ∇ · (ρ⃗v⃗v ).
X ∂
[∇ · (ρ⃗v⃗v )]j = (ρvi vj )
i
∂xi

Step 2: Product Rule Differentiate the product (ρvi ) · (vj ).

X  ∂(ρvi ) ∂vj

= vj + (ρvi )
i
∂xi ∂xi

Step 3: Reassemble Vectors Look at the two terms:

• Term 1: vj i ∂(ρv i)
P
∂xi
. The sum is exactly ∇ · (ρ⃗v ).

15
 ∂v
• Term 2: ρ
P
i vi ∂xji . The sum is exactly the convective operator (⃗v · ∇) acting
on vj .

Result:
∇ · (ρ⃗v⃗v ) = ⃗v [∇ · (ρ⃗v )] + ρ(⃗v · ∇)⃗v —#
Significance:

• The first term ⃗v [∇·(ρ⃗v )] vanishes if we use the Continuity Equation (∇·ρ⃗v =
−∂ρ/∂t).

• The second term is the standard ”Convective Acceleration” seen in Navier-

Stokes.

5. Vector-Tensor Cross Product (⃗u × τ )

(Reference: Pages 16-18 of Notes)
Your notes define this using the Levi-Civita permutation symbol ϵijk .
X
[⃗u × τ ]il = ϵijk uj τkl
j,k

Note: This operation usually results in a Tensor.

Why is this in your notes? It is used when defining Vorticity or angular
momentum conservation. If the stress tensor τ is symmetric (τij = τji ), then
certain cross-product operations vanish. The symmetry of the stress tensor is
a fundamental requirement for conservation of angular momentum in non-polar
fluids.

Part 4: Conservation Laws & Continuity

1. The General Conservation Law
(Reference: Page 25 of Notes)
Physical Setup: Consider a general physical quantity ϕ (e.g., mass, energy,
charge).

• ρϕ : Volume density of the quantity (amount per unit volume).

• J⃗ϕ : Flux density (amount flowing through unit area per unit time).

The Integral Balance: For a fixed control volume V , the rate of change of
total amount inside must equal the net flux entering across the surface S.
ZZZ
d
ρϕ dV = −S J⃗ϕ · n̂dS
dt V

16
 ⃗ · n̂dS = ⃗ ):
RRR
Using the Divergence Theorem (A ∇ · AdV
ZZZ  
∂ρϕ ⃗
+ ∇ · Jϕ dV = 0
V ∂t

Since this holds for any arbitrary volume, the integrand must be zero.
The Differential Form:
∂ρϕ
+ ∇ · J⃗ϕ = 0 — #
∂t
This is the master equation.

• For Mass: ρϕ = ρ, Flux J⃗ = ρ⃗v . Result: Continuity Equation.

• For Heat: ρϕ = ρCp T , Flux J⃗ = ⃗q. Result: Heat Equation.

2. The Substantial (Convective) Derivative

(Reference: Page 27 of Notes)
Concept: If you stand still in a river and measure temperature, you measure
the Local Derivative (∂T /∂t). If you drift on a boat with the water and measure
temperature, you measure the Substantial Derivative (DT /Dt). The difference
is that on the boat, you are moving to new locations where the water might be
hotter or colder.
Mathematical Definition: For a scalar field ψ(x, y, z, t):

∂ψ ∂ψ ∂ψ ∂ψ
dψ = dt + dx + dy + dz
∂t ∂x ∂y ∂z

Divide by dt, noting that dx/dt = vx , etc.:

 
dψ ∂ψ ∂ψ ∂ψ ∂ψ
= + vx + vy + vz
dt ∂t ∂x ∂y ∂z

Vector Form (Memorize This):

Dψ ∂ψ
= + (⃗v · ∇)ψ —#
Dt ∂t
• ∂ψ
∂t
: Local term (Unsteady term).

• (⃗v · ∇)ψ: Convective term (Change due to spatial gradients).

17
 Applying to a Vector (Acceleration): The acceleration ⃗a is the substantial
derivative of velocity ⃗v :
D⃗v ∂⃗v
⃗a = = + (⃗v · ∇)⃗v —#
Dt ∂t
Warning: The term (⃗v · ∇)⃗v is non-linear and makes fluid dynamics difficult.

3. Euler’s Equation (Inviscid Flow)

(Reference: Page 59/Concept of ”No Viscosity”)
Physical Setup: Consider a fluid where viscous forces are negligible (η = 0).
This is called an Ideal Fluid or Inviscid Flow. Starting from Navier-Stokes:
D⃗v
ρ = −∇P + 
µ∇2
⃗v + ρ⃗g
Dt
Euler’s Equation:
 
∂⃗v
ρ + ⃗v · ∇⃗v = −∇P + ρ⃗g —#
∂t

Bernoulli’s Equation (Derivation from Euler): If flow is Steady, In-

compressible, and Irrotational (∇ × ⃗v = 0):

1. Use identity (⃗v · ∇)⃗v = ∇( 12 v 2 ) − ⃗v × (∇ × ⃗v ).

2. Since irrotational, ⃗v × (∇ × ⃗v ) = 0.

3. Euler becomes: ρ∇( 12 v 2 ) = −∇P − ρ∇(gh).

 
4. Rearrange: ∇ 21 v 2 + Pρ + gh = 0.

This implies the quantity in brackets is constant everywhere:

1 2
ρv + P + ρgh = Constant — #
2

Part 5: Kinetic Theory of Viscosity (Gases)

Exam Tip: Professor-Proofing
The ”Pheel”: Viscosity in gases is caused by Momentum Transport.
Molecules from a fast-moving layer jump into a slow-moving layer, speeding
it up (and vice versa). This exchange of momentum feels like a drag force
(friction).

18
1. The Physical Setup
Consider a gas flowing in the x-direction with a velocity gradient along the y-axis
(dvx /dy).

• System: Ideal Gas of particle mass m and number density n (particles/m3 ).

• Molecular Motion: Random thermal velocity ⟨v⟩.

• Mean Free Path (λ): The average distance a particle travels before col-
liding.

• Concept: Particles crossing a plane at y = 0 come from an average vertical

distance of y = ±λ cos θ.

2. The Derivation (Flux Integration)

(Reference: Pages 28-30 of Notes)
Step 1: Particle Flux (dN ) Consider an elemental solid angle dΩ = sin θdθdϕ.
The number of particles hitting unit area per unit time from direction (θ, ϕ) is:
 
n⟨v⟩
dN = cos θ dΩ
4π

(The cos θ accounts for the projected area).

Step 2: Momentum Transport A particle coming from above (y = +λ cos θ)
carries x-momentum:
 
dvx
px (+) = m vx (0) − λ cos θ
dy

A particle coming from below (y = −λ cos θ) carries x-momentum:

 
dvx
px (−) = m vx (0) + λ cos θ
dy

Step 3: Net Momentum Flux (Shear Stress) The Shear Stress τyx is the
net momentum transfer per unit area per time.
Z
τyx = (px (+) − px (−))dN
hemisphere

Since momentum flows from fast to slow, the net transfer is proportional to the
difference:
dvx
∆px = −2mλ cos θ
dy
Integrating over the upper hemisphere (θ : 0 → π/2, ϕ : 0 → 2π):
Z 2π Z π/2   
dvx n⟨v⟩
τyx = dϕ −2mλ cos θ cos θ sin θ dθ
0 0 dy 4π

19
 Step 4: Solving the Integral The integral of the angular part:
π/2 π/2
cos3 θ
Z 
2 1
cos θ sin θdθ = − =
0 3 0 3

The 2π from dϕ cancels part of the 4π.

1 dvx
τyx = − nmλ⟨v⟩ —#
3 dy

3. The Result: Coefficient of Viscosity

Comparing this with Newton’s Law of Viscosity (τyx = −η dv
dy
x
), we get:
1
η = nmλ⟨v⟩ — #
3

4. Temperature and Pressure Dependence

This is a common conceptual question. We know from Kinetic Theory:
√
• Mean velocity: ⟨v⟩ ∝ T .

• Mean free path: λ = √ 1

2πd2 n
(where d is diameter).

• Density n appears in the numerator of η and denominator of λ, so they

cancel.

Substituting these into the viscosity formula:

√ √
 
1
η ∝n· · T =⇒ η ∝ T
n

Crucial√Conclusions: 1. Gases: Viscosity INCREASES with Tempera-

ture (η ∝ T ). (More collisions = more momentum transfer). 2. Pressure:
Viscosity is INDEPENDENT of Pressure (for ideal gases). (Higher P means
more molecules, but shorter path λ; they cancel out).

5. Thermal Conductivity (Fourier’s Law)

Your notes (Page 31) perform a similar derivation for Heat Transport. Instead of
transferring momentum (mvx ), particles transfer Energy (E = 32 kB T ).
The resulting Thermal Conductivity κ is:
1
κ = ncv λ⟨v⟩ — #
3
where cv is the specific heat capacity per particle.

20
Part 9: Exam Strategy & Master Formula Sheet
1. The ”Professor-Proof” Strategy
Your instructor penalizes students who ”jump to the math.” To maximize your
score, follow this **4-Step Protocol** for every long problem:

1. The ”Setup” Paragraph (Crucial): Before writing a single equation,

write 3 sentences defining the physical system. Example: ”Consider a steady,
laminar flow of an incompressible Newtonian fluid through a horizontal circu-
lar pipe. We utilize a cylindrical coordinate system (r, θ, z) where flow aligns
with the z-axis.”

2. List Assumptions Explicitly: Do not just assume them. List them as

bullet points:

• ρ = constant (Incompressible)
• ∂/∂t = 0 (Steady State)
• vθ = 0, vr = 0 (Unidirectional flow)

3. Define Control Volume: ”We select a cylindrical shell of length L, radius

r, and thickness ∆r.”

4. State Boundary Conditions: Write them mathematically AND in En-

glish. ”At r = R, vz = 0 (No-Slip Condition).”

2. Master Formula Sheet

A. Conservation Equations (Vector Form)
Continuity Equation (Mass):

∂ρ
+ ∇ · (ρ⃗v ) = 0 — #
∂t
For Incompressible fluid (ρ = const):

∇ · ⃗v = 0 — #

Navier-Stokes Equation (Momentum): (For Incompressible, Newtonian

fluid)  
∂⃗v
ρ + ⃗v · ∇⃗v = −∇P + µ∇2⃗v + ρ⃗g — #
∂t

21
 Creeping Flow Equation (Re ≪ 1):

0 = −∇P + µ∇2⃗v + ρ⃗g —#

B. Constitutive Laws
Newton’s Law of Viscosity (Tensor):
 
T 2
τ = −µ[∇⃗v + (∇⃗v ) ] + µ − κ (∇ · ⃗v )δ —#
3

Fourier’s Law of Heat Conduction:

⃗q = −k∇T —#

Fick’s Law of Diffusion:

J⃗A = −DAB ∇cA —#

C. Shell Balance Results (High Yield)

1. Falling Film (Inclined Plane):
ρgδ 2 cos β
• Velocity: vz (x) = 2η
[1 − (x/δ)2 ]
W ρgδ 3 cos β
• Flow Rate: Q = 3η

2. Pipe Flow (Hagen-Poiseuille):

∆P R2
• Velocity: vz (r) = 4ηL
[1 − (r/R)2 ]
π∆P R4 ρ
• Mass Flow: w = 8ηL

• Average Velocity: ⟨v⟩ = 12 vmax

D. Mathematical Identities
Substantial Derivative:
D ∂
= + ⃗v · ∇ — #
Dt ∂t
Stress Tensor Relation:

Π = Pδ + τ —#

22
E. Kinetic Theory (Gases)
Viscosity: √
1 2 mkB T
η = nmλ⟨v⟩ = —#
3 3π d2
√
(Remember: η ∝ T for gases).

Part 10: Essential Visuals for the Exam

Exam Tip: Professor-Proofing
Drawing Rule: Never draw a ”naked” diagram. Every arrow must have
a label (vz , τrz ). Every axis must have a coordinate (r, z). An unlabeled
diagram is free marks for your instructor to deduct.

1. The Falling Film (Inclined Plane)

Search Term: "falling film velocity profile derivation"
How to Draw It:

1. Draw a surface inclined at angle β to the vertical (or horizontal, matches

your derivation).

2. Draw the liquid layer of thickness δ.

3. Coordinate Axes: Draw z pointing down the slope, x pointing outward

from the surface.

4. The Shell: Draw a thin rectangle inside the fluid parallel to the slope. Label
its thickness ∆x.

5. Velocity Profile: Draw arrows starting small at the wall (x = δ) and

getting largest at the free surface (x = 0). Connect the tips with a Parabolic
Curve.

6. Forces: Draw a gravity vector ⃗g pointing straight down, and resolve it into
g cos β (driving flow) and g sin β (normal).

2. Flow in a Circular Tube (Hagen-Poiseuille)

Search Term: "hagen poiseuille flow velocity profile"
How to Draw It:

1. Draw a long horizontal cylinder (Pipe).

2. The Shell: Draw a concentric inner cylinder (a ring). Label radius r and
thickness ∆r.

23
 3. Velocity Profile: Draw a centerline (r = 0). Draw velocity vectors pointing
right.

4. Shape: The vector at r = 0 is vmax . The vector at wall r = R is 0. Connect

them with a parabola.

5. Forces: Label P0 pushing on the left face and PL pushing on the right face.

3. The Stress Tensor Cube

(Reference: Page 11-13 of Notes) Search Term: "fluid stress tensor cube
components"
How to Draw It:

• Draw a 3D Cube.

• On the Right Face (Normal to x-axis):

– Draw arrow τxx (Normal Stress) pointing out.

– Draw arrow τxy (Shear) pointing up (y-direction).
– Draw arrow τxz (Shear) pointing back (z-direction).

• Why this matters: This diagram justifies the definition of the tensor Π
used in the Navier-Stokes derivation.

4. Viscosity vs Temperature Graph

(Reference: Page 9 of Notes) Search Term: "viscosity vs temperature liquids
vs gases"
How to Draw It:

• Axes: Y-axis = Viscosity η, X-axis = Temperature T .

• Curve 1 (Liquids): Starts high, decays exponentially downwards. (Liquids

get thinner as they heat up).

• Curve 2 (Gases): Starts low, increases linearly/square-root upwards. (Gases

get thicker as they heat up).
√
• Label: ηgas ∝ T .

5. The Differential Control Volume

(For Navier-Stokes Derivation) Search Term: "Eulerian control volume mass
balance"
How to Draw It:

• Draw a small static cube in space dimensions ∆x, ∆y, ∆z.

• Inlet Face (Left): Arrow labeled ρvx |x .

24
 • Outlet Face (Right): Arrow labeled ρvx |x+∆x .

• Repeat for Y and Z directions.

• Center: Label Density ρ.

Part 11: Comprehensive ”End-Sem” Problems

Exam Tip: Professor-Proofing

Warning: Your notes (Page 37) suddenly switch from Fluids to Diffusion
(Semiconductors). This implies your course is actually Transport Phe-
nomena. Do not neglect the diffusion problem below; it is a likely candidate
since it appeared in an assignment.

Problem 1: Flow of Two Immiscible Fluids

(Reference: Pages 51–53 of Notes)
The Question: Consider two immiscible, incompressible fluids flowing be-
tween two infinite parallel stationary plates located at y = b and y = −b.

• Region I (0 < y < b): Viscosity ηI , Density ρI .

• Region II (−b < y < 0): Viscosity ηII , Density ρII .

The flow is driven by a constant pressure gradient dP dx

= PLL−P0 . Derive the velocity
profiles vxI (y) and vxII (y) and the condition for the interface velocity.

Step 1: Physical Setup & Assumptions

Setup:

• Coordinates: Cartesian (x, y, z). Flow is along x. Plates are parallel to the
xz-plane.

• Interface: The plane y = 0 separates the two fluids.

Assumptions:

• Steady state (∂/∂t = 0).

• Laminar, fully developed flow (vx depends only on y, vy = 0, vz = 0).

• Gravity effects are neglected (or included in pressure).

25
Step 2: The Equation of Motion
Start with the Navier-Stokes equation (x-component) or Shell Balance. Since flow
is fully developed, inertial terms are zero.

∂P d2 v x
0=− +η 2
∂x dy
Rearranging for each region:
d2 vx 1 dP
=
dy 2 η dx
dP
Let the pressure gradient be constant: dx
= − ∆P
L
(where ∆P > 0 favors flow).

d2 vx ∆P
2
=− —#
dy ηL

Step 3: General Solutions

Integrate twice with respect to y:
∆P 2
vx (y) = − y + C1 y + C 2
2ηL

We must apply this to both regions separately. Region I (y > 0):

∆P 2
vxI (y) = − y + C1I y + C2I
2ηI L

Region II (y < 0):

∆P 2
vxII (y) = − y + C1II y + C2II
2ηII L

Step 4: Boundary Conditions (The Crux)

We need 4 conditions to solve for 4 constants (C1I , C2I , C1II , C2II ).
BC 1 (No-Slip at Top Wall): At y = b, vxI = 0.

∆P 2
0=− b + C1I b + C2I
2ηI L

BC 2 (No-Slip at Bottom Wall): At y = −b, vxII = 0.

∆P 2
0=− b − C1II b + C2II
2ηII L

BC 3 (Velocity Continuity at Interface): At y = 0, vxI = vxII .

C2I = C2II

26
(Let’s call this interface velocity v0 ). BC 4 (Shear Stress Continuity at In-
terface): At y = 0, τyx I II
= τyx . From Newton’s Law: τyx = −η dv
dy
x
.

dvxI dvxII
−ηI = −ηII
dy y=0 dy y=0

Differentiating our velocity equations:

ηI I
ηI (C1I ) = ηII (C1II ) =⇒ C1II = C
ηII 1

Step 5: The Final Solution

Solving this system of algebra (as done in your notes Page 52): Interface Velocity
(v0 ):
∆P b2 ηI + ηII
 
v0 = × . . . (Algebraic simplification)
4L ηI ηII
(Note: For the exam, unless asked to simplify fully, setting up the 4 equations
and 4 constants clearly is worth 90% of the marks. Ensure you explicitly state BC
3 and BC 4.)

27
Problem 2: Diffusion with Recombination
(Reference: Assignment 2, Page 37 of Notes)
The Question: Consider an n-type semiconductor. The transport of minor-
ity charge carriers (holes, pn ) is governed by diffusion and recombination. The
governing equation is:
∂pn pn − pn0 ∂ 2 pn
=− + Dp 2
∂t τp ∂x
Derive the steady-state concentration profile ∆pn (x) given the boundary condi-
tions:

• ∆pn (0) = ∆pn (0) (Injection at surface).

• ∆pn (∞) = 0 (Decays far away).

Step 1: Physical Setup

Assumptions:

• Steady State: ∂pn /∂t = 0.

• 1D Flow: Variation only along x.

• Excess Carriers: Define ∆pn (x) = pn (x) − pn0 .

Step 2: Simplifying the Equation

Substitute ∆pn into the governing equation. Note that ∂ 2 (pn0 )/∂x2 = 0 since
equilibrium concentration is constant.

∆pn d2 ∆pn
0=− + Dp
τp dx2

Rearrange:
d2 ∆pn 1
2
= ∆pn
dx Dp τp

Step 3: Defining Diffusion Length

The term Dp τp has dimensions of length squared. We define the Diffusion Length
Lp : p
Lp = Dp τp
The equation becomes:

d2 ∆pn 1
2
− 2 ∆pn = 0 — #
dx Lp

28
Step 4: General Solution
This is a second-order linear differential equation.

∆pn (x) = Aex/Lp + Be−x/Lp

Step 5: Applying Boundary Conditions

BC 2 (x → ∞): The concentration cannot grow to infinity. Therefore, the
coefficient of the growing exponential (ex/Lp ) must be zero.

A=0

Equation reduces to: ∆pn (x) = Be−x/Lp .

BC 1 (x = 0): At the surface, the excess concentration is fixed at ∆pn (0).

∆pn (0) = Be0 =⇒ B = ∆pn (0)

Step 6: Final Result

∆pn (x) = ∆pn (0)e−x/Lp ∆pn (x) = ∆pn (0)e−x/Lp ∆pn (x) = ∆pn (0)e−x/Lp ∆pn (x) = ∆pn (0)e−x/Lp
Physical Interpretation: The excess carriers decay exponentially as they diffuse
into the material. The characteristic distance they travel before recombining is the
Diffusion Length Lp .

Part 12: The Conceptual ”Pheel”

1. Viscosity: It’s Not Friction, It’s Transport
The Intuition: We usually think of viscosity as ”fluid friction” or ”thickness.”
But in Transport Phenomena, this is imprecise. The Real Physics: Viscosity is
Momentum Transport.

• Imagine two trains running parallel at different speeds.

• Passengers on the fast train throw heavy sandbags to the slow train. The
slow train gains momentum (speeds up).

• Passengers on the slow train throw bags to the fast train. The fast train
loses momentum (slows down).

• Result: A ”drag force” appears between the trains, trying to equalize their
speeds.

• In Fluids: The ”sandbags” are molecules moving randomly due to thermal

energy. This is why Gases get more viscous as they get hotter (sandbags are

29
 thrown faster/harder), whereas Liquids get less viscous (molecules slide over
each other easier).

2. The ”Convective” Derivative (⃗v · ∇⃗v )

The Intuition: Why does a fluid particle accelerate even in ”Steady State”
(∂⃗v /∂t = 0)? The Real Physics: Imagine a river flowing through a nozzle
(narrowing pipe).

• You stand on the bank and watch point A. The velocity at point A is constant
5 m/s today, tomorrow, and forever. (∂v/∂t = 0).

• Now, drop a leaf in the water. As the leaf moves from the wide part to the
narrow part, it must speed up to conserve mass.

• The leaf accelerates not because time changed, but because its position
changed to a region of higher velocity.

• This ”acceleration due to moving to a new neighborhood” is the Convective

Derivative.

3. Why Parabolic? (Pipe Flow)

The Question: Why is the velocity profile in a pipe a parabola? Why not a
triangle or a flat line? The Real Physics: It is a balance of forces.

1. Pressure Force: Pushes uniformly on the whole cross-section. Since Force

∝ r2 , the cumulative push grows rapidly with radius.

2. Shear Stress (τ ): To balance this uniform pressure push, the shear stress
must increase linearly from the center to the wall (τ ∝ r).

3. Velocity (vz ): Newton’s Law says τ ∝ −dv/dr. If the slope of veloc-

ity (dv/dr) must be linear (∝ r), then the velocity function itself must be
quadratic (∝ r2 ).

Key Phrase for Exam: ”The linear distribution of shear stress required to balance
a constant pressure gradient imposes a parabolic velocity profile for Newtonian
fluids.”

4. Pressure vs. Viscous Stress

The Intuition: Both terms appear in the Navier-Stokes equation. What’s the
fundamental difference? The Real Physics:

• Pressure (P ): Acts Normal to the surface. It is associated with Reversible

Work. You can compress a spring (fluid) and get the energy back. It creates
potential energy (P/ρ).

30
 • Viscous Stress (τ ): Acts Tangential (mostly). It is associated with Irre-
versible Entropy Generation. The energy used to overcome viscosity is
lost forever as Heat.

• Consequence: This is why Bernoulli’s Equation (Conservation of Energy)

only works when Viscosity is zero. Once you add viscosity, you bleed energy.

5. The No-Slip Condition

The Intuition: Why does fluid velocity equal zero at a solid wall? The Real
Physics: It is due to molecular interactions between the fluid and the solid lattice.

• At the microscopic scale, the solid wall looks like a jagged mountain range.

• Fluid molecules get trapped in the valleys of the solid surface.

• They collide billions of times per second, thermalizing with the wall.

• Macroscopically, this looks like the fluid layer is ”stuck” to the wall.

• Note: This fails only in rarefied gases (high vacuum) where molecules act
like bouncy balls (Slip Condition).

Part 13: HOTS & Curveball Questions

1. The Annulus Trap (Boundary Conditions)
The Question: In Hagen-Poiseuille flow (circular pipe), we set the integration
constant C1 = 0 because shear stress cannot be infinite at r = 0. Consider flow
in an Annulus (the space between two coaxial cylinders, radii κR and R). Does
C1 = 0 here? Why or why not?
The Answer: NO, C1 ̸= 0.

• Reasoning: The fluid exists only in the region κR < r < R. The center
of the pipe (r = 0) is outside the fluid domain (it is inside the inner solid
cylinder).

• Since the fluid never occupies r = 0, there is no mathematical singularity to

worry about.

• Consequence: The shear stress distribution is:

∆P C1
τrz = r+ —#
2L r

• Velocity: You must apply the No-Slip condition at two walls (r = κR and
r = R) to find C1 and C2 . This results in a much more complex velocity
profile with a logarithmic term.

31
2. Couette vs. Poiseuille (Driving Forces)
The Question: Consider flow between two flat parallel plates. Case A: Both
plates stationary, pressure gradient applied. Case B: Top plate moves, zero pressure
gradient. How do the velocity profiles differ?
The Answer:

• Case A (Plane Poiseuille): Driven by Pressure (∆P ).

dτ
= const =⇒ τ ∝ y =⇒ vx ∝ y 2
dy
Result: Parabolic Profile.

• Case B (Plane Couette): Driven by Wall Motion (Shear).

dτ dvx
= 0 =⇒ τ = const =⇒ = const
dy dy

Result: Linear Profile (vx = V hy ).

Exam Tip: If both exist (moving plate + pressure drop), the solution is the linear
sum of both (Parabola + Line).

3. The ”Fake” Gravity

The Question: In the Navier-Stokes equation, we write −∇P + ρ⃗g . Can we
combine these terms? What is ”Dynamic Pressure”?
The Answer: Yes. Gravity is a conservative field: ⃗g = −∇(gz) (assuming z
is up). We can define a Modified Pressure P:

P = P + ρgz —#

The N-S equation becomes:

D⃗v
ρ = −∇P + µ∇2⃗v —#
Dt
Significance: This proves that flow is driven by the gradient of total potential
(Pressure + Gravity), not just pressure. Fluid won’t flow simply because pressure
is higher at the bottom of a tank (hydrostatic balance). It flows only if ∇P ̸= 0.

4. Non-Newtonian Fluids (The Constitutive Curveball)

The Question: If the fluid is Non-Newtonian (e.g., blood or toothpaste), which
equations remain valid and which fail?
The Answer:

1. Continuity Equation: VALID. (Mass is always conserved).

32
 2. Cauchy Momentum Equation: VALID. (ρ⃗a = ∇ · Π + ρ⃗g comes from
Newton’s 2nd Law, which applies to all matter).

3. Navier-Stokes Equation: INVALID.

4. Why? Navier-Stokes assumes τ ∝ ∇v. For Non-Newtonian fluids, the

relationship (Constitutive Law) is different (e.g., Power Law τ ∝ |∇v|n ).
You must plug the new stress law into the Cauchy equation.

5. The Incompressible Paradox

The Question: Does ∇ · ⃗v = 0 mean the density ρ is constant everywhere?
The Answer: No.

• ∇·⃗v = 0 means the Substantial Derivative of density is zero (Dρ/Dt = 0).

• This means the density of a specific fluid particle does not change as it moves.

• Example: Oil flowing next to water. Both are incompressible (∇ · ⃗v = 0).

However, if you look at the whole pipe, ρ changes from 1000 to 800 across
the interface. ρ is not uniform, but the flow is incompressible.

Part 14: Deep Dive - The ”Hidden” Physics

1. Why is the Stress Tensor Symmetric? (τij = τji )
The Question: Why must the shear stress on the x-face pointing in y (τxy ) equal
the shear stress on the y-face pointing in x (τyx )?
The Deep Explanation: This is a requirement of the **Conservation of
Angular Momentum**. Consider a tiny cubic fluid element of side length L.

• Torque (Γ): Generated by shear stresses on opposite faces.

Γ ∼ (Force) × (Arm) ∼ (τ L2 ) × L = τ L3

• Moment of Inertia (I): Mass × Radius squared.

I ∼ (ρL3 ) × L2 = ρL5

• Angular Acceleration (α): α = Γ/I.

τ L3 τ
α∼ 5
=
ρL ρL2

33
The Consequence: As the cube shrinks to a point (L → 0), the angular ac-
celeration α → ∞ unless the net torque vanishes. For the torque to vanish, the
opposing shear stresses must balance exactly:

τxy = τyx —#

2. The Reynolds Number (Re): The Ratio of Truth

The Question: What does Re = ρvD
η
actually mean? Why does it predict tur-
bulence?
The Deep Explanation: It is the ratio of Inertial Forces to Viscous
Forces.

• Inertial Force: The force required to stop a moving fluid packet (∼ ρv 2 L2 ).

• Viscous Force: The friction force trying to smooth out the flow (∼ ηvL).
ρv 2 L2
• Ratio: ηvL
= ρvL
η
= Re.

Significance:

• Low Re (< 2000): Viscosity wins. Friction dampens out any small distur-
bances instantly. Flow is Laminar.

• High Re (> 4000): Inertia wins. If a particle wiggles, its momentum

carries it further off course before friction can stop it. Disturbances grow.
Flow becomes Turbulent.

3. Kinematic Viscosity (ν): Momentum Diffusivity

The Question: Why do we have two viscosities, η (Dynamic) and ν = η/ρ
(Kinematic)?
The Deep Explanation: Compare the equations for Heat and Momentum:

• Heat Equation: ∂T
∂t
= α∇2 T (where α is Thermal Diffusivity).

• Navier-Stokes (Simplified): ∂⃗v

∂t
= ν∇2⃗v (where ν = η/ρ).

Conclusion: ν is the Momentum Diffusivity. It tells you how fast a change in

velocity (momentum) spreads through the fluid.

• Water: ν is low. If you stir it, the swirl stays local for a while.

• Honey: ν is high. If you stir one part, the whole pot moves almost instantly
because momentum diffuses rapidly.

34
4. The Stokes Assumption (2/3µ − κ)
The Question: In the Generalized Newton’s Law, we saw a term involving κ
(Bulk Viscosity). Why do we usually ignore it?
The Deep Explanation: The term is 32 µ − κ (∇ · ⃗v ).

1. For Incompressible fluids, ∇ · ⃗v = 0. The term vanishes mathematically.

2. For Monatomic Gases, Kinetic Theory predicts that Bulk Viscosity κ = 0.

3. Stokes’ Hypothesis: Stokes assumed that the thermodynamic pressure

equals the mechanical pressure (Pmech = P̄ ). This forces κ = 0.

4. Reality: κ ̸= 0 only for dense gases with complex molecules (polyatomic)

absorbing vibrational energy during rapid compression (sound waves). For
standard fluid dynamics, it is safely ignored.

Part 16: Targeted Quiz Practice

Exam Tip: Professor-Proofing
Analysis of your Quiz Image:

• Q1: Tests assumptions of Navier-Stokes (Constant Viscosity).

• Q2: Tests Divergence of a Scalar × Unit Tensor (∇ · (P δ̂)).

• Q3 & Q4: Test Divergence of Dyadic Products (∇ · (A

⃗ B)).
⃗

• Q5: Tests the Continuity Equation using the Substantial Derivative

(D/Dt).

The problems below cover these exact vector identities and physical concepts.

Topic A: Vector Tensor Identities (Like Q2, Q3, Q4)

Problem 1: Divergence of Pressure Tensor (Matches Quiz Q2) Prove that
∇ · (P δ̂) = ∇P .

• Starting Step: Write the expressionPin Index Notation. Recall that [δ̂]ij =
δij (Kronecker delta) and [∇ · T]j = i ∂i Tij .

Solution:
X ∂
[∇ · (P δ̂)]j = (P δij )
i
∂x i

35
Since δij is constant (0 or 1), it comes out of the derivative:
X ∂P
= δij
i
∂xi

The summation over i disappears because δij is non-zero only when i = j:

∂P
=
∂xj

In vector form, a vector with components ∂P/∂xj is the gradient:

∇ · (P δ̂) = ∇P - - - #

Problem 2: Divergence of a Dyadic Product (Matches Quiz Q4) Derive

the expansion for ∇ · (⃗u⃗v ).

• Starting Step: Use index notation for the divergence of a tensor: ∂i (ui vj ).
Use the Product Rule of calculus.

Solution: We look at the j-th component:

X ∂
[∇ · (⃗u⃗v )]j = (ui vj )
i
∂x i

Apply Product Rule ∂(ab) = (∂a)b + a(∂b):

X  ∂ui  
∂vj

= vj + ui
i
∂xi ∂xi

Separate the sums: !

X ∂ui X ∂vj
= vj + ui
i
∂xi i
∂xi
Translate back to Vector Notation:

• Term 1: ( ∂i ui ) is (∇ · ⃗u). So, (∇ · ⃗u)⃗v .

P

• Term 2:
P
ui ∂i is the operator (⃗u · ∇). So, (⃗u · ∇)⃗v .

Final Result:
∇ · (⃗u⃗v ) = (∇ · ⃗u)⃗v + (⃗u · ∇)⃗v - - - #
⃗ and A
*(Note: This matches Option B in your Quiz Q4, assuming ψ ⃗ are vectors).*

Problem 3: Divergence of Scalar-Dyadic (Matches Quiz Q3) Expand

∇ · (ρ⃗v⃗v ). (This is the Convective Momentum Term).

• Starting Step: Treat ρ⃗v as vector A

⃗ and the second ⃗v as vector B.
⃗ Use the
result from Problem 2.

36
 ⃗ B)
Solution: Using ∇ · (A ⃗ = (∇ · A)
⃗ B⃗ + (A
⃗ · ∇)B
⃗ with A
⃗ = ρ⃗v and B
⃗ = ⃗v :

∇ · (ρ⃗v⃗v ) = [∇ · (ρ⃗v )]⃗v + (ρ⃗v · ∇)⃗v

∇ · (ρ⃗v⃗v ) = ⃗v [∇ · (ρ⃗v )] + ρ(⃗v · ∇)⃗v - - - #

*(This derivation is key to converting Conservation of Momentum into the Sub-
stantial Derivative form).*

Topic B: Navier-Stokes Assumptions (Like Q1)

Problem 4: Deriving Navier-Stokes The Cauchy Equation is ρ D⃗ v
Dt
= −∇P +
∇ · τ + ρ⃗g . Under what condition does ∇ · τ become µ∇2⃗v ?

• Starting Step: Substitute Newton’s Law τij = µ(∂j vi + ∂i vj ). Pull µ out

of the derivative.

Solution:   
∂ ∂vi ∂vj
∇·τ = µ +
∂xj ∂xj ∂xi
If µ is constant (spatially uniform viscosity):
 
∂ ∂vi ∂ ∂vj
=µ +µ
∂xj ∂xj ∂xi ∂xj

∂
= µ∇2 vi + µ (∇ · ⃗v )
∂xi
If the fluid is Incompressible (∇ · ⃗v = 0), the second term vanishes. Conclusion:
Navier-Stokes requires: 1. Constant Viscosity (Matches Quiz Q1 Option A). 2.
Incompressible flow (usually).

Topic C: Continuity Equation Forms (Like Q5)

Problem 5: Converting Continuity to Substantial Derivative Transform
∂ρ
∂t
+ ∇ · (ρ⃗v ) = 0 into the form involving Dρ/Dt.

• Starting Step: Expand the divergence term ∇ · (ρ⃗v ) using the product rule:
⃗v · ∇ρ + ρ∇ · ⃗v .

Solution:
∂ρ
+ ⃗v · ∇ρ + ρ(∇ · ⃗v ) = 0
∂t
Identify the definition of the Substantial Derivative: DρDt
= ∂ρ
∂t
+ ⃗v · ∇ρ. Substitute
this in:
Dρ
+ ρ(∇ · ⃗v ) = 0 - - - #
Dt

37
*(This matches Quiz Q5 Option A).*

Problem 6: Incompressible Continuity Show that for an incompressible

fluid, the substantial derivative of density is zero.

• Starting Step: Incompressible means density of a specific particle doesn’t

change as it moves.
Dρ
Solution: If ρ = constant, then Dt
= 0. Using the result from Problem 5:

0 + ρ(∇ · ⃗v ) = 0 =⇒ ∇ · ⃗v = 0 - - - #

38