Número de Reynolds en Flujos No Compresibles
Número de Reynolds en Flujos No Compresibles
Flujo de fluidos no
compresibles en tuberías
y canales de conducción
En los procesos industriales, los ingenieros químicos tratan de manera frecuente con el
flujo de fluidos a través de tuberías, tubos y canales de conducción con sección transversal
no circular. Normalmente las tuberías se llenan con fluidos en movimiento, pero algunos
problemas involucran el flujo en tuberías parcialmente llenas, en capas descendentes sobre
superficies inclinadas o verticales, a través de lechos de sólidos o en tanques agitados.
En este capítulo se estudia el flujo estacionario de fluidos no compresibles a través de
tuberías cerradas y canales de conducción. En el capítulo 4 se trató el flujo en capas;
en el capítulo 6 se cubrirá el flujo de fluidos compresibles; en el capítulo 7, el flujo a
través de lechos de sólidos, y en el capítulo 9 el flujo en tanques agitados.
r
rw
FIGURA 5.1
y τ Elemento de fluido en el flujo estacionario a
dL través de una tubería.
Sa � Sb � �r 2 pa � p pa Sa � �r 2 p � �
pb Sb � �r 2 � p � dp�
El esfuerzo cortante Fs que actúa sobre el borde del elemento es igual al producto de
la tensión de corte por el área del cilindro, es decir (2πr dL)τ. [Fuerza FS igual –FW en la
ecuación (4.52).] Puesto que el canal de conducción es horizontal, Fg es igual a cero. Al
sustituir estos términos en la ecuación (4.52) se obtiene
� F � �r 2
p � �r 2 � p � dp� � �2�r dL�� � 0
FIGURA 5.2
rw Variación del esfuerzo cortante en una
r tubería.
cL de la tubería
Pared de la tubería
τw
0
Esfuerzo cortante, τ
de energía cinética se cancelan. Por otra parte, el único tipo de fricción es la fricción de
superficie entre la corriente de fluido y la pared, y la única caída de presión es la que
resulta de la fricción en la superficie. La fricción en la superficie se representa por hfs y
la caída de presión por ps. Entonces la ecuación (4.71) se convierte en
pa pa � �ps
� � h fs
� �
�ps
o � h fs (5.4)
�
Factor de fricción
Otro parámetro frecuente, especialmente útil en el estudio del flujo turbulento, es el fac-
tor de fricción de Fanning,† que se representa por f y se define como la relación entre el
–
esfuerzo cortante y el producto de la carga de velocidad V 2/2 por la densidad:
�w 2�
f� 2 � w2 (5.6)
�V /2 �V
†
Otro factor de fricción muy empleado en la literatura de la mecánica de fluidos es el llamado factor de
fricción de Blasius o de Darcy, que es 4f (cuatro veces el valor del factor de fricción de Fanning f).
2 �w �p L V2
h fs � L � s � 4f (5.7)
� rw � D 2
de la cual
�ps D (5.8)
f �
2 L�V 2
y �ps 2 f�V 2
� (5.9)
L D
La ecuación (5.9) es la ecuación que se emplea generalmente para calcular la pérdida
de fricción superficial en la tubería recta. El subíndice s se utiliza en ps y hfs para indicar
el hecho de que en las ecuaciones (5.7) a (5.9), las magnitudes mencionadas están aso-
ciadas con el factor de fricción de Fanning, y que corresponden solamente a la fricción
sobre la superficie. Si están presentes otros términos en la ecuación de Bernoulli, o si
existe fricción de forma, el valor de pa – pb difiere de ps . Si ocurre la separación de
capa límite, hf es mayor que hfs.
� D2 /4 D
rH � �
�D 4
El diámetro equivalente es 4rH, o simplemente D.
Un caso especial importante es el espacio anular entre dos tuberías concéntricas.
Aquí el radio hidráulico es
� Do2 /4 � � Di2 /4 Do � Di
rH � � (5.11)
� Di � � Do 4
donde Di y Do son los diámetros de entrada y salida de los espacios anulares, respec-
tivamente. El diámetro equivalente de un espacio anular es por lo tanto la diferencia
Las ecuaciones (5.1) a (5.9) se aplican tanto al flujo laminar como turbulento, siempre que
el fluido sea no compresible y el flujo estacionario y completamente desarrollado. Debido
a la simplicidad de la relación del esfuerzo cortante-viscosidad para el flujo laminar, se
pueden deducir fácilmente nuevas ecuaciones para este caso.
dS � 2� r dr (5.12)
u �w r
� 0
du � � � rdr
rw � rw
�
u�
2r �
w
w
�r � r �2
w
2
(5.15)
�w � w rw
� �r �
rw
2
V� w � r 2 r dr � (5.18)
rw3 � 0 4�
Distancia relativa desde el centro de la tubería
1.0
0.8
0.6
0.4
Flujo laminar
0.2
0
0.2
Flujo turbulento
0.4
0.6
0.8
1.0
0 0.10 0.20 0.30 0.40 0.50 0.60 0.70 0.80 0.90 1.00
Fracción de la velocidad máxima
FIGURA 5.3
Distribución de velocidad en una tubería, con flujo totalmente desarrollado de un fluido newto-
niano, para flujo laminar y turbulento hasta Re = 10 000.
Ecuación de Hagen-Poiseuille
Para efectuar cálculos prácticos, la ecuación (5.18) se suele transformar eliminando τw
en función de ps mediante la ecuación (5.7) y utilizando el diámetro de la tubería en
lugar del radio de la misma. El resultado es el siguiente:
�ps rw rw �p D2
V� � s
L 2 4 � 32 L �
Despejando ps se obtiene
32 LV �
�ps � (5.20)
D2
y puesto que ps = 4 τw /(DL),
8V �
�w � (5.21)
D
Sustituyendo de la ecuación (5.21) en la ecuación (5.6) se obtiene
16 � 16
f � � (5.22)
DV � Re
La ecuación (5.20) es la ecuación de Hagen-Poiseuille. Uno de sus usos es en la
determinación experimental de la viscosidad, midiendo la caída de presión y la velocidad
volumétrica de flujo a través de un tubo de longitud y diámetro conocidos. A partir de
–
la velocidad de flujo se calcula V por medio de la ecuación (4.12) y µ con la ecuación
(5.20). En la práctica es necesario hacer correcciones para tener en cuenta la energía
cinética y los efectos de entrada.
n�
� 3n� � 1� V n� (5.24)
�ps � 2 K � � � L
� n� � rwn ��1
2.8
2.6
Dilatante: n´ = 2
2.4
2.2
Newtoniano: n´ = 1
2.0
1.8
1.6
u
1.4
V̄
Seudoplástico:
1.2 n´ = 0.5
1.0
0.8
0.6
0.4
0.2
0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 FIGURA 5.4
r Perfiles de velocidad para flujo laminar de líquidos newto-
rw
nianos y no newtonianos.
Pared de la tubería
0 Velocidad, u
a)
Pared de la tubería
Flujo pistón
Eje de la tubería
rc
rw
τ0
Pared de la tubería
0 τw
Esfuerzo cortante, τ
b)
FIGURA 5.5
a) Perfil de velocidad y b) diagrama del esfuerzo cortante para el flujo de un plástico Bingham.
corte límite τ0. El diagrama del esfuerzo cortante se representa en la figura 5.5b. Para la
variación de la velocidad en el espacio anular comprendido entre la pared del tubo y el
flujo pistón, se aplica la ecuación siguiente:
1 �� � r� �
u� � rw � r �� w �1 � � � � 0 � (5.25)
K � 2 � rw � �
donde K es una constante. El límite entre el flujo pistón y el resto del fluido se obtiene
diferenciando la ecuación (5.25) e igualando a cero el gradiente de velocidad o, más
fácilmente, interpretando el valor en la figura 5.5b. El resultado es
�0
rc � r (5.26)
�w w
La velocidad en el núcleo central uc, que es la velocidad de flujo pistón, se determina
sustituyendo el valor de rc de la ecuación (5.26) en vez de r en la ecuación (5.25) y
reordenando, se obtiene
�0
uc �
2 Krc
�rw � rc �2 (5.27)
�psr 20 �� �� r ��2 1� � 2 r ��
u� ��1� �� �� � ln �� (5.28)
�� ��r0 �� ln�1/� � r0 ��
4� L �� ��
donde ro = radio de la pared exterior del espacio anular
κ = relación ri /ro
ri = radio de la pared interior del espacio anular
� Do � Di �V � (5.29)
Re �
�
1.40
φa
1.30
1.24
0.01 0.1 1.0
Di /Do
donde Di y Do son los diámetros interior y exterior del espacio anular, respectivamente.
De la ecuación (5.28) se puede calcular el factor de fricción como se hizo la deducción
de la ecuación (5.22). El flujo laminar en un espacio anular se relaciona con el número de
Reynolds por la ecuación5
16
f � �a (5.30)
Re
donde φa es la función de Di /DO mostrada en la figura 5.6. El valor de φa es 1 para secciones
transversales circulares y 1.5 para planos paralelos. Las ecuaciones para flujo laminar
a través de otras secciones transversales se proporcionan en varios textos.5 El flujo de
líquidos viscosos no newtonianos a través de canales de forma compleja, como en las
operaciones de inyección para moldeo inyección, se analiza por métodos de aproximación
numérica empleando la ecuación (4.22).
procedentes del fluido turbulento que circula en esta región. Muy cerca de la pared los
remolinos son infrecuentes, pero no existe ninguna región que esté completamente
libre de remolinos. Dentro de la subcapa viscosa solamente es importante el esfuerzo
cortante viscosa (o esfuerzo cortante), mientras que la difusión de remolino, es en todo
caso, menor.
La subcapa viscosa ocupa sólo una muy pequeña fracción de la sección transversal
total. No existe un límite nítido superior y es difícil de definir su espesor. Existe una capa
de transición inmediatamente adyacente a la subcapa viscosa, en la que hay un esfuerzo
cortante debido a la viscosidad y también un esfuerzo cortante debido a la difusión de
remolino. La capa de transición, que se llama a veces capa amortiguadora, también
es relativamente delgada. La masa global de la sección transversal de la corriente está
ocupada por flujo enteramente turbulento que recibe el nombre de núcleo turbulento.
En el núcleo turbulento el esfuerzo cortante viscoso es insignificante en comparación con
el de la viscosidad de remolino mismo.
f �w
u� � V � (5.31)
2 �
u
u� � (5.32)
u�
yu� � y
y� � � � w�
� �
rw � r � y (5.34)
Puesto que la subcapa viscosa es muy delgada, r ≈ rw, y sustituyendo –dy por dr, la
ecuación (5.14) queda de esta forma
du � w
� (5.35)
dy �
Sustituyendo los valores de u* de la ecuación (5.31), u+ de la ecuación (5.32) y y+ de la
ecuación (5.33) en la ecuación (5.35) se obtiene
du �
�1
dy �
Integrando con los límites inferiores u+ = y+ = 0, se obtiene para la distribución de velo-
cidad en la subcapa laminar,
u� � y� (5.36)
La figura 5.7 muestra una gráfica semilogarítmica de las ecuaciones (5.36), (5.37)
y (5.38). De acuerdo con las dos intersecciones de las tres líneas que representan las
ecuaciones, los intervalos que comprenden las ecuaciones son los siguientes:
25
Subcapa Capa
viscosa Núcleo turbulento
amortiguadora
20 .5
+ +5
y
In
.5
+ =2
u
15
u+
05
3.
10
−
y+
In
5
=
u+
5
u+ = y +
0
1 3 5 10 30 50 100 300 500 1 000
y+
FIGURA 5.7
Distribución universal de velocidad; flujo turbulento de un fluido newtoniano en un tubo liso.
Velocidad media. Para la línea central de la tubería, la ecuación (5.38) puede escri-
birse de la siguiente manera:
donde u+c y y+c son los valores de u+ y y+ en la línea central, respectivamente. Así como,
por las ecuaciones (5.32) y (5.33)
umáx
uc� � (5.40)
u*
rw u *
yc� � (5.41)
v
donde umáx es la velocidad máxima en la línea central y ν = µ /ρ.
Restando la ecuación (5.39) de la ecuación (5.38) resulta
y�
u � � uc� � 2.5 ln (5.42)
yc�
A partir de la ecuación (4.11), sustituyendo πr 2w por S y 2πr dr por ds, se obtiene para
–
la velocidad media V
2� rw
V�
� rw2 � ur dr
0
(5.43)
5v 2 yc� � y� � �
V�
rw2 u * �
0
�
� �
� 0.4uc � ln � � yc � y dy
� yc �
�
� (5.45)
La ecuación (5.47) proporciona valores que son demasiado altos en relación a las bajas
velocidades del fluido en las capas próximas a la pared, de tal forma que 2% del flujo
volumétrico total, no es tomado en cuenta adecuadamente.
1 � f�
� 2.5 ln � Re 8 � � 1.75 (5.50)
f /2 � �
La ecuación (5.50) concuerda bien con el experimento. En el intervalo de los nú-
meros de Reynolds 104 < Re < 106 que predicen los factores de fricción dentro de 2% de
aquellos indicados en la figura 5.10.
†
La sustitución de x = y+/y+c y el uso de las tablas estándar de integrales son suficientes para efectuar esta
–
integración. Además, para el límite inferior, x ln x = 0, y x2 ln x = 0 cuando x = 0. De la ecuación (5.31), V /u*
= 1/f/2.
las ecuaciones para estos factores de corrección se obtiene fácilmente integrando las
ecuaciones (4.50) y (4.70) y utilizando la ley logarítmica de velocidad. Las ecuaciones
que así se obtienen son
�
� � 1 � 0.78 f 15 � 15.9 f � (5.51)
� � 1 � 3.91 f (5.52)
Efecto de la rugosidad
Hasta ahora el tratamiento se ha restringido al caso de los tubos lisos. Se sabe desde
hace mucho tiempo que para un flujo turbulento, una tubería rugosa conduce a un factor
de fricción mayor que una tubería lisa para un determinado número de Reynolds. Si se
pule una tubería rugosa, se reduce el factor de fricción. Cuando un mayor pulido no da
lugar a una reducción adicional del factor de fricción para un número de Reynolds dado,
se dice que el tubo es hidrodinámicamente liso. La ecuación (5.50) se refiere a este tipo
de tubos.
La figura 5.9 muestra varios tipos idealizados de rugosidad. La altura de una unidad
simple de rugosidad se indica por medio de k y se denomina parámetro de rugosidad.
Con base en el análisis dimensional, f es tanto función de Re como de la rugosidad rela-
tiva k/D, donde D es el diámetro de la tubería. Para un tipo dado de rugosidad, tal como
el que se muestra en la figura 5.9a y b, es de esperar que se halle una curva distinta de f
versus Re para cada magnitud de la rugosidad relativa, al igual que para otros tipos, como
0.9
0.8
V̄/umáx
0.7
0.6
0.5
103 104 105 106
DV̄ ρ/µ
los que se muestran en la figura 5.9c y d, se encuentre una familia diferente de curvas
para cada tipo de rugosidad. Los experimentos realizados sobre tuberías donde se creó
la rugosidad artificialmente, han confirmado estas expectativas. También se encontró que
las tuberías comerciales nuevas, completamente limpias, parecen tener el mismo tipo de
rugosidad y que cada material de construcción tiene su propio parámetro característico
de rugosidad. En el capítulo 8 se discuten los tamaños estándar de tuberías y ductos, los
cuales se enlistan en los apartados 3 y 4.
k
k D
D
a) c)
k
k
D D
b) d)
FIGURA 5.9
Tipos de rugosidad.
Las tuberías usadas, sucias y corroídas suelen ser muy rugosas, y es muy probable
que el carácter de la rugosidad sea diferente del de una tubería limpia.
La rugosidad no ejerce un efecto apreciable sobre el factor de fricción para el flujo
laminar a menos que k sea tan grande que la medida del diámetro sea incierta.
0.125
f � 0.0014 � (5.54)
Re 0.32
Las demás líneas curvas en el intervalo de flujo turbulento representan los factores
de fricción para varios tipos de tuberías comerciales, cada uno de los cuales se caracteriza
por tener un valor diferente de k. En la figura se dan los parámetros para algunos metales
comunes. Por ejemplo, una tubería limpia de hierro o acero forjado tiene un valor de k
de 1.5 10–4 ft, sin importar el diámetro de la tubería. Las tuberías estiradas de cobre y
latón se consideran hidráulicamente lisas.
Para tuberías de acero y otro tipo de tuberías rugosas, el factor de fricción se vuelve
independiente del número de Reynolds para números de Reynolds mayores que 106. Una
ecuación empírica para esta región es
(5.55)
La figura 5.10 es útil para calcular hfs cuando se conoce el tamaño de la tubería y la
velocidad de flujo, pero no puede utilizarse directamente para determinar la velocidad
nar,
0.0125 0.02
lami
0.015
f�
0.0100
Flujo 16/Re
0.01
0.008
0.0080 Recr 0.006
Flujo turbulento
SECCIÓN II
0.0070 0.004
Región de
0.0060 transición 0.002
Rugosidad relativa, k/D
0.0050 0.001
0.0008
erí
as 0.0002
Tubo rolado Lisa lisa
0.0030 Hierro forjado, acero 0.00015 s 0.0001
Hierro galvanizado 0.0005 0.000,05
0.0025 Hierro fundido 0.00085 k/D � 0.000005
Concreto 0.01–0.001
k/D � 0.000001
0.0020 0.000,01
1 2 3 4 56 8 1 2 3 4 56 8 1 2 3 4 56 8 1 2 3 4 56 8 1 2 3 4 56 8 1
� 103 � 104 � 105 � 106 � 107
Número de Reynolds, Re
FIGURA 5.10
Gráfica del factor de fricción para tuberías circulares (Adaptado de L. W. Moody, “Friction Factors for Pipe Flow”, Trans. ASME 66: 672 [1944]).
14/12/06 22:26:00
CAPÍTULO 5 Flujo de fluidos no compresibles en tuberías y canales de conducción 123
de flujo para una caída de presión dada, en virtud de que no se conoce Re mientras no
–
esté determinado el valor de V . Sin embargo, ya que f varía sólo ligeramente con Re
para el flujo turbulento, una solución por prueba y error converge rápidamente, como se
muestra en el ejemplo 5.1.
EJEMPLO 5.1 El agua fluye a 50 °F a través de una tubería larga, plástica y horizontal, de
3 in. de diámetro interior, a una velocidad de 8 ft/s. a) Calcular la caída de presión en lbf /in2
por 100 ft de tubería. b) Si la caída de presión se limita a 2 lbf /in.2 por 100 ft de tubería, ¿cuál
es la velocidad máxima permisible del agua?
Solución
a) Utilice la ecuación (5.9), la cual está en unidades fps
�ps 2 f�V 2
� (5.56)
L Dgc
Las propiedades del agua a 50 °F, del apéndice 6, son ρ = 62.42 lb/ft3 y µ = 1.310 cP. Además
–
D = 12
–3 = 0.25 ft, V = 8 ft/s, y L = 100 ft.
El número de Reynolds es
0.25 � 8 � 62.42
Re � � 1.4 � 10 5
1.310 � 6.72 � 10 �4
De la figura 5.10, f = 0.0041. Al sustituir en la ecuación (5.56) se obtiene
Esto da un número de Reynolds de 1.17 105, de modo que si se ajusta a la ecuación (5.53).
(5.57)
A partir de esta ecuación se define el número de Reynolds Ren para fluidos no newto-
nianos, con la suposición de que para flujo laminar
16
f � (5.58)
Re n
La figura 5.11 muestra que para fluidos seudoplásticos (n < 1), el flujo laminar
persiste para números de Reynolds más altos que con fluidos newtonianos.
f = 16/Ren
Diagrama del factor
de fricción no newtoniano
Factor de fricción de Fanning, f
n´=
0.01
9 2.0
8
7 1.4
6
1.0
5
0.8
4
0.6
3
0.4
2 0.3
Regiones experimentales 0.2
Regiones extrapoladas 0.0
0.001
1 000 10 000 100 000
Número de Reynolds modificado, Ren
FIGURA 5.11
Gráfica del factor de fricción para un fluido que sigue la ley de potencia. (Según D.W. Dodge y
A.B. Metzner.2)
0.01
Agua
en flu
0.008 jo tu
rbule
nto
0.006
1 ppm
f 0.004 5 ppm
0.003 10 ppm
Flu
jo
la
0.002
m
in
ar
50 ppm
FIGURA 5.12
Factores de fricción para el flujo turbu-
0.001 lento de disoluciones diluidas de óxido
2 × 103 3 4 6 8 104 2 3 4 6 8 105 de polietileno (PM 106). [Según R.W.
Re Patterson y F.H. Abernathy.7]
Flujo no isotérmico
Cuando los líquidos fluyen a través de un tubo liso que se calienta o enfría, el factor de
fricción aumenta o disminuye debido al efecto de la temperatura sobre la viscosidad del
líquido cercano a la pared. Cuando el líquido se enfría, la viscosidad más elevada cerca
de la pared, aumenta el gradiente de velocidad, lo que conduce a un incremento de la
fricción en la pared. El efecto es especialmente pronunciado cuando la viscosidad varía
fuertemente con la temperatura. Para los gases, debido a que la viscosidad aumenta con
la temperatura, lo contrario es cierto. En la práctica de ingeniería ordinaria se justifica,
tanto para gases como para líquidos, el empleo del siguiente método empírico.
Las ecuaciones (5.61) y (5.62) están basadas en datos para valores de µ/µw compren-
didos entre 0.1 y 10, y no deberán utilizarse fuera de estos límites.
En tuberías comerciales rugosas, más que en tubos lisos, el flujo está usualmente
en “turbulencia completa”, donde el factor de fricción es independiente del número de
Reynolds. Bajo estas circunstancias puede ser que f no sea correcto para el efecto de la
temperatura en la pared.10a
Q
� c� � �T � �ps
V
Por lo tanto,
�ps
�T � (5.63)
c� �
Como se mencionó antes, la relación para flujo turbulento en tuberías puede aplicarse a
canales no circulares sustituyendo el diámetro equivalente Deq (o 4 veces el radio hidráu-
lico rH) por el diámetro D en las ecuaciones relevantes.
V a2
h fe � Ke (5.64)
2
donde Ke es un factor de proporcionalidad llamado coeficiente de pérdida por expansión
–
y V a es la velocidad media en la parte estrecha del conducto. En este caso, es posible
calcular Ke teóricamente con resultados satisfactorios. Para hacer el cálculo se utiliza
la ecuación de continuidad (4.13), la ecuación del balance de momento para flujo es-
tacionario (4.51) y la ecuación de Bernoulli (4.71). Considere un volumen de control
definido por las secciones AA y BB y la superficie interna del conducto expandido que
existe en estas dos secciones, tal como se muestra en la figura 5.13. Debido a que en
la tubería horizontal no existen fuerzas de gravedad y a que la pared es relativamente
pequeña, la fricción en ella es insignificante, de manera que prácticamente no hay gra-
diente de velocidad en la pared comprendida entre las dos secciones. Por lo tanto, las
únicas fuerzas que actúan son las fuerzas de presión sobre las secciones AA y BB. La
ecuación de momento resulta
�
pa Sa � pb Sb � m � b Vb � � a Va � (5.65)
h fe �
� Va � Vb �
2
(5.67)
2
A B
pb
pa
V̄a Sa V̄b Sb
Dirección de flujo
A B
FIGURA 5.13
Flujo en una expansión súbita de la sección transversal.
– –
De acuerdo con la ecuación (4.13), V b = V a (Sa /Sb) como ρ es constante, la ecuación
(5.67) queda
2
V 2a � Sa �
h fe � �1 � � (5.68)
2 � Sb �
Si el tipo de flujo entre las dos secciones es diferente, han de tenerse en cuenta los
factores de corrección α y β. Por ejemplo, si el flujo es laminar en la tubería expandida
y turbulento en la estrecha, se considera αb igual a 2 y βb como –43 en las ecuaciones
(5.65) y (5.66).
C
Dirección del flujo
Plano C-C de la vena contracta
FIGURA 5.14
Flujo en una contracción súbita de la sección transversal.
V b2
h fc � Kc (5.70)
2
donde Kc es el factor de proporcionalidad, que recibe el nombre de coeficiente de pér-
–
dida por contracción y V b es la velocidad media de la corriente a la salida de la sección
estrecha. En forma experimental, se encuentra que para el flujo laminar, Kc < 0.1, y la
pérdida de contracción hfc es insignificante. Para el flujo turbulento, el valor de Kc está
dado por la ecuación empírica
� S �
Kc � 0.4 �1 � b � (5.71)
� Sa �
donde Sa y Sb son las áreas de sección transversal de los ductos de las corriente al ingreso
y a la salida, respectivamente. Cuando se descarga un líquido a través de una tubería sol-
dada a la pared de un gran tanque, Sb/Sa es casi cero y Kc 0.4. El valor de Kc no puede
ser predicho con exactitud a partir de la teoría, sino que corresponde al área de una vena
contracta de aproximadamente 0.6 Sb. La velocidad en la vena contracta es entonces de
–
alrededor de Vb/0.6 y la energía cinética es (1/0.6)2 o 2.8 veces la energía cinética de la
corriente final. Sin embargo, de acuerdo con la ecuación (5.68), la fracción de pérdida de
energía cuando la corriente se expande más allá de la vena contracta es sólo (1 – 0.6)2, o
una pérdida de que corresponde a Kc = 0.44.
V a2
h ff � K f (5.72)
2
donde Kf = factor de pérdida para el accesorio
–
V a = velocidad media en la tubería que conduce al accesorio
TABLA 5.1
Coeficientes de pérdida por fricción en accesorios de tubería con cuerda estándar9b
Accesorio Kf
Codo estándar
45° 0.35
90° 0.75
T
En líneas rectas 0.4
Usados como codo 1.0
Retorno, 180° 1.5
Válvula de compuerta
Semiabierta 4.5
Completamente abierta 0.17
Válvula de ángulo, completamente abierta 2.0
Válvula de globo, completamente abierta 6.0
–
a la ecuación (5.7), es 4f(L/D)(V 2/2); la pérdida de contracción a la entrada del tubo, de
–
acuerdo con la ecuación (5.70), es Kc(V 2/2); la pérdida por expansión a la salida del tubo,
–2
según la ecuación (5.64), es Ke(V /2); y la pérdida por fricción en la válvula de
–
globo, según la ecuación (5.72), es Kf (V 2/2). Si la fricción de superficie en la entrada y
salida de los cabezales es insignificante, la fricción total es
L V2
h f � � 4 f � Kc � Ke � K f � (5.73)
� D � 2
a b
Presión = p1
Flujo
Presión = pb
Presión = pa
Velocidad = V̄ b
Velocidad = V̄ a
Velocidad en el tubo = V̄ Presión = p2
a L b
Cabezal de entrada Cabezal de salida
FIGURA 5.15
Flujo de un fluido incompresible a través de una instalación típica.
pa � pb L V2
� g� Za � Zb � � � 4 f � Kc � Ke � K f � (5.74)
� � D � 2
EJEMPLO 5.2 El petróleo crudo, que tiene una gravedad específica de 0.93 y una viscosidad
de 4 cP, se descarga por gravedad desde el fondo de un tanque. La profundidad del líquido
sobre la conexión de descarga en el tanque es de 6 m. La línea de descarga es una tubería 3 in.
norma 40 para la tubería. Su longitud es de 45 m, y contiene un codo de 90 ° y dos válvulas de
compuerta. El crudo descarga a la atmósfera en un punto situado a 9 m por debajo del punto
de conexión al tanque. ¿Qué velocidad de flujo, en métros cúbicos por hora, cabe esperar a
través de la línea de descarga?
� � 0.004 kg/m � s L � 45 m
3.068
D= � 0.256 ft � Ap. 3� � 0.078 m
12
� = 0.93 � 998 = 928 kg/m 3
–
A partir de la ecuación (4.71), suponiendo que αb = 1, y ya que pa = pb y Va = 0,
V b2
� h f � g� Za � Zb � � 9.80665�6 � 9� � 147.1 m 2 /s2 (5.75)
2
Se utiliza la ecuación (5.73). En este caso no hay pérdida de expansión final, ya que la corriente
no se expande en la descarga, y Kc = 0. Según la ecuación (5.71), como Sa es muy grande,
Kc = 0.4. Por lo tanto,
V2
�
L
h f � � 4 f � Kc � Kf � b
� D � 2
2
� 4 � 45 f �V V2
� � 0.4 � 1.09 b � �2 308 f � 1.49� b
� 0.078 � 2 2
V 2b V2
� h f � b �1 � 2 308 f � 1.49� � 147.1
2 2
147 .1 � 2 294.2
V 2b � �
2.49 � 2 308 f 2.49 � 2 308 f
– –
Vb,est., m/s Re 10–4 f (de la fig. 5.10) Vb,cal., m/s
Las velocidades de 2 a 4 m/s de los ejemplos 5.1 y 5.2 son características del flujo
de líquidos en tuberías de tamaño medio. La velocidad óptima o el tamaño óptimo de
tubería para una determinada velocidad de flujo depende de los costos de bombeo y del
capital invertido. En el capítulo 8 se da una ecuación empírica para la velocidad óptima.
Para flujo de gases, la velocidad óptima es mucho mayor que para los líquidos, debido
a la menor densidad de aquéllos.
b
a
a
b
s
A
C
FIGURA 5.16
Separación de la capa límite en un canal divergente.
accesorios, las expansiones y las contracciones son numerosas, las pérdidas por fricción
en la tubería sola resultan despreciables.
Los vórtices formados entre la pared y la corriente separada del fluido, una vez
sobrepasado el punto de separación, dan lugar a grandes pérdidas por fricción de forma.
La separación tiene lugar tanto para el flujo laminar como para el turbulento. El punto de
separación está más alejado en el conducto para el flujo turbulento que para el laminar.
Es posible evitar la separación si el ángulo formado entre la pared y el eje del conducto
es suficientemente pequeño. El ángulo máximo tolerable para una expansión cónica sin
que se produzca separación es de 7°.
SÍMBOLOS
b Ancho de un ducto cuadrado, m o ft; también espacio entre placas paralelas, m o ft
cυ Calor específico a volumen constante, J/kg · °C o Btu/lb · °F
D Diámetro, m o ft; Deq, diámetro equivalente de un canal no circular; Di, diámetro
interior del espacio anular; Do, diámetro exterior del espacio anular
F Fuerza, N o lbf; Fg, fuerza de gravedad; Fs, fuerza de corte; Fw, fuerza neta de
la pared del canal sobre el fluido
f Factor de fricción de Fanning, adimensional
gc Factor de proporcionalidad de la ley de Newton, 32.174 ft · lb/lbf · s2
hf Pérdida de fricción, J/kg o ft · lbf /lb; hfc, en una contracción súbita; hfe, de una
expansión súbita; hff, en flujo a través de accesorios o válvulas; hfs, fricción de
superficie
K Constante en las ecuaciones (5.25) y (5.27)
Kc Coeficiente de pérdida por contracción, adimensional
Ke Coeficiente de pérdida por expansión, adimensional
Kf Factor de pérdida para accesorios o válvulas, adimensional
K Índice de consistencia del flujo, kg/m · s2–n lb/ft · s2–n
k Parámetro de rugosidad, m o ft
L Longitud, m o ft
Lp Perímetro húmedo, m o ft
.
m Velocidad de flujo másico kg/s o lb/s
Letras griegas
α Factor de corrección de la energía cinética, adimensional; αa, en el punto a; αb,
en el punto b
β Factor de corrección del momento, adimensional; βa, en el punto a; βb, en el
punto b
p Pérdida de presión, pa – pb; ps, por la fricción de superficie
T Incremento de la temperatura por la disipación viscosa, °C
κ Relación ri/ro en flujo anular [ecuación (5.28)]
µ Viscosidad absoluta, P o lb/ft · s; µw, a la temperatura de la pared de la tubería
ν Viscosidad cinemática, µ/ρ, m2/s o ft2/s
ρ Densidad, kg/m3 o lb/ft3
τ Esfuerzo cortante, N/m2 o lbf /ft2; τw, en la pared de la tubería; τ0, tensión mínima
en un fluido plástico
φa Factor en la ecuación (5.30) para flujo laminar en el espacio anular, adimen-
sional
ψ Factor de corrección de la temperatura para la fricción de superficie, adimen-
sional
PROBLEMAS
5.1. Demuestre que el flujo de un líquido en régimen laminar entre placas planas paralelas infinitas
está dado por
12 �VL
pa � pb �
b2
5.3. Por un conducto rectangular, con una elevada relación entre las caras, circula un fluido
newtoniano con flujo laminar. Deduzca la relación entre la velocidad local y la máxima, y
–
determine también la relación umáx /V .
5.4. Calcule la potencia requerida por metro de ancho de la corriente, para hacer pasar un aceite
lubricante a través del espacio que existe entre dos placas planas horizontales bajo las si-
guientes condiciones:
5.5. ¿Cuál sería el incremento de la temperatura en el aceite causado por la disipación viscosa
bajo las condiciones del problema 5.4 si el calor específico del aceite es 0.5 cal/g · °C y su
densidad es de 820 kg/m3?
5.6. Un líquido con una gravedad específica de 2.6 y viscosidad de 2.0 cP fluye a través de una
tubería lisa de diámetro desconocido, dando lugar a una caída de presión de 0.183 lbf /in.2
por 1.73 mi. ¿Cuál es el diámetro de la tubería en pulgadas si la velocidad de masa del flujo
es de 7 000 lb/h?
5.7. A través de una tubería de acero de 8 in. circula agua con una velocidad promedio de 6 ft/s.
La corriente a la salida de la tubería se divide en una rama principal de 8 in. y una desviación
alterna de 2 in. La longitud equivalente de la desviación es de 22 ft; la longitud de la tubería
es de 8 in. y en la sección que comprende la desviación es de 16 ft. Despreciando las pérdidas
de entrada y salida, ¿qué fracción del flujo total del agua pasa a través de la desviación?
5.8. Una tubería de acero de 2 ft de diámetro conduce agua a 15 ft/s. Si la tubería tiene una
rugosidad de 0.0003 ft, ¿podría aumentarse la capacidad mediante un recubrimiento liso de
plástico que reduzca el diámetro interior a 1.9 ft? Calcule la variación de la caída de presión
para el mismo flujo y el cambio en la capacidad para una caída de presión dada.
5.9. Se bombea agua a 60 °F desde un recipiente hasta la cima de una montaña a través de una
tubería de 6 in. Norma 120, con una velocidad media de 12 ft/s. La tubería descarga a la
atmósfera a una altura de 3 000 ft sobre el nivel del recipiente. La longitud de la tubería de
4 500 ft. Si la eficacia total de la bomba y el motor que la acciona es de 70%, y si el costo
de la energía eléctrica es de 4 centavos por kilowatt-hora, ¿cuál es el costo de la energía por
hora para el bombeo de esta agua?
5.10. Una unidad de ósmosis inversa para purificación de agua tiene cerca de 900 000 fibras huecas
que permiten la difusión del agua pero impiden el paso de la mayor parte de la sal. Las fibras
5.19. a) ¿Cuál es la velocidad de flujo máximo en gal/min que circula a través de 120 ft de una
manguera de incendios de 2.5 in de diámetro, si la presión manométrica de descarga de la
bomba es de 125 lbf /in2? Asuma que la superficie de la manguera es dos veces más rugosa
que una tubería de acero; b) ¿cuánto se podría aumentar el flujo si el agua contuviera 50 ppm
de óxido de etileno?
5.20. Un dispositivo para el análisis químico de pequeñas porciones de solución acuosa (un
“laboratorio en un chip”) posee canales de 100 µm de ancho y 10 µm de profundidad. Si la
viscosidad de la solución es de 1.3 cP y la densidad, de 1 050 kg/m3, ¿cuáles son el número
de Reynolds y la caída de presión por unidad de longitud cuando la velocidad promedio es de
0.1 m/s? Exprese la caída de presión en Pa/m y en bar/cm.
5.21. Para el flujo en un chip con canales cuadrados de 300 µm por lado, ¿qué velocidad del fluido,
en m/s, corresponde a un número de Reynolds de 10, si el fluido es agua a 20 ºC?, ¿cuál es
si es aire a 70 ºC?
5.22. De un conducto cuadrado de 3 m de lado de un incinerador municipal fluye gas a 800 ºC y a
aproximadamente 1 atm. La viscosidad y la densidad son cercanas a las del aire a 800 ºC. a)
Para una velocidad de un gas de 4.0 m/s, ¿cuál es la caída de presión por unidad de longitud?;
b) Estime el espesor de la subcapa viscosa y la de la capa amortiguadora (buffer); c) ¿Qué
fracción del flujo total se encuentra en cada una de estas capas?
–
5.23. Utilice la distribución de velocidad universal para calcular la relación V / umax para Re =
104, 105 y 106. Use las propiedades del agua a 20 ºC para simplificar el cálculo. Compare sus
relaciones con los valores de la figura 5.8.
REFERENCIAS
1. Bird, R.B., W.E. Stewart y E.N. Lightfoot. Transport Phenomena. Nueva York: Wiley, 1960,
pp. 51-53.
2. Dodge, D.W. y A.B. Metzner. AIChE J. 5: 189 (1959).
3. Garde, R.J. Turbulent Flow. Nueva York: Wiley, 1994.
4. Kays, W.M. y A.L. London. Compact Heat Exchangers. 2a. ed. Nueva York: McGraw-Hill,
1964.
5. Knudsen, J.G. y D.L. Katz. Fluid Dynamics and Heat Transfer. Nueva York: McGraw-Hill,
1958, pp. 97, 101-105, 158-171.
6. Lumley, J.A. Physics of Fluids 20(10): S64 (1977).
7. Middleman, S. An Introduction to Fluid Dynamics. Nueva York: Wiley, 1998, p. 469.
8. Patterson, R.W. y F.H. Abernathy. J. Fluid Mech. 51: 177 (1972).
9. Pearson, J.R.A. Mechanics of Polymer Processing. Londres: Elseiver, 1958, pp. 191-194.
10. Perry, R.H. y D. W. Green, eds. Perry’s Chemical Engineers’Handbook. 7a. ed. Nueva York:
McGraw-Hill, 1997, a) p. 6-12; b) p. 6-18.
11. Prandtl, L. VDIZ. 77: 105 (1933).
12. Schlichting, H. Boundary Layer Theory. 7a. ed. Nueva York: McGraw-Hill, 1979.
13. Toms, B.A. Proc. Intern. Congr. Rheology, Holland. Amsterdam: North-Holland, 1949, p.
II-135.
14. Vinogradov, G.V., G.B. Froishteter, K.K. Trilisky y E.L. Smorodinsky. Rheologica Acta 14: 765
(1975).