Ecuacion Hazem Williams PDF

Universidad de los Andes
Facultad de Ingeniería
Departamento de Ingeniería Civil y Ambiental
Centro de Investigaciones en Acueductos y

Alcantarillados
CIACUA
Efecto del uso de la ecuación de Hazen-Williams vs la

ecuación de Darcy-Weisbach en tuberías de gran diámetro
(Redes matrices).
Proyecto de grado
Maestría Ingeniería Civil
Rafael Alejandro Flechas Hernández
Bogotá, Diciembre de 2011
I
Universidad de los Andes MIC 201120 07
Centro de Investigaciones en Acueductos y Alcantarillados - CIACUA
Efecto del uso de la ecuación de Hazen-Williams vs la ecuación de
Darcy-Weisbach en tuberías de gran diámetro (Redes matrices).
AGRADECIMIENTOS
El autor agradece la gran colaboración brindada por el ingeniero Juan Guillermo

Saldarriaga Valderrama, asesor de esta tesis, sin la cual hubiera sido posible llevar a cabo
este trabajo. De igual manera, a los ingenieros Daniel Rodriguez y Omar Lopez, jurados
de este trabajo, por sus aportes valiosos para evidenciar y mejorar las falencias del
mismo.
Igualmente, el autor agradece la colaboración prestada por los ingenieros Juan Alberto
Carval Pajaro, Felipe Hernández Cruz y María Alejandra Escovar Bernal, así como del
Centro de Investigaciones en Acueductos y Alcantarilldos (CIACUA) por la ayuda
suministrada para el desarrollo de este trabajo.
II
A Dios, a mi madre, mi padre, mis hermanos y mi abuelita quienes siempre me han
brindado su apoyo incondicional a lo largo de mi vida y son la razón de mi existir.
III
RESUMEN
En el diseño y análisis de redes de distribución de Agua potable (RDAP), las ecuaciones

más empleadas para calcular las pérdidas por fricción son: la ecuación de Darcy-
Weisbach y la ecuación de Hazen-Williams. El uso de esta última ecuación empírica se
ha popularizado alrededor del mundo por su facilidad de manejo, ya que es una ecuación
explícita para el cálculo de la velocidad. Sin embargo, el uso de esta ecuación muchas
veces se hace sin tener en cuenta sus límites de aplicabilidad. Por lo tanto, el uso de esta
no se debería justificar por este hecho sino por los efectos que tenga en el diseño y
análisis de redes matrices, parte fundamental de una RDAP, dado que estas transportan,
con tuberías de gran diámetro, el agua desde las plantas de tratamiento hasta los tanques
de almacenamiento. Por ello, en esta investigación se analizó el efecto del uso de la
ecuación de Darcy-Weisbach vs la Ecuación de Hazen-Williams en el diseño y análisis de
redes matrices. Para ello, se diseñaron y analizaron redes matrices tales como la red
matriz de Medellín y red Hanoi para distintos materiales como el CCP, GRP y Acero. Con
base en este análisis, se encontró que los coeficientes de dicha ecuación empírica,
sugeridos en la normativa colombiana, no son apropiados para los materiales CCP y
Acero dado que estos tienden a sobreestimar las pérdidas por fricción hasta en un 20 y
40%, respectivamente, cuando se comparan con los resultados obtenidos a partir de la
ecuación de Darcy-Weisbach. Adicionalmente, para números de Reynolds superiores a
400.000, independiente del tipo de material que se utilice en el diseño de redes matrices,
la ecuación de Hazen-Williams no satisface los parámetros de diseño de estas redes ya
que no tiene en cuenta la variación que presenta la rugosidad relativa en las redes de
estudio. Por lo anterior, se sugiere utilizar la ecuación de Hazen-Williams solamente si los
factores de fricción obtenidos en el diseño de las redes matrices se localizan en el Flujo
Turbulento Hidráulicamente liso (FTHL), a partir del límite definido por el autor. De lo
contrario, se debe utilizar la ecuación de Darcy-Weisbach en el diseño y análisis de redes
matrices, ya que esta es una ecuación físicamente basada que se puede utilizar para
cualquier tipo de flujo, fluido y rango de diámetros.
IV
TABLA DE CONTENIDO
Capítulo 1. INTRODUCCIÓN......................................................................................................... 1
1.1 Objetivos .............................................................................................................................. 4
Capítulo 2. MARCO TEÓRICO ...................................................................................................... 5
2.1 Ecuaciones de Resistencia Fluida ...................................................................................... 5
2.1.1 Ecuación de Darcy-Weisbach ..................................................................................... 6
2.1.2 Ecuación de Hazen-Williams ..................................................................................... 14
2.1.3 Antecedentes sobre el efecto del uso de las ecuaciones de diseño ........................ 17
2.1.3.1 Vennard (1958) .................................................................................................. 17
2.1.3.2 Diskin (1960) ..................................................................................................... 18
2.1.3.3 Barlow y Markland (1975).................................................................................. 21
2.1.3.4 Jain-Mohan-Khanna (1978) ............................................................................... 24
2.1.3.5 Kamand (1988) .................................................................................................. 26
2.1.3.6 Liou (1998) ........................................................................................................ 29
2.1.3.7 Bombardelli y García (2003).............................................................................. 31
2.2 Redes Matrices .................................................................................................................. 34
2.3 Programa REDES ............................................................................................................. 36
2.3.1 Metodología superficie óptima del Gradiente hidráulico (SOGH) ............................. 37
2.4 Programa EPANET ........................................................................................................... 43
2.5 MATERIALES DE TUBERÍAS ........................................................................................... 44
2.5.1 Rugosidad absoluta (ks) ............................................................................................ 44
2.5.2 Coeficientes de Hazen-Williams (CHw) .................................................................... 45
Capítulo 3. METODOLOGÍA PARA DETERMINAR EL EFECTO DEL USO DE LA ECUACIÓN
DE HAZEN-WILLIAMS ...................................................................................................................... 46
3.1 Selección Redes de Estudio. ............................................................................................ 47
3.1.1 Red matriz Medellín................................................................................................... 47
3.1.2 Red Hanói .................................................................................................................. 48
3.1.3 Red matriz analizada por Bombardelli y García (2003) ............................................ 48
3.2 Selección de materiales .................................................................................................... 48
3.3 Selección Rugosidad de materiales .................................................................................. 49
3.4 Procedimiento cálculo hidráulica de las redes .................................................................. 49
V
3.5 Cálculo de la variación del coeficiente de Hazen-Williams y su efecto en las pérdidas por
fricción y presiones. ....................................................................................................................... 53
3.6 Corrección Coeficientes de Hazen-Williams ..................................................................... 55
3.7 Diseño redes matrices ....................................................................................................... 57
3.7.1 Procedimiento Diseño ............................................................................................... 57
3.8 Validación diseños............................................................................................................. 61
Capítulo 4. RESULTADOS .......................................................................................................... 62
4.1 Efecto del uso de la ecuación de Hazen-Williams en las redes de estudio ...................... 62
4.1.2 Red Hanoi .................................................................................................................. 69
4.1.3 Red matriz analizada por Bombardelli y García (2003) ............................................ 75
4.2 Corrección coeficientes de Hazen-Williams ...................................................................... 80
4.2.1 Programa Corrección Coeficientes de Hazen-Williams ............................................ 81
4.2.2 Coeficientes de Hazen-Williams Corregidos ............................................................. 84
4.3 Diseño redes matrices ....................................................................................................... 89
4.3.1.1 CCP ................................................................................................................... 91
4.3.1.2 GRP ................................................................................................................... 92
4.3.1.3 Acero ................................................................................................................. 93
4.3.2 Red Hanoi .................................................................................................................. 95
4.3.2.1 CCP ................................................................................................................... 95
4.3.2.2 GRP ................................................................................................................... 96
4.3.2.3 Acero ................................................................................................................. 98
4.4 Validación diseño redes matrices ...................................................................................... 99
4.4.1 Red matriz Medellín................................................................................................. 100
4.4.1.1 CCP ................................................................................................................. 100
4.4.1.2 GRP ................................................................................................................. 103
4.4.1.3 Acero ............................................................................................................... 105
4.4.2 Red Hanoi ................................................................................................................ 108
4.4.2.1 CCP ................................................................................................................. 108
4.4.2.2 GRP ................................................................................................................. 111
4.4.2.3 Acero ............................................................................................................... 113
VI
Capítulo 5. ANÁLISIS DE RESULTADOS ................................................................................. 116

5.1 Efecto del uso de la ecuación de Hazen-Williams en las Redes de Estudio .................. 116
5.1.2 Red Hanoi ................................................................................................................ 118
5.1.3 Red matriz analizada por Bombardelli y García (2003) .......................................... 120
5.2 Corrección Coeficientes de Hazen-Williams ................................................................... 120
5.3 Validación Diseños redes matrices ................................................................................. 122
5.3.2 Red Hanoi ................................................................................................................ 129
Capítulo 6. CONCLUSIONES .................................................................................................... 136
BIBLIOGRAFÍA................................................................................................................................ 138
ANEXOS .......................................................................................................................................... 140
VII
ÍNDICE DE FIGURAS
Figura 2.1-1 Diagrama de Moody (Moody, 1940) ............................................................. 13

Figura 2.1-2 Ecuación de Hazen-Williams gráficada en el Diagrama de Moody (Diskin,
1960). .............................................................................................................................. 20
Figura 2.1-3 Nomograma que relaciona el coeficiente de la ecuación de Hazen-Williams y
el factor de fricción para diferentes velocidades y diámetros (Kamand, 1988). ................ 26
Figura 2.1-4 Errores en el cálculo de las pérdidas por fricción entre la ecuación de Hazen-
Williams y la Ecuación de Darcy-Weisbach para Hierro Fundido (Kamand, 1988). .......... 27
Figura 2.1-5 Variación del Coeficiente de la Ecuación de Hazen-Williams como función del
número de Reynolds, rugosidad relativa para agua fluyendo a 15.56ºC, y una rugosidad
absoluta de 0.0003m. ...................................................................................................... 30
Figura 2.1-6 Coeficientes de Hazen-Williams obtenidos de las firmas consultoras en la
red matriz analizada por Bombardelli y García (2003). .................................................... 33
Figura 2.5-1. Diagrama de Flujo Metodología .................................................................. 46
Figura 3.4-1 Selección ecuación de resistencia fluida en el programa EPANET. ............. 50
Figura 3.4-2 Interfaz gráfica del programa EPANET para calcular la hidráulica de la red. 50
Figura 3.4-3 Interfaz gráfica del programa REDES para calcular la hidráulica de la red.
Resultados hidráulica....................................................................................................... 51
Figura 3.4-4 Interfaz programa REDES para obtener factores hidráulicos de la red. ....... 52
Figura 3.4-5 Interfaz programa EPANET para seleccionar resultados de los nudos o
tuberías............................................................................................................................ 52
Figura 3.4-6 Interfaz programa EPANET para seleccionar los resultados de la hidráulica
de la red........................................................................................................................... 52
Figura 3.4-7. Interfaz gráfica programa de EPANET despliegue de resultados. ............... 53
Figura 3.5-1 Diagrama de Flujo para analizar la variación del coeficiente de la ecuación de
Hazen-Williams. ............................................................................................................... 54
Figura 3.6-1 Diagrama de Flujo corrección coeficientes de Hazen-Williams. ................... 56
Figura 3.7-1 Interfaz programa REDES para definir lista de diámetros, curva de costos y
ecuación de fricción en el diseño de RDAP. .................................................................... 58
Figura 3.7-2 Interfaz programa REDES para definir rugosidad absoluta o coeficiente de
Hazen-Williams en las tuberías. ....................................................................................... 59
Figura 3.7-3 Interfaz programa REDES para realizar el diseño de RDAP. ....................... 59
Figura 3.7-4 Interfaz programa Redes para obtener los factores hidráulicos de la red. .... 60
Figura 4.2-1 Interfaz programa para ingresar datos de entrada. ...................................... 82
Figura 4.2-2 Interfaz programa para visualizar resultados del coeficiente de Hazen-
Williams corregido. .......................................................................................................... 83
Figura 4.2-3 Mensaje de Advertencia del programa cuando el diámetro mínimo no cumple
con el error admisible en las pérdidas por fricción. .......................................................... 83
Figura 4.3-1 Sector Red matriz de Medellín delimitado para el Embalse Ayurá. ............. 90
VIII
ÍNDICE DE TABLAS
Tabla 2.1-1 Límites de aplicabilidad de la ecuación de Hazen-Williams (Diskin, 1960). ... 19
Tabla 2.1-2 El efecto del aumento del número de Reynolds en el valor del Coeficiente de
Hazen-Williams, para un diámetro de 1 ft (Barlow y Markland, 1975). .......................... 22
Tabla 2.1-3 Ks/d calculados por Barlow y Markland (1975) en la ecuación de Colebrook-
White, para un diámetro de 1 ft. ...................................................................................... 23
Tabla 2.1-4 Porcentaje de errores en el cálculo de la velocidad y pérdidas por fricción por
el uso de la ecuación de Hazen-Williams (Jain et al, 1978).............................................. 25
Tabla 2.1-5 Porcentaje de error en la estimación del factor de fricción de la ecuación de
Darcy-Weisbach por el uso de la ecuación de Hazen-Williams (Jain et al, 1978)............. 25
Tabla 2.1-6 Valores de Rugosidad Absoluta y coeficientes de la Ecuación de Hazen-
Williams para el PVC y Hierro Fundido (Kamand, 1988). ................................................. 27
Tabla 2.1-7 Valores del Coeficiente de Hazen-Williams promedio para PVC y Hierro
Fundido para un rango de velocidades de 0.3-1.5 m/s y diámetros de 25 a 1050 mm
(Kamand, 1988). .............................................................................................................. 28
Tabla 2.1-8 Rugosidades absolutas obtenidas por Bombardelli y García (2003) para la
tubería principal de la red matriz de estudio. .................................................................... 32
Tabla 2.5-1 Rugosidades absolutas contenidas en las normas de EPM (2009). .............. 44
Tabla 2.5-2 Coeficientes de Hazen-Williams contenidos en las normas de EPM (2009). . 45
Tabla 3.3-1. Rugosidad de los materiales seleccionados para el estudio......................... 49
Tabla 4.2-1 Coeficientes de Hazen-Williams corregidos para los materiales CCP, GRP y
Acero. .............................................................................................................................. 84
Tabla 4.3-1 Lista de diámetro utilizados en el diseño de las redes matrices. .................. 89
IX
ÍNDICE DE GRÁFICAS
Gráfica 4.1-1 Variación del Coeficiente de Hazen-Williams en la red matriz de Medellín
para el CCP. .................................................................................................................... 63
para el GRP. .................................................................................................................... 64
para el Acero. .................................................................................................................. 64
Gráfica 4.1-4 Errores en las pérdidas por fricción en la red matriz de Medellín, utilizando la
ecuación de Darcy-Weisbach y Hazen-Williams para el CCP. ......................................... 66
ecuación de Darcy-Weisbach y Hazen-Williams para el GRP. ......................................... 66
ecuación de Darcy-Weisbach y Hazen-Williams para el Acero. ....................................... 67
Gráfica 4.1-7 Errores en las presiones en los nudos a partir de los resultados obtenidos en
la red matriz de Medellín por medio de la ecuación de Darcy-Weisbach y Hazen-Williams
para el CCP. .................................................................................................................... 68
para el GRP. .................................................................................................................... 68
para el Acero. .................................................................................................................. 69
Gráfica 4.1-10 Variación del Coeficiente de Hazen-Williams en la red Hanoi para el CCP.
........................................................................................................................................ 70
Gráfica 4.1-11 Variación del Coeficiente de Hazen-Williams en la red Hanoi para el GRP.
........................................................................................................................................ 70
Gráfica 4.1-12 Variación del Coeficiente de Hazen-Williams en la red Hanoi para el Acero.
........................................................................................................................................ 71
Gráfica 4.1-13 Errores en las pérdidas por fricción en la red Hanoi, utilizando la ecuación
de Darcy-Weisbach y Hazen-Williams para el CCP. ........................................................ 72
de Darcy-Weisbach y Hazen-Williams para el GRP. ........................................................ 72
de Darcy-Weisbach y Hazen-Williams para el Acero. ...................................................... 73
Gráfica 4.1-16 Errores en las presiones en los nudos a partir de los resultados obtenidos
en la red Hanoi por medio de la ecuación de Darcy-Weisbach y Hazen-Williams para el
CCP. ................................................................................................................................ 74
GRP................................................................................................................................. 74
X
Acero. .............................................................................................................................. 75
Gráfica 4.1-19 Variación del Coeficiente de Hazen-Williams en el primer tramo de la
tubería principal de la red matriz analizada por Bombardelli y García (2003). .................. 76
Gráfica 4.1-20 Error en las pérdidas por fricción en el primer tramo de la tubería principal
de la red matriz analizada por Bombardelli y García (2003). ............................................ 77
Gráfica 4.1-21 Variación del Coeficiente de Hazen-Williams en el segundo tramo de la
Gráfica 4.1-22 Error pérdidas por fricción en el segundo tramo de la tubería principal de la
red matriz analizada por Bombardelli y García (2003). .................................................... 78
Gráfica 4.1-23 Variación del Coeficiente de Hazen-Williams en el tercer tramo de la
Gráfica 4.1-24 Error en las pérdidas por fricción en el tercer tramo de la tubería principal
de la red matriz analizada por Bombardelli y García (2003). ............................................ 80
Gráfica 4.2-1 Corrección Coeficiente de Hazen-Williams para el CCP. ............................ 86
Gráfica 4.2-2 Diagrama de Moody para el Coeficiente de Hazen-Williams corregido del
CCP. ................................................................................................................................ 86
Gráfica 4.2-3 Corrección Coeficiente de Hazen-Williams para el GRP. ........................... 87
GRP................................................................................................................................. 87
Gráfica 4.2-5 Corrección Coeficiente de Hazen-Williams para el Acero. .......................... 88
Acero. .............................................................................................................................. 88
Gráfica 4.3-1 Diámetros obtenidos en el diseño de la red matriz de Medellín para el
material CCP. .................................................................................................................. 91
Gráfica 4.3-2 Diagrama de dispersión presiones obtenidas en el diseño de la red matriz de
Medellín para el material CCP. ........................................................................................ 92
material GRP. .................................................................................................................. 93
Medellín para el material GRP. ........................................................................................ 93
material Acero.................................................................................................................. 94
Medellín para el material Acero. ...................................................................................... 94
Gráfica 4.3-7 Diámetros obtenidos en el diseño de la red Hanoi para el material CCP. ... 96
Gráfica 4.3-8 Diagrama de dispersión presiones obtenidas en el diseño de la red Hanoi
para el material CCP........................................................................................................ 96
Gráfica 4.3-9 Diámetros obtenidos en el diseño de la red Hanoi para el material GRP. ... 97
XI

para el material GRP. ...................................................................................................... 97
Gráfica 4.3-11 Diámetros obtenidos en el diseño de la red Hanoi para el material Acero. 98
para el material Acero. ..................................................................................................... 99
Gráfica 4.4-1 Diagrama de Dispersión de presiones diseño Red matriz de Medellín para
el material CCP, específicamente para los coeficientes de Hazen-Williams iguales a 144 y
120. ............................................................................................................................... 101
Gráfica 4.4-2 Errores pérdidas por fricción en el diseño de la Red matriz de Medellín, para
el material CCP (CHw: 144). .......................................................................................... 102
el material CCP (CHw: 120). .......................................................................................... 102
Gráfica 4.4-4 Diagrama de Dispersión de presiones diseño Red matriz de Medellín para el
material GRP, específicamente para los coeficientes de Hazen-Williams iguales a 147 y
150. ............................................................................................................................... 103
el material GRP (CHw: 147)........................................................................................... 104
el material GRP (CHw: 150)........................................................................................... 105
Gráfica 4.4-7 Diagrama de Dispersión de presiones diseño Red matriz de Medellín para el
material Acero, específicamente para los coeficientes de Hazen-Williams iguales a 139 y
120. ............................................................................................................................... 106
el Acero (CHw: 139)....................................................................................................... 107
el material Acero (CHw: 120). ........................................................................................ 107
Gráfica 4.4-10 Diagrama de Dispersión de presiones diseño Red Hanoi para el material
CCP, específicamente para los coeficientes de Hazen-Williams iguales a 144 y 120. ... 109
Gráfica 4.4-11 Errores pérdidas por fricción en el diseño de la Hanoi, para el material CCP
(CHw: 144). ................................................................................................................... 110
Gráfica 4.4-12 Errores pérdidas por fricción en el diseño de la Red Hanoi, para el material
CCP (CHw: 120). ........................................................................................................... 110
Gráfica 4.4-13 Diagrama de Dispersión de presiones diseño Red Hanoi para el material
GRP, específicamente para los coeficientes de Hazen-Williams iguales a 147 y 150. ... 111
GRP (CHw: 147). ........................................................................................................... 112
Gráfica 4.4-15 Errores pérdidas por fricción en el diseño de la Red matriz de Medellín,
para el material GRP (CHw: 150)................................................................................... 113
Gráfica 4.4-16 Diagrama de Dispersión de presiones diseño Red Hanoi para el Acero,
específicamente para los coeficientes de Hazen-Williams iguales a 139 y 120.............. 114
XII
Gráfica 4.4-17 Errores pérdidas por fricción en el diseño de la Red Hanoi, para el Acero
(CHw: 139). ................................................................................................................... 115
Gráfica 4.4-18 Errores pérdidas por fricción en el diseño de la Red Hanoi, para el Acero
(CHw: 120). ................................................................................................................... 115
Gráfica 5.3-1 Factores de Fricción obtenidos para el diseño de la red matriz de Medellín
para el CCP, específicamente para coeficientes de Hazen-Williams igual a 144 y 120.. 123
para el GRP, específicamente para coeficientes de Hazen-Williams iguales a 147 y 150.
...................................................................................................................................... 126
para el Acero, específicamente para coeficientes de Hazen-Williams iguales a 139 y 120.
...................................................................................................................................... 128
Gráfica 5.3-4 Factores de Fricción obtenidos para el diseño de la red Hanoi para el CCP,
específicamente para coeficientes de Hazen-Williams iguales a 144 y 120. .................. 130
Gráfica 5.3-5 Factores de Fricción obtenidos para el diseño de la red Hanoi para el GRP,
Gráfica 5.3-6 Factores de Fricción obtenidos para el diseño de la red Hanoi para el Acero,
XIII
GLOSARIO DE TÉRMINOS
- Análisis de redes de distribución de agua potable: Procedimiento que consiste en

determinar la magnitud de las presiones en los nudos y caudales que pasan por cada una
de las tuberías de una red específica, resolviendo las ecuaciones de conservación de
masa y energía.
- Coeficiente de Hazen-Williams (CHw): Coeficiente que describe la rugosidad de las

tuberías en la ecuación de Hazen-Williams, el cual depende de la rugosidad relativa, del
flujo y de las propiedades del fluido.
- Densidad (ρ): Característica del fluido que describe la relación entre la masa y el
volumen del mismo.
- Diagrama de Moody: Es un diagrama que permite encontrar de forma gráfica el factor

de fricción de la ecuación de Darcy-Weisbach.
- Diseño de Redes de Distribución de Agua potable: Procedimiento que consiste en

determinar la magnitud de los diámetros de una red especifica, teniendo en cuenta los
caudales de consumo y la presión mínima exigida por una norma local.
- Ecuación de Colebrook-White: Ecuación logarítmica e implícita que permite calcular el

factor de fricción de la Ecuación de Darcy-Weisbach de forma analítica.
- Ecuación de Darcy-Weisbach: Ecuación físicamente basada que permite determinar

las pérdidas por fricción que ocurren en una tubería debido a la fricción entre el fluido y las
paredes de la misma.
- Ecuación de Hazen-Williams: Ecuación de origen empírico que describe el cálculo de

las pérdidas por fricción en tuberías que transportan agua en condiciones normales.
- EPANET: Es un programa desarrollado por la Agencia de Protección Ambiental de los

Estados Unidos, para el estudio y análisis del comportamiento de redes de distribución de
agua potable.
XIV
- Factor de fricción (f): Es un parámetro característico de la Ecuación de Darcy-

Weisbach, el cual depende de la rugosidad relativa y del número de Reynolds.
- Fluido: Sustancia que no puede resistir esfuerzo cortante.
- Flujo: Es el movimiento de las partículas de un fluido con respecto a un sistema inercial.
- Flujo laminar: Es el tipo de flujo donde se presentan bajas velocidades y viscosidades

altas, lo cual impide el intercambio de moléculas y mezcla entre las capas de flujo.
- Flujo turbulento: Es el tipo de flujo que se caracteriza por presentar altas velocidades,
lo que permite el intercambio de moléculas y la mezcla entre las capas de flujo, anulando
los efectos viscosos.
- Flujo turbulento hidráulicamente liso (FTHL): Flujo turbulento en una tubería donde la
rugosidad absoluta es inferior al 30% del espesor de la subcapa laminar viscosa. En este
tipo de flujo, el factor de fricción depende exclusivamente de las propiedades del flujo y de
las características del fluido.
- Flujo turbulento hidráulicamente rugoso (FTHR): Flujo turbulento en una tubería

donde la rugosidad absoluta es 6.1 veces mayor que el espesor de la subcapa laminar
viscosa. En este tipo de flujo, el factor de fricción es función únicamente de la rugosidad
relativa de la tubería.
- Flujo turbulento transicional: Flujo turbulento acotado entre los flujos que se pueden
clasificar como máximo en flujo turbulento hidráulicamente liso y como mínimo en flujo
turbulento hidráulicamente rugoso.
- Método del Gradiente: Es un método de cálculo de redes de distribución de agua

potable que permite corregir de forma simultánea las presiones en los nudos y caudales
en las tuberías, utilizando matrices dispersas que garantizan la convergencia y la
eficiencia del mismo.
- Número de Reynolds (Re): Parámetro adimensional que describe la relación entre las
fuerzas gravitacionales y fuerzas viscosas. Se utiliza para clasificar los flujos, ya sea en el
régimen de flujo laminar o flujo turbulento.
XV
- Pendiente de la línea de fricción (S): Es la relación entre las pérdidas por fricción y la
longitud de una tubería.
- Pérdidas por fricción (hf): Son las pérdidas de energía que se generan en una tubería
debido a la fricción que se genera entre el fluido y las paredes de la misma.
- Radio hidráulico (R): Es la relación entre el área mojada de una tubería y el perímetro
mojado de la misma.
- REDES: Es un programa desarrollado por el Centro de Investigaciones de Acueductos y

Alcantarillados (CIACUA) de la Universidad de los Andes, el cual permite analizar y
diseñar redes de distribución de agua potable.
- Redes de distribución de agua potable (RDAP): Las redes de distribución de agua

potable (RDAP) son una serie de obras hidráulicas que permiten el suministro de agua
potable a una población específica.
- Redes matrices: Es un componente de las redes de distribución de agua potable, las

cuales se caracterizan por presentar los siguientes componentes: fuentes de
abastecimiento, obras de captación, obras de aducción y conducción para el transporte de
las aguas sin tratar hasta las plantas de tratamiento, plantas de tratamiento, tanques de
almacenamiento y tuberías de gran diámetro.
- Rugosidad relativa (ks/d): Relación entre la rugosidad absoluta de las tuberías y el

diámetro de las mismas.
-Subcapa laminar viscosa (δ`): Es una parte de la capa límite turbulenta localizada
cerca a las paredes de una tubería, donde se atenúa las vibraciones y el intercambio
molecular característico del flujo turbulento. Por tal razón, el flujo en esta zona siempre es
laminar.
-Viscosidad cinemática (ν): Es una característica del fluido, que relaciona la viscosidad
dinámica y la densidad del mismo.
XVI
Capítulo 1. INTRODUCCIÓN
Las redes de distribución de agua potable (RDAP) son una serie de obras hidráulicas que
permiten el suministro de agua potable a una población específica. Estas redes están
compuestas por dos partes, las cuales son: la red matriz, parte encargada de conducir el
agua desde la planta de tratamiento hasta los tanques de almacenamiento, y la segunda
parte son las redes de distribución secundaria o terciaria, cuya función es suministrar el
agua potable desde los tanques de almacenamiento hasta los usuarios de la red.
Teniendo en cuenta la importancia que tienen las redes de distribución de agua potable
(RDAP) como pieza fundamental del desarrollo de una sociedad, a través de la historia el
análisis y diseño de este tipo de sistemas se ha convertido en el principal objetivo de los
ingenieros. Para llevar a cabo dicho propósito, es necesario establecer con certeza la
magnitud de las pérdidas por fricción que se generan en las tuberías, dado que estas
determinan la magnitud de los diámetros, la distribución de presiones en los nudos y los
caudales en cada una de estas. Las ecuaciones más empleadas para calcular las
pérdidas por fricción son: la ecuación de Darcy-Weisbach, ecuación físicamente basada, y
la ecuación de Hazen-Williams, ecuación de origen experimental que es válida solamente
en los rangos investigados para su formulación.
En relación con la ecuación de Darcy-Weisbach, se puede decir que esta se puede aplicar
para cualquier fluido, rango de velocidades, diámetros y longitudes de tubería, debido a
que es una ecuación físicamente basada. Sin embargo, en la primera mitad del siglo XX,
el cálculo de las pérdidas por fricción por medio de esta ecuación se dificultó, ya que el
factor de fricción de esta era un parámetro difícil de cuantificar analíticamente. Esto se
refleja en la ecuación de Colebrook-White, la cual requiere algún método numérico para
calcular el factor de fricción, debido a que es una ecuación logarítmica no explícita para
esta variable hidráulica.
1
Dicho problema ocasionó la aparición de ecuaciones empíricas más sencillas de utilizar

como la ecuación de Hazen-Williams, la cual presenta la ventaja de ser explícita para la
velocidad y el caudal, pero presenta inconvenientes desde el punto de vista hidráulico,
dado que utiliza un coeficiente de rugosidad que no depende solamente de la rugosidad
absoluta de la tubería, sino de la velocidad y la magnitud del diámetro, así como de las
propiedades del fluido ( Vennard 1958; Diskin 1960; Liou 1998).
Por lo anterior, el uso de la ecuación de Hazen-Williams frente a la ecuación de Darcy-

Weisbach no se debería justificar por su facilidad de manejo, la cual es relativa debido a la
popularización del uso de las computadoras que permite resolver de forma sencilla
ecuaciones implícitas como la ecuación de Colebrook-White, sino por los efectos que
tenga el uso de esta en el diseño y análisis de RDAP desde el punto de vista hidráulico,
teniendo en cuenta el creciente aumento de la demanda de agua en las ciudades, la
continúa operación y las características de las tuberías de estas redes.
Esto podría establecer realmente como debería utilizarse dicha ecuación empírica en la
práctica de la ingeniería, ya que su aplicación inadecuada podría determinar si una RDAP
estaría en capacidad de suministrar agua a una determinada población, dado que se
podría estar perdiendo mayor cantidad de energía que la esperada, lo cual no permitiría
cumplir con las presiones mínimas exigidas en las normativas de cada uno de los países
del mundo.
En este sentido, algunos investigadores como Bombardelli y García (2003) han analizado
las posibles implicaciones que puede tener el uso de la ecuación de Hazen-Williams
posteriormente al diseño de la línea principal de una Red Matriz, pieza fundamental de
una RDAP, en cuanto a errores en el análisis hidráulico cuando se utiliza está por fuera de
sus rangos de aplicabilidad. Sin embargo, en la actualidad no se han realizado
investigaciones para establecer el efecto del uso de esta ecuación empírica en el diseño
de redes matrices en si (Walski, 2003).
2
Por esta razón, en el presente documento se analizará el efecto del uso de la ecuación de
Darcy-Weisbach frente a la ecuación de Hazen-Williams en RDAP, pero específicamente
en redes matrices, teniendo en cuenta que el creciente aumento de la demanda exige
redes que presenten tuberías de gran diámetro.
Por ello, este documento está estructurado de la siguiente manera:
El Capítulo 1 contiene los objetivos de la tesis, explicando lo que se buscó con el

desarrollo del proyecto.
El Capítulo 2 muestra la base teórica necesaria para llevar a cabo este proyecto. En este
se hace una breve reseña histórica del desarrollo de la ecuación de Darcy-Weisbach y de
la ecuación de Hazen-Williams. Igualmente, se resume la revisión bibliográfica realizada
en cuanto a los antecedentes del efecto de uso de estas ecuaciones en tuberías y redes
matrices. Adicionalmente, en este capítulo se hace referencia a los programas REDES y
EPANET, programas utilizados para analizar el efecto de estas ecuaciones. Así mismo,
se ilustran las rugosidades y coeficientes de Hazen-Williams de cada uno de los
materiales de tuberías que más se utilizan en las redes matrices en el medio colombiano.
El Capítulo 3 contiene la metodología que se aplicó para establecer el efecto del uso de
estas ecuaciones de resistencia fluida en el diseño y análisis de RDAP.
El Capítulo 4 presenta los resultados obtenidos en el análisis de estas ecuaciones en las

redes matrices seleccionadas. Posteriormente, en el Capítulo 5 se ilustra el análisis de
estos resultados.
En el Capítulo 6 se muestran algunas conclusiones de esta tesis a tener en cuenta en el

diseño de redes matrices. Por último, se incluyen las referencias bibliográficas
consultadas para el desarrollo de este trabajo.
3
1.1 Objetivos
A continuación se presentan los objetivos que motivaron el desarrollo de este proyecto.
Objetivo general
• Determinar el efecto del uso de ecuaciones empíricas vs ecuaciones físicamente

basadas en el comportamiento hidráulico de tuberías de gran diámetro (Redes
Matrices).
Objetivos específicos
Para este objetivo general, se plantearon los siguientes objetivos específicos para poder
desarrollar esta tesis:
• Establecer el efecto del uso de la ecuación de Hazen-Williams vs la ecuación de

Darcy-Weisbach en el cálculo de las pérdidas por fricción y la distribución de
presiones en una red matriz.
• Determinar si los coeficientes de la ecuación de Hazen-Williams sugeridos en la
normativa colombiana, deben ser modificados o no, con base en los efectos
hidráulicos que tienen en una red matriz.
• Establecer una metodología que permita corregir los valores sugeridos del
coeficiente de Hazen-Williams para un determinado material, que permita
minimizar los errores en el cálculo de las pérdidas por fricción al utilizar esta
ecuación.
• Verificar si esta metodología en el proceso de diseño de redes matrices, garantiza
el cumplimiento de la presión mínima en la normativa colombiana.
4
Capítulo 2. MARCO TEÓRICO

A través de la historia, el cálculo de las pérdidas por fricción en redes de distribución de
agua potable (RDAP) ha sido motivo de investigación, dado que estas pueden determinar
la magnitud de los diámetros, presiones en los nudos y caudales que pasan por las
tuberías de este tipo redes. Teniendo en cuenta lo anterior, a continuación se describen el
desarrollo de las dos ecuaciones que más se utilizan en la práctica de la ingeniería para el
cálculo de estas y los efectos del uso de las mismas que pueden tener en RDAP.
2.1 Ecuaciones de Resistencia Fluida

En el flujo de tuberías, las pérdidas de energía pueden ser cuantificadas a partir de la
ecuación de Bernoulli (Conservación de la Energía), la cual se ilustra a continuación.
𝑝1 𝑉12 𝑝2 𝑉22
+ + 𝑧1 = + + 𝑧2 + 𝐻𝑝 1−2
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔
Ecuación 2.1-1
𝑝1 𝑉12 𝑝2 𝑉22
𝐻𝑝 1−2 = � + + 𝑧1 � − � + + 𝑧2 �
𝜌𝑔 2𝑔 𝜌𝑔 2𝑔
Ecuación 2.1-2
donde:
z= Cota del fondo de la tubería.

p= Presión.
V= Velocidad de flujo.
ρ= Densidad del fluido.
g= Gravedad.
Hp1-2= Pérdidas de energía entre el punto 1 y 2 de la tubería.
Según la Ecuación 2.1-2, se puede ver que las pérdidas de energía, las cuales son
ocasionadas por la fricción con las paredes de las tuberías o accesorios que se colocan a
lo largo de esta, solamente se pueden determinar si se hacen mediciones de caudal y
presión en una red existente, lo cual no tiene ningún sentido en el diseño de tuberías.
5
Por lo anterior, en el diseño y análisis de redes de distribución de agua potable (RDAP),

es indispensable conocer la relación que existe entre las pérdidas de energía como
función del fluido, la velocidad, el diámetro y el tipo de material en cada una de las
tuberías, para poder cumplir y verificar los requerimientos de demanda y presión en una
red específica, sin hacer mediciones. Por ello, a lo largo de la historia se han planteado
diversas ecuaciones para tal fin. Estas se pueden clasificar en dos familias: las teóricas y
las empíricas. La ecuación más reconocida de la primera familia es la ecuación de Darcy
(1803-1871)-Weisbach (1806-1871). Por otro lado, la Ecuación de Hazen-Williams es una
de las ecuaciones empíricas más exitosas (desarrollada por G.S.Williams y A.H. Hazen en
1906), debido a que es explícita en relación con la velocidad y presenta ventajas de
aplicabilidad en el diseño de tuberías con respecto a las ecuaciones teóricas. Teniendo
en cuenta lo anterior, a continuación se presentará una breve reseña histórica sobre el
desarrollo de cada una de estas ecuaciones y el efecto del uso de las mismas.
2.1.1 Ecuación de Darcy-Weisbach
En 1845, Julius Weisbach (1806-1871) ingeniero, de nacionalidad sajona, propuso una

ecuación para predecir las pérdidas por fricción (sin tener en cuenta las pérdidas menores
ocasionadas por accesorios) con base en el fluido, la velocidad, el diámetro y el tipo de
material. Esta ecuación se ilustra a continuación.
𝐿 𝑉2
ℎ𝑓 = 𝑓
𝐷 2𝑔
Ecuación 2.1-3
donde:
hf= Pérdidas por fricción.

L= Longitud de la tubería.
D= Diámetro de la tubería.
g= Aceleración de la gravedad.
f= Factor de fricción.
6
Esta es la forma que actualmente se conoce de la ecuación de Darcy-Weisbach; sin

embargo para ese momento, Weisbach describió el factor de fricción de la siguiente
manera.
𝛽
𝑓=𝛼+
√𝑉
Ecuación 2.1-4
donde α y β son coeficientes de fricción que pueden variar con el diámetro y material de
la tubería. Esta última ecuación fue obtenida a partir de un pequeño conjunto de
experimentos basados en artículos de investigadores contemporáneos como Gaspard
Riche de Prony (1755-1839), ingeniero francés del siglo XIX. A pesar de esto, la
comunidad francesa, en ese entonces, se resistió a cambiar la ecuación de Weisbach por
la ecuación propuesta por Prony (Brown, 2004).
𝐿
ℎ𝑓 = (𝑎𝑉 + 𝑏𝑉 2 )
𝐷
Ecuación 2.1-5
donde a y b eran coeficientes empíricos.
La desventaja de esta ecuación era que los coeficientes empíricos dependían del sistema
de unidades, lo cual no se presentaba en la ecuación propuesta por Weisbach, ya que
esta era dimensionalmente homogénea. Por ello, el factor de fricción era un número
adimensional, consistente con el sistema de unidades.
De otro lado, en 1857, Henry Darcy (1803-1858) propuso una nueva forma de la ecuación
de Prony, utilizando experimentos de flujo en tuberías, los cuales se desarrollaron en
Chaillot cerca a Paris desde el 3 de agosto de 1849 hasta el 27 de octubre de 1850, y se
realizaron bajo las siguientes condiciones:
• Diámetros desde 0.0122 hasta 0.5 m.

• Materiales como vidrio, hierro, plomo, hierro recubierto con bitumen, hierro fundido
y latón.
• Rango de velocidades promedio entre 0.03 hasta 6.01 m/s.
7
• Longitud de tuberías superior a 100 m, sin contar las tuberías de vidrio y de plomo,
las cuales fueron de 44.8 y 50 m respectivamente.
• Utilización de piezómetros, sin tener en cuenta los 4.7 m de la longitud de las
tuberías.
A partir de lo anterior, la ecuación propuesta por Darcy para nuevas tuberías fue la
siguiente
𝐿 𝛽 𝛽´
ℎ𝑓 = ��𝛼 + 2 � 𝑉 + �𝛼´ + � 𝑉 2 �
𝐷 𝐷 𝐷
Ecuación 2.1-6
Donde α, β, α´ y β´ son coeficientes de fricción. Según Darcy, el primer término de esta

ecuación se podía omitir para tuberías de material usado. Por ello, propuso la siguiente
ecuación para estas tuberías.
𝐿 𝛽´
ℎ𝑓 = �𝛼´ + � 𝑉 2
𝐷 𝐷
Ecuación 2.1-7
donde α´ y β´ pueden ser coeficientes de fricción distintos a los valores de tuberías de

material nuevo. Es importante mencionar, que a partir de estas ecuaciones, Darcy mostró
que el factor de fricción era función de la rugosidad y el diámetro, a pesar de que esta
conclusión no se haya descrito en la publicación de estas ecuaciones (Brown, 2004).
Posteriormente, J.T Fanning (1837-1911), en 1877, fue el primero en combinar los

resultados experimentales de Henry Darcy y la ecuación de Julius Weisbach para obtener
mejores estimaciones del factor de fricción. Para ello, publicó tablas con valores de este
parámetro como función del material, diámetro y velocidad en las tuberías. Sin embargo,
Fanning en estas tablas utilizó el radio hidráulico en vez del diámetro en la ecuación de
fricción.
Teniendo en cuenta lo anterior, la ecuación de Darcy-Weisbach fue deducida por Henry

Darcy, ingeniero frances del siglo XIX que llevó a cabo numerosos experimentos en
tuberías con flujo de agua, y Julius Weisbach, ingeniero sajon de la misma época, el cual
8
propuso el uso del factor de fricción. Esta ecuación es considerada como una fórmula
racional debido a que un balance de fuerzas simple y un análisis dimensional establece
que:
ℎ𝑓 𝛼 𝐿 𝐷 −1 𝑉 2 𝑔−1
Ecuación 2.1-8
Por ello, esta ecuación se deduce mediante análisis dimensional, y por lo tanto es una
ecuación basada en la física clásica, haciendo uso de las variables encontradas por
Weisbach (ver Ecuación 2.1-3), que son relevantes para entender el comportamiento de
las pérdidas por fricción en tuberías. De este análisis, se obtiene la ecuación de Darcy-
Weisbach.
𝐿 𝑉2 𝑘𝑠
ℎ𝑓 = � � 𝐹´´ � , 𝑅𝑒�
𝐷 2𝑔 𝑑
Ecuación 2.1-9
Si los argumentos de la función F´´ se denotan como el factor de fricción, se puede ver
𝑘𝑠
que este parámetro depende de la rugosidad relativa � � y el número de Reynolds (Re),
𝑑
que se define de la siguiente manera.
𝑉𝐷
𝑅𝑒 =
𝜈
Ecuación 2.1-10
donde V es la velocidad, D el diámetro y 𝜈 la viscosidad cinemática del fluido, parámetro

que depende de la temperatura. Este número a dimensional permite establecer los tipos
de flujo que se pueden presentar en tuberías de sección circular, los cuales son:
• Flujo laminar Re<2000.

• Flujo transicional 2000<Re<4000.
• Flujo turbulento Re>4000.
9
De todos los parámetros de la ecuación de Darcy-Weisbach, el factor de fricción es el más

difícil de estudiar ya que es una función compleja de la velocidad de flujo, de la densidad y
viscosidad del fluido, del diámetro y de la rugosidad de la tubería, como se mostró
anteriormente. Esta complejidad, se reflejó en los intentos realizados por investigadores
como Ludwig Prandtl (1875-1953), Theodor Von Kármán (1881-1963), Paul Richard
Heinrich Blasius (1883-1970) y Johann Nikuradse (1894-1979) para predecir
analíticamente su magnitud. Sin embargo, en 1939, C.F. Colebrook y C.M. White lograron
encontrar una relación matemática para describir el comportamiento de este en la zona de
transición, dicha relación se describe a continuación:
1 𝑘𝑠 2.51
= −2𝑙𝑜𝑔10 � + �
�𝑓 3.7𝑑 𝑅𝑒�𝑓
Ecuación 2.1-11
donde
f= Factor de fricción.
ks= Rugosidad absoluta de la tubería.
d= Diámetro de la tubería.
Re= Número de Reynolds.
Esta ecuación probó ser válida para todo tipo de flujo en tuberías (flujo laminar, zona de
transición y flujo turbulento). Sin embargo, tiene el problema de no ser explícita para
encontrar el factor de fricción; por ello se requiere del uso de algún método numérico para
encontrar este parámetro. Por esta razón, a principios de la década de 1940, el ingeniero
Lewis Moody, desarrollo un diagrama 1 que reúne la información respecto a los tipos de
flujo que se pueden presentar en tuberías, con base en los experimentos realizados por
J. Nikuradse y C.F. Colebrook, donde es posible encontrar de una forma ágil y práctica el
factor de fricción conociendo el número de Reynolds y la rugosidad relativa.
1
(Moody, 1944)
10
En este diagrama (ver Figura 2.1-1) se puede diferenciar cuatro áreas, las cuales son:
flujo laminar, zona crítica, zona de transición y turbulencia completa. La primera región
está delimitada por números de Reynolds menores a 2000. En esta, el factor de fricción
se encuentra mediante la siguiente ecuación, la cual pudo determinar Weisbach con base
en la ecuación de Hagen (1794-1884)-Poiseuille (1799-1869) para flujo laminar.
64
𝑓=
𝑅𝑒
Ecuación 2.1-12
La segunda área, por su parte, está delimitada por números de Reynolds entre 2000 y
4000. En esta región es difícil de determinar el factor de fricción ya que no es posible
encontrar una expresión matemática debido a la inestabilidad que se presenta en dicha
zona.
Para números de Reynolds superiores a 4000, en el diagrama de Moody se pueden ver

dos regiones, la zona de transición y turbulencia completa. El paso de flujo laminar a
flujo turbulento no se hace de forma súbita sino gradual, por ello la zona que describe
dicho cambio de flujo es la zona de transición. A partir de los resultados obtenidos por
Colebrook y White, los límites de la zona de transición están definidos por el límite
superior del flujo turbulento hidráulicamente liso (FTHL) y el límite inferior del flujo
turbulento hidráulicamente rugoso (FTHR), lo cual se resumen en la siguiente expresión.
0.305𝛿´ < 𝑘𝑠 < 6.1𝛿´

Ecuación 2.1-13
donde:
ks= Rugosidad de la tubería.

δ´= Espesor de la subcapa laminar viscosa.
11
Para el caso del FTHL, se puede utilizar la ecuación de Blasius para calcular el factor de
fricción, la cual es válida para números de Reynolds situados entre 5000 y 100.000, y se
describe a continuación.
0.316
𝑓=
𝑅𝑒 0.25
Ecuación 2.1-14
Por su parte, Prandtl y Von Kármán en sus investigaciones, encontraron las siguientes
ecuaciones para calcular el factor de fricción tanto en el FTHL y FTHR.
1 2.51
= −2𝑙𝑜𝑔10 � � 𝐹𝑇𝐻𝐿
�𝑓 𝑅𝑒�𝑓
Ecuación 2.1-15
1 𝑘𝑠
= −2𝑙𝑜𝑔10 � � 𝐹𝑇𝐻𝑅
�𝑓 3.7𝑑
Ecuación 2.1-16
De igual forma, recientemente Flechas (2010) dedujo las ecuaciones para delimitar el
FTHR y FTHL, con base en la ecuación de Colebrook-White (Ver Ecuación 2.1-11) y los
límites establecidos por estos dos investigadores para la zona de transición (Ver Ecuación
2.1-13).
1 5.21
= −2𝑙𝑜𝑔10 � � 𝐹𝑇𝐻𝐿
Ecuación 2.1-17
1 56.6
= −2𝑙𝑜𝑔10 � � 𝐹𝑇𝐻𝑅
Ecuación 2.1-18
Por último, en la zona de turbulencia completa, el factor de fricción se puede determinar

con base en la ecuación de Colebrook-White (ver Ecuación 2.1-11). Es importante aclarar
que para encontrar el factor de fricción en cualquiera de las zonas del diagrama de Moody
se puede utilizar dicha ecuación, debido a que esta describe perfectamente todos los tipos
de flujos que se pueden presentar en las tuberías.
12
Efecto del uso de la ecuación de Hazen-Williams vs la ecuación de Darcy-Weisbach en tuberías de gran
diámetro (Redes matrices).
Figura 2.1-1 Diagrama de Moody (Moody, 1940)
13
Universidad de los Andes ICIV 201120 07
Efecto del uso de la ecuación de Hazen-Williams vs la ecuación de Darcy-
Weisbach en tuberías de gran diámetro (Redes matrices).
Como se describió anteriormente, el problema matemático que presentaba la ecuación de

Colebrook-White dio origen al desarrollo de métodos gráficos como el diagrama de Moody
para cuantificar el factor de fricción. Por otro lado, este problema también ocasionó la
aparición de ecuaciones empíricas más fáciles de utilizar como la ecuación de Hazen-
Williams y ecuaciones explícitas para encontrar el factor de fricción, tales como la
ecuación de Swamee y Jain (1976), ecuación explícita más exitosa para calcular dicho
parámetro hidráulico debido a su precisión cuando se compara con la ecuación de
Colebrook-White.
0.25
𝑓=
𝑘𝑠 5.74 2
�− log �3.7𝐷 + 0.9 ��
𝑅𝑒
Ecuación 2.1-19
2.1.2 Ecuación de Hazen-Williams
Una de las ecuaciones empíricas más exitosas, que no tuvieron en cuenta el análisis de la
ecuación de Darcy-Weisbach, fue la ecuación de Hazen-Williams (Hazen & Williams,
1905) propuesta por Gardner Williams (Profesor de ingeniera civil, sanitaria e hidráulica
de la Universidad de Michigan) y Allen Hazen (Ingeniero civil y sanitario) en el año de
1903. La forma original de esta ecuación se ilustra a continuación.
𝑉 = 𝑏𝐶𝐻𝑊 𝑅0.63 𝑆 0.54

Ecuación 2.1-20
donde:
CHW= Coeficiente de la ecuación de Hazen-Williams, el cual depende del material y edad
de la tubería.
R= Radio hidráulico.
S= Pendiente de la línea de fricción o la relación entre las pérdidas por fricción y la
longitud de la tubería.
b= Factor multiplicador para igualar la velocidad media calculada con base en la ecuación
de Chézy 2, la cual se muestra a continuación:
2
Ecuación propuesta por Antoine Chézy (1718-1798) en el año 1770. Esta ecuación, en el año de
1903, era la más utilizada para el cálculo de la velocidad de flujo en tuberías y canales.
14
𝑉 = 𝐶√𝑅𝑆
Ecuación 2.1-21
Teniendo en cuenta que el propósito de la ecuación de Hazen-Williams, en su momento,

era modificar la ecuación de Chézy (ver Ecuación 2.1-21), debido a que el factor C de
esta ecuación era difícil de calcular dado que dependía de la rugosidad de la tubería, el
radio hidráulico y la pendiente de fricción, la magnitud del factor b de la ecuación de
Hazen-Williams se obtuvo igualando esta con la ecuación de Chézy para lograr la misma
velocidad para un coeficiente dado, suponiendo un radio hidráulico de 1ft (0.3 m) y una
pendiente de la línea de fricción promedio de 1/1000, tal como se ilustra a continuación:
𝑏𝐶𝐻𝑊 𝑅0.63 𝑆 0.54 = 𝐶√𝑅𝑆
𝐶𝐻𝑤 = 𝐶
𝑏 = 𝑅 −0.13 𝑆 −0.04
Ecuación 2.1-22
Por lo anterior, el factor b puede adoptar el valor de 0.849 para el Sistema Internacional
de unidades y 1.318 si el sistema de unidades utilizado en el cálculo de flujo tuberías es
el Sistema Inglés. Por ello, la forma original de esta ecuación en el sistema internacional
de unidades es:
𝑉 = 0.849𝐶𝐻𝑊 𝑅 0.63 𝑆 0.54
Ecuación 2.1-23
Adicionalmente, para solucionar el problema de la variabilidad del coeficiente C de la

ecuación de Chézy, Hazen y Williams (1905) vieron la necesidad de modificar los
exponentes de esta ecuación, de tal forma de que éstos representaran las condiciones
promedio en cuanto al radio hidráulico y la pendiente de la línea de fricción, y de esa
manera determinar los valores de dicho coeficiente que dependiera exclusivamente de la
rugosidad. Para ello, estos investigadores llevaron a cabo un análisis estadístico de la
información de diversos experimentos realizados por investigadores como Darcy, Iben y
Fenkell entre otros, sobre flujo de agua en tuberías que se relacionaban con las
condiciones que más se presentaban en la práctica de la ingeniería. Con base en ello, los
15
exponentes obtenidos por estos investigadores fueron 0.63 y 0.54 para el radio hidráulico
y la pendiente de la línea de fricción respectivamente (ver Ecuación 2.1-23).
Teniendo en cuenta lo anterior, en el año de 1905, Hazen y Williams decidieron publicar el

libro titulado Hydraulic Tables, para mostrar tablas que permitieran encontrar el valor de
CHw para cada uno de los materiales de las tuberías (hierro, acero, plomo, bronce, vidrio y
estaño) en "condiciones normales", para luego reemplazar en la Ecuación 2.1-23 y
encontrar la velocidad de flujo o pérdidas de fricción en las tuberías, como se muestra a
continuación.
6.8241𝐿 1.851
ℎ𝑓 = 1.851 1.167 𝑉
𝐶𝐻𝑊 𝐷
Ecuación 2.1-24
donde:
hf= Pérdidas por fricción.
L= Longitud de las tuberías.
D=Diámetro de la tubería.
CHW= Coeficiente de la ecuación de Hazen-Williams, el cual depende del material y edad
de la tubería.
A partir de las tablas contenidas en el libro Hydraulic Tables 3, se puede ver que el valor de
este coeficiente puede variar entre 95 y 140 dependiendo del material. Por ejemplo, para
tuberías como el acero el CHw es igual a 95, pero para tuberías nuevas de vidrio puede
llegar a valores de hasta 140. Adicionalmente, es importante destacar que en dicho libro,
los autores establecen que esta ecuación solamente se puede utilizar cuando el fluido es
agua a temperaturas normales y que el coeficiente de esta ecuación empírica puede
verse disminuido por el deterioro de las tuberías o por diámetros pequeños ( inferiores a 4
pulgadas).
3
(Williams & Hazen, 1905).
16
Por lo anterior, la ecuación de Hazen-Williams presenta los siguientes rangos de validez,

ya que fue obtenida para "condiciones normales" que se presentaban en la práctica de la
ingeniería:
• El fluido debe ser agua a temperaturas normales (20oC).

• El diámetro debe superior a 100 mm o 4 pulgadas.
• La velocidad de flujo debe ser inferior a 3 m/s.
• Pendiente de la línea de fricción promedio de 1/1000.
• Radio hidráulico promedio de 0.3 m.
2.1.3 Antecedentes sobre el efecto del uso de las ecuaciones de

diseño
Teniendo en cuenta que la ecuación de Hazen-Williams solamente se puede aplicar para

un rango limitado de diámetros y velocidades, a lo largo del Siglo XX se han realizado
diversos estudios para establecer cuáles son los efectos de utilizar esta ecuación en el
cálculo de flujo de tuberías en comparación con la ecuación de Darcy-Weisbach. Los más
relevantes se ilustran a continuación en orden cronológico.
2.1.3.1 Vennard (1958)
El primero estudio conocido en relación con el efecto del uso de la ecuación de Hazen-
Williams fue el realizado por Vennard (1958), para establecer si esta ecuación empírica se
podía utilizar para cualquier tipo de flujo y si el coeficiente de esta ecuación dependía
únicamente de la rugosidad o podía verse afectado por el número de Reynolds y la
rugosidad relativa de la tubería. Para ello, este investigador reescribió la ecuación de
Hazen-Williams en la forma de la ecuación de Darcy-Weisbach, con el fin de relacionar el
coeficiente de esta ecuación empírica con el factor de fricción, como se ilustra a
continuación.
1090
𝑓= 1.85 0.015
𝐶𝐻𝑊 𝐷 𝑅𝑒 0.015
Ecuación 2.1-25
17
Aproximando el diámetro a la unidad y suponiendo una viscosidad cinemática nominal del

agua correspondiente a 15oC, en el Sistema Inglés de unidades, Vennard (1958) con base
en la Ecuación 2.1-24, graficó los factores de fricción en el diagrama de Moody para
distintos valores del coeficiente de Hazen-Williams. A partir de ello, este investigador
encontró que el coeficiente no dependía únicamente de la rugosidad de la tubería, sino
también de la rugosidad relativa. Así mismo, pudo establecer que la ecuación de Hazen-
Williams no se podía utilizar en el FTHR, teniendo en cuenta que, para un valor dado del
coeficiente de esta ecuación, el factor de fricción variaba con una pendiente constante de
0.15 a medida que aumentaba el número de Reynolds en dicho diagrama.
2.1.3.2 Diskin (1960)
Posteriormente, a diferencia de Vennard (1958), Diskin (1960) estableció claramente los

límites de aplicabilidad de la ecuación de Hazen-Williams para los distintos tipos flujo.
Para ello, este investigador, primero, reorganizó la ecuación de Hazen-Williams en la
forma de la ecuación de Darcy-Weisbach de una manera similar a lo realizado por
Vennard (1958), para obtener una relación potencial que relacionara el factor de fricción,
el coeficiente de la ecuación de Hazen-Williams, el diámetro, el número de Reynolds y la
viscosidad cinemática del agua. Dicha relación se ilustra en la siguiente ecuación.
100 1.852
0.2004 �
𝐶𝐻𝑊 � 1
𝑓= 0.019 0.148
𝐷 𝑅𝑒
Ecuación 2.1-26
De esta relación, se puede ver que el coeficiente de la ecuación de Hazen-Williams no es

un parámetro que depende exclusivamente del material y edad de la tubería, como se
había concebido inicialmente, sino es una característica del tubo (rugosidad y diámetro),
del flujo (número de Reynolds) y del fluido (viscosidad cinemática).
Por ello, para determinar los límites de aplicabilidad de la ecuación de Hazen-Williams,

Diskin graficó en el diagrama de Moody la Ecuación 2.1-23 (ver Figura 2.1-2) para
distintos valores del coeficiente de la ecuación de Hazen-Williams (entre 60 y 160),
diámetros que oscilaban entre 0.05 m y 1.0 m, a diferencia del estudio realizado por
18
Vennard (1958), y suponiendo una viscosidad cinemática del agua de 1.14x10-6 m2/s. A
partir de este diagrama, este investigador concluyó que esta ecuación es aplicable en una
parte de la zona de transición del diagrama de Moody (ver Figura 2.1-2), la cual está
delimitada por un rango de números de Reynolds, donde las líneas de este diagrama
presentan una pendiente de -0.148 en virtud de la Ecuación 2.1-25 (ver Figura 2.1-2).
Por lo anterior, Diskin estableció que esta ecuación solamente es aplicable para tuberías
cuyo coeficiente de Hazen-Williams se encuentre en el rango de 100 a 160, como se
ilustra en la siguiente tabla.
Ks/d Re min Re max CHw
2x10-2 2000 5000 100
1.5x10-2 2000 7500 110
1x10-2 2000 10.000 110
6x10-3 4000 20.000 120
4x10-3 8000 25.000 120
2x10-3 10.000 40.000 130
1x10-3 20.000 100.000 130
6x10-4 30.000 150.000 140
4x10-4 40.000 200.000 140
2x10-4 60.000 400.000 140
1x10-4 80.000 800.000 150
5x10-5 100.000 1´000.000 150
4x10-5 400.000 4´000.000 160
4x10-6 6´000.000 20´000.000 160
Tabla 2.1-1 Limites de aplicabilidad de la ecuación de Hazen-Williams (Diskin, 1960).
19
Efecto del uso de la ecuación de Hazen-Williams vs la ecuación de Darcy-Weisbach en tuberías de
gran diámetro (Redes matrices).
Figura 2.1-2 Ecuación de Hazen-Williams gráficada en el Diagrama de Moody (Diskin, 1960).
20
2.1.3.3 Barlow y Markland (1975)
Teniendo en cuenta las limitaciones que presentaba la ecuación de Hazen-Williams,

expuestas por Vennard (1958) y Diskin (1960), diversos investigadores intentaron corregir
dicha ecuación a fin de obtener una mayor exactitud en el cálculo de las pérdidas por
fricción en tuberías en comparación con la ecuación de Darcy-Weisbach. Barlow y
Markland (1975) en su artículo “Converting the Hazen-Williams equation to the Colebrook
function” ilustraron la forma de cómo corregir los errores de dicha ecuación empírica
modificando los parámetros de la misma, como el coeficiente que describe la rugosidad
de la tubería CHw (ver Ecuación 2.1-23). Para ello, estos investigadores reorganizaron la
ecuación de Hazen-Williams en la forma de la ecuación de Darcy-Weisbach, de la misma
forma como lo efectuó Vennard (1958), para obtener una ecuación que permitiera
calcular el factor de fricción, tal como se ilustra a continuación:
1055
𝑓= 1.85
𝐶𝐻𝑊 𝑅𝑒 0.15 𝐷0.015
Ecuación 2.1-27
Luego, gráficaron en el diagrama de Moody la Ecuación 2.1-27 y la ecuación de

Colebrook-White (ver Ecuación 2.1-11), con el propósito de comparar los factores de
fricción obtenidos por medio de estas dos ecuaciones. A partir de este análisis,
observaron que el coeficiente de la ecuación de Hazen-Williams, para un coeficiente
conocido previamente de una tubería, variaba con el aumento del caudal (Barlow y
Markland, 1975), debido a que los factores de fricción obtenidos por medio de la Ecuación
2.1-27 en el diagrama de Moody disminuían linealmente a medida que aumentaba el
número de Reynolds, lo cual no sucedía con los factores de fricción obtenidos a partir de
la ecuación de Colebrook-White. Lo anterior, se puede ver en la siguiente tabla, donde se
muestran los cálculos que realizaron Barlow y Markland (1975) para determinar los
coeficientes de la ecuación de Hazen-Williams aumentando el número de Reynolds en
proporciones de 2, 5 y 10 veces.
21
Tabla 2.1-2 El efecto del aumento del número de Reynolds en el valor del Coeficiente de Hazen-Williams, para un diámetro de 1 ft
(Barlow y Markland, 1975).
22
Adicionalmente, los investigadores en dicha comparación, establecieron los valores

correctos de la rugosidad relativa para un valor conocido del coeficiente de la ecuación de
Hazen-Williams, despejando de la ecuación de Colebrook-White dicho parámetro (ver
Ecuación 2.1-11) en términos del número de Reynolds y el factor de fricción calculado a
partir de la Ecuación 2.1-26.
1
𝑘 − 2.51
= 3.7 �10 2�𝑓 − �
𝐷 𝑅𝑒�𝑓
Ecuación 2.1-28
Los valores obtenidos, por Barlow y Markland (1975), en relación con la rugosidad relativa
se ilustran en la siguiente tabla:
Re 3.00E+03 1.00E+04 3.00E+04 1.00E+05 3.00E+05 1.00E+06

CHw VD (ft/s) 0.037 0.123 0.368 1.23 3.68 12.3
60 0.205951 0.160635 0.123875 0.090287 0.065761 0.044979
70 0.129138 0.098217 0.073177 0.050958 0.035402 0.022862
80 0.079538 0.059710 0.043211 0.028842 0.019147 0.011691
90 0.047187 0.035754 0.025360 0.016309 0.010375 0.006001
100 0.025852 0.020723 0.014646 0.009156 0.005606 0.003082
110 0.011594 0.011198 0.008163 0.005044 0.002997 0.001573
120 0.001905 0.005090 0.004203 0.002661 0.001558 0.000789
130 - 0.00111 0.00175 0.00127 0.00076 0.00038
140 - - 0.00022 0.00044 0.00031 0.00016
150 - - - - 0.00005 0.00004
Tabla 2.1-3 Ks/d calculados por Barlow y Markland (1975) en la ecuación de Colebrook-White, para un
diámetro de 1 ft.
A partir de lo anterior, Barlow y Markland (1975), en su momento, recomendaron utilizar la

ecuación de Hazen-Williams para calcular las pérdidas por fricción, sólo si se tiene en
cuenta las Tablas 2.1-2 y 2.1-3 para modificar los coeficientes de esta ecuación
dependiendo de la magnitud del número de Reynolds, a fin de minimizar los errores que
se puedan presentar en el cálculo de dichas pérdidas en comparación con la ecuación de
Darcy-Weisbach.
23
Darcy-Weisbach en tuberías de gran diámetro (Redes matrices)
2.1.3.4 Jain-Mohan-Khanna (1978)
A diferencia de Barlow y Markland (1975), Jain et al (1978) posteriormente plantearon

una nueva ecuación para reemplazar la ecuación de Hazen-Williams, debido a las
limitaciones inherentes que presentaba esta ecuación, las cuales son:
• El factor b que se utiliza para deducir la ecuación de Hazen-Williams se obtuvo

suponiendo un radio hidráulico igual a 1 ft (0.3 m) y una pendiente de la línea de
fricción de 1/1000. Por ello, dicho factor se puede ver alterado si estos parámetros
varían en su magnitud.
• El coeficiente de la ecuación de Hazen-Williams no es constante, varía
dependiendo de las características de la tubería y del número de Reynolds.
Las consecuencias de estas limitaciones se vieron reflejadas en los errores encontrados

por estos investigadores cuando se utilizaba dicha ecuación empírica en el cálculo de las
pérdidas por fricción y en la estimación del factor de fricción, los cuales oscilaban entre
±60% y ±39 % respectivamente, como se puede ver en las Tablas 2.1-4 y 2.1-5.
Por lo anterior, con base en una relación explícita para el diámetro en el flujo turbulento de
tuberías lisas (Swamee & Jain, 1976) y la ecuación de Colebrook-White para el cálculo del
factor de fricción, Jain et al (1978) plantearon una ecuación que minimizara estos errores
en el cálculo del flujo de tuberías, la cual se ilustra a continuación:
𝑉 = 143.534 𝐶𝑅 𝑟 0.6575 𝑆 0.5525
Ecuación 2.1-29
donde
V= Velocidad.
Cr=Coeficiente de rugosidad, cuya magnitud oscila entre 0 y 1, y depende del material,
diámetro de la tubería, velocidad de flujo y viscosidad cinemática del agua.
r= Radio hidráulico.
S= Pendiente de la línea de fricción.
24
Radio Porcentaje
Diámetro en Hidráulico, r, Magnitud Porcentaje de error en las
Número pulgadas en pies Pendiente error en la pérdidas por
Pendiente S (1) (milímetros) (2) (metros) (3) (4) velocidad (5) fricción (6)
1 4 (100) 1/12 (0.025) 0.33/1000 -30.74 -57
2 4(100) 1/12 (0.025) 1/1000 -27.6 -51.1
3 12 (300) 1/4 (0.76) 1/1000 -16.6 -30.7
4 72 (1800) 1.5 (0.46) 1/1000 5.5 10.2
5 72 (1800) 1.5 (0.46) 7.8/1000 14.32 26.5
6 72 (1800) 1.5 (0.46) 800/1000 32.45 60.2
Tabla 2.1-4 Porcentaje de errores en el cálculo de la velocidad y pérdidas por fricción por el uso de la
ecuación de Hazen-Williams (Jain et al, 1978).
Diámetro de la Velocidad en Porcentaje

Número Rugosidad Coeficiente tubería en metros por de error en
Pendiente S Material Tubería en metros de Hazen- milímetros segundo el factor de
(1) (2) (3) Williams (4) (5) (6) fricción (7)
Hierro fundido
1 nueva 0.00015 130 100 3 -13
Hierro fundido
2 nueva 0.00015 130 1500 0.3 26.19
Hierro fundido
3 vieja 0.0024 100 100 3 -39
Hierro fundido
4 vieja 0.0024 100 1500 0.3 32.52
5 Concreto nueva 0.000035 130 100 3 -12.92
6 Concreto nueva 0.000035 130 1500 0.3 34.92
7 Acero nueva 0.00006 140 100 3 -9.22
8 Acero nueva 0.00006 140 1500 0.3 17
Tabla 2.1-5 Porcentaje de error en la estimación del factor de fricción de la ecuación de Darcy-
Weisbach por el uso de la ecuación de Hazen-Williams (Jain et al, 1978).
Para calcular la velocidad con base en la Ecuación 2.1-29, adicionalmente Jain et al

(1978) elaboraron un nomograma para determinar el valor de Cr para distintos materiales
como el acero, hierro fundido y concreto, para un rango de diámetros de 100 a 2000 mm,
de velocidades de 0.3 a 6 m/s y temperaturas de 10ºC a 30ºC. Por lo anterior, a pesar de
este esfuerzo realizado por estos investigadores, el factor Cr de la Ecuación 2.1-29
presenta el mismo problema del coeficiente de la ecuación de Hazen-Williams, ya que
este también depende del material, diámetro, velocidad y viscosidad cinemática del agua,
lo cual dificulta determinar dicho factor con precisión debido a que se requieren ayudas
gráficas que tengan en cuenta las características de las tuberías, el tipo de flujo y de
fluido, las cuales se limitan a ciertos materiales como el acero, hierro fundido y concreto.
25
2.1.3.5 Kamand (1988)
Luego, Kamand (1988), miembro de la American Society of Civil Engineers (ASCE),

teniendo en cuenta la limitación que presentaba el coeficiente de la ecuación de Hazen-
Williams para representar la rugosidad de las tuberías independientemente del diámetro y
la velocidad de flujo, desarrolló un nomograma (ver Figura 2.1-3) que permitiera
determinar la diferencia en los cálculos de las pérdidas por fricción cuando se utiliza dicha
ecuación empírica en comparación con la ecuación de Darcy-Weisbach, para minimizar
dichos errores.
Figura 2.1-3 Nomograma que relaciona el coeficiente de la ecuación de Hazen-Williams y el factor de

fricción para diferentes velocidades y diámetros (Kamand, 1988).
Luego, el investigador llevó a cabo un análisis para determinar los factores que influyen
en los errores que se presentan en los cálculos de las pérdidas por fricción cuando se
utilizan las ecuaciones mencionadas. Revisando los valores sugeridos por diversos
investigadores para establecer el coeficiente de la ecuación de Hazen-Williams para el
PVC y el Hierro fundido, Kamand (1998) pudo observar que un factor a tener en cuenta en
estos errores es la incertidumbre que existe en sí para determinar dicho coeficiente para
un material dado, ya que existen diferencias en los valores sugeridos para este
coeficiente, como se puede ver en la Tabla 2.1-6.
26
Tabla 2.1-6 Valores de Rugosidad Absoluta y coeficientes de la Ecuación de Hazen-Williams para el

PVC y Hierro Fundido (Kamand, 1988).
Así mismo, Kamand (1988) calculando los errores en las pérdidas por fricción con la
ecuación de Darcy-Weisbach y la ecuación de Hazen-Williams para distintos diámetros y
velocidades, pudo determinar que estos oscilan entre -37.4% y 15.5%, y que estos
aumentan en la medida en que la velocidad es mayor y el diámetro es menor a 100 mm
(ver Figura 2.1-4), límite de aplicabilidad de la ecuación de Hazen-Williams en cuanto al
diámetro.
Figura 2.1-4 Errores en el cálculo de las pérdidas por fricción entre la ecuación de Hazen-Williams y la
Ecuación de Darcy-Weisbach para Hierro Fundido (Kamand, 1988).
27
Por esta razón, para minimizar los errores en los cálculos de las pérdidas por fricción en la
ecuación de Hazen-Williams, Kamand (1988) propuso utilizar el coeficiente de la ecuación
Hazen-Williams correspondiente a las condiciones de flujo promedio que se presentan en
las tuberías, ya que este varía dependiendo del material, velocidad y diámetro. Para ello,
este investigador elaboró una tabla con los valores promedio de estos coeficientes para el
PVC (ks: 0.03 mm) y Hierro fundido (ks: 0.259 mm) para un rango de diámetros de 25 a
1050 mm y rango de velocidades de 0.3 a 1.5 m/s, la cual se puede ver en mayor detalle
en la siguiente tabla.
Tabla 2.1-7 Valores del Coeficiente de Hazen-Williams promedio para PVC y Hierro Fundido para un
rango de velocidades de 0.3-1.5 m/s y diámetros de 25 a 1050 mm (Kamand, 1988).
A pesar de todos los esfuerzos realizados por Barlow y Markland (1975), Jain et al (1978)
y Kamand (1988) para mejorar la utilización de la ecuación de Hazen-Williams en la
práctica de la ingeniería, en la actualidad esta ecuación se sigue utilizando alrededor del
mundo sin tener en cuenta las limitaciones inherentes que presenta dicha ecuación
empírica y los efectos que podría tener su uso en las redes de distribución de agua
potable (RDAP) y específicamente en las redes con tuberías de gran diámetro. Por esta
razón, recientemente los investigadores Liou (1998), Bombardelli y García (2003)
realizaron algunos estudios para analizar dicho efecto.
28
2.1.3.6 Liou (1998)
Liou (1998) mostró cuantitativamente los efectos y las limitaciones que presenta la
utilización de la ecuación de Hazen-Williams en tuberías. Para ello, este investigador,
primero dedujo la relación que existe entre el coeficiente de esta ecuación empírica, el
número de Reynolds y la rugosidad relativa, con base en la ecuación de Darcy-Weisbach
(Ver Ecuación 2.1-3), la ecuación de Colebrook-White (Ver Ecuación 2.1-11) y la ecuación
de Hazen-Williams en el Sistema Internacional de unidades (Ver Ecuación 2.1-23), tal
como se ilustra a continuación:
𝜀 −0.01 −0.01 −0.08

𝐶𝐻𝑊 = 14.07𝑓 −0.54 𝑅𝑒 0.08 � � 𝜀 𝜈
𝐷
Ecuación 2.1-30
donde:
CHW= Coeficiente de Hazen-Williams.

f= Factor de Fricción.
Re= Número de Reynolds.
ε= Rugosidad absoluta de la tubería.
D= Diámetro.
𝜈= Viscosidad cinemática.
A partir de la Ecuación 2.1-30, Liou graficó (Ver Figura 2.1-5) el coeficiente de Hazen-
Williams en términos del número de Reynolds y la rugosidad relativa para diferentes
diámetros, calculando el factor de fricción a partir de la ecuación de Colebrook-White (ver
Ecuación 2.1-11), a fin de mostrar la variación que presenta este coeficiente, y establecer
los rangos donde se puede emplear dicha ecuación con base en la información de los 14
primeros experimentos que utilizaron Hazen y Williams para establecer su ecuación. Para
estos, Liou adoptó un valor constante para la rugosidad de 0.0003 m y un valor de
1.133x10-6 m2 /s para la viscosidad cinemática del agua (viscosidad que corresponde a
15.56oC). Adicionalmente, este investigador delimitó el límite que diferencia el flujo
turbulento hidráulicamente rugoso y el flujo turbulento transicional, para saber en qué tipo
de flujo se encontraban dichos experimentos, utilizando la Ecuación 2.1-16.
29
Figura 2.1-5 Variación del Coeficiente de la Ecuación de Hazen-Williams como función del número de
4
Reynolds, rugosidad relativa para agua fluyendo a 15.56ºC, y una rugosidad absoluta de 0.0003m.
Con base en la Figura 2.1-5, Liou coincide con Diskin en afirmar que la ecuación de
Hazen-Williams, como fue obtenida, se puede utilizar para un limitado número de
Reynolds y diámetros, dado que los coeficientes sugeridos para utilizar esta ecuación se
ubican exclusivamente en la zona de transición. Adicionalmente, a partir de la
Figura 2.1-5, Liou (1998) advierte que el coeficiente de la ecuación de Hazen-Williams
varía significativamente con el número de Reynolds y el diámetro, para una rugosidad
absoluta dada, lo cual imposibilita el uso de esta ecuación empírica para cualquier tipo de
flujo.
Por otro lado, Liou (1998) analizó cuantitativamente el efecto del uso de la ecuación de
Hazen-Williams en comparación con la ecuación de Darcy-Weisbach en el cálculo de la
pendiente de energía, utilizando la ecuación que este investigador dedujo para establecer
el valor de la rugosidad relativa de una tubería a partir de los valores sugeridos para el
coeficiente de la ecuación de Hazen-Williams (Ver Ecuación 2.1-31).
𝜀 0.926 𝑑 0.0093 (𝑅𝑒𝑣)0.074 � 5.74

= 3.7 �10−�0.0432𝐶 − 0.9 �
𝐷 𝑅𝑒
Ecuación 2.1-31
4
(Liou,1998).
30
A partir de este análisis, el investigador encontró que los errores en el cálculo de esta
pendiente son relativamente pequeños cuando se utiliza la ecuación de Hazen-Williams
dentro de su rango de validez. Sin embargo, si esta ecuación se aplica por fuera de dicho
rango se pueden presentar errores de hasta el 40% cuando se compara con la pendiente
obtenida con la ecuación de Darcy-Weisbach. Por lo anterior, Liou recomienda evitar la
utilización de la ecuación de Hazen-Williams para el cálculo de flujo de tuberías.
2.1.3.7 Bombardelli y García (2003)
Teniendo en cuenta los hallazgos encontrados por Liou (1998), los errores que se pueden
presentar en el cálculo de las pérdidas por fricción pueden resultar mayores en redes
matrices cuando se utiliza la ecuación de Hazen-Williams, en cuanto a la distribución de
presiones. En este sentido, Fabián Bombardelli y Marcelo García (2003) analizaron los
efectos que tiene el uso de la ecuación de Hazen-Williams posteriormente al diseño de
estas redes. Para ello, estos investigadores partieron del problema que se presentaba en
un sistema de distribución de agua potable de una ciudad construido con tuberías de
concreto después de 5 años de operación, el cual era la disminución significativa de los
coeficientes de Hazen-Williams calculados por firmas de ingeniería a partir de mediciones
de velocidad y presión en campo, utilizando esta ecuación empírica, en comparación a los
coeficientes utilizados inicialmente en el diseño de las tuberías del sistema (CHw:120).
Esto se vio reflejado en los coeficientes obtenidos en estas mediciones, ya que estos
oscilaban entre 85 y 95 para tuberías de 2.289 m de diámetro, valores que según las
tablas contenidas en el libro Hydraulic Tables (1905), correspondían a tuberías de hierro
fundido deterioradas. Lo anterior, indicaría en un principio que la red matriz en el futuro no
estaría en capacidad de suplir de agua a la ciudad, dado que se estaría perdiendo mayor
cantidad de energía que la esperada.
Después de verificar que las mediciones hechas por estas firmas no presentaban un error
significativo cuando no se consideraban las pérdidas menores, Bombardelli y García
plantearon que las posibles causas de este fenómeno podrían estar asociadas por un lado
con el uso de la ecuación de Hazen-Williams o por el aumento de la rugosidad en sí de
las tuberías de la red matriz de estudio. En relación con la primera causa, estos
investigadores encontraron, de acuerdo con los resultados de las mediciones realizadas
31
por las firmas consultoras, que el tipo de flujo que se presentaba en la operación del
sistema era flujo turbulento hidráulicamente rugoso (ver Figura 2.1-6 y Tabla 2.1-8),
debido al aumento de la demanda en la red, la presencia de tuberías de gran diámetro y
la continua operación del sistema, flujo donde no es posible aplicar la ecuación de
Hazen-Williams, dado que esta se puede utilizar solamente en la zona de transición,
como lo describieron claramente Diskin (1960) y Liou (1998), para obtener precisión en
los cálculos hidráulicos.
Por otro lado, para verificar si la rugosidad de las tuberías había aumentado en la red
matriz, estos investigadores con base en las mediciones realizadas por las firmas
consultoras, calcularon la rugosidad absoluta (ε) de la tubería principal de la red matriz de
estudio (Ver Tabla 2.1-8), la cual se subdivide en tres tramos, utilizando la ecuación de
Swamme y Jain (Ver Ecuación 2.1-19) para un diámetro de 2.286 m, caudal igual a 6.795
m3/s, velocidad correspondiente a 1.655 m/s y un número de Reynolds igual a 3.8x106 .
Tabla 2.1-8 Rugosidades absolutas obtenidas por Bombardelli y García (2003) para la tubería principal
de la red matriz de estudio.
Comparando estas rugosidades con los valores esperados de rugosidad absoluta para el
concreto (0.0033 y 0.0005 m) y llevando a cabo un análisis de sensibilidad de la variación
del coeficiente de Hazen-Williams (para un n de Manning de 0.014, valor esperado para
una tubería de concreto) debido a errores en la medición de presiones y velocidades,
estos investigadores lograron establecer que las rugosidades calculadas y los coeficientes
obtenidos por las firmas consultoras estaban por fuera de los rangos obtenidos. Con base
en ello, estos investigadores concluyeron que en realidad la rugosidad de las tuberías si
era alta pero no en la medida como lo indicaba la ecuación de Hazen-Williams. Por lo
anterior, estos investigadores recomiendan utilizar la ecuación de Darcy-Weisbach en el
diseño de redes matrices debido a que esta incluye todos los tipos de flujo, y de esa
manera se evita mal interpretar los resultados de la ecuación de Hazen-Williams cuando
se utiliza por fuera de su rango de aplicabilidad.
32
Efecto del uso de la ecuación de Hazen-Williams vs la ecuación de Darcy-Weisbach en tuberías
de gran diámetro (Redes matrices).
Figura 2.1-6 Coeficientes de Hazen-Williams obtenidos de las firmas consultoras en la red matriz analizada por Bombardelli y García (2003).
33
2.2 Redes Matrices

Las redes de distribución de agua potable (RDAP) son una serie de obras hidráulicas que
permiten el suministro de agua potable a una población específica. Estas redes están
compuestas por dos partes, las cuales son: la red matriz, parte encargada de conducir el
agua desde la planta de tratamiento e incluye los siguientes componentes: fuentes de
abastecimiento, obras de captación, obras de aducción y conducción para el transporte de
las aguas sin tratar hasta las plantas de tratamiento, plantas de tratamiento y tanques de
almacenamiento. La segunda parte son las redes de distribución secundaria o terciaria,
cuya función es suministrar el agua potable desde los tanques de almacenamiento hasta
los usuarios de la red.
Con base en lo anterior, las redes matrices, hidráulicamente, se clasifican como redes
abiertas, dado que estas no presentan ningún circuito cerrado en el sistema. Por ello, el
comportamiento hidráulico de estas se basa en el balance de cantidad, donde se busca
que en cada una de las uniones o nudos de la misma se conserve la masa y en cada uno
de los tanques de almacenamiento o embalses se suministre el caudal demandado para
alimentar las redes secundarias o terciarias. Por ello, las ecuaciones que se utilizan en
este tipo de redes para calcular el flujo en las tuberías son:
- Ecuación de Continuidad en las Uniones
𝑁𝑇
� 𝑄𝑖𝑢 − 𝑄𝐿𝑢 = 0
𝑖=1
Ecuación 2.2-1
donde
-Qiu= Caudal en la tubería iu (se toma positivo si llega al nodo y negativo si sale del
mismo).
-NT= Número de tuberías que llegan a una determinada unión.
- QLu= Caudal demandado en la unión.
34
- Ecuación para calcular las pérdidas de energía en las tuberías:
𝑙𝑖𝑢 𝑄𝑖𝑢2
|𝑍𝑖 − 𝑍𝑢| = (𝑓𝑖𝑢 + � 𝑘𝑚𝑖𝑢)
𝑑𝑖𝑢 2𝑔𝐴𝑖𝑢2
𝑙𝑖𝑢
(𝑓𝑖𝑢 + ∑ 𝑘𝑚𝑖𝑢)
𝑘𝑖𝑢 = 𝑑𝑖𝑢
2𝑔𝐴𝑖𝑢2
Ecuación 2.2-2
|𝑍𝑖 − 𝑍𝑢| = 𝑘𝑖𝑢𝑄𝑖𝑢2
Ecuación 2.2-3
donde
-Zu= Altura piezométrica en la unión.
- Zi=Altura piezométrica en el tanque.
- Viu= Velocidad de flujo de la tubería.
- diu= Diámetro de la tubería.
- liu= Longitud de la tubería.
-fiu= Factor de friccion de Darcy para la tubería iu.
-kmiu= Coeficiente de pérdidas menores en la tubería.
- Kiu= Parámetro característico de la tubería que incluye los factores de pérdidas por
fricción y pérdidas menores.
Con base en estas ecuaciones, para analizar la hidráulica de una red matriz se tiene
como incógnitas las alturas piezométricas en las uniones y los caudales que llegan a cada
uno de los tanques de almacenamiento. Para ello, es necesario en una primera iteración
suponer alturas en dichas uniones y luego corregirlas teniendo en cuenta la metodología
propuesta por Feaherstone en 1981, la cual se resume en la siguiente ecuación:
35
2(∑ 𝑄𝑖𝑗 − 𝑄𝐿𝑢)

∆𝑍𝑢𝑗 =
𝑄𝑖𝑗
∑
𝑍𝑖 − 𝑍𝑢𝑗
Ecuación 2.2-4
donde Δzuj es la corrección de la altura piezométrica supuesta en las uniones y Qij el

caudal que pasa por las tuberías que llegan al nudo, el cual se obtiene utilizando la
siguiente relación entre las pérdidas por fricción y el caudal.
2�2𝑔𝑑ℎ𝑓 𝑘𝑠 2.51𝑣√𝑙
𝑄=− 𝑙𝑜𝑔10 � + �
√𝑙 3.7𝑑 𝑑�2𝑔𝑑ℎ𝑓
Ecuación 2.2-5
Sin embargo, para propósitos de analizar redes matrices de gran tamaño, esta
metodología no es tan eficiente, dado que se tendrían que corregir una a una las alturas
piezométricas en las uniones de la red. Por tal razón, para esta tesis, se utilizarán
herramientas computacionales como REDES y EPANET, que se basan en el Método del
Gradiente, el cual permite corregir de forma simultánea las alturas piezométricas y los
caudales que pasan por cada una de las tuberías, sin importar la configuración que
presenten este tipo de redes.
2.3 Programa REDES

El programa REDES “es una herramienta de simulación de sistemas de tuberías con flujo
a presión” 5 creada para realizar el análisis y diseño optimizado de redes de distribución de
agua potable (RDAP). Este programa nació en la Universidad de los Andes en el año de
1998 y ha sido el resultado de diversas investigaciones y proyectos de grado del
Departamento de Ingeniería Civil y Ambiental de la Universidad de los Andes, bajo la
dirección del ingeniero Juan G. Saldarriaga. Este programa ofrece diversas opciones de
cálculo de RDAP para el usuario, entre las cuales se encuentran:
5
(Saldarriaga, 2007).
36
• Calculo hidráulico estático: Esta opción permite al usuario calcular las alturas
piezométricas y presiones, utilizando el Método del Gradiente, desarrollado por
Ezio Todini y Enda O´Conell, por medio de las ecuaciones de Hazen-Williams y
Darcy-Weisbach.
• Calculo hidráulico en período extendido: Esta opción permite al usuario visualizar
la variación de las alturas piezométricas y presiones, definiendo previamente
patrones consumo en cada uno de los nodos de la RDAP.
• Cálculo Calidad de agua.
• Calibración de RDAP.
• Análisis de Resiliencia de RDAP.
• Diseños de RDAP con métodos modernos, como diseño rápido, superficie óptima
del gradiente hidráulico (SOGH), Algoritmos Genéticos entre otros.
Teniendo en cuenta que en esta tesis se diseñaran redes matrices utilizando la

metodología de diseño SOGH, a continuación se describe en detalle dicha metodología.
2.3.1 Metodología superficie óptima del Gradiente hidráulico (SOGH)
Esta metodología de diseño fue desarrollada por Susana Ochoa en el año 2009 como
proyecto de grado de maestría de Ingeniería Civil en la Universidad de los Andes, bajo la
asesoría del ingeniero Juan G. Saldarriaga. A diferencia de la metodología de Superficie
Óptima, la cual establece que el mínimo costo de construcción de una red se obtiene
cuando la línea piezométrica forma "una curva cóncava hacia arriba con una flecha del
15% en el centro con respecto a la línea recta que une las alturas totales al inicio y al final
de la serie de tuberías" 6 , la metodología SOGH estima dicha flecha óptima como función
de la distribución de la demanda, la demanda total de RDAP, longitud total de las tuberías
de la red y la función de costos. Adicionalmente, es importante mencionar que esta última
metodología es aplicable en redes que presenten cualquier tipo de topología.
6
(Saldarriaga, 2007)
37
A priori al utilizar la metodología SOGH en el diseño optimizado de RDAP, el usuario debe

fijar algunos aspectos básicos de diseño como: la presión mínima requerida, función de
costos (depende de los materiales), diámetros comerciales disponibles, ecuación de
pérdidas por fricción (Darcy-Weisbach o Hazen-Williams), número de iteraciones y error
admisible en el cálculo de la superficie óptima del gradiente hidráulico.
Teniendo en cuenta lo anterior, el usuario debe seguir el siguiente procedimiento para

llevar a cabo el diseño de RDAP por medio de la metodología SOGH 7.
• Asignación inicial de diámetros proporcionales a la distancia topológica.
Dicha distancia corresponde a la ruta más corta que puede recorrer el agua en la
red para llegar a un nodo específico desde una fuente de abastecimiento. Para
ello, es necesario hacer una simulación hidráulica. Luego, se efectúa una
asignación de diámetros en las tuberías proporcional a la distancia topológica: a
las tuberías que presenten menor distancia se les asignará un diámetro grande y a
las tuberías con mayor distancia topológica se les asignará un diámetro pequeño.
Es importante aclarar, que la distancia topológica de cada una de las tuberías se
calcula como el promedio de las distancias topológicas del nodo inicial y final.
• Estimación de la flecha óptima de la curva línea de gradiente hidráulico (LGH) vs

Abscisa.
Para calcular dicha flecha, primero se calcula el centroide de demandas de la

RDAP (𝑥) utilizando la siguiente ecuación.
∑𝑁𝑁
𝑖=0 𝑞𝑖 ∗ 𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎𝑖
𝑄𝑡𝑜𝑡𝑎𝑙
𝑥=
𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎𝑚𝑎𝑥𝑖𝑚𝑎
Ecuación 2.3-1
donde:
NN= Corresponde al número de nodos de la RDAP.

qi= Es el caudal demandado en cada uno de los nodos.
7
(Ochoa, 2009)
38
dtopologicai= Es la distancia topológica de cada uno de los nodos con respecto a

la fuente de abastecimiento.
dtopologicamaxima= Es la distancia topológica máxima que se presenta en la red.
Posteriormente, se calcula el coeficiente de uniformidad 7 (CU7) de la RDAP.

Para ello, es necesario dividir la red en dos subconjuntos. El primer subconjunto
pertenece a los nodos cuya distancia topología es inferior al centroide de
demandas de la RDAP y el último subconjunto corresponde a los nodos con
distancia topológica mayor a dicho centroide. Teniendo en cuenta lo anterior, se
calculan los centroides de cada uno de los subconjuntos para encontrar dicho
coeficiente como se puede ver a continuación.
𝐿1 𝐿2
𝐶𝑈7 = 𝑥1 ∗ � � + 𝑥2 ∗ � �
𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎𝑚𝑎𝑥𝑖𝑚𝑎 𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎𝑚𝑎𝑥𝑖𝑚𝑎
Ecuación 2.3-2
donde:
L1= Longitud del primer subconjunto (𝐿1 = 𝑥 ∗ 𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎𝑚𝑎𝑥𝑖𝑚𝑎)
L2= Longitud del segundo subconjunto (𝐿2 = 𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎𝑚𝑎𝑥𝑖𝑚𝑎 − 𝐿1)
𝑥1 𝑦 𝑥2= Centroides de cada uno de los subconjuntos, los cuales se calculan con
la siguiente ecuación
∑𝑁𝑠
𝑖=0 𝑞𝑖 ∗ |𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎 𝑖 − 𝑥 ∗ 𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎𝑚𝑎𝑥𝑖𝑚𝑎|
∑𝑁𝑠𝑖=0 𝑞𝑖
𝑥𝑖 =
𝑑𝑡𝑜𝑝𝑜𝑙𝑜𝑔𝑖𝑐𝑎𝑚𝑎𝑥𝑖𝑚𝑎
Ecuación 2.3-3
donde Ns es el número de nodos del subconjunto.
Luego, se calcula la flecha óptima con base en el centroide de demandas y el

Coeficiente de Uniformidad 7, para una relación Q2/L3=1x10-9 y un exponente de la
función de costos de 1.46, utilizando la siguiente ecuación.
39
𝐹𝑙𝑒𝑐ℎ𝑎 = 𝑎 + 𝑏 ∗ 𝑥 + 𝑐 ∗ 𝐶𝑈7 + 𝑑 ∗ (𝐶𝑈7)2

Ecuación 2.3-4
donde:
a=0,435521465
b=-0,176612805
c=-0,977366227
d=0,906254447
Después, se determina la flecha óptima para el exponente de la función de costos

(n) que se tenga, utilizando la siguiente ecuación:
𝐹𝑙𝑒𝑐ℎ𝑎 𝑜𝑝𝑡𝑖𝑚𝑎(𝑛) = 𝛼 ∗ 𝑛2 + 𝛽 ∗ 𝑛 + 𝛾
Ecuación 2.3-5
donde:
n= Exponente de la función de costos.
𝛼 = −0,1134 ∗ 𝐹𝑙𝑒𝑐ℎ𝑎(1.46) + 0,0032
𝛽 = 0,6443 ∗ 𝐹𝑙𝑒𝑐ℎ𝑎(1.46) + 0,0043
𝛾 = 0,2835 ∗ 𝐹𝑙𝑒𝑐ℎ𝑎(1.46) − 0,0111
𝐹𝑙𝑒𝑐ℎ𝑎(1.46)= Flecha óptima para un exponente de la función de costos de 1.46.
Teniendo en cuenta esta última flecha óptima, se procede a calcular la flecha

óptima para la relación real Q2/L3, utilizando la siguiente ecuación.
𝑄2
𝐹𝑙𝑒𝑐ℎ𝑎 𝑜𝑝𝑡𝑖𝑚𝑎 = 𝑎 ∗ ln � �+𝑏
𝐿3
Ecuación 2.3-6
40
donde:
𝑎 = 0,00868 ∗ 𝐹𝑙𝑒𝑐ℎ𝑎(𝑛) + 0,00066
𝑏 = 1,18069 ∗ 𝐹𝑙𝑒𝑐ℎ𝑎(𝑛) + 0,01345
𝐹𝑙𝑒𝑐ℎ𝑎(𝑛)= Flecha óptima para una relación Q2/L3=1x10-9 y exponente de la

función de costos que se tenga.
Q= Caudal total demandado (Q), el cual es la sumatoria de todos los caudales

demandados en cada uno de los nodos de la red.
L= Longitud total de la red abierta (L), la cual es la sumatoria de las longitudes de

todos los tubos que pertenecen a la red.
• Con base en la asignación inicial de diámetros, realizar una simulación hidráulica

para determinar caudales, sentidos de flujo en tuberías y calcular las distancias
topológicas de los nodos a las fuentes de abastecimiento.
• Cálculo de la superficie óptima de gradiente hidráulico.
Para calcular dicha superficie, primero se asigna a cada nodo de la red la fuente
principal de abastecimiento, lo cual se hace inyectando un trazador en cada fuente
de alimentación y calculando su trayectoria a través de la red. Luego, se detectan
los sumideros de la red, es decir, aquellos nodos que se encuentran al final de una
o varias rutas de la red. Posteriormente, se ordenan los sumideros de forma
descendente según su distancia topológica, para recorrer cada uno de los
sumideros y así poder ajustar los coeficientes de la ecuación cuadrática que se
utilizará para calcular la LGH óptima en cada uno de los recorridos, teniendo en
cuenta que a cada sumidero se le debe asignar una LGH óptima equivalente a su
cota más la presión mínima requerida. Dicho ajuste se ilustra continuación.
𝐿𝐺𝐻(𝑑) = 𝛼 ∗ 𝑑2 + 𝛽 ∗ 𝑑 + 𝛾
Ecuación 2.3-7
41
donde:
d= Distancia topológica del nodo i.

α,β,γ= Son los coeficientes de la ecuaci
ón y se encuentran con las siguientes
ecuaciones:
𝐿𝐺𝐻𝑓𝑢𝑒𝑛𝑡𝑒 𝑎𝑏𝑎𝑠𝑡𝑒𝑐𝑖𝑚𝑖𝑒𝑛𝑡𝑜 − 𝐿𝐺𝐻𝑠𝑢𝑚𝑖𝑑𝑒𝑟𝑜𝑖

𝛼 = 4 ∗ 𝐹𝑙𝑒𝑐ℎ𝑎 𝑜𝑝𝑡𝑖𝑚𝑎 ∗
𝑑𝑠𝑢𝑚𝑖𝑑𝑒𝑟𝑜𝑖 2
Ecuación 2.3-8
𝐿𝐺𝐻𝑓𝑢𝑒𝑛𝑡𝑒 𝑎𝑏𝑎𝑠𝑡𝑒𝑐𝑖𝑚𝑖𝑒𝑛𝑡𝑜 − 𝐿𝐺𝐻𝑠𝑢𝑚𝑖𝑑𝑒𝑟𝑜𝑖

𝛽 = −(1 + 4 ∗ 𝐹𝑙𝑒𝑐ℎ𝑎 𝑜𝑝𝑡𝑖𝑚𝑎) ∗
𝑑𝑠𝑢𝑚𝑖𝑑𝑒𝑟𝑜𝑖
Ecuación 2.3-9
𝛾 = 𝐿𝐺𝐻𝑓𝑢𝑒𝑛𝑡𝑒 𝑎𝑏𝑎𝑠𝑡𝑒𝑐𝑖𝑚𝑖𝑒𝑛𝑡𝑜
Ecuación 2.3-10
• Asignación de las pérdidas de energía objetivo a cada tramo de tubería.
A partir de la LGH ideal calculada para cada uno de los nodos de la red a partir de
la Ecuación 2.3-7, se calculan las pérdidas de energía objetivo para cada tramo de
la tubería, encontrando la diferencia de la LGH ideal de los nodos inicial y final de
cada una de ellas.
• Calculo diámetros continuos de las tubería según las pérdidas objetivo.
Teniendo en cuenta los caudales de consumo de cada uno de los nodos de la red,
se calculan los caudales de diseño de cada una de las tuberías de la red. Con ello,
se procede a calcular los diámetros de cada una de ellas utilizando la ecuación de
Darcy-Weisbach en conjunto con la ecuación de Colebrook-White o la Ecuación de
Hazen-Williams.
• Realizar una simulación hidráulica, para calcular las pérdidas de energía reales
obtenidas según los nuevos diámetros calculados para cada una de las tuberías
de la red.
42
• Calcular el error entre las pérdidas reales y objetivo. Si se cumple con el error
admisible se continua con la siguiente etapa, de lo contrario se regresa hacer una
nueva simulación hidráulica (Etapa 3).
• Realizar una simulación hidráulica para la configuración de diámetros que se tiene
al final de la anterior etapa, con ello se determinan los caudales en las tuberías y
presiones en todos los nodos.
• Si alguno de los nodos de la red no cumple con la presión mínima, es necesario
corregir los nodos con cotas muy elevadas sino se continua con la siguiente etapa.
• Aproximación de los diámetros continuos a diámetros comerciales.
• Aplicación de programación por restricciones para aumentar diámetros para
cumplir con los requerimientos de presión mínima.
• Aplicación de programación por restricciones para disminuir diámetros de algunas
tuberías cumpliendo con los requisitos de presión mínima, para disminuir el costo
de las tuberías, el cual se calcula con la siguiente ecuación:
𝐶𝑇𝑢𝑏𝑜 = 𝐾 ∗ 𝐿 𝑇𝑢𝑏𝑜 ∗ 𝐷 𝑥 𝑇𝑢𝑏𝑜

Ecuación 2.3-11
donde D es el diámetro de la tubería, Ctubo es el costo de la tubería, Ltubo la

longitud de la tubería, K y x son los coeficientes obtenidos a partir de una
regresión potencial.
2.4 Programa EPANET
Es una herramienta computacional desarrollada por la Agencia de Protección Ambiental

de los Estados Unidos (Environmental Protection Agency), la cual permite analizar las
redes de distribución de agua potable, utilizando el Método del Gradiente, definiendo
claramente la topología de la red, los niveles en los tanques o embalses y las
características de las tuberías como la rugosidad, el diámetro, el coeficiente de pérdidas
menores, entre otros factores. Igualmente, permite monitorear la calidad de agua en la
red y simular el comportamiento dinámico de una RDAP asignando patrones de consumo
en cada uno de los nudos de la misma. A diferencia del programa REDES, esta
herramienta no permite diseñar, calibrar y llevar a cabo análisis de Resiliencia.
43
2.5 MATERIALES DE TUBERÍAS

Teniendo en cuenta que en las ecuaciones de resistencia fluida la rugosidad de las
tuberías es un parámetro indispensable para el cálculo de las pérdidas por fricción, se
consultaron los valores de estas para las rugosidades absolutas y los coeficientes de
Hazen-Williams en la normatividad colombiana y específicamente en las normas de
diseño de acueductos de las Empresas Públicas de Medellín (2009).
2.5.1 Rugosidad absoluta (ks)
Si en el análisis de redes matrices, se utiliza la ecuación de Darcy-Weisbach en conjunto

con la ecuación de Colebrook-White, es necesario conocer el valor de la rugosidad
absoluta ya que este parámetro en cierta medida puede afectar la hidráulica de la red,
debido a que éste puede aumentar las pérdidas de energía sobre todo por el rozamiento
entre las paredes de la tubería y el fluido. Los valores de las rugosidades absolutas para
distintos materiales se obtuvieron de la normas de las Empresas Públicas de Medellín
(EPM), los cuales se ilustran en la siguiente tabla:
Material Rugosidad absoluta ks (mm)
Acero comercial 0.45
CCP 0.12
Hierro dúctil 0.25
GRP 0.030
Polietileno 0.007
PVC 0.0015
Tabla 2.5-1 Rugosidades absolutas contenidas en las normas de EPM (2009).
44
2.5.2 Coeficientes de Hazen-Williams (CHw)
Al emplear la ecuación de Hazen-Williams en el análisis de redes matrices, es

indispensable conocer el coeficiente de esta ecuación para calcular las pérdidas por
fricción en cada una de las tuberías de la red, ya que este es un factor que depende de la
tubería (rugosidad y diámetro), el flujo (número de Reynolds) y el fluido (viscosidad
cinemática). Los valores de estos coeficientes para distintos materiales se obtuvieron de
las normas de las Empresas Públicas de Medellín (EPM), y se muestran a continuación:
Material Condición Diámetro (in) CHW

d≥12 120
Acero soldado Constante 8≤d≤10 119
4≤d≤6 118
d≥24 113
Acero bridado Constante 12≤d≤20 111
4≤d≤10 107
Formaleta de acero Todos 140

Concreto Formaleta de madera Todos 120
Centrifugado Todos 135
PVC Constante Todos 150
Con revestimiento
Hierro dúctil 140 - 150
interior de cemento
Polietileno 150
GRP
150
Tabla 2.5-2 Coeficientes de Hazen-Williams contenidos en las normas de EPM (2009).
45
Capítulo 3. METODOLOGÍA PARA DETERMINAR EL

EFECTO DEL USO DE LA ECUACIÓN DE HAZEN-
WILLIAMS
En este trabajo se analizará el efecto del uso de la ecuación de Darcy-Weisbach frente a
la ecuación de Hazen-Williams en tuberías de gran diámetro (Redes Matrices). Para ello,
se llevará a cabo la siguiente metodología, a fin de calcular la variación de los
coeficientes de Hazen-Williams en las redes de estudio, determinar su efecto en los
errores en las pérdidas por fricción y distribución de presiones, y su correspondiente
corrección para minimizar estos errores.
INICIO
1. Seleccionar Redes Matrices de Estudio
2. Definir los materiales a analizar en las ecuaciones de Resistencia Fluida.
3. Para estos materiales, establecer las rugosidades para las ecuaciones de

Resistencia fluida, según las normas.
4. Correr la hidráulica de la red estudio, tanto para la ecuación de Hazen-

Williams y Darcy-Weisbach para estos materiales.
Obtener los resultados de la red, relacionados con la presión en los nudos,

pérdidas por fricción y caudales en las tuberías.
5. Calcular errores en cuanto a las pérdidas por fricción, presiones y los coeficientes de la
Ecuación de Hazen-Williams en las redes de estudio, entre los resultados de la Ecuación de
Darcy-Weisbach y Ecuación de Hazen-Williams.
No
Errores pérdidas>5% 8. Rugosidades FIN
materiales.
Si
6. Corregir los valores de los coeficientes de Hazen-Williams para cada uno de

los materiales, para un determinado rango de número de Reynolds, velocidad
y diámetros.
7. Diseñar las redes de estudio con la Ecuación de Hazen-Williams
Figura
Figura2.5-1.
3.0-1Diagrama
Diagramade
deFlujo
FlujoMetodología
Metodología.
46
Por último, se validaron los diseños obtenidos con base en los nuevos coeficientes de
Hazen-Williams calculados a partir de la Figura 3.0-1, realizando comprobaciones de
diseño utilizando la ecuación de Darcy-Weisbach, para establecer la precisión de los
mismos en cuanto a la distribución de presiones y errores en las pérdidas por fricción.
Igualmente, se compararan los resultados de estos diseños con los obtenidos a partir de
los coeficientes de Hazen-Williams sugeridos en la normativa colombiana.
3.1 Selección Redes de Estudio.

A continuación se describen las redes matrices seleccionadas para desarrollar la
metodología propuesta.
3.1.1 Red matriz Medellín
Es una red matriz hipotética muy parecida a la red real (Ver Anexo 1), que se encuentra
ubicada como su mismo nombre lo indica, en la ciudad de Medellín, la cual presenta las
siguientes características:
• Está conformada por 685 nudos, que presentan una altura promedio de 1650
msnm.
• Presenta 768 tuberías, donde el 68% de estas presentan diámetros superiores a
los 500 mm y el 32% de estas son menores a este diámetro.
• Está conformada por 3 embalses, Ayurá, Manantiales, Villa Hermosa, los cuales
presentan una línea de gradiente hidráulico de 1760, 1758 y 1700 msnmm,
respectivamente.
• No presenta válvulas u otro tipo de accesorios.
• El funcionamiento de la red es por gravedad (no incluye bombas).
• Presenta una demanda máxima total de 8 m3/s, distribuidos de la siguiente manera
Ayurá 5 m3/s y Manantiales 3 m3/s.
• Topografía irregular cotas desde hasta 1735 msnm hasta 1760 msnm.
47
3.1.2 Red Hanói
Es una Red matriz hipotética (Ver anexo 2), que presenta las siguientes características en
cuanto a su topología.
• Está conformada por 34 tuberías, que presentan diámetros desde 12 pulgadas

hasta 40 pulgadas.
• Está conformada por 31 nudos, que presentan una demanda de 830 m3/h.
• Presenta únicamente una fuente de abastecimiento, con una línea del gradiente
hidráulico de 100 msnm.
• El funcionamiento de la red es por gravedad (no incluye bombas).
• No presenta válvulas u otro tipo de accesorios.
• No existe variación en la cota de los nodos de esta red (Topografía plana).
3.1.3 Red matriz analizada por Bombardelli y García (2003)
Con base en el estudio realizado por Bombardelli y Garcia (2003) en una red matriz de los
Estados Unidos, en este trabajo se analizará solamente el efecto de la variación del
coeficiente de Hazen-Williams en el cálculo las pérdidas por fricción en la tubería principal
de esta red, por la falta de información de esta red. Dicha tubería presenta las
características expuestas en la Tabla 2.1-8.
3.2 Selección de materiales
Teniendo en cuenta que en el medio colombiano, se utilizan materiales como el acero,

concreto (CCP) y GRP para diseñar las redes matrices, se seleccionaron estos materiales
para analizar el efecto del uso de la ecuación de Darcy-Weisbach vs. la ecuación de
Hazen-Williams en este tipo de redes.
48
3.3 Selección Rugosidad de materiales

Con base en las normas de las Empresas Públicas de Medellín para el diseño de RDAP
(2009), se seleccionó la rugosidad de los materiales escogidos para llevar a cabo este
trabajo, los cuales se ilustran en la siguiente tabla. Es importante aclarar que según estas
normas el diámetro mínimo que se puede utilizar para estos materiales es igual a 300
mm, para el diseño de redes matrices.
Material Ks (mm) CHw
CCP 0.12 120
GRP 0.03 150
Acero 0.45 120
Tabla 3.3-1. Rugosidad de los materiales seleccionados para el estudio.
3.4 Procedimiento calculo hidráulica de las redes
Teniendo en cuenta las características topológicas de cada una de las redes de estudio
seleccionadas, a continuación se describe el procedimiento que se realizó en el desarrollo
de este trabajo para calcular la hidráulica de estas, a fin de obtener la presión en los
nudos, las pérdidas por fricción y los caudales en cada una las tuberías.
• Definir la rugosidad absoluta y los coeficientes de Hazen-Williams para cada una

de las redes a analizar. Lo anterior, se llevó a cabo modificando en los archivos
.red o .inp de las redes seleccionadas, la rugosidad en todas las tuberías con base
en lo descrito en el Numeral 3.3.
• Cargar la red a analizar en la herramienta computacional que se vaya a utilizar
para calcular la hidráulica de la red. Para ello, es importante tener en cuenta que el
programa EPANET acepta solamente archivos que tengan extensión .net o en su
defecto .inp. Así mismo, para el programa REDES, el archivo debe contener la
extensión .red o .inp para cargar la red.
• Establecer la viscosidad cinemática del agua, con base en la temperatura del
agua. Para este trabajo se adoptó una viscosidad cinemática de 1.14x10-6 m2/s,
correspondiente a unta temperatura de 15ºC.
49
• Correr la hidráulica de cada una de las redes en las herramientas computacionales

para cada una de las ecuaciones de resistencia fluida, de acuerdo con los valores
de rugosidad modificados en los archivos de entrada de las redes seleccionadas.
Figura 3.4-1 Selección ecuación de resistencia fluida en el programa EPANET.
Figura 3.4-2 Interfaz gráfica del programa EPANET para calcular la hidráulica de la red.
50
Figura 3.4-3 Interfaz gráfica del programa REDES para calcular la hidráulica de la red. Resultados
hidráulica.
• Posteriormente, para obtener los resultados del análisis hidráulico de las redes
seleccionados, en cuanto a las presiones en los nudos, pérdidas por fricción (hf),
factores de fricción (f) y velocidades en cada una de las tuberías (v), es necesario
seleccionar las ventanas de despliegue de resultados de los programas EPANET y
REDES, tal como se ilustra en las Figuras 3.4-4, 3.4-5, 3.4-6 y 3.4-7.
51
Figura 3.4-4 Interfaz programa REDES para obtener factores hidráulicos de la red.
Figura 3.4-5 Interfaz programa EPANET para seleccionar resultados de los nudos o tuberías.
Figura 3.4-6 Interfaz programa EPANET para seleccionar los resultados de la hidráulica de la red.
52
Figura 3.4-7. Interfaz gráfica programa de EPANET despliegue de resultados.
3.5 Cálculo de la variación del coeficiente de Hazen-Williams y

su efecto en las pérdidas por fricción y presiones.
Luego de haber obtenido los resultados de cada una de las redes de estudio utilizando la
ecuación de Darcy-Weisbach y la ecuación de Hazen-Williams, se procede a calcular la
variación de los coeficientes de Hazen-Williams y su efecto en los errores en el cálculo de
las pérdidas por fricción y presiones, tal como se ilustra a continuación:.
• Variación coeficiente de Hazen-Williams
A partir de la Ecuaciones 2.1-26 y 2.1-30, propuestas por Diskin (1960) y Liou

(1998), se puede ver claramente que el coeficiente de la ecuación de Hazen-
Williams no es un parámetro constante, ya que este depende de las características
de la tubería (rugosidad y diámetro), del flujo (número de Reynolds) y del fluido
(viscosidad cinemática). Por lo anterior, se calcularon los coeficientes de Hazen-
Williams teóricos, llevando a cabo un procedimiento muy similar al efectuado por
Liou (1998), para analizar la variabilidad de este parámetro con base en los
resultados obtenidos en las redes de estudio. Lo anterior, se ilustra en mayor
detalle en la Figura 3-5.1 para una red específica:
53
Inicio
Leer hfi, vi, di, Li de los resultados

de la hidráulica utilizando la
Ecuación de Darcy-Weisbach.
i=i+1
NO
i>Número
tuberías
SI
Calcular CHwi en la
Ecuación 2.1-23
NO
i>Número
de tuberías i=i+1
SI
Imprimir todos CHw
Fin
Figura 3.5-1 Diagrama de Flujo para analizar la variación del coeficiente

de la ecuación de Hazen-Williams.
Posteriormente, se comparan los coeficientes de Hazen-Williams calculados con el

coeficiente sugerido para cada una de los materiales seleccionados en el Numeral
3.3, con el propósito de analizar la variabilidad que presentan estos coeficientes.
• Errores en pérdidas por fricción:
Teniendo en cuenta la variabilidad del coeficiente de Hazen-Williams,

posteriormente se calculan los errores en las pérdidas por fricción, para determinar
el efecto del uso de la ecuación de Hazen-Williams en las redes de estudio. Para
54
ello, se calculó para cada tubería el error relativo, tomando como valor teórico los
resultados obtenidos por medio de la ecuación de Darcy-Weisbach, como se
muestra a continuación:
𝑆𝐻𝑤 − 𝑆𝐷𝑤
𝐸𝑟𝑟𝑜𝑟 =
𝑆𝐷𝑤
Ecuación 3.5-1
A partir de la Ecuación 3.5-1, se busca detectar los errores que produce la

ecuación de Hazen-Williams en la sobreestimación (Error positivo) y en la
subestimación de las pérdidas de energía (Error negativo), con especial énfasis en
este último, dado que en la realidad se estaría perdiendo mayor energía y por lo
tanto la red diseñada a partir de la ecuación de Hazen-Williams no estaría en
capacidad de suministrar agua a una determinada población.
• Errores presión en los nudos:
Dado que en las redes matrices es importante cumplir con la presión mínima de
diseño exigida en las normas, se calcularon errores relativos entre los resultados
obtenidos en cuanto a la presión utilizando la ecuación de Darcy-Weisbach (DW) y
Hazen-Williams (HW), para determinar en qué magnitud no se está cumpliendo
dicha presión:
𝑃𝑟𝑒𝑠𝑖ó𝑛(𝐻𝑊) − 𝑃𝑟𝑒𝑠𝑖ó𝑛 (𝐷𝑊)
𝐸𝑟𝑟𝑜𝑟 𝑝𝑟𝑒𝑠𝑖ó𝑛 =
𝑃𝑟𝑒𝑠𝑖ó𝑛 (𝐷𝑊)
Ecuación 3.5-2
3.6 Corrección Coeficientes de Hazen-Williams
Si los errores calculados en el inciso anterior con los valores sugeridos para los
coeficientes de Hazen-Williams, no permiten cumplir con la presión mínima de diseño, y si
estos son mayores al 5% para las pérdidas por fricción, se propone el procedimiento
descrito en la Figura 3.6-1, dada una rugosidad absoluta de cualquier material, a fin de
cumplir dichos errores.
55
INICIO
Leer 𝜈, ks y CHw sugerido
Definir diámetro
mínimo y máximo
a analizar
Generar Reynolds desde

1000 hasta 30.000.000
Calcular
𝑅𝑒𝑖 ∗ 𝑑𝑖
𝑉𝑖 =
𝜈
Calcular fdw en la
ecuación (2) y fhw en
la ecuación (6)
Calcular
𝐟𝐡𝐰 − 𝐟𝐝𝐰
𝐞𝐫𝐫𝐨𝐫 𝐡𝐟 =
𝐟𝐝𝐰
NO
j>Reynolds
máximo j=j+1
SI
NO
i>diámetro
máximo i=i+1
SI
Para el diámetro mínimo

y máximo establecer el
error hf máximo
SI
NO
|𝐞𝐫𝐫𝐨𝐫 𝐡𝐟 𝐦𝐚𝐱| Imprimir CHw, y rango número de
CHw=CHw+1 < 𝟓% Reynolds para error hf max
correspondiente al diámetro mínimo
FIN
Figura 3.6-1 Diagrama de Flujo corrección coeficientes de Hazen-Williams.
56
3.7 Diseño redes matrices
Luego de haber corregido los coeficientes de la ecuación de Hazen-Williams, se procede

a diseñar las redes de estudio con base en estos nuevos valores, para verificar si en
realidad se logró minimizar el error en las pérdidas por fricción en las tuberías y en la
presión de los nudos. Para ello, se tuvieron en cuenta los siguientes parámetros de diseño
en común para estas redes:
• Presión mínima de diseño 20 m.c.a.

• Los materiales utilizados para el diseño de estas redes fueron CCP, GRP y acero.
Por ello, se adoptó la rugosidad absolutas, los coeficientes de Hazen-Williams
sugeridos en la normativa colombiana y los calculados a partir de la Figura 3.0-1.
• Se adoptó un valor para la viscosidad cinemática equivalente a una temperatura
de15°C.
• Aceleración de la gravedad de 9.8m/s2.
• El funcionamiento de las redes es por gravedad (no incluye bombas).
• No se tiene en cuenta el efecto de las pérdidas menores, teniendo en cuenta los
resultados obtenidos por Bombardelli y García (2003)
3.7.1 Procedimiento Diseño
Con base en los parámetros escogidos para el diseño de las redes matrices de estudio, a
continuación se describe el procedimiento que se realizó para llevar a cabo el diseño de
estas, a fin de obtener los diámetros de las tuberías y factores hidráulicos utilizando el
programa REDES, dado que el programa EPANET no permite diseñar RDAP.
• Definir la rugosidad absoluta y los coeficientes de Hazen-Williams para cada una

de las redes a analizar.
• Cargar la red en REDES de una forma similar a lo descrito en el Numeral 3.4
• Establecer la viscosidad cinemática del agua, con base en la temperatura del
agua. Para este trabajo se adoptó una viscosidad cinemática de 1.14x10-6 m2/s,
correspondiente a una temperatura de 15ºC.
57
• Suministrar al programa REDES la lista de diámetros, Igualmente, en este paso

es importante establecer la presión mínima y la ecuación de fricción que se va a
utilizar para el diseño de las redes, que para este caso es la ecuación de Hazen-
Williams con base en la Figura 3.7-1.
Figura 3.7-1 Interfaz programa REDES para definir lista de diámetros, curva de costos y
ecuación de fricción en el diseño de RDAP.
58
• Con base en la selección de la ecuación de fricción, se debe suministrar la

rugosidad absoluta (Darcy-Weisbach) o los coeficientes de Hazen-Williams para
los materiales de cada una de las tuberías, utilizando los valores descritos en el
Numeral 2.4
Figura 3.7-2 Interfaz programa REDES para definir rugosidad absoluta o coeficiente de Hazen-
Williams en las tuberías.
• Diseñar de la red matriz, utilizando alguna metodología que incluye el programa

REDES.
Figura 3.7-3 Interfaz programa REDES para realizar el diseño de RDAP.
59
• Obtener los factores hidráulicos calculados con el diseño de una RDAP, como
diámetro, número de Reynolds, factores de fricción, pérdidas por fricción y presión
en los nudos, para analizar el efecto de la corrección de los coeficientes de Hazen-
Williams.
Figura 3.7-4 Interfaz programa Redes para obtener los factores hidráulicos de la red.
• Luego de obtener los factores hidráulicos verificar si estos se encuentran dentro de

los rangos de número de Reynolds y velocidades establecidos para el nuevo
coeficiente de Hazen-Williams. Si se encuentran por fuera de ese rango, debe
aumentarse el diámetro en las tuberías que no cumplan este criterio.
• Posteriormente, se realiza de nuevo el procedimiento descrito anteriormente pero
utilizando los coeficientes de Hazen-Williams sugeridos en la normativa
colombiana, a fin de visualizar el efecto de los nuevos coeficientes de Hazen-
Williams sugeridos a partir de la Figura 3.7-1 en la disminución de los errores de
las pérdidas por fricción y distribución de presiones.
60
3.8 Validación diseños
Luego de haber diseñado las redes matrices de estudio con los coeficientes de Hazen-
Williams obtenidos a partir del Diagrama de Flujo de la Figura 3.7-1, y los coeficientes de
Hazen-Williams sugeridos en la normativa colombiana, se procede a validar estos diseños
comprobando los mismos con la ecuación de Darcy-Weisbach, ecuación físicamente
basada que representa el comportamiento hidráulico real de las redes matrices, llevando
a cabo el siguiente procedimiento:
• Después de diseñar una red matriz, se toman los diámetros obtenidos para cada
uno de las tuberías de la misma para un coeficiente de Hazen-Williams
correspondiente a cualquier material ya sea acero, GRP y CCP.
• Luego, tomando la misma topología de la red matriz de estudio, se ingresan los
diámetros obtenidos en el anterior paso y la rugosidad de las mismas por la
rugosidad absoluta equivalente al material que se está analizando, sugerido en la
normativa colombiana (Ver Tabla 3.3-1) .
• Con base en ello, se procede a calcular la hidráulica o realizar la comprobación de
diseño en el programa REDES utilizando la ecuación de Darcy-Weisbach,
siguiendo el mismo procedimiento descrito en el Numeral 3.4.
• Posteriormente, se obtiene la presión en los nudos y las pérdidas por fricción en
cada una de las tuberías de la red matriz.
• Después, se realiza la comprobación de diseño utilizando la ecuación de Hazen-
Williams, utilizando el programa EPANET y se comparan los resultados de esta
con los obtenidos en el anterior paso, para establecer el error en la distribución de
presiones, errores en las pérdidas por fricción y cumplimiento de la presión minima
establecida en el diseño de las redes de estudio.
• Si existen diferencias significativas entre estas se analizan con base en el
diagrama de flujo propuesto en la Figura 3.7-1 y en el diagrama de Moody.
61
Capítulo 4. RESULTADOS
En el presente capítulo se presentan los resultados obtenidos en relación con los efectos
del uso de la ecuación de Hazen-Williams vs la Ecuación de Darcy-Weisbach en redes
matrices. Para ello, este capítulo se dividirá en cuatro partes. En la primera, se analizará
el efecto del uso de la ecuación de Hazen-Williams cuando se utiliza un coeficiente de
rugosidad constante (sugerido en la normativa colombiana) en las redes de estudio, a fin
de establecer el efecto de esta suposición en la estimación de las pérdidas por fricción y
distribución de presiones de estas redes. Luego, con base en los hallazgos encontrados
en esta primera parte, en la segunda se ilustrará el programa que se desarrolló para
corregir los coeficientes de Hazen-Williams sugeridos en la normativa colombiana,
procurando minimizar los errores en las pérdidas por fricción hasta en un 5%.
Posteriormente, en la tercera parte se mostrarán los resultados obtenidos en relación con
el diseño de las redes matrices de estudio, utilizando tanto los coeficientes de Hazen-
Williams corregidos y sugeridos en la normativa colombiana. En la última parte, se
expondrá la comparación de estos diseños y su respectiva validación, utilizando la
ecuación de Darcy-Weisbach, para analizar el efecto del uso de los coeficientes de
Hazen-Williams corregidos en la disminución de los errores en las pérdidas por fricción y
distribución de presiones.
4.
4.1 Efecto del uso de la ecuación de Hazen-Williams en las

redes de estudio
En esta parte del trabajo se mostrarán los resultados obtenidos en relación con la
variabilidad del coeficiente de Hazen-Williams en las redes estudio y su efecto en la
estimación de las pérdidas por fricción y distribución de presiones, cuando se utiliza un
coeficiente de rugosidad constante, el cual es sugerido por prescripción en la normativa
colombiana. Estos resultados se obtuvieron siguiendo los pasos 2, 3, 4 y 5 de la
metodología propuesta en el Numeral 3.1.
62
• Variación Coeficiente de Hazen-Williams
En relación con los resultados obtenidos en la variación del coeficiente de

Hazen-Williams en la red matriz de Medellín, para materiales como el CCP y Acero, se
puede ver en las Gráficas 4.1-1 y 4.1-3, que los valores sugeridos en la normativa
colombiana subestiman la magnitud de éstos para cualquier diámetro y rango de números
de Reynolds. Sin embargo, para el material GRP (Ver Gráfica 4.1-2), se puede apreciar
en la Gráfica 4.1-2, que el coeficiente de Hazen-Williams sugerido para este material
representa, para la mayoría de tuberías de gran diámetro, la variabilidad que presenta
este coeficiente para un rango de Números de Reynolds entre 40000 y 2.000.000.
Adicionalmente, es importante mencionar que independiente del tipo de material, los
coeficientes de Hazen-Williams sugeridos en las normativa colombiana, no pueden
representar la variabilidad de este en el flujo laminar, dado que en este tipo de flujo no es
posible predecir la magnitud del coeficiente de Hazen-Williams, como se puede ver en las
Gráficas 4.1-1, 4.1-2 y 4.1-3.
Variación Coeficiente de Hazen-Williams (CCP, ks:0.12 mm, CHw:120)

160
150
d=2.5 m
140
d=2.0 m
CHw
130 d=1.5 m
d=1.0 m
120 CHw sugerido
CHw Teorico Red

110
100
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-1 Variación del Coeficiente de Hazen-Williams en la red matriz de Medellín para el CCP.
63
Variación Coeficiente de Hazen-Williams (GRP, ks:0.03 mm, CHw:150)

160
150
140
130 d=2.5 m
d=2.0 m
CHw
120
d=1.5 m
110 d=1.0 m
100 CHw sugerido
CHw Teorico red

90
80
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-2 Variación del Coeficiente de Hazen-Williams en la red matriz de Medellín para el GRP.
Variación Coeficiente de Hazen-Williams (Acero, ks:0.45 mm, CHw: 120)

160
150
140
d=2.5 m
130
d=2.0 m
CHw
120
d=1.5 m
110 d=1.0 m
CHw sugerido
100
CHw Teorico Red
90
80
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-3 Variación del Coeficiente de Hazen-Williams en la red matriz de Medellín para el Acero.
64
• Errores en el cálculo de las pérdidas por fricción
Con base en los resultados obtenidos en relación con la variabilidad del coeficiente de
Hazen-Williams para los materiales CCP, GRP y Acero en la red matriz de Medellín,
luego se procedió a estimar los errores en la estimación de las pérdidas por fricción,
siguiendo el procedimiento descrito en el Numeral 3.5. A partir de lo anterior, se obtuvo
que los valores sugeridos en las normativa colombiana para el Acero y CCP, los cuales
subestiman la magnitud del coeficiente de Hazen-Williams real, sobreestiman las
pérdidas por fricción hasta en un 40% y 20% respectivamente (ver Gráficas 4.1-4 y 4.1-6).
Adicionalmente, para estos materiales se puede apreciar en la Gráficas 4.1-4 y 4.1-6 que
los coeficientes de Hazen-Williams sugeridos en la normativa colombiana son válidos para
un rango de número de Reynolds superior a 200.000 y 1.000.000 para el Acero y CCP
respectivamente, dado que en estos rangos los errores en las pérdidas por fricción en la
red matriz de Medellín no superan el 10% en la subestimación o sobreestimación de las
pérdidas por fricción. Lo anterior, se refleja en el comportamiento que presenta en
general el error en las pérdidas por fricción en función del diámetro y el número de
Reynolds en las Gráficas 4.1-4 y 4.1-6, donde se llega a un valor máximo y luego este
disminuye a medida que aumenta el número de Reynolds.
Por otro lado, en relación con el efecto de la variabilidad del coeficiente de Hazen-
Williams en la red matriz de Medellín cuando se utiliza el material GRP, se puede ver en
la Gráfica 4.1-5 que el valor sugerido en la normativa colombiana es adecuado dado que
los errores en las pérdidas por fricción no superan el 10% en la subestimación en las
pérdidas por fricción. Pese a ello, se corre el riesgo de que para números de Reynolds
inferiores a 10000 o superiores a 1.000.000 se subestimen las pérdidas por fricción en la
red, de tal manera que las presiones en ésta sean inferiores a las calculadas a partir de la
ecuación de Hazen-Williams y, por lo tanto, no se pueda cumplir con las presiones
mínimas exigidas en la red de estudio.
65
Errores pérdidas por fricción (CCP, ks: 0.12 mm, CHw:120 )

50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10%
d=2.0 m
0%
-10% d=1.5 m
-20% d=1.0 m
-30% Errores Red

-40%
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-4 Errores en las pérdidas por fricción en la red matriz de Medellín, utilizando la ecuación de
Darcy-Weisbach y Hazen-Williams para el CCP.
Errores pérdidas por Fricción (GRP, ks:0.03 mm, CHw:150 )

20%
10%
0%
(Shw-Sdw)/Sdw
d=2.5 m
-10%
d=2.0 m
-20%
d=1.5 m
-30% d=1.0 m
Error Red
-40%
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Darcy-Weisbach y Hazen-Williams para el GRP.
66
(Shw-Sdw)/Sdw vs. Número de Reynolds (Acero, ks:0.45 mm, CHw:120)

50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10%
d=2.0 m
0% d=1.5 m
-10% d=1.0 m
-20% Series5
-30% Error Red
-40%
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Darcy-Weisbach y Hazen-Williams para el Acero.
• Errores en las presiones de los nudos
Con base en los resultados obtenidos en relación con el efecto de la variabilidad del
coeficiente de Hazen-Williams en los errores en las presiones de los nudos en la red
matriz de Medellín, se puede ver en las Gráficas 4.1-7, 4.1-8 y 4.1-9 que los valores
sugeridos en las normas para el CCP, GRP y Acero sobreestiman las presiones en los
nudos hasta en un 10%, 17% y 70% respectivamente, cerca al embalse de Ayurá, el cual
suministra el mayor caudal a esta red (5 m3/s). Lo anterior, es consistente con los
resultados obtenidos en cuanto a los errores en las pérdidas por fricción para estos
materiales, dado que para este sector de la red matriz de Medellín se presentan números
de Reynolds superiores a 9.000.000, rango en el cual los coeficientes de Hazen-Williams
de estos materiales, tiende a subestimar las pérdidas por fricción como se puede ver en
las Gráficas 4.1-4, 4.1-5 y 4.1-6. Adicionalmente, en las Gráficas 4.1-7 y 4.1-9 se puede
observar que para el CCP y Acero, respectivamente, en algunos nudos de esta red se
presenta subestimación de las presiones, dado que en algunas tuberías de la red matriz
de Medellín se sobreestima las pérdidas por fricción para los coeficientes de rugosidad
sugeridos en la normativa colombiana.
67
Manantiales
Ayurá
Gráfica 4.1-7 Errores en las presiones en los nudos a partir de los resultados obtenidos en la red
matriz de Medellín por medio de la ecuación de Darcy-Weisbach y Hazen-Williams para el CCP.
Ayurá
matriz de Medellín por medio de la ecuación de Darcy-Weisbach y Hazen-Williams para el GRP.
68
Ayurá
matriz de Medellín por medio de la ecuación de Darcy-Weisbach y Hazen-Williams para el Acero.
4.1.2 Red Hanoi
• Variación Coeficiente de Hazen-Williams
En relación con los resultados obtenidos en la variación del coeficiente de

Hazen-Williams en la red Hanoi, y específicamente para el material CCP, se puede ver en
la Gráfica 4.1-10 que el coeficiente de Hazen-Williams sugerido en la normativa
colombiana (CHw: 120) tiende a subestimar la magnitud real de este para cualquier
diámetro y rango de números de Reynolds. Sin embargo, los valores sugeridos en dicha
normativa para los materiales GRP y Acero (Ver Gráficas 4.1-12 y 4.1-13) tienden a
sobreestimar la magnitud real del coeficiente de Hazen-Williams, para Números de
Reynolds superiores a 100.000 y 500.000, respectivamente.
69
Variación Coeficiente de Hazen-Williams (CCP, ks:0.12 mm, CHw:120)

160
150
140
130
d=1.5 m
CHw
120 d=1.0 m
CHw sugerido
110
Chw Teorico Red
100
90
80
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-10 Variación del Coeficiente de Hazen-Williams en la red Hanoi para el CCP.
Variación Coeficiente de Hazen-Williams (GRP, ks:0.03 mm, CHw:150)

160
150
140
130
d=1.5 m
CHw
120 d=1.0 m
110 CHw sugerido
100
90
80
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-11 Variación del Coeficiente de Hazen-Williams en la red Hanoi para el GRP.
70
Variación Coeficiente de Hazen-Williams (Acero, ks: 0.45 mm, CHw: 120)

160
150
140
130
d=1.5 m
CHw
120 d=1.0 m
CHw sugerido
110
Chw teorico red
100
90
80
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-12 Variación del Coeficiente de Hazen-Williams en la red Hanoi para el Acero.
• Errores en el cálculo de las pérdidas por fricción
Con base en los resultados obtenidos en relación con la variabilidad del coeficiente de
Hazen-Williams para los materiales CCP, GRP y Acero en la red Hanoi, se procedió a
estimar los errores en la estimación de las pérdidas por fricción, siguiendo el
procedimiento descrito en el Numeral 3.5. A partir de lo anterior, se obtuvo que el valor
sugerido en las normativa colombiana para el material CCP, el cual subestiman la
magnitud del coeficiente de Hazen-Williams real, sobreestiman las pérdidas por fricción
en un rango de 10% y 30% (ver Gráfica 4.1-14). Por otro lado, en relación con el efecto de
la variabilidad del coeficiente de Hazen-Williams en la red Hanoi para los materiales GRP
y Acero, se puede ver en la Gráficas 4.1-15 y 4.1-16 que el valor sugerido en la normativa
colombiana es adecuado para estos, dado que los errores en las pérdidas por fricción en
la red Hanoi no superan el 10% en la subestimación o sobreestimación en las pérdidas
por fricción. Pese a ello, se corre el riesgo de que para números de Reynolds superiores a
1.000.000 se subestimen las pérdidas por fricción en la red, de tal manera que las
presiones en esta sean inferiores a las calculadas a partir de la ecuación de Hazen-
Williams y por lo tanto no se pueda cumplir con las presiones mínimas exigidas en la red
de estudio.
71
Errores pérdidas por fricción (CCP, ks: 0.12 mm, CHw:120 )

50%
40%
30%
20%
(Shw-Sdw)/Sdw
10%
d=1.5 m
0%
d=1.0 m
-10%
Errores red
-20%
-30%
-40%
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-13 Errores en las pérdidas por fricción en la red Hanoi, utilizando la ecuación de Darcy-
Weisbach y Hazen-Williams para el CCP.
Errores pérdidas por fricción (GRP, ks:0.03 mm, CHw:150 )

50%
40%
30%
20%
(Shw-Sdw)/Sdw
10%
d=1.5 m
0%
d=1.0 m
-10%
Errores red
-20%
-30%
-40%
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Weisbach y Hazen-Williams para el GRP.
72
Errores pérdidas por fricción (Acero, ks:0.45 mm, CHw: 120 )

50%
40%
30%
(Shw-Sdw)/Sdw
20%
10% d=1.5 m
0% d=1.0 m
-10% Errores red
-20%
-30%
-40%
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Weisbach y Hazen-Williams para el Acero.
• Errores en las presiones de los nudos
Con base en los resultados obtenidos en relación con el efecto de la variabilidad del
coeficiente de Hazen-Williams en los errores de la distribución de presiones en los nudos
de la red Hanoi, se puede ver en la Gráfica 4.1-17 que el valor sugerido en las normas
colombianas para el CCP subestima las presiones en todos los nudos de la red, hasta en
un 60%. Lo anterior, es consistente con los resultados obtenidos en cuanto a los errores
en las pérdidas por fricción para este material, dado que en la red Hanoi se presentan
números de Reynolds entre 100.000 y 6.000.000, rango en el cual el coeficiente de
Hazen-Williams de este material tiende a sobrestimar las pérdidas por fricción como se
puede ver en la Gráfica 4.1-14. En contraste, para el material GRP y Acero, el coeficiente
de Hazen-Williams sugerido para estos materiales sobreestima las presiones en todos los
nudos de la red (Ver Gráficas 4.1-15 y 4.1-16), lo cual es producto de la subestimación de
las pérdidas por fricción cuando se utiliza la ecuación de Hazen-Williams para estos
materiales (Ver Gráficas 4.1-15 y 4.1-16). Es importante aclarar que el coeficiente de
Hazen-Williams sugerido para el GRP según los errores en la distribución de presiones es
adecuado dado que estos no superan el 10%. Sin embargo, el coeficiente de Hazen-
Williams para el Acero no es adecuado dado que la distribución de presiones es
altamente sensible si en alguna tubería de la red es el número de Reynolds es superior a
5.000.000, rango donde se puede subestimar en exceso las pérdidas por fricción.
73
Hanoi por medio de la ecuación de Darcy-Weisbach y Hazen-Williams para el CCP.
Hanoi por medio de la ecuación de Darcy-Weisbach y Hazen-Williams para el GRP.
74
Hanoi por medio de la ecuación de Darcy-Weisbach y Hazen-Williams para el Acero.
Para este trabajo se analizará la red matriz estudiada por Bombardelli y García (2003)
para establecer los efectos hidráulicos que puede tener el uso inadecuado de la ecuación
de Hazen-Williams en el diseño de redes matrices. Para ello, se tomaron los cálculos
realizados por estos investigadores, en cuanto a la rugosidad absoluta, en los tres tramos
de la tubería principal de esta red (Ver Tabla 2.1-8). Luego, se efectuó el procedimiento
descrito en el Numeral 3.5 para cada uno de estos, teniendo en cuenta el coeficiente de
Hazen-Williams utilizado en el diseño original para esta tubería (Concreto, CHw: 120).
75
• Tramo 1 (ks:0.009 m, Longitud: 3959 m)
Con respecto a los resultados obtenidos en cuanto a la variación del coeficiente de

Hazen-Williams en el primer tramo de la tubería principal de la red matriz analizada
por Bombardelli y Garcia (2003), se puede apreciar en la Gráfica 4.1-19 que el
coeficiente de Hazen-Williams utilizado para el diseño original de esta red (CHw: 120)
sobreestima el comportamiento real que presenta dicho coeficiente, para las
condiciones de flujo que presenta el tramo, donde el número de Reynolds es igual
3.8x106 y el coeficiente de Hazen-Williams a 86.
Variación Coeficiente de Hazen-Williams (CCP, ks:0.009 m, CHw:120)

125
120
115 d=2.5 m
110
d=2.0 m
105
CHw
d=1.5 m
100
95 CHw teorico
90
86 d=2.286 (Bombardelli-
85 García)
80
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-19 Variación del Coeficiente de Hazen-Williams en el primer tramo de la tubería principal
de la red matriz analizada por Bombardelli y García (2003).
Teniendo en cuenta lo anterior, el coeficiente de Hazen-Williams real en este tramo es

inferior al utilizado en el diseño original, lo cual implicaría que este último subestime las
pérdidas por fricción ya que supone una tubería más lisa. Esto se refleja en la
Gráfica 4.1-20, donde se presenta un error de 46% en la subestimación de las pérdidas
por fricción, para un número de Reynolds igual a 3.8x106. Este hecho representa un
riesgo para la red matriz analizada por Bombardelli y García (2003), dado que desde el
primer tramo podrían presentarse problemas en la operación de la misma, en cuanto a las
presiones.
76
Errores pérdidas por Fricción (CCP, ks: 0.009 m, CHw:120 )

50%
40%
30%
d=2.5 m
20%
(Shw-Sdw)/Sdw
10%
d=2.0 m
0%
-10% d=1.5 m
-20%
v=0.5 m/s d=2.286 (Bombardelli-
-30%
v= 1 m/s García)
-40%
-46.118% v=2 m/s
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-20 Error en las pérdidas por fricción en el primer tramo de la tubería principal de la red
matriz analizada por Bombardelli y García (2003).
• Tramo 2 (ks: 0.004 m, Longitud: 5179 m)
Con respecto a los resultados obtenidos en cuanto a la variación del coeficiente de

Hazen-Williams en el segundo tramo de la tubería principal de la red matriz analizada por
Bombardelli y Garcia (2003), se puede apreciar en la Gráfica 4.1-21 que el coeficiente de
Hazen-Williams utilizado para el diseño original de esta red (CHw: 120) sobreestima el
comportamiento real que presenta dicho coeficiente, para las condiciones de flujo que
presenta dicho tramo, donde el número de Reynolds es igual 3.8x106 y el coeficiente de
Hazen-Williams a 97. Por lo anterior, el coeficiente de Hazen-Williams real en este tramo
es inferior al utilizado en el diseño original, lo cual implicaría que este último subestime las
pérdidas por fricción ya que supone una tubería más lisa. Esto se refleja en la Gráfica 4.1-
22, donde se presenta un error de 33% en la subestimación de las pérdidas por fricción,
para un número de Reynolds igual a 3.8x106. Este hecho representa un riesgo para la red
matriz analizada por Bombardelli y García (2003), dado que se acumularía el efecto
presentado en el primer tramo de la tubería principal de esta red, lo cual podría generar
problemas en cuanto a la operación de la misma.
77

140
130 d=2.5 m
120 d=2.0 m
CHw
110 d=1.5 m
100 CHw teorico

97
90 d=2.286
(Bombardelli-
García)
80
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-21 Variación del Coeficiente de Hazen-Williams en el segundo tramo de la tubería principal
de la red matriz analizada por Bombardelli y García (2003).

50%
40%
30%
d=2.5 m
20%
(Shw-Sdw)/Sdw
10% d=2.0 m
0%
-10% v=0.5 m/s d=1.5 m
-20% v= 1 m/s
-30% v=2 m/s d=2.286 (Bombardelli-

-32.687% García)
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-22 Error pérdidas por fricción en el segundo tramo de la tubería principal de la red matriz
analizada por Bombardelli y García (2003).
78
• Tramo 3 (ks: 0.007m, Longitud 4560 m)
En relación con los resultados obtenidos en cuanto a la variación del coeficiente de

Hazen-Williams en el tercer tramo de la tubería principal de la red matriz analizada por
Bombardelli y Garcia (2003), se puede apreciar en la Gráfica 4.1-23 que el coeficiente
de Hazen-Williams utilizado para el diseño original de esta red (CHw: 120)
sobreestima el comportamiento real que presenta dicho coeficiente, para las
condiciones de flujo que presenta dicho tramo, donde el número de Reynolds es igual
3.8x106 y el coeficiente de Hazen-Williams a 89.41. Por lo anterior, el coeficiente de
Hazen-Williams real en este tramo es inferior al utilizado en el diseño original, lo cual
implicaría que este último subestime las pérdidas por fricción ya que supone una
tubería más lisa. Esto se refleja en la Gráfica 4.1-24, donde se presenta un error de
42% en la subestimación de las pérdidas por fricción, para un número de Reynolds
igual a 3.8x106. Este hecho representa un riesgo para la red matriz analizada por
Bombardelli y García (2003), dado que se acumularía el efecto presentado en el
primer y segundo tramo de la tubería principal de esta red, lo cual podría generar
problemas aun mayores en cuanto a la operación de la misma.

130
125
d=2.5 m
120
115 d=2.0 m
110
CHw
d=1.5 m
105
100
CHw teorico
95
90 d=2.286
89,41 (Bombardelli-
85 García)
80
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-23 Variación del Coeficiente de Hazen-Williams en el tercer tramo de la tubería principal de
la red matriz analizada por Bombardelli y García (2003).
79

50%
40%
30%
20%
(Shw-Sdw)/Sdw
10% d=2.5 m
0%
d=2.0 m
-10%
v=0.5 m/s d=1.5 m
-20%
-30% v= 1 m/s
d=2.286 (Bombardelli-
-40% v=2 m/s García)
-42.107%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.1-24 Error en las pérdidas por fricción en el tercer tramo de la tubería principal de la red
matriz analizada por Bombardelli y García (2003).
4.2 Corrección coeficientes de Hazen-Williams
Teniendo en cuenta que en las redes de estudio, los coeficientes de Hazen-Williams

sugeridos por la normativa colombiana para los materiales CCP y Acero no son
adecuados con base en el efecto que tiene la variación de los coeficientes de Hazen-
Williams para estos materiales, según la ecuación de Darcy-Weisbach, en la distribución
de presiones, se llevó a cabo la metodología propuesta en el Numeral 3.6 de este trabajo
para corregir dichos coeficientes. Es importante mencionar que igualmente se corrigió el
coeficiente de Hazen-Williams para el material GRP, a fin de evitar cualquier
subestimación en las pérdidas por fricción en el diseño de redes matrices, de acuerdo con
el comportamiento visto en las Gráficas 4.1-5 y 4.1-14.
Para corregir dichos coeficientes, primero se desarrolló un programa con macros en

Excel, utilizando el lenguaje programación Visual Basic. Luego, para cada uno de los
materiales de estudio, se establecieron los valores de los coeficientes de Hazen-Williams,
con sus respectivos límites de diámetros, velocidades y rangos de números de Reynolds.
80
4.2.1 Programa Corrección Coeficientes de Hazen-Williams
Para desarrollar este programa, se tuvo en cuenta la metodología propuesta en el

Numeral 3.6 de este trabajo para corregir los coeficientes de Hazen-Williams para cada
uno de los materiales seleccionados. Con base en ello, se elaboró el código fuente del
mismo, el cual está compuesto por los siguientes módulos:
• Módulo 1: Macros que incluyen funciones básicas para calcular el factor de

fricción de la ecuación de Hazen-Williams y la ecuación de Darcy-Weisbach, y
Algoritmo de la metodología propuesta en el Numeral 3.6 de este trabajo (Ver
Anexo 7).
• Módulo 2: Macros para manejar la interfaz del programa (Ver Anexo 7)
Posteriormente, se procedió a desarrollar la interfaz de este programa a fin de que el

usuario pueda visualizar los resultados obtenidos en cuanto a la corrección de los
coeficientes de Hazen-Williams, así como el efecto de los mismos en la estimación de las
pérdidas por fricción en el diagrama de Moody. Para ello, la hoja de Excel incluye cuatro
hojas las cuales son: Programa, Resultados, Error pérdidas por fricción y Diagrama de
Moody (Ver Anexos 3, 4, 5 y 6). El procedimiento que debe seguir el usuario del
programa para obtener dichos resultados se ilustra a continuación:
• Primero el usuario debe ingresar el diámetro mínimo y máximo de la lista de

diámetros disponible para el diseño de redes matrices en la hoja electrónica
denominada “Programa” (Ver Figura 4.2-1).
• Segundo, se debe ingresar la rugosidad absoluta del material, en metros, al cual
se desea corregir el coeficiente de Hazen-Williams en la hoja electrónica
denominada “Programa” (Ver Figura 4.2-1).
• Tercero, el usuario debe ingresar la viscosidad cinemática del agua,
correspondiente a la temperatura del agua en la hoja mencionada anteriormente
(Ver Figura 4.2-1).
81
• Posteriormente, el usuario puede ingresar el error admisible en las pérdidas por

fricción en la hoja “Programa” (Ver Figura 4.2-1). Para mayor eficiencia de la
metodología propuesta de esta tesis se sugiere utilizar 5%.
• Luego, el usuario debe hacer clic sobre el botón “calcular coeficiente de Hazen-
Williams”, en la hoja “Programa”, para obtener el coeficiente de Hazen-Williams
apropiado para el material de estudio (Ver Figura 4.2-1).
Figura 4.2-1 Interfaz programa para ingresar datos de entrada.
• Posteriormente, si el diámetro mínimo de la lista de diámetros cumple con el error

admisible en las pérdidas por fricción, en la hoja “Programa”, se puede visualizar el
coeficiente de Hazen-Williams corregido para el material de estudio, el rango de
números de Reynolds válido para su aplicabilidad, así como la velocidad límite que
puede presentar las tuberías en la red matriz para cumplir dicho error (Ver Figura
4.2-2). En caso contrario el mismo programa advierte que el diámetro mínimo
debe ser aumentado tal como se ilustra en la Figura 4.2-3.
82
• Luego, el usuario tiene la opción de visualizar los resultados en la hoja

denominada con su mismo nombre. Igualmente, haciendo clic sobre los botones
“Error pérdidas por fricción” y “Diagrama de Moody”, el usuario puede visualizar el
efecto del uso de los coeficientes de Hazen-Williams en la disminución de los
errores en la estimación de las pérdidas por fricción (Ver Anexos).
Figura 4.2-2 Interfaz programa para visualizar resultados del coeficiente de Hazen-Williams corregido.
Figura 4.2-3 Mensaje de Advertencia del programa cuando el diámetro mínimo no cumple con el error
admisible en las pérdidas por fricción.
83
4.2.2 Coeficientes de Hazen-Williams Corregidos
Con base en el programa desarrollado y teniendo en cuenta que el diámetro máximo

considerado es igual a 2500 mm, así como los valores sugeridos de la rugosidad absoluta
de los materiales de estudio (Ver Numeral 3.3), se obtuvieron los siguientes coeficientes
de Hazen-Williams, para un rango especifico de números de Reynolds, diámetro mínimo y
una velocidad límite, a fin de evitar la subestimación de las pérdidas por fricción, caso
crítico en el diseño de redes matrices para cumplir con las presiones de operación de las
mismas:
Número de Número de Diámetro Velocidad

CHw
Material ks (mm) Reynolds Reynolds mínimo límite
corregido
mínimo máximo (mm) (m/s)
CCP 0.12 144 4.50E+04 4.15E+05 500 0.95
GRP 0.03 147 7.00E+04 2.48E+06 500 5.6
Acero 0.45 139 2.50E+04 2.55E+05 1000 0.3
Tabla 4.2-1 Coeficientes de Hazen-Williams corregidos para los materiales CCP, GRP y Acero.
A partir de la Tabla 4.2-1 y la Gráfica 4.2-1, se puede apreciar que para el material CCP
se obtuvo un coeficiente de Hazen-Williams de 144, el cual permite, en el diseño de
redes matrices, que las tuberías con diámetros de 2.5 metros no sobreestimen las
pérdidas por fricción en un 5% y, para tuberías de diámetro igual a 500 mm, no se
subestimen en la misma magnitud. Adicionalmente, se puede observar que este valor
sugerido permite manejar velocidades de hasta 0.95 m/s para minimizar en lo posible el
error en las pérdidas por fricción. Así mismo, se puede observar en la Gráfica 4.2-2, que
para este coeficiente presenta buena exactitud, con respecto a la ecuación de Darcy-
Weisbach, cuando los factores de fricción de este se ubican en la zona de FTHL
delimitada por la ecuación deducida por el autor (Flechas, 2010).
84
De otro lado, se puede apreciar en la Tabla 4.2-1 y la Gráfica 4.2-5, que el coeficiente de
Hazen-Williams adecuado para el Acero es igual a 139, el cual permite controlar la
sobrestimación o subestimación de las pérdidas por fricción hasta en un 5%, para
velocidades inferiores a 0.3 m/s y específicamente para diámetro superiores a 1000 mm.
Además, en relación con el rango de número de Reynolds establecido para dicho
coeficiente, se puede ver en la Gráfica 4.2-6, que en este presenta buena exactitud, con
respecto a la ecuación de Darcy-Weisbach, cuando los factores de fricción de este se
ubican en la zona de FTHL definida por el autor (Flechas, 2010).
Por último, con respecto a la corrección del coeficiente de Hazen-Williams para el

material GRP (Ver Tabla 4.2-1 y Gráfica 4.2-3), se obtuvo que el valor de este debe ser
igual a 147, a fin de evitar, en lo máximo posible, la subestimación en las pérdidas por
fricción en tuberías en un 5%, específicamente para diámetros iguales o superiores a
500 mm. Adicionalmente, a diferencia de los materiales CCP y Acero, el coeficiente de
Hazen-Williams sugerido para este material permite manejar velocidades de hasta 5.6 m/s
en las tuberías de una red matriz. Así mismo, este coeficiente permite diseñar tuberías en
un rango de número de Reynolds más amplio en comparación a los obtenidos para el
Acero y CCP. Sin embargo, se encuentra una similitud con estos materiales, dado que la
buena exactitud de este coeficiente, con respecto a la ecuación de Darcy-Weisbach, se
logra cuando los factores de fricción de éste se ubican en la zona del FTHL (Ver Gráfica
4.2-4), definida por el autor (2010).
85
Corrección Coeficiente de Hazen-Williams (CCP, ks: 0.12 mm, CHw: 144)

50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10%
d=2.0 m
0%
d=1.5 m
-10% v=0.5 m/s
v= 1 m/s d=1.0 m
-20% v=2 m/s
v= 3 m/s d=0.5 m
-30% v= 5m/s
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.2-1 Corrección Coeficiente de Hazen-Williams para el CCP.
Diagrama de Moody (CCP, ks: 0.12 mm, CHw: 144)

0.1
Factor de Fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Factor de Fricción (Darcy-Weisbach) Factor de fricción (Hazen-Williams)
Número de Reynolds Mínimo Número de Reynolds Máximo
Límite FTHR (Colebrook-White) Límite FTHR (Prantdl-Von Kármán)
Límite FTHL (Colebrook-White) Límite FTHL (Prantdl-Von Kármán)
Límite FTHL (Blasius) Flujo laminar
Gráfica 4.2-2 Diagrama de Moody para el Coeficiente de Hazen-Williams corregido del CCP.
86
Corrección Coeficiente de Hazen-Williams (GRP, ks:0.03 mm, CHw:147 )

50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10%
d=2.0 m
0%
v=2 m/s d=1.5 m
-10% v=3 m/s
v= 5m/s d=1.0 m
-20% d=0.5 m
-30%
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.2-3 Corrección Coeficiente de Hazen-Williams para el GRP.
Diagrama de Moody (GRP, ks: 0.03 mm, CHw: 147)

0.1
Factor de Fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.2-4 Diagrama de Moody para el Coeficiente de Hazen-Williams corregido del GRP.
87
Corrección Coeficiente de Hazen-Williams (Acero, ks: 0.45 mm,CHw: 139)

50%
40%
30%
20%
(Shw-Sdw)/Sdw
10% d=2.5 m
0% d=2.0 m
-10% d=1.5 m
-20% d=1.0 m
v=0.5 m/s
-30% v= 1 m/s
v=2 m/s
-40% v=3 m/s
v=5 m/s
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.2-5 Corrección Coeficiente de Hazen-Williams para el Acero.
Diagrama de Moody (Acero, ks: 0.45 mm, CHw: 139)

0.1
Factor de Fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.2-6 Diagrama de Moody para el Coeficiente de Hazen-Williams corregido del Acero.
88
4.3 Diseño redes matrices
Teniendo en cuenta los coeficientes de Hazen-Williams corregidos para los materiales

CCP, GRP y Acero, se procede a evaluar estos coeficientes diseñando las redes de
estudio con estos, para luego comparar estos resultados con los diseños obtenidos a
partir de los coeficientes de Hazen-Williams sugeridos en la normativa colombiana. Para
ello, se tuvieron en cuenta las siguientes listas de diámetros, para llevar a cabo la
metodología propuesta en el Numeral 3.7.1:
CHw sugeridos CCP (CHw: 144) y Acero

Norma GRP (CHw: 147) (CHw: 139)
Diámetro (mm) Diámetro (mm) Diámetro (mm)

300 - -
350 - -
400 - -
450 - -
500 500 -
600 600 -
750 750 -
800 800 -
1000 1000 1000
1200 1200 1200
1400 1400 1400
1500 1500 1500
1800 1800 1800
2000 2000 2000
2500 2500 2500
Tabla 4.3-1 Lista de diámetro utilizados en el diseño de las redes matrices.
89
Con base en los resultados obtenidos en el Numeral 4.1.1, se puede ver que el sector
donde se presenta mayores errores en la estimación de la distribución de presiones en la
red matriz de Medellín, cuando se utilizan los coeficientes de Hazen-Williams sugeridos
en la normativa colombiana, es el sector delimitado por el Embalse de Ayurá, debido al
caudal que aporta a esta red. Teniendo en cuenta este hecho, para este trabajo se decidió
diseñar solamente el sector comprendido por este embalse, tal como se ilustra en la
Figura 4.3-1
Figura 4.3-1 Sector Red matriz de Medellín delimitado para el Embalse Ayurá.
Es importante aclarar, que para este diseño se mantuvo la misma oferta de caudal del
embalse de Ayurá (5 m3/s). Adicionalmente, para este embalse se modificó la línea del
gradiente hidráulico a 1780 msnm, dado que los nodos cercanos a este presentaban una
línea de gradiente hidráulico similar a 1760 msnm, limitación que no permite utilizar la
metodología de diseño SOGH para cumplir con las presiones mínimas de diseño.
90
4.3.1.1 CCP
Con respecto al diseño de la red matriz de Medellín, en el sector delimitado por el

Embalse Ayurá, en la Gráfica 4.3-1 se puede apreciar que los diámetros obtenidos para el
coeficiente de Hazen-Williams corregido para este material, en su mayoría, sobreestiman
la magnitud de los diámetros cuando se compara con los obtenidos a partir del coeficiente
de Hazen-Williams sugerido en la normativa colombiana, debido a la restricción del
diámetro mínimo que tiene este para su correcta aplicación. Sin embargo, teniendo en
cuenta que este coeficiente es más liso que el recomendado en dicha normativa, en
algunos tubos se obtienen diámetros inferiores con respecto a los diámetros de este
último coeficiente.
Lo anterior, se ve reflejado en la Gráfica 4.3-2, donde se muestra las presiones obtenidas

después de realizar la comprobación del diseño correspondiente al coeficiente de Hazen-
Williams corregido, utilizando esta misma ecuación empírica, donde en algunos nudos de
la red se sobrestima o subestima las presiones con respecto a las obtenidas para el
coeficiente sugerido en la normativa colombiana, debido a la sobreestimación o
subestimación de los diámetros de esta red, respectivamente. Es importante aclarar que
en los dos diseños obtenidos para el material CCP, aparentemente la presión en los
nudos es igual o superior a la presión mínima establecida en este trabajo (20 m.c.a),
cuando se realiza la comprobación de de estos diseños utilizando la ecuación de Hazen-
Williams.
Diámetros obtenidos Diseño Red Medellín (CCP)

1400
1200
Diámetro (mm)
1000
800
600
CHw: 144
400
200 CHw:120
0
0 200 400 600 800
Número tubo
Gráfica 4.3-1 Diámetros obtenidos en el diseño de la red matriz de Medellín para el material CCP.
91
Presiones obtenidas Diseño Red Medellín (CCP)

260
240
Presion CHw: 144 (m.c.a)
220
200
180
160
140
120
100
80
60
40
20
0
0 20 40 60 80 100 120 140 160 180 200 220 240 260
Presión CHw: 120 (m.c.a)
Gráfica 4.3-2 Diagrama de dispersión presiones obtenidas en el diseño de la red matriz de Medellín
para el material CCP.
4.3.1.2 GRP
En relación con el diseño de la red matriz de Medellín utilizando el material GRP,

específicamente en el sector delimitado por el Embalse Ayurá, en la Gráfica 4.3-3 se
puede observar que los diámetros obtenidos para el coeficiente de Hazen-Williams
corregido para este material, en su mayoría, sobreestiman la magnitud de los diámetros
cuando se compara con los obtenidos a partir del coeficiente de Hazen-Williams sugerido
en la normativa colombiana, debido al aumento de la rugosidad de las tubería y la
restricción del diámetro mínimo que tiene este para su correcta aplicación. Sin embargo
existen algunos tuberías de esta red se subestiman los diámetros de la red.
Esto se ve reflejado en la Gráfica 4.3-4, donde se muestra las presiones obtenidas

Williams corregido para este material, utilizando esta misma ecuación empírica, donde en
la mayoría de los nudos se subestima o en su defecto se igualan cuando se comparan
con las obtenidas a partir del coeficiente de Hazen-Williams sugerido en la normativa
colombiana. Es importante aclarar que en los dos diseños obtenidos para el material
GRP, aparentemente la presión en los nudos es igual o superior a la presión mínima
establecida en este trabajo (20 m.c.a), cuando se realiza la comprobación de estos
diseños utilizando la ecuación de Hazen-Williams.
92
Diámetros obtenidos Diseño Red Medellín (GRP)

1400
1200
Diámetro (mm)
1000
800
600
CHw: 147
400
200 CHw:150
0
0 200 400 600 800
Número tubo
Gráfica 4.3-3 Diámetros obtenidos en el diseño de la red matriz de Medellín para el material GRP.
Presiones obtenidas Diseño Red Medellín (GRP)

260
240
220
200
180
160
140
120
100
80
60
40
20
0
0 20 40 60 80 100 120 140 160 180 200 220 240 260
para el material GRP.
4.3.1.3 Acero
Con respecto al diseño de la red matriz de Medellín utilizando el Acero como material de
diseño, específicamente en el sector delimitado por el Embalse Ayurá, en la Gráfica 4.3-5
se puede observar que los diámetros obtenidos para el coeficiente de Hazen-Williams
corregido para este material, en su mayoría, sobreestiman la magnitud de los diámetros
cuando se compara con los obtenidos a partir del coeficiente de Hazen-Williams sugerido
en la normativa colombiana, debido a la restricción del diámetro mínimo que tiene este
para su correcta aplicación.
93
Esto se ve reflejado en la Gráfica 4.3-6, donde se muestra las presiones obtenidas

la mayoría de los nudos se sobrestima las presiones en los nudos cuando se compara
con las presiones obtenidas para el coeficiente de Hazen-Williams sugerido en la
normativa colombiana. Es importante aclarar que en los dos diseños obtenidos para el
Acero, aparentemente la presión en los nudos es igual o superior a la presión mínima
Diámetros obtenidos Diseño Red Medellín (Acero)

1400
1200
Diámetro (mm)
1000
800
600 CHw: 139
400
CHw:120
200
0
0 200 400 600 800
Número tubo
Gráfica 4.3-5 Diámetros obtenidos en el diseño de la red matriz de Medellín para el material Acero.
Presiones obtenidas Diseño Red Medellín (Acero)

260
240
220
200
180
160
140
120
100
80
60
40
20
0
0 20 40 60 80 100 120 140 160 180 200 220 240 260
para el material Acero.
94
4.3.2 Red Hanoi
Luego de haber obtenido los diseños de la red matriz de Medellín, para los distintos
materiales de estudio, se procedió a diseñar la red Hanoi a fin de evaluar los coeficientes
de Hazen-Williams corregidos, comparando los diseños obtenidos por medio de estos con
los calculados a partir de los coeficientes de Hazen-Williams sugeridos en la normativa
colombiana, teniendo en cuenta los parámetros de diseño descritos en el Numeral 3.7.
4.3.2.1 CCP
En relación con el diseño de la red Hanoi para el material CCP, en la Gráfica 4.3-7 se
puede apreciar que los diámetros obtenidos para el coeficiente de Hazen-Williams
corregido para este material, en su mayoría, sobreestiman o en su defecto igualan la
magnitud de los diámetros con los obtenidos a partir del coeficiente de Hazen-Williams
sugerido en la normativa colombiana, debido a la restricción del diámetro mínimo que
tiene este para su correcta aplicación.
Lo anterior se evidencia en la Gráfica 4.3-8, donde se muestra las presiones obtenidas

la mayoría de los nudos se sobrestima o en su defecto se igualan cuando se comparan
CCP, aparentemente la presión en los nudos es igual o superior a la presión mínima
95
Diámetros obtenidos Diseño Red Hanoi (CCP)

1600
1400
1200
Diámetro (mm)
1000
800
600 CHw: 144
400 CHw: 120
200
0
0 10 20 30 40
Número tubería
Gráfica 4.3-7 Diámetros obtenidos en el diseño de la red Hanoi para el material CCP.
Presiones obtenidas Diseño Red Hanoi (CCP)

120
100
80
60
40
20
0
0 20 40 60 80 100 120 140
Gráfica 4.3-8 Diagrama de dispersión presiones obtenidas en el diseño de la red Hanoi para el material
CCP.
4.3.2.2 GRP
Con respecto al diseño de la red Hanoi para el material GRP, en la Gráfica 4.3-9 se puede
apreciar que los diámetros obtenidos para el coeficiente de Hazen-Williams corregido para
este material, en su mayoría, sobreestiman o en su defecto igualan la magnitud de los
diámetros con los obtenidos a partir del coeficiente de Hazen-Williams sugerido en la
normativa colombiana, debido a la restricción del diámetro mínimo que tiene este para su
correcta aplicación.
96
Lo anterior se evidencia en la Gráfica 4.3-10, donde se muestra las presiones obtenidas

la mayoría de los nudos se sobrestima o en su defecto se igualan cuando se comparan
GRP, aparentemente la presión en los nudos es igual o superior a la presión mínima
Diámetros obtenidos Diseño Red Hanoi (GRP)

1400
1200
Diámetro (mm)
1000
800
600 CHw: 147
400
CHw: 150
200
0
0 10 20 30 40
Número tubería
Gráfica 4.3-9 Diámetros obtenidos en el diseño de la red Hanoi para el material GRP.
Presiones obtenidas Diseño Red Hanoi (GRP)

120
100
80
60
40
20
0
0 20 40 60 80 100 120 140
Gráfica 4.3-10 Diagrama de dispersión presiones obtenidas en el diseño de la red Hanoi para el
material GRP.
97
4.3.2.3 Acero
Por último, en relación con el diseño de la red Hanoi para el Acero, en la Gráfica 4.3-11 se
puede observar que los diámetros obtenidos para el coeficiente de Hazen-Williams
corregido para este material, en su mayoría, sobreestiman o en su defecto igualan la
magnitud de los diámetros con los obtenidos a partir del coeficiente de Hazen-Williams
sugerido en la normativa colombiana, debido a la restricción del diámetro mínimo que
tiene este para su correcta aplicación.
Lo anterior, se refleja en la Gráfica 4.3-12 donde se muestra las presiones obtenidas

la mayoría de los nudos se sobrestima estas cuando se comparan con las obtenidas a
partir del coeficiente de Hazen-Williams sugerido en la normativa colombiana. Es
importante aclarar que en los dos diseños obtenidos para el Acero, aparentemente la
presión en los nudos es igual o superior a la presión mínima establecida en este trabajo
(20 m.c.a), cuando se realiza la comprobación de estos diseños utilizando la ecuación de
Hazen-Williams.
Diámetros obtenidos Diseño Red Hanoi (Acero)

1600
1400
1200
Diámetro (mm)
1000
800
600 CHw: 139
400 CHw: 120
200
0
0 10 20 30 40
Número tubería
Gráfica 4.3-11 Diámetros obtenidos en el diseño de la red Hanoi para el material Acero.
98
Presiones obtenidas Diseño Red Hanoi (Acero)

120
100
80
60
40
20
0
0 20 40 60 80 100 120 140
Gráfica 4.3-12 Diagrama de dispersión presiones obtenidas en el diseño de la red Hanoi para el
material Acero.
4.4 Validación diseño redes matrices
Luego de haber obtenido los diseños de la red matriz de Medellín y la red Hanoi, para los
distintos materiales de estudio, utilizando los coeficientes de Hazen-Williams corregidos y
los sugeridos en la normativa colombiana, se procede a validar estos diseños a fin de
verificar si estos en realidad cumplen con la presión mínima establecido en los parámetros
de diseño de estas redes matrices. Para ello, se realizaran comprobaciones de estos
diseños, utilizando la ecuación de Darcy-Weisbach, siguiendo el procedimiento descrito
en el Numeral 3.8 de este trabajo, para luego comparar estos resultados con los
obtenidos en el Numeral 4.3. Los resultados de las validaciones de estos diseños se
ilustran a continuación para cada una de las redes de estudio, teniendo en cuenta los
siguientes aspectos: error en la estimación de las pérdidas por fricción y distribución de
presiones en la red, y cumplimiento de la presión mínima de diseño establecida
(20 m.c.a).
99
Con base en los resultados obtenidos en cuanto al diseño de la red matriz de Medellín,
específicamente el sector delimitado por el Embalse Ayurá, utilizando los coeficientes de
Hazen-Williams corregidos y los sugeridos en la normativa colombiana, se procedió a
validar estos diseños para verificar el cumplimiento de los requisitos hidráulico de esta red
matriz, tal como se ilustra a continuación para cada uno de los materiales seleccionados
para este estudio.
4.4.1.1 CCP
• Dispersión presiones
De acuerdo a la metodología propuesta en el Numeral 3.8, en la Gráfica 4.4-1 se ilustran

los resultados obtenidos en cuanto a la comprobación del diseño obtenido a partir del
coeficiente de Hazen-Williams corregido en este trabajo y sugerido en la normativa
colombiana para este material. Igualmente, se muestra la comparación de estos con los
resultados descritos en el Numeral 4.3.1.1, a fin de verificar si estos diseños cumplen en
realidad con las restricciones de tipo de hidráulico de esta red en relación con los errores
en la distribución de presiones y cumplimiento de la presión mínima de diseño (20 m.c.a).
Con respecto a los errores en la distribución de presiones, en este trabajo se utilizó un

diagrama de dispersión de presiones, con el propósito de establecer si los diseños
obtenidos a partir de la ecuación de Hazen-Williams sobreestiman o subestiman las
presiones que ocurren realmente en la red matriz de Medellín, dependiendo si estas se
ubican por encima o por debajo de la línea de 45o descrita en la Gráfica 4.1-1,
respectivamente. Teniendo en cuenta lo anterior, se puede apreciar en esta gráfica que el
diseño obtenido a partir del coeficiente de Hazen-Williams corregido en este trabajo (CHw:
144) tiende a sobreestimar las presiones en la red matriz de Medellín. Así mismo, se
puede ver en la Gráfica 4.1.1 que este diseño no cumple con la presión mínima de diseño
de esta red (20 m.c.a). Por otro lado, el diseño obtenido a partir del coeficiente de Hazen-
Williams sugerido en la normativa colombiana (CHw: 120) subestima las presiones en
esta red, pero se encuentran por encima de la presión mínima de diseño.
100
Diagrama dispersión de presiones diseño Red Medellín (CCP, ks: 0.12 mm)
240
220
200
180
160
140
120
Presiones Hazen-Williams (m.c.a)
100
80
60
40
20
0
-20
-40
-60
-80
-100
-120
-140
-160
-180
-200
-220
-240
-240 -220 -200 -180 -160 -140 -120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120 140 160 180 200 220 240
Comprobación Presiones diseño Darcy-Weisbach (m.c.a)
CHw: 144 CHw: 120
Gráfica 4.4-1 Diagrama de Dispersión de presiones diseño Red matriz de Medellín para el material
CCP, específicamente para los coeficientes de Hazen-Williams iguales a 144 y 120.
• Errores pérdidas por fricción
Con base en las pérdidas por fricción obtenidas para los diseños de la red matriz de
Medellín tanto para el coeficiente de Hazen-Williams corregido para este material y el
sugerido en la normativa colombiana, utilizando la ecuación de Hazen-Williams, y las
pérdidas por fricción calculadas a partir de la comprobación de estos diseños, utilizando la
ecuación de Darcy-Weisbach, se graficaron los errores en la estimación de estas pérdidas
(Ver Gráfica 4.4-2 y Gráfica 4.4-3). A partir de lo anterior, se puede apreciar en la Gráfica
4.4.2 que para el diseño obtenido para el coeficiente de Hazen-Williams corregido para el
CCP, en su mayoría subestima las pérdidas por fricción (hasta en un 20%), dado que se
manejan velocidades superiores a 1 m/s en la red matriz de Medellín. De otro lado, para
el diseño calculado a partir del coeficiente de Hazen-Williams sugerido en la normativa
colombiana (Ver Gráfica 4.4-3), se puede observar que este en su mayoría sobreestima
las pérdidas por fricción (hasta un 30%), para todos los diámetros obtenidos.
101
Error Pérdidas por Fricción Diseño Red Medellín (CCP, ks: 0.12 mm, CHw:144 )
50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10% d=2.0 m
0% d=1.5 m
-10% d=1.0 m
v=0.5 m/s
v= 1 m/s d=0.5 m
-20% v=2 m/s
v= 3 m/s Errores Red
-30%
v= 5m/s
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.4-2 Errores pérdidas por fricción en el diseño de la Red matriz de Medellín, para el material
CCP (CHw: 144).
Error Pérdidas por Fricción Diseño Red Medellín (CCP, ks: 0.12 mm, CHw:120 )
50%
40%
30%
v=0.5 m/s
20% v= 1 m/s
(Shw-Sdw)/Sdw
d=2.5 m
v=2 m/s
10%
v= 3 m/s d=2.0 m
0% v= 5m/s d=1.5 m
-10% d=1.0 m
d=0.5 m
-20%
Errores Red
-30%
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
CCP (CHw: 120).
102
4.4.1.2 GRP
Con base en las presiones obtenidas en los diseños de la red matriz de Medellín tanto
para el coeficiente de Hazen-Williams corregido para este material y el sugerido en la
normativa colombiana, utilizando la ecuación de Hazen-Williams, y las presiones
calculadas a partir de la comprobación de estos diseños, utilizando la ecuación de Darcy-
Weisbach, se gráficaron en un diagrama de dispersión de presiones a fin de establecer si
estos sobreestiman o subestiman las presiones que ocurren realmente en la red matriz de
Medellín. A partir de lo anterior, se puede ver en la Gráfica 4.4-4 que el diseño obtenido a
partir del coeficiente de Hazen-Williams corregido en este trabajo (CHw: 147) tiende a
igualar las presiones en la red matriz de Medellín, cumpliendo con la presión mínima de
diseño de esta red (20 m.c.a). Por otro lado, el diseño obtenido a partir del coeficiente de
Hazen-Williams sugerido en la normativa colombiana (CHw: 150) para este material,
subestima las presiones en esta red y en algunos nudos no se logra cumplir con la
presión mínima de diseño.
Diagrama dispersión de presiones Diseño Red Medellín (GRP, ks: 0.03 mm)
240
220
200
180
160
140
120
100
80
60
40
20
0
-20
-40
-60
-80
-100
-120
-140
-160
-180
-200
-220
-240
-240 -220 -200 -180 -160 -140 -120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120 140 160 180 200 220 240
CHw: 147 CHw: 150
GRP, específicamente para los coeficientes de Hazen-Williams iguales a 147 y 150.
103
ecuación de Darcy-Weisbach, se gráficaron los errores en la estimación de estas pérdidas
(Ver Gráfica 4.4-5y Gráfica 4.4-6). A partir de lo anterior, se puede apreciar en la Gráfica
4.4.5 que para el diseño obtenido para el coeficiente de Hazen-Williams corregido para el
GRP, en su mayoría se ubica en el rango del número de Reynolds establecido para dicho
coeficiente, donde los errores en las pérdidas no superan el 5%. De otro lado, para el
diseño calculado a partir del coeficiente de Hazen-Williams sugerido en la normativa
colombiana (Ver Gráfica 4.4-6), se puede observar que para diámetros inferiores a 500
mm, este en su mayoría subestima las pérdidas por fricción (hasta un 14%), para todos
los diámetros obtenidos para números de Reynolds superiores a 300.000.
Errores Pérdidas por Fricción Diseño Red Medellín (GRP, ks:0.03 mm, CHw:147 )
50%
40%
30%
20% d=2.5 m
(Shw-Sdw)/Sdw
10% d=2.0 m
0% d=1.5 m
v=2 m/s v=3 m/s d=1.0 m
-10% v= 5m/s
-20% d=0.5 m
Errores red
-30%
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
GRP (CHw: 147).
104
Errores Pérdidas por Fricción Diseño Red Medellín (GRP, ks: 0.03 mm, CHw:150 )
50%
40%
30%
20% d=2.5 m
(Shw-Sdw)/Sdw
10% d=2.0 m
0% d=1.5 m
-10% v=2 m/s v=3 m/s d=1.0 m

v= 5m/s
-20% d=0.5 m
Errores Red
-30%
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
GRP (CHw: 150).
4.4.1.3 Acero
Con base en las presiones obtenidas en los diseños de la red matriz de Medellín tanto
normativa colombiana, utilizando la ecuación de Hazen-Williams, y las presiones
Weisbach, se gráficaron en un diagrama de dispersión de presiones a fin de establecer si
estos sobreestiman o subestiman las presiones que ocurren realmente en la red matriz de
Medellín. A partir de lo anterior, se puede ver en la Gráfica 4.4-7 que el diseño obtenido a
partir del coeficiente de Hazen-Williams corregido en este trabajo (CHw: 139) tiende a
diseño de esta red (20 m.c.a). Por otro lado, el diseño obtenido a partir del coeficiente de
Hazen-Williams sugerido en la normativa colombiana (CHw: 120) para este material,
sobrestima las presiones en esta red y en algunos nudos no se logra cumplir con la
presión mínima de diseño, dado que son inferiores a 20 m.c.a o en el peor de los casos
negativas.
105
Diagrama dispersión de presiones en la Red Medellín

(Acero, ks: 0.45 mm)
240
220
200
180
160
140
120
100
80
60
40
20
0
-20
-40
-60
-80
-100
-120
-140
-160
-180
-200
-220
-240
-240-220-200-180-160-140-120-100 -80 -60 -40 -20 0 20 40 60 80 100 120 140 160 180 200 220 240
CHw: 139 CHw: 120
Acero, específicamente para los coeficientes de Hazen-Williams iguales a 139 y 120.
(Ver Gráfica 4.4-8 y Gráfica 4.4-9). A partir de lo anterior, se puede apreciar en la Gráfica
4.4-8 que para el diseño obtenido para el coeficiente de Hazen-Williams corregido para el
Acero, en su mayoría se ubica en el rango del número de Reynolds establecido para
dicho coeficiente, donde los errores en las pérdidas no superan el 5%. Sin embargo, par
números de Reynolds superiores a 200.000 se presenta subestimación de las pérdidas
por fricción hasta en un 35%. De otro lado, para el diseño calculado a partir del
coeficiente de Hazen-Williams sugerido en la normativa colombiana (Ver Gráfica 4.4-9),
se puede observar que para diámetros inferiores a 1000 mm, este en su mayoría
subestima las pérdidas por fricción (hasta un 25%), para todos los diámetros obtenidos
para números de Reynolds superiores a 300.000.
106
Errores Pérdidas por Fricción Diseño Red Medellín (Acero, ks: 0.45 mm, CHw: 139)
50%
40%
30%
20%
(Shw-Sdw)/Sdw
10% d=2.5 m
d=2.0 m
0%
d=1.5 m
-10%
d=1.0 m
-20%
v=0.5 m/s Errores Red
-30% v= 1 m/s
v=2 m/s
-40% v=3 m/s
v=5 m/s
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.4-8 Errores pérdidas por fricción en el diseño de la Red matriz de Medellín, para el Acero
(CHw: 139).
Errores Pérdidas por Fricción Diseño Red Medellín (Acero, ks:0.45 mm, CHw:120)
50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10%
v=0.5 m/s d=2.0 m
0%
v= 1 m/s d=1.5 m
-10% v=2 m/s
d=1.0 m
-20% v=3 m/s
Errores Red
-30% v=5 m/s
-40%
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Acero (CHw: 120).
107
4.4.2 Red Hanoi
Con base en los resultados obtenidos en cuanto al diseño de la red Hanoi, utilizando los
coeficientes de Hazen-Williams corregidos y los sugeridos en la normativa colombiana, se
procedió a validar estos diseños para verificar el cumplimiento de los requisitos hidráulicos
de esta red matriz, tal como se ilustra a continuación para cada uno de los materiales
seleccionados para este estudio.
4.4.2.1 CCP
Con base en las presiones obtenidas en los diseños de la red matriz Hanoi tanto para el
coeficiente de Hazen-Williams corregido para este material y el sugerido en la normativa
colombiana, utilizando la ecuación de Hazen-Williams, y las presiones calculadas a partir
de la comprobación de estos diseños, utilizando la ecuación de Darcy-Weisbach, se
gráficaron en un diagrama de dispersión de presiones a fin de establecer si estos
sobreestiman o subestiman las presiones que ocurren realmente en la red matriz de
Medellín. A partir de lo anterior, se puede ver en la Gráfica 4.4-10 que el diseño obtenido
a partir del coeficiente de Hazen-Williams corregido en este trabajo (CHw: 144) tiende a
igualar las presiones en la red matriz de Medellín, pero solamente en un nodo de ésta no
se cumple con la presión mínima de diseño de esta red (20 m.c.a). Por otro lado, el diseño
obtenido a partir del coeficiente de Hazen-Williams sugerido en la normativa colombiana
(CHw: 120) para este material, subestima las presiones de la red y cumplen con la presión
minima de diseño para todos los nudos de la esta red.
108
Diagrama dispersión de presiones en la Red Hanoi (CCP, ks: 0.12 mm)

120
100
80
60
40
20
0
-20
-40
-60
-80
-100
-120
-120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120
CHw: 144 CHw: 120
Gráfica 4.4-10 Diagrama de Dispersión de presiones diseño Red Hanoi para el material CCP,
específicamente para los coeficientes de Hazen-Williams iguales a 144 y 120.
Con base en las pérdidas por fricción obtenidas para los diseños de la red Hanoi para el
colombiana, utilizando la ecuación de Hazen-Williams, y las pérdidas por fricción
Weisbach, se gráficaron los errores en la estimación de estas pérdidas (Ver Gráfica
4.4-11 y Gráfica 4.4-12). A partir de lo anterior, se puede apreciar en la Gráfica 4.4.11
que para el diseño obtenido para el coeficiente de Hazen-Williams corregido para el CCP,
en su mayoría subestiman las pérdidas por fricción (hasta un 15%) para números de
Reynolds superiores a 400.000. De otro lado, para el diseño calculado a partir del
se puede observar que para todos los diámetros obtenidos se sobrestiman las pérdidas
por fricción en un rango de 15 a 30%.
109
Error Pérdidas por Fricción Diseño Red Hanoi (CCP, ks: 0.12 mm, CHw:144 )
50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10%
d=2.0 m
0% d=1.5 m
-10% d=1.0 m
v=0.5 m/s
v= 1 m/s d=0.5 m
-20% v=2 m/s
v= 3 m/s Errores Red
v= 5m/s
-30%
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.4-11 Errores pérdidas por fricción en el diseño de la Hanoi, para el material CCP (CHw: 144).
Error Pérdidas por Fricción Diseño Red Hanoi (CCP, ks: 0.12 mm, CHw:120 )
50%
40%
30% v=0.5 m/s

20% v= 1 m/s
v=2 m/s
(Shw-Sdw)/Sdw
d=2.5 m
10% v= 3 m/s d=2.0 m
v= 5m/s
0% d=1.5 m
-10% d=1.0 m
d=0.5 m
-20%
Errores Red
-30%
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.4-12 Errores pérdidas por fricción en el diseño de la Red Hanoi, para el material CCP (CHw:
120).
110
4.4.2.2 GRP
Con base en las presiones obtenidas en los diseños de la red Hanoi tanto para el
diseño de esta red (20 m.c.a). A diferencia de este coeficiente, el diseño obtenido a partir
del coeficiente de Hazen-Williams sugerido en la normativa colombiana (CHw: 150) no
cumple con dicha presión en un nodo de esta red.
Diagrama dispersión de presiones en la Red Hanoi (GRP, ks: 0.03 mm)

120
100
80
60
40
20
0
-20
-40
-60
-80
-100
-120
-120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120
CHw: 147 CHw: 150
Gráfica 4.4-13 Diagrama de Dispersión de presiones diseño Red Hanoi para el material GRP,
específicamente para los coeficientes de Hazen-Williams iguales a 147 y 150.
111
Con base en las pérdidas por fricción obtenidas para los diseños de la red Hanoi tanto
normativa colombiana, utilizando la ecuación de Hazen-Williams, y las pérdidas por
fricción calculadas a partir de la comprobación de estos diseños, utilizando la ecuación de
Darcy-Weisbach, se gráficaron los errores en la estimación de estas pérdidas (Ver Gráfica
4.4-14y Gráfica 4.4-15). A partir de lo anterior, se puede apreciar en la Gráfica 4.4.14 que
para el diseño obtenido para el coeficiente de Hazen-Williams corregido para el GRP, en
su mayoría se ubica en el rango del número de Reynolds establecido para dicho
coeficiente, donde los errores en las pérdidas no superan el 3%. De otro lado, para el
diseño calculado a partir del coeficiente de Hazen-Williams sugerido en la normativa
colombiana (Ver Gráfica 4.4-15), se puede observar que para diámetros inferiores a 500
mm, este en su mayoría subestima las pérdidas por fricción (hasta un 5%), para números
de Reynolds entre 100.000 y 1.000.000.
Errores Pérdidas por Fricción Diseño Red Hanoi (GRP, ks:0.03 mm, CHw:147)
50%
40%
30%
20% d=2.5 m
(Shw-Sdw)/Sdw
10% d=2.0 m
0% d=1.5 m
-10% v=2 m/sv=3 m/sv= 5m/s d=1.0 m
-20% d=0.5 m
-30% Errores Red
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
GRP (CHw: 147).
112
Errores Pérdidas por Fricción Diseño Red Hanoi (GRP, ks: 0.03 mm, CHw:150 )
50%
40%
30%
20% d=2.5 m
(Shw-Sdw)/Sdw
10% d=2.0 m
0% d=1.5 m
-10% v=2 m/sv=3 m/sv= 5m/s d=1.0 m
-20% d=0.5 m
-30% Errores Red
-40%
-50%
1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
GRP (CHw: 150).
4.4.2.3 Acero
Con base en las presiones obtenidas en los diseños de la red Hanoi tanto para el
igualar sobrestimar las presiones en la red, cumpliendo con la presión mínima de diseño
de esta red (20 m.c.a). Por otro lado, el diseño obtenido a partir del coeficiente de Hazen-
Williams sugerido en la normativa colombiana (CHw: 120) para este material, iguala las
presiones de la red y en todos los nodos cumple con la presión mínima de diseño.
113
Diagrama dispersión de presiones en la Red Hanoi (Acero, ks: 0.45 mm)

120
100
80
60
40
20
0
-20
-40
-60
-80
-100
-120
-120 -100 -80 -60 -40 -20 0 20 40 60 80 100 120
CHw: 139 CHw: 120
Gráfica 4.4-16 Diagrama de Dispersión de presiones diseño Red Hanoi para el Acero, específicamente
para los coeficientes de Hazen-Williams iguales a 139 y 120.
(Ver Gráfica 4.4-17 y Gráfica 4.4-18). A partir de lo anterior, se puede apreciar en la
Gráfica 4.4-17 que para el diseño obtenido para el coeficiente de Hazen-Williams
corregido para el Acero, en su mayoría, el error en las pérdidas por fricción se ubican en
el rango del número de Reynolds establecido para dicho coeficiente, donde los errores en
las pérdidas no superan el 5%, dado que para diámetros iguales 1000 mm se manejan
velocidades de hasta 0.5 m/s. De otro lado, para el diseño calculado a partir del
se puede observar que el error en las pérdidas no supera el 5% para números de
Reynolds superiores a 100.000, teniendo en cuenta las velocidades que se presentan
para los diámetros obtenidos.
114
Errores Pérdidas por Fricción Diseño Red Hanoi (Acero, ks: 0.45 mm, CHw: 139)
50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10%
d=2.0 m
0%
d=1.5 m
-10%
d=1.0 m
-20%
v=0.5 m/s Errores Red
-30% v= 1 m/s
-40% v=2 m/s
v=3 m/s
-50%
v=5 m/s
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.4-17 Errores pérdidas por fricción en el diseño de la Red Hanoi, para el Acero (CHw: 139).
Errores Pérdidas por Fricción Diseño Red Hanoi (Acero, ks:0.45 mm, CHw:120)
50%
40%
30%
20%
(Shw-Sdw)/Sdw
d=2.5 m
10%
v=0.5 m/s d=2.0 m
0%
v= 1 m/s d=1.5 m
-10%
v=2 m/s d=1.0 m
-20% v=3 m/s
Errores Red
-30% v=5 m/s
-40%
-50%
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 4.4-18 Errores pérdidas por fricción en el diseño de la Red Hanoi, para el Acero (CHw: 120).
115
Capítulo 5. ANÁLISIS DE RESULTADOS

Los resultados obtenidos en la red matriz de Medellín, la red Hanoi y la red matriz
analizada por Bombardelli y Garcia (2003), se procede a analizar los efectos que tiene el
uso de la ecuación de Darcy-Weisbach vs. la ecuación de Hazen-Williams en estas redes
de estudio, en cuanto a la variación del coeficiente de Hazen-Williams y los efectos que
tiene dicha variación en los errores de las pérdidas por fricción y las presiones en los
nudos. Igualmente, se analizarán los resultados obtenidos en relación con los coeficientes
de Hazen-Williams corregidos para los materiales de estudio, así como su efecto en el
diseño de las redes matrices.
5.1 Efecto del uso de la ecuación de Hazen-Williams en las

Redes de Estudio
En esta parte del trabajo se mostrará el análisis de resultados obtenidos en relación con la
variabilidad del coeficiente de Hazen-Williams en las redes estudio y su efecto en la
estimación de las pérdidas por fricción y distribución de presiones, cuando se utiliza un
coeficiente de Hazen-Williams constante.
• CCP
En relación con la variación de los coeficientes de Hazen-Williams para esta red, se
puede ver en la Gráfica 4.1-1 que el coeficiente de Hazen-Williams sugerido para el
material CCP en la normativa colombiana, para un rango de número de Reynolds entre
10000 y 1.000.0000, subestima la magnitud real del coeficiente de Hazen-Williams para
los diámetros y condiciones de flujo que se presentan en ésta, lo cual representa según la
ecuación de Darcy-Weisbach, que dicho valor sobreestima las pérdidas por fricción. Este
comportamiento se puede ver en detalle en la Gráfica 4.1-4, donde se ilustran los errores
en la sobrestimación del cálculo de las pérdidas por fricción, los cuales varían desde 0
hasta el 20%. Lo anterior, implica que este coeficiente subestima las presiones de la red,
como se puede apreciar en la Gráfica 4.1-7, para algunos nudos de la misma.
116
En contraste, para números de Reynolds superiores a 9.000.000, condición de flujo que

se presenta cerca al Embalse Ayurá, el coeficiente de Hazen-Williams sugerido en la
normativa colombiana para el material CCP tiende a sobreestimar su magnitud real (Ver
Gráfica 4.1-1). Esto implica, según la ecuación de Darcy-Weisbach, que en la red se
estarían generando mayores pérdidas y, por lo tanto, dicho coeficiente estaría
subestimando dichas pérdidas (Ver Gráfica 4.1-4), situación que se debe evitar para
poder cumplir con las presiones requeridas en la misma. Este comportamiento se refleja
en la Gráfica 4.1-7 cerca al embalse de Ayurá, donde se sobreestiman las presiones en
los nudos hasta en un 10%, cuando se utiliza la ecuación de Hazen-Williams.
• GRP
Según los resultados obtenidos en relación con la variación del coeficiente de Hazen-
Williams para el material GRP, se puede apreciar en la Gráfica 4.1-2, que la variación del
coeficiente de Hazen-Williams alrededor del valor 150, sugerido en la normativa
colombiana, no es significativa para las condiciones de flujo que se presentan en la red
matriz de Medellín, dado que los errores obtenidos en cuanto a las pérdidas por fricción,
utilizando la ecuación de Darcy-Weisbach y Hazen-Williams (Ver Gráfica 4.1-5), no
sobrepasan el 10% para la mayoría de las tuberías de esta red. Igualmente, este efecto
se puede ver también en la diferencia de presiones en los nudos (Ver Gráfica 4.1-8), la
cual no sobrepasa los 15% en la sobreestimación de los nudos. Dicha sobreestimación
se explica, dado que en todas las tuberías de esta red se subestima el cálculo de las
pérdidas por fricción, según la ecuación de Darcy-Weisbach (Ver Gráfica 4.1-5), situación
riesgosa pues en la red posiblemente no se estaría cumpliendo con las presiones
requeridas en la red matriz de Medellín.
• Acero
En relación con la variación de los coeficientes de Hazen-Williams para esta red, se

puede ver en la Gráfica 4.1-3 que el coeficiente de Hazen-Williams sugerido para el Acero
en la normativa colombiana, para un rango de número de Reynolds entre 10000 y
300.000, subestima la magnitud real del coeficiente de Hazen-Williams para los diámetros
y condiciones de flujo que se presentan en esta red, lo cual representa según la ecuación
de Darcy-Weisbach, que dicho valor sobreestima las pérdidas por fricción. Este
117
comportamiento se puede ver en detalle en la Gráfica 4.1-6, donde se ilustran los errores
en la sobrestimación del cálculo de las pérdidas por fricción, los cuales varían desde 0
hasta el 40%.
En contraste, para números de Reynolds superiores a 300.000, condición de flujo que se

presenta cerca al Embalse Ayurá, el coeficiente de Hazen-Williams sugerido en la
normativa colombiana para el material CCP tiende a sobreestimar la magnitud real de
este coeficiente (Ver Gráfica 4.1-3). Esto implica, según la ecuación de Darcy-Weisbach,
que en la red se estarían generando mayores pérdidas y, por lo tanto, dicho coeficiente
estaría subestimando dichas pérdidas (Ver Gráfica 4.1-6), situación que se debe evitar
para poder cumplir con las presiones requeridas en la misma. Este comportamiento se
refleja en la Gráfica 4.1-9 cerca al embalse de Ayurá, donde se sobreestiman las
presiones en los nudos hasta en un 70%, cuando se utiliza la ecuación de Hazen-
Williams.
5.1.2 Red Hanoi
• CCP
Con respecto a la variación de los coeficientes de Hazen-Williams para esta red, se puede
ver en la Gráfica 4.1-10 que el coeficiente de Hazen-Williams sugerido para el material
CCP en la normativa colombiana, subestima la magnitud real del coeficiente de Hazen-
Williams para los diámetros y condiciones de flujo que se presentan en esta red, lo cual
representa según la ecuación de Darcy-Weisbach, que dicho valor sobreestima las
pérdidas por fricción. Este comportamiento se puede ver en detalle en la Gráfica 4.1-14,
donde se ilustran los errores en la sobrestimación del cálculo de las pérdidas por fricción,
los cuales varían desde 0 hasta el 32%. Lo anterior implica que este coeficiente subestima
las presiones de esta red, como se puede apreciar en la Gráfica 4.1-17, en todos los
nudos de la red, desde 0 hasta el 85% en los sumideros de esta red.
118
• GRP
De acuerdo con los resultados obtenidos en relación con la variación del coeficiente de
Hazen-Williams para el material GRP, se puede apreciar en la Gráfica 4.1-12, que la
variación del coeficiente de Hazen-Williams alrededor del valor 150, sugerido en la
normativa colombiana, no es significativa para las condiciones de flujo que se presentan
en la red Hanoi, dado que los errores obtenidos en cuanto a las pérdidas por fricción,
utilizando la ecuación de Darcy-Weisbach y Hazen-Williams (Ver Gráfica 4.1-15), no
sobrepasan el 10% para la mayoría de las tuberías de esta red. Igualmente, este efecto
se puede ver también en la diferencia de presiones en los nudos (Ver Gráfica 4.1-18), la
cual no sobrepasa los 10% en la sobreestimación de los nudos. Dicha sobreestimación
se explica, ya que en todas las tuberías de esta red se subestima el cálculo de las
pérdidas por fricción, según la ecuación de Darcy-Weisbach (Ver Gráfica 4.1-15),
situación riesgosa pues en la red posiblemente no se estaría cumpliendo con las
presiones requeridas.
• Acero
En tratándose de la variación de los coeficientes de Hazen-Williams para esta red, se

puede ver en la Gráfica 4.1-13 que el coeficiente de Hazen-Williams sugerido para el
material CCP en la normativa colombiana, en su mayoría, subestima la magnitud real del
coeficiente de Hazen-Williams para los diámetros y condiciones de flujo que se presentan
en ésta, lo cual representa según la ecuación de Darcy-Weisbach, que dicho valor
sobreestima las pérdidas por fricción. Este comportamiento se puede ver en detalle en la
Gráfica 4.1-16, donde se ilustran los errores en la sobrestimación del cálculo de las
pérdidas por fricción, los cuales varían desde 0 hasta el 10%. Lo anterior, implicaría que
este coeficiente subestima las presiones de esta red. Pese a ello, este comportamiento no
se presenta en cuanto a los errores en la estimación de las presiones en los nudos (Ver
Gráfica 4.19), dado que en muchas de las tuberías de esta red se subestima las pérdidas
por fricción, según la ecuación de Darcy-Weisbach, ya que en éstas los números de
Reynolds son superiores a 300.000, rango donde se subestima dichas pérdidas hasta en
un 17%.
119
Con base en los resultados obtenidos en la red matriz estudiada por Bombardelli y García
(2003) y específicamente en los tres tramos de la tubería principal de esta red (diámetro=
2.9 m), se procedió a analizar los mismos en cuanto al efecto de la variación que tiene el
coeficiente de Hazen-Williams en la estimación de las pérdidas por fricción en esta red.
En relación con la variabilidad del coeficiente de Hazen-Williams, se puede ver en las
Gráficas 4.1-19, 4.1-21 y 4.1-23, que el coeficiente de Hazen-Williams utilizado para el
diseño original de esta red (CHw: 120) sobreestima la magnitud real que presenta dicho
coeficiente (entre 85 y 97), según la ecuación de Darcy-Weisbach, para las rugosidades
absolutas y condiciones de flujo que presentan (Re=3.8x106 y V=1.655 m/s) los tres
tramos de la tubería principal de esta red. Lo anterior implicaría, según la ecuación de
Darcy-Weisbach, que el coeficiente de Hazen-Williams utilizado en el diseño original de
esta red subestima las pérdidas por fricción en la tubería principal de esta. Esto se valida
de acuerdo con las Gráficas 4.1-20, 4.1-22 y 4.1-24, donde se muestra que dicho
coeficiente subestima las pérdidas por fricción en un rango entre 32 y 42%, dependiendo
el tramo de estudio, para un número de Reynolds igual a 3.8x106. Este hecho representa
un riesgo para la red matriz analizada por Bombardelli y García (2003), dado que en esta
se podría estar perdiendo mayor cantidad de energía que la esperada, según la ecuación
de Darcy-Weisbach, lo cual no permitiría predecir con certeza las presiones que se
presentan en la misma.
5.2 Corrección Coeficientes de Hazen-Williams
Teniendo en cuenta los efectos que puede tener la variabilidad del coeficiente de Hazen-
Williams en las redes de estudio, específicamente en la estimación de las pérdidas por
fricción para los materiales CCP, GRP y Acero, se sugirió modificar los valores sugeridos
por la normativa colombiana para estos materiales, a fin de minimizar el error en las
pérdidas por fricción, hasta un 5%, en redes matrices siguiendo el Paso 7 de la
metodología propuesta en este trabajo.
120
A partir de los resultados de esta metodología, para el material CCP se obtuvo un

coeficiente de Hazen-Williams igual a 144, el cual no permite, en el diseño de redes
matrices, que las tuberías con diámetros de 2.5 metros sobreestimen las pérdidas por
fricción en un 5% y, para tuberías de diámetro igual a 500 mm, se subestimen en la
misma magnitud (Ver Gráfica 4.2-1). Adicionalmente se puede observar en la Gráfica 4.2-
1 que este valor sugerido permite manejar velocidades de hasta 0.95 m/s para minimizar
en lo posible el error en las pérdidas por fricción, ya que si en la red matriz se presentan
velocidades superiores a este límite se corre el riesgo de que la ecuación de Hazen-
Williams subestime las pérdidas por fricción en una red, lo cual no permitiría cumplir con
los requisitos de presión mínima en la misma. De igual manera, se puede observar en la
Gráfica 4.2-2, que para el rango de números de Reynolds establecido para aplicar dicho
coeficiente, este valor presenta buena exactitud, con respecto a la ecuación de Darcy-
Weisbach, cuando los factores de fricción se ubican en la zona de FTHL delimitada por la
ecuación deducida por el autor (Flechas, 2010).
De otro lado, con base en la metodología propuesta para la corrección de los coeficientes
de Hazen-Williams, para el Acero se obtuvo un coeficiente igual a 139, el cual permite
controlar la sobrestimación o subestimación de las pérdidas por fricción hasta en un 5%
para velocidades inferiores a 0.3 m/s y específicamente para diámetro superiores a 1000
mm, ya que si en la red matriz se presentan velocidades superiores a este límite se corre
el riesgo de que la ecuación de Hazen-Williams subestime las pérdidas por fricción en una
red, lo cual no permitiría cumplir con los requisitos de presión mínima en la misma.
Además, en relación con el rango de número de Reynolds establecido para dicho
coeficiente, se puede ver en la Gráfica 4.2-6, que este presenta buena exactitud, con
respecto a la ecuación de Darcy-Weisbach, cuando los factores de fricción de este se
ubican en la zona de FTHL definida por el autor (Flechas, 2010).
Por último, a pesar de que el valor sugerido para el material GRP en la normatividad
colombiana muestra un buen resultado en cuanto a los efectos en las pérdidas por fricción
en las redes de estudio, se sugirió modificar dicho valor a 147, a fin de evitar, en lo
máximo posible, la subestimación en las pérdidas por fricción en tuberías, en un 5%,
específicamente para diámetros iguales o inferiores a 500 mm. En relación con el límite
121
de velocidad de este coeficiente, se puede apreciar en la Gráfica 4.2-3 que este permite
manejar velocidades de hasta 5.6 m/s en las tuberías de una red matriz, a diferencia de
los demás materiales, lo cual se refleja en un aumento en el rango de números de
Reynolds en el cual se puede aplicar dicho coeficiente. Sin embargo, se encuentra una
similitud con estos materiales, dado que la buena exactitud de este coeficiente, con
respecto a la ecuación de Darcy-Weisbach, se logra cuando los factores de fricción de
éste se ubican en la zona del FTHL (Ver Gráfica 4.2-4), definida por el autor (2010).
5.3 Validación Diseños redes matrices
A partir de los diseños obtenidos para la red matriz de Medellín y la red Hanoi, para los
distintos materiales de estudio, utilizando los coeficientes de Hazen-Williams corregidos y
los sugeridos en la normativa colombiana, se procede en este numeral a analizar la
validación de los mismos, a fin de establecer si estos en realidad cumplen con los
requisitos hidráulicos de estas redes, en cuanto a la estimación de las pérdidas por
fricción y cumplimiento en la presión mínima de diseño.
• CCP
Teniendo en cuenta los resultados obtenidos en cuanto a la validación de los diseños de

la red Matriz de Medellín, utilizando el coeficiente de Hazen-Williams corregido y sugerido
en la normativa colombiana para este material, se puede observar en la Gráfica 4.1.1 que
el primero (CHw: 144) tiende a sobreestimar las presiones en ésta. Así mismo, se puede
observar en esta gráfica que el diseño no cumple con la presión mínima establecida en el
Numeral 3.7 (20 m.c.a). Por otro lado, el diseño obtenido a partir del coeficiente de Hazen-
Williams sugerido en la normativa colombiana (CHw: 120) subestima las presiones en
esta red, cumpliendo con la presión mínima de diseño.
122
Lo anterior se debe a que el coeficiente de Hazen-Williams corregido, utilizado en el

diseño de esta red, en su mayoría subestima las pérdidas por fricción (hasta en un 20%),
dado que se manejan velocidades superiores a 1 m/s (Ver Gráfica 4.4-3), límite de
aplicabilidad de este coeficiente. De otro lado, para el coeficiente de Hazen-Williams
sugerido en la normativa colombiana, se puede ver en la Gráfica 4.4-3, que éste en su
mayoría sobreestima las pérdidas por fricción (hasta un 30%), para todos los diámetros
obtenidos.
Para entender los errores en las pérdidas por fricción, es necesario analizar en detalle los
tipos de flujo que se presentan en la red matriz de Medellín para el material CCP,
graficando los factores de fricción obtenidos en el diagrama de Moody, delimitando la
zona de transición en dicho diagrama, zona que va desde la línea límite máxima donde se
puede presentar Flujo Turbulento Hidráulicamente Liso (FTHL) hasta la línea límite
mínima donde se puede presentar Flujo Turbulento Hidráulicamente Rugoso (FTHR), a
partir de las ecuaciones de Blasius, Prandtl-Von Kármán (Ver Ecuación 2.1-15 y Ecuación
2.1-16) y las ecuaciones deducidas por el autor (Ver Ecuación 2.1-17 y Ecuación 2.1-18)
para estos tipos de flujo. Lo anterior, se muestra en la siguiente gráfica:
Diagrama de Moody (Red Medellín, CCP, ks: 0.12 mm)

0.1
Factor de fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Hazen-Williams (CHW: 144) Comprobación Darcy-Weisbach (CHW: 144)
Hazen-Williams (CHW:120) Comprobación Darcy-Weisbach (CHW: 120)
Flujo laminar FTHL (Blasius)
Límite FTHL (Colebrook-White) Límite FTHL (Prandtl-Von Kármán)
Límite FTHR (Colebrook-White) Límite FTHR (Prandtl-Von Kármán)
Gráfica 5.3-1 Factores de Fricción obtenidos para el diseño de la red matriz de Medellín para el CCP,
específicamente para coeficientes de Hazen-Williams igual a 144 y 120.
123
A partir de la Gráfica 5.3-1, se puede ver que los factores hidráulicos obtenidos a partir del
diseño de esta red, utilizando la ecuación de Hazen-Williams y el coeficiente corregido
para este material CCP, se localizan en la zona de Flujo Turbulento Hidráulicamente Liso,
para números de Reynolds superiores a 200.000, con base en el límite definido por el
autor para este tipo de flujo. De otro lado, para este mismo rango de números de
Reynolds, los factores de fricción obtenidos en la validación de este diseño, utilizando la
ecuación de Darcy-Weisbach, se localizan en la zona de transición y tienden al límite
mínimo del FTHR. Lo anterior, explica por qué el coeficiente de Hazen-Williams sugerido
en el trabajo para este material subestima las pérdidas por fricción, dado que dicho
coeficiente obliga a que el flujo en las tuberías en esta red sea FTHL, cuando en la
realidad, para la rugosidad absoluta dada para el CCP, el flujo es Turbulento Transicional
para los diámetros obtenidos en este diseño.
Por otro lado, en relación con los factores hidráulicos obtenidos a partir del diseño de la
red matriz de Medellín, utilizando la ecuación de Hazen-Williams y el coeficiente de
Hazen-Williams sugerido en la normativa para el CCP, se puede observar en la Gráfica
5.3-1, que estos se ubican en la zona de Flujo Turbulento Transicional y tienden al límite
mínimo del FTHR para números de Reynolds superiores a 1.000.000. Por su parte, los
factores de fricción obtenidos con base en la validación de este diseño, utilizando la
ecuación de Darcy-Weisbach, se localizan en la zona de transición pero la mayoría de
estos tienden a ubicarse cerca al límite superior del FTHL, definido por el autor, para
tuberías de diámetro grande. Lo anterior, explica por qué el coeficiente de Hazen-Williams
sugerido en la normativa colombiana para este material, sobreestima las pérdidas por
fricción, dado que dicho coeficiente obliga a que el flujo en las tuberías en esta red sea
Flujo Turbulento Transicional, cuando en la realidad, para la rugosidad absoluta dada para
el CCP y los diámetros obtenidos, tiende al Flujo Turbulento Hidráulicamente Liso. Sin
embargo, es importante mencionar que para números de Reynolds superiores a
1.000.000, los factores de fricción tanto del diseño como la validación tienden a igualarse
para diámetros pequeños, como se puede ver en Gráfica 5.3-1.
124
• GRP
Con base en los resultados obtenidos en cuanto a la validación de los diseños de la red
Matriz de Medellín, utilizando el coeficiente de Hazen-Williams corregido y sugerido en la
normativa colombiana para este material, se puede ver en la Gráfica 4.4-4 que el diseño
obtenido a partir del coeficiente de Hazen-Williams corregido en este trabajo (CHw: 147)
tiende a igualar las presiones en esta red, cumpliendo con la presión mínima de diseño
(20 m.c.a), luego de realizar su respectiva validación con la ecuación de Darcy-Weisbach.
Por otro lado, posteriormente de validar el diseño obtenido a partir del coeficiente de
Hazen-Williams sugerido en la normativa colombiana (CHw: 150) para este material, se
obtuvo que este subestima las presiones en esta red y por lo tanto en algunos nudos de
ésta no se logra cumplir con la presión mínima de diseño.
Lo anterior se explica en la medida en que los números de Reynolds obtenidos se ubican

dentro del rango establecido para el coeficiente de Hazen-Williams corregido para el GRP,
donde se logra controlar en un 5% el error en las pérdidas por fricción (Ver Gráfica 4.4.5).
De otro lado, para el diseño calculado a partir del coeficiente de Hazen-Williams sugerido
en la normativa colombiana (Ver Gráfica 4.4-6), se puede observar que para diámetros
inferiores a 500 mm, este en su mayoría subestima las pérdidas por fricción (hasta un
14%), para números de Reynolds superiores a 300.000. Para entender el comportamiento
en la estimación de las pérdidas por fricción de estos coeficientes, se efectuó el mismo
análisis realizado para el material CCP, para establecer el tipo de flujo que se presenta en
los diseños mencionados anteriormente.
para este material GRP, se localizan en la zona de Flujo Turbulento Hidráulicamente Liso,
independiente de la magnitud del número de Reynolds, con base en el límite definido por
el autor para este tipo de flujo. De igual manera, los factores de fricción obtenidos en la
validación de este diseño, utilizando la ecuación de Darcy-Weisbach, se localizan en la
misma zona. Lo anterior, explica por qué el coeficiente de Hazen-Williams sugerido en
este trabajo controla las pérdidas por fricción en un 5%, dado que dicho coeficiente obliga
a que el flujo en las tuberías en esta red sea FTHL, flujo donde es posible utilizar dicho
125
coeficiente con buena precisión, con respecto a la ecuación de Darcy-Weisbach, según la

Gráfica 4.2-4.
Diagrama de Moody (Red Medellín, GRP, ks: 0.03 mm)

0.1
Factor de fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Hazen-Williams (CHw: 147) Comprobación Darcy-Weisbach (CHW: 147)
Hazen-Williams (CHw:150) Comprobación Darcy-Weisbach (CHw: 150)
Gráfica 5.3-2 Factores de Fricción obtenidos para el diseño de la red matriz de Medellín para el GRP,
específicamente para coeficientes de Hazen-Williams iguales a 147 y 150.
Hazen-Williams sugerido en la normativa para el GRP, se puede observar en la Gráfica
5.3-2, que se ubican en el FTHL. Por su parte, los factores de fricción obtenidos con base
en la validación de este diseño, utilizando la ecuación de Darcy-Weisbach, se localizan,
en su mayoría en el FTHL. Sin embargo, para números de Reynolds superiores a
1.000.000 y diámetros pequeños, estos factores tienden a ubicarse por fuera de este tipo
de flujo. Lo anterior, explica por qué el coeficiente de Hazen-Williams sugerido en la
normativa colombiana para este material, subestima las pérdidas por fricción para
diámetros pequeños, dado que dicho coeficiente obliga a que el flujo en las tuberías en
esta red sea FTHL, cuando en la realidad, para la rugosidad absoluta dada para el GRP,
el flujo es Turbulento Transicional para números Reynolds superiores a 1.000.000.
126
• Acero
A partir de los resultados obtenidos en cuanto a la validación de los diseños de la red

Matriz de Medellín, utilizando el coeficiente de Hazen-Williams corregido y sugerido en la
tiende a igualar las presiones en la red matriz de Medellín, cumpliendo con la presión
mínima de diseño de esta red (20 m.c.a). Por otro lado, el diseño obtenido a partir del
coeficiente de Hazen-Williams sugerido en la normativa colombiana (CHw: 120) para este
material, sobrestima las presiones en esta red y en algunos nudos no se logra cumplir con
la presión mínima de diseño, dado que son inferiores a 20 m.c.a o en el peor de los casos
negativas.
Lo anterior se explica en la medida en que el coeficiente de Hazen-Williams corregido

para el material GRP, como se puede apreciar en la Gráfica 4.4-8, logra controlar el error
en las pérdidas por fricción en un 5%. Sin embargo, para números de Reynolds
superiores a 200.000 (Ver Gráfica 4.4-8) se presenta subestimación de las pérdidas por
fricción hasta en un 35%. De otro lado, para el diseño calculado a partir del coeficiente de
Hazen-Williams sugerido en la normativa colombiana para este material (Ver Gráfica 4.4-
9), se puede observar que para diámetros inferiores a 1000 mm, este en su mayoría
subestima las pérdidas por fricción (hasta un 25%), para todos los diámetros obtenidos
para números de Reynolds superiores a 300.000. Para entender el comportamiento en la
estimación de las pérdidas por fricción de estos coeficientes, se efectuó el mismo análisis
realizado para el material CCP y GRP en el diagrama de Moody, para establecer el tipo
de flujo que se presenta en los diseños mencionados anteriormente.
127
Diagrama de Moody (Red Medellín, Acero, ks: 0.45 mm)

0.1
Factor de fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 5.3-3 Factores de Fricción obtenidos para el diseño de la red matriz de Medellín para el Acero,
Con base en la Gráfica 5.3-3, se puede ver que los factores hidráulicos obtenidos a partir
del diseño de esta red, utilizando la ecuación de Hazen-Williams y el coeficiente corregido
para este material, se localizan en la zona del Flujo Turbulento Transicional, para
números de Reynolds superiores a 200.000, cerca al límite superior del FTHL definido por
el autor. De otro lado, para este mismo rango de números de Reynolds, los factores de
fricción obtenidos en la validación de este diseño, utilizando la ecuación de Darcy-
Weisbach, se localizan en la zona de transición y tienden a ubicarse en el límite mínimo
del FTHR. Lo anterior, explica por qué el coeficiente de Hazen-Williams sugerido en este
trabajo para este material subestima las pérdidas por fricción en este rango del número de
Reynolds, dado que dicho coeficiente obliga a que el flujo en las tuberías en esta red sea
cercano al FTHL, cuando en la realidad, para la rugosidad absoluta dada para el Acero, el
flujo es Turbulento Transicional para los diámetros pequeños. Es importante resaltar que,
en el rango de número de Reynolds definido para utilizar dicho coeficiente, los factores de
fricción tanto del diseño como la validación del mismo, tienden a igualarse lo cual es
consistente con la Gráfica 4.4-8.
128
Hazen-Williams sugerido en la normativa para el Acero, se puede observar en la Gráfica
5.3-3, que estos se ubican en la zona de Flujo Turbulento Transicional y tienden al límite
mínimo del FTHR para números de Reynolds superiores a 200.000. Por su parte, los
factores de fricción obtenidos con base en la validación de este diseño, utilizando la
ecuación de Darcy-Weisbach, se localizan en la zona de transición pero la mayoría de
estos tienden a ubicarse cerca al límite inferior del FTHR o en su defecto en este tipo de
flujo, para tuberías de diámetro pequeño. Lo anterior, explica por qué el coeficiente de
Hazen-Williams sugerido en la normativa colombiana para este material, subestima las
pérdidas por fricción, dado que dicho coeficiente obliga a que el flujo en las tuberías en
esta red sea Flujo Turbulento Transicional, cuando en la realidad, para la rugosidad
absoluta dada para el Acero y los diámetros obtenidos, tiende al Flujo Turbulento
Hidráulicamente Rugoso.
5.3.2 Red Hanoi
• CCP
Matriz de Hanoi, utilizando el coeficiente de Hazen-Williams corregido y sugerido en la
tiende a igualar las presiones en la red matriz de Hanoi, pero solamente en un nodo de
esta no se cumple con la presión mínima de diseño de esta red (20 m.c.a). Por otro lado,
el diseño obtenido a partir del coeficiente de Hazen-Williams sugerido en la normativa
colombiana (CHw: 120) para este material, subestima las presiones de la red y cumplen
con la presión mínima de diseño para todos los nudos de la esta red.
Lo anterior se explica dado que el coeficiente de Hazen-Williams corregido para este

material (Ver Gráfica 4.4-11), en su mayoría subestiman las pérdidas por fricción (hasta
un 15%) para números de Reynolds superiores a 400.000. De otro lado, para el diseño
129
calculado a partir del coeficiente de Hazen-Williams sugerido en la normativa colombiana

(Ver Gráfica 4.4-12), se puede observar que para todos los diámetros obtenidos se
sobrestiman las pérdidas por fricción en un rango de 15 a 30%.
Para entender los errores en las pérdidas por fricción, es necesario analizar en detalle los
tipos de flujo que se presentan en la red Hanoi para el material CCP, gráficando los
factores de fricción obtenidos en el diagrama de Moody, llevando a cabo el mismo análisis
efectuado para la red matriz de Medellín.
Diagrama de Moody (Red Hanoi, CCP, ks: 0.12 mm)

0.1
Factor de fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Hazen-Williams (CHW: 144) Comprobación Darcy-Weisbach (CHW: 144)
Hazen-Williams (CHW:120) Comprobación Darcy-Weisbach (CHW: 120)
Gráfica 5.3-4 Factores de Fricción obtenidos para el diseño de la red Hanoi para el CCP,
para este material CCP, se localizan en la zona de Flujo Turbulento Hidráulicamente Liso,
el autor para este tipo de flujo. De otro lado, los factores de fricción obtenidos en la
zona de transición. Lo anterior, explica por qué el coeficiente de Hazen-Williams sugerido
en este trabajo para este material subestima las pérdidas por fricción, dado que dicho
130
coeficiente obliga a que el flujo en las tuberías en esta red sea FTHL, cuando en la
realidad, para la rugosidad absoluta dada para el CCP, el flujo es Turbulento Transicional
para los diámetros obtenidos en este diseño.
red Hanoi, utilizando la ecuación de Hazen-Williams y el coeficiente de Hazen-Williams
sugerido en la normativa para el CCP, se puede observar en la Gráfica 5.3-4, que estos
se ubican en la zona de Flujo Turbulento Transicional y tienden al límite mínimo del FTHR
para números de Reynolds superiores a 400.000. Por su parte, los factores de fricción
obtenidos con base en la validación de este diseño, utilizando la ecuación de Darcy-
Weisbach, se localizan la zona de transición pero la mayoría de estos tienden a ubicarse
cerca al límite superior del FTHL, definido por el autor, para tuberías de diámetro grande.
Lo anterior, explica por qué el coeficiente de Hazen-Williams sugerido en la normativa
colombiana para este material, sobreestima las pérdidas por fricción, dado que dicho
coeficiente obliga a que el flujo en las tuberías en esta red sea Flujo Turbulento
Transicional, cuando en la realidad, para la rugosidad absoluta dada para el CCP y los
diámetros obtenidos, tiende al Flujo Turbulento Hidráulicamente Liso. Sin embargo, es
importante mencionar que para números de Reynolds superiores a 1.000.000, los factores
de fricción tanto del diseño como la validación tienden a igualarse para diámetros
pequeños, como se puede ver en Gráfica 5.3-1.
• GRP
Matriz de Hanoi, utilizando el coeficiente de Hazen-Williams corregido y sugerido en la
tiende a igualar las presiones en la red matriz de Hanoi, cumpliendo con la presión
mínima de diseño de esta red (20 m.c.a). A diferencia de este coeficiente, el diseño
obtenido a partir del coeficiente de Hazen-Williams sugerido en la normativa colombiana
(CHw: 150) no cumple con dicha presión en un nodo de esta red.
131

material controla en un 5% el error en las pérdidas por fricción, como se puede apreciar
en la Gráfica 4.4.14, en el rango del número de Reynolds establecido para dicho
coeficiente. De otro lado, para el diseño calculado a partir del coeficiente de Hazen-
Williams sugerido en la normativa colombiana (Ver Gráfica 4.4-15), se puede observar
que para diámetros inferiores a 500 mm, este en su mayoría subestima las pérdidas por
fricción (hasta un 5%), para números de Reynolds entre 100.000 y 1.000.000. Para
entender el comportamiento en la estimación de las pérdidas por fricción de estos
coeficientes, se efectuó el mismo análisis realizado para el material CCP, para establecer
el tipo de flujo que se presenta en los diseños mencionados anteriormente.
Diagrama de Moody (Red Hanoi, GRP, ks: 0.03 mm)

0.1
Factor de fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 5.3-5 Factores de Fricción obtenidos para el diseño de la red Hanoi para el GRP,
para este material GRP, se localizan en la zona de Flujo Turbulento Hidráulicamente Liso,
el autor para este tipo de flujo. De igual manera, los factores de fricción obtenidos en la
misma zona. Lo anterior, explica por qué el coeficiente de Hazen-Williams sugerido en
132
este trabajo controla las pérdidas por fricción en un 5%, dado que dicho coeficiente obliga
a que el flujo en las tuberías en esta red sea FTHL, flujo donde es posible utilizar dicho
coeficiente con buena precisión, con respecto a la ecuación de Darcy-Weisbach, según la
Gráfica 4.2-4.
sugerido en la normativa para el GRP, se puede observar en la Gráfica 5.3-5, que estos
se ubican en el FTHL. Por su parte, los factores de fricción obtenidos con base en la
validación de este diseño, utilizando la ecuación de Darcy-Weisbach, se localizan, en su
mayoría en el FTHL. Sin embargo, para números de Reynolds superiores a 1.000.000 y
diámetros pequeños, estos factores tienden a ubicarse por fuera de este tipo de flujo. Lo
anterior, explica por qué el coeficiente de Hazen-Williams sugerido en la normativa
colombiana para este material, subestima las pérdidas por fricción para diámetros
pequeños, dado que dicho coeficiente obliga a que el flujo en las tuberías en esta red sea
FTHL, cuando en la realidad, para la rugosidad absoluta dada para el GRP, el flujo es
Turbulento Transicional para números Reynolds superiores a 1.000.000.
• Acero
En relación con los resultados obtenidos en cuanto a la validación de los diseños de la

red Hanoi, se puede ver en la Gráfica 4.4-16 que el diseño obtenido a partir del coeficiente
de Hazen-Williams corregido en este trabajo (CHw: 139) tiende a igualar o sobrestimar las
presiones en la red, cumpliendo con la presión mínima de diseño de esta red (20 m.c.a).
Por otro lado, el diseño obtenido a partir del coeficiente de Hazen-Williams sugerido en la
normativa colombiana (CHw: 120) para este material, iguala las presiones de la red y en
todos los nodos cumple con la presión mínima de diseño.

material (Ver Gráfica 4.4-17), el error en las pérdidas por fricción se ubican en el rango del
número de Reynolds establecido para dicho coeficiente, donde los errores en las pérdidas
no superan el 5%, dado que para diámetros iguales 1000 mm se manejan velocidades de
hasta 0.5 m/s, límite de aplicabilidad de este coeficiente. De otro lado, para el diseño
133
calculado a partir del coeficiente de Hazen-Williams sugerido en la normativa colombiana

(Ver Gráfica 4.4-18), se puede observar que el error en las pérdidas no supera el 5% para
números de Reynolds superiores a 100.000, teniendo en cuenta las velocidades que se
presentan para los diámetros obtenidos. Para entender el comportamiento en la
estimación de las pérdidas por fricción de estos coeficientes, se efectuó el mismo análisis
realizado para el material CCP y GRP, para establecer el tipo de flujo que se presenta en
los diseños mencionados anteriormente.
Diagrama de Moody (Red Hanoi, Acero, ks: 0.45 mm)

0.1
Factor de fricción
0.01
1.00E+03 1.00E+04 1.00E+05 1.00E+06 1.00E+07
Número de Reynolds
Gráfica 5.3-6 Factores de Fricción obtenidos para el diseño de la red Hanoi para el Acero,
Con base en la Gráfica 5.3-6, se puede ver que los factores hidráulicos obtenidos a partir
del diseño de esta red, utilizando la ecuación de Hazen-Williams y el coeficiente corregido
para este material, se localizan en la zona del FTHL, para números de Reynolds entre
25000 y 250000, cerca al límite superior del FTHL definido por el autor. De igual manera,
para este mismo rango de números de Reynolds, los factores de fricción obtenidos en la
validación de este diseño, utilizando la ecuación de Darcy-Weisbach, se localizan en esta
zona. Sin embargo, algunos de ellos tienden a ubicarse fuera de esta zona para diámetros
pequeños. Lo anterior, explica por qué el coeficiente de Hazen-Williams sugerido en este
trabajo para este material iguala o subestima las pérdidas por fricción en este rango del
número de Reynolds. Es importante resaltar que, en el rango de número de Reynolds
134
definido para utilizar dicho coeficiente, los factores de fricción tanto del diseño como la
validación del mismo, tienden a igualarse lo cual es consistente con la Gráfica 4.4-8.
sugerido en la normativa para el Acero, se puede observar en la Gráfica 5.3-6, que estos
se ubican en la zona de Flujo Turbulento Transicional y tienden al límite mínimo del FTHR
para números de Reynolds superiores a 200.000. Por su parte, los factores de fricción
obtenidos con base en la validación de este diseño, utilizando la ecuación de Darcy-
Weisbach, se localizan en la zona de transición pero la mayoría de estos tienden a
ubicarse cerca al límite inferior del FTHR, para tuberías de diámetro pequeño. Lo anterior,
explica por qué el coeficiente de Hazen-Williams sugerido en la normativa colombiana
para este material, las pérdidas por fricción no superan el 5%.
135
Capítulo 6. CONCLUSIONES
• Con base en los resultados obtenidos en relación con el efecto del uso de la
ecuación de Hazen-Williams en las redes de estudio, se puede concluir, en
principio, que los coeficientes de Hazen-Williams sugeridos en la normativa
colombiana, no son correctos para el Acero y CCP, dado que estos valores
tienden a sobrestimar las pérdidas por fricción y no tienen en cuenta la variabilidad
que presenta dicho coeficiente en el comportamiento hidráulico de estas redes.
• A pesar de que los factores de fricción obtenidos por medio de la ecuación de
Hazen-Williams se ubican en la zona de transición en el diagrama de Moody para
el Acero y GRP, como lo definió en su momento Diskin (1960) y Liou como límite
de aplicabilidad de esta ecuación, no es un condición suficiente para garantizar la
aplicabilidad de esta, dado que esta no tienen en cuenta la variación de la
rugosidad relativa que se puede presentar en el diseño de redes matrices, para
números de Reynolds superiores a 1.000.000. Lo anterior, se ve reflejado en la
validación de los diseños de las redes de estudio.
• Para la red matriz analizada por Bombardelli y García (2003), se logró establecer
que el coeficiente de Hazen-Williams utilizado en el diseño original de esta red no
es adecuado para las rugosidades absolutas presentes en la línea principal de
esta red, dado que este tiende a subestimar las pérdidas por fricción, según la
ecuación de Darcy-Weisbach, y por lo tanto no permite predecir con certeza el
comportamiento hidráulico de esta red.
• A partir de los resultados obtenidos para la corrección de los coeficientes de
Hazen-Williams, en especial para el Acero y CCP, se puede concluir que se logró
tener en cuenta la variabilidad de estos y controlar hasta en un 5% las pérdidas
por fricción para un rango específico de números de Reynolds y velocidades, para
tuberías de gran diámetro. Esto se reflejó en la validación de los diseños de las
redes de estudio, donde estos coeficientes permiten cumplir con la presión mínima
de diseño. Pese a ello, para diseñar una red matriz específica con estos
materiales, se sugiere verificar si los factores de fricción obtenidos después de
136
realizar la comprobación de los mismos con la ecuación de Darcy-Weisbach se

ubican en el FTHL, según el límite definido por el autor.
• En relación con el coeficiente de Hazen-Williams sugerido en la normativa

colombiana para el material GRP, se sugiere modificar dicho coeficiente a 147,
dado que este evita subestimar, en un 5%, las pérdidas por fricción para diámetros
superiores a 500 mm y de esa manera permite cumplir con la presión mínima de
diseño. Esto se refleja en los resultados obtenidos en cuanto a la validación de los
diseños de las redes de estudio.
• Teniendo en cuenta el creciente aumento de la demanda de agua en las ciudades,

la continúa operación y las características de las tuberías de las redes matrices,
así como los efectos que tiene la variabilidad del coeficiente de Hazen-Williams en
la estimación de las pérdidas por fricción en los materiales de estudio, se sugiere
utilizar en el diseño de las redes matrices la ecuación de Darcy-Weisbach, dado
que esta es físicamente basada, la cual se puede utilizar para cualquier fluido,
rango de velocidades, diámetros y longitudes de tubería.
• Para corregir los coeficientes de Hazen-Williams sugeridos en la normativa

colombiana fue necesario modificar dichos coeficientes, para diámetros superiores
a 500 mm. Por lo anterior, para futuras investigaciones debería evaluarse el efecto
del uso de la ecuación de Hazen-Williams en términos de costos, teniendo en
cuenta esta limitación.
• Para los materiales analizados en la red matriz de Medellín, se puede concluir que
la ecuación de Hazen-Williams no se debe utilizar cuando el flujo sea laminar dado
que en este tipo de flujo se presenta una alta variabilidad del coeficiente de Hazen-
Williams para el Acero, CCP y GRP, lo cual se refleja en los errores obtenidos en
cuanto a las pérdidas por fricción.
137
BIBLIOGRAFÍA
1. Barlow, J.F., and Markland, E. (1975). "Converting the Hazen-Williams equation to the
Colebrook function". Water Power Dam Construction , 331-334.
2. Brown, G. O. (2004). "The History of the Darcy-Weisbach equation for pipe flow
resistance". Recuperado Marzo 4, 2010, de ascelibrary: http://www.ascelibrary.org
3. Colebrook C.F. (1939). "Turbulent Flow in pipes with particular reference to the
transition region between the smooth and rough pipes laws". Proc. Institution Civil
Engrs , 393-422.
4. Committee of Pipeline Design (1975). Pipeline design for water and wastewater. New
York: ASCE.
5. Diskin, M. (1960). The limits of applicability of the Hazen-Williams formula. Houille

Blanche , 720-723.
6. EPM. (2009). Normas de diseño de sistemas de acueducto. Medellín.
7. Fabian Bombardelli & García M.(2003). "Hydraulic Design of Large-Diameter Pipes".

Recuperado 13 Marzo, 2010, de ascelibrary: http//: www.ascelibrary.com
8. Flechas R.A. (2010). Efecto del uso de la Ecuación de Darcy-Weisbach vs la Ecuación

de Hazen-Williams sobre los costos de diseños optimizados en acueductos, teniendo
en cuenta la rugosidad de las tuberías, linea del gradiente hidráulico y temperatura.
Bogotá: Trabajo de Grado: Universidad de los Andes.
9. Gardner Williams, A. H. (1905). Hydraulic tables: showing the loss of head due to the
friction of water flowing in pipes, aqueducts, sewers, etc (primera edicion ed.). New
York: Jhon Wiley & Sons.
10. Guzman, J. M. (1999). Calibraciones realizadas en los acueductos de algunas

ciudades de Colombia. Bogotá D.C.: trabajo de Grado: Universidad de los Andes.
11. Jain, A. K., Mohan, D.M, and Khanna, P. (1978). "Modified Hazen-Williams
formula". J. Environ. Eng. Div (ASCE) , 137-146.
12. Kamand, F. Z. (1988). "Hidraulic Friction Factors for Pipe Flow". ASCE , 311-323.
13. Liou, C. P. (1998). Limitations and proper use of the Hazen-Williams equation.
Retrieved Marzo 5, 2010, from ascelibrary: http://www.ascelibrary.com
14. Moody, L. (1944). Friction factors for pipe flow. Trans. ASME. 671-684.
138
15. Saldarriaga Valderrama, J. G. (2007). Hidráulica de tuberías: abastecimiento de

agua, redes, riegos . Bogotá, D.C.: Alfaomega: Uniandes.
16. Vennard, J. K. (1958). Elementary Fluid Mechanics (3rd ed.). New York: Wiley.
17. Walski, T. (2003). Discussion of “Hydraulic Design of Large-Diameter Pipes” by

Fabian A. Bombardelli and Marcelo H. García. Retrieved Marzo 5, 2010, from ascelibrary:
http://www.ascelibrary.com
139
ANEXOS
En esta parte del trabajo se ilustra la topología de las distintas redes de estudio que
fueron seleccionadas para analizar el efecto del uso de la ecuación de Hazen-Williams vs.
la ecuación de Darcy-Weisbach en tuberías de gran diámetro. Asi mismo, las hojas
electrónicas que contienen el programa desarrollado en esta tesis para corregir los
coeficientes de Hazen-Williams sugeridos en la normativa colombiana.
Manantiales
Villa Hermosa
Ayurá
Anexo 1. Red matriz Medellín.
140
Anexo 2. Red matriz Hanoi, incluyendo la magnitud de los diámetros presentes en esta red, en milímetros.
141
de gran diámetro (Redes matrices).
Anexo 3. Interfaz programa desarrollado para corregir los coeficientes de Hazen-Williams, hoja electrónica “Programa”.
142
de gran diámetro (Redes matrices)
Anexo 4. Interfaz programa desarrollado para corregir los coeficientes de Hazen-Williams, hoja electrónica “Resultados”.
143
Anexo 5. Interfaz programa desarrollado para corregir los coeficientes de Hazen-Williams, hoja electrónica “Errores pérdidas por fricción”.
144
Anexo 6. Interfaz programa desarrollado para corregir los coeficientes de Hazen-Williams, hoja electrónica “Diagrama de Moody”.
145
Anexo 7. Código fuente programa desarrollado para corregir los coeficientes de Hazen-Williams.
146

Ecuacion Hazem Williams PDF

Cargado por

Copyright:

Formatos disponibles

Ecuacion Hazem Williams PDF

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Ecuacion Hazem Williams PDF

Cargado por

Copyright:

Formatos disponibles

Universidad de los Andes

Centro de Investigaciones en Acueductos y

Efecto del uso de la ecuación de Hazen-Williams vs la

Rafael Alejandro Flechas Hernández

Bogotá, Diciembre de 2011

El autor agradece la gran colaboración brindada por el ingeniero Juan Guillermo

A Dios, a mi madre, mi padre, mis hermanos y mi abuelita quienes siempre me han

brindado su apoyo incondicional a lo largo de mi vida y son la razón de mi existir.

En el diseño y análisis de redes de distribución de Agua potable (RDAP), las ecuaciones

Capítulo 5. ANÁLISIS DE RESULTADOS ................................................................................. 116

Figura 2.1-1 Diagrama de Moody (Moody, 1940) ............................................................. 13

Gráfica 4.3-10 Diagrama de dispersión presiones obtenidas en el diseño de la red Hanoi

- Análisis de redes de distribución de agua potable: Procedimiento que consiste en

- Coeficiente de Hazen-Williams (CHw): Coeficiente que describe la rugosidad de las

- Diagrama de Moody: Es un diagrama que permite encontrar de forma gráfica el factor

- Diseño de Redes de Distribución de Agua potable: Procedimiento que consiste en

- Ecuación de Colebrook-White: Ecuación logarítmica e implícita que permite calcular el

- Ecuación de Darcy-Weisbach: Ecuación físicamente basada que permite determinar

- Ecuación de Hazen-Williams: Ecuación de origen empírico que describe el cálculo de

- EPANET: Es un programa desarrollado por la Agencia de Protección Ambiental de los

- Factor de fricción (f): Es un parámetro característico de la Ecuación de Darcy-

- Fluido: Sustancia que no puede resistir esfuerzo cortante.

- Flujo: Es el movimiento de las partículas de un fluido con respecto a un sistema inercial.

- Flujo laminar: Es el tipo de flujo donde se presentan bajas velocidades y viscosidades

- Flujo turbulento hidráulicamente rugoso (FTHR): Flujo turbulento en una tubería

- Método del Gradiente: Es un método de cálculo de redes de distribución de agua

- REDES: Es un programa desarrollado por el Centro de Investigaciones de Acueductos y

- Redes de distribución de agua potable (RDAP): Las redes de distribución de agua

- Redes matrices: Es un componente de las redes de distribución de agua potable, las

- Rugosidad relativa (ks/d): Relación entre la rugosidad absoluta de las tuberías y el

Dicho problema ocasionó la aparición de ecuaciones empíricas más sencillas de utilizar

Por lo anterior, el uso de la ecuación de Hazen-Williams frente a la ecuación de Darcy-

Por ello, este documento está estructurado de la siguiente manera:

El Capítulo 1 contiene los objetivos de la tesis, explicando lo que se buscó con el

El Capítulo 4 presenta los resultados obtenidos en el análisis de estas ecuaciones en las

En el Capítulo 6 se muestran algunas conclusiones de esta tesis a tener en cuenta en el

• Determinar el efecto del uso de ecuaciones empíricas vs ecuaciones físicamente

• Establecer el efecto del uso de la ecuación de Hazen-Williams vs la ecuación de

Capítulo 2. MARCO TEÓRICO

2.1 Ecuaciones de Resistencia Fluida

z= Cota del fondo de la tubería.

Por lo anterior, en el diseño y análisis de redes de distribución de agua potable (RDAP),

2.1.1 Ecuación de Darcy-Weisbach

En 1845, Julius Weisbach (1806-1871) ingeniero, de nacionalidad sajona, propuso una

hf= Pérdidas por fricción.

Esta es la forma que actualmente se conoce de la ecuación de Darcy-Weisbach; sin

donde a y b eran coeficientes empíricos.

• Diámetros desde 0.0122 hasta 0.5 m.

Donde α, β, α´ y β´ son coeficientes de fricción. Según Darcy, el primer término de esta

donde α´ y β´ pueden ser coeficientes de fricción distintos a los valores de tuberías de

Posteriormente, J.T Fanning (1837-1911), en 1877, fue el primero en combinar los

Teniendo en cuenta lo anterior, la ecuación de Darcy-Weisbach fue deducida por Henry

que se define de la siguiente manera.

donde V es la velocidad, D el diámetro y 𝜈 la viscosidad cinemática del fluido, parámetro

• Flujo laminar Re<2000.

De todos los parámetros de la ecuación de Darcy-Weisbach, el factor de fricción es el más

Para números de Reynolds superiores a 4000, en el diagrama de Moody se pueden ver

0.305𝛿´ < 𝑘𝑠 < 6.1𝛿´

ks= Rugosidad de la tubería.