Ley de Poiseuille

La ley de Poiseuille (ley de Hagen-Poiseuille) es una ley que permite determinar el flujo laminar estacionario ($\Phi _{V}$) de un líquido incompresible y uniformemente viscoso (fluido newtoniano) a través de un tubo cilíndrico de sección circular constante. La ley es también muy importante en hemodinámica. La ley queda formulada del siguiente modo:

$Q={\frac {dV}{dt}}=\pi \ {\bar {u}}\ R^{2}={\frac {\pi \ R^{4}}{8\ \mu }}\left(-{\frac {dp}{dz}}\right)={\frac {\pi \ R^{4}}{8\ \mu }}{\frac {\Delta p}{L}}$

donde ($V$) es el volumen del líquido que circula en la unidad de tiempo ($t$), (${\bar {u}}$) la velocidad media del fluido a lo largo del eje z del sistema de coordenadas cilíndrico, ($r$) es el radio interno del tubo, ($\Delta p$) es la caída de presión entre los dos extremos, ($\mu$) es la viscosidad dinámica (a veces representada por $\eta$) y ($L$) la longitud característica a lo largo del eje z.

La ley se puede derivar de la ecuación de Darcy-Weisbach, desarrollada en el campo de la hidráulica y que por lo demás es válida para todos los tipos de flujo. La ley de Hagen-Poiseuille se puede expresar también del siguiente modo:

$f={\frac {64}{\mathrm {Re} }}\ ,\quad \quad \mathrm {Re} ={\frac {\rho \ {\bar {u}}\ (2R)}{\mu }}$

donde ($\mathrm {Re}$) es el número de Reynolds y ($\rho$) es la densidad del fluido. En esta forma la ley aproxima el valor del factor de fricción, la energía disipada por la pérdida de carga, el factor de pérdida por fricción o el factor de fricción de Darcy ($f$) en flujo laminar a muy bajas velocidades en un tubo cilíndrico.

Etimología

Esta ecuación fue derivada experimentalmente en 1838, formulada y publicada en 1840 y 1846 por Jean Léonard Marie Poiseuille (1797-1869). Llamada también en honor a Gotthilf Heinrich Ludwig Hagen (1797-1884) por los experimentos llevados a cabo en 1839.

Historia

La derivación teórica de la fórmula original de Poiseuille fue realizada independientemente por Wiedman (1856) y Neumann y E. Hagenbach (1858, 1859, 1860). Hagenbach fue el primero que la denominó como ley de Poiseuille.

La ley de Poiseuille fue extendida en 1891 para flujo turbulento por L. R. Wilberforce, basándose en el trabajo de Hagenbach.

Simbología

Simbología

Símbolo	Nombre	Unidad	Símbolo	Nombre	Unidad
$\mathrm {Re}$	Número de Reynolds			Dimensiones
$F$	Fuerza ejercida por el fluido	N	$A$	Área transversal	m2
$f$	Factor de fricción		$A_{\rm {L}}$	Área longitudinal	m2
$g$	Gravedad	m / s2	$D$	Diámetro máximo	m
$h_{f}$	Perdida por fricción	m	$L$	Longitud	m
$p_{1}$	Presión en punto 1	Pa	$R$	Radio máximo	m
$p_{2}$	Presión en punto 2	Pa	$r$	Radio	m
$\Delta p$	Diferencia de presiones	Pa	$V$	Volumen	m3
$Q$	Caudal	m3 / s	$t$	Tiempo	s
$R_{\rm {hd}}$	Resistencia hidrodinámica			Eléctricos
$u$	Velocidad	m / s	$I$	Corriente	A
${\bar {u}}$	Velocidad media	m / s	$R_{e}$	Resistencia eléctrica	Ω
$u_{\rm {max}}$	Velocidad máxima	m / s	$V_{e}$	Voltaje	Volt
$\mu$	Viscosidad dinámica	Pa s
$\rho$	Densidad	kg / m3
$\tau$	Esfuerzo cortante	N / m2
${\frac {du}{dr}}$	Velocidad de corte	s-1

Cálculo de la fórmula

En una tubería por la cual circula un líquido se considera un trozo de ésta delimitado por los puntos ($ABDC$), de longitud ($L$) y radio interno ($R$).

En el líquido circulante se considera un cilindro coaxial interno delimitado por los puntos ($abdc$) de radio ($r$), con área transversal ($1$) y ($2$): $A=\pi r^{2}$ y área longitudinal: $A_{L}=2\pi rL$ .

Debido a la viscosidad, sobre este cilindro actúa un esfuerzo cortante ($\tau$) provocado por una fuerza cortante ($F$). Esta fuerza ($F=p_{1}A-p_{2}A$) tiene el sentido del desplazamiento del fluido, producido por un gradiente de presión ($p_{1}>p_{2}$).

De acuerdo a la segunda ley de Newton, si ($p_{1}$) y ($p_{2}$) son las presiones aplicadas en el centro de gravedad del área transversal ($1$) y ($2$):

$p_{1}A-p_{2}A+F=0$

Nota	Descripción
(*)	En un sólido, el esfuerzo de corte es proporcional a la deformación, pero un fluido se deforma continuamente mientras se aplique el esfuerzo de corte (${\frac {F}{A_{L}}}$), por lo tanto éste será proporcional a la velocidad de corte (${\frac {du}{dr}}$) por la viscosidad ($\mu$).

Deducción

	1	2	3	4	5*
Ecuaciones	$p_{1}A-p_{2}A+F=0$	$\Delta p=p_{1}-p_{2}$	$A=\pi \ r^{2}$	$A_{\rm {L}}=2\pi \ r\ L$	${\frac {F}{A_{\rm {L}}}}=\mu {\frac {du}{dr}}$
Simplificando	$(p_{1}-p_{2})A+F=0$
Despejando					$F=\mu \ A_{\rm {L}}{\frac {du}{dr}}$
Sustituyendo				$F=(2\pi \ r\ L)\ \mu \ {\frac {du}{dr}}$
Sustituyendo	$(\Delta p)(\pi \ r^{2})+{\Bigl (}2\pi \ r\ L\ \mu \ {\frac {du}{dr}}{\Bigr )}=0$
Simplificando	$\Delta p\ r+2\ L\ \mu \ {\frac {du}{dr}}=0$
Despejando	$du=-{\Bigl (}{\frac {\Delta p}{2\ L\ \mu }}{\Bigr )}\ r\ dr$
Integrando	$u=-{\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}\ r^{2}+C$
Evaluando$r=R,\ u=0$			$0=-{\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}\ (R)^{2}+C$
Despejando			$C={\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}R^{2}$
Sustituyendo	$u=-{\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}\ r^{2}+{\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}R^{2}$
Simplificando	$u={\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}\ (R^{2}-r^{2})$

Distribución de velocidades en una tubería

$u={\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}\ (R^{2}-r^{2})$

El radio elevado al cuadrado indica que se trata de un paraboloide, donde la velocidad máxima ($u_{\rm {max}}$) se obtiene en el eje de la tubería. Zona en la que los efectos del rozamiento con las paredes de la tubería es mínima.

$u_{\rm {max}}={\Bigl (}{\frac {\Delta p}{4\ \mu \ L}}{\Bigr )}R^{2}$

En la práctica es más sencillo medir la velocidad media (${\bar {u}}$) que la velocidad máxima ($u_{\rm {max}}$).

${\bar {u}}={\frac {Q}{\pi \ R^{2}}}$

Para calcular el caudal en la tubería se va a considerar un anillo diferencial de espesor ($dr$) entre dos circunferencias concéntricas con el eje de la tubería y radios ($r$) y ($r+dr$).

Deducción

	1	2
Ecuaciones	$u={\frac {Q}{\pi \ r^{2}}}$	$u={\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}\ (R^{2}-r^{2})$
Despejando	$Q=\pi \ u\ r^{2}$
Derivando	$dQ=(2\pi \ u)\ r\ dr$
Sustituyendo	$dQ=2\pi \ {\Bigl (}{\frac {\Delta p}{4\ L\ \mu }}{\Bigr )}\ (R^{2}-r^{2})\ r\ dr$
Simplificando	$dQ={\Bigl (}{\frac {\pi \ \Delta p}{2\ L\ \mu }}{\Bigr )}\ (R^{2}-r^{2})\ r\ dr$
Integrando	$\int _{0}^{Q}\ dQ=\int _{0}^{R}{\Bigl (}{\frac {\pi \ \Delta p}{2\ L\ \mu }}{\Bigr )}(R^{2}-r^{2})\ r\ dr$
	$\int _{0}^{Q}\ dQ={\Bigl (}{\frac {\pi \ \Delta p}{2\ L\ \mu }}{\Bigr )}\int _{0}^{R}(R^{2}-r^{2})\ r\ dr$
	$\int _{0}^{Q}\ dQ={\Bigl (}{\frac {\pi \ \Delta p}{2\ L\ \mu }}{\Bigr )}\int _{0}^{R}[(R^{2})\ r-r^{3}]\ dr$
	${\Bigr [}Q{\Bigr ]}_{0}^{Q}={\Bigl (}{\frac {\pi \ \Delta p}{2\ L\ \mu }}{\Bigr )}{\Bigr [}{\frac {1}{2}}R^{2}(r)^{2}-{\frac {1}{4}}(r)^{4}{\Bigr ]}_{0}^{R}$
	$Q={\Bigl (}{\frac {\pi \ \Delta p}{2\ L\ \mu }}{\Bigr )}{\Bigl (}{\frac {R^{4}}{2}}-{\frac {R^{4}}{4}}{\Bigr )}$
	$Q={\Bigl (}{\frac {\pi \ \Delta p}{8\ \mu \ L}}{\Bigr )}\ R^{4}$

Ley de Poiseuille para el caudal.

$Q={\Bigl (}{\frac {\pi \ \Delta p}{8\ \mu \ L}}{\Bigr )}\ R^{4}$

Deducción

	1	2	3	4
Ecuaciones	$R={\frac {D}{2}}$	${\bar {u}}={\frac {Q}{\pi \ R^{2}}}$	$Q={\Bigl (}{\frac {\pi \ \Delta p}{8\ \mu \ L}}{\Bigr )}\ R^{4}$	$u_{\rm {max}}={\Bigl (}{\frac {\Delta p}{4\ \mu \ L}}{\Bigr )}R^{2}$
Sustituyendo		${\bar {u}}={\Bigl (}{\frac {\pi \ \Delta p}{8\ \mu \ L}}{\Bigr )}{\frac {R^{4}}{\pi \ R^{2}}}$
Simplificando		${\bar {u}}={\Bigl (}{\frac {\Delta p}{8\ \mu \ L}}{\Bigr )}\ R^{2}$
Sustituyendo	${\bar {u}}={\Bigl (}{\frac {\Delta p}{8\ \mu \ L}}{\Bigr )}\ {\Bigl (}{\frac {D}{2}}{\Bigr )}^{2}$		${\bar {u}}={\frac {u_{\rm {max}}}{2}}$
Simplificando	${\bar {u}}={\Bigl (}{\frac {\Delta p\ D^{2}}{32\ \mu \ L}}{\Bigr )}$
Despejando	$\Delta p={\frac {32\ \mu \ {\bar {u}}\ L}{D^{2}}}$

${\bar {u}}={\frac {u_{\rm {max}}}{2}}$

Ley de Poiseuille para la pérdida de presión.

$\Delta p={\frac {32\ \mu \ {\bar {u}}\ L}{D^{2}}}$

Deducción

	1	2	3
Ecuaciones	$\Delta p={\frac {32\ \mu \ {\bar {u}}\ L}{D^{2}}}$	$h_{f}={\frac {\Delta p}{\rho \ g}}$	$\mathrm {Re} ={\frac {\rho \ {\bar {u}}\ D}{\mu }}$
Dividiendo $(2\ \rho \ g\ {\bar {v}})$	${\frac {\Delta p}{(2\ \rho \ g\ {\bar {u}})}}={\frac {32\ \mu \ {\bar {u}}\ L}{D^{2}(2\ \rho \ g\ {\bar {u}})}}$
Ordenando	${\frac {\Delta p}{\rho \ g}}={\Bigl (}{\frac {L}{D}}{\Bigr )}64{\Bigl (}{\frac {\mu }{\rho \ {\bar {u}}\ D}}{\Bigr )}{\Bigl (}{\frac {{\bar {u}}^{2}}{2\ g}}{\Bigr )}$
Sustituyendo	$h_{f}={\Bigl (}{\frac {L}{D}}{\Bigr )}{\Bigl (}{\frac {64}{\mathrm {Re} }}{\Bigr )}{\Bigl (}{\frac {{\bar {u}}^{2}}{2\ g}}{\Bigr )}$

$h_{f}={\Bigl (}{\frac {L}{D}}{\Bigr )}{\Bigl (}{\frac {64}{\mathrm {Re} }}{\Bigr )}{\Bigl (}{\frac {{\bar {u}}^{2}}{2\ g}}{\Bigr )}$

Al comparar con la ecuación de Darcy-Weisbach se deduce el factor de fricción ($f$).

$f={\frac {64}{\mathrm {Re} }}$

siendo esta otra expresión de la ecuación de Hagen-Poiseuille.

Relación con los circuitos eléctricos

La electricidad fue originalmente entendida como una clase de fluido. Esta analogía hidráulica es todavía útil en el ámbito académico con fines didácticos.

La ley de Ohm para los circuitos eléctricos se corresponde con la ley de Poiseuille.

Equivalencia

Ley de Ohm		Ley de Poiseuille
$R_{e}={\frac {V_{e}}{I}}$		$R_{\rm {hd}}={\frac {\Delta p}{Q}}={\Bigr (}{\frac {8\ \mu \ L}{\pi \ r^{4}}}{\Bigr )}$
$V_{e}$	Voltaje	$\Delta p$	Caída de presión
$I$	Corriente eléctrica	$Q$	Caudal
$R_{e}$	Resistencia eléctrica	${\Bigr (}{\frac {8\ \mu \ L}{\pi \ r^{4}}}{\Bigr )}$	Resistencia hidrodinámica

Relación con el pulmón

La ley de Poiseuille tiene aplicación en la ventilación pulmonar al describir el efecto que tiene el radio de las vías respiratorias sobre la resistencia del flujo de aire en dirección a los alveolos. De ese modo, si el radio de los bronquiolos se redujera por la mitad, la ley de Poiseuille predice que el caudal de aire que pasa por ese bronquiolo reducido tendría que oponerse a una resistencia 16 veces mayor, teniendo en cuenta que la resistencia al flujo es inversamente proporcional al radio elevado a la cuarta potencia.​

Este principio cobra importancia en el asma y otras enfermedades obstructivas del pulmón. Al reducirse el radio de las vías aéreas respiratorias, el esfuerzo de la persona se eleva a la cuarta potencia.

Referencias

• Datos: Q869854
• Multimedia: Hagen-Poiseuille equation / Q869854