Skip to main content

Derive laminar velocity in inclined pipe

Short answer

For steady gravity-driven laminar flow, the x-momentum reduces to μ (1/r) d/dr(r du/dr) + ρg sinα = 0. Integrating with u(R)=0 yields u(r) = (ρg sinα/(4μ)) (R² - r²), a parabolic profile.

Consider steady, incompressible, laminar flow of a Newtonian fluid in an infinitely long round pipe of diameter D or radius R = D/2 inclined at angle α\alpha as shown in the figure. The fluid flows down the pipe due to gravity alone. Consider the coordinate system shown, with X down the axis of pipe. Derive an expression for the X-component of velocity u as a function of radius r and the other parameters of the problem. [Figure: inclined pipe wall at angle α\alpha to the horizontal, fluid of density ρ\rho and viscosity μ\mu, radial coordinate r, x measured down the pipe axis, gravity g vertical]

Set-up. Pipe of radius RR inclined at angle α\alpha below the horizontal; the fluid flows down the pipe under gravity alone, so there is no externally imposed pressure gradient: ∂p∂x=0\dfrac{\partial p}{\partial x} = 0 (the pressure along the pipe is atmospheric/constant). Flow is steady, laminar, fully developed, so u=u(r)u = u(r) only, and vr=vθ=0v_r = v_\theta = 0.

        x (down the pipe)
     ___ \
    /     \  alpha
   / pipe  \
  /_________\ ____ horizontal
  gravity component along x: g sin(alpha)

Momentum equation along xx (Navier-Stokes in cylindrical coordinates, reduced for this case):

0=ρgsin⁡α+μ 1rddr(rdudr)0 = \rho g\sin\alpha + \mu\,\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right)

(Equivalent force balance on a cylinder of radius rr: weight component ρgsin⁡α πr2L\rho g \sin\alpha\,\pi r^2 L is balanced by the shear force −μdudr 2πrL-\mu\dfrac{du}{dr}\,2\pi r L.)

Integrate.

ddr(rdudr)=−ρgsin⁡αμ r\frac{d}{dr}\left(r\frac{du}{dr}\right) = -\frac{\rho g\sin\alpha}{\mu}\,r rdudr=−ρgsin⁡α2μ r2+C1r\frac{du}{dr} = -\frac{\rho g\sin\alpha}{2\mu}\,r^2 + C_1

C1=0C_1 = 0, since du/drdu/dr must be finite at r=0r = 0. Then

dudr=−ρgsin⁡α2μ r⇒u=−ρgsin⁡α4μ r2+C2\frac{du}{dr} = -\frac{\rho g\sin\alpha}{2\mu}\,r \quad\Rightarrow\quad u = -\frac{\rho g\sin\alpha}{4\mu}\,r^2 + C_2

Boundary condition. No slip at the wall: u=0u = 0 at r=Rr = R, so C2=ρgsin⁡α4μR2C_2 = \dfrac{\rho g\sin\alpha}{4\mu}R^2.

u(r)=ρgsin⁡α4μ(R2−r2)\boxed{u(r) = \frac{\rho g\sin\alpha}{4\mu}\left(R^2 - r^2\right)}

The profile is parabolic with umax=ρgsin⁡α4μR2u_{max} = \dfrac{\rho g\sin\alpha}{4\mu}R^2 at the axis. With R=D/2R = D/2, u=ρgsin⁡α16μ(D2−4r2)u = \dfrac{\rho g\sin\alpha}{16\mu}\left(D^2 - 4r^2\right). The mean velocity is V=umax/2=ρgD2sin⁡α32μV = u_{max}/2 = \dfrac{\rho g D^2\sin\alpha}{32\mu}.

✓ Answer checked by Claude Sonnet 5.5

All Hydraulics chapter 1 questions →