Hagen-Poiseuille Viscous Laminar Flow & Pressure Drop Calculator

📁 Physics🔒 100% Free & Private⚡ Instant Output

Calculate viscous fluid dynamics in circular pipes and microfluidic channels. Determine pressure drop ($Delta P$), volumetric discharge ($Q$), parabolic velocity profile ($u(r)$), wall shear stress ($tau_w$), pumping power, and verify laminar flow conditions via Reynolds number ($Re$).

Select dependent variable to solve
Standard fluid dynamic viscosity and density
$text{mPa}cdottext{s}$
Fluid shear viscosity ($1text{ mPa}cdottext{s} = 0.001text{ Pa}cdottext{s} = 1text{ cP}$)
$text{kg/m}^3$
Mass density for Reynolds number computation
mm
Internal conduit bore diameter
m
Total axial length of straight pipe
mL/min
Flow discharge through pipe
Total Pressure Drop ($Delta P$):
-- Pa
-- bar (-- psi)
Volumetric Flow Rate ($Q$):
-- mL/min
-- m³/s (-- L/h)
Flow Velocity (Mean & Peak):
-- m/s
Peak centerline u_max: -- m/s
Reynolds Number ($Re$):
--
Regime: Laminar (Re < 2300)
Wall Shear Stress ($tau_w$):
-- Pa
Darcy friction f: --
Hydraulic Resistance & Power:
-- Pa·s/m³
Pumping Power: -- mW

Parabolic Radial Velocity Profile $u(r)$

Axial Pressure Gradient & Head Loss along Pipe $L$

Quick Overview & Direct Answer

Hagen-Poiseuille Laminar Flow &amp; Pressure Drop Calculator: Free Hagen-Poiseuille Laminar Flow & Pressure Drop Calculator tool: calculate viscous laminar pipe flow pressure drop, volumetric flow rate Q, parabolic.

Primary Use Physical formula evaluation, mechanics, thermodynamics calculations
Input Formats Physical values, mass, velocity, acceleration, energy variables
Output Format Standardized SI units, converted physics quantities
Processing Mode 100% Client-Side Mathematical Engine
Access Tier Free, unlimited, no sign-up needed

Comprehensive Fluid Mechanics Guide to Hagen-Poiseuille Flow

Independently derived by German civil engineer Gotthilf Hagen in 1839 and French physician Jean Léonard Marie Poiseuille in 1840 (who was investigating human cardiovascular hemodynamics), the **Hagen-Poiseuille equation** describes the exact analytical solution of the Navier-Stokes equations for steady, laminar, axisymmetric viscous flow of an incompressible Newtonian fluid through a constant-diameter circular pipe.

1. Exact Derivation from the Navier-Stokes Equations

Consider a cylindrical coordinate system $(r, theta, x)$ aligned with the pipe centerline. Under steady, laminar, fully-developed flow conditions, the radial and azimuthal velocity components vanish ($v_r = 0, v_theta = 0$), and the axial velocity $u(r)$ depends exclusively on radius $r$.

The axial momentum Navier-Stokes equation reduces to a balance between the pressure gradient and viscous shear forces:

$$frac{1}{r} frac{d}{dr} left( r frac{du}{dr} right) = frac{1}{mu} frac{dP}{dx} = -frac{Delta P}{mu L}$$

Integrating twice with respect to radius $r$:

$$r frac{du}{dr} = -frac{Delta P}{2 mu L} r^2 + C_1$$ $$frac{du}{dr} = -frac{Delta P}{2 mu L} r + frac{C_1}{r}$$ $$u(r) = -frac{Delta P}{4 mu L} r^2 + C_1 ln(r) + C_2$$

Applying physical boundary conditions:

  1. Finite Centerline Velocity ($r = 0$): Because $ln(0) to -infty$, constant $C_1 = 0$ to prevent singular velocity on the pipe axis.
  2. No-Slip Boundary Condition at Wall ($r = R$): Viscous fluid adheres to the pipe boundary, yielding $u(R) = 0$:
    $$0 = -frac{Delta P}{4 mu L} R^2 + C_2 implies C_2 = frac{Delta P R^2}{4 mu L}$$

Substituting $C_2$ gives the iconic **parabolic velocity profile**:

$$u(r) = frac{Delta P}{4 mu L} (R^2 - r^2) = u_{max} left( 1 - frac{r^2}{R^2} right)$$

where the peak velocity on the centerline is $u_{max} = frac{Delta P R^2}{4 mu L}$.

2. Volumetric Flow Rate and the Fourth-Power Law

The total volumetric discharge $Q$ is computed by integrating the parabolic velocity across differential concentric annular rings $dA = 2pi r, dr$:

$$Q = int_0^R u(r) cdot 2pi r, dr = 2pi frac{Delta P}{4 mu L} int_0^R (R^2 r - r^3), dr = frac{pi Delta P}{2 mu L} left[ frac{R^4}{2} - frac{R^4}{4} right]$$ $$Q = frac{pi R^4 Delta P}{8 mu L} = frac{pi D^4 Delta P}{128 mu L}quad [text{m}^3/text{s}]$$

Solving directly for the pressure drop yields the **Hagen-Poiseuille equation**:

$$Delta P = frac{8 mu L Q}{pi R^4} = frac{128 mu L Q}{pi D^4}quad [text{Pascals}]$$

The mean cross-sectional flow velocity is exactly half the centerline maximum:

$$bar{u} = frac{Q}{pi R^2} = frac{Delta P R^2}{8 mu L} = frac{u_{max}}{2}$$

3. Wall Shear Stress & Hydraulic Resistance Analogy

According to Newton's law of viscosity, fluid shear stress is proportional to the radial velocity gradient:

$$tau(r) = -mu frac{du}{dr} = -mu left( -frac{Delta P}{2 mu L} r right) = frac{Delta P cdot r}{2 L}$$

At the inner conduit boundary ($r = R$), the wall shear stress reaches its maximum:

$$tau_w = frac{Delta P cdot R}{2 L} = frac{4 mu Q}{pi R^3} = frac{32 mu Q}{pi D^3}quad [text{Pa}]$$

In microfluidics and hemodynamic circuit modeling, the ratio of pressure drop to volumetric flow rate defines the **hydraulic resistance** $R_{hyd}$, serving as the exact fluidic equivalent of electrical resistance in Ohm's Law ($Delta P = Q cdot R_{hyd} leftrightarrow Delta V = I cdot R$):

$$R_{hyd} = frac{Delta P}{Q} = frac{8 mu L}{pi R^4} = frac{128 mu L}{pi D^4}quad [text{Pa}cdottext{s/m}^3]$$

4. Validity Limits: Reynolds Number & Hydrodynamic Entrance Length

The Hagen-Poiseuille equation is strictly applicable only when flow remains in the stable laminar regime:

$$Re = frac{rho bar{u} D}{mu} = frac{4 rho Q}{pi mu D} le 2300$$

If $Re > 2300$, turbulent eddies emerge, causing viscous dissipation to scale with $bar{u}^2$ rather than $bar{u}$, rendering the linear Hagen-Poiseuille law invalid (requiring the Colebrook-White or Darcy-Weisbach formulations).

Furthermore, near the pipe entrance, boundary layers grow along the walls until meeting at the centerline. The entrance length $L_e$ required to establish the fully-developed parabolic profile is:

$$L_e approx 0.05 cdot Re cdot D$$

The conduit length $L$ must be significantly larger than $L_e$ ($L gg L_e$) for Hagen-Poiseuille friction to dominate over entrance pressure losses.

5. Step-by-Step Worked Laminar Flow Calculation

Consider water flowing through a precision stainless steel capillary conduit:

  • Fluid: Pure water at $20^circtext{C}$ ($mu = 1.002text{ mPa}cdottext{s} = 1.002 times 10^{-3}text{ Pa}cdottext{s}$, $rho = 998.2text{ kg/m}^3$)
  • Conduit internal diameter: $D = 4.0text{ mm} = 4.0 times 10^{-3}text{ m}$ (radius $R = 2.0text{ mm} = 2.0 times 10^{-3}text{ m}$)
  • Conduit length: $L = 2.50text{ m}$
  • Volumetric discharge: $Q = 120.0text{ mL/min}$

Step 1: Convert Units and Calculate Mean & Peak Velocity

$$Q = 120.0text{ mL/min} = frac{120.0 times 10^{-6}text{ m}^3}{60text{ s}} = 2.0 times 10^{-6}text{ m}^3/text{s}$$ $$A = frac{pi D^2}{4} = frac{pi (4.0 times 10^{-3})^2}{4} = frac{pi cdot 1.6 times 10^{-5}}{4} approx 1.2566 times 10^{-5}text{ m}^2$$ $$bar{u} = frac{Q}{A} = frac{2.0 times 10^{-6}}{1.2566 times 10^{-5}} approx 0.15915text{ m/s}$$ $$u_{max} = 2 cdot bar{u} = 2 cdot 0.15915 approx 0.3183text{ m/s}$$

Step 2: Verify Laminar Flow via Reynolds Number ($Re$)

$$Re = frac{rho bar{u} D}{mu} = frac{998.2 cdot 0.15915 cdot (4.0 times 10^{-3})}{1.002 times 10^{-3}} = frac{0.6355}{0.001002} approx 634.2$$ $$text{Since } Re = 634.2 < 2300, text{ the flow is unconditionally laminar.}$$ $$L_e = 0.05 cdot (634.2) cdot (0.004text{ m}) approx 0.127text{ m}quad (12.7text{ cm} ll 2.50text{ m})$$

Step 3: Calculate Pressure Drop ($Delta P$)

$$R^4 = (2.0 times 10^{-3})^4 = 1.6 times 10^{-11}text{ m}^4$$ $$Delta P = frac{8 cdot (1.002 times 10^{-3}) cdot 2.50 cdot (2.0 times 10^{-6})}{pi cdot (1.6 times 10^{-11})} = frac{4.008 times 10^{-8}}{5.0265 times 10^{-11}} approx 797.37text{ Pa}$$ $$Delta P approx 0.00797text{ bar}quad (0.1156text{ psi})$$

Step 4: Compute Wall Shear Stress & Pumping Power

$$tau_w = frac{Delta P cdot R}{2 L} = frac{797.37 cdot (2.0 times 10^{-3})}{2 cdot 2.50} = frac{1.5947}{5.0} approx 0.3189text{ Pa}$$ $$f_{Darcy} = frac{64}{Re} = frac{64}{634.2} approx 0.1009$$ $$R_{hyd} = frac{Delta P}{Q} = frac{797.37}{2.0 times 10^{-6}} approx 3.987 times 10^8text{ Pa}cdottext{s/m}^3$$ $$P_{pump} = Q cdot Delta P = (2.0 times 10^{-6}text{ m}^3/text{s}) cdot 797.37text{ Pa} approx 0.001595text{ W} = 1.595text{ mW}$$

The fluid sustains an exact laminar pressure drop of $797.4text{ Pa}$ with wall shear stress of $0.32text{ Pa}$, demanding only $1.6text{ mW}$ of mechanical pumping power.

Frequently Asked Questions

The Hagen-Poiseuille law assumes: (1) laminar flow (Reynolds number Re < 2300), (2) incompressible Newtonian fluid with constant dynamic viscosity, (3) steady, fully-developed, axisymmetric flow through a rigid cylindrical pipe of uniform circular cross-section, and (4) no-slip boundary conditions at the pipe walls.

In laminar cylindrical flow, cross-sectional area scales as R^2, and the mean flow velocity scales as R^2 because viscous drag scales inversely with distance from the boundary. Multiplying velocity by area yields the iconic fourth-power dependency Q ∝ R^4, meaning halving pipe radius increases required driving pressure by sixteen times (16x) for identical flow.

The hydrodynamic entrance length Le is the development distance required from the pipe inlet for the boundary layer to merge and establish the permanent parabolic velocity profile. For laminar flow, it is approximated as Le ≈ 0.05 * Re * D.

Explore other useful tools related to this one.

Reynolds Number, Pipe Flow Regime & Pressure Drop Calculator

Calculate Reynolds number (Re), laminar vs turbulent flow regime, Darcy-Weisbach friction factor, head loss, and pipe pressure drop using the Colebrook-White and Swamee-Jain equations.

Fanno Flow Calculator (Duct Friction & Choking)

Calculate compressible Fanno flow in constant-area adiabatic ducts with wall friction. Determine choking length L*, downstream Mach number, pressure drops, and T-s Fanno lines.

Mach Number, Compressible Flow & Isentropic Stagnation Calculator

Calculate Mach number, speed of sound, isentropic stagnation temperature/pressure ratios (T0/T, P0/P), normal shock waves, and Mach cone angles across gases and flight regimes.

Orbital Velocity, Period & Hohmann Transfer Calculator

Calculate orbital speed, Keplerian period, escape velocity, Vis-Viva equation, and Hohmann transfer orbit delta-v across LEO, GEO, Moon, Mars, and custom.

Stress & Strain Calculator – Tensile, Compressive & Shear Solver

Calculate engineering stress, normal strain, shear modulus, deformation, and factor of safety. Fast, free, and secure client-side execution in your browser.

Tsiolkovsky Rocket Equation & Multi-Stage Delta-V Calculator

Calculate rocket delta-v (Δv), mass ratios, specific impulse (Isp), multi-stage velocity budgets, and payload fractions using the Tsiolkovsky Rocket Equation.