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, km, QD, H d
, 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.
l v2
hf =f (1.38)
d 2g
1 ks 2.51
= 2log10 +
f 3.7d Re f (1.69)
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 = + z 2 + + 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 = 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 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:
2
l v 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
= 2log10 + (1.69)
f 3.7d Re f
l v K 2.51v l
= 2log10 s
+
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= log10 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= Alog10 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, km, , , 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.
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
3
m
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
= 2log10 +
f 3.7d Re f (1.69)
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
ks
[
f i1= − 2 lo g10
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=200000
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 á lc u lo d e l fa c to r f
M é to d o d e N e w to n -R a p h s o n
0 .0 6
0 .0 5
0 .0 4
f, g (f)
0 .0 3
0 .0 2
0 .0 1
0
0 .0 1 0 .0 2 0 .0 3 0 .0 4 0 .0 5 0 .0 6
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
ks
[
F x i = −2 log10
2. 51
3 .7d Re f ] i=i+1
2. 51 f = 1/x2i+1
F ' x i = [ −2
Ln 10 ] [ ks
Re
2 .51 x i
3 . 7d R e
] xi-+1 = xi
FIN
NO
F x i −x i ?
x i1 =x i − │Xi+1 – Xi │≤ Ɛ
F ' x i −1
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, km, , , , 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
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 – z 2
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).
Siguiendo la metodología de Darcy-Weisbach, la segunda ecuació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,200 0,150 1,429 0,018 0,0252 no 0,343
2,200 0,200 1,717 0,031 0,0539 no 0,496
2,200 0,250 1,977 0,049 0,0971 no 0,658
2,200 0,300 2,217 0,071 0,1567 si 0,827*
1,373 0,300 1,741 0,071 0,1231 si 0,510
1,690 0,300 1,937 0,071 0,1369 si 0,631
1,569 0,300 1,865 0,071 0,1318 si 0,585
1,615 0,300 1,893 0,071 0,1338 si 0,603
1,597 0,300 1,882 0,071 0,1330 si 0,596
1,604 0,300 1,886 0,071 0,1333 si 0,598
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 (km = 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 km = 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
v 2pH
=
2g k m
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.
El procedimiento se esquematiza en el Diagrama de Flujo No.5 el cual es
más general que el Diagrama de Flujo No.4 ya que también sirve para el
caso de tuberías con pérdidas menores bajas. Una vez se ha calculado la
primera velocidad de pérdida, en las demás iteraciones esta velocidad se
calcula de acuerdo con la siguiente ecuación:
2g(H h f )
vp =
k m
(2.6)
INICIO
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
? Vp
Vi <
SI
Procedimiento pérdidas
menores altas A
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 – z 2
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
v p =
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 Qd ∑hm vp
(m) (m) (m/s) (m2) (m3/s) (si o 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 = vp
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
= 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
hf =0 . 0115
0 .3 2×9 . 81
hf =10 . 406 m
Luego las pérdidas menores son:
∑ hm=H−hf
∑ hm=33.2m−10.406m
∑ hm=22.794 m
El coeficiente global de pérdidas menores es:
v2
∑ h m=∑ k m 2g
2g ∑ hm
k
∑ m v2
=
2gπ2 d 4 ∑ hm
∑ k m=
16Q 2
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:
2gdhf
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:
ks
1
f
=−2log
2. 51
3. 7d Re f
ks
−
1
2f
=log10
2. 51
3. 7d Re f
1
−
2f
ks2. 51
10 =
3 .7d Re f
Y, finalmente:
k s =3.7d 10 −
2f
−
2.51
Re f (2.8)
(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
En la red matriz del sistema de abastecimiento de agua de una ciudad se
tiene una tubería de concreto con una longitud de 2.8 km, un diámetro de
1200 mm y un coeficiente global de pérdidas menores de 16.4. En una
determinada condición de operación se mide un caudal de 3.72 m3/s y
una caída en la altura piezométrica de 32 metros a lo largo de toda la
longitud. Calcular la rugosidad absoluta de la tubería. El agua se
encuentra a una temperatura de 14°C.
km = 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.22 m2=1.13m2
4
Q 3.72 m3 / s
v= = 2
=3.29m/s
A 1.13 m
2 2
v 3.29
hm=∑ k m =16.4× m=9.05m
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:
2 gdh f 2×9 .81×1 . 2×22 .95
f= 2 =
l v 2800×3 . 29 2
f =0 . 0178
Finalmente, utilizando la ecuación 2.8 se calcula la rugosidad absoluta de
la tubería:
1
k s =3.7d 10
−
2f
−
2.51
Re f
1
k s =3.7×1.2 10 −
2 0.0178
−
2.51
3374359× 0.0178
k s =0.78 mm