Chapter 2

Fluid Mechanics for Chemical Engineers

Viscosity, flow regimes, and pumping systems.

Fluid mechanics is fundamental to chemical engineering as it governs how fluids (liquids and gases) move through process equipment. Understanding fluid behavior is essential for designing pumps, pipes, reactors, and separation equipment.

Fluid Properties

Viscosity

Viscosity (μ\mu) measures a fluid's resistance to flow:

  • Dynamic viscosity: Resistance to shear stress
  • Kinematic viscosity: ν=μ/ρ\nu = \mu/\rho

Density and Specific Gravity

  • Density (ρ\rho): Mass per unit volume
  • Specific gravity: Ratio of fluid density to water density at 4°C

Flow Regimes

Reynolds Number

The Reynolds number determines flow regime:

Re=ρvDμ=vDνRe = \frac{\rho v D}{\mu} = \frac{v D}{\nu}

Where:

  • vv = average velocity
  • DD = characteristic length (pipe diameter)
  • μ\mu = dynamic viscosity
  • ν\nu = kinematic viscosity

Flow Classification

  • Laminar flow: Re<2300Re < 2300 (smooth, orderly)
  • Transitional flow: 2300<Re<40002300 < Re < 4000 (unstable)
  • Turbulent flow: Re>4000Re > 4000 (chaotic, mixed)

Bernoulli Equation

For incompressible, frictionless flow:

P1+12ρv12+ρgz1=P2+12ρv22+ρgz2P_1 + \frac{1}{2}\rho v_1^2 + \rho g z_1 = P_2 + \frac{1}{2}\rho v_2^2 + \rho g z_2

Where:

  • PP = pressure
  • vv = velocity
  • zz = elevation

Pressure Drop Calculations

Friction Factor

For laminar flow (Re<2300Re < 2300):

f=64Ref = \frac{64}{Re}

For turbulent flow, use the Colebrook equation:

1f=2.0log(ϵ/D3.7+2.51Ref)\frac{1}{\sqrt{f}} = -2.0 \log\left(\frac{\epsilon/D}{3.7} + \frac{2.51}{Re\sqrt{f}}\right)

Darcy-Weisbach Equation

Pressure drop in pipes:

ΔP=fLDρv22\Delta P = f \frac{L}{D} \frac{\rho v^2}{2}

Where:

  • LL = pipe length
  • DD = pipe diameter
  • ϵ\epsilon = pipe roughness

Pumping Systems

Pump Types

  • Centrifugal pumps: Most common in chemical plants
  • Positive displacement pumps: For high-viscosity fluids
  • Diaphragm pumps: For corrosive or abrasive fluids

Pump Performance

  • Head: Energy per unit weight
  • Capacity: Flow rate
  • Efficiency: Power output / Power input

Net Positive Suction Head (NPSH)

Prevents cavitation:

NPSHavailable=PsuctionPvaporNPSH_{available} = P_{suction} - P_{vapor}

Advanced Fluid Mechanics Concepts

Navier-Stokes Equations

The fundamental equations governing fluid motion are the Navier-Stokes equations:

Conservation of Mass (Continuity):

ρt+(ρv)=0\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0

For incompressible flow, this simplifies to:

v=0\nabla \cdot \mathbf{v} = 0

Conservation of Momentum:

ρ(vt+vv)=p+μ2v+ρg\rho \left( \frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v} \right) = -\nabla p + \mu \nabla^2 \mathbf{v} + \rho \mathbf{g}

Where:

  • ρ\rho = fluid density
  • v\mathbf{v} = velocity vector
  • pp = pressure
  • μ\mu = dynamic viscosity
  • g\mathbf{g} = gravitational acceleration

Boundary Layer Theory

The boundary layer is the thin region near a solid surface where viscous effects are significant:

Boundary Layer Thickness (δ\delta): For laminar flow over a flat plate:

δx=5.0Rex\frac{\delta}{x} = \frac{5.0}{\sqrt{Re_x}}

Where Rex=ρvxμRe_x = \frac{\rho v_\infty x}{\mu} is the local Reynolds number.

Displacement Thickness (δ\delta^*): Represents the distance by which the external streamlines are displaced due to the boundary layer:

δ=0(1uU)dy\delta^* = \int_0^\infty \left(1 - \frac{u}{U}\right) dy

Turbulence Modeling

For turbulent flow, we use time-averaged equations (Reynolds-averaged Navier-Stokes):

Reynolds Decomposition:

u=uˉ+uu = \bar{u} + u'

Where uˉ\bar{u} is the mean velocity and uu' is the fluctuating component.

Turbulent Kinetic Energy (kk):

k=12(u2+v2+w2)k = \frac{1}{2} \left( \overline{u'^2} + \overline{v'^2} + \overline{w'^2} \right)

Non-Newtonian Fluids

Many chemical process fluids exhibit non-Newtonian behavior:

Power Law Model:

τ=K(dudy)n\tau = K \left( \frac{du}{dy} \right)^n

Where:

  • KK = consistency index
  • nn = flow behavior index
  • n<1n < 1: Pseudoplastic (shear thinning)
  • n>1n > 1: Dilatant (shear thickening)
  • n=1n = 1: Newtonian

Bingham Plastic Model:

τ=τ0+μpdudy\tau = \tau_0 + \mu_p \frac{du}{dy}

Where τ0\tau_0 is the yield stress.

Real-World Application: Pipeline Design

Designing a chemical pipeline requires careful fluid mechanics calculations:

Design Considerations

  1. Flow rate requirements
  2. Pressure limitations
  3. Material compatibility
  4. Energy efficiency
  5. Safety factors

Advanced Pipeline Analysis

For complex pipeline systems, consider:

Minor Losses:

hminor=KLv22gh_{minor} = K_L \frac{v^2}{2g}

Where KLK_L is the loss coefficient for fittings, valves, and expansions/contractions.

Network Analysis: For branched or looped systems, use the Hardy-Cross method or equivalent pipe length methods.

Example: Crude Oil Pipeline

Calculate pressure drop for crude oil flowing through a pipeline:

# Pipeline parameters
pipe_diameter = 0.3  # m
pipe_length = 5000   # m
pipe_roughness = 0.000045  # m (commercial steel)
flow_rate = 0.1      # m³/s

# Fluid properties (crude oil)
density = 850        # kg/m³
viscosity = 0.005    # Pa·s

# TODO: Calculate pressure drop
# Steps:
# 1. Calculate velocity from flow rate and area
# 2. Calculate Reynolds number
# 3. Determine friction factor
# 4. Calculate pressure drop using Darcy-Weisbach

velocity = 0
Re = 0
friction_factor = 0
pressure_drop = 0

print(f"Velocity: {velocity:.2f} m/s")
print(f"Reynolds number: {Re:.0f}")
print(f"Friction factor: {friction_factor:.4f}")
print(f"Pressure drop: {pressure_drop:.0f} Pa")

# Additional analysis: Check for turbulent flow transition
if Re < 2300:
    print("Flow is laminar")
elif Re < 4000:
    print("Flow is in transition region")
else:
    print("Flow is turbulent")
    
# Calculate hydraulic power required
power = pressure_drop * flow_rate  # W
print(f"Hydraulic power: {power/1000:.1f} kW")

Your Challenge: Pump Selection Analysis

In this exercise, you'll analyze pump performance and select an appropriate pump for a chemical process.

Goal: Calculate pump power requirements and evaluate system curve.

System Description

A centrifugal pump moves water from a storage tank to a reactor. The system has:

  • Elevation difference: 15 m
  • Pipe length: 200 m
  • Pipe diameter: 0.1 m
  • Pipe roughness: 0.000046 m
  • Required flow rate: 0.02 m³/s
  • Fluid: Water at 20°C (ρ=998\rho = 998 kg/m³, μ=0.001\mu = 0.001 Pa·s)

Pump Performance Data

The pump has the following characteristics:

  • Maximum head: 50 m
  • Maximum flow: 0.05 m³/s
  • Efficiency: 70% at design point
# System parameters
elevation_diff = 15  # m
pipe_length = 200    # m
pipe_diameter = 0.1  # m
pipe_roughness = 0.000046  # m
flow_rate = 0.02     # m³/s

# Fluid properties
density = 998        # kg/m³
viscosity = 0.001    # Pa·s
gravity = 9.81       # m/s²

# TODO: Calculate system head requirement
# System head = elevation head + friction head
velocity = 0
Re = 0
friction_factor = 0
friction_head = 0
system_head = 0

# TODO: Calculate pump power requirement
pump_power = 0

print(f"System head requirement: {system_head:.1f} m")
print(f"Pump power requirement: {pump_power:.1f} kW")

# Check if pump is suitable
if system_head <= 50:
    print("Pump is suitable for this application")
else:
    print("Pump cannot provide sufficient head")

What modifications could reduce the power requirement? How would changing the pipe diameter affect the system?

ELI10 Explanation

Simple analogy for better understanding

Fluid mechanics is the study of how liquids and gases flow through pipes and equipment. Think of it like understanding how water flows through a garden hose. When you turn the faucet just a little, water flows smoothly and quietly - that's laminar flow. When you turn it all the way, water rushes out noisily and chaotically - that's turbulent flow. Chemical engineers need to understand fluid mechanics to design pumps that move chemicals through factories, predict pressure drops in pipes, and ensure processes run efficiently without wasting energy.

Self-Examination

Q1.

What is the difference between laminar and turbulent flow regimes?

Q2.

How do chemical engineers calculate pressure drops in piping systems?

Q3.

Why is viscosity important in chemical process design?