Definitions
Nearby Words

# Hagen-Poiseuille equation

The Hagen-Poiseuille equation is a physical law that describes slow viscous incompressible flow through a constant circular cross-section. It is also known as the Hagen-Poiseuille law, Poiseuille law and Poiseuille equation.

## Equation

### Standard fluid dynamics notation

In standard fluid dynamics notation:
$Delta P = frac\left\{8 mu L Q\right\}\left\{ pi r^4\right\}$
or
$Delta P = frac\left\{128 mu L Q\right\}\left\{ pi d^4\right\}$

Where:

$Delta P$ is the pressure drop
L is the length of pipe
$mu$ is the dynamic viscosity
Q is the volumetric flow rate
d is the diameter
$pi$ is the mathematical constant (approximately 3.1416).

### Physics notation

$Phi = frac\left\{dV\right\}\left\{dt\right\} = v pi R^\left\{2\right\} = frac\left\{pi R^\left\{4\right\}\right\}\left\{8 eta\right\} left\left(frac\left\{- Delta P\right\}\left\{Delta x\right\}right\right) = frac\left\{pi R^\left\{4\right\}\right\}\left\{8 eta\right\} frac\left\{ |Delta P|\right\}\left\{L\right\}$

Where:

V is a volume of the liquid poured (cubic meters)
t is the time (seconds)
v is mean fluid velocity along the length of the tube (meters/second)
x is a distance in direction of flow (meters)
R is the internal radius of the tube (meters)
ΔP is the pressure difference between the two ends (pascals)
η is the dynamic fluid viscosity (pascal-second (Pa·s)),
L is the total length of the tube in the x direction (meters).

## Relation to Darcy-Weisbach

This result is also a solution to the phenomenological Darcy-Weisbach equation in the field of hydraulics, given a relationship for the friction factor in terms of the Reynolds number:

$Lambda = \left\{64over \left\{it Re\right\}\right\} ; , quadquad Re = \left\{2rho v rover eta\right\} ; ,$

where Re is the Reynolds number and ρ fluid density. In this form the law approximates the Darcy friction factor, the energy (head) loss factor, friction loss factor or Darcy (friction) factor Λ in the laminar flow at very low velocities in cylindrical tube. The theoretical derivation of a slightly different form of the law was made independently by Wiedman in 1856 and Neumann and E. Hagenbach in 1858 (1859, 1860). Hagenbach was the first who called this law the Poiseuille's law.

The law is also very important specially in hemorheology and hemodynamics, both fields of physiology.

The Poiseuilles' law was later in 1891 extended to turbulent flow by L. R. Wilberforce, based on Hagenbach's work.

## Derivation

The Hagen Poiseuille equation can be derived from the Navier-Stokes equations.

### Viscosity

The derivation of Poiseuille's Law is surprisingly simple, but it requires an understanding of Viscosity. When two layers of liquid in contact with each other move at different speeds, there will be a force between them. This force is proportional to the area of contact A, the velocity difference in the direction of flow Δvxy, and a proportionality constant η and is given by

$F_\left\{text\left\{viscosity, top\right\}\right\} = - eta A frac\left\{Delta v_x\right\}\left\{Delta y\right\}$

The negative sign is in there because we are concerned with the faster moving liquid (top in figure), which is being slowed by the slower liquid (bottom in figure). By Newton's third law of motion, the force on the slower liquid is equal and opposite (no negative sign) to the force on the faster liquid. This equation assumes that the area of contact is so large that we can ignore any effects from the edges and that the fluids behave as Newtonian fluids.

### Liquid flow through a pipe

In a tube we make a basic assumption: the liquid in the center is moving fastest while the liquid touching the walls of the tube is stationary (due to friction). To simplify the situation, let's assume that there are a bunch of circular layers (lamina) of liquid, each having a velocity determined only by their radial distance from the center of the tube.

To figure out the motion of the liquid, we need to know all forces acting on each lamina:

1. The force pushing the liquid through the tube is the change in pressure multiplied by the area: F = -ΔPA. This force is in the direction of the motion of the liquid - the negative sign comes from the conventional way we define $Delta P = P_\left\{end\right\}-P_\left\{top\right\} < 0$.
2. The pull from the faster lamina immediately closer to the center of the tube
3. The drag from the slower lamina immediately closer to the walls of the tube

The first of these forces comes from the definition of pressure. The other two forces require us to modify the equations above that we have for viscosity. In fact, we are not modifying the equations, instead merely plugging in values specific to our problem. Let's focus on the pull from the faster lamina (#2) first.

### Faster lamina

Assume that we are figuring out the force on the lamina with radius s. From the equation above, we need to know the area of contact and the velocity gradient. Think of the lamina as a cylinder of radius s and thickness ds. The area of contact between the lamina and the faster one is simply the area of the inside of the cylinder: A = 2πsΔx. We don't know the exact form for the velocity of the liquid within the tube yet, but we do know (from our assumption above) that it is dependent on the radius. Therefore, the velocity gradient is the change of the velocity with respect to the change in the radius at the intersection of these two laminae. That intersection is at a radius of s. So, considering that this force will be positive with respect to the movement of the liquid (but the derivative of the velocity is negative), the final form of the equation becomes

$F_\left\{text\left\{viscosity, fast\right\}\right\} = - eta 2 pi s Delta x left . frac\left\{dv\right\}\left\{dr\right\} right vert_s$

where the vertical bar and subscript s following the derivative indicates that it should be taken at a radius of s.

### Slower lamina

Next let's find the force of drag from the slower lamina. We need to calculate the same values that we did for the force from the faster lamina. In this case, the area of contact is at s+ds instead of s. Also, we need to remember that this force opposes the direction of movement of the liquid and will therefore be negative (and that the derivative of the velocity is negative).

$F_\left\{text\left\{viscosity, slow\right\}\right\} = eta 2 pi \left(s+ds\right) Delta x left . frac\left\{dv\right\}\left\{dr\right\} right vert_\left\{s+ds\right\}$

### Putting it all together

To find the solution for the flow of liquid through a tube, we need to make one last assumption. There is no acceleration of liquid in the pipe, and by Newton's first law, there is no net force. If there is no net force then we can add all of the forces together to get zero

$0 = F_\left\{text\left\{pressure\right\}\right\} + F_\left\{text\left\{viscosity, fast\right\}\right\} + F_\left\{text\left\{viscosity, slow\right\}\right\}$

or

$0 = - Delta P2 pi sds - eta 2 pi s Delta x left . frac\left\{dv\right\}\left\{dr\right\} right vert_s + eta 2 pi \left(s+ds\right) Delta x left . frac\left\{dv\right\}\left\{dr\right\} right vert_\left\{s+ds\right\}$

Before we move further, we need to simplify this ugly equation. First, to get everything happening at the same point, we need to do a Taylor series expansion of the velocity gradient, keeping only the linear and quadratic terms (a standard mathematical trick).

$left . frac\left\{dv\right\}\left\{dr\right\} right vert_\left\{r+dr\right\} = left . frac\left\{dv\right\}\left\{dr\right\} right vert_r + left . frac\left\{d^2 v\right\}\left\{dr^2\right\} right vert_r dr$

Let's use this relation in our equation. Also, let's use r instead of s since the lamina we chose was arbitrary and we want our expression to be valid for all laminae. Grouping like terms and dropping the vertical bar since all derivatives are assumed to be at radius r,

$0 = - Delta P2 pi rdr + eta 2 pi dr Delta x frac\left\{dv\right\}\left\{dr\right\} + eta 2 pi r dr Delta x frac\left\{d^2 v\right\}\left\{dr^2\right\} + eta 2 pi \left(dr\right)^2 Delta x frac\left\{d^2 v\right\}\left\{dr^2\right\}$

Finally, let's get this in the form of a differential equation, moving some terms around to make it easier to solve later, and neglecting the term quadratic in dr since this will be really small compared to the rest (another standard mathematical trick).

$frac\left\{1\right\}\left\{eta\right\} frac\left\{Delta P\right\}\left\{Delta x\right\} = frac\left\{d^2 v\right\}\left\{dr^2\right\} + frac\left\{1\right\}\left\{r\right\} frac\left\{dv\right\}\left\{dr\right\}$

It can be seen that both sides of the equations are negative: there is a drop of pressure along the tube (left side) and both first and second derivatives of the velocity are negative (velocity has a maximum value of the center of the tube). The equation may be re-arranged to:

$frac\left\{1\right\}\left\{eta\right\} frac\left\{Delta P\right\}\left\{Delta x\right\} = frac\left\{1\right\}\left\{r\right\} frac\left\{d\right\}\left\{dr\right\} r frac\left\{dv\right\} \left\{dr\right\}$

This differential equation is subject to the following boundary conditions:

v(r) = 0 at r = R -- "No-slip" Boundary Condition at the Wall
$frac\left\{dv\right\} \left\{dr\right\} = 0$ at r = 0 -- Axial symmetry

Axial symmetry means that the velocity v(r) is maximum at the center of the tube, therefore the first derivative $frac\left\{dv\right\}\left\{dr\right\}$ is zero at r = 0.

The differential equation can be integrated to:

$v\left(r\right) = frac\left\{1\right\}\left\{4 eta\right\}r^2frac\left\{Delta P\right\}\left\{Delta x\right\} + A ln\left(r\right) + B$

To find A and B, we use the boundary conditions.

First, the symmetry boundary condition indicates:

$frac\left\{dv\right\}\left\{dr\right\} = frac\left\{1\right\}\left\{2 eta\right\} r frac\left\{Delta P\right\}\left\{Delta x\right\} + A frac\left\{1\right\}\left\{r\right\} = 0$ at r = 0

A solution possible only if A = 0. Next the no-slip boundary condition is applied to the remaining equation:

$v\left(R\right) = frac\left\{1\right\}\left\{4 eta\right\} R^2 frac\left\{Delta P\right\}\left\{Delta x\right\} + B = 0$

so therefore

$B = - frac\left\{1\right\}\left\{4 eta\right\} R^2 frac\left\{Delta P\right\}\left\{Delta x\right\}$

Now we have a formula for the velocity of liquid moving through the tube as a function of the distance from the center of the tube

$v = - frac\left\{1\right\}\left\{4 eta\right\} frac\left\{Delta P\right\}\left\{Delta x\right\} \left(R^2 - r^2\right)$

or, at the center of the tube where the liquid is moving fastest (r = 0) with R being the radius of the tube,

$v_\left\{max\right\} = - frac\left\{1\right\}\left\{4 eta\right\} frac\left\{Delta P\right\}\left\{Delta x\right\}R^2$

### Poiseuille's Law

To get the total volume that flows through the tube, we need to add up the contributions from each lamina. To calculate the flow through each lamina, we multiply the velocity (from above) and the area of the lamina.

$Phi \left(r\right) = frac\left\{1\right\}\left\{4 eta\right\} frac$
{Delta x} (R^2 - r^2) 2 pi rdr = frac{pi}{2 eta} frac{|Delta P{Delta x} (rR^2 - r^3)dr >

Finally, we integrate over all lamina via the radius variable r.

$Phi = frac\left\{pi\right\}\left\{2 eta\right\} frac$>{Delta x} int_{0}^{R} (rR^2 - r^3), dr = frac{|Delta P| pi R^4}{8 eta Delta x}

## Poiseuille's equation for compressible fluids

For a compressible fluid in a tube the volumetric flow rate and the linear velocity is not constant along the tube. The flow is usually expressed at outlet pressure. As fluid is compressed or expands work is done and the fluid is heated and cooled. This meaning that the flow rate is dependant upon heat transfer to and from the fluid. For an ideal gas in the isothermal case, where the temperature of the fluid is permitted to equilibrate with its surroundings, and when the pressure difference between ends of the pipe is small, the volumetric flow rate at the pipe outlet is given by

$Phi = frac\left\{dV\right\}\left\{dt\right\} = v pi R^\left\{2\right\} = frac\left\{pi R^\left\{4\right\} left\left(P_\left\{i\right\}-P_\left\{o\right\} right\right)\right\}\left\{8 eta L\right\} times frac\left\{ P_\left\{i\right\}+P_\left\{o\right\}\right\}\left\{2 P_\left\{o\right\}\right\} = frac\left\{pi R^\left\{4\right\}\right\}\left\{16 eta L\right\} left\left(frac\left\{ P_\left\{i\right\}^\left\{2\right\}-P_\left\{o\right\}^\left\{2\right\}\right\}\left\{P_\left\{o\right\}\right\} right\right)$

Where:

$P_\left\{i\right\}$ inlet pressure
$P_\left\{o\right\}$ outlet pressure
$L$ is the length of tube
$eta$ is the viscosity
$R$ is the radius
$V$ is the volume of the fluid at outlet pressure
$v$ is the velocity of the fluid at outlet pressure

This is usually a good approximation when the flow velocity is less than mach 0.3

This equation can be seen as Poiseuille's law with an extra correction factor $frac\left\{P_\left\{i\right\}+P_\left\{o\right\}\right\}\left\{2\right\} times frac\left\{1\right\}\left\{P_\left\{o\right\}\right\}$ expressing the average pressure relative to the outlet pressure.

## Electrical Circuits analogy

Electricity was originally understood to be a kind of fluid. This hydraulic analogy is still conceptually useful.

Poiseuille's law corresponds to Ohm's law for electrical circuits (V = IR), where the pressure drop ΔP is analogous to the voltage V and voluminal flow rate Φ is analogous to the current I. Then the resistance

$R = frac\left\{ 8 eta Delta x\right\}\left\{pi r^4\right\}$
This concept is useful because the effective resistance in a tube is inversely proportional to the fourth power of the radius. This means that halfing the size of the tube increases the resistance to fluid movement by 16 times.

Both Ohm's law and Poiseuille's law illustrate transport phenomena.

## History

It was developed independently by Gotthilf Heinrich Ludwig Hagen (1797-1884) and Jean Louis Marie Poiseuille.

Poiseuille's law was experimentally derived in 1838 and formulated and published in 1840 and 1846 by Jean Louis Marie Poiseuille (1797-1869). Hagen did his experiments in 1839.

## References

• S. P. Sutera, R. Skalak, "The history of Poiseuille's law," Annual Review of Fluid Mechanics, Vol. 25, 1993, pp. 1-19