Poiseuille flow

ComplémentA video (in french) about Poiseuille's flow (with some fun science experiments)

La physique animée : Écoulement de Poiseuille d'un fluide visqueux

("Unisciel" and "Culture Sciences Physiques")

Pour lire la vidéo, cliquer ici :

FondamentalVelocity Fields and pressures

In 1835 a French physicist Poiseuille made a series of experiments to determine how a viscous fluid flows in a straight pipe.

His goal was to understand the dynamics of blood flow in humans knowing that the blood plasma behaves like a Newtonian fluid.

A density of incompressible viscous fluid \(\rho\) flows through a cylindrical tube of length \(L\) and radius \(R\).

The pressure at the input of the tube (\(z=0\)) is \(P_E\).

The pressure at the output of the tube (\(z=L\)) is \(P_S\).

We will calculate the velocity and pressure fields within the tube being placed in the steady state.

Laminar flow is assumed :

\(\vec v =v(r)\; \vec u_z\)

It is assumed that \(P=P(r,z)\) (invariance by rotation about the axis Oz) and the gravity is neglected.

Applying the second law of Newton to a system consisting of the fluid contained in the cylinder of radius \(r<R\) at time \(t\) and the mass \(dm_1=dm\) that enters the cylinder between \(t\) and \(t+dt\).

At the moment \(t+dt\), this system is made up of fluid in the cylinder of radius \(r\) at \(t+dt\) and the same mass  \(dm_2=dm\) that exits between \(t\) and \(t+dt\) (incompressible flow).

In the case of steady state, the change in momentum of the system is then simply zero :

\((P_e \pi r^2 \vec u_z - P_s \pi r^2 \vec u_z ) + \left( {\eta \frac{{dv}}{{dr}}2\pi rL\vec u_z } \right) = \vec 0\)

Where :

\(\frac{{dv}}{{dr}} = - \frac{{P_e - P_s }}{{2\eta L}}r\)

By integration, taking account of \(v(R)=0\) :

\(v(r)=\frac {\Delta P}{4 \eta L}(R^2-r^2)\)

This gives the velocity profile in a Poiseuille flow.

The graph above shows the parabolic profile of the velocity field and the variation of the pressure.

Called "charge loss" the amount \(\Delta P=P_E-P_S\).

FondamentalVolumetric flow rate

The volumetric flow rate is here :

\(D_v = \int_0^R {v(r)2\pi rdr} = \int_0^R {\frac{{\Delta P}}{{4\eta L}}(R^2 - r^2 )2\pi rdr}\)

Either :

\(D_v = \frac{{\Delta P}}{{\eta L}}\frac{\pi }{8}R^4\)

We define a "hydraulic resistance", by analogy with an ohmic resistor (\(\Delta U=V_1-V_2=RI\)) :

\(\Delta P=P_E-P_S=R_{hyd}D_v\)

Either :

\(R_{hyd}=\frac {8\eta L}{\pi R^4}\)

The hydraulic resistance is as much higher as the viscosity of the fluid is large and the radius of the tube is small.

ExempleCharge loss and hydraulic resistance

Placing a succession of vertical tubes along the pipe.

The measurement of water levels used to determine the pressure drop :

\(\Delta P=\mu _{eau} g (h_1-h_2)=2.10^3\;Pa\)

If the flow rate is \(D_v=10\;cm^3.s^{-1}=10^{-5}\;m^3.s^{-1}\), the hydraulic resistance can be deduced as :

\(R_{hyd}=\frac {\Delta P}{D_v}=2.10^8\;Pa.s.m^{-3}\)

SimulationJAVA animation by JJ.Rousseau

  • Flow in a tube : click here