Physics:Hagen–Poiseuille flow from the Navier–Stokes equations

From HandWiki

In fluid dynamics, the derivation of the Hagen–Poiseuille flow from the Navier–Stokes equations shows how this flow is an exact solution to the Navier–Stokes equations.[1][2]

Derivation

The laminar flow through a pipe of uniform (circular) cross-section is known as Hagen–Poiseuille flow. The equations governing the Hagen–Poiseuille flow can be derived directly from the Navier–Stokes momentum equations in 3D cylindrical coordinates by making the following set of assumptions:

  1. The flow is steady ( [math]\displaystyle{ \partial(...)/\partial t = 0 }[/math] ).
  2. The radial and swirl components of the fluid velocity are zero ( [math]\displaystyle{ u_r = u_\theta = 0 }[/math] ).
  3. The flow is axisymmetric ( [math]\displaystyle{ \partial(...)/\partial \theta = 0 }[/math] ).
  4. The flow is fully developed ([math]\displaystyle{ \partial u_z/\partial z = 0 }[/math] ). However, this can be proved via mass conservation, and the above assumptions.

Then the angular equation in the momentum equations and the continuity equation are identically satisfied. The first momentum equation reduces to [math]\displaystyle{ \partial p/\partial r = 0 }[/math], i.e., the pressure [math]\displaystyle{ p }[/math] is a function of the axial coordinate [math]\displaystyle{ z }[/math] only. The third momentum equation reduces to:

[math]\displaystyle{ \frac{1}{r}\frac{\partial}{\partial r}\left(r \frac{\partial u_z}{\partial r}\right)= \frac{1}{\mu} \frac{\partial p}{\partial z} }[/math] where [math]\displaystyle{ \mu }[/math] is the dynamic viscosity of the fluid.
The solution is
[math]\displaystyle{ u_z = \frac{1}{4\mu} \frac{\partial p}{\partial z}r^2 + c_1 \ln r + c_2 }[/math]

Since [math]\displaystyle{ u_z }[/math] needs to be finite at [math]\displaystyle{ r = 0 }[/math], [math]\displaystyle{ c_1 = 0 }[/math]. The no slip boundary condition at the pipe wall requires that [math]\displaystyle{ u_z = 0 }[/math] at [math]\displaystyle{ r = R }[/math] (radius of the pipe), which yields

[math]\displaystyle{ c_2 = -\frac{1}{4\mu} \frac{\partial p}{\partial z}R^2. }[/math]

Thus we have finally the following parabolic velocity profile:

[math]\displaystyle{ u_z = -\frac{1}{4\mu} \frac{\partial p}{\partial z} (R^2 - r^2). }[/math]

The maximum velocity occurs at the pipe centerline ([math]\displaystyle{ r=0 }[/math]):

[math]\displaystyle{ {u_z}_{max}=\frac{R^2}{4\mu} \left(-\frac{\partial p}{\partial z}\right). }[/math]

The average velocity can be obtained by integrating over the pipe cross section:

[math]\displaystyle{ {u_z}_\mathrm{avg}=\frac{1}{\pi R^2} \int_0^R u_z \cdot 2\pi r dr = 0.5 {u_z}_\mathrm{max}. }[/math]

The Hagen–Poiseuille equation relates the pressure drop [math]\displaystyle{ \Delta p }[/math] across a circular pipe of length [math]\displaystyle{ L }[/math] to the average flow velocity in the pipe [math]\displaystyle{ {u_z}_\mathrm{avg} }[/math] and other parameters. Assuming that the pressure decreases linearly across the length of the pipe, we have [math]\displaystyle{ - \frac{\partial p}{\partial z} = \frac{\Delta p}{L} }[/math] (constant). Substituting this and the expression for [math]\displaystyle{ {u_z}_\mathrm{max} }[/math] into the expression for [math]\displaystyle{ {u_z}_\mathrm{avg} }[/math], and noting that the pipe diameter [math]\displaystyle{ D = 2R }[/math], we get:

[math]\displaystyle{ {u_z}_{avg} = \frac{D^2}{32 \mu} \frac{\Delta p}{L}. }[/math]

Rearrangement of this gives the Hagen–Poiseuille equation:

[math]\displaystyle{ \Delta p = \frac{32 \mu L ~{u_z}_\mathrm{avg}}{D^2}. }[/math]

References

  1. White, Frank M. (2003). "6". Fluid Mechanics (5 ed.). McGraw-Hill. 
  2. Bird, Stewart, Lightfoot (1960). Transport Phenomena. 

See also