Hidráulica de Tuberías

Capítulo 1

Juan Saldarriaga
Flujo uniforme en tuberías

En flujo uniforme las características del flujo (presión y velocidad

media) permanecen constantes en el espacio y en el tiempo. Por
consiguiente, es el tipo de flujo más fácil de analizar y sus ecuaciones
se utilizan para el diseño de sistemas de tuberías. Como la velocidad no
está cambiando, el fluido no está siendo acelerado. Luego, según la
segunda ley de Newton:

∑ Fx=0
Es decir, existe un equilibrio de fuerzas.

- fuerzas de presión
En el flujo en tuberías actúan tres fuerzas: - fuerzas gravitacionales
- fuerzas de fricción.
Distribución de Esfuerzos en Tuberías Circulares

dx

0

dz

1 0 D ir e c c ió n
d e l f lu jo
z 2
X

Volumen de control para el flujo en una tubería. Se muestran todas las fuerzas
que actúan en un fluido contenido en un volumen de control.
 Si el fluido contenido en el volumen de control mostrado en la
figura no se está acelerando , entonces ∑ Fx = 0. Teniendo en
cuenta ésto, se obtiene la siguiente ecuación:

pA  ( p  dp ) A  W cos   Pdx 0 0

donde:
P = perímetro mojado ( Perímetro interno del tubo ).
W = peso del fluido en el volumen de control.
El peso del volumen de control es:
W  gAdx
Luego:
 dpA  gAdx cos    o Pdx 0
Pero, de acuerdo con la figura:
dx cos  dz
Por consiguiente:
−dpA−ρ gAdz−τ o Pdx=0

τ o Pdx=A (dp+ ρ gdz)

En esta última ecuación dz incluye el signo negativo. Como (p + gz) =
p*, la presión piezométrica, entonces:

τ o Pdx=A (dp∗)

de donde finalmente se obtiene:

A dp*
τ o= (1.1)
P dx
Despejando el esfuerzo cortante en la pared de la tubería se obtiene:
A d p
τ o=
A d
P dx
( p+ ρ gz ) o τ o = ρg
P dx ρg
+z( )
Esta ecuación representa la forma más general de las ecuaciones de
pérdidas por fricción en el flujo a través de una tubería. En esta ecuación
el término:
d p
( )
dx ρg
+z

Representa la pendiente de fricción del flujo en la tubería, Sf .

Reemplazando A/P por el radio hidráulico R, y el producto de la densidad

y la gravedad por el peso específico del fluido se tiene:

τ o =γ RS f (1.2)
 Experimento de Reynolds

Tinta

Boquilla

Válvula

Aparato utilizado por O. Reynolds para establecer el régimen del flujo

en tuberías. Reynolds varió tanto el diámetro de las tuberías como el
líquido que fluía por éstas.
 Filamento de tinta Tinta

a) Caudales bajos: La b) Caudales intermedios: El

tinta no se mezcla filamento de tinta
comenzaba a hacerse
inestable
T in ta m e z c la d a
T in ta m e z c la d a

d) Caudales más altos: c) Caudales altos: Mezcla agua-

Mezcla agua-tinta tinta.
Reynolds define los tipos de flujo de la siguiente forma:

Flujo laminar: Cuando la tinta no se mezcla. El flujo se mueve en

capas sin intercambio de “paquetes” de fluido entre ellas.

Flujo turbulento: Cuando la tinta se mezcla completamente. Se

presenta intercambio de “paquetes” de fluido entre las capas que se
mueven a diferente velocidad. Las partículas no tienen un vector de
velocidad muy definido. Velocidad promedio (flujo cuasi-
permanente).

Flujo en transición: Cuando el filamento de la tinta comienza a

hacerse inestable, con un patrón de oscilación manifiesto.
 NÚMERO DE REYNOLDS

Reynolds fue el primero en introducir el

concepto de similaridad dinámica.

Su experimento lo realizó con diferentes

diámetros de tubería y fluidos, y encontró que
los resultados eran similares. Por lo cual pensó
que las leyes físicas que gobernaban éste
fenómeno eran similares.

Para reproducir su experimento se debían

tener condiciones de velocidad y geometría
similares, medidas en sitios similares, también
las líneas de corriente debían ser similares
como se muestra en la siguiente figura:
 2
2’
1 1’

Flujos similares alrededor de esferas de diferentes tamaños. Los puntos

homólogos 1 y 1´ y 2 y 2´ deben estar sometidos a fuerzas homólogas.

F I= m a
F I= m a
2’ F
2
V
F V L í n e a d e c o r r ie n t e
F P
L ín e a d e c o r r ie n te
F P

Triángulos de fuerza para los puntos homólogos 2 y 2' de los flujos alrededor
de esferas similares.
Teniendo en cuenta los puntos 2 y 2’ de la figura anterior, por semejanza
de triángulos
FI 2 FI 2 '
=
FV 2 FV 2'
Donde las fuerzas inerciales:
dv ds ds
FI = ma ; a= ; v= entonces dt=
dt dt v
2
dv v
a=v α
ds L
3
m=Vρ α ρL
2 2
ma=FI α ρv L (1.3)
y las fuerzas viscosas:

dv
FV = τA en donde τ=μ
dy

dv
FV =μ A α μvL (1.4)
dy
FI ρv 2 L2 vLρ vL
= = = (1.5)
FV μvL μ υ

vLρ vLρ
( ) ( )
μ 2
=
μ 2'
 ρdv
Número de
Reynolds μ

Reynolds encontró que:

Re < 2000 Flujo Laminar

2000 < Re < 5000 Flujo en Transición

Re > 5000 Flujo Turbulento

Explicó numéricamente el cambio entre los diferentes regímenes de

flujo.
Segundo Experimento de Reynolds
C o n tr o l ( v á lv u la )
L>250d l a g u a s a b a jo

d V

1

h

2

Aparato diseñado por Reynolds para estudiar la caída de presión

por unidad de longitud a lo largo de una tubería en función del tipo
de flujo.
Resultados del segundo experimento de Reynolds:

T u b e ría ru g o s a
Log p
T u b e r ía lis a
l 2 4
1
1 .7 5
1

2
1
1
1

F lu jo la m in a r F . T r a n s ic ió n F lu jo tu r b u le n to Log v

Gráfica logarítmica de los resultados del segundo experimento de

Reynolds
 De la figura anterior se pueden obtener las siguientes
conclusiones:

a) La línea 0-1 presentaba una pendiente de 1 a 1, lo cual implicaba una

variación directa del logaritmo de las pérdidas de presión por unidad
de longitud con respecto al logaritmo de la velocidad. Esto era
válido para el flujo laminar.

b) Si el experimento se hacía abriendo la válvula se alcanzaba a tener

flujo laminar hasta el punto 2 (Re=2500-4000).
Continuación conclusiones:

c) Si el experimento se hacía cerrando la válvula, el flujo laminar se

restablecía en el punto 1 (Re=2200).

d) En el punto 3 (Re>5000) la variación log(∆p/l) vs. log(v) volvía

a ser aproximadamente lineal con pendientes desde 1.75 para
tubos muy lisos, hasta 2.0 para tubos muy rugosos.

e) La zona de transición se obtenía para 2000≤Re≤5000. En ésta la

variación log(∆p/l) vs. log(v) era muy compleja.
 Viscosidad de remolino. Longitud de mezcla
Boussines introdujo el concepto de viscosidad de remolino, teniendo
en cuenta la siguiente figura:
Y
C a m p o d e f lu jo

 v
v x x
  y d y
dy v x





X
In te r c a m b io d e p a q u e te s
Z d e m ó le c u la s
Dos placas de fluido moviéndose a diferente velocidad
dentro de un campo de flujo en una tubería.

Supuso que el intercambio de paquetes de moléculas entre capas de flujo

Turbulento añade momentum acelerando las capas o resta momentum
frenando las capas, esto produce un efecto similar al de viscosidad dinámica µ.
El intercambio de momentum producía un nuevo
esfuerzo cortante, causado por la turbulencia del
flujo:
τyxT = Esfuerzo Cortante Turbulento
δv x
τ yxT =η
δy
δv x δv x
τ yx =μ +η
δy δy
δv x
τ yx =( μ+η ) (1.6)
δy
 → propiedad del flujo
η → propiedad del fluido
Por similaridad con la viscosidad cinemática:
η
ε= (1.7)
ρ
-Para flujo laminar:  = 0.
-Para flujo turbulento: 0 ≤  ≤ 10000.
 Esfuerzo de Reynolds

Reynolds se basó es sus observaciones de flujo turbulento en sus

experimentos y encontró que la trayectoria seguida por una
partícula individual era aleatoria y que ninguna otra partícula
seguiría una trayectoria exactamente igual a la de otra, como se
muestra en la siguiente figura:

Dirección
2 3 5
del flujo
1
VX
Y  4

Trayectoria de la partícula
X
individual

Trayectoria seguida por una partícula de fluido

en un campo de flujo turbulento
 Esfuerzo Cortante de Reynolds
Con el fin de obtener una mejor ecuación que definiera el aumento del
esfuerzo cortante cuando el flujo era turbulento, Reynolds supuso lo
siguiente:
1. La velocidad en el sentido principal del flujo (eje x) está compuesta por
una velocidad media y una variación aleatoria de ésta:
'
v x =v x + v x
Trayectoria
seguida por
una partícula
de un fluido
En un
campo de
Y flujo
turbulento

X
 Esfuerzo de Reynolds
2. La velocidad en el sentido vertical del flujo (eje y) está compuesta
por:
v y=v y+v ' y

X
Trayectoria seguida por una partícula de
fluido en un campo de flujo turbulento
 Interacción flujo-pared sólida
Siempre que un fluido en movimiento interactúa con una pared sólida,
el esfuerzo cortante que se genera afecta una zona de dicho flujo. Esta
zona recibe el nombre de la capa límite, la cual puede ser laminar o
turbulenta.

A´ B´
Dirección
del flujo Capa límite

A B

Pared sólida

La pared sólida afecta la distribución de velocidades del flujo. La velocidad

es cero en la superficie sólida y crece a medida que el flujo se aleja de ésta.
En el flujo turbulento la superficie sólida impide que cerca a ella ocurran
las vibraciones de v'y en forma libre, razón por la cual se genera una zona
de flujo laminar. Esta se conoce con el nombre de la subcapa laminar
viscosa.

C a p a lím it e la m in a r C a p a lím it e tu r b u le n ta

D ir e c c ió n
d e l flu jo 
S u b c a p a la m in a r v is c o s a
´

P a r e d s ó lid a

Desarrollo de una capa límite turbulenta mostrando la subcapa

laminar que se genera cerca a la superficie.
El espesor de la subcapa laminar es mucho menor que el de la capa límite
(’). 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).

k s < ´
´ F lu jo h id r á u lic a m e n t e lis o
B ks

S u b c a p a la m in a r R u g o s id a d m e d ia

k s > ´
ks
F lu jo h id r á u lic a m e n t e r u g o s o
´

Flujos hidráulicamente lisos e hidráulicamente rugosos.

 Distribución de Esfuerzos en Tuberías Circulares
Para conocer el esfuerzo cortante en las paredes de la tubería (τo) hay que
conocer la caída en la presión piezométrica en una determinada longitud
de la misma tubería y la geometría de ésta.
Para tuberías circulares:
A=πr P=2πr o
o2
Se sabe: A dp*
τ 0= (1.11)
P dx

Luego: πr 2
o Δp* r o Δp*
τ o= =
2 πr o Δx 2 Δx

Adicionalmente, si se aplica un razonamiento similar a un "tubo" de fluido

de radio r<ro se obtendría:
r Δp*
τr =
2 Δx (1.9’)
Dividiendo esta última ecuación por la ecuación para τo se obtiene lo
siguiente:
τr r
=
τ o ro

Esta última ecuación indica una variación lineal de τ con respecto a r, tal
como se muestra en la siguiente figura:
τo

Distribución lineal del esfuerzo cortante en tuberías

circulares
Utilizando la ecuación (1.11) se obtiene:

A ρ gh f
τ o= (1.12)
P L
donde:
hf = pérdidas por fricción
L = longitud del tramo de tubería
Nuevamente, para tubos circulares:
ρ gh f r o
τ o=
2L
Esta última ecuación establece una relación directa entre el esfuerzo cortante
en la pared de una tubería y la caída en la altura piezométrica hf . Luego, es
un primer paso para la deducción de una ecuación de diseño.
Distribución de Velocidades
(Tuberías Circulares)

Distribución de velocidades
en secciones transversales de
tuberías circulares.

Diferencias entre los tipos de

flujo, para encontrar las
ecuaciones de resistencia
fluida.

Aplicación: Diseño de
sistemas de Tuberías.
1. Flujo Laminar
En el caso de flujo laminar en tuberías el
esfuerzo cortante está definido por la
ecuación de Newton para fluidos viscosos:

dv
τ o =μ
dr
Para tuberías circulares, de acuerdo con la
distribución de esfuerzos, se tiene lo
siguiente:
dv r
τ r =μ =τ o
dr ro
Donde: r < ro
De esta última ecuación se obtiene:

dv τ o r
=
dr μ r o
Integrando con respecto al radio r se obtiene la distribución de velocidades:

τo r
v =∫ dv= ∫0 rdr
μr o

τo r 2
v=
μr o 2 ( )
+cte

Para evaluar la constante de integración se utiliza el hecho de que cuando

r = ro, entonces v = 0; luego:
τo
cte=− ro
2μ
Finalmente:
τ o r2
v=
(
2 μ ro
−r o
) (1.13)

La cual es una distribución parabólica de la velocidad. En flujo laminar la

velocidad sigue esta distribución con su máximo en el centro y su mínimo
(v = 0) en las paredes internas de la tubería.

0
D ir e c c ió n
d e l flu jo r0 r

0

Distribución de velocidades para flujo laminar en

una tubería de sección circular.
2. Flujo Turbulento
La presencia de esfuerzos cortantes en las fronteras fluidos-sólidos y entre
las diferentes capas del fluido afecta la distribución de velocidades que en
principio, debería ser uniforme. En flujo turbulento la presencia de la
subcapa laminar viscosa modifica aún más dicha distribución.
 Eje de la tubería


 Vx  Zona turbulenta
Dirección  
del flujo   Distribución potencial
 
 
 
  Zona
  de transición
 
xy
 
  Distribución logarítmica
 
 
´


0 

Distribución lineal Zona Laminar
Pared de
la tuberia

Distribución de esfuerzos y velocidades para flujo turbulento en

una tubería de sección circular. Se muestra media tubería.
Subcapa laminar viscosa (flujo laminar):
El esfuerzo debe estar gobernado por la ecuación de viscosidad de Newton:

dv x
τ o =μ
dy
lo cual implica que los gradientes de velocidad en esta zona deben ser muy
grandes. Igualmente, por ser de espesor tan pequeño, en la subcapa laminar
se cumple que:
dv x vx
≈
dy y
Luego: vx
τ o =μ
y
Ahora, dividiendo por la densidad:

τo μ vx vx
= =υ
ρ ρ y y
Donde: υ = viscosidad cinemática

La raíz cuadrada del término τo /ρ, localizado en la parte izquierda de

esta última ecuación, tiene las siguientes dimensiones:

[√ ]τo
ρ
=
m
s
(dimensiones de velocidad )
Por definición:

Luego,
√
2
τo
ρ
=v *
vx
v =υ ¿

y
De donde se obtiene la siguiente ecuación:

vx v* y
= (1.14)
v* υ

La ecuación (1.14) indica que la velocidad sigue una distribución lineal con
respecto a y, siempre y cuando y ≤ δ’, el espesor de la subcapa laminar
viscosa. Esta ecuación ha probado ser válida hasta el siguiente límite:
vx
=11. 6
v*
lo cual implica que:

v¿ y
≤11. 6
υ

Este último resultado permite establecer una ecuación para el cálculo del
espesor de la subcapa laminar viscosa y así establecer claramente la
diferencia entre flujos hidráulicamente lisos e hidráulicamente rugosos. En
el límite y = ’, lo cual implica que:

11.6υ (1.15)
δ '=
v*
Simplificando la expresión, se obtiene:
dv x 1 dy
=
v* k y

Con el fin de encontrar la distribución de velocidad, se debe integrar la

ecuación anterior:
dv x 1 dy
∫ v =k∫ y
*
Lo que resulta en:
vx 1
= ln y +c (1.16)
v¿ k
utilizando la ecuación (1.16): vl 1
= ln δ ' +c
v* k

Se llega a:
vl 1 11. 6 υ
= ln +c (1.18)
v* k v*

Luego, reemplazando la ecuación ( 1.15 ) en la ecuación ( 1.16 ):

1 11. 6 υ
11.6= ln +c (1.19)
k v*
Despejando la constante de integración se obtiene:

1 11.6 ν
c=11.6− ln
k v*

Luego, reemplazando la ecuación ( 1.16 ) se obtiene:

vx 1 1 11. 6 υ
= ln y +11. 6− ln
v¿ k k v*
Esta última ecuación se puede transformar en:

vx 1 v* y
= ln +A
v* k υ
Donde: 1
A=11.6− ln ( 11. 6 )
k
A=5 . 47
Luego, la distribución de velocidades en la zona de transición está definida
por la siguiente ecuación:

vx 1 v* y
= ln +5 . 47 (1.20)
v * 0. 4 υ
La ecuación (1.20) puede ser graficada en forma lineal
escogiendo unos ejes de coordenadas apropiados. Esto se hace
de la siguiente forma:
vx 1 v* y
= ln +5 . 47
v * 0. 4 υ
vx
=Y
v*
1
=m
0.4 Y =mX +b
v* y
ln =X
υ
5 . 47=b
La gráfica de la línea recta sería la siguiente:

vx
Y
v¿
1
m 0.4

b 5.47
v¿ y
X ln
υ
 Efecto de la rugosidad de la pared interna de la
tubería en la distribución de velocidades.
Vx Vx
V* C a p a v is c o s a Z o n a t u r b u le n ta V*
d e p a re d

B Corrimiento hacia la
vx 1 vy derecha. (Nikuradse)
 5 .5  ln
v k v

R a n g o d e d a to s e x p e r im e n ta le s v ks
ln  7 .5 2 5
v
5 .4 7 v x
v 

v  y
v 5 .4 7

vy  vy 
ln   ln  
 v   v 

a) Flujo turbulento hidráulicamente b) Flujo turbulento hidráulicamente

Liso. (Ecuación 1.19) rugoso comparado con el liso.
Para el FHR J. Nikuradse demostró que la distribución de
velocidades seguía siendo logarítmica pero que dependía de
(ks) de la tubería.

Como se muestra en la figura anterior (b), la relación

entre ln v * y y vx/v* se corría hacia la derecha para lo que
υ
Nikuradse encontró:

1 ks v*
ΔB= ln −3 . 01 (1.21)
0.4 υ
Es claro que el valor de la relación vx/v* para un FHR debe
ser un B menor que el valor de la misma relación en un
FHL. Es decir:
vx 1 v* y
= ln +5 . 47−ΔB
v * 0. 4 υ

vx 1 y
= ln +8 . 48 (1.22)
v* 0. 4 k s
 Zona turbulenta
Algunos autores afirman que la distribución de velocidad en esta zona es
similar a la de la zona de transición, especialmente en el flujo en tuberías.
Por lo cual la distribución logarítmica encontrada antes es aplicable a esta
zona.

Otros autores dicen que la distribución es potencial, regida por la siguiente

ecuación:
Donde:
vx y n
v̄ Q
___
vx
=
( )
ro
x = velocidad media =
A
ro = radio total de la tubería

n = 1/7 si Re < 100.000

Si el número de Reynolds ( Re ) aumenta, el exponente n empieza a disminuir. Esta
distribución de velocidades se conoce como la Ley de la Potencia 1/7.
La diferencia en la presión piezométrica puede ser reemplazada por la
diferencia en altura piezométrica con lo cual se obtiene la siguiente
expresión:
4
πd
Q= ρgΔh
128 μL
de donde se obtiene la siguiente ecuación:

4
πd Δh (1.33)
Q= ρg
128 μ L

La relación entre la caída piezométrica y la longitud es la pendiente de

fricción:
Δh
=S f
L