La ley de Poiseuille (ley de Hagen-Poiseuille) es una ley que permite determinar el flujo laminar estacionario ( Φ V {displaystyle 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 = d V d t = π   u ¯   R 2 = π   R 4 8   μ ( − d p d z ) = π   R 4 8   μ Δ p L {displaystyle 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 {displaystyle V} ) es el volumen del líquido que circula en la unidad de tiempo ( t {displaystyle t} ), ( u ¯ {displaystyle {bar {u}}} ) la velocidad media del fluido a lo largo del eje z del sistema de coordenadas cilíndrico, ( r {displaystyle r} ) es el radio interno del tubo, ( Δ p {displaystyle Delta p} ) es la caída de presión entre los dos extremos, ( μ {displaystyle mu } ) es la viscosidad dinámica y ( L {displaystyle 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 = 64 R e   , R e = ρ   u ¯   ( 2 R ) μ {displaystyle f={frac {64}{mathrm {Re} }} ,quad quad mathrm {Re} ={frac {rho {bar {u}} (2R)}{mu }}}

donde ( R e {displaystyle mathrm {Re} } ) es el número de Reynolds y ( ρ {displaystyle 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 {displaystyle f} ) en flujo laminar a muy bajas velocidades en un tubo cilíndrico.

Etimología

[

editar

]

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

[

editar

]

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

[

editar

]

SimbologíaSímboloNombreUnidadSímboloNombreUnidad

R e {displaystyle mathrm {Re} }

Número de ReynoldsDimensiones

F {displaystyle F}

Fuerza ejercida por el fluidoN

A {displaystyle A}

Área transversalm2

f {displaystyle f}

Factor de fricción

A L {displaystyle A_{rm {L}}}

Área longitudinalm2

g {displaystyle g}

Gravedadm / s2

D {displaystyle D}

Diámetro máximom

h f {displaystyle h_{f}}

Perdida por fricciónm

L {displaystyle L}

Longitudm

p 1 {displaystyle p_{1}}

Presión en punto 1Pa

R {displaystyle R}

Radio máximom

p 2 {displaystyle p_{2}}

Presión en punto 2Pa

r {displaystyle r}

Radiom

Δ p {displaystyle Delta p}

Diferencia de presionesPa

V {displaystyle V}

Volumenm3

Q {displaystyle Q}

Caudalm3 / s

t {displaystyle t}

Tiempos

R h d {displaystyle R_{rm {hd}}}

Resistencia hidrodinámicaEléctricos

u {displaystyle u}

Velocidadm / s

I {displaystyle I}

CorrienteA

u ¯ {displaystyle {bar {u}}}

Velocidad mediam / s

R e {displaystyle R_{e}}

Resistencia eléctricaΩ

u m a x {displaystyle u_{rm {max}}}

Velocidad máximam / s

V e {displaystyle V_{e}}

VoltajeVolt

μ {displaystyle mu }

Viscosidad dinámicaPa s

ρ {displaystyle rho }

Densidadkg / m3

τ {displaystyle tau }

Esfuerzo cortanteN / m2

d u d r {displaystyle {frac {du}{dr}}}

Velocidad de cortes-1

Cálculo de la fórmula

[

editar

]

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

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

Debido a la viscosidad, sobre este cilindro actúa un esfuerzo cortante ( τ {displaystyle tau } ) provocado por una fuerza cortante ( F {displaystyle F} ). Esta fuerza ( F = p 1 A − p 2 A {displaystyle 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 {displaystyle p_{1}>p_{2}} ).

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

p 1 A − p 2 A + F = 0 {displaystyle p_{1}A-p_{2}A+F=0}

NotaDescripció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 (

F A L {displaystyle {frac {F}{A_{L}}}}

velocidad de corte (

d u d r {displaystyle {frac {du}{dr}}}

viscosidad (

μ {displaystyle mu }

Deducción12345*Ecuaciones

p 1 A − p 2 A + F = 0 {displaystyle p_{1}A-p_{2}A+F=0}

Δ p = p 1 − p 2 {displaystyle Delta p=p_{1}-p_{2}}

A = π   r 2 {displaystyle A=pi r^{2}}

A L = 2 π   r   L {displaystyle A_{rm {L}}=2pi r L}

F A L = μ d u d r {displaystyle {frac {F}{A_{rm {L}}}}=mu {frac {du}{dr}}}

Simplificando

( p 1 − p 2 ) A + F = 0 {displaystyle (p_{1}-p_{2})A+F=0}

Despejando

F = μ   A L d u d r {displaystyle F=mu A_{rm {L}}{frac {du}{dr}}}

Sustituyendo

F = ( 2 π   r   L )   μ   d u d r {displaystyle F=(2pi r L) mu {frac {du}{dr}}}

Sustituyendo

( Δ p ) ( π   r 2 ) + ( 2 π   r   L   μ   d u d r ) = 0 {displaystyle (Delta p)(pi r^{2})+{Bigl (}2pi r L mu {frac {du}{dr}}{Bigr )}=0}

Simplificando

Δ p   r + 2   L   μ   d u d r = 0 {displaystyle Delta p r+2 L mu {frac {du}{dr}}=0}

Despejando

d u = − ( Δ p 2   L   μ )   r   d r {displaystyle du=-{Bigl (}{frac {Delta p}{2 L mu }}{Bigr )} r dr}

Integrando

u = − ( Δ p 4   L   μ )   r 2 + C {displaystyle u=-{Bigl (}{frac {Delta p}{4 L mu }}{Bigr )} r^{2}+C}

Evaluando

r = R ,   u = 0 {displaystyle r=R, u=0}

0 = − ( Δ p 4   L   μ )   ( R ) 2 + C {displaystyle 0=-{Bigl (}{frac {Delta p}{4 L mu }}{Bigr )} (R)^{2}+C}

Despejando

C = ( Δ p 4   L   μ ) R 2 {displaystyle C={Bigl (}{frac {Delta p}{4 L mu }}{Bigr )}R^{2}}

Sustituyendo

u = − ( Δ p 4   L   μ )   r 2 + ( Δ p 4   L   μ ) R 2 {displaystyle u=-{Bigl (}{frac {Delta p}{4 L mu }}{Bigr )} r^{2}+{Bigl (}{frac {Delta p}{4 L mu }}{Bigr )}R^{2}}

Simplificando

u = ( Δ p 4   L   μ )   ( R 2 − r 2 ) {displaystyle u={Bigl (}{frac {Delta p}{4 L mu }}{Bigr )} (R^{2}-r^{2})}

Distribución de velocidades en una tubería

u = ( Δ p 4   L   μ )   ( R 2 − r 2 ) {displaystyle 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 m a x {displaystyle 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 m a x = ( Δ p 4   μ   L ) R 2 {displaystyle 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 ( u ¯ {displaystyle {bar {u}}} ) que la velocidad máxima ( u m a x {displaystyle u_{rm {max}}} ).

u ¯ = Q π   R 2 {displaystyle {bar {u}}={frac {Q}{pi R^{2}}}}

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

Deducción12Ecuaciones

u = Q π   r 2 {displaystyle u={frac {Q}{pi r^{2}}}}

u = ( Δ p 4   L   μ )   ( R 2 − r 2 ) {displaystyle u={Bigl (}{frac {Delta p}{4 L mu }}{Bigr )} (R^{2}-r^{2})}

Despejando

Q = π   u   r 2 {displaystyle Q=pi u r^{2}}

Derivando

d Q = ( 2 π   u )   r   d r {displaystyle dQ=(2pi u) r dr}

Sustituyendo

d Q = 2 π   ( Δ p 4   L   μ )   ( R 2 − r 2 )   r   d r {displaystyle dQ=2pi {Bigl (}{frac {Delta p}{4 L mu }}{Bigr )} (R^{2}-r^{2}) r dr}

Simplificando

d Q = ( π   Δ p 2   L   μ )   ( R 2 − r 2 )   r   d r {displaystyle dQ={Bigl (}{frac {pi Delta p}{2 L mu }}{Bigr )} (R^{2}-r^{2}) r dr}

Integrando

∫ 0 Q   d Q = ∫ 0 R ( π   Δ p 2   L   μ ) ( R 2 − r 2 )   r   d r {displaystyle int _{0}^{Q} dQ=int _{0}^{R}{Bigl (}{frac {pi Delta p}{2 L mu }}{Bigr )}(R^{2}-r^{2}) r dr}

∫ 0 Q   d Q = ( π   Δ p 2   L   μ ) ∫ 0 R ( R 2 − r 2 )   r   d r {displaystyle int _{0}^{Q} dQ={Bigl (}{frac {pi Delta p}{2 L mu }}{Bigr )}int _{0}^{R}(R^{2}-r^{2}) r dr}

∫ 0 Q   d Q = ( π   Δ p 2   L   μ ) ∫ 0 R [ ( R 2 )   r − r 3 ]   d r {displaystyle int _{0}^{Q} dQ={Bigl (}{frac {pi Delta p}{2 L mu }}{Bigr )}int _{0}^{R}[(R^{2}) r-r^{3}] dr}

[ Q ] 0 Q = ( π   Δ p 2   L   μ ) [ 1 2 R 2 ( r ) 2 − 1 4 ( r ) 4 ] 0 R {displaystyle {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 = ( π   Δ p 2   L   μ ) ( R 4 2 − R 4 4 ) {displaystyle Q={Bigl (}{frac {pi Delta p}{2 L mu }}{Bigr )}{Bigl (}{frac {R^{4}}{2}}-{frac {R^{4}}{4}}{Bigr )}}

Q = ( π   Δ p 8   μ   L )   R 4 {displaystyle Q={Bigl (}{frac {pi Delta p}{8 mu L}}{Bigr )} R^{4}}

Ley de Poiseuille para el caudal.

Q = ( π   Δ p 8   μ   L )   R 4 {displaystyle Q={Bigl (}{frac {pi Delta p}{8 mu L}}{Bigr )} R^{4}}

Deducción1234Ecuaciones

R = D 2 {displaystyle R={frac {D}{2}}}

u ¯ = Q π   R 2 {displaystyle {bar {u}}={frac {Q}{pi R^{2}}}}

Q = ( π   Δ p 8   μ   L )   R 4 {displaystyle Q={Bigl (}{frac {pi Delta p}{8 mu L}}{Bigr )} R^{4}}

u m a x = ( Δ p 4   μ   L ) R 2 {displaystyle u_{rm {max}}={Bigl (}{frac {Delta p}{4 mu L}}{Bigr )}R^{2}}

Sustituyendo

u ¯ = ( π   Δ p 8   μ   L ) R 4 π   R 2 {displaystyle {bar {u}}={Bigl (}{frac {pi Delta p}{8 mu L}}{Bigr )}{frac {R^{4}}{pi R^{2}}}}

Simplificando

u ¯ = ( Δ p 8   μ   L )   R 2 {displaystyle {bar {u}}={Bigl (}{frac {Delta p}{8 mu L}}{Bigr )} R^{2}}

Sustituyendo

u ¯ = ( Δ p 8   μ   L )   ( D 2 ) 2 {displaystyle {bar {u}}={Bigl (}{frac {Delta p}{8 mu L}}{Bigr )} {Bigl (}{frac {D}{2}}{Bigr )}^{2}}

u ¯ = u m a x 2 {displaystyle {bar {u}}={frac {u_{rm {max}}}{2}}}

Simplificando

u ¯ = ( Δ p   D 2 32   μ   L ) {displaystyle {bar {u}}={Bigl (}{frac {Delta p D^{2}}{32 mu L}}{Bigr )}}

Despejando

Δ p = 32   μ   u ¯   L D 2 {displaystyle Delta p={frac {32 mu {bar {u}} L}{D^{2}}}}

u ¯ = u m a x 2 {displaystyle {bar {u}}={frac {u_{rm {max}}}{2}}}

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

Δ p = 32   μ   u ¯   L D 2 {displaystyle Delta p={frac {32 mu {bar {u}} L}{D^{2}}}}

Deducción123Ecuaciones

Δ p = 32   μ   u ¯   L D 2 {displaystyle Delta p={frac {32 mu {bar {u}} L}{D^{2}}}}

h f = Δ p ρ   g {displaystyle h_{f}={frac {Delta p}{rho g}}}

R e = ρ   u ¯   D μ {displaystyle mathrm {Re} ={frac {rho {bar {u}} D}{mu }}}

Dividiendo

( 2   ρ   g   v ¯ ) {displaystyle (2 rho g {bar {v}})}

Δ p ( 2   ρ   g   u ¯ ) = 32   μ   u ¯   L D 2 ( 2   ρ   g   u ¯ ) {displaystyle {frac {Delta p}{(2 rho g {bar {u}})}}={frac {32 mu {bar {u}} L}{D^{2}(2 rho g {bar {u}})}}}

Ordenando

Δ p ρ   g = ( L D ) 64 ( μ ρ   u ¯   D ) ( u ¯ 2 2   g ) {displaystyle {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 = ( L D ) ( 64 R e ) ( u ¯ 2 2   g ) {displaystyle h_{f}={Bigl (}{frac {L}{D}}{Bigr )}{Bigl (}{frac {64}{mathrm {Re} }}{Bigr )}{Bigl (}{frac {{bar {u}}^{2}}{2 g}}{Bigr )}}

h f = ( L D ) ( 64 R e ) ( u ¯ 2 2   g ) {displaystyle 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 {displaystyle f} ).

f = 64 R e {displaystyle f={frac {64}{mathrm {Re} }}}

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

Relación con los circuitos eléctricos

[

editar

]

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.

EquivalenciaLey de OhmLey de Poiseuille

R e = V e I {displaystyle R_{e}={frac {V_{e}}{I}}}

R h d = Δ p Q = ( 8   μ   L π   r 4 ) {displaystyle R_{rm {hd}}={frac {Delta p}{Q}}={Bigr (}{frac {8 mu L}{pi r^{4}}}{Bigr )}}

V e {displaystyle V_{e}}

Voltaje

Δ p {displaystyle Delta p}

Caída de presión

I {displaystyle I}

Corriente eléctrica

Q {displaystyle Q}

Caudal

R e {displaystyle R_{e}}

Resistencia eléctrica

( 8   μ   L π   r 4 ) {displaystyle {Bigr (}{frac {8 mu L}{pi r^{4}}}{Bigr )}}

Resistencia hidrodinámica

Relación con el pulmón

[

editar

]

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.[1]​

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

[

editar

]