HIDRÁULICA

______________________________________________________________________________________

PROFESOR: CARLOS VALENZUELA ING.CIVIL

INGENIERÍA CIVIL
Flujo en Tuberías
Pérdida de presión:

P1/γ

P2/γ

Z1
Z2
p1 v12 p2 v22
 z1    z2   hf
 2g  2g
Caudal constante y sección constante V1=V2 La pérdida de energía
debida a la fricción, se
 p1   p2  evidencia en el flujo como
   z1 
 
   z2 
  hf una pérdida de altura
    piezométrica (presión)

p  p hf hf
(  z )  h f (  z)  d p
(  z) 
g x g x dx g l
Relación del esfuerzo cortante con las pérdidas por fricción

pA 0
dx
 Fx  0
El flujo no se acelera ni se frena
dz
θ
0 (p+dp)A

z
W
datum

pA  ( p  dp) A  W cos   0 Pdx  0    g A d ( p  z )

P dx g
0
 dpA  Agdx cos   0 Pdx  0
A hf
dpA  0 Pdx 0  
  Adz  0 P l
g g
dp  0 Pdx P Perímetro
 A(  dz )  h f Pérdidas por fricción
g g l Longitud del tramo de tubería
Relación del esfuerzo cortante con las pérdidas por fricción

A hf
0  
P l
Ecuación que relaciona las pérdidas debidas a la fricción con el
esfuerzo cortante en la pared de la tubería.

Los científicos de la época (1800-1900) buscaron relacionar las pérdidas

con la velocidad, y de esta forma obtener una relación entre la velocidad y
las perdidas que se pudiera utilizar en el diseño: Ecuación de resistencia al
flujo

1. Determinar las grandes diferencias de comportamiento entre los

diferentes tipos de flujo experimento de Reynolds
Número de Reynolds

L=Longitud Característica
d  diámetro en tuberías
Rh  conductos de sección no circular (A/P)

Para flujo en tuberías

Reynolds encontró que:

Re < 2000 Flujo Laminar
2000 < Re < 5000 Flujo en Transición
Re > 5000 Flujo Turbulento
Relación entre el tipo de flujo y la pérdida de presión

A partir de un segundo experimento, Reynolds obtuvo

Tubería rugosa
Log p
Tubería lisa
l 2 4
1
1.75
1
Variación Potencial entre las
3
perdidas de presión por
2 unidad de longitud y V
1
1
1

Flujo laminar F. Transición Flujo turbulento Log v

Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

Re>5000
Variación lineal de la perdida de presión por unidad de longitud y V
A medida que velocidad aumenta, aumentan las perdidas pero con un
comportamiento diferente dependiendo del tipo de flujo
Si las perdidas aumentaban esfuerzo cortante aumentaba ¿cómo explicarlo?
Relación entre el esfuerzo cortante y la velocidad del
fluido
dv
Flujo laminar:    Ley de Newton de la Viscosidad
dy
(Boussinesq)
Flujo Turbulento:
 yx  (   )
dv =0 Laminar
Esfuerzo Cortante
Turbulento
dy 0 < η < 10000µ Turbulento
Efectos
Efectos debidos a la
debidos a la
viscosidad laminar
viscosidad
turbulenta
La forma mas exitosa fue la teoría de Prandtl- Von Karman, que relaciono  con
el cuadrado de la velocidad

 yx   k 2
v x / y 4
 v
2
x / y 
2 2
Prandtl-Von Karman k  0.4
2
 dv 
 yx   l 2 
 
 l  ky Longitud de mezcla
 dy 
Interacción flujo-pared sólida
Flujo Laminar
A´ B´
Zona del flujo afectada
Dirección
del flujo por el esfuerzo cortante
Capa límite
producido por el
d movimiento del flujo
A B
interactuando con una
pared
Pared sólida

Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

Flujo Turbulento
Capa límite laminar Capa límite turbulenta

Dirección
del flujo d
Subcapa laminar viscosa
d´
Tamaño de esta subcapa
Pared sólida laminar viscosa permite
clasificar el flujo turbulento en:
Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes
Tipo de flujo turbulento

ks < d´
d´ Flujo hidráulicamente liso
B ks

Subcapa laminar Rugosidad media

ks > d´
ks
Flujo hidráulicamente rugoso
d´
Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

El tipo de flujo depende del tamaño relativo entre el espesor de la

subcapa laminar viscosa (δ') y el tamaño de la rugosidad media Ks

Ks Rugosidad [L]

Espesor de subcapa laminar viscosa es

V* = Velocidad de corte
d '
11.6 v o
 v* Magnitud relativa del esfuerzo
v*  cortante en la superficie sòlido-fluido
v = Viscosidad cinemática
Distribución de Esfuerzos

A gh f
En tuberías circulares:
A  ro
2
P  2ro Reemplazando en o 
P l
gh f ro
o  (1) o en las paredes de la tubería
2l

gh f r
Si el análisis se hace para un tubo de fluido de radio r r  (2)
2l

Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

Relación   0 rDistribución lineal de 

(2)/(1) r0
Gráfico Distribución de Velocidades

1. Flujo Laminar
dv
•Partiendo de la ecuación para el esfuerzo cortante en F. L   
dy
•Conociendo la distribución de esfuerzos en tuberías   0 r
r0
•Igualando estas dos ecuaciones y tendiendo en cuenta que y=r
dv  0 r o  r2 
 integrando v   ro  V = 0 implica r = r0
2  ro 
dr  r0

0
Dirección
del flujo r0 r
Distribución
parabólica de la
velocidad

0
Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes
 2. Distribución de Velocidades

Flujo turbulento
Eje de la tubería La presencia

Vx
 de la subcapa
  Zona turbulenta
Dirección   laminar
del flujo   Distribución potencial
  viscosa

modifica la

 
 Zona
distribución de

  de transición


xy
  velocidades
  Distribución logarítmica
 
 

´ Zona Laminar

0
  Distribución lineal
Pared de
la tuberia
Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

11.6 o
d '  v*
v* 

a. Zona Laminar
dv
Se parte de 0   vx v* y Relación lineal
dy Se obtiene 
dv v x v*  respecto a y
Asumiendo 
dy y
2. Flujo turbulento

b.1) Zona de transición - Flujo turbulento hidráulicamente liso:

Eje de la tubería


 Vx  Zona turbulenta
Dirección  
del flujo   Distribución potencial
 
 
 
  Zona
  de transición
 
xy

  Distribución logarítmica
 
 

´ Zona Laminar

0
  Distribución lineal
Pared de
la tuberia
Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

Se parte del esfuerzo cortante para F. turbulento

 dv 
2 l  ky
 yx   l2
 dy 

 
Asumiendo  yx   0 Verificado experimentalmente

v 1 v* y
Desarrollando  ln( ) 5.47
v * 0.4 
2. Flujo turbulento

b.2) Zona de transición - Flujo turbulento hidráulicamente rugoso:

Eje de la tubería


 Vx  Zona turbulenta
Dirección  
del flujo   Distribución potencial
 
 
 
  Zona
  de transición
 
xy
 
  Distribución logarítmica
 
 

´ Zona Laminar

0
  Distribución lineal
Pared de
la tuberia
Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

Cuando la rugosidad afecta la subcapa laminar viscosa,

Nikuradse demostrò:
vx 1 y
 ln( ) 8.48
v * 0.4 k s
 2. Flujo turbulento

c) Zona de turbulenta: Eje de la tubería



 Vx  Zona turbulenta
Dirección  
del flujo   Distribución potencial
 
 
 
  Zona
  de transición
 
xy
 
  Distribución logarítmica
 
 

´ Zona Laminar

0
  Distribución lineal
Pared de
Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes
la tuberia

Algunos autores establecen para esta zona una relación potencial:

n
vx  y 
   n  1/ 7 para Re < 100 000
v x  r0 

Las relaciones para la zona de transición describen muy bien el perfil de velocidades

Distribución de velocidad logarítmica para una tubería con F. Turbulento

Comparación de distribución de Velocidades

Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

Usando las tres ecuaciones anteriores se obtiene el gráfico mostrado arriba.

Conocidas

Ecuaciones para la distribución de esfuerzos

Ecuaciones para la distribución de velocidades

Interacción entre el flujo y la pared sólida

Establecer las ecuaciones de resistencia

para el diseño de tuberías
Ecuaciones para el diseño de tuberías simples

Ecuaciones para la Pérdida de Carga basadas en la física clásica:

1.) Flujo Laminar

dv
A partir del esfuerzo cortante    dr

gh f r
Distribución de esfuerzos r 
2l

 
(Ver desarrollo en
1 ( P1  P2 ) 2
V r0  r Saldarriaga y en
Igualando y desarrollando 2

4 L Apèndice B de Sotelo-
Àvila)

d 4
h h
Q g  Sf
En términos de Q 128 L L

Ecuación de Hagen-Poiseuille
2.) Flujo Turbulento
Algunas de estas ecuaciones incluyen el caso de flujo laminar

Ecuación de Darcy-Weisbach (1845)

Desarrollo por análisis Dimensional y utilizando la física clásica

2
l v
hf  f
d 2g

Varios trabajos tanto teóricos

 Ks 
f  F  Re,  como en laboratorio para
 d  desarrollar expresiones para
el factor f.
f para flujo laminar
También estudiaron el factor f, y su relación con el tipo de flujo

A partir de la ec. de D-W y teniendo en cuenta la relación de Hagen-

Poiseuille:
2
l v
hf  f
d 2g

d 4 hf
Q gS f  Sf
128 l

Se obtiene
64
f 
Re
 f para flujo turbulento hidráulicamente liso
Blassius siguió experimentando y encontró que el tipo de fricción cambiaba si era
flujo hidráulicamente liso o flujo hidráulicamente rugoso

Blassius, encontró empíricamente que para números de Reynolds entre

5000 y 100000 el factor de fricción se podía calcular de acuerdo a la
siguiente ecuación:

0.316
Ecuación de Blassius para flujo hid. liso: f  0.25
Re

Utilizando la ecuación de Darcy-Weisbach y reemplazando en ella la ecuación de

Blassius se llega a:

p 1.75  4.75
 0.241   Q d
0.25

l
La mejor forma de reducir las pérdidas de presión por unidad de longitud es
aumentando el diámetro de la tubería.
 f para flujo turbulento hidráulicamente liso y rugoso

Nikuradse estudió la relación para flujo turbulento hidráulicamente liso y rugoso

1 Re f
Para flujo hidráulicamente liso  2 log 10 ( )
f 2.51

1 3.71D
Para flujo hidráulicamente rugoso  2 log 10( )
f 

Los resultados de estas investigaciones se resumen en el diagrama de

Nikuradse
Diagrama de Nikuradse

f no depende Re y si depende
de la rugosidad del tubo

f depende
solo de Re y
no de la f depende Re y de
rugosidad del la rugosidad del
tubo tubo

Tomado de Hidráulica de tuberías, Juan Saldarriaga. U de los Andes

 c. Flujo turbulento –Ecuaciones de Prandtl – Von Karman

A partir de estos estudios Prandtl y Von Karman completaron las fórmulas

para tuberías de rugosidad comercial…
k s  0.305 d '
a. Flujo hidráulicamente liso:
11.6
Partiendo de Espesor de la subcapa laminar viscosa: d '
v*
vx 1 v y
 ln *  5.47
Con distribución de velocidades en flujo turbulento: v * 0.4 
gh f ro
Esfuerzo cortante en la pared de la tubería:  o 
8 0
2l
Relación entre el factor de fricción y el esfuerzo cortante: f 
v 2
Tubo de Pitot

Dirección
del flujo r Integrando para el diferencial de caudal
x
dy que pasa por el diferencial de área
y

(Perfil de velocidad)
vx 1 y
 ln  8 .48
v k ks

1
Ecuación de Prandtl-Von Kármán:  2 log10 (Re f )  0.8
f
c. Flujo turbulento –Ecuaciones de Prandtl – Von Karman

b. Flujo hidráulicamente rugoso:

-Tamaño de la rugosidad: k s  6.10 d '

De igual manera que para el flujo anterior ahora con la distribución de
velocidades
vx 1 y
 ln( ) 8.48
v * 0.4 k s
Se tiene la Ecuación de : Prandtl-Von Kármán para este tipo de flujo

1 d
 2 log 10( ) 1.14
f ks

A partir de estos estudios Moody realizó el diagrama universal para tuberías

de rugosidad comercial
 DIAGRAMA DE MOODY

Relaciona el diámetro, el factor de fricción y la rugosidad relativa para cada uno

de los materiales con un determinado ks.
Otras investigaciones relacionadas con f

Posterior a estos estudios otros siguieron investigando….

c. Flujo transicional:

-Tamaño de la rugosidad: 0.305d '  k s  6.10 d '

-Ecuación de Colebrook-White:

1  ks 2.51 
 2 log 10   
 
f  3.7d Re f 

Válida para todos los flujos turbulentos en tuberías

Otras investigaciones relacionadas con f

ECUACION DE MOODY:

  k 10 6
 3
1

f  0.0055 1   20000 s
 
  d Re  
 

f = factor de fricción de Darcy

Factor de fricción es una función explícita de la rugosidad relativa y del número de

Reynolds.
Otras investigaciones relacionadas con f

ECUACIÓN DE WOOD:

f  a  b Re c

f = factor de fricción de Darcy

0.225
 ks   ks 
a  0.094   0.53 
d d
0.44
 ks 
0.134
b  88  k 
c  1.62  s 
d d
Factor de fricción de Darcy es función explícita de la rugosidad relativa y del
número de Reynolds. Su forma es radicalmente diferente a la de la ecuación de
Colebrook-White.
 Otras investigaciones relacionadas con f

ECUACIÓN DE BARR

1  A   C 
 2 log10   B    D E  
f   1    2  
A  0.325 * 0.027 log10  1
B  0.93 * 0.0068 log10  1
C  2.95 * 0.29 log10  2
D  0.914 * 0.0052 log10  2

E  1  2
0.44
11 1
0.54

2
1
 hf  3 5
1  0 . 958866 Q 5
.
132786 Q  g  5
1
2   l   hf 5
v  g  k s
 l 
 Otras investigaciones relacionadas con f

ECUACIÓN DE SWAMEE – JAIN

La última ecuación explícita, y por consiguiente la más exitosa, fue desarrollada por los
investigadores Prabhata K. Swamee y Akalank K. Jain en el año de 1976.

A partir de Colebrook – White calculando el factor de fricción para los dos

extremos de flujo turbulento con altos números de Reynolds

Flujo Turbulento
Hidráulicamente Rugoso Flujo Turbulento
Hidráulicamente Liso
0.25
f  2 f 
0.25
  k s    5.74 
2

log 3.7d  log Re 0.9 

     
 ECUACIÓN DE SWAMEE – JAIN

0.25
f  2
  ks 5.74 
log10  3.7 d  Re 0.9 
  

Utilizada para todo el rango de flujo turbulento.

Comparada con la ecuación de Colebrook – White con el fin de establecer
su exactitud. Se encontró que para los rangos:

ks
10   102
6
y 5 103  Re  108
d
Los errores involucrados en el factor de fricción siempre fueron menores que el 1%.
Varios Programas comerciales para el cálculo de redes de distribución de agua potable
utilizan la ecuación de Darcy – Weisbach en conjunto con la ecuación de Swamee – Jain.
Por fuera de la deducción física de la ecuación de Darcy Weisbach se
desarrollaron otras ecuaciones a partir de ensayos prácticos….

Ecuaciones empíricas para el

cálculo de las pérdidas por
fricción
ECUACIÓN DE HAZEN-WILLIAMS (1)

v  0.849C HW R 0.63 0.54

Sf SI

v  1.318C HW R 0.63 0.54

Sf Sistema
Inglés

Donde:
v = Velocidad media de la tubería
R = Radio Hidráulico
Sf = Pérdida de energía por unidad de peso (altura) por unidad de longitud
CHW = Coeficiente de rugosidad de la tubería

Ecuación explícita para la velocidad y por consiguiente para el caudal.

ECUACIÓN DE HAZEN-WILLIAMS (2)

Reemplazando en la ec anterior
d hf
Rh  y Sf 
4 l
Se tiene
0.63 0.54
d hf
v  0.849CHW
2.395 l 0.54
Despejando hf se obtiene:

6.8241l (1’)
hf  1.851 1.167
v1.851

CHW d
 ECUACIÓN DE HAZEN-WILLIAMS (3)

hf
 v1.851
l
Esta relación entre las perdidas y la velocidad coincide con:
Tubería rugosa
Log p
Tubería lisa
l 2 4
1
1.75
1

2
1
1
1

Flujo laminar F. Transición Flujo turbulento Log v

 ECUACIÓN DE HAZEN-WILLIAMS (4)

Comparando H-W en la forma 1’ con D–W se obtiene

14.09
CHW  0.54 0.009 0.081 0.081
f d v Re

‘rugosidad relativa’, depende del tubo, del flujo y del fluido

Límites para el uso de H-W

a) El fluido debe ser agua a temperaturas normales.

b) El diámetro debe ser superior o igual a 3 pulgadas (75mm).
c) La velocidad en las tuberías se debe limitar a 10 pies/s (3m/s).
ECUACIÓN DE HAZEN-WILLIAMS (4)

Volviendo a H-W espresada en terminos de hf

6.8241l
hf  1.851 1.167
v1.851

CHW d

Ventaja sobre la ecuación de Darcy-Weisbach en conjunto con la de

Colebrook White, es que ecuación de Hazen-Williams es explícita
para las pérdidas por fricción. Por lo cual dicha ecuación se
popularizó.

El coeficiente CHW se encuentra tabulado en muchos libros de

ingeniería.
OTRAS FÓRMULAS EMPÍRICAS

V  D x S f
y
Relaciones del tipo

Tubos rugosos en zona turbulenta Agua a T<24ºC

1 1 8g
CHEZY V  CRh S f
2 2 C
f
1/ 6
Rh
al considerar
C
n
MANNING
1 2 / 3 1/ 2 1.486 2 / 3 1 / 2
V  Rh S f V Rh S f
n n
SI Sist. Inglés
OTRAS FÓRMULAS EMPÍRICAS

DARCY-WEISBACH Y COLEBROOK-WHITE
Partiendo de D-W 2
hf  f l v
h f 2 gd d 2g
2
 f
lv Despejando de D-W
lv 1

h f 2 gd f
1 Vd
Reemplazando en C-W y Re en función de la velocidad Re 
f 
 
lv  k 2.51 lv 
 2 log 10  s
 
h f 2 gd  3.7d Vd h f 2 gd 
 
  
Despejando la velocidad se obtiene:
 OTRAS FÓRMULAS EMPÍRICAS

h f 2 gd  k 2.51  l 
V  2 log 10  s
 
l  3.7d d h 2 gd 
 f 
Ecuación explícita para la velocidad en función del
diámetro y las pérdidas, la cual permite obtener d
iterando.
4Q Y reemplazando en la anterior
Despejando V de Q=VA V
d 2 ecuación se obtiene

 2 2 gdh f  k 2.51 l  Ecuación implícita para d, en

4q
 log 10  s
  función del caudal.
d 2
l  3.7d d 2 gdh f 
 