FLUID MECHANICS

Chapter 10 — Approximate Solutions of the Navier–

Stokes Equation

COMPLETE STUDY NOTES | Articles 10.1 – 10.6

Topics Covered:
Creeping Flow | Inviscid Flow | Irrotational Flow |
Boundary Layer Flow
 TABLE OF CONTENTS

Section Topic
Introduction — Exact vs. Approximate Solutions
Nondimensionalized Equations of Motion
Creeping Flow (Stokes Flow / Low Re Flow)
Inviscid Regions of Flow — Euler Equation
Irrotational Flow — Potential Flow Analysis
Boundary Layer Approximation (Prandtl)

NOTE: Numerical examples (10-1 through 10-10) are excluded. Focus is entirely on
concepts, theory, and exam-ready understanding.

QUICK REFERENCE — THE FOUR FLOW APPROXIMATIONS

CREEPING INVISCID IRROTATIONAL BOUNDARY

FLOW FLOW FLOW LAYER
Re << 1 (Very Re >> 1 (Very Vorticity = 0 Near wall region
Low) High)
Viscous forces Net viscous Inviscid AND no Thin region, full
dominate forces negligible rotation N-S needed
Article 10.1 — Introduction
What is the Navier–Stokes (N–S) Equation?
The Navier–Stokes equation is the fundamental equation of fluid mechanics that describes
the motion of a viscous, incompressible, Newtonian fluid. It is derived from Newton's
second law applied to a fluid element, accounting for pressure forces, viscous (friction)
forces, gravity, and inertial effects.
Important: The N–S equation is itself a MODEL — it assumes Newtonian fluid behaviour,
constant properties, and incompressibility. It is not exact in an absolute sense, but is an
excellent engineering model.

Exact Solutions vs. Approximate Solutions

Type Definition
Start with the FULL N–S equation. Some terms drop out
naturally due to the geometry or symmetry of the problem
Exact Solution
(e.g., infinite parallel plates). No prior simplification of the
equation.
Certain terms in the N–S equation are eliminated BEFORE
Approximate solving, based on the class of problem. The simplification is
Solution made a priori (beforehand). Different regions may use
different approximations.

KEY IDEA: Very few exact analytical solutions of the N–S equation exist. For the vast
majority of engineering problems, we must use approximations or computers (CFD). This
chapter focuses on approximations.

Different Approximations for Different Regions

In any practical flow (e.g., liquid flowing from one tank to another), different
approximations are valid in different parts of the flow field:
• Fluid Statics: used in a region far from the inlet where fluid is nearly stationary
• Irrotational (Potential) Flow: used near the pipe inlet where viscous effects are
absent
• Boundary Layer: used very close to the pipe wall where velocity gradients are
large
 • Full N–S: required in complex regions like near a pipe exit

DANGER: If you apply the wrong approximation to a region, your mathematics may be
perfectly correct but your answer will be physically WRONG. Always verify that your
approximation is appropriate.

How Do We Know Which Approximation to Use?

By comparing the ORDER OF MAGNITUDE of the various terms in the equations of
motion. If one term is much smaller than another, it can be neglected. This is formalised
using nondimensionalization — covered in Article 10.2.
Article 10.2 — Nondimensionalized Equations of Motion

Why Nondimensionalise?
Nondimensionalization converts the equations of motion into a dimensionless form. This
allows us to compare the relative magnitudes of different physical effects (inertia,
viscosity, gravity, unsteadiness) and identify which terms are negligible.

Scaling Parameters
We introduce characteristic (reference) values to make variables dimensionless:
Scaling Parameter Physical Meaning
L (Characteristic
Typical size of the problem domain
Length)
V (Characteristic
Typical fluid velocity
Speed)
f (Characteristic
Relevant frequency for unsteady flows
Frequency)
(P0 − P∞)
(Reference Pressure Typical pressure difference driving the flow
Difference)
g (Gravitational
Used for normalising gravity term
Acceleration)

The Four Dimensionless Numbers in the N–S Equation

After nondimensionalizing the full N–S equation, four dimensionless groups appear. The
relative importance of each term is governed by these numbers:

Strouhal (St) St = fL/V Importance of Flow is highly

unsteadiness unsteady
Euler (Eu) Eu = Pressure vs. inertia Pressure effects
(P0−P∞)/ρV² forces dominate
 Froude (Fr) Fr = V/√(gL) Inertia vs. gravity Gravity effects
forces negligible
Reynolds Re = ρVL/μ Inertia vs. viscous Viscous effects
(Re) forces negligible

Key Conclusions from Nondimensionalisation

• The continuity equation has NO dimensionless parameters — all terms are equal in
magnitude, so nothing can be dropped from it.
• Nondimensional variables (starred quantities) are all of order ONE if proper
scaling is used. This means we can directly compare the dimensionless groups.
• For STEADY flow: Strouhal number drops out (St = 0). No unsteady term.
• For flows WITHOUT a free surface: Froude number drops out. Gravity only
causes hydrostatic pressure, which is absorbed into Modified Pressure.

Modified Pressure — Removing Gravity from N–S

For flows without a free surface, gravity only adds a hydrostatic pressure distribution on
top of the dynamic pressure. We define:
Modified Pressure (P'): P' = P + ρgz (where z is measured upward from a datum)
Using modified pressure P' instead of actual pressure P, the gravity term disappears from
the N–S equation. After solving, we simply add back the hydrostatic term. This approach
is commonly used in Computational Fluid Dynamics (CFD).

WARNING: Do NOT use modified pressure in flows with free-surface effects (open
channels, waves, etc.) because gravity plays a physical role there.
Article 10.3 — Creeping Flow Approximation

Also Called: Stokes Flow | Low Reynolds Number Flow

Definition
Creeping flow occurs when the Reynolds number is VERY SMALL (Re << 1). The flow
is so slow (or the fluid is so viscous, or the object is so tiny) that VISCOUS FORCES
completely dominate over INERTIAL FORCES.

Physical Conditions for Creeping Flow

Condition Example
Very high viscosity Pouring honey, syrup, molasses
Very small velocity Slow seeping of oil
Very small length Microorganisms (bacteria ~1 micron) swimming in
scale water
Small gaps + high
Lubricating oil in bearing gaps (tens of microns)
viscosity

What Happens to the N–S Equation in Creeping Flow?

Because Re << 1, the inertial (acceleration) terms on the LEFT side of the N–S equation
are orders of magnitude smaller than the viscous terms on the RIGHT. So the ENTIRE
LEFT SIDE IS DROPPED.

Creeping Flow N–S (dimensional form): Pressure gradient = Viscosity × (Laplacian of

velocity) i.e., Grad(P) = μ × Laplacian(V)

This is also called the Stokes Equation. Notice: inertia is gone. The equation is LINEAR
and much easier to solve.
Pressure Scale in Creeping Flow
Flow Type Pressure Scales With
Normal (inertia-
ρV² (dynamic pressure — from Bernoulli)
dominated) flow
Creeping flow
μV/L (viscous pressure scale)
(viscosity-dominated)

Critical Property: Density Disappears!

In creeping flow, density (ρ) does NOT appear in the momentum equation at all. This leads
to an important result:
• Aerodynamic drag on an object in creeping flow does NOT depend on fluid
density
• Drag depends only on: velocity (V), object size (L), and viscosity (μ)
• This is fundamentally different from high-Re flow where drag strongly depends on
ρ

Drag Force in Creeping Flow — General Result

General Creeping Flow Drag: Drag Force = Constant × μ × V × L The constant
depends only on the SHAPE and ORIENTATION of the object. It cannot be found from
dimensional analysis alone.

Drag on a Sphere — The Famous Stokes Law

Stokes Law (Drag on a Sphere in Creeping Flow): Drag Force = 3π × μ × V × D
where D = diameter of the sphere Note: 2/3 of this drag comes from VISCOUS forces,
1/3 from PRESSURE forces.

Terminal Velocity of a Particle

When a small particle falls at constant speed (terminal velocity), its weight downward
equals drag plus buoyancy upward. Using Stokes Law, terminal velocity depends on:
• Particle diameter squared (D²) — smaller particle falls MUCH slower
• Density difference between particle and fluid
• Fluid viscosity
• NOT significantly on fluid density (since density drops from momentum equation)
Practical result: The terminal velocity of a dense, small particle in creeping flow is
essentially INDEPENDENT of fluid density, but highly sensitive to fluid viscosity. A
volcanic ash particle settles at nearly the same speed regardless of altitude (despite air
density varying 10× from 15,000 m altitude to sea level).

Absence of Inertia — Physical Consequences

• No gliding is possible: The instant a microorganism stops moving its tail, it STOPS
moving. There is no coasting.
• Microorganism propulsion: Flagella (corkscrew or whip motion) work because
they exploit viscosity, not inertia. A simple flapping tail (like a dolphin) would not
propel a microorganism.
• Reversibility: Creeping flow is time-reversible. If you reverse the forces, the flow
reverses exactly.
Article 10.4 — Approximation for Inviscid Regions of Flow

Important Clarification — What Does 'Inviscid' Mean?

COMMON MISCONCEPTION: 'Inviscid flow' does NOT mean the fluid has no
viscosity. All engineering fluids have viscosity. 'Inviscid region of flow' means a region
where the NET VISCOUS FORCES are negligible compared to pressure and/or inertial
forces.

A better term (used by some authors) is 'regions of flow with negligible net viscous forces'.
Viscous stresses still exist, but they cancel each other out — leaving no net force on fluid
elements.

When Does an Inviscid Region Exist?

When the Reynolds number Re is VERY LARGE (Re >> 1), the viscous term (1/Re) in
the nondimensional N–S equation becomes negligibly small compared to the other terms.
This is exactly the opposite of creeping flow.
Flow Regime Re Dominant Forces
Creeping (Stokes) Re << 1 Viscous forces dominate
Inviscid region Re >> 1 Inertial and pressure forces dominate

The Euler Equation

In inviscid regions, the viscous term is dropped from the N–S equation. The result is called
the EULER EQUATION:
Euler Equation: ρ × [DV/Dt] = −Grad(P) + ρg (This is the N–S equation without the
viscous term μ × Laplacian(V))

Where is the Euler Equation Valid?

• In the BULK of a high-Re flow, far away from solid walls
• Away from wakes behind bodies
• In the outer flow region (outside boundary layers)
Where is the Euler Equation NOT Valid?
• NEAR solid walls — because the no-slip condition requires large velocity
gradients, making viscous terms important regardless of Re
• In the WAKE of a body — velocity gradients are large there too
• In regions of FLOW SEPARATION

KEY LIMITATION: When we drop the viscous term, we lose the highest-order
derivative of velocity. Mathematically, this means we can NO LONGER impose the no-
slip condition at walls. The Euler equation allows SLIP at walls — which is physically
unrealistic near walls.

Bernoulli Equation from the Euler Equation

For steady, incompressible flow in inviscid regions, the Euler equation can be integrated
along a streamline to yield the Bernoulli equation:
Bernoulli Equation (Inviscid Region — Along a Streamline): P/ρ + V²/2 + gz = C
(constant ALONG each streamline) Note: The constant C may differ from one streamline
to another.

Solid Body Rotation — An Inviscid but Rotational Region

Solid body rotation (like fluid in a spinning container) is an example of a flow that is
INVISCID (no net viscous force) but also ROTATIONAL (the vorticity is not zero). In
this case:
• The Bernoulli equation holds along each streamline
• BUT the Bernoulli constant C is DIFFERENT on each streamline
• Pressure increases parabolically with radius (centripetal effect)
Article 10.5 — The Irrotational Flow Approximation

Definition of Irrotational Flow

A region of flow is irrotational if fluid particles have NO NET ROTATION.
Mathematically, this means the VORTICITY VECTOR is zero (or negligibly small)
everywhere in that region.
Irrotational Condition: Vorticity (ζ) = Curl of Velocity (∇ × V) = 0 (This means no
spinning of fluid elements)

Velocity Potential Function (φ — phi)

From vector mathematics: if the curl of a vector is zero, that vector can be expressed as the
gradient of a scalar function. For irrotational flow:
Velocity = Gradient of Potential Function: V = ∇φ (phi is called the velocity
potential function) This combines 3 unknown velocity components into ONE unknown
scalar function!
This is a huge mathematical simplification — 3 equations reduced to 1.

The Laplace Equation — Heart of Irrotational Flow

Substituting the velocity potential into the continuity equation gives:
Laplace Equation (for velocity potential φ): ∇²φ = 0 (The Laplacian of phi equals
zero) This is valid ONLY in irrotational regions of flow.

The Laplace equation is a well-known, well-studied equation in mathematics and physics

(electrostatics, heat conduction, etc.). Many analytical and numerical solution techniques
exist.

Key Properties of the Laplace Equation

• Mass (density) does NOT appear — solutions are valid for ANY incompressible
fluid
• Time does NOT appear — irrotational flow instantly adjusts to boundary
conditions at any moment
 • It is LINEAR — superposition of solutions is valid
• Solutions are dominated by GEOMETRY (boundary conditions)

Why Does the Viscous Term Vanish in Irrotational Flow?

In the N–S equation, the viscous term equals μ × Laplacian(V). For irrotational flow, V =
∇φ, so:
μ × Laplacian(V) = μ × Laplacian(∇φ) = μ × ∇(Laplacian(φ)) = μ × ∇(0) = 0
The viscous term vanishes NOT because the fluid has no viscosity, but because
irrotationality forces it to zero. All irrotational regions are automatically also inviscid
regions.

Bernoulli Equation for Irrotational Flow

Bernoulli Equation (Irrotational Region): P/ρ + V²/2 + gz = C (constant
EVERYWHERE in the irrotational region) This is more powerful than the inviscid
version — C is the SAME for ALL streamlines!

Inviscid vs. Irrotational — Important Distinction

Property Inviscid Region Irrotational Region
Net viscous
Negligible (Re >> 1) Zero (due to zero vorticity)
forces
Vorticity May or may not be zero Always zero (by definition)
Bernoulli Constant ALONG a
Constant EVERYWHERE
constant C streamline
More restrictive? Less restrictive More restrictive
Euler eq. applies? Yes Yes (from different reason)

Stream Function (ψ — psi)

For TWO-DIMENSIONAL irrotational flow, we can also define a stream function ψ. Both
φ and ψ satisfy the Laplace equation. They are harmonic conjugates:
• Stream function ψ: defined to automatically satisfy CONTINUITY. Laplace
equation comes from irrotationality.
• Velocity potential φ: defined to enforce IRROTATIONALITY. Laplace equation
comes from continuity.
 • Lines of constant ψ = streamlines
• Lines of constant φ = equipotential lines
• Streamlines and equipotential lines are MUTUALLY ORTHOGONAL (intersect
at 90°)

CAUTION — Axisymmetric Flow: For axisymmetric irrotational flow, the Laplace

equation is valid for φ but NOT for ψ. The stream function equation in axisymmetric
flow is a different (but still linear) PDE.

Superposition of Irrotational Flows

Because the Laplace equation is LINEAR and HOMOGENEOUS, any sum of solutions is
also a solution. This allows us to build complex flow fields by adding together simple
'building block' flows.

Elementary Building Block Flows

Building
Description Stream Function ψ
Block
1. Uniform Fluid flowing at constant speed V in
ψ=V×y
Stream x-direction
2. Line Fluid emerging (source) or entering
Source / (sink) uniformly from a line parallel ψ = (V̇/L)/(2π) × θ
Sink to z-axis
3. Line Fluid spinning around an axis —
ψ = −(Γ/2π) × ln(r)
Vortex circular streamlines
Source + Sink of equal strength at
4. Doublet ψ = −K × sinθ/r
same point — used to model bodies

Important Superpositions
• Uniform Stream + Doublet: Models potential flow over a CIRCULAR
CYLINDER. One streamline (c = 0) forms a perfect circle which acts as the
cylinder wall.
• Line Sink + Line Vortex: Models a spiralling drain (water spiralling into a
bathtub drain).
• Uniform Stream + Source + Sink: Can model flow around elongated body shapes
(Rankine ovals).
Potential Flow over a Circular Cylinder
By superposing a uniform stream and a doublet, we get a flow that has a circular streamline.
This circle is treated as the cylinder surface. Key results:
• Flow is symmetric fore-and-aft (front-to-back) and top-to-bottom
• At the top of the cylinder: fluid speed = 2V (twice free stream). Pressure is lowest
here.
• At stagnation points (nose and tail): speed = 0, pressure is highest.
• No-slip condition CANNOT be satisfied — irrotational flow allows slip at walls

D'Alembert's Paradox (1752): With the irrotational flow approximation, the

aerodynamic drag on ANY non-lifting body of ANY shape immersed in a uniform stream
is ZERO. This is because the fore-aft pressure symmetry means front and rear pressures
exactly cancel.
In reality, drag is non-zero because: (1) pressure on the rear of a body never fully recovers
(pressure drag), and (2) the no-slip condition causes viscous drag. Both effects are missed
by irrotational flow.

Fish Eye Location — Real-World Application

The zero-pressure point on a body is where the surface pressure equals the ambient free-
stream pressure (Cp = 0). Fish eyes are located near this point so that the eye pressure does
not change as the fish swims faster or slower — preventing eye distortion and blurred
vision.

Two-Region Tornado Model

A horizontal slice of a tornado is modelled with two concentric regions:
• Inner region (r < R): Solid body rotation — inviscid but ROTATIONAL.
Pressure increases parabolically outward.
• Outer region (r > R): Irrotational flow — inviscid AND irrotational. Bernoulli
applies with the SAME constant everywhere. Pressure increases as speed decreases
outward.
The lowest pressure is at the centre (eye) of the tornado. This is why a tornado can be
extremely destructive — the very large pressure difference between the eye and
surroundings can cause buildings to explode outward.
Article 10.6 — The Boundary Layer Approximation

The Problem with the Euler Equation

The Euler equation (inviscid/irrotational approximation) cannot satisfy the no-slip
condition at solid walls. This leads to:
• Zero shear stress (friction) at walls — physically wrong
• Zero aerodynamic drag (D'Alembert's paradox) — physically wrong
• No flow separation prediction — physically wrong for bluff bodies

Prandtl's Breakthrough — 1904

Ludwig Prandtl (1875–1953): Introduced the boundary layer concept in 1904. He
divided the flow into two regions: 1. OUTER FLOW — treated as inviscid and/or
irrotational (Euler/potential flow applies) 2. BOUNDARY LAYER — a thin region near
the wall where viscous effects and vorticity cannot be ignored This bridged the gap
between the Euler equation (slip at wall) and the N–S equation (no-slip at wall).

What is a Boundary Layer?

A boundary layer is the very thin region of fluid immediately adjacent to a solid surface.
Within this region:
• The fluid velocity transitions from zero at the wall (no-slip) to the free-stream
velocity some distance away
• Velocity gradients normal to the wall are LARGE
• Viscous forces are NOT negligible, even at high Re
• The flow is ROTATIONAL (vorticity is generated at the wall)

Boundary Layer Thickness (δ — delta)

By convention, boundary layer thickness δ is defined as the distance from the wall at which
the velocity component parallel to the wall reaches 99% of the free-stream value.
• δ is NOT constant — it grows with downstream distance x
• δ is much smaller than the body length L (that is the whole point — the layer is
THIN)
 • The curve δ(x) is NOT a streamline — streamlines cross it

Reynolds Number and Boundary Layer Thickness

Key Result: At a given streamwise location x: δ/x ~ 1/√(Re_x) where Re_x = ρVx/μ
The HIGHER the Reynolds number, the THINNER the boundary layer at that location.

Laminar vs. Turbulent Boundary Layer

Property Laminar Turbulent
Fuller profile — more
Velocity Smooth, parabolic-like
uniform away from wall,
profile (Blasius profile)
steep near wall
Shear stress at
Lower Higher (more mixing)
wall
δ thickness Thinner Thicker
Resistance to Less resistant — separates More resistant — stays
separation earlier attached longer
Re_x range
Re_x < ~10⁵ Re_x > ~3 × 10⁶
(flat plate)

Transition from Laminar to Turbulent

• Critical Re (smooth flat plate): Re_x ≈ 10⁵ — disturbances begin to grow
• Transition Re: Re_x ≈ 3 × 10⁶ — fully turbulent boundary layer
• Engineering Critical Re: Re_x ≈ 5 × 10⁵ — used in practice because real surfaces
are not smooth
• Trip wires: Rough wires placed on the surface to force EARLY transition to
turbulence at a desired location (useful in experiments and design)

Key Assumptions of the Boundary Layer Approximation

• Boundary layer is VERY THIN: δ << x (or δ << L)
• The outer flow is unaffected by the boundary layer (thin layer has negligible effect)
• Pressure in the y-direction (normal to wall) is approximately CONSTANT through
the boundary layer
 • The velocity component normal to the wall (v) is MUCH SMALLER than the
component parallel to the wall (u)

Most Important Result — Pressure Across Boundary Layer

Pressure Gradient Through Boundary Layer (y-direction): ∂P/∂y ≈ 0 (pressure
does NOT change across the boundary layer) Consequence: The pressure at the wall =
The pressure at the edge of the boundary layer Practical use: Static pressure taps in the
wall give the same pressure as the outer flow!

The Boundary Layer Equations (Laminar, 2D, Steady, Incompressible)

After order-of-magnitude analysis, the N–S equation in the boundary layer reduces to TWO
equations:
1. Continuity: ∂u/∂x + ∂v/∂y = 0 (same as incompressible continuity)
2. x-Momentum: u(∂u/∂x) + v(∂u/∂y) = U(dU/dx) + ν(∂²u/∂y²) [Inertia on left] =
[Pressure term from outer flow] + [Viscous term]

The pressure term is replaced using the Bernoulli equation applied to the outer flow: dP/dx
= −ρU(dU/dx)

Nature of Boundary Layer Equations — Parabolic

Property Full N–S Equation Boundary Layer Equation
Mathematical
Elliptic Parabolic
type
Information In ALL directions (upstream
Only DOWNSTREAM
propagation + downstream)
Boundary
Only 3 sides (upstream, top,
conditions On ALL sides of domain
bottom — NOT downstream)
needed
Solution March downstream — solve
Iterative (coupled)
method step by step

Boundary Layer Procedure (Step by Step)

3. Solve the OUTER FLOW first (Euler or potential flow) to get U(x) and P(x)
4. ASSUME the boundary layer is very thin — it does not disturb the outer flow
 5. Solve the boundary layer equations with boundary conditions: u = v = 0 at wall, u
→ U(x) far from wall, known starting profile at x = x_start
6. Calculate quantities of interest: δ(x), wall shear stress, total skin friction drag
7. VERIFY that δ << x — if not, the approximation is invalid

Flat Plate Boundary Layer — Blasius Solution

For a uniform stream flowing parallel to a flat plate, there is no pressure gradient (dU/dx =
0, called 'zero pressure gradient' boundary layer). The boundary layer equations were first
solved numerically by Heinrich Blasius (1908, Prandtl's PhD student).
• The solution uses a similarity variable — combining x and y into a single
dimensionless variable η (eta)
• The Blasius velocity profile is unique — it applies at ANY x-location when plotted
using the similarity variable
• Result: δ grows proportional to √x along the plate
• Excellent agreement with experimental data

Limitations of the Boundary Layer Approximation

• Re too low: If Re < ~1000, δ is not thin enough (δ/L > 3%), approximation is poor
• High wall curvature: If the radius of curvature R ≈ δ, centripetal effects matter
and ∂P/∂y ≠ 0
• Turbulent flow: The laminar BL equations (above) do not apply. Turbulent BL
equations are much more complex
• Flow separation: When adverse pressure gradient is too strong, the boundary
layer separates from the wall. The parabolic nature of BL equations breaks down
— REVERSE FLOW cannot be handled

Boundary Layer on Surfaces Other Than Flat Plates

The boundary layer concept applies not just to flat walls but also to:
• Jets: Thin shear region around a high-speed jet mixing with surrounding fluid
• Wakes: Region of disturbed flow behind a body
• Mixing Layers: Between two streams of different speed
In all cases, the boundary layer is thin, the Reynolds number is high, and the same equations
apply (with appropriate coordinate systems).
 CONCEPTUAL QUESTIONS & MODEL ANSWERS
These questions are designed to test deep understanding. Read both the question AND the
full answer. Each answer explains the CONCEPT — not just a formula.

Q1. What is the difference between an 'exact solution' and an 'approximate

solution' of the Navier–Stokes equation?
★ Answer: Exact vs. Approximate
✦ An EXACT solution starts with the FULL N–S equation. Terms may drop out
naturally due to the specific geometry or symmetry (e.g., fully developed channel flow),
but no terms are eliminated beforehand.
✦ An APPROXIMATE solution eliminates one or more terms BEFORE solving,
based on the class of problem. The simplification is chosen depending on which physical
effects are dominant.
✦ Example: In creeping flow, the inertial terms are dropped a priori because Re <<
1. In boundary layer theory, the y-momentum equation and second-order x-derivative of
u are dropped because the layer is thin.
✦ Risk: If the approximation is wrong from the start, the mathematics may be correct
but the answer is physically meaningless.

Q2. What does 'inviscid flow' actually mean? Does it mean the fluid has no
viscosity?
★ Answer: Inviscid = Negligible NET Viscous Forces
✦ NO — 'inviscid flow' does NOT mean the fluid has no viscosity. All engineering
fluids have viscosity.
✦ 'Inviscid region of flow' means a region where NET viscous forces on fluid
elements are negligible compared to pressure and inertial forces.
✦ Viscous stresses still exist in the fluid, but they cancel each other out, leaving no
significant NET force.
✦ This occurs in HIGH Reynolds number regions far from walls and wakes.
✦ The preferred term is 'inviscid region of flow' rather than 'inviscid flow' or
'frictionless flow', both of which are misleading.
Q3. Why does density disappear from the momentum equation in creeping flow?
What is the physical significance?
★ Answer: Density Vanishes in Creeping Flow
✦ In creeping flow (Re << 1), the inertial terms on the left side of the N–S equation
are much smaller than viscous terms. The ENTIRE inertial side is dropped.
✦ Since density appears ONLY in the inertial terms (ρ × DV/Dt), when these terms
are eliminated, density disappears from the equation.
✦ Physical significance: Aerodynamic drag on an object in creeping flow is
INDEPENDENT of fluid density. It depends only on velocity, object size, and viscosity.
✦ Practical example: A volcanic ash particle settling at 50,000 ft altitude settles at
essentially the same speed as at sea level — even though air density increases 10× —
because in creeping flow, density doesn't matter.

Q4. What is the velocity potential function φ (phi)? Why is it so useful?

★ Answer: Velocity Potential Function
✦ In an irrotational region of flow, the velocity vector can be written as the
GRADIENT of a scalar function φ (phi), called the velocity potential function: V = ∇φ
✦ This is possible because: if vorticity = 0, then curl(V) = 0, and any vector whose
curl is zero equals the gradient of some scalar.
✦ Usefulness: Instead of solving for 3 unknown velocity components (u, v, w), we
solve for ONE unknown scalar φ. This reduces the problem enormously.
✦ Substituting V = ∇φ into the continuity equation gives the Laplace equation (∇²φ =
0) — a well-studied linear PDE with many known solution techniques.
✦ Once φ is known, all velocity components follow by differentiation.

Q5. Why is the Bernoulli constant the SAME everywhere in an irrotational region,
but only constant ALONG a streamline in an inviscid region?
★ Answer: Bernoulli Constant — Irrotational vs. Inviscid
✦ In an INVISCID region (not necessarily irrotational): Bernoulli is derived by
integrating the Euler equation ALONG a streamline. The integral holds only along that
path — C can vary from streamline to streamline.
✦ In an IRROTATIONAL region: When vorticity = 0, the cross product term V × ζ
in the Euler equation vanishes. The remaining equation says the gradient of (P/ρ + V²/2 +
gz) = 0 everywhere. If a gradient is zero everywhere, the quantity itself must be constant
everywhere — not just along streamlines.
✦ Therefore, irrotational flow is a STRICTER condition than inviscid flow. All
irrotational regions are also inviscid, but not all inviscid regions are irrotational (example:
solid body rotation).

Q6. What is D'Alembert's Paradox? Why does it occur, and what does it tell us
about irrotational flow?
★ Answer: D'Alembert's Paradox
✦ D'Alembert's Paradox (1752) states: with the irrotational flow approximation, the
aerodynamic drag on ANY non-lifting body of any shape in a uniform stream is ZERO.
✦ Why? Because irrotational flow predicts fore-aft symmetry of pressure. The
pressure at the front stagnation point exactly matches the pressure at the rear — so
pressure forces cancel. Additionally, since no-slip cannot be applied, there is no viscous
drag.
✦ In reality, drag is non-zero because: (1) Real pressure on the rear of a body never
recovers fully (separation occurs), causing pressure drag. (2) No-slip condition causes
viscous (friction) drag.
✦ What this tells us: Irrotational flow is excellent for predicting the FRONT of a
body (near the nose), but fails badly at the REAR where separation and wake effects
dominate.

Q7. What is Prandtl's boundary layer concept? Why was it a major breakthrough?
★ Answer: Prandtl's Boundary Layer (1904)
✦ Before Prandtl, fluid mechanics had two unsatisfactory extremes: the Navier–
Stokes equation (correct but unsolvable for complex geometry) and the Euler equation
(solvable but gives zero drag — D'Alembert's paradox).
✦ Prandtl's idea: Divide the flow into TWO regions. In the OUTER FLOW (away
from walls), treat the flow as inviscid/irrotational — Euler/potential flow applies. Near
walls, insert a thin BOUNDARY LAYER where viscous effects are fully included.
✦ This bridged the Euler equation (slip allowed) and N–S equation (no-slip
enforced) through the boundary layer — enforcing no-slip at the wall while simplifying
the equations in the outer region.
✦ The breakthrough: For the first time, drag could be predicted, the no-slip condition
could be applied, and separation could be understood — all with equations that were
actually solvable before computers.

Q8. Why is the pressure in the y-direction (normal to wall) approximately constant
across a boundary layer? What is the practical consequence?
★ Answer: Zero Pressure Gradient Across Boundary Layer
✦ The boundary layer is VERY THIN (δ << L). This means streamlines within the
boundary layer have negligible curvature when viewed at the scale of δ.
✦ Curved streamlines require centripetal acceleration, which comes from a pressure
gradient across the radius of curvature. Since streamlines are nearly straight (not curved)
in a thin BL, there is no centripetal force needed, so ∂P/∂y ≈ 0.
✦ Consequence: The pressure at any y-location across the boundary layer equals the
pressure at the outer edge of the boundary layer (which is obtained from the outer flow
analysis using Bernoulli).
✦ Practical use: Engineers measure the static pressure through tiny holes (pressure
taps) drilled into the wall surface. Because ∂P/∂y = 0, this wall pressure equals the outer
flow pressure — so wall pressure taps directly give the outer inviscid flow pressure!

Q9. What is superposition in irrotational flow? Give examples of what can be built.
★ Answer: Superposition of Irrotational Flows
✦ Because the Laplace equation is LINEAR, any sum of two or more solutions is
also a solution. This means we can ADD simple building block flows to create complex
ones.
✦ The velocity vectors also add directly: V_total = V_1 + V_2 at every point.
✦ Examples of superposition: (1) Uniform stream + Doublet = Flow over a circular
cylinder. (2) Line Sink + Line Vortex = Spiralling drain flow. (3) Multiple sources/sinks
= Complex body shapes (Rankine ovals). (4) Method of Images: A source + its mirror
image across a wall = Flow with a solid wall.
✦ Important limit: Superposition only works for IRROTATIONAL flows. You
cannot add a turbulent jet (rotational) with an irrotational outer flow.

Q10. When does a boundary layer separate from the wall? What is the
consequence?
★ Answer: Boundary Layer Separation
✦ Flow separation occurs when the pressure INCREASES in the flow direction
(adverse pressure gradient — dP/dx > 0, or equivalently dU/dx < 0). The fluid near the
wall, which already has low momentum, cannot push against this rising pressure.
✦ The slow near-wall fluid decelerates, stops, and then reverses direction. The
boundary layer detaches from the wall — this is separation.
✦ Consequence for boundary layer equations: The parabolic BL equations require
flow information to march DOWNSTREAM only. Reverse flow (going upstream) breaks
this parabolic nature. The equations become invalid once separation occurs.
 ✦ Engineering consequence: After separation, a large recirculating wake forms
behind the body. Pressure on the rear of the body drops significantly — causing large
PRESSURE DRAG. This is why bluff bodies (cylinders, flat plates perpendicular to
flow) have much higher drag than streamlined bodies.

Q11. Compare and contrast creeping flow and boundary layer flow — when does
each apply?
★ Answer: Creeping vs. Boundary Layer — Comparison
✦ CREEPING FLOW: Applies when Re << 1 (entire flow). Viscous forces dominate
everywhere. Inertial terms are ZERO. Density drops from momentum equation. Example:
bacteria swimming, honey flow, bearing lubrication.
✦ BOUNDARY LAYER: Applies near a wall when Re >> 1 globally. ONLY in the
thin near-wall region are viscous forces important. In the outer flow, inertial forces
dominate. Density does NOT drop from equations.
✦ Key contrast: In creeping flow, viscosity rules EVERYWHERE. In boundary
layer flow, viscosity only matters in the thin near-wall strip — the rest of the flow is
effectively inviscid.
✦ Both are approximations of the full N–S equation — but applied in completely
opposite Reynolds number regimes.

Q12. Why can't the boundary layer approximation handle flow separation?
★ Answer: BL Equations and Flow Separation
✦ The boundary layer equations are PARABOLIC — meaning information travels
only in the downstream direction. This allows us to solve the equations by marching
forward in x.
✦ When separation occurs, there is REVERSE FLOW near the wall — fluid moving
UPSTREAM. This means information must travel backwards — which is impossible for
a parabolic equation.
✦ Mathematically, the parabolic character is lost and the boundary layer equations
break down at and beyond the separation point.
✦ Physical consequence: We cannot use the BL approximation in separated flow
regions. The full N–S equation (or CFD) must be used there.
 MASTER SUMMARY — DECISION FLOWCHART FOR FLOW
APPROXIMATIONS

Use this flowchart thinking to choose the correct approximation in an exam:

STEP 1: What is the Reynolds number?

Re << 1 (Very Small) Re >> 1 (Very Large)

→ CREEPING FLOW → INVISCID or IRROTATIONAL
• Drop all inertial terms • Stokes FLOW (far from walls)
equation applies • Drag = C × μ × V × • Drop viscous term → Euler equation
L • Density is irrelevant in momentum • Use Bernoulli (along streamline) •
Near walls: need Boundary Layer

STEP 2 (for Re >> 1): Is the flow also irrotational?

Irrotational (Vorticity = 0) Rotational (but still inviscid)

→ POTENTIAL FLOW → EULER EQUATION
• Use velocity potential φ • Solve • Use Euler equation • Bernoulli only
Laplace equation ∇²φ = 0 • Use along streamlines • Example: solid
superposition of building blocks • body rotation • C varies from
Bernoulli constant is SAME streamline to streamline
everywhere

STEP 3: Are we near a solid wall in a high-Re flow?

YES — Near the wall NO — Far from wall

→ BOUNDARY LAYER → OUTER FLOW REGION
• Use BL equations (thin layer) • No- • Use Euler or Potential Flow •
slip enforced at wall • Pressure ∂P/∂y Bernoulli equation applies • No
= 0 across BL • March downstream to viscous effects • Superposition may be
solve • Check if laminar or turbulent used