Manual Del Sistema TransDocencia 2.1

Manual del sistema de simulación
de transitorios hidráulicos en
conductos a presión TRANS, versión
docencia
Libia Georgina Carmona Paredes

Gabriel Castillo Hernández
Rafael Carmona Paredes
Informe final del proyecto PE104017

Abril, 2019
MANUAL DEL SISTEMA DE SIMULACIÓN DE TRANSITORIOS HIDRÁULICOS EN CONDUCTOS A PRESIÓN TRANS, VERSIÓN DOCENCIA
CONTENIDO
Contenido
NOMENCLATURA iii
CAPÍTULO 1. Intoducción 1
1.1 Introducción 1
CAPÍTULO 2. Transitorios hidráulicos en sistemas de bombeo 4

2.1 Introducción 4
2.2 Oscilación de masa 5
2.2.1 Modelo matemático de oscilación de masa 6
2.3 Golpe de ariete 17
2.3.1 Modelo matemático de golpe de ariete 18
CAPÍTULO 3 Sistema TransDocencia2-0 31

3.1 Introducción 31
3.2 Descripción del sistema TransDocencia2-0 31
3.3 Pantalla principal del sistema TransDocencia2-0 34
3.4 Creación de un proyecto 40
3.5 Creación de una propuesta 42
3.6 Inspector de propiedades 44
3.7 Simulación y visualización de resultados 46
CAPÍTULO 4 Elementos de un sistema de bombeo que no 50

producen transitorios
4.2 Tuberías 51
4.2.1 Simulación en el sistema TransDocencia2-0 52
4.3 Cambio de tubería o tramo 53
4.4 Cámara de aire 55
4.4.1 Modelo matemático de cámara de aire 58
4.4.2 Predimensionamiento de cámaras de aire 62
4.5 Tanques de carga constante 73
4.5.1 Tanque en la frontera aguas arriba 74
i
CONTENIDO

4.5.3 Tanque en la frontera aguas abajo 75
4.6 Tanque unidireccional 77
4.6.1 Modelo matemático de tanque unidireccional 79
4.6.2 Predimensionamiento de tanques 81
unidireccionales
4.7 Torre de oscilación 84
4.7.1 Modelo matemático de torre de oscilación 86
4.7.2 Predimensionamiento de torres de oscilación 87
4.8 Válvula de admisión y expulsión de aire 89
4.8.1 Modelo matemático de válvulas de admisión 90
y expulsión de aire
4.8.3 Dimensionamiento de válvulas de admisión 94
y expulsión de aire
CAPÍTULO 5 Elementos de sistemas de bombeo que generan 96

transitorios
5.2 Válvulas de movimiento programado o controlado 96
5.2.1 Modelo matemático de válvulas de 98
movimiento controlado
5.3 Disparo de equipos de bombeo 104
5.3.1 Modelo matemático 112
5.3.2.1 Bombas en la frontera aguas arriba 116
5.3.2.2 Bombas a lo largo de la conducción 120
CAPÍTULO 6 Ejemplos 122

6.1 Ejemplo 1 122
6.1.1 Archivos de resultados 138
6.2 Ejemplo 2 145
BIBLIOGRAFÍA 149
ii
NOMENCLATURA
Nomenclatura
𝐴 sección transversal tuberías, área transversal de tanques o torres
𝑎 celeridad
𝑐𝑎 constante del método de las características
𝑐𝑛 parámetro de la característica negativa
𝑐𝑝 parámetro de la característica positiva
𝐷 diámetro de la tubería
𝐸𝑉 módulo de elasticidad volumétrica del agua
𝐸𝑡 módulo de elasticidad del material de la tubería
𝑒 espesor de la tubería
𝐹𝐴 valor de la variable 𝐹 en la posición 𝑥 − Δ𝑥 en el tiempo 𝑡
𝐹𝐵 valor de la variable 𝐹 en la posición 𝑥 + Δ𝑥 en el tiempo 𝑡
𝐹𝑃 valor de la variable 𝐹 en la posición 𝑥 en el tiempo 𝑡 + ∆𝑡
𝑓 factor de fricción de Darcy-Weisbach de la tubería
𝑓ℎ función característica de carga, adimensional
𝑓𝛽 función característica de par o momento, adimensional
𝑔 aceleración de la gravedad
𝐻 carga piezométrica
ℎ carga, adimensional
𝐼 término inercial
𝐾 constante de pérdida de carga
𝑙 longitud de la tubería
𝑁 velocidad de la bomba
𝑁𝑅 velocidad de diseño de la bomba
𝑛 velocidad de la bomba, adimensional
𝑃 potencia mecánica en la operación con máxima eficiencia
𝑃(𝜏) polinomio de sexto grado
𝑝 presión absoluta del aire en la cámara
𝑝𝑎 presión atmosférica
𝑄 gasto en la tubería
𝑞 gasto, adimensional
𝑅 radio de la tubería
𝑅𝑒 número de Reynolds
iii
NOMENCLATURA
SLA superficie libre del agua

𝑇 periodo de la onda de presión, torque, momento angular o par aplicado al eje de la
bomba
𝑇𝑅 torque, momento angular o par aplicado al eje de la bomba en operación con máxima
eficiencia
𝑡 tiempo
𝑉 volumen de aire
𝑊𝑅 2 momento angular de inercia total (bomba, motor y agua que pasa por la bomba, en
función del radio y la masa)
𝑍 elevación del agua en cámaras de aire, tirante de agua, altura de tanques de oscilación
𝛼 coeficiente politrópico del aire
𝛽 torque, momento angular o par, adimensional
𝛾 peso específico del agua
∆𝑡 intervalo de integración
𝜀 promedio de la rugosidad de la tubería
𝜂 eficiencia
𝜃 ángulo en grados
𝑣 velocidad del flujo
𝜗 viscosidad del fluido
𝜌 densidad del agua
𝜏 variable adimensional
𝜔 velocidad angular de la bomba
iv
INTRODUCCIÓN
1
Introducción
1.1 Introducción
TransDocencia2.0 es una herramienta de apoyo para analizar el comportamiento de

sistemas de bombeo durante la operación en flujo transitorio, de ninguna manera sustituye
la lectura de libros especializados en el tema. El objetivo de su desarrollo es contribuir a
formación de profesionistas y especialistas con amplio conocimiento en el manejo de flujo
transitorio en conducciones de agua a presión, en particular en sistemas de bombeo o
acueductos.
El software que se presenta, TransDocencia2-0, es la versión docencia del Sistema de

Simulación de Transitorios Hidráulicos en Conductos a Presión: TRANS desarrollado en el
Instituto de Ingeniería de la UNAM.
Al fenómeno que ocurre durante el intervalo de tiempo que le toma a un sistema hidráulico
pasar de una condición de operación en flujo permanente que fue perturbado, a otra
condición de operación en flujo permanente se le llama transitorio hidráulico.
Algunas perturbaciones ocurren súbitamente como consecuencia de eventos imprevistos

otras, son consecuencia de las maniobras que deben realizarse de manera programada y
1
INTRODUCCIÓN
controlada para cambiar el modo de operación del sistema, como arranques y paros de
equipos de bombeo o apertura y cierre de válvulas para modificar caudales. Durante el flujo
transitorio que se desarrolla después haber ocurrido un evento imprevisto se presentan las
máximas variaciones en la presión interna de las tuberías, en la velocidad del agua y en la
velocidad de los elementos rodantes de bombas y turbinas, a éste se le llama transitorio
rápido o golpe de ariete.
Cambios súbitos en el consumo de energía eléctrica en la red alimentada por un sistema

hidroeléctrico y fallas en el suministro de energía eléctrica en plantas de bombeo, generan
golpes de ariete. Un golpe de ariete puede causar daños con costos muy altos, tanto
económicos como sociales. Los primeros, debido a las reparaciones que se tuvieran que
realizar en la infraestructura comunitaria y los segundos, por la suspensión del servicio de
entrega de agua o de generación de energía eléctrica.
El fenómeno del golpe de ariete se analiza a partir de las ecuaciones dinámica y de

continuidad, las que en muchas ocasiones apenas son presentadas en los cursos de
licenciatura y sólo se dan soluciones simplificadas. Para diseñar sistemas hidráulicos con
conductos a presión, no es suficiente considerar las soluciones que se obtienen con los
métodos simplificados, debido a ello el desarrollo del presente trabajo, el cual está
estructurado de la siguiente manera.
En el capítulo 2 se presentan los modelos de los transitorios por oscilación de masa y por
golpe de ariete, debido al tipo de ecuaciones que los determinan, ambos modelos se
resuelven a través de métodos numéricos. La solución del primero se plantea con un
esquema de integración de Euler y el segundo se resuelve utilizando el método de las
características. En este mismo capítulo se hace una descripción espacio-temporal de la
variación de la presión y el gasto en el interior de las tuberías durante un transitorio
producido por golpe de ariete.
En el capítulo 3 se describe el manejo del sistema TransDocencia2-0, formado por dos

sistemas de cómputo. El primero, es la interface de comunicación tipo CAD entre el usuario
y el sistema TransMoMaDoc2-0.exe, en este segundo se resuelven los modelos
matemáticos del flujo transitorio y de los elementos hidráulicos que integran los
acueductos. Al ejecutar el programa TransDocencia2-0.exe se accede a la interface de
comunicación, cuando se pide realizar la simulación de un transitorio se ejecuta
automáticamente TransMoMaDoc2-0.exe. Ambos sistemas están desarrollados en la
plataforma Delphi XE2.
En el capítulo 4 se presentan los elementos y estructuras hidráulicas que pueden

modelarse en el sistema TransDocencia2-0 que no generan transitorios, los modelos
matemáticos de cada uno de ellos, cómo es su operación, procedimientos para obtener
sus dimensiones preliminares y cómo se simulan en el sistema TransDocencia2-0.
2
INTRODUCCIÓN
En el capítulo 5 se presentan los elementos que tienen la capacidad de producir transitorios

hidráulicos (bombas y válvulas), sus modelos matemáticos y cómo se simulan en el sistema
TransDocencia2-0.
En el capítulo 6 se presentan dos ejemplos, a través de ellos, se muestra cómo se usa el

sistema TransDocencia2-0 para describir un sistema de bombeo, de qué manera se
obtiene la información para simular la operación de equipos de bombeo, se determinan las
dimensiones preliminares de algunos dispositivos de control de transitorios, se indica cómo
interpretar los resultados de las simulaciones y su uso para optimizar las dimensiones de
los dispositivos de control de transitorios.
3
TRANSITORIOS HIDRÁULICOS EN SISTEMAS DE BOMBEO
2
Transitorios hidráulicos en
sistemas de bombeo
2.1 Introducción
En los sistemas de bombeo se presentan transitorios hidráulicos cuando se realizan las
maniobras necesarias para cambiar su forma de operación. Dependiendo de las
estructuras hidráulicas que forman parte de un sistema, además de los transitorios por
golpe de ariete pueden presentarse transitorios por oscilación de masa. En comparación
de la rapidez con la que se amortiguan los transitorios producidos por golpe de ariete, a los
transitorios por oscilación de masa se les conoce como transitorios lentos.
El estudio numérico de estos dos tipos de transitorios se realiza con enfoques distintos,
para los de golpe de ariete el modelo de estudio está basado en la solución numérica de
las ecuaciones diferenciales hiperbólicas de un fenómeno ondulatorio en un medio elástico,
se considera el fluido compresible y el material elástico con el que está construida la
tubería. El análisis de los transitorios por oscilación de masa se basa en la solución de las
ecuaciones de continuidad y de energía en una tubería que une dos tanques abiertos (fluido
incompresible y tubería rígida), la ecuación de la energía no es lineal en la velocidad del
flujo, por lo que la solución de éstas ecuaciones debe plantearse también en forma
numérica, comúnmente basado en un esquema de integración de Euler.
4
2.2 Oscilación de masa

La determinación del comportamiento transitorio por oscilación de masa es de gran
importancia en el análisis del funcionamiento de un sistema de conducción, a través de él
se determinan algunas características de su obra civil como el uso de vasos de regulación,
el diámetro de las torres de oscilación, sus alturas para que no derramen agua o en caso
de que se permita se construyan las obras civiles necesarios para encausar el agua
derramada. En la figura 2.1 se muestra el tramo de conducción de un sistema de bombeo
en el que después de un cambio en la operación del sistema se presentará un transitorio
por oscilación de masa entre el tanque 1 y el tanque 2. En la figura 2.2 se muestra un
tanque de oscilación derramando agua.
Tanque 1
Tanque 2
𝐻1
𝐻2
𝑄0
Línea de referencia
Figura 2.1 Después de un cambio de operación del sistema se presentará un transitorio

de oscilación de masa entre los tanques o torres con cargas 𝐻1 y 𝐻2
5
Figura 2.2 Torre de oscilación con agua derramando (Grupo de Hidromecánica, IIUNAM)
2.2.1 Modelo matemático de oscilación de masa. El modelo matemático está

formado por las ecuaciones que representan la dinámica del agua en la tubería y la de la
superficie libre del agua en las torres de oscilación.
Dinámica del agua en la tubería. El estudio de los transitorios de oscilación de masa se

realiza considerando el agua incompresible y la tubería rígida, bajo estas condiciones, la
dinámica del agua en la tubería se describe con la ecuación de la energía, ecuación 2.1.
𝑙 𝑑𝑄 𝑓𝑙 𝑄|𝑄|
= 𝐻1 − 𝐻2 − (2.1)
𝑔𝐴 𝑑𝑡 2𝑔𝐷 𝐴2
donde
𝐻1 carga piezométrica en el tanque aguas arriba, en m
𝐻2 carga piezométrica en el tanque aguas abajo, en m
𝑄 gasto en la tubería, en m3/s
𝐴 sección transversal de la tubería, en m²
𝐷 diámetro de la tubería que separa al tanque 1 del tanque 2, en m
𝑙 longitud de la tubería que separa al tanque 1 del tanque 2, en m
𝑓 factor de fricción de Darcy-Weisbach de la tubería que separa al tanque 1
del tanque 2
𝑔 aceleración de la gravedad, en m/s²
𝑡 tiempo, en s
6
En la ecuación 2.1, cuando la conducción está formada por tuberías con diferentes
características los términos 𝐼 = 𝑙 ⁄(𝑔𝐴), conocido como término inercial y 𝐾 = 𝑓𝑙⁄(2𝑔𝐷𝐴2 )
llamado constante de pérdida de carga, deben calcularse como indican las ecuaciones 2.2
y 2.3, donde 𝑙𝑖 , 𝐷𝑖 , 𝐴𝑖 y 𝑓𝑖 son la longitud, diámetro, área y factor de fricción de la tubería 𝑖.
𝑛
1 𝑙𝑖
𝐼= ∑ (2.2)
𝑔 𝐴𝑖
𝑖=1
𝑛
1 𝑓𝑖 𝑙𝑖
𝐾= ∑ (2.3)
2𝑔 𝐷𝑖 𝐴𝑖2
𝑖=1
donde
𝐼 término inercial, en s2/m2
𝐾 constante de pérdida de carga, en s2/m5
𝑛 número de tuberías
Si es necesario, las pérdidas que introducen las válvulas, orificios, cambios de dirección,
ampliaciones, reducciones, etc., pueden incluirse en la constante de pérdida de carga 𝐾.
Dinámica de la superficie libre del agua en un tanque. Considerando que puede

despreciarse la inercia del agua en el tanque o torre, la dinámica de la superficie libre del
agua (SLA) en una torre está dada por la ecuación 2.4.
𝑑𝐻𝑖 𝑄𝑇𝑖
= (2.4)
𝑑𝑡 𝐴 𝑇𝑖
donde
𝐻𝑖 elevación de la superficie libre del agua en el tanque 𝑖, en m
𝑄𝑇𝑖 gasto que ingresa al tanque 𝑖 (positivo) o que sale (negativo), en m³/s
𝐴 𝑇𝑖 área transversal del tanque 𝑖, en m2
La solución de las ecuaciones 2.1 y 2.4 se obtiene a través de un esquema numérico, para
ello las derivadas temporales de la carga y el gasto se escriben en diferencias finitas
(ecuación 2.5), de tal forma que éstas toman la forma de las ecuaciones 2.6 y 2.7,
respectivamente.
𝑑𝑥 𝑥(𝑡 + ∆𝑡) − 𝑥(𝑡)

= (2.5)
𝑑𝑡 ∆𝑡
∆𝑡
𝑄(𝑡 + ∆𝑡) = 𝑄(𝑡) + [𝐻1 (𝑡 + ∆𝑡) − 𝐻2 (𝑡 + ∆𝑡) − 𝐾𝑄(𝑡)|𝑄(𝑡)|] (2.6)
𝐼
7
∆𝑡
𝐻𝑖 (𝑡 + ∆𝑡) = 𝐻𝑖 (𝑡) + 𝑄 (𝑡 + ∆𝑡) (2.7)
𝐴 𝑇𝑖 𝑇𝑖
donde
∆𝑡 intervalo de integración, en s
Ejemplo 1. En la figura 2.3 se muestra un esquema de la conexión en serie de plantas

de bombeo. El objetivo del ejemplo es determinar las variaciones de las elevaciones de los
niveles de agua en las torres TO y TS cuando se suspende simultáneamente el bombeo
en las plantas 1 y 2 (pb1 y pb2).
TO TS
1906 msnm 1912 msnm
1900.63 msnm
1895.11 msnm
1855 msnm
𝑄3
1862 msnm
Planta de
𝑄2 bombeo 2
𝑄1 𝑄2
Planta de bombeo 1
Figura 2.3 Cuando se suspenden los bombeos se presenta un transitorio de oscilación

de masa entre las torres TO y TS
La conducción que une a las torres tiene es de longitud 𝑙 = 3540 m, diámetro 𝐷 = 2.51 𝑚
y factor de fricción de Darcy-Weisbach 𝑓 = 0.012. Ambas torres tienen un diámetro de 10
m, sus desplantes y coronamientos se muestran en la figura 2.3. Durante la operación en
flujo estable 𝑄1 = 𝑄2 = 𝑄3 = 12.52 m3/s y las elevaciones de las SLA en la torre TO (𝐻1 =
1900.63 m) y en la torre TS (𝐻2 = 1895.11 m) permanecen constantes. Si se considera que
el paro del bombeo ocurre en 𝑡 = 0, entonces para 𝑡 > 0, 𝑄1 = 𝑄3 = 0.
Cuando las plantas de bombeo pb1 y pb2 suspenden la operación, las válvulas en las
descargas ambas bombas se cierran de tal forma que 𝑄1 = 𝑄3 = 0. Debido a la diferencia
de energías entre las torres TO y TS y siendo mayor en la TO que en la TS, habrá un flujo
8
desde la torre TO hacia la torre TS, el agua saldrá de la torre TO y entrará a la torre TS.
Dependiendo de las características de la tubería de la conducción que une a las torres y
las del flujo en estado permanente, el nivel del agua en la torre TS puede llegar a ser
superior que en la torre TO, en ese caso cambiará el sentido del flujo, el agua saldrá de la
torre TS y entrará a la torre TO. Mientras haya un desbalance de energía entre las torres
el proceso descrito volverá a repetirse, los niveles de agua en ambas subirán y bajaran
siguiendo un movimiento oscilatorio hasta que, debido a la fuerza de fricción entre la tubería
y el agua, la energía del sistema se disipe, a partir de entonces las elevaciones de los
niveles de agua en las torres serán iguales y no habrá flujo en la conducción.
A través de las ecuaciones 2.2 y 2.3 se obtienen el término inercial y la constante de pérdida
de carga de la conducción,
3540 𝑠2 0.012(3540) 𝑠2
𝐼= = 72.93 2 𝐾= = 0.035
9.81(4.95) 𝑚 2(9.81)(2.51)(4.95)2 𝑚5
Al aplicar la ecuación 2.7 para obtener los niveles del agua en las torres y la ecuación 2.6
para determinar el gasto en la conducción (𝑄2 , ver figura 2.3) se obtiene el sistema de
ecuaciones 2.8 a 2.10.
∆𝑡
𝐻TO (𝑡 + ∆𝑡) = 𝐻TO (𝑡) + (−𝑄2 (𝑡)) (2.8)
78.54
∆𝑡
𝐻TS (𝑡 + ∆𝑡) = 𝐻TS (𝑡) + 𝑄 (𝑡) (2.9)
78.54 2
∆𝑡
𝑄2 (𝑡 + ∆𝑡) = 𝑄2 (𝑡) + [𝐻 (𝑡 + ∆𝑡) − 𝐻TS (𝑡 + ∆𝑡) − 0.035𝑄2 (𝑡)|𝑄2 (𝑡)|] (2.10)
72.93 TO
El método iterativo para resolver el sistema de ecuaciones puede programarse en una hoja
de cálculo. En la Tabla 2.1 se muestran, para un incremento de integración de 0.1 s, los
resultados obtenidos para el primer segundo de cálculo, observa que para 𝑡 = 0 las
condiciones del flujo son las del estado estacionario.
9
Tabla 2.1 Niveles del agua en las torres y gasto en la conducción
𝑡 𝐻TO 𝐻TS 𝑄2
(s) (m) (m) (m3/s)
0 1900.63 1895.11 12.520
0.1 1900.61 1895.12 12.520
0.2 1900.60 1895.14 12.520
0.3 1900.58 1895.16 12.520
0.4 1900.57 1895.17 12.520
0.5 1900.55 1895.19 12.519
0.6 1900.53 1895.20 12.519
0.7 1900.52 1895.22 12.519
0.8 1900.50 1895.23 12.518
0.9 1900.49 1895.25 12.518
1 1900.47 1895.27 12.518
En las gráficas de la figura 2.4 se presentan las variaciones de los niveles del agua en las
torres, observa que cuando el nivel del agua en la torre TO disminuye en la torre TS
aumenta y viceversa. La simetría que se observa en la gráfica, cuando en una torre el nivel
es máximo en la otra es mínimo, se debe a que los diámetros de las torres son iguales.
En la Tabla 2.2 se indican los desplantes de las torres, las elevaciones de los bordos y los
niveles máximos y mínimos que se presentan durante el transitorio de oscilación de masa,
se concluye que para este transitorio las torres no se vacían y tampoco se desborda el
agua.
1906
1902
Elevación de la SLA (msnm)
1898
1894
TO TS
1890
0 1 2 3 4 5 6 7 8 9 10
Tiempo (min)
Figura 2.4 Variación de los niveles de agua en las torres de oscilación
10
Tabla 2.2 Elevaciones de las torres de oscilación
Nivel
Desplante Corona máximo Nivel mínimo
(msnm) (msnm) (msnm) (msnm)
TO 1855 1906 1902.91 1891.07
TS 1862 1912 1904.67 1892.82
Existen varias expresiones para determinar el factor de fricción para flujos

turbulentos, como lo son en los sistemas de bombeo, entre ellas la propuesta por
Swamee y Jain (Swamee, 1993) expresada en la ecuación 2.11.
−2
𝜀 5.74
𝑓 = 1.325 [𝑙𝑛 ( + )] para 𝑅𝑒 > 4000 (2.11)
3.7 𝐷 𝑅𝑒 0.9
donde
𝑅𝑒 = 𝑣𝐷⁄𝜗 número de Reynolds, adimensional
𝑣 velocidad del flujo, en m/s
𝜗 viscosidad del fluido, en m2/s
𝜀 promedio de la rugosidad de la tubería, en m
𝑅 radio de la tubería, en m
11
Ejemplo 2. El objetivo del ejemplo es determinar los niveles de agua máximo y mínimo
que se presentan en la torre de sumergencia cuando se pone en marcha la bomba de la
figura 2.5. Se considera que el área del tanque de entrega es lo suficientemente grande
que puede suponerse que el nivel del agua permanece en la elevación 1900 msnm durante
el transitorio.
Tanque de
entrega
1900 msnm
𝐷
Torre de Tramo 4
sumergencia
1775 msnm 𝐻𝑇
𝑄1 Tramo 2 𝑄3
Presa 𝑇
𝑐
𝑄2
𝑃
Bomba
Tramo 1
𝑎 𝑏
Tramo 3
Figura 2.5 Al iniciar el bombeo se presenta una oscilación de masa entre la presa y la
torre de sumergencia
En la Tabla 2.3 se presentan las características de las tuberías que forman la conducción,
las cargas piezométricas con gasto nulo se indican en la figura 2.5. En la figura 2.6 se
muestra la curva carga-gasto de la bomba, la cual puede representarse con el polinomio
de la ecuación 2.12, donde 𝐻𝐵 es la carga que entrega la bomba y 𝑄 el gasto bombeado.
Tabla 2.3 Característica de las tuberías
𝑙 𝐷
Tramo 𝑓(D-W)
(m) (m)
100.00 2.90 0.012
Sección 𝑃 - sección 𝑎 1 512.45 3.10 0.012
1345.00 3.55 0.012
Sección 𝑇 - sección 𝑎 2 220.52 3.10 0.012
Sección 𝑎 – sección 𝑏 3 95.00 3.10 0.012
Sección 𝑐 – sección 𝐷 4 4000.00 3.10 0.012
𝐻𝐵 = 169.94 − 6.7692 𝑄 + 0.8904 𝑄 2 − 0.4942 𝑄 3 (2.12)
12
180
170
160
150
Carga HB, en m
140
130
120
110
100
90
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Gasto Q, en m /s
3
Figura 2.6 Curva la bomba
Los términos inerciales y las constantes de pérdida de carga de cada tramo se obtienen a
través de las ecuaciones 2.2 y 2.3.
100 512.45 1345 𝑠2

𝐼1 = + + = 22.32
9.81(6.10) 9.81(7.55) 9.81(9.90) 𝑚2
0.012(100) 0.012(512.45) 0.012(1345) 𝑠2

𝐾1 = + + = 0.005 5
2(9.81)(2.90)(6.61)2 2(9.81)(3.10)(7.55)2 2(9.81)(3.55)(9.90)2 𝑚
220.52 𝑠2
𝐼2 = = 2.98
9.81(7.55) 𝑚2
0.012(220.52) 𝑠2
𝐾2 = = 0.001
2(9.81)(3.10)(7.55)2 𝑚5
95 𝑠2
𝐼3 = = 1.283
9.81(7.55) 𝑚2
0.012(95) 𝑠2
𝐾3 = = 0.0003
2(9.81)(3.10)(7.55)2 𝑚5
4000 𝑠2
𝐼4 = = 54.006
9.81(7.55) 𝑚2
13
0.012(4000) 𝑠2
𝐾4 = = 0.0138
2(9.81)(3.10)(7.55)2 𝑚5
El nivel del agua en la torre de sumergencia se obtiene resolviendo el sistema de

ecuaciones determinado por la ecuación de continuidad (ecuación 2.13), las ecuaciones
de la energía entre la presa y la sección 𝑎 (ecuación 2.14), entre la torre y la sección 𝑎
(ecuación 2.15), entre las secciones 𝑎 y 𝑏 (ecuación 2.16), entre la sección 𝑐 y la descarga
𝐷 (ecuación 2.17), la variación del nivel del agua en la torre (ecuación 2.18) y la relación
de carga entre las secciones 𝑏 y 𝑐 (ecuación 2.19), ver figura 2.5, y la curva de operación
de la bomba (ecuación 2.12).
𝑄1 + 𝑄2 = 𝑄3 (2.13)
𝑑𝑄1
𝐻𝑃 = 𝐻𝑎 + 𝐾1 𝑄1 |𝑄1 | + 𝐼1 (2.14)
𝑑𝑡
𝑑𝑄2
𝐻𝑇 = 𝐻𝑎 + 𝐾2 𝑄2 |𝑄2 | + 𝐼2 (2.15)
𝑑𝑡
𝑑𝑄3
𝐻𝑎 = 𝐻𝑏 + 𝐾3 𝑄3 |𝑄3 | + 𝐼3 (2.16)
𝑑𝑡
𝑑𝑄3
𝐻𝑐 = 𝐻𝐷 + 𝐾4 𝑄3 |𝑄3 | + 𝐼4 (2.17)
𝑑𝑡
𝑑𝐻𝑇 𝑄2
= − (2.18)
𝑑𝑡 𝐴𝑇
𝐻𝑏 = 𝐻𝑐 − 𝐻𝐵 (2.19)
Restando la ecuación 2.14 de la 2.15, sustituyendo la ecuación 2.13 y despejando 𝑑𝑄2 ⁄𝑑𝑡
se obtiene la ecuación 2.20.
𝑑𝑄2 1 𝑑𝑄3
= [𝐻𝑇 − 𝐻𝑃 + 𝐾1 𝑄1 |𝑄1 | − 𝐾2 𝑄2 |𝑄2 | + 𝐼1 ] (2.20)
𝑑𝑡 𝐼1 + 𝐼2 𝑑𝑡
Sustituyendo la ecuación 2.17 en la 2.19 y la que resulta en la ecuación 2.16, se obtiene la

ecuación 2.21.
𝑑𝑄3 1
= [𝐻 − 𝐻𝐷 + 𝐻𝐵 − (𝐾3 + 𝐾4 )𝑄3 |𝑄3 |] (2.21)
𝑑𝑡 𝐼3 + 𝐼4 𝑎
14
El sistema determinado por las ecuaciones 2.20 y 2.21 puede resolverse con el método de
Euler modificado, para ello la ecuación 2.21 se escribe como la ecuación 2.22 y la ecuación
2.20 como la ecuación 2.23.
∆𝑡
𝑄3 (𝑡 + ∆𝑡) = 𝑄3 (𝑡) + [𝐻 (𝑡) − 𝐻𝐷 + 𝐻𝐵 (𝑡) − (𝐾3 + 𝐾4 )𝑄3 (𝑡)|𝑄3 (𝑡)|] (2.22)
𝐼3 + 𝐼4 𝑎
∆𝑡
𝑄2 (𝑡 + ∆𝑡) = 𝑄2 (𝑡) + [𝐻 (𝑡) − 𝐻𝑃 + 𝐾1 𝑄1 (𝑡)|𝑄1 (𝑡)| − 𝐾2 𝑄2 (𝑡)|𝑄2 (𝑡)|] +
𝐼1 + 𝐼2 𝑇
𝐼1
+ (𝑄 (𝑡 + ∆𝑡) − 𝑄3 (𝑡)) (2.23)
𝐼1 + 𝐼2 3
El nivel del agua en la torre se obtiene a través de la ecuación 2.24, la carga que entrega
la bomba (𝐻𝐵 ) con la ecuación 2.12 y el gasto 𝑄1 a través de la ecuación de continuidad,
ecuación 2.14.
∆𝑡
𝐻𝑇 (𝑡 + ∆𝑡) = 𝐻𝑇 (𝑡) − 𝑄 (𝑡 + ∆𝑡) (2.24)
A𝑇 2
Antes de poner en marcha a la bomba los gastos 𝑄1 = 𝑄2 = 𝑄3 = 0, los niveles de agua

en la torre (𝐻𝑇 ) y la presa (𝐻𝑃 ) están en la elevación 1755 msnm y la carga aguas abajo de
la bomba (𝐻𝑐 ), determinada por la elevación del agua en el tanque de entrega (𝐻𝐷 ) es 1900
msnm. Tanto 𝐻𝑃 como 𝐻𝐷 permanecen constantes durante todo el transitorio.
El método iterativo para resolver el sistema de ecuaciones puede programarse en una hoja
de cálculo. En la Tabla 2.4 se muestran, para un incremento de integración de 0.1 s, los
resultados obtenidos para el primer segundo de cálculo, observa que para 𝑡 = 0 las
condiciones del flujo son las del estado estacionario.
En la figura 2.7 se muestran los gastos que entregan y reciben la presa (𝑄1 ) y la torre de
sumergencia (𝑄2 ) y el gasto que se establece en la conducción, que es igual al gasto
bombeado (𝑄3 ). Durante el primer minuto después del arranque de la bomba la aportación
más importante la entrega la torre, esto se debe a que la distancia entre la torre y la bomba
es menor que la distancia entre la presa y la bomba. Conforme transcurre el tiempo el gasto
en la torre tiende a cero y el gasto en la conducción es aportado por la presa.
15
Tabla 2.4 Gastos, nivel del agua en la torre y carga de bombeo para el primer segundo de
cálculo
𝑡 𝑄3 𝑄2 𝑄1 𝐻𝑇 𝐻𝐵
(s) (m3/s) (m3/s) (m3/s) (msnm) (m)
0.0 0.000 0.000 0.000 1775.00 169.40
0.1 0.080 0.071 0.009 1775.00 168.86
0.2 0.156 0.137 0.018 1775.00 168.36
0.3 0.231 0.203 0.027 1775.00 167.88
0.4 0.305 0.269 0.036 1775.00 167.41
0.5 0.378 0.333 0.045 1775.00 166.94
0.6 0.450 0.397 0.053 1775.00 166.49
0.7 0.522 0.460 0.062 1775.00 166.04
0.8 0.593 0.523 0.070 1775.00 165.60
0.9 0.663 0.584 0.078 1775.00 165.16
1.0 0.732 0.646 0.086 1775.00 164.73
6
Gasto Q, en m3/s
-2
Q1 Q2 Q3
-4
0 5 10 15
Tiempo t, en min
Figura 2.7 Gasto entregado por la presa (𝑄1 ), la torre (𝑄2 ) y bombeado (𝑄3 )
En la figura 2.8 se muestra la variación del nivel del agua en la torre y el gasto que entrega
o recibe. El nivel mínimo que alcanza el agua está en la elevación 1773.03 msnm, para
que la torre no se vacíe y se evite que ingrese aire a la conducción durante el arranque de
la bomba la torre debe estar desplantada por debajo de esa elevación. El nivel máximo que
se presenta está en la elevación 1776.49 msnm, su corona debe estar arriba de esta
elevación para que durante el arranque de la bomba no se desborde el agua. La amplitud
16
de la oscilación del tirante de agua disminuye conforme transcurre el tiempo debido a la

pérdida de energía por fricción.
Nivel del agua en la torre HT, en msnm 1777 4
1776 2
Gasto Q, en m3/s
1
1775 0
-1
1774 -2
-3
HT Q2
1773 -4
0 5 10 15
Tiempo t, en min
Figura 2.8 Nivel del agua en la torre de sumergencia y gasto (𝑄2 ) que entrega o recibe
2.3 Golpe de ariete
Los transitorios rápidos o de golpe de ariete se caracterizan porque los cambios en la

presión y velocidad del fluido ocurren de manera abrupta, en este estado de operación es
cuando se presentan los valores extremos de la presión capaces de colapsar o reventar
las tuberías. El arranque de equipos de bombeo y el paro del bombeo por corte de la
energía eléctrica, así como la apertura o cierre súbito de válvulas durante la operación
estacionaria en sistemas de bombeo, son ejemplos de eventos que generan este tipo de
transitorios.
Como consecuencia de un golpe de ariete, en el interior de algunas secciones de la

conducción de un sistema de bombeo pueden presentarse:
 Sobrepresiones, presión mayor que la de operación en régimen estacionario. Si ésta es

superior a la resistencia nominal de la tubería, la conducción está en riesgo de falla o
rotura por sobrepresión.
 Depresiones, presión menor que la presión atmosférica. Las depresiones favorecen el

ingreso de aire a la conducción, pueden causar tanto que las tuberías colapsen como
que se presente el fenómeno de separación de columna.
17
Las depresiones favorecen que el aire mezclado en el agua se separe del agua, se
aglomere y forme cavidades de aire, las cuales pueden contener vapor de agua en caso
de que se hubiera alcanzado la presión de vapor del agua, la columna de agua pierde
su continuidad presentándose la llamada separación de columna.
Cuando el origen del transitorio es un paro súbito del bombeo, la dinámica de la bomba
queda bajo la acción del flujo transitorio y la inercia de las masas en rotación del grupo
motor-bomba. En acueductos con pendiente ascendente (la elevación del tanque de
entrega es mayor que la de la planta de bombeo) a los pocos segundos de haberse
suspendido el bombeo el sentido del flujo se invierte. Mientras no cierren las válvulas que
lo eviten, el agua entrará a las bombas haciéndolas girar en sentido inverso pudiendo
alcanzar velocidades mayores que las de diseño. Las partes que giran serán sometidas a
esfuerzos excesivos debidos a la fuerza centrífuga, la integridad estructural de estos
elementos se verá seriamente comprometida.
Por lo señalado anteriormente es evidente la importancia que tiene conocer los valores
máximos y mínimos de la presión en el interior de las tuberías y la velocidad de la bomba
como consecuencia de un golpe de ariete. Las presiones se determinan a través del modelo
de golpe de ariete y la velocidad de la bomba a través del modelo que representa su
dinámica.
2.3.1 Modelo matemático del golpe de ariete. En el modelo de golpe de ariete

se considera que el fluido es compresible y el material con el que están fabricadas las
tuberías es elástico.
Los cambios en la presión del fluido que se presentan después de perturbar rápidamente
la velocidad del flujo se transmiten como ondas de presión con una velocidad finita a la que
se le llama celeridad, la cual se determina como se indica a continuación.
La celeridad depende principalmente de las características elásticas del fluido y del material
con el que está fabricada la tubería, así como, la relación espesor-diámetro y la forma como
está sujeta, debido a que de ello depende la libertad para expandirse longitudinal y
axialmente. Para la forma más común de la instalación de la tubería, totalmente anclada y
con juntas de expansión, la celeridad se determina con la ecuación 2.25, para otros tipos
de sujeción se puede consultar a Cafaggi et al (2016), Chaudhry (2014) y Wylie y Streeter
(1978).
𝐸𝑉
√
𝜌
𝑎= (2.25)
𝐷
√1 + 𝐸𝑉
𝑒𝐸𝑡
18
donde
𝑎 celeridad, en m/s
𝐷 diámetro de la tubería, en m
𝐸𝑡 módulo de elasticidad del material de la tubería, en N/m2
𝑒 espesor de la tubería, en m
𝐸𝑉 módulo de elasticidad volumétrica del agua, en N/m2
𝜌 densidad del agua, en kg/m3
Mientras más rígidos son los materiales el módulo de elasticidad es mayor, entonces, si se
considera que el módulo de elasticidad de tuberías totalmente rígidas es infinito, el
denominador de la ecuación 2.25 es igual a uno, por lo que para tuberías rígidas la
celeridad está dada por la ecuación 2.26.
𝐸𝑉
𝑎=√ (2.26)
𝜌
La celeridad de la onda de presión en una tubería rígida que conduce agua a 20°C (𝐸𝑉 =
2.2 (109 ) N/m2 y 𝜌 = 998.2 kg/m3) es 1484.5 m/s, cuando el agua está a 50°C (𝐸𝑉 =
2.3 (109 ) N/m2 y 𝜌 = 988.1 kg/m3) es 1525.68 m/s, en el sistema TransDocencia2-0 la
celeridad máxima que se permite es 1530 m/s.
El hecho de que el cambio en la presión del fluido se transmita a través de una onda que
se propaga dentro de la conducción significa que la presión en el interior de la tubería va
cambiando conforme avanza esa onda de presión. Cuando en un sistema de bombeo como
el de la figura 2.9 deja de operar súbitamente la bomba, inmediatamente aguas abajo de
la bomba la presión disminuye, en ese mismo instante la presión en la mitad de la
conducción o en cualquier otro sitio, no ha ocurrido ningún cambio.
Con la finalidad de imaginarnos lo que sucede en el interior de la conducción cuando se

suspende súbitamente el bombeo en un sistema, consideremos el de la figura 2.9 en el
que se considera el área del tanque de descarga lo suficientemente grande de tal forma
que el nivel del agua permanece constante durante un transitorio. La válvula en la descarga
de la bomba cierra inmediatamente después del paro de la bomba, la perturbación de la
presión se transmite con celeridad 𝑎 y se desprecian las pérdidas de energía debidas a la
fricción. En flujo permanente la carga es 𝐻𝑜 y el gasto 𝑄𝑜 .
19
Tanque de descarga
𝐻0
𝐻𝑑𝑒𝑠
𝑄0
𝑥
Figura 2.9 Después de un paro súbito del bombeo se genera una onda de baja presión
que se transmite a través del agua hacia toda la conducción
En las figuras 2.10 a 2.17 se muestran en el plano posición-tiempo los cambios de presión
en el interior de la tubería conforme avanza la onda de presión, en la parte inferior se indica
la presión en el interior de la conducción.
En 𝑡 = 0 se suspende súbitamente el bombeo y se cierra la válvula de descarga.

Inmediatamente aguas abajo de la bomba (𝑥 = 0) se genera una onda de baja presión que
a través del agua se propaga con celeridad 𝑎 a toda la conducción, figura 2.10.
En 𝑡 = 𝑙 ⁄𝑎, la onda de baja presión habrá recorrido toda la tubería de longitud 𝑙, figura
2.11. Al incidir en el tanque con carga constante se impone la carga del tanque (se dice
que hay un nodo de presión), la onda de presión se refleja con una amplitud tal que la
presión en el interior de la conducción vuelve a la que había antes de presentarse el
transitorio, como se muestra en la figura 2.12.
20
4𝑙
0.8
𝑎
3𝑙
0.6
Tiempo 𝑡, en s 𝑎
2𝑙
0.4
𝑎
𝑙
0.2
𝑎
0 H0+DH
0 𝑙1 𝑙
2 3𝑙
3 𝑙
4
4 2 4
Posición 𝑥, en m
𝐷− 𝐷 H0-DH
Carga = H0 - DH Carga = H0 𝐷
Figura 2.10 En 𝑡 = 0 se suspende súbitamente el bombeo, se genera una onda de baja

presión que avanza en sentido aguas arriba
4𝑙
0.8
𝑎
3𝑙
0.6
𝑎
Tiempo 𝑡, en s
2𝑙
0.4
𝑎
𝑙
0.2
𝑎
0
0 𝑙1 2
𝑙 3
3𝑙 4
𝑙
4 2 4
𝐷− 𝐷 Carga = H0 - DH
Figura 2.11 En 𝑡 = 𝑙 ⁄𝑎, la onda de baja presión ha recorrido toda la tubería de longitud 𝑙
21
4𝑙
0.8
𝑎
3𝑙
0.6
𝑎
Tiempo 𝑡 , en s
2𝑙
0.4
𝑎
𝑙
0.2
𝑎
0
0 𝑙1 𝑙
2 3𝑙
3 𝑙
4
4 2 4
𝐷− 𝐷 Carga = H0 - DH Carga = H0 𝐷
Figura 2.12 Entre 𝑡 = 𝑙 ⁄𝑎 y 𝑡 = 2𝑙 ⁄𝑎, se restablece la presión que había en flujo

permanente, la onda de presión viaja en sentido aguas abajo
En 𝑡 = 2𝑙 ⁄𝑎, la presión en el interior de la tubería es igual a la que había en flujo

permanente, figura 2.13.
4𝑙
0.8
𝑎
3𝑙
0.6
𝑎
Tiempo 𝑡 , en s
2𝑙
0.4
𝑎
𝑙
0.2
𝑎
0
0 𝑙1 𝑙
2 3𝑙
3 𝑙
4
4 2 4
𝐷 Carga = H0
Figura 2.13 En 𝑡 = 2𝑙 ⁄𝑎 la presión en toda la conducción es igual a la de flujo

permanente
22
En 𝑡 = 2𝑙 ⁄𝑎 la onda de presión llega a la válvula de descarga completamente cerrada,

inmediatamente aguas abajo de la válvula la presión aumenta. Se genera una onda de alta
presión que se transmite hacia el tanque de descarga aumentando la presión como se
muestra en la figura 2.14.
4𝑙
0.8
𝑎
3𝑙
0.6
𝑎
Tiempo 𝑡 , en s
2𝑙
0.4
𝑎
𝑙
0.2
𝑎
0
0 𝑙1 𝑙
2 3𝑙
3 𝑙
4
4 2 4
𝐷+ 𝐷 Carga = H0 + DH Carga = H0 𝐷
Figura 2.14 Entre 𝑡 = 2𝑙 ⁄𝑎 y 𝑡 = 3𝑙 ⁄𝑎, una onda de alta presión se transmite en sentido
aguas abajo
En 𝑡 = 3𝑙 ⁄𝑎 la onda de alta presión habrá recorrido toda la conducción, figura 2.15.
4𝑙
0.8
𝑎
3𝑙
0.6
𝑎
Tiempo 𝑡 , en s
2𝑙
0.4
𝑎
𝑙
0.2
𝑎
0
0 𝑙1 𝑙
2 3𝑙
3 𝑙
4
4 2 4
𝐷+ 𝐷 Carga = H0 + DH
Figura 2.15 En 𝑡 = 3𝑙 ⁄𝑎 la onda de alta presión termina de recorrer la tubería
23
Cuando la onda de alta presión incide en el tanque de descarga con carga constante vuelve
a imponerse la carga del tanque, conforme avanza la onda reflejada en el tanque se
restablece la carga del flujo que había en el estado permanente, figura 2.16.
4𝑙
0.8
𝑎
3𝑙
0.6
𝑎
Tiempo 𝑡 , en s
2𝑙
0.4
𝑎
𝑙
0.2
𝑎
0
0 𝑙1 𝑙
2 3𝑙
3 𝑙
4
4 2 4
𝐷+ 𝐷 Carga = H0 + DH Carga = H0 𝐷
Figura 2.16 Entre 𝑡 = 3𝑙 ⁄𝑎 y 𝑡 = 4𝑙 ⁄𝑎, se restablece la presión del flujo permanente
En 𝑡 = 4𝑙 ⁄𝑎, la onda reflejada en el tanque de descarga llega a 𝑥 = 0, se concluye el

primer periodo o ciclo de la onda de presión, el cual se determina con la ecuación 2.27, y
comienza otro igual, figura 2.17.
4𝑙
𝑇= (2.27)
𝑎
donde
𝑇 periodo de la onda de presión, en s
𝑙 longitud de la tubería, en m
24
4𝑙
0.8
𝑎
3𝑙
0.6
𝑎
Tiempo 𝑡 , en s 2𝑙
0.4
𝑎
𝑙
0.2
𝑎
0
0 𝑙1 𝑙
2 3𝑙
3 𝑙
4
4 2 4
𝐷 Carga = H0
Figura 2.17 En 𝑡 = 4𝑙 ⁄𝑎 se termina el primer ciclo de la onda de presión
En los comentarios anteriores únicamente se ha hecho referencia a los cambios de la

presión, sin embargo, al haber cambios en la presión también los hay en la velocidad del
agua y por tanto en el gasto, como se muestra en las figuras 2.18a y 2.18b
.
Figura
G4 2.18a Cambio de la presión en función del tiempo y la posición
25
G3 Figura 2.18b Cambio del gasto en función del tiempo y la posición
Las relaciones entre la carga y el gasto en función de la posición y el tiempo se obtienen a

través del modelo de golpe de ariete o modelo de tubería elástica y fluido compresible, las
ecuaciones que lo determinan son la ecuación dinámica (2.28) y la de continuidad (2.29),
el desarrollo para obtenerlas puede consultarse en Abreu et al., (1995) o Guarga et al
(1985) algunos autores le llaman modelo simplificado de golpe de ariete debido a las
simplificaciones que se realizan para obtenerlo.
En el sistema de ecuaciones 2.28 y 2.29 hay dos variables independientes, la posición (𝑥)
y el tiempo (𝑡), y dos variables dependientes, la carga (𝐻) y el gasto (𝑄).
𝜕𝐻 1 𝜕𝑄 𝑓 𝑄|𝑄|
+ + =0 (2.28)
𝜕𝑥 𝑔𝐴 𝜕𝑡 2g𝐷 𝐴2
𝜕𝑄 𝑔𝐴 𝜕𝐻
+ 2 =0 (2.29)
𝜕𝑥 𝑎 𝜕𝑡
El área transversal de la conducción (𝐴) y el diámetro (𝐷) son características de la

conducción, la celeridad (𝑎) además depende de las características del fluido (ecuación
2.25), el factor de fricción (𝑓) varía con el número de Reynolds sin embargo, se puede
considerar constante porque los efectos de su variación en el estado transitorio son
despreciables. Todos los parámetros mencionados son invariantes en el tiempo, pero
pueden cambiar a lo largo de la conducción porque, en general, las conducciones de los
sistemas de bombeo están formadas por tuberías o tramos de distintas características.
26
Este es un sistema de ecuaciones hiperbólicas de primer orden sin solución analítica, la

solución se obtiene a través de métodos numéricos, en el sistema TransDocencia2-0 se
utiliza el método de las características.
Método de las características. A través del método de las características, el sistema de

ecuaciones 2.28 y 2.29 se transforma en el sistema de ecuaciones desacoplado con
derivadas totales en el tiempo determinado por las ecuaciones 2.30 y 2.31. En Abreu et al,
(1995) y Guarga et al, (1985) se presenta el desarrollo para pasar de un sistema al otro.
𝑑𝑄 𝑔𝐴 𝑑𝐻 𝑓 𝑑𝑥
+ + 𝑄|𝑄| = 0 cuando =𝑎 (2.30)
𝑑𝑡 𝑎 𝑑𝑡 2𝐷𝐴 𝑑𝑡
𝑑𝑄 𝑔𝐴 𝑑𝐻 𝑓 𝑑𝑥
− + 𝑄|𝑄| = 0 cuando = −𝑎 (2.31)
𝑑𝑡 𝑎 𝑑𝑡 2𝐷𝐴 𝑑𝑡
Las ecuaciones 2.28 y 2.29 son válidas en todo el espacio y el tiempo, las ecuaciones 2.30
y 2.31 únicamente son válidas en las regiones del espacio y el tiempo que satisfacen
𝑑𝑥 ⁄𝑑𝑡 = ±𝑎. Estas regiones son el conjunto de rectas en el plano 𝑥, 𝑡 señaladas en la
figura 2.18 que dividen al plano en las zonas perturbadas por la onda de presión de las
zonas que aún no han sido perturbadas, a estas rectas se les llama rectas características,
la recta con pendiente positiva es la característica positiva y la de pendiente negativa la
característica negativa, ver figuras 2.10 a 2.18. En las figuras 2.10 a 2.17, las
características positivas son las rectas de color rojo y las características negativas son las
de color azul.
Al definir las constantes 2.32 y 2.33 y escribir las ecuaciones 2.30 y 2.31 en diferencias
finitas (ecuación 2.5 y recordando que la carga y el gasto dependen de la posición y el
tiempo) se obtienen las ecuaciones 2.34 y 2.35.
𝑔𝐴
𝑐𝑎 = (2.32)
𝑎
𝑓∆𝑡
𝑐𝑓 = (2.33)
2𝐷𝐴
𝑄(𝑥, 𝑡 + ∆𝑡) − 𝑄(𝑥 − ∆𝑥, 𝑡) 𝐻(𝑥, 𝑡 + ∆𝑡) − 𝐻(𝑥 − ∆𝑥, 𝑡)

+ 𝑐𝑎 + 𝑐𝑓 𝑄(𝑥 − ∆𝑥, 𝑡)|𝑄(𝑥 − ∆𝑥, 𝑡)| = 0
∆𝑡 ∆𝑡
(2.34)
27
𝑄(𝑥, 𝑡 + ∆𝑡) − 𝑄(𝑥 + ∆𝑥, 𝑡) 𝐻(𝑥, 𝑡 + ∆𝑡) − 𝐻(𝑥 + ∆𝑥, 𝑡)

− 𝑐𝑎 + 𝑐𝑓 𝑄(𝑥 + ∆𝑥, 𝑡)|𝑄(𝑥 + ∆𝑥, 𝑡)| = 0
∆𝑡 ∆𝑡
(2.35)
Tomando en consideración la notación definida a través de las ecuaciones 2.36 a 2.38

donde 𝐹 representa a las variables 𝑄 o 𝐻, las definiciones de los parámetros 𝑐𝑝 y 𝑐𝑛
(ecuaciones 2.39 y 2.40) y sustituyendo en las ecuaciones 2.34 y 2.35, se obtienen la
ecuación de la característica positiva (ecuación 2.41) y de la característica negativa
(ecuación 2.42).
𝐹𝑃 = 𝐹𝑥𝑡+∆𝑡 = 𝐹(𝑥 , 𝑡 + ∆𝑡) (2.36)
𝐹𝐴 = 𝐹(𝑥 − ∆𝑥 , 𝑡) (2.37)
𝐹𝐵 = 𝐹(𝑥 + ∆𝑥 , 𝑡) (2.38)
𝑐𝑝 = 𝑄𝐴 + 𝑐𝑎 𝐻𝐴 − 𝑐𝑓 𝑄𝐴 |𝑄𝐴 | (2.39)
𝑐𝑛 = 𝑄𝐵 − 𝑐𝑎 𝐻𝐴 − 𝑐𝑓 𝑄𝐵 |𝑄𝐵 | (2.40)
𝑄𝑃 = 𝑐𝑝 − 𝑐𝑎 𝐻𝑃 (2.41)
𝑄𝑃 = 𝑐𝑛 + 𝑐𝑎 𝐻𝑃 (2.42)
donde
𝐹𝑃 valor de la variable 𝐹 en la posición 𝑥 en el tiempo 𝑡 + ∆𝑡
𝐹𝐴 valor de la variable 𝐹 en la posición 𝑥 − Δ𝑥 en el tiempo 𝑡
𝐹𝐵 valor de la variable 𝐹 en la posición 𝑥 + Δ𝑥 en el tiempo 𝑡
𝑐𝑝 parámetro de la característica positiva, en m3/s
𝑐𝑛 parámetro de la característica negativa, en m3/s
El análisis de los transitorios se realiza a través de un proceso iterativo en el tiempo con un

incremento de integración ∆𝑡 y discretizando a la conducción en secciones de cálculo
distantes ∆𝑥 entre sí, la primera sección está en 𝑥 = 0 y la última en 𝑥 = 𝑙 donde 𝑙 es la
longitud de la conducción, ∆𝑥⁄∆𝑡 = ±𝑎 (ecuaciones 2.30 y 2.31). La malla de cálculo del
método de las características se muestra en la figura 2.19.
Las condiciones iniciales del proceso iterativo son las condiciones de operación en flujo
permanente. A través de la figura 2.19 se observa que en la frontera aguas arriba (𝑥 =0)
sólo existe la característica negativa, ecuación 2.42, una ecuación con dos incógnitas la
28
otra ecuación que se necesita para completar el sistema se obtiene de la condición

hidráulica que impone la frontera. Lo mismo sucede en la frontera aguas abajo (𝑥 =𝑙) en la
que sólo existe la característica positiva (ecuación 2.41). En los capítulos 4 y 5 se presentan
algunas de las condiciones de frontera que existen en los sistemas de bombeo.
Característica negativa Característica positiva
𝑡
Tiempo
𝑡 + ∆𝑡 𝑃
𝑡 𝐴 𝐵
𝑥 − ∆𝑥 𝑥 𝑥 + ∆𝑥
Posición 𝑥
Condiciones de operación en flujo permanente

Secciones de cálculo o puntos interiores
Frontera aguas arriba
Frontera aguas abajo
Figura 2.19 Malla de cálculo del método de las características
La presión en el interior de la tubería que se presenta durante los transitorios por golpe de
ariete puede alcanzar valores lo suficientemente bajos como para generar la separación
de la columna líquida (figura 2.20) o colapsar tuberías (figura 2.21) o lo suficientemente
grandes como para reventarlas, figura 2.22.
29
Figura 2.20 Fenómeno de separación de la columna líquida causado por un golpe de

ariete (Grupo de Hidromecánica, IIUNAM)
Figura 2.21 Colapso de una tubería causado por la depresión originada por un golpe de
ariete (Grupo de Hidromecánica, IIUNAM)
30
Figura 2.22 Falla en una tubería causada por un golpe de ariete (Grupo de Hidromecánica,
IIUNAM)
31
SISTEMA TRANSDOCENCIA2-0
3
Sistema TransDocencia2-0
3.1 Introducción
El sistema TransDocencia2-0 permite simular la operación en flujo transitorio de sistemas
hidráulicos con tuberías a presión, en particular sistemas de bombeo o acueductos, a los
que de aquí en adelante se hará referencia. Para ello se cuenta con una interfaz tipo CAD
a través de la cual, mediante operaciones de arrastre de componentes gráficos se describe
el sistema que se desea analizar.
Una vez descritas las estructuras físicas que forman al sistema, la condición de operación
en flujo permanente y la maniobra que produce el flujo transitorio se pueden realizar los
cálculos numéricos que determinan los valores de las variables hidráulicas en flujo
transitorio. Con la finalidad de que el usuario utilice cómodamente el sistema
TransDocencia2-0, en este capítulo se presenta el uso general del mismo.
3.2 Descripción del sistema TransDocencia2-0

El sistema TransDocencia2-0 es un sistema tipo CAD, que ofrece la posibilidad de crear
gráficamente el arreglo físico de acueductos seleccionando los elementos hidráulicos que
los forman a través de íconos que los representan, los cuales se “conectan” a la tubería de
conducción.
El análisis del flujo transitorio de un sistema de bombeo se realiza a través de lo que se ha

llamado proyecto. Un proyecto, en el sistema TransDocencia2-0 es el estudio que se hace
32
para resolver el problema de sobrepresiones o depresiones, generado por golpe de ariete,

que se presenta en un acueducto particular, para ello se analizan alternativas de solución
a las que se les ha llamado propuestas, cada propuesta se analiza por separado.
En la primera propuesta de un proyecto, debe definirse el sistema de bombeo a través de

sus características físicas y de la condición de operación en flujo permanente. Una vez que
se ha creado la primera propuesta se realiza la simulación de la ejecución de una maniobra
que genera un flujo transitorio, con base en el análisis de los resultados obtenidos se crean
nuevas propuestas para proponer, en caso de ser necesarios, dispositivos de control de
transitorios o realizar simulaciones para determinar las dimensiones óptimas de los
dispositivos.
El sistema TransDocencia2-0 está constituido por dos sistemas de cómputo (dos archivos
ejecutables) que interactúan entre sí. Al ejecutar el programa TransDocencia2-0.exe el
primer sistema que se ejecuta es la interface de comunicación con el usuario, a través de
la cual se generan las propuestas. Para simular la ejecución de la maniobra que genera el
transitorio y realizar los cálculos numéricos que modelan al flujo y la operación de los
elementos que forman parte del acueducto, se ejecuta automáticamente el sistema
TransMoMaDoc2-0.exe (figura 3.1).
33
Interfaz del usuario
Módulo de captura Acueducto.trs
Componente
TransMoMaDoc2-0
Propuesta.txt
Módulo de simulación
Propuesta.txt
Módulo de
cálculos
numéricos
Módulo presentación
de resultados Archivos de
resultados *
* Los archivos de resultados poseen sufijos: _H, _HEX, _IMP, _Q, _RPM
Figura 3.1 Elementos del sistema TransDocencia2-0
Como puede observarse en la figura 3.1, la interfaz de usuario está constituida por tres
módulos.
 Captura. Permite al usuario utilizar un ambiente gráfico para dibujar el sistema e

introducir la información de cada componente, como resultado de ello por cada proyecto
se genera un archivo con extensión “.trs”, que constituye la base sobre la cual se puede
realizar el proceso de simulación y presentación de resultados. Este archivo contiene la
información de todas las propuestas del proyecto, no debe borrarse, de lo contrario no
podrán recuperarse los esquemas del acueducto ni generarse los archivos para la
simular la operación transitoria del sistema.
34
 Simulación. Una vez que se ha creado una propuesta y se quiere realizar una
simulación, a partir del archivo “.trs” se genera un archivo tipo texto (*.txt) para esa
propuesta, éste contiene la información y el formato adecuado para ser procesada por
el módulo de cálculos numéricos (Carmona y Carmona, 2014) el cual genera los
archivos de resultados de la simulación o reporta errores de ejecución.
 Resultados. Si la simulación tuvo éxito se generan los archivos que contienen los
resultados de los cálculos numéricos, en el capítulo 6 se presenta cómo se analiza la
información que contienen los archivos.
En las secciones siguientes se describe el uso del sistema TransDocencia2-0 desde el

punto de vista de interfaz de usuario, capacitando al lector en el uso cómodo del sistema.
3.3 Pantalla principal del sistema TransDocencia2-0

La pantalla principal del sistema se muestra en la figura 3.2, en ésta se distinguen las seis
áreas principales de trabajo.
1. Menús. En esta área se especifican las acciones que pueden realizarse con los
proyectos, con el sistema en general o elementos administrativos y de información.
2. Elementos de una propuesta. En esta área se localizan todos los elementos que
pueden utilizarse durante el proceso de esquematización o dibujo del acueducto, el
usuario toma un elemento y lo coloca en el área de dibujo, uniendo los diferentes
elementos por líneas que representan las tuberías que los conectan.
3. Administrador del proyecto. Aquí se puede visualizar las diferentes propuestas que
constituyen el proyecto, en su versión más simple un proyecto posee una sola
propuesta.
4. Área de dibujo y resultados. Aquí se realiza el proceso de esquematización, además

cuando se simula una propuesta los resultados de dicha simulación también se
presentarán en esta área.
5. Inspector de propiedades. Cada elemento que se emplea para esquematizar una

propuesta posee atributos que lo describen, como su color para mostrarlos
gráficamente, sus dimensiones físicas y según sea el caso, las condiciones de
operación en flujo permanente o la descripción de la maniobra que se realiza con ese
componente. Se emplea esta área para alimentar la propuesta con los valores
adecuados dependiendo del elemento seleccionado.
35
6. Mensajes. Esta área la emplea normalmente el sistema para notificar comentarios,

advertencias, errores y en general cualquier información que pueda ser útil al usuario
respecto del sistema.
3 4
5
1. Menús
2. Elementos de una propuesta
3. Administrador del proyecto
4. Área de dibujo y resultados
5. Inspector de propiedades
6. Mensajes
6
Figura 3.2 Pantalla principal del sistema TransDocencia2-0
Elementos del menú principal. A continuación, se describen las opciones del menú
principal.
Proyecto. En este apartado se localizan las opciones para administrar un proyecto y sus
propuestas. A través de las opciones que hay se pueden crear nuevos proyectos, abrir
proyectos existentes y guardar proyectos, además de imprimir el esquema con el que se
está trabajando.
36
Figura 3.3 Opción principal: Proyecto
Editar. Aquí están las opciones típicas para realizar operaciones de cortar, copiar, pegar y
eliminar selecciones dentro del área de dibujo, además de los comandos de deshacer y
rehacer.
Figura 3.4 Opción principal: Editar
Catálogos. Actualmente hay dos catálogos a saber: catálogo de curvas de bombas y

catálogo de curvas de válvulas. En esta versión del sistema TransDocencia2-0 los
37
catálogos no pueden editarse, únicamente pueden simularse las bombas y válvulas que
están en ellos.
Figura 3.5 Opción principal: Catálogos
El catálogo de curvas de bombas contiene curvas de operación de bombas en flujo

transitorio para las velocidades específicas 𝑁𝑆 = 29, 35, 38, 44, 46, 147 y 261 (rpm
m1.5 ⁄s0.5 m-0.75 ) equivalentes a 𝑛𝑆 = 1498, 1800, 1950, 2270, 2400, 7600 y 13500 (rpm
gpm0.5 pies-0.75 ). Al elegir esta opción se abre una ventana secundaria mostrando la
información correspondiente para cada tipo de bomba disponible en el sistema, como se
muestra en la figura 3.6
Figura 3.6 Ventana secundaria que se abre en la opción Catálogo de bombas.
38
Como se mencionó, también se cuenta con un catálogo de curvas de válvulas en el cual

se indica la curva de pérdida de carga en función de la apertura de la válvula. Al igual que
el catálogo de bombas, cuando se elige esta opción se muestra la información
correspondiente en una ventana secundaria, figura 3.7.
Figura 3.7 Ventana secundaria que se abre en la opción Catálogo de válvulas.
Simulación. Se ejecuta el módulo en el que se realizan los cálculos numéricos para simular
la operación en flujo transitorio del sistema de bombeo que tiene las características de la
propuesta que en ese momento se encuentra activa y visible en el “área de dibujo y
resultados”.
39
Figura 3.8 Opción principal: Simulación
Una vez que se han terminado los cálculos de la simulación, los resultados se muestran en
la misma área, “área de dibujo y resultados”. Antes de realizar una simulación, se verifica
que la información proporcionada en la propuesta sea correcta reportando en el área de
mensajes los problemas detectados.
Ventana. Esta opción del menú permite habilitar la visualización del inspector del
administrador de proyectos y de los paneles o barras de herramientas, con la finalidad de
que el usuario configure a su gusto el sistema.
Figura 3.9 Opción principal: Ventana
Acerca de. Muestra la información correspondiente a los créditos y propiedad intelectual

del sistema
40
3.4 Creación de un proyecto

Al iniciar el programa TransDocencia2-0 lo primero que se genera es un proyecto
conteniendo una propuesta, ambos si información. El nombre por omisión del proyecto es
“Sin titulo” y el de la propuesta es “Propuesta 1”.
Figura 3.10 Proyecto creado por omisión cuando se abre el TransDocencia2-0
El nombre del proyecto se modifica al momento de salvarlo. Por otro lado, el de la propuesta
se puede modificar en cualquier momento, para ello se selecciona la propuesta en el
administrador del proyecto y al dar un nuevo clic en el nombre de la propuesta
inmediatamente se podrá modificar el nombre (figura 3,11), una alternativa es el menú que
se despliega mediante el botón derecho colocando el menú sobre la propuesta (figura
3.12).
Figura 3.11 Cambio del nombre de una propuesta
41
Figura 3.12 Opciones de menú para una propuesta
Por medio del administrador de proyectos se pueden agregar nuevas propuestas, para ello
se coloca el ratón sobre el nombre del proyecto y se presiona el botón derecho del ratón,
se desplegará un menú con las opciones disponibles, se selecciona nueva propuesta.
Figura 3.13 Agregando una propuesta al proyecto
De esta manera se crea una nueva propuesta. Como puede observarse en la figura 3.14
se agrega una nueva pestaña en el área de dibujo.
Figura 3.14 Elementos para navegar entre propuestas
42
Para desplazarse entre propuestas de manera indistinta puede darse clic en la propuesta
deseada ya sea en el administrador del proyecto o en las pestañas dentro de la zona de
dibujo.
Una vez desarrollada la o las propuestas, se puede salvar el proyecto (y con él todas las
propuestas) para ello se emplean las opciones correspondientes en el menú principal
Proyecto.
3.5 Creación de una propuesta

Una propuesta es una representación esquemática de un acueducto. Esta representación
está formada por una secuencia de elementos hidráulicos unidos por líneas que
representan a las tuberías, la información de cada elemento se proporciona al sistema a
través del Inspector.
En la figura 3.15 se muestra el esquema de un acueducto, todos los elementos están

interconectados. En una propuesta solo puede haber un diagrama, es decir no se permite
más de un acueducto, en esta versión no se admiten derivaciones ni más de una
aportación.
Figura 3.15 Esquema de un acueducto desarrollado en TransDocencia2-0
43
Los elementos que pueden formar una propuesta se encuentran agrupados en las cuatro
pestañas que están dentro de la zona de elementos de una propuesta (figura 3.16), estas
son:
 Fronteras y maniobras aguas arriba
 Tubos y elementos entre tramos
 Maniobras entre tramos
 Fronteras y maniobras aguas abajo
Figura 3.16 Pestañas que agrupan los elementos que pueden emplearse para
esquematizar un acueducto
En cada pestaña se muestran los diferentes elementos disponibles. Para elegir un

elemento se da un clic con el ratón sobre el ícono correspondiente, a continuación, se sitúa
el cursor del ratón sobre el área de dibujo en la posición deseada y se da otro clic. El
sistema colocará el elemento seleccionado en la posición del cursor. Las propiedades del
elemento colocado se presentan en el Inspector de propiedades de manera automática,
figura 3.17.
Figura 3.17 Agregar un elemento al esquema
44
Todos los elementos del sistema, excepto las tuberías, poseen “puntos de anclaje”,
representados por un pequeño cuadro color amarillo dentro del dibujo (figura 3.18). Estos
elementos sirven para unir por medio de una tubería un elemento con otro. Si se coloca el
cursor sobre uno de estos puntos ancla el cursor cambia de una flecha a una mano, lo que
indica que ahí se puede conectar una tubería, al dar clic en ese punto y sin dejar de
presionar el ratón se puede trazar una línea hasta otro punto ancla de otro elemento o
simplemente dejar la tubería sin conectar, con la idea de conectarlo más tarde. Todas las
tuberías del esquema deben estar necesariamente conectadas en sus dos extremos a
algún elemento dentro del diagrama.
Figura 3.18 Puntos de anclaje (indicados por las flechas) y dos estructuras conectadas
por una tubería
3.6 Inspector de Propiedades

Cada elemento que se agrega al esquema de la propuesta requiere que el usuario
proporcione los datos que lo caracterizan, físicas o de operación. La información requerida
para cada elemento se muestra en el Inspector. Por ejemplo, tratándose de un paro súbito
de bombeo en la frontera aguas arriba del acueducto, pestaña Fronteras y maniobras
aguas arriba, la información requerida incluye el tipo de válvula en la descarga de la bomba,
características de la bomba, etc. (figura 3.19).
Como una propuesta tiene más de un elemento, cada vez que se selecciona un
componente en el área de dibujo, el Inspector cambia para mostrar la información del
elemento seleccionado.
45
Figura 3.19 Inspector mostrando la información requerida para el “Paro súbito de

bombeo” en la frontera aguas arriba
El Inspector consta de dos columnas, la primera es la descripción de la información solicitada y

la segunda corresponde a los valores proporcionados por el usuario. Si la cantidad de datos
requeridos por el elemento es tal que no cabe en el espacio del Inspector, entonces se presentará
a la derecha una barra de desplazamiento a fin de que el usuario la emplee para desplazarse
por toda la información.
Valores
Barra de
desplazamiento
Propiedades
sobre las
propiedades
Figura 3.20 Las tres principales zonas de trabajo del Inspector
46
Dependiendo del valor que se proporcione y el rango de posibles valores válidos para una
propiedad en particular se realiza una verificación, en su caso se emite un error en el área
de mensajes de la aplicación y se colorea de rojo la propiedad que generó dicho error.
3.7 Simulación y visualización de resultados

Una vez que se ha creado una propuesta se puede proceder a realizar una simulación.
Para ello se emplea el menú principal en la opción Simulación. Como primer paso del
proceso el sistema valida la información, en caso de encontrar algún error, éstos se
reportan en la zona de mensajes del programa (figura 3.21).
Figura 3.21 Errores generados al intentar simular una propuesta incorrecta (hay un tubo
que no está conectado)
Cuando el sistema no detecta errores en la propuesta se procede a la simulación, para lo

cual se solicita al usuario, mediante una ventana emergente (figura 3.22) el tiempo de
simulación y la frecuencia de grabación de resultados, en el capítulo 6 se trata este punto.
47
Figura 3.22 Ventana para capturar el tiempo de simulación y la frecuencia de grabación

de resultados una vez que se ha verificado la información de la propuesta
Una vez realizada la simulación, se presentan los resultados en el “área de dibujo y

resultados” mostrándose en la parte inferior una nueva pestaña denominada Resultados
Simulación esta pestaña permite visualizar a su vez las pestañas Reporte general y
Envolvente de presiones. La información mostrada en la pestaña Reporte general
corresponde en su totalidad a la información plasmada en el archivo con sufijo _IMP.
Figura 3.23 Resultados de la simulación, si la simulación es correcta, se activan dos

pestañas para visualizar el reporte general de la simulación y la envolvente de presiones
48
La pestaña Envolvente de presiones muestra las presiones máximas y mínimas a lo largo

de la conducción además de la piezométrica de operación. Para visualizarla se selecciona
esta pestaña (figura 3.24)
Figura 3.24 Resultados de la simulación, envolvente de presiones resultante
El perfil del terreno puede incluirse en la gráfica, para ello se presiona el botón Archivo del
perfil, eligiendo el archivo que contenga el perfil del terreno. Al seleccionar el archivo
correspondiente, se dibuja en la gráfica el perfil y se habilita la casilla Perfil de terreno
permitiendo mostrarlo u ocultarlo de la imagen (figura 3.25).
Para regresar al esquema de la propuesta, se selecciona la pestaña Diagrama, lo que

ocultará las pestañas Reporte general y Envolvente de presiones dejando visible la pestaña
Resultados simulación en la parte inferior del área de trabajo.
Si se realiza cualquier modificación en el esquema de la propuesta, invalidará los

resultados de la simulación, ocultando la pestaña Resultados simulación con la
consecuente necesidad de volver a ejecutar el proceso de simulación para visualizar
nuevamente los resultados.
49
Figura 3.25 Resultados de la simulación, envolvente de presiones resultante incluyendo

el perfil del terreno
50
ELEMENTOS DE UN SISTEMA DE BOMBEO QUE NO PRODUCEN TRANSITORIOS
4
Elementos de un sistema de
bombeo que no producen
transitorios
4.1 Introducción
En general, los sistemas de bombeo están formados por conducciones construidas con
tuberías de diferentes características, elementos hidráulicos instalados a lo largo de la
conducción y estructuras para proteger a la conducción de las presiones que se presentan
en su interior después de un golpe de ariete o suspensión súbita del bombeo.
En este capítulo se presentan los modelos matemáticos y cómo se simulan en el sistema

TransDocencia2-0, las tuberías, los cambios de tuberías, cámaras de aire, tanques de
carga constante, tanques o torres de oscilación, tanques unidireccionales y las válvulas de
admisión y expulsión de aire.
Cabe señalar que los modelos matemáticos que se resuelven en los softwares TRANS y
TransDocencia2-0 son los mismos.
El sistema TransDocencia2-0 valida los datos que el usuario le proporciona, los intervalos
de validez que se indican en las secciones de este capítulo se aplican a la versión del
sistema TRANS que se presenta en este manual.
51
Las estructuras y elementos que se presentan en este capítulo no generan transitorios, con
ellos no se pueden realizar maniobras para cambiar la forma de operar de los acueductos,
a excepción de las tuberías y cambios de tuberías, los demás están diseñados para
proteger a la conducción de las sobrepresiones y depresiones que se presentan en su
interior debido al golpe de ariete que se produce después de que se suspende súbitamente
el bombeo o después de un cierre o apertura de válvulas.
4.2 Tuberías
En el sistema de simulación de transitorios TransDocencia2-0 se consideran tuberías
elásticas y flujo compresible, la dinámica del flujo se representa con las ecuaciones de
golpe de ariete. A través de las ecuaciones características 2.41 y 2.42, que por
conveniencia en este capítulo se escribirán con la notación de las ecuaciones 4.1 y 4.2, se
obtienen la carga y el gasto en cada punto o sección de cálculo de las tuberías (ver figuras
2.19 y 4.1) en el instante 𝑡 + ∆𝑡.
𝑄1𝑡+∆𝑡 = 𝑐𝑝 − 𝑐𝑎 𝑖 𝐻1𝑡+∆𝑡 (4.1)
𝑄1𝑡+∆𝑡 = 𝑐𝑛 + 𝑐𝑎 𝑖 𝐻1𝑡+∆𝑡 (4.2)
Al sumar las ecuaciones anteriores se obtiene el gasto en el punto 1 señalado en la figura

4.1, ecuación 4.3, la carga se calcula sustituyéndolo en cualquiera de las ecuaciones 4.1 o
4.2.
𝑄1𝑡+∆𝑡 = 0.5(𝑐𝑝 + 𝑐𝑛 ) (4.3)
donde
𝐻1𝑡+∆𝑡 carga en el punto 1 en el tiempo 𝑡 + ∆𝑡, en m
𝑄1𝑡+∆𝑡 gasto en el punto 1 en el tiempo 𝑡 + ∆𝑡, en m3/s
𝑐𝑎𝑖 constante del método de las características del tubo 𝑖, en m2/s
Tramo 𝑖
𝑄0
Punto 1
Figura 4.1 Las secciones de la conducción se identifican a través de puntos de cálculos
52
4.2.1 Simulación en el sistema TransDocencia2-0. Para dar de alta un tramo

de tubería se selecciona la pestaña Tubos y elementos entre tramos y se elige el ícono
señalado en la figura 4.2. A través del Inspector se indican las propiedades de la tubería.
Figura 4.2 Pantalla para dar de alta a los tramos de las tuberías
Propiedades generales
 Nombre, con éste se identifica al tramo
Propiedades físicas
 Longitud [m], la longitud mínima que acepta el sistema TransDocencia2-0 es de 2 m.
 Diámetro [m], el diámetro mínimo permitido en TransDocencia2-0 es de 0.0508 m (2").
53
 Celeridad [m/s], su valor puede determinarse a través de la ecuación 2.25, el intervalo

permitido es 100 < 𝑎 < 1530 m/s.
 Coeficiente de fricción [D-W], es el factor de fricción de Darcy-Weisbach, puede

determinarse a través de la ecuación 2.11. El intervalo permitido es 0.008 ≤ 𝑓 ≤ 0.07.
 Gasto [m3/s], es el gasto en flujo permanente a través del tramo.
 Carga piezométrica al inicio del tramo [m], cuando en el tramo el flujo es cero debe
indicarse la carga en la frontera aguas arriba del tramo, cuando el flujo es diferente de
cero este campo se desactiva.
4.3 Cambio de tubería o tramo

Despreciando las pérdidas por velocidad en una región muy cercana a la unión de dos
tubos y la pérdida de carga en la junta, las cargas y los gastos en los puntos 1 y 2 de la
conducción de la figura 4.3 son iguales. Éstos se obtienen a partir de las ecuaciones
características 4.4 y 4.5, en donde los subíndices 𝑖 y 𝑖 + 1 se refieren a los tramos aguas
arriba y aguas abajo de la unión, los subíndices 1 y 2 identifican los puntos de la conducción
inmediatos (aguas arriba y aguas abajo) a la unión.
Observa que la constante 𝑐𝑎 𝑖 de la ecuación 4.4 debe calcularse con las propiedades del
tramo 𝑖 y la de la ecuación 4.5 con las propiedades del tramo 𝑖 + 1, las que a su vez se
determinan con la ecuación 2.32.
Al restar la primera ecuación de la segunda y despejar la carga, se obtiene la carga en la

unión de los tramos a través de la ecuación 4.6, los gastos se obtienen sustituyendo la
carga en cualquiera de las ecuaciones 4.4 o 4.5, los valores obtenidos de las variables
corresponden al tiempo 𝑡 + ∆𝑡, ver figura 4.3.
Tramo 𝑖
Tramo 𝑖 + 1
𝑄0
Punto 1 Punto 2
Figura 4.3 Cambio de tramo o de tubería
𝑄1𝑡+∆𝑡 = 𝑐𝑝 − 𝑐𝑎 𝑖 𝐻1𝑡+∆𝑡 (4.4)
𝑄2𝑡+∆𝑡 = 𝑐𝑛 + 𝑐𝑎 𝑖+1 𝐻2𝑡+∆𝑡 (4.5)
54
𝑐𝑝 − 𝑐𝑛
𝐻1𝑡+∆𝑡 = (4.6)
𝑐𝑎𝑖 + 𝑐𝑎𝑖+1
4.3.1 Simulación en el sistema TransDocencia2-0. Para dar de alta un cambio

de tubería se selecciona la pestaña Tubos y elementos entre tramos y se elige el ícono
señalado en la figura 4.4.
Figura 4.4 Cambio de tramo o de tubería
Cuando en el Inspector se palomea ¿Registro de carga y gasto en la unión de los tubos?,

al terminar una simulación, en el archivo de cargas, se graba la carga en mca y en el de
gastos, el gasto en m3/s en la conexión de los tubos, ambas variables en función del tiempo,
en s.
55
4.4 Cámara de aire

Las cámaras de aire son tanques a presión que contienen en su parte inferior agua y en la
superior un gas comprimido, comúnmente aire. En éstas se almacena la energía que
permite reducir las fluctuaciones de presión transitorias en el interior de las tuberías
producidas por el disparo de equipos de bombeo en los acueductos.
Durante la operación en gasto establecido la piezométrica de trabajo pasa por arriba de la

cámara, no hay flujo entre la cámara y la conducción, figura 4.5.
aire
agua
Conducción
𝐻𝑐 𝐻1
𝑄𝑐 = 0
Nivel de referencia
Figura 4.5 En flujo estable la piezométrica de trabajo pasa por arriba de la cámara
Cuando una o todas las bombas de una planta dejan de operar disminuye la presión del
agua en la conexión de la cámara con la tubería principal, por lo que el aire comprimido en
la cámara se expande y expulsa el agua que contiene hacia la tubería, figura 4.6.
Comúnmente el sistema de bombeo entrega el agua en un tanque, durante el tiempo en el

que se desarrolla un transitorio de golpe de ariete puede considerarse que la carga en el
tanque permanece constante, lo que impone en el extremo aguas abajo de la conducción
una sección con carga constante. La onda de baja presión producida por la suspensión de
bombeo se refleja en este tanque y produce un aumento de la presión en toda la
conducción conforme la recorre en dirección aguas arriba (ver figuras 2.10 a 2.18a). Al
llegar la onda de presión a la cámara, la diferencia de presiones que se establece entre la
línea principal y el interior de la cámara induce un flujo de agua hacia el interior de la
cámara; como resultado, el aire se comprime nuevamente (figura 4.7) pudiendo alcanzar
presiones suficientemente altas (superiores a la carga en el tanque de entrega) para volver
a empujar el agua contenida en la cámara hacia la línea principal.
56
Piezométrica en la
cámara
Piezométrica en
la conducción
aire
agua Conducción
𝐻𝑐 𝐻1
𝑄𝑐
Nivel de referencia
Figura 4.6 Cuando la presión en la conducción en el sitio de conexión con la cámara es

menor que la presión en el interior de la cámara, el aire en la cámara se expande y ésta
le entrega agua a la conducción
Piezométrica en
la conducción
aire Piezométrica
en la cámara
agua
𝐻1
Conducción
𝐻𝑐
𝑄𝑐
Nivel de referencia
Figura 4.7 Cuando la presión en la conducción en el sitio de conexión con la cámara es

mayor que la presión en el interior de la cámara, la cámara de aire recibe agua de la
conducción y el volumen de aire se comprime
57
El aire se comprime y se expande como consecuencia del viaje de ida y vuelta de la onda
de presión entre los extremos de la conducción hasta que se disipa la energía y se alcanza
una condición de equilibrio de cargas en toda la tubería.
Las cámaras de aire protegen bien a la tubería cuando contienen el volumen de aire
determinado en su diseño. El volumen de aire se mide con mirillas que se instalan a lo largo
de la cámara y que permiten registrar la altura de la interfaz agua-aire en el interior de la
cámara. El volumen de aire no puede medirse con manómetros que registren la presión en
el interior de la cámara puesto que la presión queda establecida por las condiciones de
flujo en la tubería principal, en la figura 4.7 se señala la mirilla de una cámara de aire.
mirilla
Figura 4.8 Cámara de aire con su mirilla para verificar el volumen de aire que contiene
(Grupo de Hidromecánica, IIUNAM)
El aire de las cámaras no se pierde cuando operan para controlar un transitorio por corte
de bombeo, sólo se pierde por fugas a través de conexiones mal realizadas y, muy
lentamente, por mezclarse con el agua del interior de la cámara, el volumen de aire que se
pierde se suministra con compresores. La capacidad de los compresores para mantener el
volumen de aire especificado en el diseño sólo tomará en cuenta esas pequeñas pérdidas,
mientras que, para el llenado de las cámaras por primera vez, o después de un
mantenimiento o reparación que haya requerido su vaciado, podrán usarse compresores
grandes que no requieren instalación permanente junto a las cámaras.
58
Un nivel de agua por arriba del indicado, significa falta de aire; el cual debe inyectarse con
el compresor. Cuando el nivel del agua esté por debajo de la referencia determinada se
tendrá exceso de aire que deberá expulsarse a través de una válvula de purga
especialmente instalada para ese propósito, con el fin de evitar que éste entre a la tubería
en caso de presentarse un corte de bombeo.
La conexión de las cámaras de aire con la conducción se realiza a través de arreglos que
contienen válvulas de seccionamiento necesarias para realizar mantenimiento y
reparaciones en las cámaras y en muchos casos, orificios diferenciales con pérdida de
carga mayores para el flujo que regresan de la tubería principal a la cámara y menores
para los que se salen de la cámara hacia la tubería.
En la figura 4.9 se muestra una de las cámaras y el arreglo de válvulas y placas orificio en
la tubería de conexión con la línea principal del acueducto Chapala Guadalajara. El arreglo
de la conexión no siempre es como el este acueducto, para definirlo se necesitan realizar
simulaciones del comportamiento de la cámara en flujo transitorio, más adelante se trata
este tema.
De la cámara hacia la conducción con

válvulas de no retorno
De la conducción hacia la cámara de

aire con orificios diferenciales
(a) (b)
Figura 4.9 (a) Cámara de aire (b) arreglo de conexión con la línea principal (Grupo de
Hidromecánica, IIUNAM)
4.4.1 Modelo matemático de cámara de aire. El modelo matemático está

formado por las ecuaciones de continuidad en la intersección de la cámara con la tubería,
la dinámica del volumen de agua en la cámara, del volumen de aire, de la tubería que
conecta a la cámara con la línea principal y las ecuaciones características en la tubería.
59
En la unión de la cámara con la tubería se debe satisfacer la ecuación de continuidad 4.7,

el sentido de los flujos se indica en la figura 4.10, considerando que la pérdida de carga es
muy pequeña en la interconexión de la cámara con la conducción, se tiene la ecuación 4.8.
𝑄1 + 𝑄𝐶 = 𝑄2 (4.7)
𝐻1 = 𝐻2 (4.8)
donde
𝑄𝐶 gasto que entrega o recibe la cámara, en m³/s
𝑄1 gasto en la conducción aguas arriba de la cámara, en m³/s
𝑄2 gasto en la conducción aguas abajo de la cámara, en m³/s
𝐻1 carga en la conducción inmediatamente aguas arriba del sitio de conexión de la
cámara, en m
𝐻2 carga en la conducción inmediatamente aguas abajo del sitio de conexión de la
cámara, en m
Conducción
𝑄2
Tramo 𝑖 + 1
Cámara de aire
𝑄𝐶 > 0
Punto 2
𝐴𝐶
𝑙
Punto 1
𝑄1 Tramo 𝑖
Figura 4.10 Conexión de la cámara de aire con la tubería principal
Considerando que el flujo entre la cámara y la tubería es como el de una columna rígida su
ecuación dinámica es la 4.9, ver figura 4.11.
𝑙 𝑑𝑄𝐶 𝑙 𝑄𝐶 |𝑄𝐶 |
𝐻𝑐 − 𝐻1 = +𝑓 + 𝐾𝑉 𝑄𝐶 |𝑄𝐶 | (4.9)
60
donde
𝐻𝑐 carga en la cámara de aire, en m
𝐾𝑉 coeficiente de pérdida de entrada o salida del agua de la cámara, en s²/m5
𝐴 área de la tubería que conecta a la cámara con la conducción, en m²
𝐷 diámetro de la tubería que conecta a la cámara con la conducción, en m
𝑙 longitud de la tubería que conecta a la cámara con la conducción, en m
𝑓 coeficiente de fricción de Darcy-Weisbach de la tubería que conecta a la cámara con
la conducción
𝑔 aceleración de la gravedad, en m/s²
Piezométrica en la
cámara
Piezométrica en
la conducción
aire
agua
Conducción
𝐻𝑐 𝑍𝑐 Z 𝐻1
𝑄𝑐 > 0 𝐷
Nivel de referencia
Figura 4.11 En la tubería de conexión de la cámara de aire con la conducción se

considera el modelo de columna rígida
La variación del nivel del agua en la cámara, 𝑍, depende de su área transversal, 𝐴𝐶 , y del
gasto de agua que entra o sale de la cámara, como indica la ecuación 4.10.
𝑑𝑍 𝑄𝐶
= − (4.10)
𝑑𝑡 𝐴𝐶
La presión absoluta del aire está dada por la ecuación 4.11.
𝑝 = 𝑝𝑎 + 𝛾(𝐻𝑐 − 𝑍) (4.11)
donde
𝑝 presión absoluta del aire en la cámara, en Pa
𝑝𝑎 presión atmosférica, en Pa
𝛾 peso específico del agua, en N/m³
𝑍 elevación del agua en la cámara, en msnm
61
Se ha visto que el proceso de expansión y compresión del aire en el interior de la cámara,

se representa con muy buena aproximación considerando al aire como un gas perfecto,
por lo que, 𝑝𝑉 𝛼 = 𝐶 ∗ donde 𝑝 es la presión absoluta del aire, 𝑉 su volumen, 𝐶 ∗ una
constante y 𝛼 puede considerarse constante con un valor entre 1.0 y 1.4. El volumen de
aire en la cámara es 𝑉 = 𝐴𝐶 (𝑍𝐶 − 𝑍) debido a que el área es constante (ver figura 4.11),
definiendo 𝐶 = 𝐴𝐶 −𝛼 𝐶 ∗, la ecuación constitutiva del aire puede escribirse como la ecuación
4.12. En el sistema TransDocencia2-0 se considera 𝛼 = 1.2 debido a los buenos resultados
que se han obtenido al comparar datos obtenidos en campañas de medición de transitorios
en prototipo con resultados de simulaciones numéricas.
𝑝(𝑍𝐶 − 𝑍)𝛼 = 𝐶 (4.12)

donde
𝑍𝑐 elevación del límite superior de la cámara de aire, en msnm
Las incógnitas del sistema formado por las ecuaciones 4.7 a 4.12 son 𝑄𝐶 , 𝑄1 , 𝑄2 , 𝐻1 , 𝐻2 ,
𝐻𝑐 , 𝑝 y 𝑍, se tienen seis ecuaciones con ocho incógnitas. El sistema se completa con el
modelo de golpe de ariete en la conexión de la cámara con la tubería a través de las
ecuaciones características 4.4 y 4.5.
El sistema se resuelve numéricamente, para ello, las ecuaciones 4.7 a 4.12 se escribe en
diferencias finitas. El sistema completo está formado por las ecuaciones 4.4, 4.5, 4.13 a
4.18. Ahora la nueva nomenclatura de las incógnitas es 𝑄𝐶𝑡+∆𝑡 , 𝑄1𝑡+∆𝑡 , 𝑄2𝑡+∆𝑡 , 𝐻1𝑡+∆𝑡 , 𝐻2𝑡+∆𝑡 ,
𝐻𝐶𝑡+∆𝑡 , 𝑝𝑡+∆𝑡 y 𝑍𝑡+∆𝑡 . Los subíndices 𝑡 y 𝑡 + ∆𝑡 indican el valor de la variable
correspondiente en el tiempo 𝑡 y en el tiempo 𝑡 + ∆𝑡, respectivamente.
𝑄2 𝑡+∆𝑡 = 𝑄1 𝑡+∆𝑡 + 𝑄𝐶𝑡+∆𝑡 (4.13)
𝐻1 𝑡+∆𝑡 = 𝐻2 𝑡+∆𝑡 (4.14)
𝑙 𝑄𝐶𝑡+∆𝑡 − 𝑄𝐶𝑡 𝑓𝑙
𝐻𝐶𝑡+∆𝑡 − 𝐻1 𝑡+∆𝑡 = + ( + 𝐾𝑉 ) 𝑄𝐶𝑡+∆𝑡 |𝑄𝐶𝑡+∆𝑡 | (4.15)
𝑔𝐴 ∆𝑡 2𝑔𝐷𝐴2
𝑍𝑡+∆𝑡 − 𝑍𝑡 1 𝑄𝐶𝑡+∆𝑡 + 𝑄𝐶𝑡

= − (4.16)
∆𝑡 𝐴𝐶 2
𝑝𝑡+ ∆𝑡 = 𝑝𝑎 + 𝛾(𝐻𝐶𝑡+∆𝑡 − 𝑍𝑡+∆𝑡 ) (4.17)
𝑝𝑡+ ∆𝑡 (𝑍𝐶 − 𝑍𝑡+∆𝑡 )1.2 = 𝐶 (4.18)
A través del método de sustitución, el sistema anterior puede reducirse a una ecuación no
lineal en 𝑄𝐶𝑡+∆𝑡 , la cual puede resolverse con el método de Newton-Raphson.
62
Cuando 𝑄𝐶𝑡+∆𝑡 > 0 significa que la cámara se está vaciando y debe considerarse el
coeficiente de pérdida del orificio de salida de agua. Si 𝑄𝐶𝑡+∆𝑡 < 0, significa que está
entrando agua a la cámara y debe tomarse en cuenta el coeficiente de pérdida 𝐾𝑉 de
entrada a la cámara.
Una vez obtenido el gasto 𝑄𝐶 𝑡+∆𝑡 , la elevación del agua en la cámara 𝑍𝑡+∆𝑡 se obtiene con
la ecuación 4.16, la presión en el interior de la cámara 𝑝𝑡+∆𝑡 con la ecuación 4.18 (la
constante 𝐶 se determina a partir de las condiciones de operación en flujo estable), la carga
piezométrica en la cámara 𝐻𝐶 𝑡+∆𝑡 con la ecuación 4.17, las cargas en los puntos 1 y 2 (ver
figura 4.10) con las ecuaciones 4.15 y 4.14, los gastos en esas secciones de la conducción
con las ecuaciones características 4.4 y 4.5.
4.4.2 Predimensionamiento de cámaras de aire. El volumen de aire, la

elevación del agua, el área transversal de la cámara y los coeficientes de pérdida de carga
del flujo que entra o sale de la cámara son los parámetros que caracterizan una cámara de
aire, para predeterminarlos se han propuestos varios métodos, como los presentados por
Guarga et al. (1985), Thorley (1991) y Carmona et al. (2002).
A continuación, se presenta el algoritmo propuesto en (Carmona et al., 2002), con el cual

se obtiene un valor preliminar del volumen de aire. El método se basa en la hipótesis de
que la presencia de cámaras de aire reduce significativamente la rapidez de variación de
las condiciones de frontera producidas por corte de bombeo, de tal manera que el
comportamiento transitorio puede representarse con el modelo de oscilación de masa o de
transitorio lento. Para su aplicación se debe proporcionar la siguiente información, ver
figura 4.12.
𝐻1 𝑍0 𝐻𝑚í𝑛 𝐻2
Nivel de referencia
𝐿
Figura 4.12 Datos para calcular un volumen de aire preliminar, Carmona et al. (2002)
63
1. Carga piezométrica en la descarga de las bombas en flujo estable, 𝐻1 en m.
2. Carga piezométrica en la descarga del tramo de conducción, 𝐻2 en m.
3. Carga piezométrica mínima permitida en algún punto alto del tramo de la conducción
que protegerá la cámara de aire, 𝐻𝑚í𝑛 en m.
4. Nivel inicial del agua en la cámara de aire, 𝑍0 en msnm. De manera preliminar, podría
considerarse un tirante de 2.0 a 3.0 m.
5. Longitud del tramo de la conducción que protegerá la cámara, 𝐿 en m.
6. Área transversal de la tubería principal, 𝐴𝑡 en m².
7. Gasto en flujo estable, 𝑄0 en m³/s.
Con los datos anteriores se calculan los siguientes parámetros
𝑝1 = (𝐻1 − 𝑍0 + 10.33)(9 810)
𝑝2 = (𝐻2 − 𝑍0 + 10.33)(9 810)
𝐻𝑚í𝑛 − 𝐻2
𝑧𝑚í𝑛 =
𝐻1 − 𝐻2
9 810 𝐿 2
𝐸𝐾 = 𝑄0
2𝑔𝐴𝑡
𝑝1
𝑟= ⁄𝑝2
Calcular T ∗ con la ecuación

∗ 𝜋
𝑇 ∗𝑒 𝑇 + =0
2 𝑧𝑚í𝑛
se puede utilizar la función solver o función objetivo de Excel.
Una vez determinado el valor de T ∗, calcular:
𝜋 2
1 + (2𝑇 ∗ ) 1
𝐾= 1⁄ −1
1 − 𝑧𝑚í𝑛 1 − 𝑧𝑚í𝑛 1.2
{(𝑧𝑚í𝑛 + 𝑟 ) }
64
Una vez realizados los cálculos indicados, el volumen preliminar de aire en la cámara se
obtiene con la ecuación 4.19.
2 𝐸𝐾
𝑉0 = (4.19)
1
𝑝1 𝐾 (1 − )
𝑟
Se debe poner atención en la variación del nivel del agua en la cámara, pues es
recomendable que ésta no se vacíe para evitar que entre aire a la tubería. La variación del
nivel del agua en la cámara depende de su área, para determinarla se han propuesto varios
métodos, uno de ellos se presenta en Guarga et al. (1985) en el que se propone que la
altura de la columna de aire, antes de que la cámara comience a operar sea igual al
diámetro de la cámara (𝐷𝐶 ). Para esa geometría, el diámetro de la cámara se obtiene con
la ecuación 4.20.
1⁄
4𝑉0 3
𝐷𝐶 = ( ) (4.20)
𝜋
donde
𝐷𝐶 diámetro preliminar de la cámara, en m
𝑉0 volumen inicial de aire preliminar, en m3
4.4.3 Simulación en el sistema TransDocencia2-0. Para dar de alta una

cámara de aire se selecciona la pestaña Tubos y elementos entre tramos y se elige el ícono
señalado en la figura 4.13, a través del Inspector se proporcionan los parámetros que
caracterizan a la cámara, el cual contiene las secciones Propiedades generales y
Propiedades físicas.
 Nombre, con éste se identifica a la cámara de aire
 Coeficiente de pérdida a la salida de la cámara [s2/m5], en la etapa de

predimensionamiento puede considerarse un coeficiente de pérdida de carga del flujo
que sale de la cámara igual a 0.01 s2/m5.
 Volumen de aire inicial en la cámara [m3], el volumen de aire preliminar se obtiene con
la ecuación 4.19, si este volumen es insuficiente para controlar el transitorio deberá
duplicarse o triplicarse y realizar nuevas simulaciones.
 Área de la cámara [m2], el área de la cámara es el que corresponde al diámetro que se

obtiene con la ecuación 4.20, cuando se opta por la geometría propuesta con la que se
obtuvo esa ecuación.
65
 Cota del agua en la cámara [msnm], la elevación de la superficie libre del agua en la
cámara se determina a partir del sitio en donde se construirá. Un valor preliminar puede
obtenerse considerando que la base de la cámara (figura 4.14) este entre 1.5 y 2.0 m
arriba de la elevación del terreno en donde se construirá, a partir de la elevación de la
base de la cámara pueden considerarse de 2.0 a 3.0 m de columna de agua, como
resultado se obtiene una elevación preliminar de la superficie libre del agua.
Figura 4.13 Pantalla para dar de alta una cámara de aire
 Elevación base la cámara [msnm], la elevación de la base de la cámara se determina a

partir de la elevación del terreno del sitio en donde se construirá, un valor preliminar
puede ser entre 1.5 y 2.0 m arriba de esa elevación, figura 4.14.
 Coeficiente de pérdida a la entrada de la cámara [s2/m5], en la etapa de

predimensionamiento puede considerarse un coeficiente de pérdida de carga del flujo
que retorna a la cámara igual a 0.01 s2/m5, más adelante se presenta cómo puede
determinarse el valor definitivo de este coeficiente.
66
 Longitud de la tubería de conexión cám-conducción [m], (ver figuras 4.10 y 4.11) en etapas
de pre diseño puede considerase de unos cinco metros, la longitud definitiva depende
del sitio de construcción.
 Diámetro de la tubería de conexión cám-conducción [m], (ver figuras 4.10 y 4.11) en etapas
de pre diseño puede considerarse igual o menor al de la conducción.
Aire
Aire
Agua
Agua
Base de la cámara
Figura 4.14 Durante la operación de la cámara, se recomienda que el nivel del agua
siempre esté por arriba de la base de la cámara
Para etapas posteriores a la de pre-diseño o factibilidad, la longitud de la tubería de

conexión se determina a partir del sitio en el que se instalarán las cámaras y del arreglo de
válvulas necesario (figuras 4.9 y 4.20) y su diámetro de tal forma que la velocidad del flujo
en las tuberías de conexión entre las cámaras y la línea de conducción sea menor a 4 m/s.
Al terminar una simulación, en el archivo de cargas se graba la carga en mca en la conexión

de la cámara con la línea de conducción; en el de gastos, el gasto en m 3/s que entrega o
recibe la cámara, positivo cuando la cámara entrega agua a la conducción y negativo
cuando el agua regresa a la cámara; en el archivo con sufijo rpm, se graba la elevación del
agua en la cámara, en msnm; todas las variables están en función del tiempo, en s.
En la figura 4.15 se presentan las líneas de carga piezométricas en flujo estable y las
máximas y mínimas obtenidas al simular dos transitorios generados por la suspensión
súbita del bombeo y el cierre de las válvulas en la descarga de las bombas en un sistema
similar al de la figura 4.12, la planta de bombeo se localiza en el cadenamiento 0+000. En
una simulación no hay cámara de aire y en la otro opera una cámara de aire considerando
ambos coeficientes de pérdida de carga del flujo que entra o sale de la cámara iguales a
0.01 s2/m5. Se observa que cuando se instala una cámara de aire con un volumen de aire
suficiente se evitan las subpresiones pero, en general, las sobrepresiones son mayores
que cuando no hay cámara de aire.
67
240
Piezométrica y elevación, en m 180
120
60
0
0+000 0+200 0+400 0+600 0+800 1+000 1+200 1+400
Cadenamiento, en km
Terreno Flujo estable Sin cámara Con cámara, Ke=0.01
Figura 4.15 Piezométrica en flujo estable y piezométricas máximas y mínimas en flujo

transitorio con y sin cámara de aire, 𝐾𝑒 en s2/m5
Lo anterior ocurre en los sistemas de bombeo con tanques de entrega desplantados en

sitios más elevados que en los que se construyen las cámaras de aire (generalmente las
cámaras de aire se instalan muy cerca de las plantas de bombeo, las que a su vez, en
general están desplantadas más abajo que los tanques de entrega), esto provoca que al
suspenderse el bombeo el agua se regrese del tanque de entrega hacia las cámaras, el
flujo de agua que entra a las cámaras comprime el aire en la cámara produciendo un
aumento de la presión en su interior y en la conducción, como se muestra en la figura 4.15.
Las sobrepresiones se evitan colocando disipadores de energía, como placas orificio, en

la línea que lleva el agua de regreso hacia la cámara, ver figura 4.20. En la figura 4.16 se
comparan las líneas de carga piezométrica que resultan al considerar el coeficiente de
pérdida de carga (𝐾𝑒 ) de 0.01, 10 y 100 s2/m5 para el flujo que regresa a la cámara. Se
observa que cuando 𝐾𝑒 = 100 s2/m5, la piezométrica de cargas máximas se aproxima
mucho a la de operación en flujo estable.
Al observar las gráficas de la figura 4.16 podría pensarse que a medida que se aumenta el
coeficiente de pérdida de carga del flujo que llena la cámara la sobrepresión disminuye,
pero no es así, como se observa en la figura 4.17. Lo que sucede es que, con pérdidas
muy grandes, cuando el agua regresa de la conducción hacia la cámara su paso está
obstruido, la velocidad del agua disminuye y la presión aumenta. Una vez determinado el
volumen de aire en la cámara, existe un coeficiente de pérdida óptimo para el cual la
sobrepresión se acerca mucho a la piezométrica de flujo estable, este puede determinarse
realizando simulaciones de transitorios con diferentes coeficientes.
68
240
Piezométrica y elevación, en m 180
120
60
0
0+000 0+200 0+400 0+600 0+800 1+000 1+200 1+400
Cadenamiento, en km
Terreno Flujo estable Ke=0.01 Ke=10 Ke=100
Figura 4.16 Sobrepresión ocasionada por la compresión del aire en la cámara para
diferentes coeficientes de pérdida de carga del flujo que ingresa a la cámara, 𝐾𝑒 en s2/m5
240
Piezométrica y elevación, en m
180
120
60
0
0+000 0+200 0+400 0+600 0+800 1+000 1+200 1+400
Cadenamiento, en km
Terreno Flujo estable Ke=100 Ke=10000
Figura 4.17 Con coeficientes de pérdida de carga muy grandes para el flujo que regresa
a la cámara se producen sobrepresiones por golpe de ariete, 𝐾𝑒 en s2/m5
69
Para determina el coeficiente óptimo se recomienda realizar simulaciones con coeficientes

mayores al que se obtiene con la ecuación 4.21.
𝐻𝑚á𝑥 − 𝐻0
𝐾𝑒 = 2 (4.21)
(𝑄𝑒𝑚á𝑥 )
donde
𝐾𝑒 coeficiente de pérdida de carga del flujo que regresa a la cámara, en s2/m5
𝐻0 carga en flujo estable en la conexión de la cámara con la conducción, en m
𝐻𝑚á𝑥 carga máxima en la conexión de la cámara con la conducción obtenida de la
simulación en la que 𝐾𝑒 = 0.01 s2/m5, en m
𝑄𝑒𝑚á𝑥 gasto máximo que entra a la cámara obtenido de la simulación en la que 𝐾𝑒 = 0.01
s2/m5, en m3/s
Ejemplo. En la simulación de un paro súbito del bombeo en un acueducto, considerando

una cámara de aire con un coeficiente de pérdida de carga del flujo que entra a la cámara
𝐾𝑒 = 0.01 s2/m5, la carga en el sitio de conexión de la cámara con la conducción y el gasto
que entrega la cámara se comportan como se muestra en la figura 4.18. Se observa una
carga piezométrica máxima de 𝐻𝑚á𝑥 = 201.476 m y un gasto máximo que entra a la cámara
(negativo en la figura 4.18) de |𝑄𝑒𝑚á𝑥 | = 1.7 m3/s.
240 2
180 1
Piezométrica, en mca
Gasto, en m3/s
120 0
60 -1
0 -2
0 10 20 30 40 50
Tiempo, en s
Carga Gasto
Figura 4.18 Carga piezométrica en la conexión de la cámara con la conducción y gasto

que entrega (positivo) y recibe la cámara (negativo) cuando 𝐾𝑒 = 0.01 s2/m5
70
Al aplicar la ecuación 4.21 se obtiene:
𝐾𝑒 = (201.476 − 85.675)⁄(1.7)2 = 40 s2 /m5
Entonces, para obtener el coeficiente óptimo de pérdida de carga del flujo que retorna a la
cámara se realizan simulaciones con coeficientes 𝐾𝑒 > 40 s2/m5 en las que tanto el volumen
de aire inicial como el coeficiente de pérdida de carga del flujo que sale de la cámara no
deben cambiarse.
En la figura 4.19 se muestra la carga piezométrica máxima en la conexión de la cámara

con la conducción en función del coeficiente de pérdida de carga del flujo que entra a la
cámara. En este ejemplo el coeficiente óptimo puede considerarse igual a 100 s2/m5.
240
Carga piezométrica, en m
180
120
60
0
0 100 200 300 400 500 600 700 800 900 1000
Coeficiente de pérdida de carga, en s2/m5
Figura 4.19 En este caso 𝐾𝑒 = 100 s2/m5 es el valor óptimo
Para el diseño final de una cámara de aire, se deben tomar en cuenta los siguientes puntos.
La longitud y el diámetro de la tubería de conexión de la cámara con la línea de conducción

y el desplante de la cámara dependerán del sitio donde se puede construir la cámara. Éstos
deberán considerarse al realizar el análisis definitivo del funcionamiento de la cámara.
Para evitar flujo de agua de la línea de principal de conducción hacia la cámara durante la
operación a gasto establecido, deben instalarse válvulas check (retención o no retorno) en
la conexión de la cámara con la línea principal. Estas válvulas producen pérdidas a la salida
de la cámara que deben tomarse en cuenta en el momento de realizar análisis de
transitorios hidráulicos, figura 4.20.
71
Mirilla
Válvula de retención
𝑄𝐶
𝑄𝐶
Placas orificio
Figura 4.20 Esquema de la conexión de dos cámaras de aire con la línea principal
(Carmona et al., 1999)
La altura de la cámara la determinan las longitudes de las columnas de agua y aire,

dependiendo del diámetro de la cámara y la geometría de su tapa, ésta última puede
contener un volumen considerable de aire que debe tomarse en cuenta al determinar la
altura total de la cámara.
Debido a que puede haber pérdida del volumen de aire por fugas o disolución del aire en
el agua, el nivel del agua aumenta. Si al reponer el aire perdido se hubiera inyectado con
un compresor más aire del necesario, el nivel del agua quedará por debajo del proyectado.
Lo anterior indica que el volumen de aire y la columna de agua varían, por lo que deberá
determinarse un intervalo de la columna de agua (𝑍𝑚á𝑥 , 𝑍𝑚í𝑛 ), ver figura 4.21, que garantice
que cuando la interface agua-aire esté dentro de ese intervalo la cámara protegerá al
acueducto de las presiones generadas durante un golpe de ariete y al expandirse el aire
éste no ingrese a la conducción. Generalmente, la elevación de diseño de la interface agua-
aire (𝑍0 ) está a la mitad de ese intervalo. El nivel del agua en la cámara se revisa través de
una mirilla, figuras 4.8 y 4.21.
72
𝑍𝑚á𝑥
𝑍0
𝑍𝑚í𝑛
Mirilla
Figura 4.21 Se debe determinar un intervalo de la elevación de la interface agua-aire que

asegure la protección del sistema de bombeo
Cuando se necesita un volumen de aire muy grande para controlar los transitorios de tal
forma que la cámara resulta ser muy grande, éste se puede repartir en varias cámaras, con
las ventajas y desventajas que conlleva.
Ventajas: se tienen varias cámaras más pequeñas que pudieran proporcionar mayor
versatilidad a la operación del sistema.
Desventajas: más instalaciones de válvulas, más compresores, estar al cuidado y dar

mantenimiento a más estructuras hidráulicas.
Al repartir el volumen de aire en partes iguales en 𝑛 cámaras, el diámetro de cada una se

obtiene con la ecuación 4.22.
𝐷𝐶
𝐷𝑖 = (4.22)
√𝑛
donde
𝐷𝐶 diámetro de la cámara grande
𝐷𝑖 diámetro de las cámaras pequeñas
𝑛 número de cámaras
Generalmente las cámaras de aire se construyen en el mismo predio de la planta de

bombeo, las mismas personas que dan mantenimiento a la planta también se lo dan a las
cámaras.
En la figura 4.22 se muestran algunas cámaras de aire construidas en acueductos

mexicanos.
73
Figura 4.22 El tamaño y número de las cámaras de aire depende de las características
de cada sistema de bombeo (Grupo de Hidromecánica, IIUNAM)
4.5 Tanques de carga constante

Es muy común que en los extremos de las conducciones de los sistemas de bombeo
existan tanques cuyos niveles de agua pueden considerarse constantes durante el tiempo
en que se desarrolla un transitorio rápido, lo que impone fronteras de carga constante.
A continuación, se presentan los modelos matemáticos para los tanques de carga

constante localizados en la frontera aguas arriba del sistema y para los ubicados en la
frontera aguas abajo. En TransDocencia2-0 se considera que la carga de velocidad en la
entrada de la tubería, para el primer caso y a la salida de la tubería para el segundo, es
muy pequeña comparada con la carga de presión impuesta por la elevación de la superficie
74
libre del agua en el tanque, por lo que las pérdidas en la entrada y salida de la tubería no
se toman en consideración.
4.5.1 Tanque en la frontera aguas arriba. En la figura 4.23 se muestra un tanque

en la frontera aguas arriba de la conducción. La carga en la sección de la conducción
inmediata aguas abajo del tanque (punto 1 de la figura 4.32) es igual a la carga de presión
(𝐻𝑆 ) impuesta por el tanque, ecuación 4.23. El gasto se obtiene sustituyendo la carga en la
característica negativa, ecuación 4.2.
𝐻1𝑡+∆𝑡 = 𝐻𝑆 (4.23)
𝐻𝑆 Tramo 𝑖
𝑄0
Punto 1
Figura 4.23 Tanque de carga constante en el extremo aguas arriba de la conducción
4.5.2 Simulación en el sistema TransDocencia2-0. Para dar de alta un tanque

de carga constante en la frontera aguas arriba se elige la pestaña Fronteras y maniobras
aguas arriba (figura 3.16) y se selecciona el ícono del tanque señalado en la figura 4.24. A
través del Inspector se piden las propiedades generales del tanque.
 Nombre, con éste se identifica al tanque o embalse.
 Carga [msnm], carga de presión que impone el tanque o elevación de la superficie libre
del agua.
 Cadenamiento del embalse [m], es la ubicación del tanque respecto al cadenamiento

topográfico.
75
Al terminar una simulación, en el archivo de cargas se graba la carga inmediatamente

aguas abajo del tanque en mca, que en este caso es un valor constante y en el de gastos,
el gasto que entrega el tanque en m3/s, ambas variables en función del tiempo, en s.
Figura 4.24 Pantalla para dar de alta un tanque de carga constante en la frontera aguas
arriba
4.5.3 Tanque en la frontera aguas abajo. En la figura 4.25 se muestra un tanque

en la frontera aguas abajo de la conducción, se considera descarga ahogada. La carga en
la sección de la conducción inmediata aguas arriba del tanque (punto 𝑛 en la figura 4.25)
es igual a la carga de presión (𝐻𝐸 ) impuesta por el tanque, ecuación 4.24. El gasto se
obtiene sustituyendo la carga en la característica positiva, ecuación 4.1.
𝐻𝑛𝑡+∆𝑡 = 𝐻𝐸 (4.24)
76
𝐻𝐸
Tramo 𝑖
𝑄0
Punto 𝑛
Figura 4.25 Tanque de carga constante en el extremo aguas abajo de la conducción

de carga constante en la frontera aguas abajo de un sistema de bombeo o tanque de
entrega se elige la pestaña Frontera y maniobras aguas abajo (figura 3.16) y se selecciona
el ícono del tanque que se muestra en la figura 4.26. A través del Inspector se pide como
propiedad general su nombre.
 Nombre, con éste se identifica al tanque o embalse.
La carga de presión que impone el tanque se calcula a través de la ecuación de la energía

a partir de la condición de operación en flujo estable y las características de la conducción.
Al terminar una simulación, en el archivo de cargas se graba la carga de la sección

inmediatamente aguas arriba del tanque, en mca, que en este caso es un valor constante.
En el archivo de gastos, el gasto en el tanque en m3/s, ambas variables en función del
tiempo, en s.
77
Figura 4.26 Pantalla para dar de alta un tanque de carga constante en el extremo aguas
abajo de la conducción
4.6 Tanque unidireccional

Un tanque unidireccional es una estructura que se utiliza para limitar las depresiones
durante un transitorio. Este dispositivo está formado por un tanque abierto en su parte
superior que se conecta a la conducción a través de una tubería en la que se instalan, por
un lado, válvulas de retención (no retorno o check) y por otro, válvulas que permiten el
llenado del tanque después haber operado.
Los tanques unidireccionales operan cuando la piezométrica en la tubería de conducción

es menor que la elevación del agua en el tanque. Cuando la situación anterior se presenta,
las válvulas de retención se abren y el agua contenida en el tanque fluye hacia la
conducción; manteniendo en una región cercana al tanque, la piezométrica en la tubería
principal con un valor similar al nivel del agua en el tanque (siempre y cuando las pérdidas
en la conexión entre el tanque y la línea principal sean pequeñas). Un tanque con las
78
dimensiones y ubicación adecuadas, evita que se presenten depresiones en la conducción

en la cercanía de la conexión del tanque con la tubería principal.
El sentido correcto de la instalación de las válvulas de retención es el que impide que

durante la operación en flujo permanente el agua fluya de la conducción hacia el tanque
unidireccional, de tal forma que el nivel del agua en el tanque puede estar por debajo de la
piezométrica de operación del sitio en donde se conecta el tanque con la conducción.
Una vez que el tanque ha operado, el nivel del agua debe recuperarse para que el tanque
pueda volver a funcionar en cualquier momento, para ello se instalan las válvulas de
llenado. Éstas trabajan mecánicamente a través de un flotador o con un presóstato fijado
a la pared del tanque que abre la válvula cuando la presión en la pared está por debajo de
la correspondiente a la de tanque lleno. Cuando el nivel del agua en el tanque ha bajado
hasta el programado para abrir las válvulas de llenado y la carga de presión en el interior
de la tubería es mayor que ese nivel, se crea un flujo de agua de la línea de conducción
hacia el tanque que volverá a llenarlo. En la figura 4.27 se muestra un tanque unidireccional
con su arreglo de válvulas y la línea que lo conecta con la conducción principal.
Válvulas de no retorno
Válvulas de llenado
Conexión con la línea principal
Figura 4.27 Tanque unidireccional y el arreglo de válvulas para conectarlo a la línea

principal (Grupo de Hidromecánica, IIUNAM)
79
4.6.1 Modelo matemático de tanque unidireccional. El modelo matemático

está formado por la ecuación de continuidad en la intersección de la conducción principal
con la tubería que conecta al tanque, las ecuaciones dinámicas del volumen de agua en el
tanque y de la tubería que lo conecta con la línea principal y las ecuaciones características
en la tubería.
Considerando el sentido del flujo del tanque hacia la línea principal como positivo y
despreciable la pérdida de carga en la conexión entre la conducción principal y la tubería
que conecta al tanque, se tienen las ecuaciones 4.25 y 4.26, ver figura 4.28.
Tanque
unidireccional
Válvula de
retención
𝑄0 𝑄𝑇 > 0
𝑄𝑇 < 0
𝑙
Válvula de llenado
Figura 4.28 Conexión de tanques unidireccionales con la conducción
𝑄1 + 𝑄𝑇 = 𝑄2 (4.25)
𝐻1 = 𝐻2 (4.26)
Considerando el flujo en la tubería que conecta al tanque con la conducción como el de

una columna rígida, su ecuación dinámica es la 4.27, ver figura 4.28.
𝑙 𝑑𝑄𝑇 𝑙 𝑄𝑇 |𝑄𝑇 |
𝑍 − 𝐻1 = +𝑓 + 𝐾𝑣 𝑄𝑇 |𝑄𝑇 | (4.27)
80
Despreciando la inercia del agua en el tanque y constante la sección transversal del tanque,
la dinámica del agua en el mismo está representada por la ecuación 4.28.
𝑑𝑍 𝑄𝑇
= − (4.28)
𝑑𝑡 𝐴𝑇
donde
𝑄𝑇 gasto que entrega (positivo) o recibe el tanque (negativo), en m³/s
𝐴 𝑇 área transversal del tanque, en m2
𝑍 tirante de agua en el tanque, en m
𝐴 área de la tubería que conecta al tanque con la conducción, en m²
𝐷 diámetro de la tubería que conecta al tanque con la conducción, en m
𝑙 longitud de la tubería que conecta al tanque con la conducción, en m
Las incógnitas del sistema formado por las ecuaciones 4.25 a 4.28 son 𝑄𝑇 , 𝑄1 , 𝑄2 , 𝐻1 , 𝐻2
y 𝑍, se tienen cuatro ecuaciones con seis incógnitas, el sistema se completa con las
ecuaciones características 2.41 y 2.42.
El sistema de ecuaciones se resuelve numéricamente transformándolo en el sistema de

ecuaciones en diferencias finitas determinado por las ecuaciones 4.4, 4.5, 4.29 a 4.32. Los
subíndices 𝑡 y 𝑡 + ∆𝑡 indican el valor de la variable correspondiente en el instante de tiempo
𝑡 y 𝑡 + ∆𝑡, respectivamente. Ahora la nueva nomenclatura de las incógnitas es 𝑄𝑇𝑡+∆𝑡 ,
𝑄1𝑡+∆𝑡 , 𝑄2𝑡+∆𝑡 , 𝐻1𝑡+∆𝑡 , 𝐻2𝑡+∆𝑡 y 𝑍𝑡+∆𝑡 .
𝑄2 𝑡+∆𝑡 = 𝑄1 𝑡+∆𝑡 + 𝑄𝑇𝑡+∆𝑡 (4.29)
𝐻1 𝑡+∆𝑡 = 𝐻2 𝑡+∆𝑡 (4.30)
𝑙 𝑄𝑇𝑡+∆𝑡 − 𝑄𝑇𝑡 𝑓𝑙
𝑍𝑡+∆𝑡 − 𝐻1 𝑡+∆𝑡 = + ( + 𝐾𝑉 ) 𝑄𝑇𝑡+∆𝑡 |𝑄𝑇𝑡+∆𝑡 | (4.31)
𝑔𝐴 ∆𝑡 2𝑔𝐷𝐴2
𝑍𝑡+∆𝑡 − 𝑍𝑡 1 𝑄𝑇𝑡+∆𝑡 + 𝑄𝑇𝑡

= − (4.32)
∆𝑡 𝐴𝑇 2
Por el método de sustitución el sistema se reduce a la ecuación 4.33, ecuación de segundo

grado.
(𝑈 + 𝐾𝑉 ) 𝑄𝑇𝑡+∆𝑡 |𝑄𝑇𝑡+∆𝑡 | + 𝑉 𝑄𝑇𝑡+∆𝑡 + 𝑊 = 0 (4.33)
donde
𝑓𝑙
𝑈=
2𝑔𝐷𝐴2
81
𝑙 1 ∆𝑡
𝑉= + +
𝑔𝐴∆𝑡 𝑐𝑎1 + 𝑐𝑎2 𝐴 𝑇
𝑐𝑝 − 𝑐𝑛 𝑙
𝑊= − 𝑄 − 𝑍𝑡
𝑐𝑎1 + 𝑐𝑎2 𝑔𝐴∆𝑡 𝑇𝑡
Cuando 𝑄𝑇 𝑡+∆𝑡 > 0 el tanque entrega agua a la línea y 𝐾𝑉 corresponde al coeficiente de

pérdida de las válvulas de retención, cuando 𝑄𝑇 𝑡+∆𝑡 < 0 el agua entra al tanque, 𝐾𝑉 es el
coeficiente de pérdida de las válvulas de llenado.
Una vez obtenido el gasto 𝑄𝑇 𝑡+∆𝑡 , el tirante de agua en el tanque se obtiene con la ecuación
4.32, las cargas en los puntos 1 y 2 de la figura 4.28 con las ecuaciones 4.31 y 4.30, los
gastos en esas secciones de la conducción con las ecuaciones características 4.4 y 4.5.
4.6.2 Predimensionamiento de tanques unidireccionales. Las dimensiones

del tanque y el volumen de agua que contiene deben ser tales que durante su
funcionamiento no se vacíe para evitar la entrada de aire a la tubería, además, deberá estar
preparado para volver a operar inmediatamente previendo que el sistema se vuelva a poner
en operación antes de que recupere el volumen de agua a través de las válvulas de llenado.
El área transversal preliminar puede obtenerse con la ecuación 4.34

4𝐿 𝑄0
𝐴𝑇 = (4.34)
𝑎 𝐷
donde
𝐴𝑇 área transversal preliminar del tanque unidireccional, en m2
𝑄0 gasto en la conducción en flujo establecido, en m3/s
𝐷 diámetro de la conducción, en m
𝐿 longitud de la conducción, aguas abajo del tanque, que protegerá el mismo tanque
unidireccional, en m
𝑎 celeridad, en m/s

unidireccional se selecciona la pestaña Tubos y elementos entre tramos (figura 3.16) y se
elige el ícono señalado en la figura 4.29, a través del Inspector se proporcionan las
Propiedades generales y las Propiedades físicas de los tanques unidireccionales.
 Nombre, con éste se identifica al tanque unidireccional.
82
 Número de válvulas de retención, en la etapa de predimensionamiento puede

considerarse una sola válvula.
 Número de válvulas de llenado, en la etapa de predimensionamiento puede
considerarse una sola válvula.
Figura 4.29 Pantalla para dar de alta un tanque unidireccional
 Coeficiente de pérdida de carga de una válvula de retención [s2/m5], durante la etapa de
pre-dimensionamiento puede considerarse igual a 0.001.
 Coeficiente de pérdida de carga de una válvula de llenado [s2/m5], durante la etapa de

pre-dimensionamiento puede considerarse igual a 0.1.
 Desplante del tanque [msnm], lo determina la topografía del sitio en donde se construirá.
Se recomienda que los tanques unidireccionales se ubiquen en zonas altas de la
topografía de la conducción.
83
 Área del tanque [m2], el área preliminar del tanque puede calcularse con la ecuación
4.34.
 Tirante del agua en el tanque [m], para que un tanque unidireccional opere durante un
transitorio, la carga piezométrica debe alcanzar un valor menor que la suma de la cota
de desplante del tanque más el tirante de agua, en el sitio en donde se conecta el tanque
con la conducción. A partir de esta condición puede proponerse un valor para el tirante
de agua en el tanque.
 Longitud de la tubería de conexión tanque-conducción [m], (ver figura 4.28) la determina

el sitio en donde se construirá, cuando aún no se ha determinado puede considerarse
de 5 m.
 Diámetro de la tubería de conexión tanque-conducción [m], (ver figura 4.28) en la etapa

de predimensionamiento puede considerarse igual al diámetro de la conducción. En
sistemas de bombeo se recomienda que la velocidad del flujo en las tuberías de
conexión entre los tanques unidireccionales y la línea de conducción sea menor que 4
m/s.
 Bordo libre del agua [m], cuando la distancia entre la superficie del agua y la corona del
tanque es igual al bordo libre del agua se cierran las válvulas de llenado.
En las simulaciones posteriores, los coeficientes de pérdida de carga deberán

corresponder al diseño del arreglo funcional de las válvulas. En la figura 4.30 se presenta
el arreglo funcional de las válvulas de los tanques unidireccionales del acueducto Linares-
Monterrey, México (Guarga et al., 1984).
Al final de la simulación se graba en el archivo de cargas, en mca, la carga en el punto 2

de la figura 4.28; en el de gastos, el gasto que recibe o entrega el tanque en m3/s, positivo
cuando el tanque entrega agua a la conducción y negativo cuando el tanque recibe agua;
en el archivo con sufijo rpm, se graban la elevación del agua en el tanque en msnm y el
volumen de agua acumulado que ha entregado el tanque, en m3, todas las variables en
función del tiempo, en s.
Para que el tanque esté preparado para volver a operar antes de que recupere su volumen
inicial de agua, se recomienda que su área transversal se obtenga con la ecuación 4.35,
3𝑉
𝐴𝑇 = (4.35)
𝑍𝑡𝑜𝑡𝑎𝑙 − 1
donde 𝑉 es volumen de agua entregado por el tanque, en m3 y 𝑍𝑡𝑜𝑡𝑎𝑙 altura total del tanque,
en m, la cual se determina considerando un bordo libre de un metro por debajo de su
84
corona. La altura máxima del tanque depende de las condiciones del suelo y la zona de
construcción.
Tanque unidireccional Línea de

Válvulas de
llenado conducción
Válvulas de
retención
Línea de
Arreglo de fontanería conducción
Figura 4.30 Arreglo funcional de las válvulas de los tanques unidireccionales del
acueducto Linares-Monterrey, México (Guarga et al., 1984)
4.7 Torre de oscilación

Una torre o tanque de oscilación es un tanque abierto en su parte superior que se conecta
directamente a la línea de conducción. Por tanto, para evitar que el agua se derrame por
su coronamiento, la cota de la corona debe ser mayor que la piezométrica que se establece
en flujo establecido en la conexión de la línea que conecta al tanque con la conducción.
Cuando en un acueducto se suspende el bombeo la presión en el interior de la tubería

disminuye y en el sitio en donde se localiza una torre se produce un flujo de la torre hacia
la conducción, debido a la diferencia de presiones que se establece entre la elevación del
agua en la torre y la piezométrica en la conducción, como resultado, se limita la disminución
de presión en la conducción. En presencia de una sobrepresión en la tubería, al ser mayor
la carga en la conducción que en la torre se establecerá un flujo de la conducción hacia la
torre, el ingreso de agua en la torre evita que la presión en la tubería siga aumentando.
Las torres de oscilación dividen a la conducción en dos tramos, uno aguas arriba de la torre
y el otro aguas abajo de ésta. Durante un transitorio por golpe de ariete operan como un
reflector de las ondas de presión, de tal forma que las ondas que se generan en uno de los
85
tramos inciden en la torre y no se transmiten o se transmiten con menor amplitud hacia el

otro tramo, como se muestra en la figura 4.31.
Carga (m)
Figura 4.31 La variación de la presión permanece entre el sitio en donde se produjo el
transitorio y las torres de oscilación
En la figura 4.32 se muestran dos torres de oscilación del acueducto Linares Monterrey,
Nuevo León, México, es común que a las torres localizadas aguas arriba de las plantas de
bombeo de les llame torre de sumergencia.
Figura 4.32 Torres de oscilación, una localizada aguas arriba de la planta de bombeo y la
otra aguas abajo (Grupo de Hidromecánica, IIUNAM)
86
4.7.1 Modelo matemático de torre de oscilación. El modelo matemático está

formado por la ecuación de continuidad en la conexión de la línea principal con la tubería
que va hacia la torre, la ecuación dinámica del agua en la torre y la de la tubería que la
conecta con la línea principal.
El sistema de ecuaciones que modelan el comportamiento de las torres de oscilación es

similar al de los tanques unidireccionales, excepto la ecuación de la tubería que la conecta
con la conducción, debido a que en este caso no hay válvulas. Ésta cambia de la ecuación
4.31 a la 4.36, ver figura 4.33.
Torre de
oscilación
𝑙 𝑄𝑇 < 0
𝑄0
𝑄𝑇 > 0
Punto 1
Punto 2
Figura 4.33 Conexión de torres de oscilación con la conducción
𝑙 𝑄𝑇𝑡+∆𝑡 − 𝑄𝑇𝑡 𝑓𝑙
𝑍𝑡+∆𝑡 − 𝐻1 𝑡+∆𝑡 = + ( )𝑄 |𝑄 | (4.36)
𝑔𝐴 ∆𝑡 2𝑔𝐷𝐴2 𝑇𝑡+∆𝑡 𝑇𝑡+∆𝑡
Entonces, el sistema de ecuaciones en diferencias finitas que modelan a las torres de

oscilación está formado por las ecuaciones 4.4, 4.5, 4.29, 4.30, 4.32 y 4.36, las incógnitas
son 𝑄𝑇𝑡+∆𝑡 , 𝑄1𝑡+∆𝑡 , 𝑄2𝑡+∆𝑡 , 𝐻1𝑡+∆𝑡 , 𝐻2𝑡+∆𝑡 y 𝑍𝑡+∆𝑡 , se tiene un sistema compatible de seis
ecuaciones con seis incógnitas. Por el método de sustitución éste se reduce a la ecuación
4.37 de segundo grado.
𝑈 𝑄𝑇𝑡+∆𝑡 |𝑄𝑇𝑡+∆𝑡 | + 𝑉 𝑄𝑇𝑡+∆𝑡 + 𝑊 = 0 (4.37)
donde
𝑓𝑙
𝑈=
2𝑔𝐷𝐴2
87
𝑙 1 ∆𝑡
𝑉= + +
𝑔𝐴∆𝑡 𝑐𝑎1 + 𝑐𝑎2 𝐴 𝑇
𝑐𝑝 − 𝑐𝑛 𝑙
𝑊= − 𝑄 − 𝑍𝑡
𝑐𝑎1 + 𝑐𝑎2 𝑔𝐴∆𝑡 𝑇𝑡
Cuando 𝑄𝑇 𝑡+∆𝑡 > 0 la torre entrega agua a la línea, cuando 𝑄𝑇 𝑡+∆𝑡 < 0 la torre recibe agua.
Una vez que se ha obtenido 𝑄𝑇 𝑡+∆𝑡 el tirante de agua en la torre (𝑍𝑡+∆𝑡 ) se obtiene con la
ecuación 4.32, las cargas en los puntos 1 y 2 con las ecuación 4.36 y 4.30 y sus gastos
con las ecuaciones características 4.4 y 4.5.
4.7.2 Predimensionamiento de torres de oscilación. La altura de las torres de

oscilación debe ser mayor que la piezométrica que se establece en flujo permanente en el
punto de conexión de las torres con la conducción. Para evitar que el aire ingrese a la
conducción, el diámetro de las torres debe ser tal que durante operaciones en flujo
transitorio no se vacíen.
Una vez elegido el sitio en donde se construirá una torre de oscilación se determina la
elevación de desplante, a través de ésta y la piezométrica en flujo establecido en el punto
de conexión de la torre con la conducción, la primera propuesta para la altura de la torre se
puede calcular con la ecuación 4.38.
𝑍𝑇 = 𝐻1 − 𝑍0 + 20 (4.38)
donde
𝑍𝑇 altura preliminar de la torre de oscilación, en m

𝐻1 piezométrica en flujo estable en la conexión de la torre con la conducción, en m
𝑍0 elevación de desplante de la torre, en msmn
La variación de la superficie libre del agua en la torre depende de su área transversal (𝐴 𝑇 ),

como primera propuesta se sugiere que el diámetro de la torre sea 10 veces el diámetro
de la tubería de la conducción, ecuación 4.39.
𝐷𝑇 = 10 𝐷 (4.39)
donde
𝐷𝑇 diámetro preliminar de la torre, en m
𝐷 diámetro de la conducción principal, en m
88
4.7.3 Simulación en el sistema TransDocencia2-0. Para dar de alta una torre

de oscilación se selecciona la pestaña Tubos y elementos entre tramos (figura 3.16) y se
elige el ícono señalado en la figura 4.34, a través del Inspector se proporcionan las
Propiedades generales y las Propiedades físicas de las torres de oscilación.
Figura 4.34 Pantalla para dar de alta una torre de oscilación
 Nombre, con éste se identifica a la torre de oscilación
 Desplante [msnm], se obtiene de los datos topográficos y la ubicación de la torre. Se
recomienda que las torres de oscilación se ubiquen en zonas altas de la topografía para
que la altura de la torre sea menor que si se construyera en sitios bajos.
 Área de la torre [m2], el área preliminar de una torre puede obtenerse a partir de la
ecuación 4.39.
 Altura [m], la altura preliminar puede calcularse con la ecuación 4.38. Las dimensiones
óptimas de la torre se obtienen a través de analizar los resultados de las simulaciones
realizadas con diferentes áreas y alturas de la torre.
89
 Longitud de la tubería de conexión torre-conducción [m], (ver figura 4.33) la determinan

el trazo del acueducto y la ubicación de la torre.
 Diámetro de la tubería de conexión torre-conducción [m], (ver figura 4.33) en las

primeras simulaciones puede considerarse igual al diámetro de la conducción. En
sistemas de bombeo se recomienda que en las tuberías de conexión entre las torres de
oscilación y la línea de conducción la velocidad del flujo sea máxima de 4 m/s. A partir
del gasto que entrega la torre y la condición mencionada se determina el diámetro de la
tubería de conexión, una vez determinado se debe realizar otra simulación para revisar
el comportamiento de la torre.
 Factor de fricción de la tubería de conexión [D-W], en las primeras simulaciones puede

considerarse igual al de la conducción.
Cuando la longitud, el diámetro o el factor de fricción de la tubería que conecta a la torre

con la conducción es cero, en la simulación del transitorio se considera que no existe esa
tubería.
Al final de la simulación en el archivo de cargas se graba la carga en mca en la conexión

de la línea principal de la conducción con la línea que llega a la torre (punto 1 de la figura
4.33); en el de gastos, el gasto en m3/s que entra o sale de la torre, positivo cuando el agua
sale y negativo cuando la torre recibe agua; en el archivo con sufijo rpm, se graba la
elevación de la superficie libre del agua en la torre en msnm, todas las variables están en
función del tiempo, en s.
4.8 Válvula de admisión y expulsión de aire

A través de las válvulas de admisión y expulsión de aire (VAEA) o válvulas rompedoras de
vacío se crea una conexión entre la atmósfera y el interior de la conducción para evitar
depresiones. Cuando la presión en el interior de la tubería alcanza un valor determinado
por los operadores del sistema se abre la válvula para que ingrese aire, en ese momento
se establece la presión atmosférica en el interior de la tubería.
El flujo de aire que ingresa a la tubería es función del diámetro de la válvula y la presión de
calibración determinada por el operador a la que ésta abrirá, es decir, de la diferencia de
presión entre la atmósfera y el interior de la tubería al momento de abrir la válvula. En la
figura 4.35 se nuestra una VAEA.
90
Figura 4.35 Válvula de admisión y expulsión de aire (VAEA) (Grupo de Hidromecánica.

IIUNAM)
4.8.1 Modelo matemático de válvulas de admisión y expulsión de aire.

En el modelo matemático de la operación de las válvulas de admisión y expulsión de aire
deben considerarse las etapas de admisión y compresión del aire.
Modelo de admisión de aire. Cuando las cargas de presión en las secciones de la

tubería indicadas con los puntos 1 y 2 de la figura 4.36 se igualan al valor de la carga que
comanda la apertura de la válvula y se considera despreciable la pérdida de carga entre
esas dos secciones de la tubería, se cumplen las ecuaciones 4.40 y 4.41. El volumen de
aire que entra a la tubería interfiere con el flujo de agua en los puntos 1 y 2, de tal manera
que el gasto de aire que ingresa a la tubería está relacionado con los gastos de agua en
los puntos 1 y 2 a través de la ecuación 4.42.
𝐻1 = 𝑍 (4.40)
𝐻2 = 𝑍 (4.41)
𝑑𝑉
= 𝑄2 − 𝑄1 (4.42)
𝑑𝑡
donde
𝑍 carga piezométrica que comanda la apertura de la válvula, en msmn
𝑉 volumen de aire que entra a la tubería, en m3
91
Orificio de entrada y salida del aire
Flotador arriba,
Flotador abajo, válvula cerrada
válvula abierta
Q0
Tramo i Tramo i+1
Punto 1 Punto 2
Figura 4.36 Válvula de admisión y expulsión de aire
El sistema formado por las tres ecuaciones 4.40 a 4.42 tiene cinco incógnitas (𝑄1 , 𝑄2 , 𝐻1 ,
𝐻2 y 𝑉), éste se completa con las ecuaciones características 2.41 y 2.42, el cual se
resuelve numéricamente escribiéndolo en diferencias finitas (ecuaciones 4.4, 4.5 y 4.43 a
4.45). En esta nueva nomenclatura las incógnitas son 𝑄1𝑡+∆𝑡 , 𝑄2𝑡+∆𝑡 , 𝐻1𝑡+∆𝑡 , 𝐻2𝑡+∆𝑡 y 𝑉𝑡+∆𝑡 ,
los subíndices 𝑡 y 𝑡 + ∆𝑡 indican el valor de la variable correspondiente en el tiempo 𝑡 y en
el tiempo 𝑡 + ∆𝑡, respectivamente. A través del método de sustitución el sistema se reduce
a la ecuación lineal 4.46.
𝐻1 𝑡+∆𝑡 = 𝑍 (4.43)
𝐻2 𝑡+∆𝑡 = 𝑍 (4.44)
𝑉𝑡+∆𝑡 − 𝑉𝑡 𝑄2𝑡+∆𝑡 + 𝑄2𝑡 𝑄1𝑡+∆𝑡 + 𝑄1𝑡

= − (4.45)
∆𝑡 2 2
∆𝑡
𝑉𝑡+∆𝑡 = 𝑉𝑡 + [𝑄 − 𝑄1𝑡 + (𝑐𝑎1 + 𝑐𝑎2 )𝑍 + 𝑐𝑛 − 𝑐𝑝 ] (4.46)
2 2𝑡
92
Una vez que se determina el volumen 𝑉𝑡+∆𝑡 , se obtiene la diferencia de gastos en 𝑡 + ∆𝑡

con la ecuación 4.45, a través de las ecuaciones 4.4, 4.5, 4.43 y 4.44 se obtienen las cargas
en los puntos 1 y 2, por último con las ecuaciones características se obtienen los gastos.
Modelo de compresión del aire. Existen varias condiciones para que el volumen de
aire que ingresó a la tubería se comprima, entre ellas, cuando el sentido del flujo aguas
abajo de la VAEA se invierte (ahora se dirige del tanque de entrega hacia el bombeo) y su
magnitud es mayor que el flujo en la sección aguas arriba de la VAEA. El proceso de
compresión del aire, escrito en diferencias finitas, se representa con el sistema de
ecuaciones 4.4, 4.5, 4.43, 4.44, 4.47 y 4.48. Las incógnitas son 𝑄1𝑡+∆𝑡 , 𝑄2𝑡+∆𝑡 , 𝐻1𝑡+∆𝑡 , 𝐻2𝑡+∆𝑡 ,
𝑉𝑡+∆𝑡 y 𝑝𝑡+∆𝑡 .
1
𝐻1𝑡+∆𝑡 = (𝑝 − 𝑝𝑎 ) + 𝑍 (4.47)
𝛾 𝑡
𝛼
𝑝𝑡+∆𝑡 𝑉𝑡+∆𝑡 =𝐶 (4.48)
donde
𝑝𝑡 presión absoluta en el interior de la conducción en el tiempo 𝑡, en N/m2
𝑝𝑡+∆𝑡 presión absoluta en el interior de la conducción en el tiempo 𝑡 + ∆𝑡, en N/m2
𝑝𝑎 presión atmosférica, en N/m2
𝛾 peso específico del agua, en N/m3
𝛼 coeficiente politrópico del aire
A través de comparaciones de medidas en prototipo con resultados de simulaciones

numéricas se ha visto que 𝛼 = 1.2 representa adecuadamente el proceso de compresión
del aire que entró a la tubería, la constante 𝐶 de la ecuación 4.48 se determina con la
condición de operación inmediata anterior al inicio de la compresión.
El sistema anterior de ecuaciones se reduce a la ecuación 4.49 cuya incógnita (𝐻1𝑡+∆𝑡 ) es

la carga en la sección de la tubería en donde se localiza la VAEA, ésta puede resolverse
con el método iterativo de Newton-Raphson. Una vez calculada 𝐻1𝑡+∆𝑡 se obtiene 𝐻2𝑡+∆𝑡 con
las ecuaciones 4.43 y 4.44. Con las ecuaciones características 4.4 y 4.5 se determinan
𝑄1𝑡+∆𝑡 y 𝑄2𝑡+∆𝑡 , con la ecuación 4.45 se calcula 𝑉𝑡+∆𝑡 y finalmente con la ecuación 4.48 se
determina 𝑝𝑡+∆𝑡 .
𝛼
∆𝑡
(𝐻1𝑡+∆𝑡 + 𝐻1𝑡 − 𝑍) [𝑉𝑡 + (𝑄2𝑡 − 𝑄1𝑡 + 𝑐𝑛 − 𝑐𝑝 + (𝑐𝑎1 + 𝑐𝑎2 )𝐻1𝑡 )] − 𝐶 = 0 (4.49)
2
93
4.8.2 Simulación en el sistema TransDocencia2-0. Para dar de alta una

válvula de admisión y expulsión de aire se selecciona la pestaña Tubos y elementos entre
tramos (figura 3.16) y se elige el ícono señalado en la figura 4.37.
 Nombre, con éste se identifica a la válvula de admisión y expulsión de aire.
Figura 4.37 Pantalla para dar de alta una válvula de admisión y expulsión de aire
Para simular la operación de una VAEA en el sistema TransDocencia2-0 no se necesita
proporcionar sus dimensiones, en este caso, los resultados que se obtienen después de
simular una VAEA se utilizan para determinar su diámetro. La información que debe
proporcionarse es la siguiente.
94
 Cota de la entrada de aire [msnm], valor de la carga piezométrica a partir del cual inicia
el ingreso de aire a la conducción.
 Presión atmosférica [mca], presión atmosférica en el sitio en donde se localiza la VAEA.
 Coeficiente politrópico, se ha visto que la relación entre el volumen del aire y la presión
está bien representada a través de la ley de compresión adiabática de un gas perfecto
(ecuación 4.48) con un coeficiente politrópico igual a 1.2, sin embargo, éste puede
cambiarse.
Durante la simulación, se considera que la masa de aire que entró a la conducción

permanece en el sitio en donde está la válvula, pero su volumen varía debido a las
variaciones de presión en el interior de la tubería. Como consecuencia del viaje de la onda
de presión pueden presentarse depresiones que ocasionen que la válvula vuelva a abrir e
ingrese más aire aumentando la masa y el volumen de aire.
Al final de la simulación en el archivo de cargas se graba la carga en la conexión de la

válvula con la conducción en mca, en el de gastos el flujo de aire que ingresa a la
conducción, en m3/s y en el archivo con sufijo rpm el volumen de aire, en m3; todas las
variables en función del tiempo, en s.
4.8.3 Dimensionamiento de válvulas de admisión y expulsión de aire. El

dimensionamiento de una VAEA consiste en determina el diámetro de la válvula que
permita el flujo de aire necesario para evitar depresiones en el interior de la tubería. El flujo
de aire que ingresa a la tubería es función del diámetro de la válvula y de la diferencia de
presión entre la atmósfera y el interior de la tubería.
La presión en el interior de la tubería es la diferencia de la Cota de la entrada de aire [mca]

(figura 4.37) menos la elevación del sitio en donde se instala la VAEA (dato que se obtiene
de la topografía del terreno). Por ejemplo, si una VAEA se instala en la elevación 250 msnm
y se decide que la válvula abra cuando la carga piezométrica sea de 253 mca en el sitio de
conexión de la VAEA, la válvula abrirá cuando la presión en el interior de la tubería sea de
3 mca (presión diferencial de 3 mca).
El diámetro de la VAEA se determina a través de los nomogramas que proporcionan los

fabricantes de válvulas y de la gráfica del flujo de aire que ingresa a la conducción que se
obtiene como resultado de simulaciones de transitorios.
Como ejemplo, supongamos que para una presión diferencial de 3 mca (4.26 PSI) el flujo
máximo de aire que ingresó a la tubería es de 2.35 m3/s (83 ft3/s) como se muestra la figura
4.38. A partir del nomograma de la figura 4.39 se observa que el punto (83, 4.26) está entre
las curvas correspondientes a las válvulas con diámetros 6” y 8”. Para garantizar el flujo de
83 ft3/s debe elegirse la válvula de 8” de diámetro.
95
2.5
Gasto de aire que entra a la tubería, en

2.0
1.5
m3/s
1.0
0.5
0.0
0 50 100 150 200
Tiempo, en s
Figura 4.38 El flujo máximo que ingresa a la tubería es de 2.35 m3/s
Figura 4.39 Un flujo de 83 ft3/s a una presión diferencial de 4.26 PSI se obtiene con una
válvula de 8” de diámetro
En los casos en que las válvulas resulten ser demasiado grandes y por tanto de costo muy
elevado, está la alternativa de instalar dos o más válvulas de menor diámetro que
proporcionen entre todas ellas el gasto de aire necesario que debe ingresar a la tubería.
96
ELEMENTOS DE UN SISTEMA DE BOMBEO QUE GENERAN TRANSITORIOS
5
Elementos de sistemas de
bombeo que generan
transitorios
5.1 Introducción
Un flujo transitorio se desarrolla en el intervalo de tiempo en el que un sistema pasa de una
condición de operación en flujo estacionario a otra condición de operación de flujo
estacionario. En sistemas de bombeo los cambios en la operación se producen de manera
programada o súbita, los primeros se realizan para incrementar o disminuir la demanda del
gasto de agua y para dar mantenimiento al sistema, los segundos se presentan cuando las
bombas dejan de operar súbitamente como consecuencia de fallas en el suministro de
energía eléctrica. Los elementos de los sistemas de bombeo que tienen la capacidad de
cambiar la forma de operación son las válvulas y las bombas, en este capítulo se presentan
los modelos matemáticos que los representan y cómo se simulan en el sistema
TransDocencia2-0.
5.2 Válvulas de movimiento programado o controlado

La operación de las válvulas de movimiento programado o controlado no depende del flujo.
El estado de apertura es independiente del flujo, su cierre y apertura se operan
manualmente o por medio de un motor, el tiempo en que se realizan las maniobras es
programable.
97
La pérdida de carga que producen las válvulas depende de su posición o porcentaje de

apertura con respecto al eje de la tubería, como se muestra en la figura 5.1, se usan para
controlar el flujo.
abierta semi abierta cerrada

(a) (b) (c)
Figura 5.1 Válvula esférica, el flujo que la atraviesa depende del grado de apertura
En sistemas de bombeo, se instalan en puntos de la conducción en donde se necesite

controlar el flujo y en la descarga de algunas bombas de eje horizontal, figura 5.2.
98
Figura 5.2 Válvula esférica, en la descarga de una bomba de eje horizontal (Grupo de
Hidromecánica, IIUNAM)
En conducciones ascendentes cuando las bombas dejan de operar, limitan el caudal del
agua que retorna hacia las bombas con la finalidad de evitar que giren en sentido inverso
a altas velocidades.
5.2.1 Modelo matemático de válvulas de movimiento controlado. El

modelo de la operación de la válvula se plantea para dos casos
 Cuando la válvula está completamente cerrada.
 Cuando la válvula está parcialmente abierta o en movimiento.
Válvula completamente cerrada. El gasto a través de la válvula y en las secciones de

la tubería inmediatas aguas arriba y aguas abajo de la válvula (secciones indicadas como
punto 1 y punto 2 de la figura 5.3) es cero, ecuación 5.1.
𝑄1𝑡+∆𝑡 = 𝑄2𝑡+∆𝑡 = 0 (5.1)
99
Las cargas en las secciones de la tubería inmediatas aguas arriba y aguas abajo de la
válvula se obtienen a través de las ecuaciones 5.2 y 5.3, las que a su vez se obtienen de
las ecuaciones características 4.4 y 4.5.
𝑐𝑝
𝐻1𝑡+∆𝑡 = (5.2)
𝑐𝑎𝑖
𝑐𝑛
𝐻2𝑡+∆𝑡 = − (5.3)
𝑐𝑎𝑖+1
Tramo 𝑖 Tramo 𝑖 + 1
𝜃
𝑄0
Punto 1 Punto 2
Figura 5.3 Secciones aguas arriba y abajo de la válvula
Válvula parcialmente abierta o en movimiento. La válvula introduce una diferencia

de carga en las secciones inmediatas aguas arriba y aguas abajo de la válvula que depende
de su grado de apertura, la relación entre las cargas en el tiempo 𝑡 + ∆𝑡 está dada por la
ecuación 5.4. Cuando no cambia la apertura de la válvula el coeficiente de pérdida es
constante, no es función del tiempo.
𝐻1𝑡+∆𝑡 − 𝐻2𝑡+∆𝑡 = 𝐾𝑉𝑡+∆𝑡 𝑄𝑉𝑡+∆𝑡 |𝑄𝑉𝑡+∆𝑡 | (5.4)
Los gastos en esas secciones son iguales entre sí e iguales al gasto que pasa a través de
la válvula, ecuaciones 5.5 y 5.6
𝑄1𝑡+∆𝑡 = 𝑄2𝑡+∆𝑡 (5.5)
𝑄1𝑡+∆𝑡 = 𝑄𝑉𝑡+∆𝑡 (5.6)
donde
𝐻𝑖𝑡+∆𝑡 carga en el punto 𝑖 en el tiempo 𝑡 + ∆𝑡, en m

𝑄𝑖𝑡+∆𝑡 gasto en el punto 𝑖 en el tiempo 𝑡 + ∆𝑡, en m3/s
𝑄𝑉𝑡+∆𝑡 gasto a través de la válvula en el tiempo 𝑡 + ∆𝑡, en m3/s
𝐾𝑉𝑡+∆𝑡 coeficiente de pérdida de carga que introduce la válvula en el tiempo 𝑡 + ∆𝑡, en
s2/m5
100
El sistema determinado por las tres ecuaciones 5.4 a 5.6 contiene a las cinco incógnitas
𝑄𝑉𝑡+∆𝑡 , 𝑄1𝑡+∆𝑡 , 𝑄2𝑡+∆𝑡 , 𝐻1𝑡+∆𝑡 , 𝐻2𝑡+∆𝑡 , éste se completa con las ecuaciones características 4.4
y 4.5. A través del método de sustitución el sistema de reduce a la ecuación de segundo
grado 5.7.
1 1 𝑐𝑛 𝑐𝑝
𝐾𝑉𝑡+∆𝑡 |𝑄𝑉𝑡+∆𝑡 |𝑄𝑉𝑡+∆𝑡 + ( + ) 𝑄𝑉𝑡+∆𝑡 − − =0 (5.7)
𝑐𝑎i 𝑐𝑎i+1 𝑐𝑎i+1 𝑐𝑎i
donde
𝑐𝑎𝑗 constante del método de las características del tramo 𝑗, en m2/s
Cuando la válvula está en movimiento, para determinar el gasto en la ecuación 5.7 primero
debe conocerse el coeficiente de pérdida de carga 𝐾𝑉𝑡+∆𝑡 , el cual depende del tiempo
debido a que la velocidad de cierre o apertura la determina el operador del sistema. Este
coeficiente se obtiene a partir de la información que proporcionan los fabricantes de
válvulas a través de gráficas como la de la figura 5.4. En el sistema TransDocencia2-0 se
considera el movimiento de las válvulas con velocidad constante.
Figura 5.4 Coeficiente de pérdida de carga en función del ángulo de apertura
101
5.2.2 Simulación en el sistema TransDocencia2-0. En el sistema

TransDocencia2-0, a las válvulas de movimiento programado o controlado se le llama
válvulas de control, están presentes en los siguientes arreglos.
 En la descarga de las bombas, ya sea que la planta de bombeo determine la frontera

aguas arriba del sistema o que esté a lo largo de la conducción, ambos casos se
presentan en la sección 5.3.
 La válvula se localiza entre dos tramos de la conducción, pueden simularse el cierre

total o la apertura total. En ambos casos la apertura en flujo estable puede ser parcial.
 La válvula se localiza en la frontera aguas abajo del sistema. Únicamente puede

simularse el cierre total de la válvula, la apertura en flujo estable puede ser parcial.
A continuación, se indica la información necesaria para simular la operación de una válvula

de control localizada entre dos tramos de la conducción.
En la figura 5.5 se muestra la pantalla para dar de alta una válvula que se localiza entre
dos tramos de la conducción. Esta pantalla se presenta seleccionando primero la pestaña
Maniobras entre tamos (figura 3.16) y después el ícono Válvula de control, la información
necesaria para simular su operación se pide en los siguientes campos de captura.
Figura 5.5 Información para simular una válvula de movimiento programado
102
 Nombre, con éste se identifica a la válvula.
 Maniobra, se indica si se trata de una apertura o cierre.
 Tipo de válvula de control, se elige un tipo de válvula a través de la lista que se

despliega, figura 5.6. La información detallada de las válvulas puede verse
seleccionando la pestaña Catálogos y después Curvas de válvulas (figura 3.5).
Figura 5.6 Se despliega una lista para seleccionar el tipo de válvula de control
 Tiempo en que se realiza la maniobra [s], es el tiempo en el que la apertura de la válvula

pasa de su posición en flujo estable a su posición final.
 Instante en el que inicia la maniobra [s], el movimiento de la válvula puede iniciar en

cualquier instante de la simulación. Para simular que la válvula no se abre ni se cierra,
basta con indicar que este instante es mayor que el tiempo de simulación.
103
 Posición de la válvula en flujo estable, durante la operación en flujo estable la válvula

puede estar parcialmente abierta.
El coeficiente de pérdida de carga que introducen las válvulas depende del estado de
apertura de la válvula. Algunos fabricantes presentan el coeficiente de pérdida en
función de la posición en grados del elemento que obstruye el paso del agua como es
el caso de las válvulas de mariposa y esférica, figura 5.4, otros lo presentan en función
del porcentaje de apertura, para ambos casos, la posición igual a cero significa que la
válvula está completamente cerrada. Cuando la válvula está completamente abierta el
dato de la posición depende del tipo de válvula, 90° o 100% (el catálogo de válvulas
únicamente contiene información de válvulas de mariposa y esférica). En la figura 5.4
se presenta la correspondiente a la válvula de mariposa, para la que:
𝑝𝑜𝑠𝑖𝑐𝑖ó𝑛 = 0 𝑠𝑖𝑔𝑛𝑖𝑓𝑖𝑐𝑎 𝑣á𝑙𝑣𝑢𝑙𝑎 𝑐𝑜𝑚𝑝𝑙𝑒𝑡𝑎𝑚𝑒𝑛𝑡𝑒 𝑐𝑒𝑟𝑟𝑎𝑑𝑎
𝑝𝑜𝑠𝑖𝑐𝑖ó𝑛 = 90 𝑠𝑖𝑔𝑛𝑖𝑓𝑖𝑐𝑎 𝑣á𝑙𝑣𝑢𝑙𝑎 𝑐𝑜𝑚𝑝𝑙𝑒𝑡𝑎𝑚𝑒𝑛𝑡𝑒 𝑎𝑏𝑖𝑒𝑟𝑡𝑎
 Diámetro de la válvula [m], es el diámetro de la válvula.
Cuando la válvula está en la descarga de una bomba o en la frontera aguas abajo

únicamente puede simularse el cierre, por lo que no existe el campo Maniobra.
En el sistema TransDocencia2-0, el coeficiente de pérdida de la válvula se obtiene a través

de la ecuación 5.8
𝐾𝑉 = 𝐾𝑚í𝑛 10𝑃( 𝜏) (5.8)
donde
𝐾𝑉 coeficiente de pérdida de carga que introduce la válvula, en s2/m5
𝐾𝑚í𝑛 coeficiente de pérdida de carga que introduce la válvula cuando está totalmente
abierta, en s2/m5
𝜏 variable adimensional
𝑃(𝜏) polinomio de sexto grado
Cuando el coeficiente de pérdida de la válvula está en función de la posición (𝜃, en

grados) del elemento que obstruye el paso del agua 𝜏 = 𝜃⁄90, cuando está en función del
porcentaje (𝑦) de apertura 𝜏 = 𝑦⁄100. De tal forma que cuando la válvula está
completamente abierta 𝜏 = 1 (𝑘𝑉 = 𝑘𝑚í𝑛 ) y cuando está totalmente cerrada 𝜏 = 0 (𝑘𝑉 =
𝑘𝑚á𝑥 ), ver figura 5.4.
El polinomio de sexto grado 𝑃(𝜏) = 𝑎0 + 𝑎1 𝜏 + 𝑎2 𝜏 2 + 𝑎3 𝜏 3 + 𝑎4 𝜏 4 + 𝑎5 𝜏 5 + 𝑎6 𝜏 6 se

obtiene a partir de un ajuste por mínimos cuadrados de la curva determinada por los puntos
(𝜏, log(𝑘⁄𝑘𝑚í𝑛 )), éste se presenta en la parte inferior de la figura 5.4. Hay dos condiciones
104
necesarias para que un polinomio pueda representar el movimiento de una válvula de

control.
1. La suma de sus coeficientes debe ser igual a cero, la tolerancia que permite el
sistema TransDocencia2-0 está dada por la ecuación 5.9
|𝑎0 + 𝑎1 + 𝑎2 + 𝑎3 + 𝑎4 + 𝑎5 + 𝑎6 | ≤ 0.05 (5.9)
2. El término independiente debe satisfacer la condición dada por la ecuación 5.10
𝑎0 > 1 (5.10)
Al final de la simulación, en el archivo de cargas se graba la carga inmediatamente aguas

arriba y aguas abajo de la válvula, en mca y en el de gastos, el flujo a través de la válvula
en m3/s; todas las variables en función del tiempo, en s.
En el archivo que resume los datos del sistema de bombeo (*_IMP.txt) la posición de la
válvula en flujo estable siempre se reporta en porcentaje de apertura.
5.3 Disparo de equipos de bombeo

La función de las bombas hidráulicas es transformar energía mecánica en energía
hidráulica. La energía mecánica proviene de otro tipo de energía que, en los sistemas de
conducción de agua es muy común que sea energía eléctrica, razón por la que en este
trabajo se hará referencia a la energía eléctrica.
Mientras exista suministro de energía eléctrica, la velocidad de giro de las bombas es

constante y el flujo del fluido de trabajo (en este caso agua) es estable o permanente. Los
fabricantes de bombas representan la operación de las mismas a través de curvas carga
vs gasto, potencia vs gasto, eficiencia vs gasto y el nivel mínimo de succión requerido
(NPSH, siglas en inglés) para evitar cavitación, también en función del gasto. Las unidades
de las variables dependen del lugar de fabricación de las bombas, es común que el gasto
se presente en galones por minuto (gpm), litros por segundo (l/s) o metros cúbicos por
segundo (m3/s), la carga y el NPSH en pies (ft) o metros (m) y la potencia en caballos de
fuerza (hp) o kilo Watts (kW), como se muestra en las figuras 5.7 y 5.8.
Las gráficas de la figura 5.7 presentan la información para un único diámetro de impulsor,
a diferencia de la figura 5.8 en la que se presenta para tres diámetros del impulsor (12”,
14” y 16”); en ambos casos, el tipo de bomba se identifica a través de su velocidad
específica que se calcula con la ecuación 5.11. Observa que la velocidad específica es la
misma para los tres impulsores de la figura 5.8.
105
√𝑄𝑅
𝑛
𝑁𝑠 = 𝑁𝑅 (5.11)
𝐻 0.75
( 𝑅)
𝑚
donde
𝑁𝑆 velocidad específica de la bomba, en rpm (m3/s)0.5 m-0.75 o rpm (gpm)0.5 ft-0.75
𝑁𝑅 velocidad de rotación de diseño o nominal de la bomba, en rpm
𝑄𝑅 gasto que entrega la bomba en operación a máxima eficiencia, en m3/s o gpm
𝐻𝑅 carga que entrega la bomba en operación a máxima eficiencia, en m o ft
𝑛=1 cuando la bomba es de una succión
𝑛=2 cuando la bomba es de doble succión
𝑚 número de impulsores
Cada una de las curvas de las figuras 5.7 y 5.8, pueden representarse a través de
polinomios. Es muy común que con polinomios de tercer grado se obtenga una buena
aproximación.
180
170
160
150
Carga (m)
140
130
120
110
100
6000
Potencia (kW)
5000
4000
3000
90
80
Eficiencia (%)
70
60
50
40
30
20
10
0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5
Gasto (m 3/s)
Figura 5.7 Curvas características de operación en flujo permanente de una bomba de eje
horizontal, doble succión, de velocidad específica 𝑁𝑆 = 46 (rpm (m3/s)0.5 m-0.75)
106
300
250
200
Carga (pies)
150
100
50
0 500 1000 1500 2000 2500 3000 3500 4000
Gasto (gpm)
Figura 5.8 Curvas características de operación en flujo permanente de una bomba con
velocidad específica 𝑁𝑆 = 93.6 (rpm (m3/s)0.5 m-0.75)
Cuando se suspende el suministro de energía eléctrica la velocidad de la bomba no puede

mantenerse constante, ésta variará en función del flujo de agua que pasa a través de ella
y de los momentos de inercia de la misma bomba y del motor eléctrico. La operación de la
bomba pasa de ser permanente a transitoria, esta última forma de operación se representa
a través de gráficas como las de las figuras 5.9 y 5.10 a las que se les conoce como curvas
de cuatro cuadrantes.
Observa que cuando la válvula en la descarga de la bomba es de no retorno, tanto el flujo

como la velocidad de la bomba no pueden ser negativos, por lo que la operación de la
bomba en flujo transitorio se desarrolla en el primer cuadrante de las gráficas de las figuras
5.9 y 5.10.
107
Velocidad, n
100
-200 -100 100 200 Gasto, q
-100
Figura 5.9 Curva de carga para la operación en flujo transitorio de una bomba de
velocidad específica Ns = 49 (rpm m1.5/s0.5 m-0.75)
108
Velocidad, n
220 % 
180 % 
140 % 
100 % 
100
60 % 
20 % 
-200 -100 100 200 Gasto, q
-100
-200
Figura 5.10 Curva de momento o par para la operación en flujo transitorio de una
bomba de velocidad específica Ns = 49 (rpm m1.5/s0.5 m-0.75)
Las curvas de cuatro cuadrantes relacionan a las variables adimensionales de carga (ℎ),
gasto (𝑞), velocidad (𝑛) y par o momento (𝛽) definidas a través de las ecuaciones 5.12 a
5.15.
109
𝐻
ℎ= (5.12)
𝐻𝑅
𝑄
𝑞= (5.13)
𝑄𝑅
𝑁
𝑛= (5.14)
𝑁𝑅
𝑇
𝛽= (5.15)
𝑇𝑅
donde
𝑁𝑅 velocidad de diseño de la bomba, en rpm
𝑇𝑅 torque, momento angular o par aplicado al eje de la bomba en operación con máxima
eficiencia, en Nm o lb-ft
𝐻 carga que entrega la bomba, en m o ft
𝑄 gasto en la bomba, en m3/s o gpm
𝑁 velocidad de la bomba, en rpm
𝑇 torque, momento angular o par aplicado al eje de la bomba, en Nm o lb-ft
ℎ carga, adimensional
𝑞 gasto, adimensional
𝑛 velocidad de la bomba, adimensional
𝛽 torque, momento angular o par, adimensional
La información que proporcionan todas las curvas de las figuras 5.9 y 5.10 puede
representarse a través de las funciones definidas en las ecuaciones 5.16 y 5.17,
respectivamente, las que dependen del ángulo definido por la ecuación 5.18, en Carmona
y Carmona (2013) se presenta en procedimiento para obtener las funciones. A la función
𝑓ℎ se le conocen como función característica de carga y a 𝑓𝛽 como función característica
de momento, par o torque.
ℎ
𝑓ℎ (𝜃) = (5.16)
𝑛2 + 𝑞2
𝛽
𝑓𝛽 (𝜃) = (5.17)
𝑛2 + 𝑞2
𝑛
𝜃 = tan−1 ( ) (5.18)
𝑞
donde
𝑓ℎ función característica de carga, adimensional
𝑓𝛽 función característica de par o momento, adimensional
𝜃 ángulo en grados
110
En el sistema TransDocencia2-0, la operación de las bombas en flujo transitorio se

representa con las funciones características de carga (𝑓ℎ ) y momento o par (𝑓𝛽 ) definidas
por las ecuaciones 5.16 y 5.17, dadas a través de tablas de datos como la Tabla 5.1
Tabla 5.1 Funciones características de carga y momento (Grupo de Hidromecánica,

IIUNAM)
𝑁𝑆 = 38 𝑁𝑆 = 44 𝑁𝑆 = 46
𝜃° 𝑓ℎ 𝑓𝛽 𝑓ℎ 𝑓𝛽 𝑓ℎ 𝑓𝛽
0 -0.530 -0.350 -0.574 -1.045 -0.530 -0.350
5 -0.461 -0.265 -0.550 -0.825 -0.476 -0.265
10 -0.392 -0.180 -0.500 -0.583 -0.392 -0.180
15 -0.271 -0.072 -0.403 -0.354 -0.291 -0.062
20 -0.150 0.037 -0.276 -0.134 -0.150 0.037
25 -0.038 0.133 -0.124 0.045 -0.037 0.135
30 0.075 0.228 0.023 0.176 0.075 0.228
35 0.199 0.309 0.160 0.309 0.200 0.320
40 0.322 0.390 0.316 0.405 0.320 0.410
45 0.500 0.500 0.505 0.501 0.500 0.500
50 0.678 0.579 0.669 0.577 0.650 0.555
55 0.834 0.619 0.839 0.640 0.810 0.604
60 0.989 0.658 0.968 0.680 0.930 0.625
65 1.112 0.657 1.084 0.700 1.050 0.640
70 1.234 0.656 1.220 0.713 1.170 0.643
75 1.307 0.611 1.272 0.735 1.250 0.645
80 1.380 0.566 1.373 0.750 1.330 0.640
85 1.427 0.518 1.408 0.767 1.370 0.630
90 1.473 0.470 1.472 0.778 1.410 0.620
95 1.472 0.501 1.502 0.786 1.420 0.635
100 1.470 0.532 1.513 0.824 1.440 0.690
105 1.464 0.574 1.462 0.835 1.445 0.760
110 1.457 0.615 1.427 0.861 1.450 0.820
115 1.429 0.718 1.425 0.876 1.440 0.890
120 1.400 0.820 1.418 0.923 1.420 0.960
125 1.376 0.883 1.411 1.011 1.400 1.020
130 1.351 0.946 1.380 1.099 1.350 1.080
135 1.301 0.997 1.366 1.200 1.310 1.143
140 1.251 1.047 1.318 1.300 1.245 1.180
145 1.178 1.069 1.310 1.400 1.180 1.178
150 1.105 1.090 1.257 1.420 1.100 1.165
155 0.990 1.069 1.215 1.407 1.030 1.140
160 0.875 1.048 1.154 1.427 0.940 1.110
𝑁𝑆 en rpm m1.5/s0.5 m-0.75 (continua)
111
Tabla 5.1 Funciones características de carga y momento (continuación)
𝑁𝑆 = 38 𝑁𝑆 = 44 𝑁𝑆 = 46
𝜃° 𝑓ℎ 𝑓𝛽 𝑓ℎ 𝑓𝛽 𝑓ℎ 𝑓𝛽
165 0.780 0.987 1.097 1.403 0.850 1.070
170 0.685 0.925 1.065 1.340 0.750 1.030
175 0.664 0.904 1.033 1.288 0.700 0.990
180 0.642 0.882 0.972 1.262 0.667 0.950
185 0.620 0.861 0.919 1.216 0.665 0.910
190 0.598 0.840 0.849 1.157 0.660 0.860
195 0.587 0.794 0.778 1.067 0.650 0.810
200 0.576 0.748 0.709 0.967 0.630 0.750
205 0.551 0.684 0.624 0.830 0.600 0.670
210 0.525 0.620 0.555 0.687 0.550 0.600
215 0.503 0.535 0.485 0.547 0.510 0.500
220 0.481 0.450 0.454 0.382 0.480 0.400
225 0.487 0.344 0.434 0.287 0.460 0.300
230 0.493 0.238 0.431 0.168 0.470 0.190
235 0.529 0.134 0.422 0.062 0.490 0.100
240 0.565 0.030 0.392 -0.104 0.530 0.010
245 0.607 -0.084 0.354 -0.292 0.570 -0.090
250 0.649 -0.198 0.313 -0.476 0.620 -0.200
255 0.676 -0.283 0.271 -0.643 0.660 -0.320
260 0.702 -0.368 0.239 -0.863 0.690 -0.450
265 0.711 -0.506 0.199 -1.030 0.715 -0.575
270 0.719 -0.644 0.192 -1.208 0.740 -0.700
𝑁𝑆 en rpm m1.5/s0.5 m-0.75
En la figura 5.11 se muestran las curvas características de carga y momento

correspondientes a la bomba de velocidad específica 𝑁𝑆 = 44 (rpm m1.5/s0.5 m-0.75).
112
Curvas características 𝑓ℎ , 𝑓𝛽
1
-1
-2
0 90 180 270 360
carga momento
Ángulo 𝜃°
Figura 5.11 Curvas características de carga y momento, 𝑁𝑆 = 44 (rpm m1.5/s0.5 m-0.75)
5.3.1 Modelo matemático. El disparo, paro súbito o paro de emergencia de equipos

de bombeo ocurre cuando se suspende el suministro de energía a las bombas.
Cuando las bombas se desconectan de la energía eléctrica su velocidad depende del flujo
que las atraviesa y de los momentos de inercia del motor eléctrico y de la misma bomba.
La velocidad de la bomba se determina a través de la ecuación dinámica 5.19 también
conocida como de las masas rodantes.
𝑑𝜔
𝑇 = −𝑊𝑅2 (5.19)
𝑑𝑡
donde
𝑊𝑅2 momento angular de inercia total (bomba, motor y agua que pasa por la bomba,
determinado con el radio del impulsor de la bomba y masa), en kg m2
𝜔 velocidad angular de la bomba, en rad/s
Al escribir la velocidad angular de la bomba en revoluciones por minuto, la ecuación 5.19

se escribe como la ecuación 5.20.
2𝜋 𝑑𝑁
𝑇 = −𝑊𝑅2 (5.20)
60 𝑑𝑡
Al sustituir las definiciones de la velocidad y el torque adimensionales (ecuaciones 5.14 y

5.15), escribir la derivada como el cociente del cambio de la velocidad adimensional entre
113
el intervalo de tiempo en el que ocurre el cambio y se considera el promedio del torque

adimensional en ese intervalo de tiempo, la ecuación 5.20 toma la forma de la ecuación
5.21.
𝑛𝑡+∆𝑡 − 𝑛𝑡 60 𝑇𝑅 𝛽𝑡+∆𝑡 + 𝛽𝑡
= (5.21)
∆𝑡 2𝜋 𝑊𝑅2 𝑁𝑅 2
donde
𝑛𝑖 velocidad de la bomba en el tiempo 𝑖, adimensional
𝛽𝑖 torque en el tiempo 𝑖, adimensional
El momento de inercia del conjunto motor-bomba el sistema Trans lo determina con la

ecuación 5.22 tomando en consideración la ecuación 5.23 (Thorley, 1991).
𝑃 0.9556 𝑃 1.48
𝑊𝑅 2 = 0.03768 ( ) + 0.0043 ( ) (5.22)
𝑁3 𝑁
𝜌𝑔𝑄𝑅 𝐻𝑅
𝑃= (5.23)
𝜂
donde
𝑃 potencia mecánica en la operación con máxima eficiencia, en kW
𝑁 velocidad de diseño de la bomba, en miles de rpm
𝜂 eficiencia máxima
𝜌 densidad del fluido, en kg/m3
𝑔 aceleración de la gravedad, en m/s2
Como se mencionó anteriormente, ante la falta de energía eléctrica, la velocidad de la

bomba variará en función del flujo de agua que pasa a través de ella y de los momentos de
inercia de la misma bomba y del motor eléctrico, como resultado la carga que entrega
también estará cambiando. La relación que existe entre la carga (ℎ), el gasto (𝑞), la
velocidad (𝑛) y el par (𝛽), está dada por las funciones características de carga y momento
(ecuaciones 5.16 a 5.18) cuyos valores se obtienen a través de los datos de la Tabla 5.1.
Cuando la velocidad de la bomba y el gasto que pasa por ella determinan un ángulo 𝜃
(ecuación 5.18) cuyo valor no está en la Tabla 5.1, las funciones características de carga
(𝑓ℎ ) y momento (𝑓𝛽 ) se obtienen a través de las interpolaciones lineales dadas por las
ecuaciones 5.24 y 5.25 respectivamente.
𝑛
𝑓ℎ (𝜃) = 𝑎1 + 𝑎2 𝑡𝑎𝑛−1 ( ) (5.24)
𝑞
𝑛
𝑓𝛽 (θ) = 𝑎3 + 𝑎4 𝑡𝑎𝑛−1 ( ) (5.25)
𝑞
Las constantes 𝑎1 , 𝑎2 , 𝑎3 y 𝑎4 se obtienen con las ecuaciones 5.26 a 5.29.
114
𝑓ℎ (𝜃2 ) − 𝑓ℎ (𝜃1 )
𝑎1 = 𝑓ℎ (𝜃1 ) − ( ) 𝜃1 (5.26)
𝜃2 − 𝜃1
𝑓ℎ (𝜃2 ) − 𝑓ℎ (𝜃1 )
𝑎2 = (5.27)
𝜃2 − 𝜃1
𝑓𝛽 (𝜃2 ) − 𝑓𝛽 (𝜃1 )
𝑎3 = 𝑓𝛽 (𝜃1 ) − ( ) 𝜃1 (5.28)
𝜃2 − 𝜃1
𝑓𝛽 (𝜃2 ) − 𝑓𝛽 (𝜃1 )
𝑎4 = (5.29)
𝜃2 − 𝜃1
Donde los ángulos 𝜃1 y 𝜃2 para los que sí hay información, determinan el intervalo que
contiene al ángulo 𝜃, para el que se está buscando el valor de las funciones características.
En el tiempo 𝑡 + ∆𝑡, la carga aguas abajo de la válvula de descarga de la bomba, punto 2

de la figura 5.12, se determina a través de la ecuación 5.30 y por continuidad se satisface
la ecuación 5.31. El gasto que entrega la bomba está dado por la ecuación 5.32.
Para cada sección, aguas arriba y abajo del conjunto bomba-válvula de descarga, la
relación entre la carga y el gasto está dada por las ecuaciones características de los puntos
1 y 2 de la figura 5.12, ecuaciones 4.4 y 4.5.
Válvula de descarga
Tramo 2
𝐻𝑆 Tramo 1
Punto 1 Punto 2
Figura 5.12 Bomba entre dos tramos de la conducción (con línea de succión)
𝐻2𝑡+∆𝑡 = 𝐻1𝑡+∆𝑡 + 𝐻𝑅 ℎ𝑡+∆𝑡 − 𝐾𝑉𝑡+∆𝑡 𝑄𝑅2 |𝑞𝑡+∆𝑡 |𝑞𝑡+∆𝑡 (5.30)
𝑄1𝑡+∆𝑡 = 𝑄2𝑡+∆𝑡 (5.31)
115
𝑄2𝑡+∆𝑡 = 𝑄𝑅 𝑞𝑡+∆𝑡 (5.32)

donde
𝐾𝑉𝑡+∆𝑡 coeficiente de pérdida de carga que introduce la válvula de descarga, en s2/m5
El sistema de ecuaciones que determinan el modelo de la operación en flujo transitorio de

una bomba está formado por las ecuaciones 5.21, 5.24, 5.25, 5.30 a 5.32, 4.4 y 4.5, es un
sistema de ocho incógnitas con las variables 𝑄1𝑡+∆𝑡 , 𝑄2𝑡+∆𝑡 , 𝐻1𝑡+∆𝑡 , 𝐻2𝑡+∆𝑡 , 𝑞𝑡+∆𝑡 , ℎ𝑡+∆𝑡 , 𝑛𝑡+∆𝑡 y
𝛽𝑡+∆𝑡 .
Para válvulas de cierre controlado el coeficiente 𝐾𝑉𝑡+∆𝑡 se obtiene a través de la ecuación

5.8. Para válvulas de no retorno, cuando el sentido del flujo es igual al del flujo en estado
permanente 𝐾𝑉𝑡+∆𝑡 es igual al coeficiente de pérdida mínimo de la válvula, en s2/m5, cuando
el agua se regresa hacia la bomba el sistema TransDocencia2-0 le asigna el valor 1012 .
A través de operaciones algebraicas el sistema puede reducirse al sistema de ecuaciones

5.33 y 5.34 cuyas variables son 𝑛𝑃 y 𝑞𝑃 , en el que además de haber sustituido el subíndice
𝑡 + ∆𝑡 por 𝑃, se ha considerado el caso en el que hay 𝑁𝐵 número de bombas operando en
paralelo, las cuales todas dejan de operar simultáneamente (la ecuación 5.32 cambia por
𝑄2𝑡+∆𝑡 = 𝑁𝐵 𝑄𝑅 𝑞𝑡+∆𝑡 ). Este sistema puede resolverse con el método iterativo de Newton-
Raphson.
1 1 𝑛𝑃 𝐶𝑝 𝐶𝑛
𝐾𝑉𝑃 𝑄𝑅2 |𝑞𝑃 |𝑞𝑃 + ( + ) 𝑁 𝑄 𝑞 − 𝐻𝑅 (𝑛𝑃2 + 𝑞𝑃2 ) [𝑎1 + 𝑎2 𝑡𝑎𝑛−1 ( )] − − = 0 (5.33)
𝐶𝑎1 𝐶𝑎2 𝐵 𝑅 𝑃 𝑞𝑃 𝐶𝑎1 𝐶𝑎2
𝜋𝑁𝑅 𝑊𝑅2 𝑛𝑃
{(𝑛𝑃2 + 𝑞𝑃2 ) [𝑎3 + 𝑎4 𝑡𝑎𝑛−1 ( )] + 𝛽} + 𝑛𝑃 − 𝑛 = 0 (5.34)
15 𝑇𝑅 ∆𝑡 𝑞𝑃
donde
𝑁𝐵 número de bombas
Una vez que se han calculado los valores de 𝑛𝑃 y 𝑞𝑃 , las demás variables se calculan
realizando sustituciones en las demás ecuaciones.
Para el caso en el que la planta de bombeo determina la frontera aguas arriba del
acueducto, el sistema de ecuaciones está formado por la característica negativa, ecuación
4.5 y las ecuaciones 5.21, 5.24, 5.25, 5.35 y 5.36, ver figura 5.13.
116
Válvula de descarga
Tramo 1
𝐻𝑆
Punto 1
Figura 5.13 Bomba en la frontera aguas arriba (sin línea de succión)
𝐻1𝑡+∆𝑡 = 𝐻𝑆 + 𝐻𝑅 ℎ𝑡+∆𝑡 − 𝐾𝑉𝑡+∆𝑡 𝑄𝑅2 |𝑞𝑡+∆𝑡 |𝑞𝑡+∆𝑡 (5.35)
𝑄1𝑡+∆𝑡 = 𝑄𝑅 𝑞𝑡+∆𝑡 (5.36)
Se tiene un sistema de seis ecuaciones con las seis incógnitas 𝑄1𝑡+∆𝑡 , , 𝐻1𝑡+∆𝑡 , 𝑞𝑡+∆𝑡 , ℎ𝑡+∆𝑡 ,
𝑛𝑡+∆𝑡 y 𝛽𝑡+∆𝑡 , el cual puede reducirse al sistema de ecuaciones 5.34 y 5.37 y resolverse
con el método iterativo de Newton-Raphson. El valor del coeficiente 𝐾𝑉𝑡+∆𝑡 se determina
como ya se indicó.
𝑛𝑃 𝑁𝐵 𝑄𝑅 𝐶𝑛
𝐻𝑅 (𝑛𝑃2 + 𝑞𝑃2 ) [𝑎1 + 𝑎2 𝑡𝑎𝑛−1 ( )] − 𝐾𝑉𝑃 𝑄𝑅2 |𝑞𝑃 |𝑞𝑃 − 𝑞𝑃 + + 𝐻𝑆 = 0 (5.37)
𝑞𝑃 𝑐𝑎 𝐶𝑎
Una vez determinados los valores de 𝑛𝑃 y 𝑞𝑃 , las demás variables se calculan realizando
sustituciones en las demás ecuaciones.
5.3.2 Simulación en el sistema TransDocencia2-0. Las plantas de bombeo

pueden determinar la frontera aguas arriba del acueducto o localizarse a lo largo de la
conducción, para cada caso existe una pantalla para simularlas.
5.3.2.1 Bombas en la frontera aguas arriba

En la figura 5.14 se muestra la pantalla para dar de alta el paro de emergencia de bombas
localizadas en la frontera aguas arriba de un sistema de bombeo. Esta pantalla se presenta
seleccionando el ícono Paro súbito del bombeo, que se muestra al iniciar la ejecución de
TransDocencia2-0 (figura 3.2) o primero seleccionando la pestaña Fronteras y maniobras
aguas arriba (figura 3.16). En la sección Inspector se capturan los datos para simular la
operación de la bomba en flujo transitorio, el cual contiene las secciones Propiedades
generales, Características de las bombas y Datos de las válvulas descarga.
117
Figura 5.14 Información para simular el paro súbito de bombas localizadas en la frontera
aguas arriba
118
 Nombre, con éste se identifica a la planta de bombeo.
 Número de bombas, es el número de bombas que operan en paralelo en una planta de

bombeo para suministrar el gasto de operación o diseño del acueducto.
 Instante en el que paran las bombas [s], el paro de las bombas puede ser posterior al
inicio de la simulación de la operación del sistema de bombeo, cuando es mayor al
tiempo de simulación las bombas nunca paran.
 Elevación del nivel de succión [msnm], es la elevación de la superficie libre del agua en
el cárcamo de succión o pozos.
 Cadenamiento de la planta de bombeo [m], es la ubicación de la planta de bombeo

respecto a la topografía del terreno con la que se deberán comparar las piezométricas
máximas y mínimas que se presentan durante el tiempo de simulación del transitorio.
 Tipo de válvula en la descarga de las bombas, las válvulas en la descarga de las bombas
pueden ser de cierre controlado o de no retorno (retención o check).
Características de las bombas
 Curvas fh y fb, del catálogo de bombas se elige la velocidad específica más cercana a
la velocidad específica de las bombas del proyecto, la cual se determina a través de la
ecuación 5.11. En el catálogo, las unidades de la velocidad específica denotada con 𝑁𝑆
son (rpm m1.5 ⁄s0.5 m-0.75 ) y las de la velocidad específica denotada con 𝑛𝑆 son (rpm
gpm0.5 pies-0.75 ).
 Gasto (máxima eficiencia) [m3/s], es el gasto con el que opera la bomba con máxima
eficiencia.
 Carga (máxima eficiencia) [m], es la carga que entrega la bomba cuando opera con
máxima eficiencia.
 Velocidad nominal [rpm], es la velocidad de diseño de la bomba.
 Eficiencia máxima, es la eficiencia máxima con la que puede operar la bomba.
 Velocidad específica [Ns], es la velocidad específica de la bomba, se calcula con la

ecuación 5.11.
Con la información de la operación de la bomba con máxima eficiencia se determina el

momento de inercia del conjunto motor-bomba, ecuaciones 5.22 y 5.23. Tanto el momento
119
de inercia, como la velocidad específica calculador por TransDocencia, se despliegan en

la sección de mensajes (figura 3.2).
Bombas en flujo permanente
 Gasto por bomba [m3/s], es el gasto de operación del acueducto entre el número de
bombas.
 Carga de bombeo [m], es la carga que entrega cada bomba en flujo permanente.
 Eficiencia, es la eficiencia con la que trabaja la bomba en flujo permanente.
Datos de las válvulas de descarga
 Cuando la válvula en la descarga de la bomba es de retención se pide.
 Coeficiente mínimo de pérdida [s2/m5], el valor prestablecido es 0.001, sin embargo, éste
puede cambiarse.
 Cuando la válvula de descarga es de cierre programado.
 Tipo de válvula de control, del catálogo de válvulas se selecciona el tipo de válvula de

cierre programado, su acceso es a través de la pestaña Catálogos y después Curvas de
válvulas, figura 3.5.
 Tiempo en el que cierra la válvula [s], es el tiempo que tarda la válvula en pasar de su
posición en flujo estable a totalmente cerrada.
 Instante en el que inicia el cierre de la válvula [s], es el tiempo que transcurre entre el
inicio de la simulación y el inicio del cierre de la válvula. Para simular que la apertura de
la válvula permanece igual a la de flujo permanente durante todo el transitorio (la válvula
no cierra) este dato debe ser mayor que el tiempo de simulación.
 Posición de las válvulas en flujo estable, las válvulas de descarga pueden estar
parcialmente abiertas en flujo permanente. Posición igual a cero significa que la válvula
está completamente cerrada. El valor de la posición cuando está parcialmente abierta
depende del porcentaje o grado de apertura, esta información se obtiene de las gráficas
del comportamiento de cada válvula, misma que se obtiene a través del catálogo de
válvulas.
120
 Diámetro de la válvula [m], es el diámetro de la válvula.
Al término de la simulación, en el archivo de presiones se graba, en mca, la carga

piezométrica aguas abajo de las válvulas de descarga, en el archivo de gastos, el gasto
por bomba en m3/s y en el archivo de velocidades, la velocidad de las bombas en rpm,
todas las variables como función del tiempo en s.
5.3.2.2 Bombas a lo largo de la conducción
En la figura 5.15 se muestra la pantalla para dar de alta el paro de emergencia de bombas
localizadas a lo largo de la conducción.
121
Figura 5.15 Información para simular el paro súbito de bombas localizadas a lo largo de
la conducción
El acceso a la pantalla de la figura 5.15 es a través del ícono Disparo de bombas que se
presenta al seleccionar la pestaña Maniobras entre tramos (figura 3.16). La información
necesaria para simular la operación de la bomba en flujo transitorio es similar al caso en el
que la bomba determina la frontera aguas arriba, excepto que, para las actuales plantas de
bombeo, el paro súbito siempre ocurre al inicio de la simulación, el cadenamiento y la carga
de succión los determina el sistema TransDocencia2-0.
122
Al final de la simulación en el archivo de cargas se graba la carga, en mca, inmediatamente

aguas arriba de las bombas y aguas abajo de las válvulas de descarga; en el de gastos,
en m3/s, el gasto por bomba y en el archivo con sufijo rpm la velocidad de las bombas en
rpm, todas las variables en función del tiempo, en s.
123
EJEMPLOS
6
Ejemplos
6.1 Ejemplo 1
La presentación de este ejemplo tiene varios objetivos.
 Mostrar cómo puede realizarse la simulación de un transitorio producido por el paro

súbito del bombeo cuando se desconoce la curva de operación en flujo transitorio de
la bomba. Esto puede realizarse determinando de manera preliminar la velocidad
específica de la bomba para elegir del catálogo de bombas alguna de ellas.
 Una vez elegida la bomba del catálogo de bombas presentar cómo se describe un
sistema de bombeo en el sistema TransDocencia2-0 y simular el transitorio generado
por la suspensión súbita del bombeo.
 Proponer el tipo y ubicación de los dispositivos de control de transitorios y determinar

sus dimensiones preliminares, para el sistema del ejemplo.
 Indicar cómo puede utilizarse la información que entrega el sistema TransDocencia2-

0 en los archivos de resultados que se crean al final de cada simulación para optimizar
las dimensiones de los dispositivos de control de transitorios.
124
EJEMPLOS
Consideremos que se necesita llevar agua a razón de 0.1 m3/s desde un tanque de succión
a un tanque de entrega a través de una tubería de acero de diámetro de 12”, espesor 0.25”
y rugosidad absoluta 0.05 mm que sigue la topografía que se muestra en la figura 6.1.
El tanque de succión se localiza en el cadenamiento 0+000, el tanque se descarga está en

el cadenamiento 6+221.8 desplantado en la elevación 183.57 msnm, la planta de bombeo
está a 250 m del tanque de succión en la elevación 145.40 msnm. Las cargas de succión
y descarga son 152.58 y 185 m, respectivamente. En la Tabla 6.1 se indican los puntos
que determinan la topografía.
190
Tanque de entrega
180
Elevación, en msnm
170
160
150
Tanque de succión
140
0 800 1600 2400 3200 4000 4800 5600 6400
Distancia, en m
Figura 6.1 Topografía que sigue la conducción para llevar agua del tanque de succión al
tanque de entrega
Tabla 6.1 Topografía
Cadenamiento Elevación Cadenamiento Elevación Cadenamiento Elevación

(m) (msnm) (m) (msnm) (m) (msnm)
0 142.90 1681.8 155.78 3901.8 159.65
164 142.90 1781.8 155.94 4131.8 157.30
164 145.40 2013.3 159.49 4391.8 162.30
250 145.40 2181.8 157.40 4481.8 160.70
384 145.40 2381.8 162.00 4581.8 160.56
384 147.56 2552.8 170.10 4881.8 163.90
440 150.11 2561.8 164.80 4981.8 166.80
641.8 154.15 2781.8 157.97 5181.8 167.91
781.8 152.62 2943.8 162.48 5591.8 181.22
902.8 157.88 3161.8 157.38 5681.8 180.72
1108.8 153.15 3381.8 162.47 5891.8 179.08
1281.8 155.70 3581.8 156.38 5981.8 179.30
1481.8 157.23 3681.8 156.65 6081.8 182.68
6221.8 183.57
125
EJEMPLOS
Para determinar la carga de bombeo necesaria para llevar el agua del tanque de succión
al de descarga, primero determinemos con la ecuación 2.11 el factor de fricción y con la
ecuación 2.3 la constante de pérdida de carga considerando que en la descarga de la
bomba hay una válvula de retención con coeficiente de pérdida de carga igual a 0.001
s2/m5. Previa conversión de unidades al sistema internacional, el cálculo de la velocidad
del agua y del número de Reynolds, se obtiene:
𝑣𝐷 1.37(0.305)
𝑅𝑒 = = = 417850 (6.1)
𝜗 10−6
−2
0.000164 5.74
𝑓 = 1.325 [𝑙𝑛 ( + )] = 0.015 (6.2)
3.7 4178500.9
1 0.015(6221.8) s2
𝐾= + 0.001 = 2926.60 5 (6.3)
2(9.81) 0.305(0.073)2 m
La carga de bombeo (𝐻𝐵 ) necesaria para transportar 0.01 m3/s es la que se necesita para
vencer la pérdida de energía debido a la fricción más la carga estática, en la figura 6.2 se
muestra la curva del sistema.
𝐻𝐵 = 2926.60(0.1)(0.1) + (185 − 152.58) = 61.69 m (6.4)
120
100
80
Carga H, en m
60
40
20
0
0.00 0.05 0.10 0.15 0.20
Gasto Q, en m3/s
Figura 6.2 Curva del sistema
126
EJEMPLOS
La velocidad específica de la bomba se determina a través de la ecuación 5.11, para

obtenerla deben conocerse el punto de operación de la bomba cuando trabaja con la
eficiencia máxima, la velocidad nominal de la bomba, el número de succiones y pasos o
impulsores. Consideremos que el punto de operación con eficiencia máxima coincide con
las condiciones de operación del sistema, 𝑄 = 0.1 m3/s y 𝐻 = 61.69 m y que la bomba es
de succión sencilla, debido a que el gasto no es muy grande.
Para aplicar la ecuación 5.11 falta conocer la velocidad de la bomba y el número de

impulsores, propongamos el número de pasos y velocidades teóricas de giro que se indican
en la Tabla 6.2. Aplicando la ecuación 5.11, como se indica en la ecuación 6.5 para el caso
de una bomba con velocidad 1200 rpm, de 2 pasos y succión simple, se obtienen las
velocidades específicas de la Tabla 6.2.
Tabla 6.2 Velocidades específicas en función del número de pasos y velocidad de la

bomba para 𝑄𝑅 = 0.1 m3/s y 𝐻𝑅 = 61.69 m
𝑁 𝑁𝑆
# pasos
(rpm) (rpm, (m3/s)0.5(m)-0.75)
1 3600 51.72
1 1800 25.86
1 1200 17.24
2 3600 86.98
2 1800 43.49
2 1200 28.99
3 3600 117.90
3 1800 58.95
3 1200 39.30
√0.1
1
𝑁𝑠 = 1200 = 28.99 (6.5)
61.69 0.75
( )
2
El catálogo de bombas contiene curvas características de la operación de bombas en flujo

transitorio para las velocidades específicas 𝑁𝑆 = 29, 35, 38, 44, 46, 147 y 261 (rpm
(m3/s)0.5 (m)-0.75), las velocidades específicas más cercanas de la Tabla 6.2 y del catálogo
son 28.99 y 29, bomba de dos pasos con velocidad 1200 rpm. Para la simulación del
transitorio se usará 𝑁𝑆 = 29.
Otro dato que pide TransDocencia2-0 para simular la operación de las bombas es la
eficiencia máxima (ver figura 5.15) ésta puede estimarse a partir de la figura 6.3 (KSB,
2020).
Para un gasto de 100 l/s y una velocidad específica de 29 (rpm (m3/s)0.5 (m)-0.75), la eficiencia
máxima es aproximadamente 83 %.
127
EJEMPLOS
100
90 más de
630 𝑙/𝑠
80 63 𝑙/𝑠 189 𝑙/𝑠
Eficiencia h
31.5 𝑙/𝑠 630 𝑙/𝑠

70 12.6 𝑙/𝑠
60 6.3 𝑙/𝑠
50
40
10 20 30 40 50 60 80 100 200 300
Velocidad específica, Ns en rpm-(m3/s)0.5-m-0.75
Figura 6.3 Eficiencia máxima esperada de una bomba en función de su velocidad

específica (KSB, 2002)
Un dato que se pide en la descripción de las tuberías es la celeridad, ésta puede

determinarse con la ecuación 2.25, se obtiene:
9
√2.2(10 )
1000
𝑎= = 1199.91 m/s (6.6)
12 2.2(109 )
√1 +
0.25 2(1011 )
A través de las pantallas de las figuras 6.4 a 6.8 se da la información del sistema de
bombeo, la secuencia no necesariamente debe ser la que aquí se presenta.
128
EJEMPLOS
Figura 6.4 Se da de alta el embalse
Figura 6.5 Se da de alta la tubería de succión
129
EJEMPLOS
Figura 6.6 Se da de alta a la bomba
Figura 6.7 Se da de alta a la tubería de descarga
130
EJEMPLOS
Figura 6.8 Se da de alta al tanque de descarga
Una vez que la información del sistema de bombeo está completa, puede realizarse una
simulación seleccionando la pestaña Simulación, figura 3.8. En seguida se despliega la
pantalla de la figura 6.9 para indicar el tiempo de simulación, cuando no hay tanques de
oscilación ni cámaras de aire que generen transitorios de oscilación de masa (transitorios
lentos), un tiempo de simulación de 5 a 10 veces el periodo de la onda de presión (ecuación
2.27) es suficiente, para este caso se propuso 5.8 veces lo que resulta 120 s,
aproximadamente, como se muestra en la figura 6.9.
En esa misma pantalla se muestra el intervalo de integración del método de las

características (∆𝑡) y se indica a TransDocencia2-0 la frecuencia o cada cuántos intervalos
de integración se grabarán los resultados de la simulación.
Cuando el intervalo de integración es muy pequeño y el tiempo de simulación grande, los

archivos de resultados son muy grandes, dependiendo de la capacidad del equipo de
cómputo son difíciles de analizar. Al indicar que la frecuencia de grabación es 2 o 3, el
tamaño de los archivos de resultados son de la mitad o un tercio de los de los archivos
para los que la frecuencia de grabación es igual a uno.
131
EJEMPLOS
Figura 6.9 Se determina el tiempo de simulación
Cuando termina la simulación sin errores de ejecución se presenta la pantalla de la figura

6.10, al seleccionar el ícono OK se presenta el reporte de la figura 6.11.
Este reporte contiene las características físicas del sistema de bombeo, información de la
operación en flujo permanente y las cargas máximas y mínimas en los extremos de las
tuberías que se obtuvieron durante el tiempo de simulación.
132
EJEMPLOS
Figura 6.10 Se informa que terminó la simulación
Figura 6.11 Reporte del sistema de bombeo

133
EJEMPLOS
A través de la pestaña superior del reporte general que dice Envolvente de presiones,
figuras 3.23 y 6.11 se despliega la pantalla de la figura 6.12.
Las gráficas de la figura 6.12 son la línea de cargas piezométricas de operación en flujo
permanente y las líneas de carga piezométricas máximas y mínimas que se obtuvieron en
la simulación. En esta gráfica también puede mostrarse el perfil topográfico, las
coordenadas de la topografía deben estar en un archivo tipo texto que contiene dos
columnas de datos, la primera es el cadenamiento o distancia y la segunda las elevaciones
del terreno. Las columnas deben estar separadas por espacios y no puede haber
encabezados.
Figura 6.12 Piezométricas de cargas máximas, mínimas y operación en flujo permanente
Debido a que la topografía del terreno no es información que utiliza el sistema

TransDocencia2-0, el valor de las cargas mínimas no está limitado por la presión de vapor
del agua. La interpretación que debe darse a las cargas que están por debajo de la
topografía, más de allá de la carga de vaporización es simplemente que habrá problemas
de depresión, la tubería puede colapsar o puede presentarse el fenómeno de separación
de columna.
134
EJEMPLOS
El hecho de que las cargas mínimas estén por debajo del perfil del terreno, ver figura 6.12,
indica que el sistema de bombeo necesita de dispositivos para controlar o limitar las
presiones que se presentan durante el transitorio. A continuación, se presenta el
procedimiento que se siguió para determinar el tipo, ubicación y dimensiones preliminares
de las estructuras hidráulicas que controlan el transitorio.
De la figura 6.1 y la Tabla 6.1 se identifica el punto de coordenadas (2552.8, 170.10) en

éste se propone ubicar un tanque de oscilación, la ubicación de las estructuras hidráulica
se da a través de las longitudes de las tuberías. Aplicando la ecuación 4.38 para determinar
una altura preliminar (𝑍𝑇 ) y la ecuación 4.39 para obtener el área preliminar (𝐴 𝑇 ), se
obtiene:
𝑍𝑇 = 202.82 − 170.10 + 20 = 52.72 m (6.7)
3.1416(3.05)2
𝐴𝑇 = = 7.31 m2 (6.8)
4
La carga piezométrica en flujo permanente en el cadenamiento 2+552.8 puede obtenerse

a través de la ecuación de la energía o aproximar su valor a través de la gráfica
correspondiente de la figura 6.12.
La propuesta de la figura 6.8, puede duplicarse (figura 3.12) para que a partir de ella se
genere una que contenga un tanque de oscilación, en la figura 6.13 se muestra el esquema
de la nueva propuesta, en la misma figura se indican las dimensiones de la torre.
Figura 6.13 Sistema de bombeo con una torre de oscilación
135
EJEMPLOS
Una vez que se realiza la simulación del transitorio, a través de Resultados simulación y
Envolventes de presiones (figura 3.23) se muestran las cargas piezométricas extremas que
se presentan durante 3600 s de simulación, se observa que aguas abajo de la torre de
oscilación, no hay problemas de depresión, figura 6.14.
Figura 6.14 Piezométricas con una torre de oscilación en el cadenamiento 2+552.8
Para evitar las depresiones entre la planta de bombeo y la torre de oscilación se propondrá
un tanque unidireccional en el punto de coordenadas (902.8, 157.88). Al aplicar la
ecuación 4.34 para proponer un área preliminar, se obtiene:
4(2552.8 − 902.8) 0.1

𝐴𝑇 = = 1.803 m2 (6.9)
1200 0.305
A partir de la figura 6.14, observando la línea de cargas piezométricas mínimas y

considerando el sitio de desplante del tanque unidireccional, se propone que la cota de la
elevación de la SLA en el tanque en flujo permanente sea 165 msnm, esto lleva a que el
tirante de agua sea de 165 – 157.88 = 7.12 m.
La propuesta de la figura 6.13 puede duplicarse (figura 3.12) para generar una nueva
propuesta que contenga ambos tanques, el de oscilación y el unidireccional, como se
muestra en la figura 6.15. En la misma figura se muestran las dimensiones del tanque
unidireccional.
136
EJEMPLOS
Figura 6.15 Sistema de bombeo con un tanque unidireccional y una torre de oscilación
Una vez que se realiza la simulación del transitorio pueden observarse las cargas
piezométricas, figura 6.16. En este caso el tiempo de simulación fue de 3600 s, se observan
depresiones en la línea de succión.
Para evitar las depresiones en la línea de succión se propone instalar un tanque de

oscilación, a éste se le llamará tanque de succión, en el punto de coordenadas (200, 145.4).
Aplicando la ecuación 4.38 para determinar la altura preliminar del tanque (𝑍𝑇 ), se obtiene:
𝑍𝑇 = 151.64 − 145.4 + 20 = 26.24 m (6.10)
El área preliminar del tanque de sumergencia es 7.31 m2, según la ecuación 6.8. En la
figura 6.17 se muestra el esquema del sistema de bombeo con la torre de sumergencia, el
tanque unidireccional y la torre de oscilación, en ésta se muestran las dimensiones de la
torre de sumergencia. Ésta nueva propuesta de solución puede crearse duplicando la
propuesta de la figura 6.15.
137
EJEMPLOS
Figura 6.16 Piezométricas con el tanque unidireccional y la torre de oscilación en los

cadenamientos 0+902.8 y 2+552.8, respectivamente
Figura 6.17 Sistema de bombeo con la torre de sumergencia, el tanque unidireccional y

la torre de oscilación
138
EJEMPLOS
En la figura 6.18 se presentan las líneas de carga piezométrica que se obtienen durante
3600 s de simulación, se observa que la torre se sumergencia evita que se presenten
depresiones en la línea de succión.
Figura 6.18 Piezométricas con la torre de sumergencia, el tanque unidireccional y la torre

de oscilación
El sistema TransDocencia2-0 determina la malla de cálculo del método de las

características, figura 2.19, a partir de las celeridades y longitudes de las tuberías. Para el
arreglo de la figura 6.17 el intervalo de integración es de 0.0416 s, como se indica en la
figura 6.19, a diferencia del resto de los arreglos en los que ∆𝑡 = 0.1 s (figura 6.9). Con la
finalidad de reducir el tamaño de los archivos de resultados, los cuales se muestran más
adelante, se eligió una frecuencia de grabación igual a dos.
139
EJEMPLOS
Figura 6.19 Intervalo de integración y frecuencia de grabación de resultados en los

archivos de resultados
Con el arreglo de la figura 6.17 las presiones que se presentan durante el transitorio no
causarán daño a la tubería. Lo que sigue es revisar la operación de cada uno de los
dispositivos de control, la variación de los niveles de agua, los gastos que entregan, la
velocidad de la bomba (aunque en este ejemplo no puede girar en sentido inverso debido
a que la válvula en su descarga es de retención), en las cámaras de aire determinar la
presión máxima a la que estarán sometidas, etc.
6.1.1 Archivos de resultados. Al terminar la simulación, TransDocencia2-0 genera

cinco archivos de resultados, en tres de ellos se graban las variaciones temporales de la
presión, el gasto, la velocidad de las bombas, los volúmenes de agua entregados por los
tanques y los volúmenes de aire admitidos por las válvulas de admisión y expulsión de aire.
En otro se registra las variaciones máximas de la presión a lo largo de la conducción
durante el periodo de simulación, estos cuatro archivos son tipo ASCII con valores
separados por coma, con extensión csv que pueden ser leídos por Excel. En el archivo
restante, tipo texto, se graba la descripción física del sistema de bombeo e información de
las operaciones en estado establecido y transitorio.
Los nombres de los archivos con extensión csv se forman con el nombre de la propuesta
seguidos por los sufijos _HEX, _H, _Q y _RPM, el nombre del archivo tipo texto se forma
con el nombre de la propuesta seguido de _IMP.
El archivo tipo texto contiene la misma información del reporte que se presenta al finalizar
una simulación (figura 6.11).
El archivo _HEX contiene, en función de la distancia, las cargas piezométricas de operación

en flujo permanente y las extremas máximas y mínimas que se obtienen en la simulación.
140
EJEMPLOS
En la Tabla 6.3 se presenta el archivo *_HEX.csv que se genera al terminar la simulación

del transitorio hidráulico producido por la suspensión súbita del bombeo del sistema de la
figura 6.17, la información correspondiente al flujo permanente está completa, la del flujo
transitorio se presenta sólo un segmento. Las gráficas de la figura 6.18, excepto el perfil
del terreno, se obtienen a partir de este archivo.
Para determinar dimensiones preliminares de algunos dispositivos de control de transitorios

es necesario conocer la carga piezométrica en flujo permanente en los sitios en donde se
propone ubicarlos, esta información puede obtenerse a través del archivo *_HEX.csv.
Tabla 6.3 Archivo tipo *_HEX.csv
Flujo transitorio Flujo permanente

Cadenamiento H (máx) H (mín) Cdenamiento H (op. est.)
(m) (mca) (mca) (m) (mca)
0.00 152.580 152.580 0.0 152.580
25.00 152.606 152.463 200.0 151.641
50.00 152.631 152.345 200.0 151.641
75.00 152.656 152.228 250.0 151.406
100.00 152.681 152.110 250.0 213.096
125.00 152.707 151.993 902.8 210.031
150.00 152.732 151.876 902.8 210.031
175.00 152.757 151.758 2552.8 202.282
200.00 152.782 151.641 2552.8 202.282
200.00 152.782 151.641 6221.8 185.053
225.00 155.693 151.311
250.00 159.138 151.157
250.00 250.442 149.798
275.11 250.237 149.995
300.22 249.974 150.210
325.32 249.725 150.444
350.43 249.420 150.695
375.54 249.050 150.966
400.65 248.663 151.259
425.75 248.254 151.575
450.86 247.834 151.918
475.97 247.419 152.290
501.08 247.032 152.689
526.18 246.624 153.120
551.29 246.204 153.588
576.4 245.788 154.089
601.51 245.121 154.581
626.62 244.464 155.113
651.72 243.877 155.636
676.83 243.305 156.207
701.94 242.758 156.831
141
EJEMPLOS
Los archivos _H.csv, _Q.csv y _RPM.csv están formados por columnas de datos, la primera
columna es el tiempo. El número de columnas depende del tipo y número de estructuras
que hay en el sistema de bombeo. Los encabezados de las columnas contienen el nombre
de cada estructura, como ejemplo, en la Tabla 6.4 se muestra un segmento del archivo de
gastos para el arreglo de la figura 6.17.
Tabla 6.4 Segmento del archivo de resultados de gastos (*_Q.csv)
Tiempo Q (Embalse) Q(TS) Q/bomba (bomba) Q (TU) Q (TO) Q (Descarga)

(s) (m3/s) (m3/s) (m3/s) (m3/s) (m3/s) (m3/s)
0.000 0.1000 0.0000 0.1000 0.0000 0.0000 0.1
0.083 0.1000 -0.0050 0.0954 0.0000 0.0000 0.1
0.167 0.1000 -0.0122 0.0880 0.0000 0.0000 0.1
0.250 0.1000 -0.0172 0.0830 0.0000 0.0000 0.1
0.333 0.1000 -0.0211 0.0790 0.0000 0.0000 0.1
0.417 0.1000 -0.0240 0.0761 0.0000 0.0000 0.1
0.500 0.1000 -0.0261 0.0739 0.0000 0.0000 0.1
0.583 0.1000 -0.0278 0.0722 0.0000 0.0000 0.1
0.667 0.1000 -0.0292 0.0708 0.0000 0.0000 0.1
0.750 0.1000 -0.0304 0.0696 0.0000 0.0000 0.1
0.833 0.1000 -0.0314 0.0686 0.0000 0.0000 0.1
0.917 0.1000 -0.0323 0.0677 0.0000 0.0000 0.1
1.000 0.1000 -0.0330 0.0670 0.0000 0.0000 0.1
1.083 0.1000 -0.0336 0.0664 0.0000 0.0000 0.1
1.167 0.1000 -0.0342 0.0658 0.0027 0.0000 0.1
1.250 0.1000 -0.0348 0.0652 0.0055 0.0000 0.1
En la figura 6.20(a) se muestra la variación temporal del nivel del agua en la torre de
sumergencia y en la figura 6.20(b) el gasto que recibe la torre y entrega la misma a la
conducción.
En la figura 6.20(a) se observa que la elevación máxima del nivel del agua no llega a 153
msnm. Debido a que la cota de desplante es 145.4, la altura de la torre puede ser de
153-145.4 = 7.6 m, la altura propuesta en la simulación fue de 26.24 m, ecuación 6.10, el
análisis de estas gráficas permite optimizar las dimensiones de los dispositivos de control
de transitorios.
A través de la información que se muestra en la gráfica de la variación temporal del gasto,

figura 6.20(b) se determina el diámetro de la tubería de conexión de la torre con la
conducción.
142
EJEMPLOS
154
Elevación de la SLA en la TS, en 153

msnm
152
151
150
0 100 200 300 400 500 600
Tiempo, en s
(a)
0.03
0.00
Gasto en la TS, en m3
-0.03
-0.06
-0.09
-0.12
0 100 200 300 400 500 600
Tiempo, en s
(b)
Figura 6.20 (a) Nivel del agua y (b) gasto en la torre de sumergencia
En la primera gráfica de la figura 6.21 se observa cómo el tanque unidireccional entrega

agua a la línea cuando la carga en la conexión disminuye hasta la carga del agua en el
tanque, en la segunda gráfica se observa que el tirante de agua apenas disminuye unos
cuantos centímetros. En la figura 6.18 se observa a las piezométricas mínimas muy cerca
de la topografía, para elevarlas habría que aumentar el tirante de agua en el tanque
unidireccional.
143
EJEMPLOS
240 0.05
Carga en la conexión del TU con la

220 0.04
conducción, en mca
Gasto en el TU, en m3/s

200 0.03
180 0.02
160 0.01
H( TU) Q(TU)
140 0
0 10 20 30
Tiempo, en s
(a)
165.02
Elevación de la SLA en el TU, en
165.00
164.98
msnm
164.96
164.94
164.92
164.90
0 600 1200 1800 2400
Tiempo, en s
(b)
Figura 6.21 (a) Carga en la conexión del tanque unidireccional con la conducción, (b)
gasto que entrega el tanque y variación de la elevación de la SLA
144
EJEMPLOS
En la gráfica de la figura 6.22 se muestra la variación del nivel del agua en la torre de
oscilación, se observa que el nivel del agua en la torre no supera al que tenía en flujo
permanente, 𝑡 = 0, por lo que, la altura de la torre puede ser igual a la carga piezométrica
en flujo permanente menos la cota de desplante de la torre 𝑍𝑇 = 202.8 − 170.10 = 32.70 m
más un bordo, la altura preliminar fue de 52.72 m. El diámetro preliminar fue de 3.05 m,
ésta sería una torre muy esbelta, podría aumentarse el diámetro, pero no disminuir la altura
de 32.70 m más un bordo.
204
Elevación de la SLA en la TO, en msnm
200
196
192
188
184
180
0 600 1200 1800 2400 3000 3600
Tiempo, en s
Figura 6.22 Variación de la SLA en la torre de oscilación y gasto que entrega o recibe de
la conducción
En la figura 6.23(a) se muestra el gasto en la torre de oscilación obtenido durante 3600 s

de simulación, en la figura 6.23(b) se muestra con mayor detalle el gasto durante los
primeros 600 s de la simulación.
En la gráfica de la figura 6.23(b) se observan cambios muy rápidos del gasto, estos se
deben al transitorio por golpe de ariete que se desarrolla entre la planta de bombeo y la
torre de oscilación. La forma de la variación lenta del gasto de la gráfica de la figura 6.23(a)
se debe al transitorio de oscilación de masa que se desarrolla entre la torre de oscilación y
el tanque de entrega.
145
EJEMPLOS
0.15
0.12
Gasto en la TO, en m3 0.09
0.06
0.03
0.00
-0.03
0 600 1200 1800 2400 3000 3600
Tiempo, en s
(a)
0.15
0.12
Gasto en la TO, en m3
0.09
0.06
0.03
0.00
-0.03
0 100 200 300 400 500 600
Tiempo, en s
(b)
Figura 6.23 Variación de la SLA en la torre de oscilación y gasto que entrega o recibe de
la conducción
Cada vez que se modifique cualquier parámetro de los dispositivos de control de

transitorios, su ubicación y tiempos de cierre o aperturas de válvulas, se tiene que volver a
revisar que el diseño del sistema protección de la conducción es adecuado.
146
EJEMPLOS
6.2 Ejemplo 2
La presentación de este ejemplo tiene como objetivos:
 Mostrar cómo elegir del catálogo de bombas la bomba adecuada cuando se conoce la
curva de operación de la bomba en flujo permanente.
 Indicar cómo se debe dar la información de la bomba.
El sistema de bombeo a considerar es el mismo del ejemplo 1, para el cual el gasto a

transportar es 𝑄0 = 0.1 m3/s y la carga de bombeo necesaria es 𝐻𝐵 = 61.68 m.
En la figura 6.24 se muestra la curva de la bomba en la que se ha identificado el punto de

operación del sistema de bombeo (1585.04 gpm, 202.38 pies) y el punto de operación con
máxima eficiencia, 𝜂 = 81 %, es 𝑄𝑅 =1760 gpm = 0.11 m3/s y 𝐻𝑅 = 210 pies = 64.01 m.
Esta bomba es de succión sencilla, de un paso con velocidad de diseño 𝑁𝑅 =1770 rpm.
300
15.00” 40 55
65 70
75
78
13.50” 81.0 %
80
200
78
12.00”
Head-ft
100 Hp
50 Hp
75 Hp
100
0
20
NPSH-ft
0
0 1000 2000
Capacity-USgpm
Figura 6.24 Curva de operación de la bomba en flujo estacionario
A través de la ecuación 5.11 se obtiene la velocidad específica de la bomba, ecuación 6.12.
147
EJEMPLOS
1/2
√0.11 m3
𝑁𝑆 = 1770 = 26.06 rpm ( ) m-0.75 (6.11)
64.010.75 s
La velocidad específica más cercana 26.06 (rpm m1.5 ⁄s0.5 m-0.75 ) del catálogo de bombas
es 𝑁𝑆 = 29, entonces se eligen esas curvas características para realizar la simulación, en
la figura 6.25 se muestra la pantalla con la información de la bomba.
En la sección Características de las bombas debe darse el punto de operación de la bomba

con máxima eficiencia, el campo “Velocidad específica” se refiere a la velocidad específica
de la bomba que operará en el sistema de bombeo, en este caso 26.06, este dato sólo es
informativo. Con esta información se determina el momento de inercia del conjunto motor-
bomba, ecuación 5.22.
En la sección Bombas en flujo permanente debe darse el punto de operación o punto de

diseño del sistema de bombeo y la eficiencia con la opera la bomba en ese punto. Para
este ejemplo, 𝑄 = 0.1 m3/s, 𝐻 = 61.68 m y 𝜂 = 0.8.
Figura 6.25 Información para simular la operación de la bomba
148
BIBLIOGRAFÍA
Bibliografía
Abreu J. M., Guarga F. R., Izquierdo J., et al (1995) Transitorios y oscilaciones en sistemas
hidráulicos a presión, Abreu, Guarga e Izquierdo, España.
Baker A. J. (1985) Finit Element Computational Fluid Mechanics, Mc Graw Hill Book Company,
Singapur.
Cafaggi F. A., Rodal C. E. y Sánchez H. A. (2016) Sistemas de bombeo, Facultad de Ingeniería,

UNAM.
Carmona P. R. (1990) “Transitorios por oscilación de masa en acueductos a presión”, Ingeniería

Hidráulica en México, enero-abril, pp 40-51.
Carmona P. L., Castillo H. G. y Manzano R. S. (2015) Software generador de archivos de datos

para el sistema Trans, Instituto de Ingeniería, UNAM.
Carmona P. L. y Carmona P. R. (2014) Sistema de simulación de transitorios hidráulicos en

conductos a presión: TRANS, Instituto de Ingeniería, UNAM.
Carmona P. L. y Carmona P. R. (2013) Manual del sistema de simulación de transitorios

hidráulicos en tuberías a presión: TRANS, versión 2010.1.0, Serie Manuales, SM14, Instituto de
Ingeniería, UNAM.
Carmona P. R. y Carmona P. L. (1988) “Hydraulic Transients Caused by Pumping Failure.

Models and Comparison Between Numerical Results and Real Measurements”, Symposium,
Trondheim.
Carmona P. R. y Carmona P. L. (1987) “Transitorios por corte de bombeo”, Ingeniería Hidráulica

en México, septiembre-diciembre, pp 53-61.
Carmona P. R., Carmona P. L. y Zarco H. J. A. (2002) “Modelo paramétrico para calcular

dimensiones de tanques de oscilación y de cámaras de aire en sistemas a bombeo”, XX
Congreso Latinoamericano de Hidráulica.
Chaudhry M. H. (2014) Applied hydraulic transients, Van Nostrand Reinhold Company.
149
BIBLIOGRAFÍA
Donsky B. (1961) “Complete Pump Characteristics and the Effects of Specific Speeds on
Hydraulic Transients”, Journal of Basic Engineering, pp 685-699.
Fox J. A. (1977) Hydraulic analysis of unsteady flow in pipe networks, The MacMillan Press Ltd.
Fox P. (1960) “The Solution of Hyperbolic Partial Differential Equations by Difference Methods”.
Ralston, A. y Wilf, H. S. editors, Mathematical Methods for Digital Computers, John Wiley & Sons,
Inc., Nueva York.
Guarga F. R., Sánchez B. J. L., Carmona P. R. y Aguilar M. L. (1985) Diseño y operación

hidráulicos de conducciones de agua a presión. Fascículo II, Llenado, vaciado y operación
estacionaria de las conducciones de agua a presión, Instituto de Ingeniería, UNAM.

hidráulicos de conducciones de agua a presión. Fascículo III, Transitorios hidráulicos en
conductos a presión, Instituto de Ingeniería, UNAM.

hidráulicos de conducciones de agua a presión. Fascículo IV, Control de los transitorios
hidráulicos, Instituto de Ingeniería, UNAM.

hidráulicos de conducciones de agua a presión. Fascículo VII, Programas numéricos, Instituto de
Ingeniería, UNAM.
KSB Compañía Sudamericana de Bombas S.A. (2002) Manual de Entrenamiento, primera

edición.
Sotelo A. G. (2002) Hidráulica General, Limusa, México.
Streeter V. L. y Wylie E. B. (1993) Mecánica de los fluidos, Mc Graw Hill, México.
Swamee P.K. (1993) “Design of a Submarine Oil Pipeline”, Journal of Transportation Engineering,
119(1), pp. 159-170.
Thorley A. (1991) Fluid Transients in Pipeline Systems, D. & L. George Ltd., Inglaterra.
Wylie. E. B. y Streeter V. L. (1978) Fluid Transients, McGraw Hill.
150

Manual Del Sistema TransDocencia 2.1

Cargado por

Información del documento

Descripción 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

Manual Del Sistema TransDocencia 2.1

Cargado por

Copyright:

Formatos disponibles

Manual del sistema de simulación

Libia Georgina Carmona Paredes

Informe final del proyecto PE104017

CAPÍTULO 2. Transitorios hidráulicos en sistemas de bombeo 4

CAPÍTULO 3 Sistema TransDocencia2-0 31

CAPÍTULO 4 Elementos de un sistema de bombeo que no 50

4.5.2 Simulación en el sistema TransDocencia2-0 74

CAPÍTULO 5 Elementos de sistemas de bombeo que generan 96

CAPÍTULO 6 Ejemplos 122

SLA superficie libre del agua

TransDocencia2.0 es una herramienta de apoyo para analizar el comportamiento de

El software que se presenta, TransDocencia2-0, es la versión docencia del Sistema de

Algunas perturbaciones ocurren súbitamente como consecuencia de eventos imprevistos

Cambios súbitos en el consumo de energía eléctrica en la red alimentada por un sistema

El fenómeno del golpe de ariete se analiza a partir de las ecuaciones dinámica y de

En el capítulo 3 se describe el manejo del sistema TransDocencia2-0, formado por dos

En el capítulo 4 se presentan los elementos y estructuras hidráulicas que pueden

En el capítulo 5 se presentan los elementos que tienen la capacidad de producir transitorios

En el capítulo 6 se presentan dos ejemplos, a través de ellos, se muestra cómo se usa el

2.2 Oscilación de masa

Figura 2.1 Después de un cambio de operación del sistema se presentará un transitorio

2.2.1 Modelo matemático de oscilación de masa. El modelo matemático está

Dinámica del agua en la tubería. El estudio de los transitorios de oscilación de masa se

Dinámica de la superficie libre del agua en un tanque. Considerando que puede

𝑑𝑥 𝑥(𝑡 + ∆𝑡) − 𝑥(𝑡)

Ejemplo 1. En la figura 2.3 se muestra un esquema de la conexión en serie de plantas

Figura 2.3 Cuando se suspenden los bombeos se presenta un transitorio de oscilación

Tabla 2.1 Niveles del agua en las torres y gasto en la conducción

Figura 2.4 Variación de los niveles de agua en las torres de oscilación

Tabla 2.2 Elevaciones de las torres de oscilación

Existen varias expresiones para determinar el factor de fricción para flujos

Tabla 2.3 Característica de las tuberías

𝐻𝐵 = 169.94 − 6.7692 𝑄 + 0.8904 𝑄 2 − 0.4942 𝑄 3 (2.12)

Figura 2.6 Curva la bomba

100 512.45 1345 𝑠2

0.012(100) 0.012(512.45) 0.012(1345) 𝑠2

El nivel del agua en la torre de sumergencia se obtiene resolviendo el sistema de

Sustituyendo la ecuación 2.17 en la 2.19 y la que resulta en la ecuación 2.16, se obtiene la

Antes de poner en marcha a la bomba los gastos 𝑄1 = 𝑄2 = 𝑄3 = 0, los niveles de agua

de la oscilación del tirante de agua disminuye conforme transcurre el tiempo debido a la

Nivel del agua en la torre HT, en msnm 1777 4

2.3 Golpe de ariete

Los transitorios rápidos o de golpe de ariete se caracterizan porque los cambios en la

Como consecuencia de un golpe de ariete, en el interior de algunas secciones de la

 Sobrepresiones, presión mayor que la de operación en régimen estacionario. Si ésta es

 Depresiones, presión menor que la presión atmosférica. Las depresiones favorecen el

2.3.1 Modelo matemático del golpe de ariete. En el modelo de golpe de ariete

Con la finalidad de imaginarnos lo que sucede en el interior de la conducción cuando se

En 𝑡 = 0 se suspende súbitamente el bombeo y se cierra la válvula de descarga.

Figura 2.10 En 𝑡 = 0 se suspende súbitamente el bombeo, se genera una onda de baja

Figura 2.12 Entre 𝑡 = 𝑙 ⁄𝑎 y 𝑡 = 2𝑙 ⁄𝑎, se restablece la presión que había en flujo

En 𝑡 = 2𝑙 ⁄𝑎, la presión en el interior de la tubería es igual a la que había en flujo

Figura 2.13 En 𝑡 = 2𝑙 ⁄𝑎 la presión en toda la conducción es igual a la de flujo

En 𝑡 = 2𝑙 ⁄𝑎 la onda de presión llega a la válvula de descarga completamente cerrada,

En 𝑡 = 3𝑙 ⁄𝑎 la onda de alta presión habrá recorrido toda la conducción, figura 2.15.

Figura 2.15 En 𝑡 = 3𝑙 ⁄𝑎 la onda de alta presión termina de recorrer la tubería

Figura 2.16 Entre 𝑡 = 3𝑙 ⁄𝑎 y 𝑡 = 4𝑙 ⁄𝑎, se restablece la presión del flujo permanente

En 𝑡 = 4𝑙 ⁄𝑎, la onda reflejada en el tanque de descarga llega a 𝑥 = 0, se concluye el

Figura 2.17 En 𝑡 = 4𝑙 ⁄𝑎 se termina el primer ciclo de la onda de presión

En los comentarios anteriores únicamente se ha hecho referencia a los cambios de la

G3 Figura 2.18b Cambio del gasto en función del tiempo y la posición

Las relaciones entre la carga y el gasto en función de la posición y el tiempo se obtienen a

El área transversal de la conducción (𝐴) y el diámetro (𝐷) son características de la