Trabajo 4

Universidad Nacional
de Huancavelica
(Creada por Ley Nro. 25265)
Facultad de Ciencias de
Ingeniería
Escuela Profesional de
Ingeniería Civil-Hvca
CURSO:
ABASTECIMIENTO DE AGUA Y
ALCANTARILLADO
TEMA:
APLICACIÓN DE FENÓMENOS
TRANSITORIOS EN EL LENGUAJE DE
PROGRAMACIÓN C++
DOCENTE:
Ing. AYALA BIZARRO Iván
ESTUDIANTE:
SEDANO SALAZAR Bacilio
Huancavelica-Agosto 2020
‹‹‹‹‹
Índice General
1. Introducción 2
2. Objetivos 3
1. Objetivos General . . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
2. Objetivos Específicos . . . . . . . . . . . . . . . . . . . . . . . . . . . 3
3. Golpe de Ariete 4
1. Causas del Fenómeno . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
2. Efectos del Fenómeno . . . . . . . . . . . . . . . . . . . . . . . . . . . 6
3. Métodos de Solución . . . . . . . . . . . . . . . . . . . . . . . . . . . 7
4. Método de las Características . . . . . . . . . . . . . . . . . . . . . . 7
5. Solución del Sistema . . . . . . . . . . . . . . . . . . . . . . . . . . . 8
6. Condiciones de Borde . . . . . . . . . . . . . . . . . . . . . . . . . . . 9
1. Reservorio Aguas Arriba . . . . . . . . . . . . . . . . . . . . . 10
2. Válvulas aguas abajo . . . . . . . . . . . . . . . . . . . . . . . 10
3. Unión de Tuberías . . . . . . . . . . . . . . . . . . . . . . . . 11
4. Bifurcación . . . . . . . . . . . . . . . . . . . . . . . . . . . . 12
4. Metodología de Programación 13
1. Ejercicio N 01 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 13
2. Solución . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 15
3. Ejercicio N° 02 . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 25
4. Solución . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 26
5. Conclusiones 35
1
1 Introducción
El golpe de ariete se refiere a las fluctuaciones causadas por un repentino incremento

o disminución de la velocidad de flujo. Estas fluctuaciones de presión pueden ser lo
suficientemente severas para romper la tubería de agua.
Las ondas de presión generadas en este transitorio se propagarán por el sistema a

una velocidad determinada, de tal modo que el fenómeno sería una similitud al de
la propagación y reflexión de las ondas acústicas, en el cual estas ondas viajan a la
conocida velocidad del sonido.
Los problemas potenciales del golpe de ariete pueden ser considerados al evaluarse el
diseño de las tuberías y cuando se realiza un análisis detallado de las oscilaciones de
presión, en muchos casos para evitar malos funcionamientos costosos en el sistema
de distribución. Cualquier cambio mayor en el diseño del sistema principal o cambio
de la operación tales como el aumento en la demanda de los niveles de flujo deben
incluir la consideración de los problemas potenciales del golpe de ariete.
2
2
Objetivos
2.1 Objetivos General

1 Realizar una aplicación en lenguaje C++ para el fenómeno de Golpe de Ariete
en tubería a presión.
2.2 Objetivos Específicos

1 Realizar la selección de tuberías diseñados y calculados como flujo permanente
en condiciones iniciales.
2 Desarrollar un procedimiento de cálculo para el análisis del fenómeno el cual

permita la programación del mismo en lenguaje de programación.
3 Comparar los parámetros hidráulicos en condiciones iniciales (Flujo Perma-

nente) e incorporando el análisis de los fenómenos transitorios (Flujo no Per-
manente).
3
3
Golpe de Ariete
El golpe de ariete (choque hidráulico) es el incremento momentáneo en presión, el

cual ocurre en un sistema de agua cuando hay un cambio repentino de dirección
o velocidad del agua. Cuando una válvula de rápido cierre cierra repentinamente,
detiene el paso del agua que está fluyendo en las tuberías, y la energía de presión es
transferida a la válvula y a la pared de la tubería. Las ondas expansivas se activan
dentro del sistema. Las ondas de presión viajan hacia atrás hasta que encuentran
el siguiente obstáculo sólido, luego continúan hacia adelante, luego regresan otra
vez. La velocidad de las ondas de presión es igual a la velocidad del sonido; por lo
tanto, su ?explosión? a medida que viaja hacia adelante y hacia atrás, hasta que se
disipa por la pérdida de fricción. Cualquiera que haya vivido en una casa antigua
está familiarizado con la ?explosión? que resuena a través de las tuberías cuando
una llave de agua es cerrada repentinamente. Esto es un efecto del golpe de ariete.
Los siguientes son algunos casos en que se puede presentar golpe de ariete:
1 Cambios en la abertura de la válvula, accidental o planeado.
2 Arranque o interrupción de bombas.
3 Cambios en la demanda de potencia de turbinas.
4 Cambios de elevación del embalse.
5 Ondas en el embalse.
6 Vibración de impulsores en bombas, ventiladores o turbinas.
7 Variaciones en la apertura o cierre del gobernador o regulador de una turbina
causadas por cambios en la carga de los sistemas eléctricos.
8 Vibración de accesorios deformables tales como válvulas.
El golpe de ariete puede ser positivo o negativo de acuerdo a la forma en que se
produzca. Cuando se cierra súbitamente una válvula se presenta un golpe de ariete
positivo o sobrepresió. El golpe de ariete negativo o subpresión ocurre al efectuarse
la apertura brusca de una válvula.
En la siguiente figura 3.1 representa la variación de la presión en la válvula de
cierre rápido sin considerar pérdidas por fricción debido al desplazamiento de la
onda de sobrepresión a lo largo del conducto. El ciclo de las ondas de sobrepresión y
subpresión se repetiría indefinidamente pero las pérdidas de energía hacen que vaya
atenuándose hasta anularse por completo tal como se ilustra en la Figura 3.2.
4
Ingeniería Civil 3 Golpe de Ariete
Figura 3.1: Variación real de la presión en la válvula, considerando pérdidas de

fricción.
Figura 3.2: Variación teórica de la presión en la válvula, sin considerar pérdidas de

fricción.
C : velocidad de propagación de la onda de sobrepresión
L : longitud de la tubería
h’ : sobrepresión en la chimenea de equilibrio.
-h’ : subpresión en la chimenea de equilibrio.
hf : perdidas por fricción en la tubería
V : velocidad de la tubería.
Pág. 5
Abatecimiento de Agua y Alcantarillado
3.1 Causas del Fenómeno

1 El fenómeno comúnmente ocurre cuando una válvula está cerrada rápidamente
y repentinamente detiene el flujo de líquido en una tubería.
2 El arranque de la bomba puede inducir un colapso rápido del espacio vacío
que existe aguas abajo de la bomba.
3 El cierre rápido de la bomba puede crear un cambio rápido en el flujo, lo cual
causa un despliegue violento de presión en el lado succionante y una presión
abajo de oleada en el lado de descarga.
4 Ondas en el embalse.
Figura 3.3: Cierre brusco de la Válvula
3.2 Efectos del Fenómeno

1 La ruptura de la tubería
2 El daño para bombas, fundaciones, interiores de la tubería y válvulas.
3 Los aditamentos dañados de la tubería.
4 Un acontecimiento externo a la instalación.
Pág. 6
Figura 3.4: Efectos dañinos en la tubería
3.3 Métodos de Solución

1 Método de las Características
2 Método The Lax-Wendroff One-Step Method
3 Método de diferencias finitas.
4 Método de Elementos Finitos (MEF)
5 Métodos Híbridos (MH.s)
3.4 Método de las Características

El Método de las Características ha sido el método numérico explícito mas utilizado
para el modelo de fenómeno de propagación de ondas en las redes de tubería, debido
a su facilidad para introducir diferentes dispositivos y condiciones limite (bombas,
válvulas, tanques, etc.). Por esta razón, conducto cerrado, canal abierto y las corrien-
tes de aguas subterráneas se han analizado utilizando esta técnica. La popularidad
del método se deriva en sus principales atributos (Goldberg y Wylie,1983)
Se aplican dos ecuaciones básicas de la mecánica a un elemento de fluido en una

tubería para obtener las ecuaciones diferenciales del flujo transitorio: la ecuación de
conservación de momento lineal y la ecuación de conservación de la masa.
Pág. 7
Ecuación de Conservación de Momento Lineal:
BQ g.A.BH f
L1 “
` ` Q.|Q| (3.1)
Bt Bx 2.D.A
Ecuación de Conservación de la Masa.
BQ g.A.BH
L2 “ a2 ` “0 (3.2)
Bx Bt
Donde la a en la ecuación 3.2 es la velocidad de la onda de presión y depende de la
compresibilidad del fluido, la rigidez de la cañería y las propiedades mecánicas del
material. La misma puede calcularse como:
K
ρ
2
a “ ` K.D ˘C1 (3.3)
1` E.e
Donde:
K = Módulo de compresibilidad.
ρ = Densidad del Fluido
E = Módulo de elasticidad
e = Espesor de la Tubería
3.5 Solución del Sistema

El sistema de ecuaciones considerando este formato con dos ecuaciones en derivadas
parciales en donde las variables dependientes son H y Q. Estas ecuaciones pueden
combinarse linealmente mediante un multiplicador aún desconocido λ:
L “ λL1 ` L2 “ 0 (3.4)
Donde L1 y L2 son las ecuaciones de equilibrio dinámico 3.1 y de continuidad 3.2,
respectivamente. Cualquier par de valores de distintos λ va a dar un sistema de
ecuaciones equivalente a la combinación de 3.1 y 3.2. La solución por el método de las
características consiste en escoger un par de valores para λ que permitan transformar
las ecuaciones en derivadas parciales en un par de ecuaciones diferenciales totales,
que puedan ser resueltas numéricamente Al reemplazar las ecuaciones 3.1 y 3.2 en ,
se obtiene:
ˆ ˙ ˆ ˙
BQ BQ BH BH f
` λ.a2 . ` λ.g.A ` ` Q. |Q| “ 0 (3.5)
Bt Bx Bt λBx 2.D.A
Pág. 8
La ecuación característica, puede ser expresada en un esquema de diferencias finitas

como:
g.A f.∆t
pQp ´ QA q ` pHP ´ HA q ´ QA .|QA | “ 0 (3.6)
a 2.D.A
g.A f.∆t
pQp ´ QB q ´ pHP ´ HB q ´ QB . |QB | “ 0 (3.7)
a 2.D.A
La ecuación puede ser expresada de la siguiente manera:
Qp “ CP ´ Ca HP ^ Qp “ Cn ´ Ca HP (3.8)
Donde:
CP “ QA ` Ca HA ´ Cf QA |QA | (3.9)
Cn “ QB ` Ca HB ´ Cf QB |QB | (3.10)
Los valores Ca y Cf son valores constantes.
g.A f.∆t
Ca “ ^ Cf “ (3.11)
a 2D.A
1
QP “ pCP ` Cn q (3.12)
2
CP ´ QP
HP “ (3.13)
Ca
3.6 Condiciones de Borde

La solución completa de las ecuaciones transitoria gobernantes requiere calcular
algunas de las condiciones iniciales, así como algunas condiciones límite. En gene-
ral, las condiciones iniciales se especifican a partir de una condición conocida, por
ejemplo: estado de equilibrio de flujo. Por otro lado, las condiciones límite debe cal-
cularse por las expresiones que pueden establecer alguna relación entre las variables
de estado H y Q presente en la frontera.
Pág. 9
3.6.1 Reservorio Aguas Arriba
Cn “ QB ´ Ca HB ´ Cf QB |QB | (3.14)
la presión en el nudo.
HPi,1 “ Hres (3.15)
El caudal que circula en el instante.
QPi,1 “ Cni ` Cai Hres (3.16)
3.6.2 Válvulas aguas abajo
Pág. 10
CP “ QA ` Ca HA ´ Cf QA |QA (3.17)
El coeficiente de válvula para cada intervalo de tiempo.
pτ Qo q2
Cv “ (3.18)
Ca Ho
1´ a ¯
QPi,n`1 “ Ćv ` Cv2 ` 4CP .Cv (3.19)
2
la presión en el nudo.
CP ´ QPi,n`1
HPi,n`1 “ (3.20)
Ca
3.6.3 Unión de Tuberías
Esta ecuación característica positiva corresponde a la tubería aguas arriba, la cual

es.
CPi “ QA ` Ca HA ´ Cf QA |QA (3.21)

Esta ecuación característica negativa corresponde a la tubería aguas abajo, la cual
es.
Cni`1 “ QB ´ Ca HB ´ Cf QB |QB | (3.22)
La presión en el nudo.
CPi ´ Cni`1
HPi,n`1 “ (3.23)
Cai ` Cai`1
QPi,n`1 “ CPi ´ Cai HPi,n` (3.24)
Pág. 11
3.6.4 Bifurcación
CPi “ QA ` Cai HA ´ Cfi QA |QA (3.25)
Cni`1 “ QB ´ Cai`1 HB ´ Cfi`1 QB |QB | (3.26)
Cni`2 “ QC ´ Cai`2 HC ´ Cfi`2 QC |QC | (3.27)

La presión en el nudo.
CPi ´ Cni`1 ´ Cni`2

HP “ (3.28)
Cai ` Cai`1 ` Cai`2
QPi “ CPi ´ Ca HPi (3.29)
QPi`1 “ Cni`1 ` Cai`1 HPi`1,1 (3.30)
QPi`2 “ Cni`2 ` Cai`2 HPi`2,2 (3.31)
Pág. 12
4
Metodología de Programación
La programación esta realizado en el lenguaje de programación C++, para lo cual

se tendrá dos sistema a analizar, la primera será el mismo ejercicio desarrollado en la
clase del libro de Chaudhry y la segunda será un ejercicio propuesta del estudiante
que incluirá la bifurcación.
La programación esta realizada para i tuberías y n tramos.
4.1 Ejercicio N 01
Figura 4.1: Sistema de tubería
13
Ingeniería Civil 4 Metodología de Programación
Figura 4.2: Curva de cierre de válvula
t [s] τ
0 1.0
1 0.9
2 0.7
3 0.5
4 0.3
5 0.1
6 0
Tabla 4.1: Datos del Tiempo vs Tao
Pág. 14
4.2 Solución
La solución del ejercicio será aplicando por el método de las características y utilizan-
do la aplicación desarrollado por el alumno en lenguaje c++. Para el primer ejercicio
se utilizará la Primera aplicación desarrollado con nombre Analisis ´ Transtorioy
el segundo Analisis ´ Transtorio ´ Bifurcacion que se adjuntara junto con el
informe.
1 Datos Generales
1 // I n g r e s o de Datos G e n e r a l e s
2 f l o a t Hres , Qo , Tmax , Tc ;
3
4 Hres = 6 7 . 7 ; // A l t u r a de R e s e r v o r i o [m]
5 Qo = 1 ; // Caudal [ m3/ s ]
6 Tmax = 5 ; // Tiempo de C a l c u l o [ s ]
7 Tc = 6 ; // Tiempo de c i e r r e de l a V a l v u l a [ s ]
8 f l o a t L [ ] = {550 ,450}; // Longitud de l a s T u b e r i a s [m]
9 float D[ ] = {.75 , .60}; // Diametro de l a Tuberia [m
10 float f [ ] = {.010 ,.012}; // C o e f i c i e n t e de f r i c c i o n
11 f l o a t a [ ] = {1150 ,900}; // C e l e r i d a d [m/ s ]
2 Datos de Tiempo de operación de la Válvula vs Toa

1 //
´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´
2 f l o a t Tm[ ] = { 0 , 1 , 2 , 3 , 4 , 5 , 6 } ; // Tiempo de o p e r a c i o n
de V a l v u l a
3 f l o a t Tao [ ] = { 1 , 0 . 9 , 0 . 7 , 0 . 5 , 0 . 3 , 0 . 1 , 0 } ; // P r o c e s o de C i e r r e y
Apertura de V a l v u l a
3 Cálculo de valores constantes Los valores Ca y Cf se calculará mediante la

ecuación 3.11 y la variación del tiempo se calculara con la siguiente formula:
L
∆t “ (4.1)
N tramos.a
1 // C a l c u l o de V a l o r e s C o n s t a n t e s
2
3 f l o a t Dt [ l (L) ] ;
4 f l o a t Dx [ l (L) ] ;
5 f l o a t A [ l (D) ] ;
6
7 for ( int i =0 ; i <l (L) ; i ++){

8 Dt [ i ] = (L [ i ] ) / (N∗a [ i ] ) ;
9 Dx [ i ] = L [ i ] /N;
10 A [i ] = ( p i / 4 ) ∗D[ i ] ∗D[ i ] ;
11 }
12
13
Pág. 15
4 Resultados de los valores Constantes
El ejercicio se calculará en tres tramos como se muestra en la figura 4.3.
Resultados
Figura 4.4: Resultados de los valores constantes
Pág. 16
5 Cálculo de Cargas Hidráulicas en cada punto
Para calcular las perdidas de carga se utilizará la ecuación de Darcy-Weisbach4.2.
8∆xQo2
Hf “ f (4.2)
gπ 2 D5
1 // C a l c u l o de Carga H i d r a u l i c a en cada punto

2
3 i n t f i l a s = l (L) ;
4 i n t columnas = N+1;
5 f l o a t Ht=Hres ;
6 f l o a t hh [ f i l a s ∗ columnas ] ;
7 f l o a t c_f [ f i l a s ∗ columnas ] ;
8 f l o a t c_a [ f i l a s ∗ columnas ] ;
9
10 f o r ( i n t i =0; i < f i l a s ; i ++){

11 f o r ( i n t j =0; j <columnas ; j ++){
12 Hf [ i ] [ j ] = Ht ´ Hf [ i ] [ j ] ;
13 }
14 Ht = Hf [ i ] [ columnas ´ 1 ] ;
15 }
Resultados de carga hidráulica en estado permanente en cada nodo.
Figura 4.5: Carga Hidráulica en estado permanente
Pág. 17
6 Proceso de Interpolación con Lagrange
La interpolación de polinomios de Lagrange es una reformulación del polinomio

de Newton que evita el cálculo de la tabla de diferencias, el polinomio de
Lagrange se expresa como:
ÿ
fn pxq “ Li pxqf pxi q (4.3)
i“0
Donde: n
ź x ´ xi
Li pxq “ (4.4)
j“0
pxi ´ xj q
j‰i
El polinomio de Lagrange de primer grado se calcula con la siguiente fórmula

4.5.
x ´ x1 x ´ x0
f1 pxq “ f px0 q ` f px1 q (4.5)
x0 ´ x 1 x1 ´ x0
Algoritmo del Método de interpolación de Lagrange en C++
1 // P r o c e s o para I n t e r p o l a r Lagrange
2 float x i n t e r p o l [ cc ] , yinterpolada [ cc ] ;
3 f l o a t productosNum , productosDen ;
4 yinterpolada [ 0 ] = 0;
5
6 f o r ( i n t i =0; i <c c ; i ++){

7 xinterpol [ i ] = x [ i ] ;
8 }
9
10 f o r ( i n t i i =0; i i <c c ; i i ++){

11 f o r ( i n t i =0; i <l (Tm) ; i ++){
12 productosNum =1;
13 productosDen = 1 ;
14 f o r ( i n t j =0; j <l ( Tao ) ; j ++){
15 i f ( i != j ) {
16 productosNum = productosNum ∗ ( x i n t e r p o l [ i i ]´Tm[ j ] ) ;
17 productosDen = productosDen ∗ (Tm[ i ] ´ Tm[ j ] ) ;
18 }
19 }
20 y i n t e r p o l a d a [ i i ] = y i n t e r p o l a d a [ i i ] + ( productosNum /
productosDen ) ∗Tao [ i ] ;
21 }
22 }
23 p r i n t f ( " \n" ) ;
Pág. 18
Resultados compilados Tiempo y los valores interpolados de Tao.
Figura 4.6: Resultados interpolados del Tao
7 Cálculo en las condiciones de borde del reservorio.

1 i f ( j ==0){
2 p r i n t f ( " \n" ) ;
3 p r i n t f ( " \ t Borde [ %i , %i ] " , i , j ) ;
4 p r i n t f ( " \n" ) ;
5
6 Cn=Qt [ i ´ 1 ] [ j +1]ća [ j ] ∗ Pt [ i ´ 1 ] [ j +1]´ c f [ j ] ∗ Qt [ i ´ 1 ] [ j +1]∗ abs ( Qt

[ i ´ 1 ] [ j +1]) ;
7 Pt [ i ] [ j ]=Hh [ j ] ;
8 Qt [ i ] [ j ]=Cn+ca [ j ] ∗ Pt [ i ´ 1 ] [ j + 1 ] ;
9 }
8 Cálculo en las condiciones de válvula .

1 e l s e i f ( j ==( l (Hh) ´1) ) {
2
3 p r i n t f ( " \ t V a l v u l a [ %i , %i ] " , i , j ) ;
4 p r i n t f ( " \n" ) ;
5 i f ( menor∗ i <Tc ) {
Pág. 19
6 tao = yi n te r po l ad a [ i ] ;
7 }
8 else {
9 tao = 0 ;
10 }
11 p r i n t f ( " \ t Tao %f 0 . 5 " , t a o ) ;
12 p r i n t f ( " \n" ) ;
13
14 Cp=Qt [ i ´ 1 ] [ j ´1]+ca [ j ] ∗ Pt [ i ´ 1 ] [ j ´1]´ c f [ j ] ∗ Qt [ i ´ 1 ] [ j ´1]∗ abs ( Qt [ i

´ 1 ] [ j ´1]) ;
15 Cv=pow ( t a o ∗Qt [ i ´ 1 ] [ j ´ 1 ] , 2 ) / ( ( ca [ j ] ) ∗Pt [ i ´ 1 ] [ j ´1]) ;
16 Pt [ i ] [ j ]=(Cp´Qt [ i ´ 1 ] [ j ´1]) / ( ca [ j ] ) ;
17 Qt [ i ] [ j ]=.5∗( Ćv+pow ( pow (Cv , 2 ) +4∗Cv∗Cp , . 5 ) ) ;
18 p r i n t f ( " \ t ´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´\n" ) ;
19 p r i n t f ( " \ t i e s : %i j %i ", i ,j);
20 p r i n t f ( " \n" ) ;
21 p r i n t f ( " \ t ´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´\n" ) ;
9 Cálculo en las condiciones de empalme o unión .

1 e l s e i f ( j %(N)==0 ) {
2 p r i n t f ( " \ t Empalme [ %i , %i ] " , i , j ) ;
3 p r i n t f ( " \n" ) ;
4
5 Cp=Qt [ i ´ 1 ] [ j ´1]+ca [ j ´1]∗ Pt [ i ´ 1 ] [ j ´1]´ c f [ j ´1]∗Qt [ i ´ 1 ] [ j ´1]∗ abs (

Qt [ i ´ 1 ] [ j ´1]) ;
6
7 Cn=Qt [ i ´ 1 ] [ j +1]ća [ j +1]∗ Pt [ i ´ 1 ] [ j +1]´ c f [ j +1]∗Qt [ i ´ 1 ] [ j +1]∗ abs (

Qt [ i ´ 1 ] [ j +1]) ;
8
9 Pt [ i ] [ j ]=(CpĆn) / ( ca [ j ´1]+ca [ j +1]) ;

10
11 Qt [ i ] [ j ]=Cpća [ j ] ∗ Pt [ i ] [ j ] ;
10 Cálculo en las condiciones de tramos interiores .

1 else {
2 p r i n t f ( " \ t Tramos I n t e r i o r e s [ %i , %i ] " , i , j ) ;
3 p r i n t f ( " \n" ) ;
4

´ 1 ] [ j ´1]) ;
6
7 Cn=Qt [ i ´ 1 ] [ j +1]ća [ j ] ∗ Pt [ i ´ 1 ] [ j +1]´ c f [ j ] ∗ Qt [ i ´ 1 ] [ j +1]∗ abs ( Qt [ i

´ 1 ] [ j +1]) ;
8
9 Qt [ i ] [ j ] = . 5 ∗ (Cp+Cn) ;
10
11 Pt [ i ] [ j ]=(Cp´Qt [ i ] [ j ] ) / ca [ j ] ;
12 }
Pág. 20
11 Reconocimiento de la condiciones de cada Iteración.
Figura 4.7: Reconocimiento de la condiciones de cada Iteración
Pág. 21
Pág. 22
12 Resultados Obtenidos de las Presiones en cada Nodo.
Figura 4.9: Resultados de Presión en cada nodo
Pág. 23
13 Resultados Obtenidos de los Caudales en cada Nodo.
Figura 4.10: Resultados de Caudal en cada tramo
Pág. 24
4.3 Ejercicio N° 02
Figura 4.11: Sistema de Tuberías

Datos del ejercicio:
Tramo [s] Longitud [m] Diámetro [m] Fricción Celeridad [m/s]

0 450 0.70 0.012 1400
1 600 0.80 0.010 1300
2 500 0.65 0.011 1200
3 650 0.55 0.013 900
Tabla 4.2: Datos del Ejercicio
Para interpolar del Tao se utilizará mismo datos del ejercicio N° 01.
Figura 4.12: Curva de cierre de válvula
Pág. 25
t [s] τ
0 1.0
1 0.9
2 0.7
3 0.5
4 0.3
5 0.1
6 0
Tabla 4.3: Datos del Tiempo vs Tao
4.4 Solución
La solución del ejercicio será aplicando el método de las características y utilizando
la aplicación desarrollado por el alumno en lenguaje c++.
1 Datos Generales
EL Caudal de diseñará el primer tramo con 1 [L/s] y los dos tramos con la
mitad que equivale 0.5 [L/s].
1 // I n g r e s o de Datos G e n e r a l e s
2
3 f l o a t Hres , Qo , Tmax , Tc ;
4
5 Hres = 6 7 . 7 ; // Altura de R e s e r v o r i o [m]

6 Qo = 1 ; // Caudal [ m3/ s ]
7 Tmax = 5 ; // Tiempo de C a l c u l o [ s ]
8 Tc = 6 ; // Tiempo de c i e r r e de l a V a l v u l a [ s ]
9
10 float L[] = {450 ,600 ,500 ,650}; // Longitud de l a s T u b e r i a s [m]

11 float D[ ] = {.7 , .80 , .65 ,.55}; // Diametro de l a Tuberia [m]
12 float f [] = { . 0 1 2 , . 0 1 0 , . 0 1 1 , . 0 1 3 } ; // C o e f i c i e n t e de f r i c c i o n
13 float a[] = {1400 ,1300 ,1200 ,900}; // C e l e r i d a d [m/ s ]
2 Datos de Tiempo de operación de la Válvula vs Toa

1 //
´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´´
2 f l o a t Tm[ ] = { 0 , 1 , 2 , 3 , 4 , 5 , 6 } ; // Tiempo de o p e r a c i o n
de V a l v u l a
3 f l o a t Tao [ ] = { 1 , 0 . 9 , 0 . 7 , 0 . 5 , 0 . 3 , 0 . 1 , 0 } ; // P r o c e s o de C i e r r e y
Apertura de V a l v u l a
3 Cálculo de valores constantes
Pág. 26
Los valores Ca y Cf se calculará mediante la ecuación 3.11 y la variación del

tiempo se calculara con la siguiente formula:
L
∆t “ (4.6)
N tramos.a
1 // C a l c u l o de V a l o r e s C o n s t a n t e s
2
3 f l o a t Dt [ l (L) ] ;
4 f l o a t Dx [ l (L) ] ;
5 f l o a t A [ l (D) ] ;
6
7 for ( int i =0 ; i <l (L) ; i ++){

8 Dt [ i ] = (L [ i ] ) / (N∗a [ i ] ) ;
9 Dx [ i ] = L [ i ] /N;
10 A [i ] = ( p i / 4 ) ∗D[ i ] ∗D[ i ] ;
11 }
4 Resultados de los valores Constantes
El ejercicio se calculará en dos tramos como se muestra en la figura 4.13.
Pág. 27
Resultados
Figura 4.14: Resultados de los valores constantes
5 Cálculo de Cargas Hidráulicas en cada nodo
Para calcular las perdidas de carga se utilizará la ecuación de Darcy-Weisbach4.7.
8∆xQo2
Hf “ f (4.7)
gπ 2 D5
1 // C a l c u l o de Carga H i d r a u l i c a en cada nodo

2
4 i n t f i l a s = l (L) ;
5 i n t columnas = N+1;
6 f l o a t Ht=Hres ;
7 f l o a t hh [ f i l a s ∗ columnas ] ;
8 f l o a t c_f [ f i l a s ∗ columnas ] ;
9 f l o a t c_a [ f i l a s ∗ columnas ] ;
10
Pág. 28
11 f o r ( i n t i =0; i < f i l a s ; i ++){

12 f o r ( i n t j =0; j <columnas ; j ++){
13 Hf [ i ] [ j ] = Ht ´ Hf [ i ] [ j ] ;
14 }
15 i f ( i==columnas ) {
16 Ht = Hf [ 0 ] [ columnas ´ 1 ] ;
17 else {
18 Ht = Hf [ 0 ] [ columnas ´ 1 ] ;
19 }
20 }
Resultados de carga hidráulica en estado permanente en cada nodo.
Figura 4.15: Carga Hidráulica en estado permanente
6 Proceso de Interpolación con Lagrange
La interpolación de polinomios de Lagrange,al igual que se calculó en primer

ejercicio, se calculará también para este ejercicio 4.5:
Algoritmo del Método de interpolación de Lagrange en C++
1 // P r o c e s o para I n t e r p o l a r Lagrange
2 float x i n t e r p o l [ cc ] , yinterpolada [ cc ] ;
3 f l o a t productosNum , productosDen ;
4 yinterpolada [ 0 ] = 0;
5
6 f o r ( i n t i =0; i <c c ; i ++){

Pág. 29
7 xinterpol [ i ] = x [ i ] ;
8 }
9
10 f o r ( i n t i i =0; i i <c c ; i i ++){

11 f o r ( i n t i =0; i <l (Tm) ; i ++){
12 productosNum =1;
13 productosDen = 1 ;
14 f o r ( i n t j =0; j <l ( Tao ) ; j ++){
15 i f ( i != j ) {
16 productosNum = productosNum ∗ ( x i n t e r p o l [ i i ]´Tm[ j ] ) ;
17 productosDen = productosDen ∗ (Tm[ i ] ´ Tm[ j ] ) ;
18 }
19 }
20 y i n t e r p o l a d a [ i i ] = y i n t e r p o l a d a [ i i ] + ( productosNum /
productosDen ) ∗Tao [ i ] ;
21 }
22 }
23 p r i n t f ( " \n" ) ;
Resultados compilados, Tiempo y los valores interpolados de Tao.
7 Cálculo de las condiciones de borde Para este ejercicio la aplicación con-

tiene las mismas pseudo codigos del ejercició N° 01, solo se agrego para calcular
condición del borde para bifurcación.
Cálculo de las condiciones de borde de bifurcación

1 else i f ( j== N) {
2 p r i n t f ( " \ t b i f u r c a c i o n [ %i , %i ] " , i , j ) ;
3 p r i n t f ( " \n" ) ;
4

´ 1 ] [ j ´1]) ;
6
7 Cn1=.5∗Qt [ i ´ 1 ] [ j +1]ća [ j +1]∗ Pt [ i ´ 1 ] [ j +1]´ c f [ j +1]∗Qt [ i ´ 1 ] [ j +1]∗

abs ( Qt [ i ´ 1 ] [ j +1]) ∗ . 5 ;
8
9 Cn2=.5∗Qt [ i ´ 1 ] [ j +1]ća [ j +1]∗ Pt [ i ´ 1 ] [ j +1]´ c f [ j +1]∗Qt [ i ´ 1 ] [ j +1]∗

abs ( Qt [ i ´ 1 ] [ j +1]) ∗ . 5 ;
10
11 Pt [ i ] [ j ]=(CpĆn1 Ćn2 ) / ( ca [ j ]+ ca [ j +1] + ca [ j +2]) ;

12 Qt2 [ i ] [ j ]=Cn1ća [ j +1]∗ Pt [ i ] [ j ] ;
13 Qt3 [ i ] [ j ]=Cn2ća [ j +2]∗ Pt [ i ] [ j ] ;
14 Qt [ i ] [ j ]= Cpća [ j ] ∗ Pt [ i ] [ j ] ;
15 }
8 Reconocimiento de la condiciones de cada Iteración.
Pág. 30
Figura 4.16: Resultados interpolados del Tao
9 Resultados Obtenidos de las Presiones en cada Nodo.
10 Resultados Obtenidos de los Caudales en cada Nodo.
Pág. 31
Pág. 32
Figura 4.18: Resultados de Presión en cada nodo
Pág. 33
Figura 4.19: Resultados de Caudal en cada tramo
Pág. 34
5 Conclusiones
1 El trabajo realizado sintetiza las principales ecuaciones que rigen el fenómeno

del golpe de ariete, orientándose a la programación en un lenguaje programa-
ble y notación del fenómeno en los reservorios por la aplicación del cierre o
apertura de las válvulas.
2 La metodología utilizada actualmente para el diseño de redes de agua potable

no proporciona resultados óptimos, ya que solo se realiza la modelación del flujo
permanente, creando un incertidumbre en el diseño del diámetro de la tubería
que no cumplan los parámetros hidráulicos (Caudal, Velocidad y Presión).
3 Se realizó una aplicación en lenguaje C++ para el cálculo de flujo transitorio

en los sistemas de tuberías que calcula la presión y caudal en cada punto de
aplicación.
35
Bibliografía
[1] M. Hanif Chaudhry. Applied Hydraulic Transients. 2014.
[2] Bruce E. Larock, Roland W. Jeppson and Gary Z. Watters. Hydraulics of Pi-
peline Systems.
[3] Mancebo Del Castillo. Teoría del golpe de Ariete y sus aplicaciones en Ingenie-
ría Hidráulica.
[4] E. Benjamin Wylie and Victor I. Streeter. Fluid Transients.
Tabla 5.1: Add caption

76.2 0 100 0 0 100
63.5 0 100 0 0 100
50.8 0 100 0 0 100 100 31.9 68.1 100
38.1 0 100 6.2 6.2 93.8 86.603 27.626 63.878 91.504
25.4 0 100 13.5 19.7 80.3 70.711 22.557 54.684 77.241
19.05 0 100 29.2 48.9 51.1 61.237 19.535 34.799 54.334
12.7 0 100 26.6 75.5 24.5 50 15.95 16.685 32.635
9.525 0 100 23.1 98.6 1.4 43.301 13.813 0.953 14.766
4.76 4 4 96 1.4 100 0 30.611 9.765 0 9.765
2.36 10.4 14.4 85.6 100 0 21.554 6.876 0 6.876
1.18 22.8 37.2 62.8 100 0 15.241 4.862 0 4.862
0.6 17.6 54.8 45.2 100 0 10.868 3.467 0 3.467
0.3 13.5 68.3 31.7 100 0 7.685 2.452 0 2.452
0.149 12.8 81.1 18.9 100 0 5.416 1.728 0 1.728
0.074 12.5 93.6 6.4 100 0 3.817 1.218 0 1.218
0 0 0
36

Trabajo 4

Cargado por

Copyright:

Formatos disponibles

Trabajo 4

Cargado por

Información del documento

Título original

Derechos de autor

Formatos disponibles

Compartir este documento

Compartir o incrustar documentos

Opciones para compartir

¿Le pareció útil este documento?

¿Este contenido es inapropiado?

Copyright:

Formatos disponibles

Trabajo 4

Cargado por

Copyright:

Formatos disponibles

Universidad Nacional

El golpe de ariete se refiere a las fluctuaciones causadas por un repentino incremento

Las ondas de presión generadas en este transitorio se propagarán por el sistema a

2.1 Objetivos General

2.2 Objetivos Específicos

2 Desarrollar un procedimiento de cálculo para el análisis del fenómeno el cual

3 Comparar los parámetros hidráulicos en condiciones iniciales (Flujo Perma-

El golpe de ariete (choque hidráulico) es el incremento momentáneo en presión, el

Figura 3.1: Variación real de la presión en la válvula, considerando pérdidas de

Figura 3.2: Variación teórica de la presión en la válvula, sin considerar pérdidas de

C : velocidad de propagación de la onda de sobrepresión

h’ : sobrepresión en la chimenea de equilibrio.

-h’ : subpresión en la chimenea de equilibrio.

hf : perdidas por fricción en la tubería

3.1 Causas del Fenómeno

Figura 3.3: Cierre brusco de la Válvula

3.2 Efectos del Fenómeno

Figura 3.4: Efectos dañinos en la tubería

3.3 Métodos de Solución

2 Método The Lax-Wendroff One-Step Method

3 Método de diferencias finitas.

4 Método de Elementos Finitos (MEF)

5 Métodos Híbridos (MH.s)

3.4 Método de las Características

Se aplican dos ecuaciones básicas de la mecánica a un elemento de fluido en una

Ecuación de Conservación de Momento Lineal:

ρ = Densidad del Fluido

3.5 Solución del Sistema

La ecuación característica, puede ser expresada en un esquema de diferencias finitas

3.6 Condiciones de Borde

3.6.1 Reservorio Aguas Arriba

QPi,1 “ Cni ` Cai Hres (3.16)

3.6.2 Válvulas aguas abajo

3.6.3 Unión de Tuberías

Esta ecuación característica positiva corresponde a la tubería aguas arriba, la cual

CPi “ QA ` Ca HA ´ Cf QA |QA (3.21)

CPi “ QA ` Cai HA ´ Cfi QA |QA (3.25)

Cni`1 “ QB ´ Cai`1 HB ´ Cfi`1 QB |QB | (3.26)

Cni`2 “ QC ´ Cai`2 HC ´ Cfi`2 QC |QC | (3.27)

CPi ´ Cni`1 ´ Cni`2

QPi “ CPi ´ Ca HPi (3.29)

QPi`1 “ Cni`1 ` Cai`1 HPi`1,1 (3.30)

QPi`2 “ Cni`2 ` Cai`2 HPi`2,2 (3.31)

La programación esta realizado en el lenguaje de programación C++, para lo cual

La programación esta realizada para i tuberías y n tramos.

Figura 4.1: Sistema de tubería

Figura 4.2: Curva de cierre de válvula

Tabla 4.1: Datos del Tiempo vs Tao

2 Datos de Tiempo de operación de la Válvula vs Toa

3 Cálculo de valores constantes Los valores Ca y Cf se calculará mediante la

7 for ( int i =0 ; i <l (L) ; i ++){

4 Resultados de los valores Constantes

El ejercicio se calculará en tres tramos como se muestra en la figura 4.3.

Figura 4.3: Sistema de tubería

Figura 4.4: Resultados de los valores constantes

5 Cálculo de Cargas Hidráulicas en cada punto

Para calcular las perdidas de carga se utilizará la ecuación de Darcy-Weisbach4.2.

1 // C a l c u l o de Carga H i d r a u l i c a en cada punto

10 f o r ( i n t i =0; i < f i l a s ; i ++){