Tema 4

Flujo Laminar
Flujos unidireccionales con viscosidad
dominante

1
 Movimiento laminar con viscosidad
dominante
4.1 Introducción
4.2 Movimiento laminar permanente
4.2.1 Corrientes de Hagen-Poiseuille
4.2.2 En conductos de sección circular
4.2.3 Otras secciones
4.4 Pérdida de carga
4.4.1 Coeficiente de fricción
4.3 Efecto la longitud finita del tubo
4.5 Estabilidad de corriente laminar

2
4.1 Introducción
Laminar vs Turbulento

Fuente: You Tube. URL: http://youtu.be/3na9oezJ_Yk

Autor: Terry Jonhson
3
 Laminar vs Turbulento

Fuente: Frank M. White. Fluid Mechanics.

4.1 Introducción 4
 Diferencias
Flujo Laminar Flujo turbulento
• El fluido se mueve en capas o • Fluctuaciones rápidas tanto
láminas. espaciales como temporales de
• Las inestabilidades en el flujo son las magnitudes físicas.
atenuadas por las fuerzas • La perturbaciones se amplifican
viscosas. haciendo el flujo muy irregular.
• Las fuerzas viscosas predominan • Pierde al carácter estacionario y
sobre las fuerzas de inercia. unidireccional, tomando carácter
aleatorio.
Fuente: Frank M. White. Fluid Mechanics.

4.1 Introducción 5
 Referencias históricas
G. H. L. Hagen Osborne Reynolds
• Ingeniero alemán (1797-1884) • Ingeniero británico (1842-
• En 1839: 1912)
– Indicó la existencia de 2 • En 1883 estableció que el
regímenes de flujo viscoso. cambio de flujo dependía del
– Midió la caída de presión en nº de Reynolds.
tubos largos de latón. 𝜌𝑉𝐿

Fuente: Frank M. White. Fluid Mechanics.

𝑅𝑒 =
𝐿𝑄 𝜇
Δ𝑝 = 𝑐𝑜𝑛𝑠𝑡 + 𝑒𝑓𝑒𝑐𝑡𝑜𝑠 𝑑𝑒 𝑒𝑛𝑡𝑟𝑎𝑑𝑎
𝑟4
– Observó que la ley no era
válida cuando Q superaba
cierto límite.

4.1 Introducción 6
Transición
• Re<Red, crit → Flujo
laminar
• Re>Red, crit → Flujo
turbulento
• Re≈Red, crit → Zona de
Fuente: Frank M. White. Fluid Mechanics.
transición
• En tubo comercial el
valor aceptado es:
𝑅𝑒𝑑,𝑐𝑟𝑖𝑡 ≈ 2300

4.1 Introducción 7
 𝜌𝑉𝐿 𝑉𝐿
Parámetros 𝑅𝑒 =
𝜇
=
𝜈

Aceite - Repsol
Fluido Agua (20º) Agua (20º) Agua (20º) Aire (1 atm, 50º)
Telex E - 32
ρ [Kg/m3] 998 998 998 870 1,09

μ [N·s/m2] 1,003E-03 1,003E-03 1,003E-03 2,784E-02 1,950E-05

μ [cP] 1,003 1,003 1,003 27,840 0,0195

ν [m2/s] 1,005E-06 1,005E-06 1,005E-06 3,200E-05 1,789E-05

ν [cSt] 1,005 1,005 1,005 32 17,89

Re d,crit 2300 2300 2300 2300 2300

V [m/s] 1 1 0,05 1 1

D [mm] 25,4 1 25,4 25,4 25,4

Re 25273 995 1264 794 1420

Tipo de flujo Turbulento Laminar Laminar Laminar Laminar

4.1 Introducción 8
 Ecuaciones del movimiento
• Flujo con viscosidad dominante – Flujo laminar.
• Flujo incompresible.
𝜕𝜌
𝜌 = 𝑐𝑡𝑒 ⇒ = 0 ⇒ 𝛻 · 𝑣Ԧ = 0
𝜕𝑡
• Flujo unidireccional.
𝜕𝑣𝑥
𝑣Ԧ = 𝑣𝑥 𝑖Ԧ ⇒ 𝛻 · 𝑣Ԧ = = 0 ⇒ 𝑣Ԧ = 𝑣𝑥 𝑦, 𝑧, 𝑡 𝑖Ԧ
𝜕𝑥
• Las variaciones en la velocidad del flujo tendrán
lugar en el plano normal a la dirección del
movimiento.

4.2 Movimiento laminar permanente 9

 Ecuaciones del movimiento
• Ecuación de Navier-Stokes para fluidos incompresibles
𝜕𝑣Ԧ
𝜌 + 𝜌𝑣𝛻
Ԧ 𝑣Ԧ = −𝛻 𝑝 + 𝜌𝑈 + 𝜇∆𝑣Ԧ
𝜕𝑡
𝜕𝑣𝑖 𝜕𝑣𝑖 𝜕 𝑝 + 𝜌𝑈 𝜕 2 𝑣𝑖
𝜌 + 𝑣𝑗 =− +𝜇
𝜕𝑡 𝜕𝑥𝑗 𝜕𝑥𝑖 𝜕𝑥𝑗2
• Eje X
𝑣Ԧ = 𝑣𝑥 𝑦, 𝑧, 𝑡 𝑖Ԧ
𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥 𝜕𝑣𝑥
𝜌 + 𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧 =
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕 𝑝 + 𝜌𝑈 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
=− +𝜇 2 + 2 +
𝜕𝑥 𝜕𝑥 𝜕𝑦 𝜕𝑧 2
𝜕𝑣𝑥 𝜕 𝑝 + 𝜌𝑈 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
𝜌 =− +𝜇 2
+
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 2

4.2 Movimiento laminar permanente 10

 Ecuaciones del movimiento
• Eje Y
𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑦 𝜕𝑣𝑦
𝜌 + 𝑣𝑥 + 𝑣𝑦 + 𝑣𝑧 =
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕 𝑝 + 𝜌𝑈 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑦 𝜕 2 𝑣𝑦
=− +𝜇 2
+ 2
+
𝜕𝑦 𝜕𝑥 𝜕𝑦 𝜕𝑧 2
𝜕 𝑝 + 𝜌𝑈
0=−
𝜕𝑦
• Eje Z
𝜕 𝑝 + 𝜌𝑈
0=−
𝜕𝑧
4.2 Movimiento laminar permanente 11
 Ecuaciones del movimiento
• Ecuaciones de Navier-Stokes
𝜕𝑣𝑥 𝜕 𝑝 + 𝜌𝑈 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
𝑥 𝜌 =− +𝜇 2 + 2
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧
𝜕 𝑝 + 𝜌𝑈
𝑦 0=−
𝜕𝑦
𝜕 𝑝 + 𝜌𝑈
𝑧 0=−
𝜕𝑧
• Sea Pr la presión reducida. 𝑃𝑟 = 𝑝 − 𝑝𝑎𝑡𝑚 + 𝜌𝑈 = 𝜌𝑔𝐻𝑝𝑖𝑒𝑧𝑜𝑚é𝑡𝑟𝑖𝑐𝑎
𝜕𝑣𝑥 𝜕𝑃𝑟 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
𝑥 𝜌 =− +𝜇 +
𝜕𝑡 𝜕𝑥 𝜕𝑦 2 𝜕𝑧 2
𝜕𝑃𝑟
𝑦 0=−
𝜕𝑦
𝜕𝑃𝑟
𝑧 0=−
𝜕𝑧

4.2 Movimiento laminar permanente 12

 Ecuaciones del movimiento
𝜕𝑃𝑟 𝜕𝑃𝑟
• Como =0y = 0, la presión reducida no depende de
𝜕𝑦 𝜕𝑧
la variable y ni de la variable z.
𝜕𝑣𝑥 𝜕𝑃𝑟 𝜕2 𝑣𝑥 𝜕2 𝑣𝑥
• Como 𝑣𝑥 𝑦, 𝑧, 𝑡 y 𝜌 = − +𝜇 + ,
𝜕𝑡 𝜕𝑥 𝜕𝑦 2 𝜕𝑧 2
𝜕𝑃𝑟
entonces =𝑓 𝑡 .
𝜕𝑥
• La presión reducida es lineal respecto a la variable x.
• Movimiento estacionario, 𝑣Ԧ = 𝑣𝑥 𝑦, 𝑧 𝑖Ԧ.
𝜕𝑣𝑥 𝜕𝑃𝑟 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
𝜌 =− +𝜇 2
+
𝜕𝑡 𝜕𝑥 𝜕𝑦 𝜕𝑧 2
𝜕𝑃𝑟 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
0=− +𝜇 2
+
𝜕𝑥 𝜕𝑦 𝜕𝑧 2

4.2 Movimiento laminar permanente 13

 Corriente de Couette
• Corriente entre 2 placas
paralelas infinitas. 𝑉0
• Sin gradiente de presión 𝑣𝑥 = 𝑦
𝜕𝑃𝑟
reducida. =0
ℎ
𝜕𝑥
• La placa superior se
desplaza a una velocidad
constante V.
𝜕 2 𝑣𝑥
0=
𝜕𝑦 2
• Condiciones de contorno:
– vx(0)=0
– vx(h)=V0

4.2 Movimiento laminar permanente 14

 Corriente de Couette
𝜕 2 𝑣𝑥 𝜕𝑣𝑥
2
=0⇒ = 𝑐1 ⇒ 𝑉0
𝜕𝑦 𝜕𝑦 𝑣𝑥 = 𝑦
𝑣𝑥 = 𝑐1 𝑦 + 𝑐2 ℎ
• En y=0, v=0 y en y=h, v=V0:
𝑣𝑥 0 = 𝑐1 0 + 𝑐2 = 0 ⇒
𝑐2 = 0

𝑣𝑥 ℎ = 𝑐1 ℎ = 𝑉0 ⇒
𝑉0
𝑐1 =
ℎ

4.2 Movimiento laminar permanente 15

 Corriente de Hagen-Poiseuille
• Corriente entre 2 placas 1 𝜕𝑃𝑟 2
paralelas infinitas. 𝑣𝑥 = 𝑦 − ℎ𝑦
2𝜇 𝜕𝑥
𝜕𝑃𝑟 𝜕 2 𝑣𝑥
0=− +𝜇
𝜕𝑥 𝜕𝑦 2
𝜕𝑃𝑟
• no depende de t (flujo
𝜕𝑥
estacionario), es
constante.
• Condiciones de contorno: 𝜕𝑃𝑟 𝑃𝑟2 − Pr1
– vx(0)=0 =
𝜕𝑥 L
𝑃𝑟2 −Pr1
– vx(h)=0 𝑣𝑥 = 2𝜇𝐿 𝑦 2 − ℎ𝑦

4.2 Movimiento laminar permanente 16

 Corriente de Hagen-Poiseuille
𝜕𝑃𝑟 𝜕 2 𝑣𝑥 1 𝜕𝑃𝑟 𝜕 2 𝑣𝑥 1 𝜕𝑃𝑟 2
0=− +𝜇 2 ⇒ = ⇒
𝜕𝑥 𝜕𝑦 𝜇 𝜕𝑥 𝜕𝑦 2 𝑣𝑥 = 𝑦 − ℎ𝑦
𝜕𝑣𝑥 1 𝜕𝑃𝑟
= y + c1 ⇒
2𝜇 𝜕𝑥
𝜕𝑦 𝜇 𝜕𝑥
1 𝜕𝑃𝑟 y 2
𝑣𝑥 = + c1 y + c2
𝜇 𝜕𝑥 2
• En y=0 y en y=h, v=0.
1 𝜕𝑃𝑟 2
𝑣𝑥 0 = 0 + c1 0 + c2 = 0 ⇒
2𝜇 𝜕𝑥
𝑐2 = 0

1 𝜕𝑃𝑟 2
𝑣𝑥 ℎ = ℎ + 𝑐1 ℎ ⇒
2𝜇 𝜕𝑥 𝜕𝑃𝑟 𝑃𝑟2 − Pr1
1 𝜕𝑃𝑟 =
𝑐1 = − ℎ 𝜕𝑥 L
𝑃𝑟2 −Pr1
2𝜇 𝜕𝑥 𝑣𝑥 = 2𝜇𝐿 𝑦 2 − ℎ𝑦

4.2 Movimiento laminar permanente 17

 Corriente de Hagen-Poiseuille
1 𝜕𝑃𝑟 2
• Caudal: 𝑣𝑥 = 𝑦 − ℎ𝑦
ℎ ℎ
1 𝜕𝑃𝑟 2
2𝜇 𝜕𝑥
𝑞 = න 𝑣𝑥 𝑑𝑦 = න 𝑦 − ℎ𝑦 𝑑𝑦
0 0 2𝜇 𝜕𝑥
ℎ 3 3
1 𝜕𝑃𝑟 1 𝜕𝑃𝑟 ℎ ℎ
= න 𝑦 2 − ℎ𝑦 𝑑𝑦 = −
2𝜇 𝜕𝑥 0 2𝜇 𝜕𝑥 3 2
1 𝜕𝑃𝑟 3
=− ℎ
12𝜇 𝜕𝑥
• Velocidad máxima:
𝜕𝑣𝑥 1 𝜕𝑃𝑟 ℎ
=0⇒ 2𝑦 − ℎ = 0 ⇒ 𝑦 = ⇒
𝜕𝑦 2𝜇 𝜕𝑥 2
1 𝜕𝑃𝑟 ℎ 2 ℎ 2 1 𝜕𝑃𝑟 2
𝑣𝑥 𝑚𝑎𝑥 = − =− ℎ
2𝜇 𝜕𝑥 4 2 8𝜇 𝜕𝑥
• Velocidad media:
𝑞 1 𝜕𝑃𝑟 2 2
𝑣𝑥 𝑚𝑒𝑑 = = − ℎ = 𝑣𝑥 𝑚𝑎𝑥
ℎ 12𝜇 𝜕𝑥 3

4.2 Movimiento laminar permanente 18

 Corriente de Hagen-Poiseuille
• El esfuerzo cortante tiene un 1 𝜕𝑃𝑟 2
ℎ
valor nulo en 𝑦 = y el 𝑣𝑥 = 𝑦 − ℎ𝑦
2
módulo toma el valor máximo 2𝜇 𝜕𝑥
en 𝑦 = 0 y en 𝑦 = ℎ.

∗
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑥
𝜏𝑥𝑦 =𝜇 + =𝜇
𝜕𝑦 𝜕𝑥 𝜕𝑦
1 𝜕𝑃𝑟
= 2𝑦 − ℎ
2 𝜕𝑥
∗
1 𝜕𝑃𝑟 4𝜇𝑣𝑥 𝑚𝑎𝑥
𝜏𝑥𝑦 0 = − ℎ=
2 𝜕𝑥 ℎ
∗
1 𝜕𝑃𝑟 4𝜇𝑣𝑥 𝑚𝑎𝑥 1 𝜕𝑃𝑟 2 1 1 𝜕𝑃𝑟
𝜏𝑥𝑦 h = ℎ=− 𝑣𝑥 𝑚𝑎𝑥 = −
8𝜇 𝜕𝑥
ℎ =− ℎ
4𝜇 2 𝜕𝑥
ℎ
2 𝜕𝑥 ℎ

4.2 Movimiento laminar permanente 19

 F. Laminar en conductos de sección
circular
• En coordenadas
cilíndricas:
𝜕𝑃𝑟 𝜕 2 𝑣𝑥 1 𝜕𝑣𝑥
=𝜇 2
+ ⇒
𝜕𝑥 𝜕𝑟 𝑟 𝜕𝑟
𝜕𝑃𝑟 1 𝜕 𝜕𝑣𝑥
=𝜇 𝑟
𝜕𝑥 𝑟 𝜕𝑟 𝜕𝑟
• Condiciones de contorno:
– vx(R)=0
– vx(0) es finita

4.2 Movimiento laminar permanente 20

 F. Laminar en conductos de sección
circular
• Solución de la ecuación diferencial:
𝜕𝑃𝑟 1𝜕 𝜕𝑣𝑥 1 𝜕𝑃𝑟 𝜕𝑣𝑥
=𝜇 𝑟 ⇒න 𝑟𝑑𝑟 = න 𝑑 𝑟 ⇒
𝜕𝑥 𝑟 𝜕𝑟 𝜕𝑟 𝜇 𝜕𝑥 𝜕𝑟
1
2 ×
1 𝜕𝑃𝑟 𝑟 𝜕𝑣𝑥 𝜕𝑣𝑥 1 𝜕𝑃𝑟 𝑟 𝐶1
𝑟
+ 𝐶1 = 𝑟 = +
𝜇 𝜕𝑥 2 𝜕𝑟 𝜕𝑟 𝜇 𝜕𝑥 2 𝑟
• En r=0, debe tomar un valor finito ⇒ 𝐶1 = 0
1 𝜕𝑃𝑟 𝑟 1 𝜕𝑃𝑟 𝑟 2
𝑣𝑥 = න 𝑑𝑟 = + 𝐶2
𝜇 𝜕𝑥 2 𝜇 𝜕𝑥 4
1 𝜕𝑃𝑟 𝑅2
• En r=R, vx=0 ⇒ 𝐶2 = −
𝜇 𝜕𝑥 4
1 𝜕𝑃𝑟 2
𝑣𝑥 = 𝑟 − 𝑅2
4𝜇 𝜕𝑥

4.2 Movimiento laminar permanente 21

 F. Laminar en conductos de sección
circular
• Caudal
1 𝜕𝑃𝑟 2
𝑄 = න 𝑣𝑥 𝑑𝑆 = න 𝑣𝑥 2𝜋𝑟𝑑𝑟 = න 𝑟 − 𝑅2 2𝜋𝑟𝑑𝑟
4𝜇 𝜕𝑥
𝜋 𝜕𝑃𝑟 𝑅 3 2
𝜋 𝜕𝑃 𝑟 R 4 R 4 𝜋 𝜕𝑃𝑟 4
= න 𝑟 − 𝑅 𝑟 𝑑𝑟 = − =− R
2𝜇 𝜕𝑥 0 2𝜇 𝜕𝑥 4 2 8𝜇 𝜕𝑥
• Velocidad máxima.
1 𝜕𝑃𝑟 2
r = 0 ⇒ 𝑣𝑥 𝑚𝑎𝑥 = − 𝑅
4𝜇 𝜕𝑥
• Velocidad media.
𝜋 𝜕𝑃𝑟 4
𝑄 − R 1 𝜕𝑃𝑟 2 𝑣𝑥 𝑚𝑎𝑥
8𝜇 𝜕𝑥
𝑣𝑥 𝑚𝑒𝑑 = = 2 =− 𝑅 =
𝑆 𝜋𝑅 8𝜇 𝜕𝑥 2

4.2 Movimiento laminar permanente 22

 F. Laminar en conductos de sección
circular
• El esfuerzo cortante tiene valor 1 𝜕𝑃𝑟 2
nulo en r=0 y su valor máximo en 𝑣𝑥 = 𝑟 − 𝑅2
4𝜇 𝜕𝑥
r=R.
∗
𝜕𝑣𝑟 𝜕𝑣𝑥 𝜕𝑣𝑥
𝜏𝑟𝑥 = 𝜇 + =𝜇
𝜕x 𝜕𝑟 𝜕𝑟
1 𝜕𝑃𝑟 1 𝜕𝑃𝑟
= 2r = 𝑟⇒
4 𝜕𝑥 2 𝜕𝑥
∗
1 𝜕𝑃𝑟
𝜏𝑤 = 𝜏𝑝 = 𝜏𝑟𝑥 𝑅 = R
2 𝜕𝑥

2𝜇𝑣𝑥 𝑚𝑎𝑥 4𝜇𝑣𝑥 𝑚𝑒𝑑

𝜏𝑝 = − =−
R 𝑅
4𝜇𝑄
=− 3
𝜋𝑅 1 𝜕𝑃𝑟 2 R 1 𝜕𝑃𝑟
𝑣𝑥 𝑚𝑎𝑥 = − 𝑅 =− R
4𝜇 𝜕𝑥 2𝜇 2 𝜕𝑥

4.2 Movimiento laminar permanente 23

 Pérdida de carga
• Ec. Bernouilli generalizado entre 1
y2
𝑝1 𝑣12 𝑝2 𝑣22
+ + 𝑧1 = + + 𝑧2 + ℎ𝑓
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔
• Ec. Continuidad
𝑄1 = 𝑄2 ⇒ 𝑣1 𝐴1 = 𝑣2 𝐴2 ⇒ 𝑣1 = 𝑣2
𝑝1 𝑝2
+ 𝑧1 = + 𝑧2 + ℎ𝑓 ⇒
𝜌𝑔 𝜌𝑔
∆𝑝
Δ𝐻 = + ∆𝑧 = ℎ𝑓
𝜌𝑔
La pérdida de carga por fricción en una
tubería de sección constante es igual a
la variación de altura piezométrica.
Fuente: Frank M. White. Fluid Mechanics.

4.3 Pérdida de carga 24

 Pérdida de carga
• Para un volumen diferencial
𝑝 + ∆𝑝 𝜋𝑅2 + 𝜌 𝜋𝑅2 ∆𝐿 𝑔 𝑠𝑒𝑛𝜙
− 𝜏𝑤 2𝜋𝑅 ∆𝐿 − 𝑝𝜋𝑅2 = 0 ⇒

∆𝑝 𝜋𝑅2 + 𝜌𝑔 𝜋𝑅2 ∆𝐿 𝑠𝑒𝑛𝜙

= 𝜏𝑤 2𝜋𝑅 ∆𝐿 ⇒
• Las pérdidas de carga son proporcionales
al esfuerzo cortante en la pared y la
longitud del conducto e inversamente
proporcional al radio.
∆𝑝 2𝜏𝑤 ∆𝐿
+ ∆𝑧 = ℎ𝑓 =
𝜌𝑔 𝜌𝑔 𝑅

4.3 Pérdida de carga 25

 Ecuación de Darcy-Weisbach
𝜏𝑤 = 𝐹 𝜌, 𝜇, 𝑉, 𝐷, 𝜖 𝑓𝜌𝑉 2
• Mediante análisis 2𝜏𝑤 ∆𝐿 2 8 ∆𝐿
dimensional establecemos ℎ𝑓 = =
𝜌𝑔 𝑅 𝜌𝑔 𝑅
𝜏𝑤 𝜖
= 𝐹 𝑅𝑒, 𝑓 ∆L 𝑉 2 ∆𝐿 𝑉 2
𝜌𝑉 2 𝑑 = =𝑓
4 𝐷 𝑔 𝐷 2𝑔
• Darcy definió el factor de 2
fricción f como: • Ecuación de Darcy-
8𝜏𝑤 𝜖 Weisbach
𝑓= 2
= 𝐹 𝑅𝑒,
𝜌𝑉 𝑑 𝐿 𝑉2
• Henry Darcy (1803-1858) ℎ𝑓 = 𝑓
𝐷 2𝑔
introdujo el efecto la
rugosidad en la pérdida de • Propuesta en 1850 Julius
carga en conductos Weisbach.

4.3 Pérdida de carga 26

 Factor de fricción en flujo laminar
• Flujo en conducto de sección circular:
1 𝜕𝑃𝑟 2
𝑣𝑥 = 𝑟 − 𝑅2
4𝜇 𝜕𝑥
1 𝜕𝑃𝑟 2 𝑣𝑥 𝑚𝑎𝑥
𝑣𝑥 𝑚𝑒𝑑 = 𝑅 =
8𝜇 𝜕𝑥 2
𝜕𝑣𝑥 1 𝜕𝑃𝑟 1 𝜕𝑃𝑟
𝜏𝑤 = 𝜇 = 𝑟 = 𝑅
𝜕𝑟 𝑟=𝑅 2 𝜕𝑥 𝑟=𝑅 2 𝜕𝑥
4𝜇𝑣𝑥 𝑚𝑒𝑑
=
𝑅
4𝜇𝑣𝑥 𝑚𝑒𝑑 8𝜇𝑣𝑥 𝑚𝑒𝑑 Fuente: Frank M. White. Fluid Mechanics.
8𝜏𝑤 8 8
𝑅 𝐷
𝑓= 2 = =
𝜌𝑣𝑥 𝑚𝑒𝑑 𝜌𝑣𝑥2 𝑚𝑒𝑑 𝜌𝑣𝑥2 𝑚𝑒𝑑 64
=
64𝜇
=
64 𝑓𝑙𝑎𝑚𝑖𝑛𝑎𝑟 =
𝜌𝑣𝑥 𝑚𝑒𝑑 𝐷 𝑅𝑒 𝑅𝑒
4.3 Pérdida de carga 27
 Factor de fricción en flujo laminar
• Flujo bidimensional plano:
1 𝜕𝑃𝑟 2
𝑣𝑥 = 𝑦 − ℎ𝑦
2𝜇 𝜕𝑥
1 𝜕𝑃𝑟 2 2
𝑣𝑥 𝑚𝑒𝑑 = ℎ = 𝑣𝑥 𝑚𝑎𝑥
12𝜇 𝜕𝑥 3
𝜕𝑣𝑥 1 𝜕𝑃𝑟
𝜏𝑤 = 𝜇 = (2𝑦 − ℎ)
𝜕𝑦 𝑦=ℎ 2 𝜕𝑥 𝑦=ℎ
1 𝜕𝑃𝑟 6𝜇𝑣𝑥 𝑚𝑒𝑑
= ℎ =
2 𝜕𝑥 ℎ

8𝜏𝑤 8
6𝜇𝑣𝑥 𝑚𝑒𝑑
ℎ
8
6𝜇𝑣𝑥 𝑚𝑒𝑑
ℎ
48
𝑓= 2 = = 𝑓𝑙𝑎𝑚𝑖𝑛𝑎𝑟 =
𝜌𝑣𝑥 𝑚𝑒𝑑 𝜌𝑣𝑥2 𝑚𝑒𝑑 𝜌𝑣𝑥2 𝑚𝑒𝑑 𝑅𝑒
48𝜇 48
= =
𝜌𝑣𝑥 𝑚𝑒𝑑 ℎ 𝑅𝑒

4.3 Pérdida de carga 28

Fuente: Frank M. White. Fluid Mechanics.

4.3 Pérdida de carga 29

 Efecto de longitud finita del tubo
• Las soluciones
obtenidas son válidas
para conductos
infinitamente largos. Fuente: Antonio Crespo Martínez. Mecánica de Fluidos.

Condición de
unidireccionalidad.
• En caso de conductos
finitos, aparece
velocidades 𝑣𝑦 = 𝑂 𝑣𝑇 , 𝑣𝑧 = 𝑂(𝑣𝑇 )
transversales de orden
de magnitud vT.

4.4 Efecto de longitud finita del tubo 30

 Efecto de longitud finita del tubo
𝑣𝑥 = 𝑂 𝑢 , 𝑣𝑦 = 𝑂 𝑣𝑇 , 𝑣𝑧 = 𝑂 𝑣𝑇 • Componente en eje X del término
convectivo
𝑥 = 𝑂 𝐿 ,𝑦 = 𝑂 𝐷 ,𝑧 = 𝑂 𝐷
𝑢2 𝑣𝑇
𝜌 𝑣Ԧ ⋅ 𝛻 𝑢 = 𝑂 𝜌 + 𝑂 𝜌𝑢
𝜕𝑣𝑥 𝜕𝑣𝑦 𝜕𝑣𝑧 𝐿 𝐷
𝛻 ⋅ 𝑣Ԧ = + + =0 𝑢2
𝜕𝑥 𝜕𝑦 𝜕𝑧 =𝑂 𝜌
𝑢 𝑣𝑇 𝐿
𝛻 ⋅ 𝑣Ԧ = 𝑂 +𝑂 =0
𝐿 𝐷 • Término viscoso
𝜕2𝑢 𝜕2𝑢 𝑢
𝑣𝑇 𝑢 𝑣𝑇 𝐷 𝜇 + = 𝑂 𝜇
𝑂 =𝑂 ⇒𝑂 =𝑂 𝜕𝑦 2 𝜕𝑧 2 𝐷2
𝐷 𝐿 𝑢 𝐿
• Para despreciar el término
convectivo, se debe cumplir
• La velocidad transversal es tanto
menor cuanto mayor sea la 𝑢2
𝜌
longitud del conducto. 𝐿 ≪ 1 ⇒ 𝜌𝑢𝐷 𝐷 = 𝑅𝑒 𝐷 ≪ 1
𝑢 𝜇 𝐿 𝐿
𝜇 2
𝐷

4.4 Efecto de longitud finita del tubo 31

 Efecto de longitud finita del tubo
𝜌𝑢𝐷 𝐷 𝑢𝐷2
≪1⇒ ≪𝐿
𝜇 𝐿 𝜈
𝑢𝐷 2
𝐿𝑒 ≈
𝜈
• La longitud de entrada es la distancia
en la entrada de una conducción, en
la que el efecto de las paredes se
difunde debido a los efectos viscosos
hasta alcanzar el flujo desarrollado.
• Longitudes de entradas
características:
𝐿𝑒
– ≈ 0,06𝑅𝑒 para flujo laminar.
𝑑
1
𝐿𝑒
– 𝑑
≈ 4,4 ⋅ 𝑅𝑒 para flujo
6

turbulento. Fuente: Frank M. White. Fluid Mechanics.

4.4 Efecto de longitud finita del tubo 32

 Efecto de longitud finita en la pérdida
de carga
• En un conducto de sección circular, despreciando los efectos de la entrada:
𝐿
∆𝑝 + 𝜌𝑔∆𝑧 = 𝜌𝑔ℎ𝑓 = ∆𝑃𝑟 = 2𝜏𝑤
𝑅
4𝜇𝑣𝑥 𝑚𝑒𝑑
𝜏𝑤 =
𝑅
4𝜇𝑣𝑥 𝑚𝑒𝑑 𝐿 8𝜇𝑣𝑥 𝑚𝑒𝑑 𝐿
∆𝑃𝑟 = 2 =
𝑅 𝑅 𝑅2
• De acuerdo con los resultado experimentales, los efectos de entrada influyen de la siguiente
forma:
8𝜇𝑣𝑥 𝑚𝑒𝑑 𝐿 2 ∆𝑃𝑟 8𝜇𝑣𝑥 𝑚𝑒𝑑 𝐿 0,582𝜌𝑣𝑥2 𝑚𝑒𝑑
∆𝑃𝑟 = + 0,582𝜌𝑣𝑥 𝑚𝑒𝑑 ⇒ = +
𝑅2 1 2 1 2 1 2
𝜌𝑣𝑥 𝑚𝑒𝑑 𝜌𝑣𝑥 𝑚𝑒𝑑 𝑅 2 𝜌𝑣
2 2 2 𝑥 𝑚𝑒𝑑
64𝜇𝐿 64 𝐿
= + 1,165 = + 1,165
𝜌𝑣𝑥 𝑚𝑒𝑑 𝐷 2 𝑅𝑒 𝐷
• Para que se cumpla la hipótesis de unidireccionalidad debe cumplirse:
64 𝐿 1 𝜌𝑣𝑥 𝑚𝑒𝑑 𝑅 2 1 𝜌𝑢𝐷 𝐷
≫ 1,165 ⇒ ≪ 𝐿, Re D ≪ L; menos restrictiva que: ≪1
𝑅𝑒 𝐷 13 μ 52 𝜇 L

4.4 Efecto de longitud finita del tubo 33

 Flujos cuasi-estacionarios
• Para considerar un flujo cuasi estacionario, debe
cumplirse:
𝜕𝑣𝑥 𝜕 2 𝑣𝑥 𝜕 2 𝑣𝑥
𝜌 ≪𝜇 2
+
𝜕𝑡 𝜕𝑦 𝜕𝑧 2
• Si definimos tc como el tiempo característico de la
variación de propiedades, en órdenes de
magnitud podríamos expresar la condición de
flujo cuasi-estacionario como:
𝜌𝑢 𝜇𝑢 𝜌 𝜇 𝜌𝐷2 𝐷2
𝑂 ≪ 𝑂 2 ⇒ ≪ 2 ⇒ 𝑡𝑐 ≫ =
𝑡𝑐 𝐷 𝑡𝑐 𝐷 𝜇 𝜈
4.4 Efecto de longitud finita del tubo 34
Flujos cuasi-estacionarios: Descarga de
un depósito
¿Caudal que descarga el depósito?
𝜕𝑃𝑟 𝑝𝑠 − 𝑝𝑒 + 𝜌𝑔𝐻
= Depósito cilíndrico A=20cm2
𝜕𝑥 𝐿
𝑝𝑎𝑡𝑚 − 𝑝𝑎𝑡𝑚 + 𝜌𝑔ℎ0 + 𝜌𝑔𝐻
=
𝐿
𝜌𝑔 𝐻 + ℎ0
=−
𝐿

𝜋𝑅4 𝜕𝑃𝑟 𝜋𝑅 4 𝜌𝑔 𝐻 + ℎ0
𝑄=− = ⇒
8𝜇 𝜕𝑥 8𝜇 𝐿
10−6 𝑚3
𝑄 = 1,54 ⋅
𝑠

4Q −2 𝑚/𝑠
𝑣med = = 49 ⋅ 10
𝜋𝐷2

4.4 Efecto de longitud finita del tubo 35

Flujos cuasi-estacionarios: Descarga de
un depósito
Comprobamos las hipótesis de flujo Integrado con la condición t=0, h=h0:
laminar y unidireccionalidad −
𝜋𝑅 4 𝜌𝑔
𝑡
𝜌𝑣𝑚𝑒𝑑 𝐷 ℎ = ℎ0 + 𝐻 𝑒 8𝜇𝐿𝐴 −𝐻
𝑅𝑒 = = 980 < 2300 𝑡
𝜇 −
= ℎ0 + 𝐻 𝑒 𝑡𝑐 − 𝐻
⇒ 𝐹𝑙𝑢𝑗𝑜 𝐿𝑎𝑚𝑖𝑛𝑎𝑟
Donde el tiempo característico de
1 𝜌𝑣𝑥 𝑚𝑒𝑑 𝑅 𝑅 descarga sería:
= 0,075 ≪ 1
13 μ 𝐿 8𝜇𝐿𝐴
⇒ 𝑈𝑛𝑖𝑑𝑖𝑟𝑒𝑐𝑐𝑖𝑜𝑛𝑎𝑙 𝑡𝑐 = = 520𝑠
𝜋𝑅4 𝜌𝑔

Aplicando continuidad Cumple la condición de flujo cuasi-

𝑑ℎ 𝜋𝑅4 𝜌𝑔 𝐻 + ℎ estacionario
−𝐴 =
𝑑𝑡 8𝜇 𝐿 𝜌𝐷 2
520𝑠 ≫ = 4𝑠
𝜇

4.4 Efecto de longitud finita del tubo 36

Transición
• Re<Red, crit → Flujo
laminar
• Re>Red, crit → Flujo
turbulento
• Re≈Red, crit → Zona de
Fuente: Frank M. White. Fluid Mechanics.
transición
• En tubo comercial el
valor aceptado es:
𝑅𝑒𝑑,𝑐𝑟𝑖𝑡 ≈ 2300

4.5 Estabilidad de la corriente laminar 37

 Estabilidad de la corriente laminar
• Experimentalmente se
observa que el flujo laminar
existe cuando Re<2300:
– V es lo suficientemente baja
– D es lo suficientemente
pequeño
– ρ es lo suficientemente baja
– μ es lo suficientemente
elevada
• Las perturbaciones que
aparecen sobre el flujo son
amortiguadas rápidamente.
• Transición a la turbulencia:
2000 < 𝑅𝑒 < 4000

4.5 Estabilidad de la corriente laminar 38

Estabilidad de la corriente laminar

4.5 Estabilidad de la corriente laminar 39