Capitulo 2
Capitulo 2
a) DISEÑO EN SÍ DE LA TUBERÍA:
Se conoce el caudal demandado, la potencia
disponible y algunas características de la tubería
(longitud, accesorios). Se desconoce el diámetro necesario.
En cuanto al material de la tubería usualmente se tienen sólo
2 o 3 alternativas.
Variables conocidas Incógnita
l, Skm, QD, H d
r, m, g, (ks)
Tipos de Problemas en Hidráulica de
Conductos a Presión
a) CALIBRACIÓN DE LA TUBERÍA:
Se conoce el caudal demandado, la caída en la altura piezométrica que ocurre
entre la entrada y la salida de la tubería, algunas de las características de la
tubería y las propiedades del fluido. Se desconoce la rugosidad absoluta de la
tubería que produce esa caída en la presión piezométrica para el caudal medido.
En campo el caudal usualmente se mide a través de caudalímetros acústicos y la
caída de presión utilizando transductores de presión.
1. Ecuación de Fricción:
Utilizando la ecuación de Colebrook-White en conjunto
con la ecuación de Darcy-Weisbach:
l v2
hf = f (1.38)
d 2g
1 ks 2.51
= 2log 10 +
3.7d Re f (1.69)
f
2. Ecuación de Conservación de la Energía
Si se plantea una ecuación de energía entre un embalse (bomba) y un
punto en la tubería se obtendrá lo siguiente:
Luego: v p2
h1 + z1 = + z2 + + h f + hm
22
2g ρg
Si el punto 2 es la salida:
p2 v2
=0 ; 2 =0
ρg 2g
Es claro que para el punto 1, h1 + z1 = H, donde H es la altura del nivel
de la superficie del tanque con respecto al Datum. Luego:
H = z2 + h f + hm
De la última expresión se puede obtener la siguiente ecuación, la cual
describe las pérdidas por fricción en función de las otras variables:
v
h f = H z2 km 22
(2.1)
2g
Utilizando la ecuación de Darcy-Weisbach Ecuación (1.38), que también
predice las pérdidas por fricción, se puede despejar el factor de fricción f:
l v 2
h f d 2g
hf = f f =
d 2g l v 2
h f d 2g 1 l v
f = =
l v f h f d 2g
(2.2)
Igualando la ecuación (2.2) con la ecuación (1.69) (Colebrook-White)
se obtiene:
1 l v
= (2.2)
f h f d 2g
1 Ks 2.51
= 2log 10 + (1.69)
3.7d Re f
f
l v K
= 2log 10 s
+
2.51v l
h f d 2g 3.7d Re h f d2g
El número de Reynolds en esta última ecuación puede reemplazarse
por:
vd
Re =
υ
Por consiguiente, se obtiene la siguiente expresión:
v l k 2.51υv l
= 2log 10 s
+
2 gdh f 3.7d vd 2 gdh f
Finalmente, despejando la velocidad se encuentra una ecuación
explícita para esa variable:
2 2 gdh f k 2.51 υ l
v= log 10 s +
l 3.7d d 2 gdh f (2.3)
Esta última ecuación es la base para la solución de los tres tipos de
problemas relacionados con tuberías simples mencionados anteriormente.
Para el caudal, se multiplica la ecuación 2.3 y se obtiene la siguiente
expresión:
2 2 gdh f k 2.51υ L
Q= Alog 10 s
+
L 3.7d d 2 gdh f
a) Comprobación de Diseño en Tuberías Simples
Variables conocidas
• Características de la tubería: longitud
(l), diámetro (d), rugosidad absoluta
(ks) y coeficientes de pérdidas
menores (km).
• Propiedades del fluido: densidad (ρ)
y viscosidad dinámica (µ).
• potencia (P) y eficiencia de la
bomba(η).
Variables desconocidas
• velocidad (v) y caudal (Q).
INICIO
Q = Vi A
Leer d, ks, H, E, Skm, r, m, z2, l
Imprima Q
Suponer hf1 = H- z2
PARE
Calcular ks/d
Diagrama de flujo 1:
│hfi – hfi-1 │ ≤ E Comprobación de diseño de tuberías
simples
V V1 = Vc = Vconvergencia
V3
V5 V7 V9
Vc
V6 V8
V4
V2
0 1 2 3 4 5 6 7 8 9 10
Número de iteraciones
Ejemplo 1.
Se desea calcular el caudal de agua que puede ser movido a través de una
tubería de PVC, de 300 mm de diámetro nominal y 730 m de longitud, que
conecta dos tanques de abastecimiento de agua potable con una diferencia
de nivel de 43.5 m. El diámetro real de la tubería es de 293 mm y su
rugosidad absoluta es de 1.5 x 10-6 m. Todos los accesorios que forman
parte del sistema, incluyendo la entrada y la salida, implican un coeficiente
global de pérdidas menores km de 11.8. El agua se encuentra a 20ºC.
ρ= 998 . 2 kg /m 3
−3 −6 2
μ= 1 .005 × 10 Pa⋅ s ⇒ v = 1 . 007× 10 m /s
Con los datos anteriores se puede seguir el procedimiento del Diagrama de
Flujo 1; los resultados de las iteraciones se muestran en la siguiente tabla.
Este procedimiento es fácilmente programable
H ks/d hf v hf i+1 ∑hm
(m) (-) (m) (m/s) (m)
18,928849
43,5 5,12E-06 43,5 5,6101 24,5711 1
10,186575
43,5 5,12E-06 24,5711 4,1155 33,3135 7
14,175669
43,5 5,12E-06 33,3135 4,8549 29,3241 6
12,343460
43,5 5,12E-06 29,3241 4,5303 31,1564 4
43,5 5,12E-06 31,1564 4,6818 30,3171 13,182832
12,797675
43,5 5,12E-06 30,3171 4,6129 30,7021 2
12,974171
43,5 5,12E-06 30,7021 4,6446 30,5256 7
12,893289
43,5 5,12E-06 30,5256 4,6301 30,6065 9
12,930631
43,5 5,12E-06 30,6065 4,6368 30,5694 4
12,913347
43,5 5,12E-06 30,5694 4,6337 30,5864 3
12,921151
43,5 5,12E-06 30,5864 4,6351 30,5787 6
12,917806
43,5 5,12E-06 30,5787 4,6345 30,5822 6
43,5 5,12E-06 30,5822 4,6348 30,5806 12,919479
12,918921
43,5 5,12E-06 30,5806 4,6347 30,5813 5
12,918921
43,5 5,12E-06 30,5813 4,6347 30,581 5
12,918921
43,5 5,12E-06 30,581 4,6347 30,5812 5
12,918921
43,5 5,12E-06 30,5812 4,6347 30,5811 5
Los resultados de este ejemplo se pueden resumir así:
f = 0.011211
h f = 30.58
hm = H h f = 43.5m 30.58 m
hm = 12.92m
v = 4.635m / s
Finalmente se calcula el caudal que pasa a través de la tubería:
Q= vA
m π 2
Q= 4 . 635 0 .293 m
s4
m3
Q= 0 . 3124
s
l
Q= 312 . 4
s
b) Cálculo de la Potencia Requerida
En este caso todas las características de la tubería son conocidas, al igual
que las del fluido. Se conoce un caudal demandado y se pregunta por la
potencia requerida, ya sea de origen gravitacional o mecánico.
Para poder resolver el problema de la potencia requerida es necesario
utilizar un método numérico con el fin de poder averiguar el valor del factor
de fricción f de Darcy en la ecuación no explícita de Colebrook-White:
1 ks 2.51
= 2log 10 +
3.7d Re f (1.69)
f
Para que este método pueda ser aplicado la función no explícita debe ser
de la siguiente forma:
x = g(x)
SI Imprima fi+1
Re < 2200 f = 64/Re
NO
PARE
PARE
f1 = semilla de f
i=1
−2
[
f i 1 = − 2 lo g 10
ks 2 . 51
3 . 7d Re f ] Diagrama de flujo 2a.
i=i+1
Cálculo del factor de
?
fricción f por el método
│fi+1 - fi │≤ Ɛ NO fi-+1 = fi de iteración de un punto.
SI
f x g(x) f f x g(x) f
Re=20000
Re=20000 0
31.622.77 479.681.74 0.04346039 316.227.76 674.548.77 0.0219772
0.001 7 1 1 0.001 6 7 5
0.0434603 47.968.17 640.266.05 0.02439377 0.0219772 674.548.72 0.0160067
9 4 7 7 5 5 79.040.267 5
0.0243937 64.026.60 616.125.70 0.02634276 0.0160067 790.402.70 779.772.59 0.0164461
8 1 6 4 5 7 9 4
0.0263427 61.612.57 619.353.25 0.02606892 0.0164461 779.772.66 780.695.51 0.0164072
6 5 2 7 4 8 2 8
0.0260689 61.935.32 618.914.73 0.02610588 0.0164072 780.695.55 780.614.99 0.0164106
3 2 3 1 8 2 6 7
0.0261058 61.891.47 618.974.17 0.02610086 0.0164106 780.614.91 780.622.02 0.0164103
8 5 8 7 7 3 9 7
0.0261008 61.897.41 618.966.12 0.02610154 0.0164103 780.622.04 780.621.40
7 4 3 7 7 8 6 0.0164104
0.0261015 61.896.60 618.967.21 0.02610145 780.621.33 780.621.46 0.0164103
5 8 6 4 0.0164104 5 9 9
0.0261014 61.896.72 618.967.05 0.02610146 0.0164103 780.621.57 780.621.44
5 7 5 8 9 2 8 0.0164104
0.0261014 61.896.70 618.967.08 0.02610146 780.621.33 780.621.46 0.0164103
7 3 8 5 0.0164104 5 9 9
x = g(x)
debe cumplir 3 condiciones especiales para que exista convergencia.
• La primera condición es que exista un intervalo I = (a,b) tal que para
todo x perteneciente a I, la función g(x) esté definida y pertenezca a I, lo
cual significa que g(x) se aplica a sí misma.
• La segunda condición es que la función de iteración g(x) sea continua
en I.
• La tercera condición de convergencia es que g(x) sea diferenciable en
I y que la pendiente de g(x) sea siempre menor que 1 y mayor que -1.
La ecuación de Colebrook-White cumple con las tres.
Cálculodel factor f
MétododeNewton-Raphson
0.06
0.05
0.04
f, g(f)
0.03
0.02
0.01
0
0.01 0.02 0.03 0.04 0.05 0.06
f
f =f Re =20000
En la figura anterior se esquematizó el proceso de convergencia de este
método. En este caso se utilizó una tubería con rugosidad relativa (ks/d)
de 0.0001 y un número de Reynolds de 20.000. El valor semilla para f fue
0.002. A pesar de que este valor estaba bastante lejos del valor real (f =
0.0261) el método convergió muy rápidamente; sus ventajas sobre el
método anterior saltan a la vista.
f1 = semilla de f
? f = 64/Re
Re < 2200 SI
NO
xi = 1/f11/2
FIN
i=1
−2
[
F x i = − 2 lo g 10
ks 2 . 51
3 .7d Re f ] i=i+1
f = 1/x2i+1
[ ]
2 . 51
−2
F ' xi = [ Ln 10 ] ks
Re
2 . 51 x i
xi-+1 = xi
3 . 7d Re
NO FIN
F xi − x i ?
x i 1= x i−
F' x i − 1 │Xi+1 – Xi │≤ Ɛ
SI
Una vez se pueda calcular el
valor del factor de fricción de
Darcy f en la ecuación de
Colebrook-White el cálculo de la
potencia requerida es bastante
sencillo. En el Diagrama de
Flujo No.3 se esquematiza dicho
procedimiento de cálculo.
INICIO
Calcular H total
Leer Q, d, ks, Skm, r, m, , l, z2
1
Pot = ρ QgH
Calcular V = Q/A η
FIN
Calcular Re y ks/d
L = 970 m
Q = 0.042 m3/s
z2 = 16 m
d = 150mm
= 1.14 x 10-6 m2/s
Para una tubería de PVC de 150mm de diámetro el área es:
A = 1.77x10-2m2
v2
• Cálculo de las pérdidas menores: ∑ hm= ∑ k m 2g
= 2 . 706 m
vd
Re= = 312725
• Cálculo del Re y la rugosidad relativa: υ
ks
= 1 × 10 − 5
d
• Cálculo del factor de fricción mediante el método de Newton
(Diagrama de flujo No. 2a):
f x g(x) f
Luego f = 0.01446
Cálculo de las pérdidas por fricción utilizando la ecuación de Darcy-
Weisbach:
2
l v
hf = f = 26.92m
d 2g
Cálculo de la altura total que debe ser producida por la bomba:
H = z2 + h f + hm = 45.62
Cálculo de la potencia bomba:
Se supone:
η = 75
ρQgH
Pot = = 18.78kW
η
c) Diseño de Tuberías Simples
∑ h m ≤ 0 . 3h f
(2.4)
Esta última restricción en la práctica resulta ser irrelevante ya que en
la gran mayoría de los sistemas de tuberías esto se cumple con
facilidad. Para diseñar un sistema con altas pérdidas, menores, como
en el caso de la tubería de succión de una bomba, se debe seguir un
algoritmo diferente.
INICIO
di+1 = di+d
Q = vA
NO
? NO SI
Q ≥ Qd ? Siguiente d
Diám. comercial
comercial
SI
Calcular hf en la ecuación 2.1’
NO
? Calcular V en la
│hfi – hfi-1 │ ≤ E Q = VA
ecuación 2.3
SI
? NO ? SI
Q ≥ Qd Diám. comercial Siguiente d
comercial
SI NO
Imprimir di+1 di+1 = di+d hf = H – z2
FIN
Ejemplo 3
La tubería de descarga de la planta de tratamiento de aguas residuales del
municipio de Ubaté tiene una longitud de 150 m y por ella debe pasar un caudal
máximo de 120 l/s. La altura mínima de operación es 2.2 m y en la tubería se
tienen pérdidas menores por entrada (km = 0.5), por un codo (km = 0.8), por
uniones ( km = 10 x 0.1), y por salida (km = 1.0).
Ejemplo 4
Suponiendo que la planta de Ubaté se localiza a sólo 15 m del río Suta,
sitio de descarga, la tubería tendría un total de 17 m de longitud. Si las
uniones fueran roscadas, las pérdidas menores serían: entrada (km = 0.5),
por un codo (km = 0.8), por uniones (Skm = 4 x 0.5), y por salida
(km = 1.0). Calcular el diámetro de la tubería comercial en PVC requerido
para la descarga.
Los datos del problema son:
l = 17 m H = 2.2 m
ks = 0.00015 m Skm = 0.5 + 0.8 + 4x0.5 + 1.0 = 4.30
QD = 0.12 m3/s (14ºC) = 1.17x10-6 m2/s
Siguiendo la metodología de Darcy-Weisbach, la cuación de Colebrook-
White y el diagrama de flujo 4, se tienen los siguientes resultados:
hf d v A Q Q Qd ∑hm
(m) (m) (m/s) (m2) (m3/s) (si o no) (m)
2,20 0,100 3,360 0,008 0,026 no 2,474
2,20 0,150 4,345 0,018 0,077 no 4,137
2,20 0,200 5,204 0,031 0,163 si 5,934
-3,73
El proceso que permite tener en cuenta sistemas con pérdidas menores altas
fue desarrollado por Saldarriaga y Ferrer (1989) y modificado por Camacho
(1990). Consiste en definir una “velocidad de pérdida”, la cual, en esencia,
es la velocidad que haría que la sumatoria de las pérdidas menores fuera igual
a la altura disponible:
∑ hm= H
v2
Además: ∑ h m = ∑ k m 2g
Mediante las dos ecuaciones anteriores se obtiene el siguiente resultado
para la “velocidad” de pérdida:
v 2pH
=
2g km
2 2 gH
v =
km
p
2 gH
vp = (2.5)
km
Si en alguna iteración la vi es mayor que la vp, ésto quiere decir que la
velocidad vi implica unas pérdidas menores mayores a la altura disponible, lo
cual es físicamente imposible. Si esto sucede, se debe limitar la altura
disponible para ser perdida por fricción, dentro del procedimiento de diseño.
2g(H h f )
vp =
k m
(2.6)
INICIO
Leer Qd, ks, d, , z2, E, Ev, L, h Diagrama de flujo 5: Diseño de tuberías
Suponer hf = H – z2 simples con altas pérdidas menores
Suponer di “pequeño”
Q = vA di+1 = di+d
NO NO SI
? ? Siguiente d
Q ≥ Qd d comercial comercial
SI
Calcular Vp en la ecuación 2.5
NO Procedimiento pérdidas
Vi <? Vp menores altas A
SI
Calcular hf en la ecuación 2.1’
NO Calcular v en la Q = vA
│hfi – h?fi-1 │ ≤ E ecuación 2.3
SI
NO ? SI Siguiente d
Q ≥? Qd d comercial comercial
SI NO
Imprimir di+1 di+1 = di+d hf = H – z2
FIN
A
Asignar a hf un valor pequeño
di+1 = di+d
Q = vA
NO
? NO SI Siguiente d
Q ≥ Qd ?
Diám. comercial comercial
SI
Calcular vp en la ecuación 2.6
? SI
|vi - vp| < Ev Imprimir d FIN
NO
? SI
vi < vp
NO
hfi = hfi -1 - h hfi = hfi -1 + h
Suponer di “pequeño”
Ejemplo 4 (Continuación)
La velocidad de pérdida inicial se calcula como:
2 gH
vp =
Como: km
H= 2 . 2m ∑ k m= 4 . 3
Luego:
2 9.81 2.2
vp = m/ s
4.3
v p = 3 . 168 m /s
Con ayuda del Diagrama de Flujo No. 5 se obtienen los siguientes resultados para
el diseño de la tubería de descarga de la planta de tratamiento del municipio de
Ubaté, cuando su longitud se reduce a 17 metros en total.
hf d v A Q Q
Q(si od ∑hm vp
(m) (m) (m/s) (m2) (m3/s) no) (m) (m/s)
2,200 0,100 3,360 0,008 0,026 no 2,474 3,168
2,200 0,150 4,345 0,018 0,077 no 4,137 3,168
2,200 0,200 5,204 0,031 0,163 si 5,934 3,168
0,500 0,200 2,452 0,031 0,077 no 1,318 2,785
0,500 0,250 2,821 0,049 0,138 si 1,745 2,785
0,480 0,200 2,402 0,031 0,075 no 1,264 2,801
0,480 0,250 2,763 0,049 0,136 si 1,674 2,801
0,490 0,200 2,427 0,031 0,076 no 1,291 2,793
0,490 0,250 2,793 0,049 0,137 si 1,709 2,793
En la última iteración se tiene lo siguiente:
vi = v p
2.793 m/s = 2.799 m/s
H = hf + ∑hm
d = 250mm
hf = 0.49 metros H = 0.49 m + 1.709 m
∑hm = 1.709 metros H = 2.199 m ≈ 2.20 m
Q = 137 Lt/s
hf d real v A Q Q
Q(si od ∑hm
(m) (m) (m/s) (m2) (m3/s) no) (m)
15,20 0,08042 1,200 0,005 0,006 no 0,873
15,20 0,10342 1,421 0,008 0,012 no 1,225
15,20 0,15222 1,837 0,018 0,033 no 2,046
15,20 0,19848 2,185 0,031 0,068 si 2,897
12,30 0,19848 1,946 0,031 0,060 si 2,296
12,90 0,19848 1,997 0,031 0,062 si 2,419
12,78 0,19848 1,987 0,031 0,061 si 2,394
12,81 0,19848 1,989 0,031 0,062 si 2,399
12,80 0,19848 1,988 0,031 0,062 si 2,398
12,80 0,19848 1,989 0,031 0,062 si 2,398
Los anteriores resultados indican que en el caso del PVC es necesario
colocar una tubería de 200 mm de diámetro nominal.
= 33.2 m ks = 0.0000015 m
= 1.007 x 10-6 m2/s Q = 270 l/s
d = 0.3048 m l = 365 m
Re = 1120026,85 Ks/d=0.000005
Cálculo del Factor de Fricción
f ks d Q Re f
(m) (m) (m3/s)
0,02 0,0000015 0,3048 0,27 1120026,85 0,01101037
0,01101037 0,0000015 0,3048 0,27 1120026,85 0,01159084
0,01159084 0,0000015 0,3048 0,27 1120026,85 0,01153863
0,01153863 0,0000015 0,3048 0,27 1120026,85 0,0115432
0,0115432 0,0000015 0,3048 0,27 1120026,85 0,0115428
0,0115428 0,0000015 0,3048 0,27 1120026,85 0,01154284
0,01154284 0,0000015 0,3048 0,27 1120026,85 0,01154283
f = 0.0115
Las pérdidas por fricción son:
l v2
hf = f
d 2g
365 3 . 820 2
h f = 0 . 0115
0 .3 2 × 9 . 81
h f = 10 . 406 m
Luego las pérdidas menores son:
∑ h m = H− hf
∑ h m = 33 .2m − 10 . 406 m
∑ h m = 22 .794 m
El coeficiente global de pérdidas menores es:
v2
∑ h m = ∑ k m 2g
2g ∑ h m
∑ k m= v 2
2gπ2 d 4 ∑ h m
∑ k m=
16Q2
4
2∗ 9 .81∗ π2 0 .3048 ∗ 22 . 794
∑ k m=
16 0 . 27 2
¿32 . 661
∑ k mv = 32 . 661− 7 . 40
∑ k mv = 25 . 261
• Cálculo de la potencia de la bomba con el fin de duplicar el caudal:
En este caso, los datos del problema son:
Se considera que la válvula del literal anterior se elimina. Con estos datos se
calcula el número de Reynolds, la rugosidad relativa y el factor de fricción de
Darcy:
f ks d Q Re f
(m) (m) (m3/s)
f = 0.0103
l v2
hf = f
d 2g
Se puede despejar el factor de fricción de Darcy para obtener:
2 gdh f
f= 2 (2.7)
l v
Ahora, utilizando la ecuación de Colebrook-White se puede despejar la
rugosidad absoluta de la tubería tal como se muestra en las siguientes
ecuaciones:
1 ks 2 . 51
= − 2 lo g
f 3 . 7d Re f
1 ks 2 . 51
− = lo g 10
2 f 3 . 7d Re f
1
− ks
2 f 2 . 51
10 =
3 .7d Re f
Y, finalmente:
1
−
2 f 2 .51
k s = 3 .7d 10 − (2.8)
Re f
(a)
(b)
(b)
a) Manómetro tipo transductor de presión con data logger, con capacidad de 900 registros.
b) Mediciones de presión en una tubería de 300 mm . Medidas de 24 horas, durante 21 días, de
lunes a viernes cada 5 minutos.
Ejemplo 7
Skm = 16.4 H = 32 m
= 1.17 x 10-6 m2/s Q = 3.72 m3/s
d = 1200 mm l = 2800 m
Siguiendo el Diagrama de Flujo No. 6 se calcula en primer lugar el área, la
velocidad, las pérdidas menores y el número de Reynolds:
π
A= d 2 = π× 1 . 2 2 m 2 = 1. 13 m 2
4
3
Q 3 .72 m / s
v= = 2
= 3 . 29 m /s
A 1 .13 m
2 2
v 3 . 29
h m = ∑ k m = 16 . 4× m = 9 . 05 m
2g 2× 9 . 81
vd 3 .29 × 1 . 2 6
Re= = = 3 . 374359 x 10
v 1 .17 × 10 − 6
h f = H− h m = 32 m − 9 . 05 m = 22 . 95 m
Ahora se calcula el factor de fricción de Darcy utilizando la ecuación 2.7:
f = 0 . 0178
Finalmente, utilizando la ecuación 2.8 se calcula la rugosidad absoluta de
la tubería:
1
−
2 f 2 . 51
k s = 3 .7d 10 −
Re f
1
−
2 0 .0178 2 . 51
k s = 3 .7× 1 . 2 10 −
3374359 × 0 . 0178
k s = 0 .78 mm