Hagen-Poiseuille Flow: Viscous Fluid Dynamics & Pipe Flow

01/ Historical Context & Fundamentals

In the mid-19th century, German engineer **Gotthilf Hagen** (1839) and French physiologist **Jean Léonard Marie Poiseuille** (1840) independently investigated the behavior of water flowing through narrow tubes. While Hagen was interested in general hydraulics and civil engineering, Poiseuille’s research was driven by a biological mystery: the physical laws governing blood circulation through the human cardiovascular network (hemodynamics).

Their combined work resulted in the **Hagen-Poiseuille Law**, which describes the relationship between pressure drop ($\Delta p$), fluid viscosity ($\mu$), pipe length ($L$), and volumetric flow rate ($Q$) for an incompressible, Newtonian fluid undergoing laminar flow through a uniform cylinder. This equation remains a cornerstone of mechanical engineering, clinical medicine, and chemical transport design today.

Interactive Pipe Flow Arena

Live WebAssembly Render (Rust + Macroquad + Egui)

Initializing hagen-poiseuille.wasm...

WASM Hardware Accelerated 2D/3D Pipe Flow Simulator

Figure 1: Live Hagen-Poiseuille simulation. Drag to orbit the 3D velocity profile and boundary shear stress vectors.

02/ Mathematical Derivations & Profiles

To formulate the velocity distribution of a fluid layer, we consider a cylinder of length $L$ and radius $R$. For a steady, laminar, and uniform flow, we balance the pressure force pushing the fluid forward against the viscous shear force resisting the motion:

$$\Delta p \cdot (\pi r^2) = \tau(r) \cdot (2\pi r L)$$

Here, $\tau(r)$ represents the shear stress at radius $r$. Substituting Newton’s Law of Viscosity, $\tau(r) = -\mu \frac{dv}{dr}$, we get:

$$\frac{dv}{dr} = -\frac{\Delta p \cdot r}{2 \mu L}$$

Integrating this differential equation from the pipe boundary (where $v(R) = 0$ due to the **No-Slip Condition**) to radius $r$, we derive the classic **parabolic velocity profile**:

$$v(r) = \frac{\Delta p}{4 \mu L} \left( R^2 - r^2 \right) = v_{\text{max}} \left( 1 - \frac{r^2}{R^2} \right)$$

Where the maximum centerline velocity is $v_{\text{max}} = \frac{\Delta p \, R^2}{4 \mu L}$. Integrating the velocity profile across the cross-sectional area of the cylinder gives the volumetric flow rate ($Q$):

$$Q = \int_{0}^{R} v(r) \cdot 2\pi r \, dr = \frac{\pi \Delta p R^4}{8 \mu L}$$

03/ The Fourth-Power Relationship

Expressing the volumetric flow rate equation in terms of pipe diameter $D = 2R$, we obtain the standard **Hagen-Poiseuille Equation**:

$$\Delta p = \frac{128 \mu L Q}{\pi D^4}$$

By defining a fluidic resistance variable $R_{\text{flow}}$ analogous to electrical resistance ($R = V/I$ where pressure difference $\Delta p$ acts as voltage and flow rate $Q$ acts as current):

$$R_{\text{flow}} = \frac{\Delta p}{Q} = \frac{128 \mu L}{\pi D^4}$$

This formulation highlights the extreme sensitivity of flow dynamics to the channel width ($D$):

Narrowing (Constriction): Reducing the diameter of a tube by half ($D \to 0.5D$) increases the fluidic resistance by a factor of 16 ($2^4 = 16$). Under a constant flow rate, the pressure drop must increase by **1,600%** to compensate.
Dilation (Expansion): Increasing the pipe diameter by a factor of two ($D \to 2D$) drops resistance to just 6.25% of its initial value, allowing 16 times more fluid volume to flow under the same pressure head.

04/ Viscosity and Temperature dependency

Viscosity ($\mu$) represents the internal friction coefficient of a fluid, dictating how easily its concentric layers slide over one another. In liquid Newtonian fluids, heating the substance breaks the intermolecular attraction forces, decreasing viscosity exponentially. The simulator models this temperature dependence via the following exponential function:

$$\mu(T) = \mu_0 \cdot e^{-0.025 (T - 20)}$$

Where $\mu_0$ is the baseline viscosity at room temperature ($20^\circ\text{C}$). This model captures how heating thick industrial crude oil or human blood increases their flow velocity under identical pumping pressures.