Mechanical Engineer's Handbook

Viscosity and Turbulence

Up until now, fluid dynamics models assumed an ideal fluid with zero friction. Real fluids exhibit internal friction that resists flow, known as viscosity.

Viscosity Definition and Shear Force

Viscosity (η\eta) measures a fluid's resistance to flow.

Consider a fluid sandwiched between two flat plates of area AA separated by a distance LL. The bottom plate is fixed, and a horizontal force FF moves the top plate at constant speed vv:

F=ηvAL    η=FLvAF = \frac{\eta v A}{L} \implies \eta = \frac{F L}{v A}
(1)

The SI unit for viscosity is Pas\text{Pa} \cdot \text{s} or Ns/m2\text{N} \cdot \text{s/m}^2.

Liquid viscosity decreases as temperature increases because molecules gain thermal energy to overcome cohesive forces. Gas viscosity increases with temperature because faster-moving gas molecules collide more frequently across fluid layers.

Poiseuille’s Law and Flow Resistance

In a viscous fluid flowing through a tube, speed is highest at the center and zero at the walls due to drag. The pressure difference p2p1p_2 - p_1 drives flow against resistance RR:

Q=p2p1RQ = \frac{p_2 - p_1}{R}
(2)

For laminar flow of an incompressible fluid through a horizontal tube of radius rr and length ll, the fluid resistance is:

R=8ηlπr4R = \frac{8 \eta l}{\pi r^4}
(3)

Substituting RR yields Poiseuille’s Law:

Q=(p2p1)πr48ηlQ = \frac{(p_2 - p_1) \pi r^4}{8 \eta l}
(4)

Flow rate depends on the fourth power of the radius (r4r^4). Cutting the radius of a tube in half (rr/2r \to r/2) increases resistance by a factor of 24=162^4 = 16, reducing flow rate to 6.25%6.25\% of its original value for a constant pressure difference.

Rearranging Equation 2 shows that flow and resistance cause pressure drops throughout a system:

p2p1=RQp_2 - p_1 = R Q
(5)

Reynolds Number and Turbulence

Turbulence occurs when fluid layers mix via eddies and swirls, caused by high fluid speeds or sharp obstructions. Turbulence greatly increases flow resistance.

The Reynolds number (NRN_R) is a dimensionless parameter used to predict whether flow in a tube of radius rr is laminar or turbulent:

NR=2ρvrηN_R = \frac{2 \rho v r}{\eta}
(6)

where ρ\rho is fluid density, vv is average fluid speed, and η\eta is viscosity.

  • NR<2000N_R < 2000: Flow is laminar.
  • 2000NR30002000 \le N_R \le 3000: Flow is unstable (chaotic) and oscillates between laminar and turbulent states.
  • NR>3000N_R > 3000: Flow is turbulent.

Example 1

An air conditioning duct has a diameter of 18.0 cm18.0 \text{ cm} (r=0.090 mr = 0.090 \text{ m}) and a length of 20.0 m20.0 \text{ m}. It delivers air (ρ=1.23 kg/m3\rho = 1.23 \text{ kg/m}^3, η=1.81105 Pas\eta = 1.81 \cdot 10^{-5} \text{ Pa} \cdot \text{s}) at a volume flow rate of Q=3.84103 m3/sQ = 3.84 \cdot 10^{-3} \text{ m}^3/\text{s}. Determine if the air flow is laminar or turbulent.

First calculate fluid speed v=Qπr2v = \frac{Q}{\pi r^2}, then evaluate the Reynolds number using Equation 6.

Calculate the fluid velocity inside the duct:

v=Qπr2=3.84103 m3/sπ(0.090 m)20.151 m/sv = \frac{Q}{\pi r^2} = \frac{3.84 \cdot 10^{-3} \text{ m}^3/\text{s}}{\pi (0.090 \text{ m})^2} \approx 0.151 \text{ m/s}

Calculate the Reynolds number:

NR=2ρvrη=2(1.23 kg/m3)(0.151 m/s)(0.090 m)1.81105 Pas1840N_R = \frac{2 \rho v r}{\eta} = \frac{2 (1.23 \text{ kg/m}^3)(0.151 \text{ m/s})(0.090 \text{ m})}{1.81 \cdot 10^{-5} \text{ Pa} \cdot \text{s}} \approx 1840

Since NR=1840<2000N_R = 1840 < 2000, the flow is laminar.