DAIMF - UNMSM 5 - 19

FLUJO TURBULENTO

5.6 ECUACION DE DISTRIBUCION DE VELOCIDAD DE FLUJO TURBULENTO EN

TUBERÍAS RUGOSAS
Igual que para tuberías lisas, veamos primero cuáles son las variables adimensionales
que intervienen en la ecuación de distribución de velocidad de flujo turbulento en
tuberías rugosas. En este caso, se puede suponer que la velocidad promedio temporal
( Vz ) depende de:
a).- Propiedades del fluido:  ,  ;
b).- Distancia del punto, medida desde la pared interior del tubo: y ;
c).- Condiciones de la pared del tubo, representadas por el esfuerzo cortante en la
pared (  0 ) y la rugosidad promedio de las paredes (  ).
Así, se puede plantear la siguiente relación funcional:
Vz  F (  ,  , y,  0 ,  ) ; (5.57)
y por el método de análisis dimensional se encuentra la siguiente posible relación
funcional adimensional:
Vz V y y
 F1 ( * , ) (5.58)
V*  

La condición del flujo en tubos rugosos depende de la magnitud relativa de la rugosidad

de la pared (  ) y del espesor de la subcapa viscosa (  ' ), de la siguiente manera:
 Si  se vuelve grande comparado con  ' , la rugosidad de la pared interfiere la
formación de la subcapa viscosa y de la zona de transición. Tubos con este tipo de
pared son denominados completamente rugosos (ver figura 5.8b) por que el perfil
y
de la velocidad es una función solamente de ( ) como muestra la figura 5.9. (PAO,

1967, p, 332).
 Por otro lado, si    ' , la pared del tubo es cubierta por la subcapa viscosa y los
tubos con pared de este tipo se denominan hidráulicamente lisos (ver figura 5.8a).
V y
En estas condiciones el perfil de velocidad es una función únicamente de ( * ) .

 Por último, como no hay un cambio brusco entre una condición y otra, entre estos
dos extremos habrá una pared del tubo que es clasificada como transición entre
lisa y rugosa.

(a) Pared hidráulicamente lisa (b) Pared rugosa

Figura 5.8: Tipos de pared de conductos
(Fuente: POTTER & W IGGERT, 1998, p, 286)
DAIMF - UNMSM 5 - 20
FLUJO TURBULENTO

5.6.1 TUBERÍA DE PARED COMPLETAMENTE RUGOSA

La ecuación algebraica que mejor se ajusta a los datos experimentales para flujo
turbulento en tubos completamente rugosos se muestra en la figura 5.9, y puede ser
derivada desde la ecuación (5.38) [perfil de velocidad para núcleo turbulento de tubos
lisos] de la siguiente manera:
V z V z max r /
En (5.38):   2.5 ln( 0 ) ,
V* V* y /
de donde,
V z V z max r y
  2.5 ln( 0 )  2.5 ln( ) .
V* V*  
   
B
Representando por B los dos términos del lado derecho sobre el corchete, esta
ecuación queda generalizada como
Vz y
 B  2.5 ln( ) (5.59)
V* 

Figura 5.9:
Verificación experimental de la
ecuación de distribución de
velocidad para flujo turbulento
en tuberías rugosas, en escala
semi-logarítmica.
(Fuente: PAO, 1967, p, 331)
DAIMF - UNMSM 5 - 21
FLUJO TURBULENTO

V z max r
donde B   2.5 ln( 0 ) toma un valor constante para cualquier flujo turbulento en
V* 
tubo completamente rugoso. Experimentalmente se ha encontrado que este valor
constante es igual a 8.50. Por tanto, el perfil de velocidad del flujo turbulento en un tubo
completamente rugoso puede ser escrito como (con ln  2.3 log ):
Vz y y
 8.50  2.5 ln( )  8.50  5.75 log( ) (5.60a)
V*  
A esta ecuación, SOTELO (2009, p, 542) denomina “ley universal de distribución de
velocidades, para flujo turbulento en tubos rugosos”.
Por otro lado, Nikuradse, por sus estudios experimentales realizados con rugosidad
artificial y homogénea, llegó a la conclusión que la validez de la ecuación (5.36) puede
extenderse hasta (ROCHA, 1979, p, 85)

y0  (5.60b)
30
Al reemplazar 5.60(b) en (5.36) se obtiene la siguiente ecuación de distribución de
velocidades en un contorno rugoso (tubería o canal)
V 30 y
V z  * ln( ), (5.60c)
k1 
que para k1  0.4 resulta la ecuación 5.60(a).

5.6.2 TUBERÍA DE PARED QUE NO ES COMPLETAMENTE RUGOSA

Si la pared del tubo no es completamente rugosa, la constante B es dependiente de la
rugosidad de la pared (HANSEN, 1967, p, 425). Se ha hallado experimentalmente que el
valor numérico de B es función de  , tal como muestra la figura 5.10, en la cual se han
graficado resultados experimentales obtenidos en tubos con rugosidad artificial.

Figura 5.10: Relación entre el valor numérico de B y de  [este último expresado

en forma de “número de Reynolds de rugosidad”].
(Fuente: PAO, 1967, p, 332)

Los datos experimentales indican los siguientes límites para definir las condiciones de
rugosidad de la pared de un conducto (ver también HANSEN, 1967, p, 462):
DAIMF - UNMSM 5 - 22
FLUJO TURBULENTO

V* 
 Hidráulicamente lisa: cuando 5

V* 
 Transición entre lisa y rugosa: cuando 5   70 1

V* 
 Completamente rugosa: cuando  70

Estos límites son numéricamente similares a aquellas encontradas para las diferentes
zonas en el interior de una tubería lisa, excepto que y ha sido reemplazado por  . A
partir de estas evidencias es importante notar que:
V
 Un tubo rugoso puede ser considerado hidráulicamente liso en tanto que * 5;

 La velocidad adimensional para flujo turbulento en tubería completamente rugosa
(ecuación 5.60a o 5.60b), es independiente de la viscosidad del fluido.

5.6.2.1 Tubería de pared hidráulicamente lisa

En la región hidráulicamente lisa, la curva de la figura 5.10 se ajusta a la siguiente recta
(HANSEN, 1967, p, 426):
V  V 
B  5.5  5.75 log( * )  5.5  2.5 ln( * ) (5.61a)
 
Al sustituir esta expresión de B en la ecuación generalizada (5.59), resulta
Vz V  y V y
 5.5  2.5 ln( * )  2.5 ln( )  5.5  2.5 ln( * ) , es decir se obtiene que
V*   
Vz V y V y
 5.5  2.5 ln( * )  5.5  5.75 log( * ) (5.61b)
V*  
La ecuación 5.61(b) es el perfil de velocidad para flujo turbulento en tubos hidráu-
licamente lisos y, tal como se esperaba, es la misma que la ecuación (5.55), esto es, el
perfil de velocidad en el núcleo turbulento de tubos lisos. Adicionalmente, tal como se
presentó al final de la Sección 5.5.2, para conductos hidráulicamente lisos (canal o
tubería) la distribución de velocidades también está dada por la ecuación (5.51).
V z max r
Por otra parte, si la expresión de B   2.5 ln( 0 ) se reordena de la siguiente forma:
V* 
V z max 1 V* r0
B  2.5 ln( ),
V* V  
( * )

V* V z max V r
para  1 esta expresión reordenada se transforma en B   2.5 ln( * 0 ) y
 V* 
toma el valor 5.5 en la gráfica 5.10; es decir, B se transforma en la constante (Const ) de
la ecuación (5.55), tal como se explica en la forma alternativa de hallar esta ecuación en
la Sección 5.5.3. Por lo tanto, se demuestra que 5.5 es el valor límite inferior de la curva
para pared hidráulicamente lisa y, a su vez, es el valor de la constante de la ecuación
(5.55) para tubos lisos.

1
Algunos autores consideran este límite como 80, ver por ejemplo, PAO(1967, p, 332).
DAIMF - UNMSM 5 - 23
FLUJO TURBULENTO

5.6.2.2 Tubería de pared en transición entre lisa y rugosa

[Cuando la pared de la tubería no es completamente rugosa ni es hidráulicamente lisa, el
valor de B , tal como puede observarse en la figura 5.10, varía entre 8.5 y la relación
lineal en escala semilogorítmica dada por la ecuación 5.61(a). En estas condiciones, el
V y
perfil de velocidad del flujo depende tanto de ( * ) como de ( ), y su expresión estará
 
dada por la ecuación generalizada (5.59) con un valor de B obtenido de la figura 5.10.
V V z max r
Por ejemplo, para *  10 , en la figura 5.10, B   2.5 ln( 0 )  9.6 . Entonces, el
 V* 
Vz y y
perfil de velocidad sería:  9.6  2.5 ln( )  9.6  5.75 log( ) .
V*  
DAIMF - UNMSM 5 - 24
FLUJO TURBULENTO

5.7 ECUACION DE VELOCIDAD MEDIA EN TUBERIAS CIRCULARES

5.7.1 TUBERÍAS DE PARED LISA

Por definición, la velocidad media U del flujo en un conducto de sección transversal A ,
que conduce un caudal volumétrico Q , es: U  Q / A . Si A es la sección del tubo de
radio r0 , y el caudal se calcula con una fórmula integral tal como Q   V z dA , en el
A
presente caso, la integración debería efectuarse para las ecuaciones de perfil de
velocidad promedio temporal halladas para las 3 zonas dentro del tubo liso. Sin
embargo, es permisible hacer algunas aproximaciones sin afectar la exactitud del
resultado. Así por ejemplo, se puede extender la ecuación del núcleo turbulento hasta el
límite entre la subcapa viscosa y la zona de transición, y plantear la siguiente fórmula
integral para hallar la velocidad media
r  r0 

Q
 V z (2 r ) dr
U  r 0 (5.62a)
A  r02
Según los límites de la integración, al caudal volumétrico calculado con la fórmula
(5.62a), faltaría agregar el caudal a través de la subcapa viscosa de espesor nominal  .
La integral de la fórmula (5.62a) está expresada para coordenadas r  z . Pero, como
las ecuaciones de perfil de velocidad se han hallado para coordenadas y  z , será
necesario transformar (5.62a) para estas coordenadas, de la siguiente manera:
Según (5.21): y  r0  r . Entonces, diferenciando: dy   dr . Asimismo, en (5.21),
para r  0 ,  y  r0 ; y para r  r0   ,  y   . Con estas relaciones, la
fórmula 5.62(a) en coordenadas y  z , queda expresada como:
y   r0

 V z 2 (r0  y )( dy)   2 V z (r0  y ) dy 2  V z (r0  y) dy

y  r0 r0 
U  
 r02  r02  r02
es decir,
r0
2
U
r02  V z (r0  y) dy (5.62b)

Entonces, para el perfil de velocidad dada por la ecuación (5.55) (esto es, con la
ecuación para el núcleo turbulento), según 5.62(b) se tiene la siguiente expresión para la
velocidad media:
V*  y 
r0
2  2
U
r02  V* 5.5  2.5 ln(  
)(r0  y) dy  2 I ,
r0
(5.63)


donde I es la siguiente integral:

V y 
r
0

I   V* 5.5  2.5 ln( * ) (r0  y ) dy , (5.64)
   
que puede ser efectuada de la siguiente manera:
r0
 V y V y 
I  V*  5.5 r0  2.5 r0 ln( * )  5.5 y  2.5 y ln( * ) dy
  
DAIMF - UNMSM 5 - 25
FLUJO TURBULENTO

 
r0
 V*   V*
I  V*  
 5 . 5 r0  2 . 5 r0 ln(
 

)  ln y 

 5 . 5 y  2 . 5 y ln(
 

)  ln y  dy

r0
 V V 
I  V*  5.5 r0  2.5 r0 ln( * )  2.5 r0 ln y  5.5 y  2.5 y ln( * )  2.5 y ln y dy
    
r
 V* 5.5 2 2.5 2 V*  0
I  V* 5.5 r0 y  2.5 r0 ln( ) y  y  y ln( ) 
  2 2  


r0 r0


V* 2.5r0  ln y dy  2.5 y ln y dy (5.65)

   

En esta ecuación, las integrales de los dos últimos términos del lado derecho se puede
efectuar empleando, por ejemplo, el método de integración por partes. Así,
1
sean: u  ln y  du  dy ,
y
dv  dy  v  y ;
entonces,
1
 ln y dy   udv  uv   vdu  y ln y   y ( y d y)
 ln y dy  y ln y   dy
 ln y dy  y ln y  y , (5.66)

Similarmente,
1
sean: u  ln y  du  dy
y
y2
dv  y dy  v ;
2
entonces,
y2 y2 1
 y ln y dy   udv  uv   vdu  2
ln y  
2 y
( d y)

y2 y y2 y2
 ln y dy  2
ln y   dy 
2 2
ln y 
4
, es decir,

y2 y2
 y ln y dy  2
ln y 
4
, (5.67)

Sustituyendo (5.66) y (5.67) en (5.65), resulta

r0
 V 5.5 2 2.5 2 V* 
I  V* 5.5 r0 y  2.5 r0 ln( * ) y  y  y ln( ) 
  2 2  
r0
 y2 y2  (5.68)
V* 2.5r0 ( y ln y  y )  2.5( ln y  ) 
 2 4 

Aquí se hace una simplificación matemática haciendo   0 , con el cual se está

despreciando el flujo dentro de la subcapa viscosa. Entonces,
DAIMF - UNMSM 5 - 26
FLUJO TURBULENTO

r0
 V 5.5 2 2.5 2 V* 
I  V* 5.5 r0 y  2.5 r0 ln( * ) y  y  y ln( ) 
  2 2  0
r0
 y2 y2 
V* 2.5r0 ( y ln y  y )  2.5( ln y  ) 
 2 4 0
que evaluando para los límites, resulta
 V 5.5 2 2.5 2 V* 
I  V* 5.5 r02  2.5 r02 ln( * )  r0  r0 ln( ) 
  2 2  

 r2 r2 
V* 2.5r02 (ln r0  1)  2.5( 0 ln y  0 )  V* 0

 2 4 
 V 5.5 2.5 V* 1 1 
I  V* r02 5.5  2.5 ln( * )   ln( )  2.5(ln r0  1)  2.5( ln r0  )
  2 2  2 4 
 5.5 2.5 V* 2.5 2.5 
I  V* r02   ln( )  2.5 ln r0  2.5  ln r0  
 2 2  2 4 
 5.5 2.5 V* 2.5 2.5 
I  V* r02   ln( )  ln r0  2.5  
 2 2  2 4 

 2.5 V* 2.5   2.5 V 

I  V* r02  ln( )  ln r0  0.875  V* r02  (ln( * )  ln r0 )  0.875
 2  2   2  

 2 .5 V r 
I  V* r02  (ln( * 0 )  0.875  (5.69)
 2  
Finalmente, reemplazando (5.69) en (5.63), resulta,
2  2.5 V r 
U V* r02  (ln( * 0 )  0.875 
r02  2  

 V r   V r 
U  V* 1.75  2.5 ln( * 0 )   V* 1.75  5.75 log( * 0 )  (5.70)
     

5.7.2 TUBERÍAS DE PARED COMPLETAMENTE RUGOSA

Para flujo turbulento en tubería rugosa, la velocidad del flujo está dada por la ecuación
(5.60a), es decir
 y 
V z  V* 8.50  2.5 ln( ) (5.71)
  
Entonces, de acuerdo a (5.62b) se tiene la siguiente expresión para velocidad media:
r0
2  y  2
U
r02  V* 8.50  2.5 ln(  ) z (r0  y) dy  r02 I , (5.72)

donde la integral I en este caso es

r
 y 
0

I   V* 8.5  2.5 ln( ) (r0  y ) dy (5.73)

   
Efectuando la integración,
DAIMF - UNMSM 5 - 27
FLUJO TURBULENTO

 y 
r0
y
I  V*   8.5 r0  2.5 r0 ln( )  8.5 y  2.5 y ln( )  dy
    
r0

I  V*  8.5 r 0
 2.5 r0 ln y  ln    8.5 y  2.5 yln y  ln   dy

r0

I  V*  8.5 r

0
 2.5 r0 ln y  2.5 r0 ln   8.5 y  2.5 y ln y  2.5 y ln  dy

 
r0
 
r r
8.5 2 2.5 2 0 0

I  V* 8.5 r0 y  2.5 r0 y ln   y  y ln    V* 2.5r0  ln y dy  2.5 y ln y dy 

 2 2     
En esta ecuación, sustituyendo (5.66) y (5.67), resulta
r
 
0
8.5 2 2.5 2
I  V* 8.5 r0 y  2.5 r0 y ln   y  y ln   
 2 2 
r0
 y2 y2  (5.74)
V* 2.5r0 ( y ln y  y )  2.5( ln y  ) 
 2 4 

Aquí nuevamente se hace la simplificación matemática haciendo que   0 , que en este

caso tiene sustento físico dado que, tratándose de pared rugosa, no se forma la subcapa
viscosa. Entonces,
r0
 y2 
r0
 8.5 2 2.5 2  y2
I  V* 8.5 r0 y  2.5 r0 y ln   y  y ln    V* 2.5r0 ( y ln y  y )  2.5( ln y  )
 2 2 0  2 4 0

que evaluando para los límites, resulta

 8.5 2 2.5 2 
I  V* 8.5 r02  2.5 r02 ln   r0  r0 ln   
 2 2 
 1 
r0 (ln r0  )  0
2.5 2
V* 2.5r02 (ln r0  1) 
 2 2
8.5 2.5 2 2.5 2 7.5 2 
I  V*  r02  r0 ln   r0 ln r0  r0 
2 2 2 4 
r02  7.5  r02  7.5 
I  V* 8.5  2.5 ln   2.5 ln r0    V* 8.5   2.5(ln r0  ln  )
2  2  2  2 
r02  r 
I  V*  4.75  2.5 ln( 0 )  (5.75)
2   
Finalmente, sustituyendo (5.75) en (5.72) resulta
2 r02  r 
U V*  4.75  2.5 ln( 0 ) 
r02 2   
 r   r 
U  V*  4.75  2.5 ln( 0 )   V*  4.75  5.75 log( 0 )  (5.76)
     
DAIMF - UNMSM 5 - 28
FLUJO TURBULENTO

5.8 ECUACION DE DISTRIBUCION DE VELOCIDAD DE FLUJO TURBULENTO EN

TUBERÍAS DE PARED LISA O RUGOSA
Tratándose de flujo turbulento en tuberías lisas, si se resta la ecuación (5.70) de la
ecuación (5.55), resulta:
Vz V y
(5.55) (núcleo turbulento, tubería lisa)  5.50  5.75 log( * )
V* 

U V r
(5.70) (núcleo turbulento, tubería lisa)  1.75  5.75 log( * 0 )
V* 
Vz  U V y V r
 3.75  5.75 [log( * )  log( * 0 )] ,
V*  
reordenando resulta:
Vz U y
  3.75  5.75 log( ) (5.77)
V* V* r0
Del mismo modo, para el caso de flujo turbulento en tuberías rugosas, si se resta la
ecuación (5.76) de la ecuación (5.60), resulta:
Vz y
(5.60) (completamente rugoso)  8.50  5.75 log( )
V* 

U r
(5.76) (completamente rugoso)  4.75  5.75 log( 0 )
V* 
Vz  U y r
 3.75  5.75 [log( )  log( 0 )]
V*  
reordenando, se obtiene:
Vz U y
  3.75  5.75 log( ) (5.78)
V* V* r0
Comparando (5.77) y (5.78) se observa que resulta una ecuación idéntica para ambos
casos (lisos y rugosos). Realmente, no es sorprendente encontrar que la ecuación de
distribución de velocidad para el núcleo turbulento sea idéntica para tubos lisos y
rugosos, puesto que el mecanismo del movimiento turbulento en la región lejos de la
pared del tubo es independiente de las condiciones en la pared.
Es importante enfatizar que, cuando la ecuación (5.78) es usada para graficar el perfil de
velocidad en tubos rugosos, no es necesario conocer las rugosidades de la pared del
tubo. Solamente se necesita determinar la velocidad media del flujo ( U ) en el tubo. En
general, es más fácil determinar la velocidad media que las rugosidades de la pared del
tubo. Así, la ecuación (5.78) tiene una ventaja definitiva sobre la ecuación especial
5.60(a) que es solamente para tubos completamente rugosos.
DAIMF - UNMSM 5 - 29
FLUJO TURBULENTO

5.9 PERDIDAS POR FRICCION EN FLUJO TURBULENTO EN TUBERÍAS RECTAS DE

SECCION CIRCULAR
Otro problema de flujo en tuberías de importancia en la aplicación técnica, es la
determinación de pérdidas por fricción en flujo turbulento. A diferencia del flujo laminar
en tuberías, para el cual se puede obtener una solución matemática (del sistema cerrado
de ecuaciones), la pérdida por fricción en flujo turbulento en tuberías debe ser
determinada con la ayuda de datos experimentales.
En general, en problemas de flujo en tuberías de interés práctico, se supone que el flujo
es unidimensional y unidireccional, y bajo esta hipótesis dichos problemas se pueden
resolver por la aplicación de:
a).- La ecuación de pérdida de carga por fricción, y
b).- La ecuación unidimensional de energía.
En esta Sección hallaremos la ecuación de pérdida de carga por fricción, y en la Sección
5.14 se desarrollará la ecuación unidimensional de energía.

5.9.1 ECUACIÓN DE PÉRDIDA DE CARGA POR FRICCIÓN: EUACIÓN DE DARCY

Si se tiene una tubería recta de diámetro constante a través del cual ocurre flujo
turbulento “permanente”, unidireccional y unidimensional, el componente de la ecuación
de Reynolds en la dirección del eje z , se simplifica a la ecuación (5.9) (ver Sección 5.3):
p 1  1
0  (r rzl )  ( r rzt )  g z (5.9)
z r r r r
Esta ecuación se puede expresar como:
p 1 
0  ( r rz )   g sen (a)
z r r

donde,  rz   rzl   rzt , y g z  g sen (ver figura 5.2).

Asimismo, al reemplazar la expresión de  rz (ecuación 5.18) en (a) se obtiene:
p 1  r
0  (r 0 )   g sen ,
z r r r0
de donde, realizando las operaciones algebraicas, efectuando la derivación y admitiendo
p  p ( p 2  p1 ) ( p 2  p1 )
la aproximación    , se obtiene:
z  z ( z 2  z1 ) L
( p1  p 2 ) 2 0
0    g sen (b)
L r0
Multiplicando (b) por ( A L ), y teniendo en cuenta que Pm  2 r0 y W   AL , resulta:
( p1  p 2 ) A  Wsen   0 ( Pm L)  0 , (5.79)
donde,
p1 , p 2 : presión media sobre la sección transversal del tubo en los puntos 1 y 2,
W : peso del fluido contenido en el tramo del tubo (volumen de control),
A, Pm : área del flujo y el perímetro mojado de la sección, respectivamente,
0 : esfuerzo cortante en la pared lateral del volumen de control,
L : longitud del tramo considerado (longitud del volumen de control), y
 : ángulo de inclinación del tubo.
DAIMF - UNMSM 5 - 30
FLUJO TURBULENTO

La ecuación (5.79) establece que las fuerzas que actúan sobre el tramo del tubo en
análisis se mantienen en equilibrio, es decir,  Fz  0 (ver figura 5.11a).

Figura 5.11(a):
Esquema de fuerzas sobre un tramo del tubo
(volumen de control interior al tubo).

p* r0
Pero, por la ecuación (5.15), se tiene:  0  ( ) ; que se puede expresar como:
L 2
p* /  r0 r
0 ( )   Sw 0   SwR (5.80)
L 2 2
donde:
r0 D A
R   : radio hidráulico de la sección circular, siendo D el diámetro,
2 4 Pm
p* / 
Sw   sen w : pendiente de la línea piezométrica2(ver figura 5.11b), y
L
 : peso específico del fluido contenido en el tramo del tubo (volumen de control).
Para flujo uniforme, como es el presente caso, S w  sen E  S , siendo S la “pendiente
de la línea de energía”. Entonces, en (5.80):
0  SR (5.81)

Figura 5.11(b):
Esquema de pérdida de
carga en flujo uniforme
en una tubería.

2
Con el mismo criterio para hallar la ecuación (5.10), sen    ( p * /  )  S : “pendiente de la línea piezométrica (L.P.)”.
z
w W

Pero,   ( p * /  )    ( p * /  )   ( p2 /   p1 /  )  ( p1 /   p2 /  )  p * /  . Por lo tanto, se tiene:

* * * *
p * /  .
SW 
z z z2  z1 z2  z1 L L
Por otro lado, para flujo uniforme como en el presente caso, S w representa también a la “pendiente (disminución) de la línea
de energía (L.E.)” S , puesto que  E   W (ver figura 5.11b).
DAIMF - UNMSM 5 - 31
FLUJO TURBULENTO

Por otro lado, tratándose de flujo turbulento, la velocidad media se puede calcular con la
siguiente fórmula empírica de Chézy (ver Sección 6.7):
U C R S (5.82)
donde C es denominado “coeficiente de Chézy”. De esta ecuación se tiene:
U2
RS 2 (c)
C
Sustituyendo (c) en la ecuación (5.81), resulta:
U2
0  2 (5.83)
C
Entonces, reemplazando (5.83) en (5.79), y dividiendo entre  A , se obtiene:
( p1  p 2 ) U 2 Pm
 L sen  ( )L  0 (d)
 C2 A

Pero como L sen  h1  h2 (ver figura 5.2), la ecuación anterior se transforma en:

p1 p2 U 2 Pm
(  h1 )  (  h2 )  ( )L (e)
  C2 A
En (e) el término del lado izquierdo representa la pérdida de carga piezométrica3 debido
a la fricción entre las secciones 1 y 2 en un flujo permanente unidimensional e
incompresible, y se designa por h f 1 2 , es decir,
p1 p2 p1  h1 p 2  h2 p1* p 2* p* (5.84)
(  h1 )  (  h2 )  ( )( ) [  ]  h f 1 2
      
Entonces, sustituyendo (5.84) en (e), y teniendo en cuenta que A D , se obtiene:
R 
Pm 4
U 2 Pm U2 1 U2 4
h f 1 2  2
( )L  2 ( )L  2 ( )L (f)
C A C R C D
Esta ecuación multiplicando y dividiendo por 2 g puede tomar la siguiente forma:
8g L U 2
h f 1 2  ( ) (g)
C 2 D 2g

Por último, si se define

8g
f ( ) (5.85)
C2
al que se denomina coeficiente, o factor, de fricción de Darcy, la ecuación (g) queda
expresada como:
LU2
h f 1 2  f (5.86)
D 2g

Esta ecuación es a veces llamada fórmula de Darcy (o fórmula de Darcy-Weisbach en

literatura americana), y el término h f 1 2 es igual a la diferencia de cotas piezométricas
entre las secciones 1 y 2 cuando el flujo es uniforme. En cualquier otro caso, será igual
a la diferencia de carga total entre dichas secciones.

3
El término carga se emplea para designar a la “energía por unidad de peso”. La carga total, como es sabido, es la suma de
la carga de posición, carga de presión y la carga de velocidad. La suma de la carga de posición y la carga de presión
constituye la “carga piezométrica” (ver figuras 5.2 y 5.11b).
DAIMF - UNMSM 5 - 32
FLUJO TURBULENTO

Henry Philibert Gaspard Darcy (Dijon, Francia,

10/06/1803 - París, 02/01/1858), más conocido como
Darcy, fue un hidráulico francés. Graduado como
ingeniero de puentes y caminos, fue uno de los
pioneros modernos en el abastecimiento de agua
potable. Tuvo un papel importante en el desarrollo de
su ciudad natal.
Darcy, entre 1834 y 1840, diseñó y construyó el
sistema de abastecimiento de agua de la ciudad de
Dijon, con una línea de aducción subterránea de 12
km. En 1847, el agua entubada llegó a todos los
pisos de los edificios de de Dijon, transformando así a
esta ciudad en la segunda ciudad europea con Henry P. G. Darcy
abastecimiento de agua, después de Roma. (1803 – 1858)

5.9.2 FACTOR DE FRICCIÓN DE DARCY Y FACTOR DE FRICCIÓN DE FANNING

PARA TUBERÍAS CIRCULARES
El factor de fricción de Darcy ( f ) definida en (5.85), puede ser conceptuado como una
medida relativa de las pérdidas por fricción o rozamiento en el sistema. Una expresión
empírica para este factor se puede obtener a partir de la ecuación (5.83). Así,
reordenando (5.83), y teniendo en cuenta (5.85), se obtiene:
U2 8g  2  (5.87)
 0  g 2
( 2
) U  f U2
C C 8 8
de donde,
0
f  (5.88)
U 2 1
8
Por otra parte, principalmente en la literatura europea, se denomina factor de fricción de
Fanning ( f f ) a la siguiente expresión:
0
ff  (5.89)
U 2 1
2
Al dividir (5.88) entre (5.89), resulta la siguiente relación:
f 4 ff (5.90)

5.9.3 SIGNIFICADO FÍSICO DEL FACTOR DE FRICCIÓN DE DARCY

Para flujo turbulento en tuberías circulares, mediante análisis dimensional se encuentra
la siguiente relación funcional adimensional: f   ( Re ,  / D ) , donde, Re es el número de
Reynolds y  / D se denomina rugosidad relativa, siendo  la rugosidad promedio
introducida en la Sección 5.6, denominada también “rugosidad absoluta” de las paredes
de la tubería.
La rugosidad absoluta  depende de la calidad de las paredes expresada por:
a).- Altura media de las irregularidades de la superficie,
b).- Variación de la altura con respecto a la media,
c).- Forma de las irregularidades de la superficie,
d).- Separación entre irregularidades adyacentes.
DAIMF - UNMSM 5 - 33
FLUJO TURBULENTO

Si bien es cierto que en flujo turbulento el factor de fricción de Darcy f , en el caso más

general, es función de Re y de , también lo es que puede ser función de sólo uno de
D
ellos. A seguir, se presenta las ecuaciones para tales casos.

5.9.4 ECUACIONES PARA CÁLCULO DE FACTOR DE FRICCIÓN DE DARCY

5.9.4.1 Tubería de pared hidráulicamente lisa

Como ya fue expuesto, en este tipo de superficies se desarrolla la subcapa viscosa de
espesor mayor que la rugosidad promedio y por lo tanto ésta no tiene significado en el
cálculo de f, por lo que simplemente f   ( Re ) . Así, tenemos ecuaciones para el cálculo
de f , tales como:
 Ecuación experimental de Blasius
0.316
f  Válida para 3  10 3  Re  10 5 , (5.91)
Re1 / 4

 Ecuación semiempírica de von Kármán-Prandtl (llamada también ecuación univer-

sal de factor de fricción para flujo turbulento en tuberías lisas)
1
 2.0 log( Re f )  0.8 Válida para Re  4  10 3 , (5.92)
f

 Ecuación empírica de Nikuradse

0.221
f  0.0032  (5.93)
Re0.237
Esta ecuación da los mismos resultados que (5.92) para Re entre 105 y 107.

 Ecuación de Konakov
1
f  Válida para Re  2,300 , (5.94)
(1.81 log Re  1.5) 2

En las ecuaciones presentadas R e está referido al diámetro del tubo circular.

Graficando las ecuaciones (5.91) y (5.92) en un papel logarítmico (log-log), se obtiene un
gráfico f vs. Re , llamado diagrama de Stanton, similar a la figura 5.12. En este gráfico
aparece también la ecuación analítica para flujo laminar. Obviamente, esta última
ecuación corresponde a una línea recta.

5.9.4.2 Tubería de pared completamente rugosa

En este tipo de superficies, como ya fue referido en la Sección 5.5, no se desarrolla la
subcapa viscosa. Entonces, el coeficiente de fricción es función exclusiva de rugosidad

relativa, es decir, es apenas de la forma f   ( ) , tal como se muestra a seguir:
D
 Ecuación analítica de von Kármán: Obtenida combinando las ecuaciones (5.76) y
(5.88)
1 3.35 D
 2.03 log( ) (5.95a)
f 
DAIMF - UNMSM 5 - 34
FLUJO TURBULENTO

 Ecuación semiempírica de von Kármán

Ajustando los coeficientes de (5.95a) de acuerdo a los resultados experimentales de
Nikuradse, von Kármán encontró esta otra ecuación:
1 3.71D
 2.0 log( ) (5.95b)
f 
Si quisiéramos hacer un gráfico similar a la figura 5.12 de las ecuaciones para
superficies hidráulicamente rugosas, se obtendría una familia de rectas paralelas al eje
horizontal. Para cada valor de D se obtiene el de f , tal como ilustra la figura 5.13.


Figura 5.12: Relación f vs. Re para tuberías lisas

(Fuente: ROCHA, 1980, p, 115)

Figura 5.13: Relación f vs. Re para tuberías rugosas

(Fuente: ROCHA, 1980, p, 116)
DAIMF - UNMSM 5 - 35
FLUJO TURBULENTO

Ejemplo 5.3: Demostración

Obtener la ecuación (5.95a) a partir de la combinación de: U  V* [4.75  5.75 log( r0 /  )] y
0
f  , ecuaciones (5.76) y (5.88), respectivamente.
1
8
U 2
Solución
Para facilitar la demostración expresemos la ecuación (5.95a) en la siguiente forma
equivalente:
1 3.35 D r
 2.03 log( )  1.68  2.03 log( 0 ) (a)
f  
Por su parte, las ecuaciones (5.76) y (5.88) se pueden expresar como:
U
 4.75  5.75 log( r0 /  ) (b)
V*
0 1 0 1
U2  1  U 8  8 V* (c)
8
 f f  f
Reemplazando (c) en (b), resulta:
1 4.75 5.75
  log( r0 /  ) ;
f 8 8

realizando las operaciones algebraicas indicadas y redondeando a 2 decimales, se

obtiene:
1 r
 1.68  2.03 log( 0 ) (d) = (a) L.q.q.d.
f 

Comentario:
Una comparación de la ecuación (d) con datos experimentales de flujo turbulento en
tubos completamente rugosos, muestra que una ligera modificación de las constantes
mejora el ajuste a los datos probados. La ecuación universal de factor de fricción
modificada para flujo turbulento en tubos completamente rugosos es:
1 r
 1.74  2.0 log( 0 ) (5.96)
f 
r0
que indica un valor constante de f para un valor de ; f es independiente de la

viscosidad del fluido en tubos completamente rugosos (PAO, 1967, p, 339).
DAIMF - UNMSM 5 - 36
FLUJO TURBULENTO

5.9.4.3 Tubería de pared en transición entre lisa y rugosa

En 1930, Nikuradse estudió experimentalmente el comportamiento de las tuberías lisas y
rugosas introduciendo algunos ligeros ajustes en los coeficientes de las expresiones
analíticas. Pero también estudió experimentalmente la fase que corresponde a la
transición entre paredes lisas y rugosas. De este modo, el estudio de Nikuradse
representa en conjunto el comportamiento de las tuberías lisas, rugosas y la transición
entre ambas.
La figura 5.14, conocido como diagrama de Nikuradse, es una síntesis de las figuras
5.12 y 5.13. Debe tenerse presente que este diagrama corresponde a tuberías de
rugosidad artificial y uniforme (obtenida pegando arena tamizada en la superficie interior
de 6 tubos -LOMAX & SAUL, 1979, p, 82). Los valores de  / D para estas tuberías varía
de 0.033 a 0.000986).

Figura 5.14: Gráfico de Nikuradse [escala log-log]

(Fuente: STREETER et al., 2000, p, 291)

Ejemplo 5.4: Cálculo de coeficiente de Chézy

Si ocurre flujo turbulento a través de una tubería de rugosidad absoluta   0.001 m
(concreto sin acabado especial), de 0.20 m de diámetro, calcular el coeficiente de
fricción de Chézy.
Solución
Suponiendo que el flujo turbulento ocurre en un tubo completamente rugoso, el valor de
del coeficiente de fricción de Darcy, f , se puede calcular con la fórmula (5.96). Es
decir,
1 r 0.10
 1.74  2.0 log( 0 )  1.74  2.0 log( )
f  0.001

de donde,
DAIMF - UNMSM 5 - 37
FLUJO TURBULENTO

f  0.0304
Luego, el coeficiente C de Chézy para flujo turbulento se puede calcular a partir de la
ecuación (5.85). En efecto,

8g 8  9.8 m / s 2 m
1/ 2
C   50 .78
f 0.0304 s
1/ 2
m
C  51
s
Alternativamente, C se puede calcular con la siguiente fórmula empírica (ver Sección
6.7.2, fórmula 6.32):
6R
C  18 log( ) (a)
 /2 /7
donde,
R : radio hidráulico;  : rugosidad absoluta;  : espesor de subcapa viscosa.
Se debe tener en cuenta que, para superficies rugosas   0 y la ecuación (a) se
transforma en:
6R
C  18 log( ) (b)
 /2
En cambio, para superficies lisas   0 y la ecuación (a) se transforma en:
6R
C  18 log( ) (c)
 /7
Para el presente ejemplo, ya que se supone que el flujo ocurre en una tubería
completamente rugosa (concreto sin acabado especial), el coeficiente de Chézy C se
6R
puede calcular con la fórmula (b), es decir: C  18 log( )
 /2
D
Reemplazando datos, teniendo en cuenta que R  , resulta:
4
6 D/4 3D 3  0.20
C  18 log( )  18 log( )  18 log( )
 /2  0.001
1/ 2
m
C  50 (similar al primer procedimiento).
s
DAIMF - UNMSM 5 - 38
FLUJO TURBULENTO

5.10 FACTOR DE FRICCION PARA FLUJO EN TUBERIAS COMERCIALES

Si se pretendiera aplicar el diagrama de Nukuradse a tuberías comerciales, cuya
rugosidad no es artificial sino natural, entonces en la zona de transición se encontrarían
fuertes diferencias. Para tuberías comerciales se utiliza otra carta conocida como
Diagrama de Moody, elaborada en 1944 por Lewis F. Moody, que muestra la figura 5.15.
Como ya se dijo, la mayor diferencia entre estos dos diagramas está en la zona de
transición, ya que para el caso de tuberías de rugosidad natural el fenómeno de la
transición es diferente. Esta diferencia se debe a que en una superficie con rugosidad
natural las irregularidades son de diferente tamaño; entonces, basta la presencia de
algunas protuberancias mayores que la media para alterar la subcapa viscosa.
En el diagrama de Moody los valores de f para superficies en transición entre
hidráulicamente lisa y completamente rugosa, se obtienen por medio de la siguiente
ecuación implícita, formulada por Cyril F. Colebrook en 1939 (PAO, 1967, p, 339; ver
también, DAUGHERTY & FRANZINI, 1977, p, 212; LOMAX & SAUL, 1979, p, 83):
1  / D 2.51 
 2.0 log    (5.97)
f  3.71 R e f 

Para calcular f con esta ecuación, llamada fórmula de transición de Colebrook & White,
se emplea proceso iterativo considerando los términos del lado derecho como supuesto
o valor anterior.
Según GERHART et al. (1995, p 536), la ecuación implícita de Colebrook combina
(acopla) la ecuación de factor de fricción de tubos lisos (ecuación 5.92) con la ecuación
para tubos completamente rugosos (ecuación 5.96), y que cubre también el intervalo de
transición entre liso y rugoso del diagrama de Moody (ver también, POTTER & W IGGERT,
1998, p, 295). Observe que la ecuación (5.92) es la ecuación de Colebrook con   0 , y
la ecuación (5.96) es la ecuación de Colebrook con Re   .
Existen ecuaciones explícitas aproximadas de (5.97). Por ejemplo, la siguiente fórmula
de HAALAND (1983)
1  6.9   / D  1.11 
 1.8 log     (5.98)
f  R e  3.71  

Los resultados obtenidos empleando esta relación se encuentran dentro de 2% [de

diferencia] de los obtenidos a partir de la ecuación de Colebrook & White. Así, si se
desean resultados más precisos, se puede aplicar la ecuación (5.98) como una primera
aproximación en los programas de cómputo para hallar f con la ecuación (5.97) (ver
ÇENGEL & CIMBALA, 2006, p, 324).
Otra fórmula explícita para hallar f es la siguiente (ver De NEVERS, 2005, p, 187):
   10 6 
1/ 3

f  0.0055 1   20,000    (5.98a)
  D Re  
que da muy buena aproximación de las curvas en las regiones turbulenta y transición del
diagrama de Moody.

Ejemplo 5.5: Demostración

Hallar la ecuación (5.97) a partir de la siguiente expresión (ver PAO, 1967, p, 339):
1 r0  (r /  ) 
 2.0 log ( )  1.74  2.0 log 1  18 .7 0
f   f R e 

DAIMF - UNMSM 5 - 39
FLUJO TURBULENTO

Figura 5.15: Diagrama de Moody

DAIMF - UNMSM 5 - 40
FLUJO TURBULENTO

5.11 FACTOR DE FRICCION PARA FLUJO EN TUBERÍAS NO CIRCULARES

Para el caso de tuberías de sección diferente a la circular, como por ejemplo cuadrado,
rectangular, oval, etc. (ver figura 5.16), la relación funcional f   ( Re ,  / D ) presentada
en la Sección 5.8.3 evidentemente tendría que ser ampliada a fin de incluir el factor
“forma de sección”. Así, se tendría:

f   ( R e , , forma ) (5.99)
D
Sin embargo, según los estudiosos, los errores que se pueden cometer en la
determinación de la rugosidad tienen una influencia mayor que la que resulta de ignorar
el factor de forma.

Figura 5.16: Secciones no circulares típicas de conductos

De este modo, aceptaremos que en tuberías no circulares la pérdida de carga puede

calcularse con la fórmula de Darcy. Para esto, en las fórmulas se debe emplear ( 4 R ) en
lugar de D .
En efecto, como sabemos, para una sección circular el radio hidráulico (R) es: R  D / 4 .
En consecuencia, D  4R , y la ecuación de Darcy se transforma en:
L U2
hf  f ( ) (5.100)
4R 2 g

Para usar esta fórmula, es importante remarcar lo siguiente: al hallar f en los

diagramas f vs. Re se debe considerar:
U 4R  
Re  y 
 D 4R

Algunos, autores entre ellos FOX & McDONALD (1973), ÇENGEL & CIMBALA (2006, p, 324),
al término ( 4 R ) denominan “diámetro hidráulico” Dh . Así, el diámetro hidráulico se
define de modo que se reduce a diámetro común ( D ) para tuberías circulares, o sea,
4A
Dh   4R (5.101)
Pm
De este modo, para hallar f desde diagramas para tuberías circulares (diagrama de
Nikuradse o de Moody, por ejemplo), el número de Reynolds y la rugosidad relativa se
deben calcular considerando el diámetro hidráulico ( Dh ) como longitud característica.

Ejemplo 5.6: Espesor de la subcapa viscosa en el interior de un tubo liso

Como un ejemplo de las dimensiones de las zonas del flujo en el interior de un tubo liso,
supongamos que fluye agua a 20 ºC por una tubería lisa de 2.5 cm de diámetro y que el
número de Reynolds del flujo es 106. Hallar: velocidad de corte, espesor de la subcapa
viscosa (  ' ) y el espesor de la “zona de transición + núcleo turbulento” (tnt).
DAIMF - UNMSM 5 - 41
FLUJO TURBULENTO

Solución
Para los siguientes datos:
Agua a 20 ºC:   1.0  10 6 m2/s
Dimensiones del tubo: D  0.025 m; r0  0.0125 m
Número de Reynolds del flujo: Re  10 6

Velocidad de corte:
UD   10 6 1.0  10 6  10 6
Como Re   10 6 , se tiene: U   , es decir,
 D 0.025
U  40 m/s (a)

Asimismo, del diagrama de Moody, para Re  10 6 y tuberías lisas, se tiene

f  0.0115 (b)
0
Por otro lado, de la ecuación (5.88): f  1 . Entonces, teniendo en cuenta (a) y (b),
8
U2
y la definición de velocidad de corte:
0 0 0.0115  U 2
0.0115  , de donde,  V 2
 . Por lo tanto,
U 2 
1 *
8
8

V*  U 0.0115 / 8  40 0.001438

V*  1.52 m/s (c)

Espesor de la subcapa viscosa (  ' ):

De la ecuación (5.46):
5
' (d)
V*
Por lo tanto, teniendo en cuente (c) y el valor de  , resulta:
5   5  1.0  10 6
 ' 
V* 1.52
 '  3.29  10 6 m (e)
Así, vemos que el espesor de la subcapa viscosa es extremadamente pequeño.

Espesor de la “zona de transición + núcleo turbulento” (tnt) en la mitad del tubo:

En la mitad del interior del tubo se tendrá:
tnt  r0   '  0.0125  3.29  10 6
tnt  0.0125 m  r0 (f)
En consecuencia, vemos que el espesor de la “zona de transición + núcleo turbulento”
en la mitad del tubo comprende casi la totalidad del radio interior del tubo, es decir, la
zona tnt es el predominante en el flujo.
DAIMF - UNMSM 5 - 42
FLUJO TURBULENTO

5.12 PROBLEMAS TIPICOS DE FLUJOS SIMPLES EN TUBERIAS

Los problemas que se presentan en cálculo de flujos simples en tuberías, que son
básicos para solución de problemas más complejos son:
1).- Cálculo de pérdida de carga ( h f ),
2).- Cálculo de caudal ( Q ), y
3).- Cálculo de diámetro del tubo ( D ).
En cada uno de estos casos, para determinar la magnitud desconocida se utiliza la
fórmula de Darcy, la ecuación de continuidad y el diagrama de Moody (o la ecuación de
Colebrook & White).

5.12.1 PROBLEMAS DE CÁLCULO DE PÉRDIDA DE CARGA

Para este tipo de problemas se requieren los siguientes datos:
 Caudal ( Q ),
 Longitud de tubería ( L ),
 Diámetro del tubo ( D ),
 Viscosidad cinemática del fluido ( ), y
 Rugosidad absoluta del tubo (  ).
Con estos datos se determinan R e y  / D , y h f se calcula con la fórmula de Darcy
determinando f por el diagrama de Moody. El siguiente pseudocódigo, y el diagrama
de flujo de la figura 5.17(a), muestran el procedimiento de cálculo:
Leer: Q , L , D ,  y ;
Calcular: U , R e y  / D ;
Con R e y  / D hallar f en el diagrama de Moody o de Nikuradse, según sea el caso;
Con el valor f , hallar la pérdida de carga ( h f ) con la fórmula de Darcy.

Inicio

Leer: Q, L, D,  , 

4Q UD 
U  ; Re  ;
D 2  D

Hallar f en Moody

L U2
hf  f
D 2g

Fin Figura 5.17(a)

Ejemplo 5.7: Cálculo de pérdida de carga

Se tiene una tubería de 1,000 m de largo, 0.20 m de diámetro, rugosidad absoluta
  0.001 m, a través del cual fluye un líquido de viscosidad   10 6 m2/s a razón de
0.12566 m3/s. Calcular la pérdida de carga por fricción.
DAIMF - UNMSM 5 - 43
FLUJO TURBULENTO

Solución
Para los siguientes datos: L, D,  , n y Q, se desea calcular h f . En este caso, según el
diagrama de flujo de la figura 5.17(a), se tiene (ver también planilla electrónica):
0.12566
U 2  4 m/s
  0.20
4

4  0.20
Re   800 ,000
10  6
 0.001
  0.005
D 0.20

Con ( Re y ), en el diagrama de Moody se obtiene: f  0.0305 . Por lo tanto, la pérdida
D
de carga por fricción será:
L U2 1000 42
hf  f ( )  0.0305 ( )
D 2  9.81 0.20 2  9.81
h f  124 .363 m

Por otra parte, dado que ( Re y ) se ubica en la zona completamente rugosa del
D
diagrama de Moody, el valor de f , alternativamente, se puede calcular con la fórmula
(5.96), ecuación semiempírica de von Kármán, es decir,
1 3.71D 3.71  0.20
 2.0 log( )  2.0 log( )
f  0.001

f  0.0303

5.12.2 PROBLEMAS DE CÁLCULO DE CAUDAL

Para este tipo de problemas se requieren los siguientes datos:
 Longitud de tubería ( L ),
 Diámetro del tubo ( D ),
 Viscosidad cinemática del fluido ( ),
 Rugosidad absoluta del tubo (  ), y
 Pérdida de carga en un tramo de la tubería ( h f ).
En este caso, U y f son desconocidos y una de las maneras de resolver es usando
simultáneamente la fórmula de Darcy y el diagrama de Moody para encontrar sus
valores4. El siguiente pseudocódigo, y el diagrama de flujo de la figura 5.17(b), muestran
el procedimiento de cálculo:
Leer: L , D ,  ,  y h f ;
Calcular  / D ,
Suponer un valor de f en el diagrama de Moody (podría ser, por ejemplo, el que
corresponde a superficie hidráulicamente rugosa, esto es, al valor de  / D );
Calcular un valor tentativo de la velocidad con la siguiente ecuación:

4
Otras dos maneras de resolver problemas de este tipo ver, por ejemplo, en ÇENGEL & CIMBALA (2006, p, 345).
DAIMF - UNMSM 5 - 44
FLUJO TURBULENTO

2 gD h f 2 gD S E (a partir de la fórmula de Darcy, (5.86), donde

U   S E  h f / L );
f L f
Con el valor de U , calcular Re ;
Con R e y  / D , hallar f m en el diagrama de Moody o de Nikuradse, según sea el caso;
Comparar el valor de f m con el f supuesto;
Si f m hallado es diferente de f supuesto, repetir el proceso (tomando el valor de f m ,
hasta conseguir igualdad en las dos primeras cifras significativas);
Verificar que los valores finales de f y de U satisfacen la fórmula de Darcy (esto es,
conseguida la “igualdad” de f m y f , efectuar la verificación);
Con el valor correcto de U , calcular el caudal Q .

Inicio

Leer: L, D,  ,  , h f


Con , en Moody suponer f0
D

f  f0

2g Dhf
U
fL

UD
Re 


Con ( Re ,
 ) hallar f en Moody
m
D

No
f  fm fm  f

Si

L U2
hf  f
D 2g

Q=UA

Fin Figura 5.17(b)

Ejemplo 5.8: Cálculo de caudal

Se tiene una tubería nueva de hierro fundido (   0.00025 m) de 10 in de diámetro. La
longitud es de 1,000 m. Conduce agua de viscosidad   10 6 m2/s. La perdida de carga
en el tramo considerado es 10 m. Calcular el caudal a través de la tubería.
Solución
Para los siguientes datos: L , D ,  ,  y h f , se desea calcular el caudal ( Q ).

En este caso, según el diagrama de flujo de la figura 5.17(b), se tiene:

DAIMF - UNMSM 5 - 45
FLUJO TURBULENTO

 0.00025
  0.00098
D 0.254
Con este valor de rugosidad relativa, se supone f  0.0195 (que corresponde a la zona
completamente rugosa del diagrama de Moody).
2 gD h f 2  9.81  0.254  10
U 
f L 0.0195  1000

U  1.599 m/s
UD 1.599  0.254
Re    4.06  10 5
 10 6

Con ( Re y ), en el diagrama de Moody se obtiene: f  0.0205 . Como este valor
D
difiere del valor supuesto, se repite el proceso con f  0.0205 (ver planilla electrónica).
Finalmente, se obtiene la igualdad en tres cifras significativas del valor de f , resultando:
f  0.020
U  1.579 m/s
Luego se verifica el valor de h f :

1000 1.579 2
h f  0.020 ( )
0.254 2  9.81
h f  10 .006 m  10 m

Finalmente,
  0.254 2
Q  UA  1.579 ( )  0.080 m3/s
4
5.12.3 PROBLEMAS DE CÁLCULO DE DIÁMETRO DE TUBO
Para este tipo de problemas se requiere los siguientes datos:
 Longitud de tubería ( L ),
 Viscosidad cinemática del fluido ( ),
 Rugosidad absoluta del tubo (  ),
 Pérdida de carga en un tramo de la tubería ( h f ), y
 Caudal ( Q ).
En este caso, con D desconocido, hay tres incógnitas de la fórmula (5.86), f , U , D ;
dos en la ecuación de continuidad, U , D ; y tres en la expresión del número de
Reynolds, U , D , Re . La rugosidad relativa (  / D ) también es desconocida. [Es decir,
se tiene 4 incógnitas y 3 ecuaciones, por lo que se emplea un método iterativo]. Usando
la ecuación de continuidad para eliminar U en (5.86), se tiene:
L Q2 8LQ 2
hf  f ; de donde, D 
5
f ,o
D 2g(  D ) 2 hf  2g
2

Q2
D 5  0.0827 f (a) (para g  9.8 m/s2)
SE
DAIMF - UNMSM 5 - 46
FLUJO TURBULENTO

Similarmente, usando ecuación de continuidad en la expresión de Re , se tiene:

UD 4Q 1
Re   (b)
 D
El procedimiento de cálculo se explica con el siguiente pseudocódigo y el diagrama de
flujo de la figura 5.17(c):
Leer: L , Q ,  ,
 y hf ;
Suponer un valor de f ;
Calcular un valor tentativo de D con la fórmula (a), (siempre que g  9.8 m/s )
2

Con el valor de D , hallado en el paso anterior, calcular Re con la fórmula (b),

Calcular  / D ;
Con valores de Re y  / D , hallar f m en el diagrama de Moody (o de Nikuradse);
Comparar el valor de f m con el f supuesto;
Si f m hallado es diferente de f supuesto, repetir el proceso (tomando el valor de f m
hallado, hasta conseguir igualdad en las dos primeras cifras significativas);
Verificar que los valores finales de D y de f satisfacen la fórmula de Darcy (esto es,
conseguida la “igualdad” de f m y f , efectuar la verificación);
“Escoger” el tamaño comercial de D (si h f está entre los valores que corresponden a 2
diámetros comerciales sucesivos, tomar el diámetro mayor).

Inicio

Leer: L,  ,  , h f , Q

Suponer f 0

f  f0

Q2
D 5  0.0827 f
SE

4Q 
Re  ,
D D

Con ( Re ,
 ) hallar f en Moody
m
D

No
f  fm fm  f

Si

L U2
hf  f
D 2g

Escoger valor de D

Fin Figura 5.17(c)

DAIMF - UNMSM 5 - 47
FLUJO TURBULENTO

Ejemplo 5.9: Cálculo de diámetro de tubería

Calcular el diámetro que debe tener una tubería nueva de cemento enlucido
(   0.0004 m) para conducir 2 m3/s de agua de viscosidad   1.2  10 6 m2/s. La longitud
de la tubería es de 1,000 m. La pérdida de carga admisible es de 25 m.
Solución
Para los siguientes datos: L , Q ,  ,  y h f , se desea calcular el diámetro del tubo ( D ).

Esta vez, según el diagrama de flujo de la figura 5.17(c), se tiene

f  0.020 (valor inicial supuesto)

Q2 22
D 5  0.0827 f  0.0827 0.020
(h f / L) ( 25 / 1000 )

D  0.767 m
4Q 4 2
Re  
 D   1.2  10  6  0.767

Re  2.77  10 6
 0.0004
  0.00052
D 0.767

Con ( Re y ), en el diagrama de Moody se obtiene: f  0.0168 . Como este valor
D
difiere del valor supuesto, se repite el proceso con f  0.0168 (ver planilla electrónica).
Finalmente, se obtiene la igualdad en tres cifras significativas del valor de f , resultando:
f  0.0171
D  0.743 m
Q4 2 4
U   4.613 m/s
 D   0.743 2
2

Luego se verifica el valor de h f :

1000 4.613 2
h f  0.0171 ( )  24 .962 m  25 m.
0.743 2  9.81
Finalmente,
D  0.743 m  29.134 in. En consecuencia:
D  29 .5 in o 30 in (diámetro comercial).