Skip to main content

Hagen-Poiseuille flow in a circular tube

Short answer

For laminar pipe flow, velocity u = (Δp/(4μL))(R^2 - r^2). Maximum u_max = Δp R^2/(4μL). Average V = u_max/2 = Δp R^2/(8μL). Discharge Q = πΔp R^4/(8μL) (Hagen-Poiseuille). Wall shear stress τ_w = Δp R/(2L). Head loss h_f = 32μVL/(ρgD^2). Friction factor f = 64/Re.

Derive the velocity distribution for steady laminar flow of a Newtonian fluid in a circular pipe of radius RR. Hence obtain the expressions for maximum velocity, average velocity, discharge (Hagen-Poiseuille equation), wall shear stress, head loss and the friction factor in terms of Reynolds number.

Assumptions

Steady, fully developed, laminar, incompressible flow of a Newtonian fluid in a horizontal circular pipe of radius RR and length LL. Velocity is axial only and depends on rr only.

        r
        ^   _____________________________
   R    |  /      ->  ->   ->
        | /   ->  ->  ->   ->   -> u(r)
   0 ---+--------------------------------->  x
        | \   ->  ->  ->   ->
        |  \_______________________________
           p1                       p2

Force balance on a cylinder of radius r

Pressure force =(p1−p2)πr2=Δp πr2= (p_1 - p_2)\pi r^2 = \Delta p\,\pi r^2. Shear force on the cylindrical surface =τ 2πrL= \tau\,2\pi r L. For equilibrium:

τ=Δp r2L\tau = \frac{\Delta p\, r}{2L}

This is linear in rr: zero at the axis, maximum at the wall.

Velocity distribution

With τ=−μ du/dr\tau = -\mu\,du/dr (velocity decreases as rr increases):

dudr=−Δp2μLr\frac{du}{dr} = -\frac{\Delta p}{2\mu L}r

Integrating, with u=0u = 0 at r=Rr = R:

u=Δp4μL(R2−r2)u = \frac{\Delta p}{4\mu L}\left(R^2 - r^2\right)

The profile is a paraboloid.

Maximum velocity (at r=0r = 0)

umax=ΔpR24μLu_{max} = \frac{\Delta p R^2}{4\mu L}

Discharge

Q=∫0Ru 2πr dr=πΔp2μL[R42−R44]=πΔpR48μL=πΔpD4128μLQ = \int_0^R u\,2\pi r\,dr = \frac{\pi \Delta p}{2\mu L}\left[\frac{R^4}{2} - \frac{R^4}{4}\right] = \frac{\pi \Delta p R^4}{8\mu L} = \frac{\pi \Delta p D^4}{128\mu L}

This is the Hagen-Poiseuille equation.

Average velocity

Vˉ=QπR2=ΔpR28μL=umax2\bar V = \frac{Q}{\pi R^2} = \frac{\Delta p R^2}{8\mu L} = \frac{u_{max}}{2}

Wall shear stress

τw=ΔpR2L=ΔpD4L=8μVˉD\tau_w = \frac{\Delta p R}{2L} = \frac{\Delta p D}{4L} = \frac{8\mu\bar V}{D}

Head loss

hf=Δpρg=32μVˉLρgD2h_f = \frac{\Delta p}{\rho g} = \frac{32\mu \bar V L}{\rho g D^2}

Friction factor

Equating with the Darcy-Weisbach equation hf=fLVˉ22gDh_f = \dfrac{f L \bar V^2}{2 g D}:

f=64μρVˉD=64Ref = \frac{64\mu}{\rho \bar V D} = \frac{64}{Re}

The result is valid only for laminar flow, Re<2000Re < 2000. The head loss is directly proportional to the mean velocity and independent of the pipe roughness.

✓ Answer checked by Claude Opus 5.5 (high effort)

All Fluid Mechanics chapter 6 questions →