UT7
UT7
ESCURRIMIENTO SUPERFICIAL 1
UT7 de
FLUJO EN CONDUCCIONES
Rev. A 45
MODELOS MATEMÁTICOS
CATEDRA DE HIDROLOGIA
Y
OBRAS HIDRAULICAS
El análisis de crecidas en una cuenca lleva implícito dos fenómenos que se dan en la
misma, escurrimiento ó flujo superficial y flujo en canales apreciables. Por lo tanto es preciso
conocer los procesos teóricos del fluido escurriendo tanto en una superficie como en
conducciones.
En el presente apunte se incluirá un resumen de ambos fenómenos con la finalidad de darle
mayor comprensión a estos procesos.
Cabe aclarar la conveniencia de conocer la teoría sobre los procesos hidrológicos basados
ampliamente en el TEOREMA DE TRANSPORTE DE REYNOLDS ampliamente detallado
en el capítulo 2 del libro “Hidrologia Aplicada – (Chow-Maidment-Mays / 1994). También
llamado Ecuación general del volumen de control. Este teorema se utiliza para desarrollar
las ecuaciones de continuidad, momentum y energía de los procesos hidrológicos.
Ec. 1
Donde:
B: Propiedad extensiva (valores que dependen de la cantidad de masa presente)
: Propiedad intensiva (cantidad de B por unidad de masa presente)
V: Velocidad del fluido
: Densidad del fluido
: Diferencial de volumen
dA : Diferencial de área
Siendo :
Este teorema aplica a múltiples casos de procesos hidrológicos dado que tratan leyes físicas
que se aplican a masas discretas de una sustancia o bien aplicadas a un fluido que circula
continuamente a través de un volumen de control.
Volumen almacenado
Ec.2
Ec. 3
En este caso se dice que el sistema es CERRADO, caso contrario será ABIERTO. Tal es así
que el ciclo hidrológico es un sistema cerrado para el agua, pero el proceso LLUVIA-
ESCORRENTÍA en una cuenca conforma un sistema abierto porque no toda la lluvia se
convierte en escorrentía, una parte de ella vuelve a la atmosfera mediante la evaporación y
otra parte se transforma en una escorrentía subsuperficial, aunque finalmente el agua
llegará al final de la cuenca, esta lo hará en un plazo tardío debido a las características
geomorfológicas de la cuenca.
PROCESO DE LLUVIA-ESCORRENTÍA
COMO PROCESO HIDROLÓGICO EN UNA CUENCA
Donde k es una cte. Que tiene dimensiones de tiempo. Por consiguiente, la tasa de cambio
del almacenamiento con respecto al tiempo dS/dt es igual a la diferencia entre la entrada y
la salida
Luego:
Donde D es el operador diferencial d/dt. Si la ecuación de transformación ha sido
determinada y puede ser resuelta, se encuentra la salida como función de la entrada.
La ecuación anterior describe un sistema lineal si k es constante. Si k es función de la
entrada o la salida, entonces esta ecuación describirá un sistema no lineal cuya solución es
más compleja.
Realizado por el Ing. DIEGO TAUS 2020 UT Nº7
REV. Colaborador Ing DAVID HARDOY AÑO
Archivo: [Link]
UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
ESCURRIMIENTO SUPERFICIAL 6
UT7 de
FLUJO EN CONDUCCIONES
Rev. A 45
MODELOS MATEMÁTICOS
La función f, está determinada por la naturaleza del sistema hidrológico que se está
estudiando. Por ejemplo, el modelo del embase lineal para flujo base en corrientes,
relaciona el almacenamiento y el flujo de salida según la siguiente expresión:
Donde a1, a2, an,….b1, b2, bm son constantes y las derivadas de orden superior a aquellas
mostradas se desprecian.
Realizado por el Ing. DIEGO TAUS 2020 UT Nº7
REV. Colaborador Ing DAVID HARDOY AÑO
Archivo: [Link]
UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
ESCURRIMIENTO SUPERFICIAL 7
UT7 de
FLUJO EN CONDUCCIONES
Rev. A 45
MODELOS MATEMÁTICOS
Esta ecuación fue presentada por Chow y Kulandaiswamy (1971) como el MODELO
GENERAL DEL SISTEMA HIDROLÓGICO (MGSH). La que podemos expresar en forma
mas compacta como:
A su vez, considerando la ecuación del modelo general del sistema hidrológico aplicado a un
embalse lineal de constante a1 = k, podemos expresar la misma, de la siguiente forma:
FLUJO SUPERFICIAL
Este proceso ocurre en la parte superior de las superficies antes de concentrarse en canales
reconocibles. Este proceso tiene lugar cuando transcurre tiempo suficientes desde el inicio
de la precipitación. Dos condiciones que se tienen en cuenta en al análisis es que el flujo es
permanente y la densidad del fluido es constante. Por lo tanto podemos escribir la ecuación
de continuidad de la siguiente manera:
De la siguiente imagen:
Figura 1
Donde :
i : Intensidad de precipitación
f : Velocidad de infiltración
So : Pendiente
: angulo de la pendiente
y : Tirante
v : Velocidad del flujo
f . Lo . cos + v . y – i . Lo . cos = 0
De la ecuación de continuidad anterior podemos obtener el caudal por unidad de ancho (qo)
Según Roberson y Crowe (1985) llegaron a la conclusión que la velocidad promedio para
un flujo laminar en un plano inclinado viene dada por la siguiente expresión:
Esta expresión es similar a la propuesta por DARCY para describir la resistencia de un flujo
R: (Radio hidráulico)
Donde corresponde a la fricción = 96 / Re donde Re = [Link] /
Con R = y, Sf = So para flujo uniforme y qo = v.y esta ecuación puede resolverse para y con
la siguiente expresión:
FLUJO EN CANALES
Ec. 4
Luego, se puede ver que para calcular la profundidad y la velocidad en varios puntos a lo
largo del canal es necesaria una solución iterativa de la ecuación de Manning:
Ec. 5
METODO DE NEWTON
El método de Newton puede aplicarse iterativamente para dar una solución numérica.
Partiendo de la base de la selección de la profundidad yi en la iteración j lo que permite
calcular el flujo Qj con la ecuación 5 utilizando el área y el radio hidráulico correspondientes
a yi. Este nuevo Qj, se compara con el flujo real Q (calculado con la ecuación 4), esto
permitirá obtener un nuevo yi que permita calcular un nuevo Qj iterando tantas veces como
sea necesario para reducir el error entre dos iteraciones consecutivas.
Donde el subíndice j por fuera de los paréntesis indica que sus contenidos se evalúan para y
= yi. La expresión para gradiente es útil en el método de Newton, donde, dada una
alternativa de yi, yi+1, se escoge para satisfacer
Este valor yi+1 es el valor de y en una gráfica de f vs y, donde la tangente a la curva en y=yi
interseca al eje horizontal, tal como se ilustra en la figura 2. Resolviendo entonces la
ecuación anterior para yi+1:
Ec. 6
Figura 2
Se adjunta una tabla donde dan algunas adaptaciones de este para las figuras de canales
mas conocidas.
El análisis del flujo, tanto en tuberías (conductos cerrados) como en canales abiertos, se
basa en la aplicación de tres ecuaciones básicas de la mecánica de los fluidos:
Continuidad
Energía
Por Bernoulli:
Cantidad de movimiento
A diferencia de las ecuaciones de continuidad y de la energía que son relaciones escalares
la ecuación de la cantidad de movimiento es una relación vectorial, es decir en las que
intervienen, tanto la magnitud como la dirección de las fuerzas y velocidades.
La ley de conservación de la cantidad movimiento puede enunciarse del siguiente modo: La
variación en el tiempo de la cantidad de movimiento (que se define como el caudal másico
.AV multiplicado por la velocidad V) a lo largo de un tubo de corriente ocasiona una fuerza
llamada fuerza de impulso. La aplicación de la ecuación de la cantidad de movimiento al
análisis de flujo en el tramo de tubería que se muestra es como sigue.
La cantidad de movimiento que atraviesa el área A1 durante el tiempo dt es igual a:
(2A2 . V2 dt) . V2
La fuerza neta de impulso del fluido causada por la variación de la cantidad de movimiento
entre las secciones 1 y 2 es:
Flujo laminar
En el flujo laminar las partículas de fluido no se mezclan más que a escala molecular de
modo que durante el movimiento dichas partículas se desplazan según trayectorias
paralelas bajo la acción de la viscosidad. En la práctica, el flujo laminar se produce cuando
el numero de Reynolds no excede de 1500 a 2000.
Velocidad crítica
La velocidad crítica de interés práctico para el ingeniero es aquella velocidad por debajo de
la cual toda turbulencia es amortiguada por la acción de la viscosidad del fluido. La
experiencia demuestra que un límite superior para el régimen laminar, en tuberías, viene
fijado por un valor del numero de Reynold alrededor de 2000, en la mayoría de los casos
prácticos.
Flujo turbulento
En el flujo turbulento las partículas del fluido se mezclan a escala molar, de modo que
durante el movimiento se produce un intercambio de cantidad de movimiento entre
partículas adyacentes , ocasionando una rápida y continua agitación y mezcla en el seno del
fluido. En la práctica, el flujo turbulento se produce para números de Reynolds por encima
de valores entre 6000 a lo 10.000
Pérdida de carga
Pérdida de carga es la pérdida de energía que experimentan los líquidos que fluyen en
tuberías y canales abiertos. La energía necesaria para vencer los efectos del rozamiento en
el flujo turbulento es pérdida de carga. Las pérdidas de energía localizadas en las
turbulencias inducidas por las piezas especiales y accesorios que se utilizan en tuberías y
canales , son también pérdidas de carga. La pérdida de carga se representa usualmente
con el símbolo hL.
Línea piezométrica
La línea piezométrica, como muestra también la figura, es la línea que une los puntos hasta
los que el liquido podría ascender si se insertasen tubos (denominados tubos piezometricos)
en distintos lugares a lo largo de la tubería o canal abierto. Es una medida de la altura de
presión hidrostática disponible en dichos puntos. Cuándo el agua fluye en un canal abierto,
la línea piezometrica coincide con la superficie libre del liquido.
Energía especifica
La energía especifica E también llamada carga especifica es la suma de la altura
piezométrica y la altura cinética o de presión dinámica U2/2g. El concepto de energía
especifica es especialmente útil en el análisis de flujo en canales abiertos. La relación entre
tirante y descarga por unidad de ancho del canal (caudal unitario), para energía especifica
constante, se conoce como curva q.
Flujo permanente
El flujo permanente se produce cuando la descarga o caudal en cualquier seccion
transversal permanece constante.
Flujo variado
El flujo en un canal se considera variado cuando e tirante cambia a lo largo del canal. En
general el flujo puede ser gradualmente variado (FGV) o rapidamente variado (FRV), como
se muestra en la figura el flujo rápidamente variado tiene lugar cuando él tirante varia
bruscamente.
NUMERO DE REYNOLDS
Para tuberías circulares a sección llena, el numero de Reynolds viene dado por la siguiente
expresión:
Donde:
V = Velocidad media del fluido en m/s
d = Radio de la tuvería en m
= Viscosidad cinemática del fluido en m2/seg.
= Densidad del fluido en UTM/m3 ó kg.seg2/m4.
= Viscosidad absoluta del fluido en [Link]/m2.
(1)
(2)
Luego por Chezy y Maning se tiene:
1
1
U c R i Sabemos que c Rh 6
n
2 1
1
Entonces U Rh 3 i 2 (3) Donde Rh
n Pm
a su vez
(4)
Trabajando con la ecuación (3) podemos llegar a ecuación para hallar la velocidad en
función del tirante, el caudal, la pendiente y el coeficiente de Manning.
2 1
1
U Rh 3 i 2
n
Llamando S = Rh /
2
podemos obtener los siguientes valores en función de
distintos tirantes:
Relacion Radio Area m S S1/4
h/Ø hidáulico Ø
cm cm2
Podemos decir
(5)
(6)
Podemos decir que para que se cumpla la velocidad de autolimpieza, la fuerza tractiva debe
estar entre:
(7)
DIÁMETRO DE CONDUCTOS
1
1
U c R i Sabemos que c Rh 6
n
Entonces
2 1
1
U Rh 3 i 2 Donde el Rh
n Pm
0.1
0.2
0.3
….
……
0.94 X XX
Qn Qn
Awp 8
Awp ' 8
Ø i
3
h i
3
8
0.3352 Ø i
3
Luego Q 1
n
Caño de Hº n = 0.015
Pendiente i = 1‰
Para poner en numero un numero enter en la fórmula se realiza el siguiente análisis en el
que 1‰ = 1/1000 → i / 1000 = 0.0316 . i
Entonces de 1, tenemos:
8
Q 0.706 Ø i
3
O lo que es lo mismo:
Podemos desarrollar una formulación genérica que nos permita obtener el diámetro teórico
de un conducto sin definir previamente n y/o i
Donde:
Todos y cada uno de los softwares disponibles para el modelado de flujos de agua tanto
superficial como canalizada, disponen de distintos procesos capaces de resolver la ecuación
de continuidad tal es así que se pueden tener 3 procesos de análisis, ya sean:
- Flujo Uniforme.
- Flujo según Onda Cinemática.
- Flujo según Onda Dinámica.
Estros procesos suponen la utilización de las ecuaciones de continuidad y de cantidad de
movimiento. De manera simplificada se comentan cada uno de ellos teniendo en cuenta la
siguiente imagen del flujo del agua en superficie y canalizacion:
FLUJO UNIFORME
En este tipo de proceso representa la forma mas simple de simular el comportamiento del
agua en superficies o en el interior de los conductos mediante la expresión:
Para ello se asume que en cada uno de los incrementos de tiempo de cálculo considerados,
el flujo es uniforme. De esta forma el modelo simplemente “traslada” los hidrogramas en el
nudo aguas arriba del conducto hacia el nudo final aguas abajo del mismo, con un cierto
retardo y cambio en el aspecto del mismo. Para ello se utilizarán las características de las
conducciones como ser, coeficiente de Manning, pendiente, dimensiones, etc.
Este caso requiere que la simulación se realice a cañería llena cumpliendo con los
parámetros de Manning. Por tal razón no arroja buenos resultados si se tiene un sistema de
cuencas encadenadas donde la descarga de cuencas aguas arriba puede permanecer
estancada conforme se vayan descargando las cuencas aguas abajo.
La ventaja del proceso es que permite que tanto el caudal como el área mojada, varíen tanto
espacial como temporalmente en el interior de los conductos, esto genera una cierta
amortiguación y/o retraso de los hidrogramas de salida respecto a los de entrada.
No obstante este tipo de modelación no puede tener en cuenta; el almacenamiento de agua
que se produce en los conductos, los fenómenos de resalto hidráulico, las perdidas en las
entradas o salidas de los pozos de registro, el flujo inverso ni el flujo presurizado.
Es útil en los casos con incrementos de tiempos largos de 5 a 15 minutos.
Donde:
: Componente de la velocidad de descarga lateral qo
g : Gravedad
Io , If : Pendiente del fondo y de fricción respectivamente
Las inundaciones ocurren en el sistema cuando la profundidad (calado) del agua en los
nudos excede el valor máximo disponible en las mismas, este exceso de caudal, bien puede
perderse o bien puede generar un estancamiento en la parte superior del nudo y entrar
nuevamente al sistema posteriormente.
Se tiene los sucesivos embales como modelo de las cuencas para la modelación de
crecidas según (HUIN). Se adjuntan entonces la expresión de la función de la función para
la obtención del hidrograma de descarga
Chutha y Dooge (1990) propusieron una interpretación sencilla del HUI de Rodríguez-Iturbe
y Valdés (1979) que describía los caminos o recorridos del agua como una cascada de
depósitos lineales distintos, añadiendo al final un depósito. La representación de la cuenca
de Sorogain según este modelo se presenta en la Figura anterior.
donde aCD y bCD son las fracciones de la cuenca cuyos recorridos del agua generada en
ellas atraviesan tres o dos depósitos respectivamente, y kCD1 y kCD2 son los parámetros del
modelo. La fracción de la cuenca que pasa a través de los depósitos lineales de primer
orden, I y II en la Figura anterior, es AI,II=0.527; y para el resto, III, es
AIII=0.473.
CONCLUSION
Uno de los HUI mas representativos corresponde al HUIGD dado que utiliza el criterio de
Nash y Clark
EJEMPLO DE APLICACIÓN
En el siguiente ejemplo se realiza una modelación para el impacto de una lluvia sobre una
cuenca. En la que se dimensionarán los conductos utilizando el método SCS y para la
obtención del hidrograma total se utilizarán todos los métodos de hidrogramas vistos en la
cursada.
SCA
19 has
CN 78
Ti = 10'
SCB
28.5 has
CN 89
Ti = 10'
SCC
31 has
CN 80
Ti = 15'
Depuración s/SCS
Precipitación total 50 mm
CN ponderado 83
Ia 3 mm
S 53
Pe 22 mm
Tormenta I1 33
(t = 10min) I2 I2 53
I3 26
I1 I4 20
I3
5.5 8.8 4.4 I4
3.3 P tot 22.0
10 20 30 40
Dt HU P1 P2 P3 P4 P5 P6 P7 P8 Q total
0 0.00 0 0.00
5 0.05 0.14 0 0.14
10 0.11 0.29 0.14 0 0.43
15 0.16 0.43 0.29 0.23 0 0.96
20 0.12 0.33 0.43 0.46 0.23 0 1.46
25 0.08 0.22 0.33 0.70 0.46 0.12 0 1.82
30 0.04 0.11 0.22 0.52 0.70 0.23 0.12 0 1.89
35 0.00 0.00 0.11 0.35 0.52 0.35 0.23 0.09 0 1.64
40 0.00 0.17 0.35 0.26 0.35 0.17 0.09 1.39
45 0.00 0.17 0.17 0.26 0.26 0.17 1.04
50 0.00 0.09 0.17 0.20 0.26 0.72
55 0.00 0.09 0.13 0.20 0.41
60 0.00 0.07 0.13 0.20
65 0.00 0.07 0.07
70 0.00 0.00
Dt HU P1 P2 P3 P4 P5 P6 P7 P8 Q total
0 0.00 0 0.00
5 0.05 0.15 0 0.15
10 0.11 0.299 0.15 0 0.45
15 0.16 0.449 0.299 0.24 0 0.99
20 0.22 0.599 0.449 0.479 0.24 0 1.77
25 0.27 0.748 0.599 0.719 0.479 0.12 0 2.66
30 0.33 0.898 0.748 0.958 0.719 0.24 0.12 0 3.68
35 0.29 0.798 0.898 1.198 0.958 0.359 0.24 0.09 0 4.54
40 0.25 0.699 0.798 1.437 1.198 0.479 0.359 0.18 0.0898 5.24
45 0.22 0.599 0.699 1.277 1.437 0.599 0.479 0.269 0.1796 5.54
continúa
Realizado por el Ing. DIEGO TAUS 2020 UT Nº7
REV. Colaborador Ing DAVID HARDOY AÑO
Archivo: [Link]
UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
ESCURRIMIENTO SUPERFICIAL 37
UT7 de
FLUJO EN CONDUCCIONES
Rev. A 45
MODELOS MATEMÁTICOS
50 0.18 0.499 0.599 1.118 1.277 0.719 0.599 0.359 0.2695 5.44
55 0.15 0.399 0.499 0.958 1.118 0.639 0.719 0.449 0.3593 5.14
60 0.11 0.299 0.399 0.798 0.958 0.559 0.639 0.539 0.4491 4.64
65 0.07 0.2 0.299 0.639 0.798 0.479 0.559 0.479 0.5389 3.99
70 0.04 0.1 0.2 0.479 0.639 0.399 0.479 0.419 0.479 3.19
75 0.00 0 0.1 0.319 0.479 0.319 0.399 0.359 0.4192 2.40
80 0 0.16 0.319 0.24 0.319 0.299 0.3593 1.70
85 0 0.16 0.16 0.24 0.24 0.2994 1.10
90 0 0.08 0.16 0.18 0.2395 0.66
95 0 0.08 0.12 0.1796 0.38
100 0 0.06 0.1198 0.18
105 0 0.0599 0.06
110 0 0.00
Para el análisis se tuvo en cuenta la siguiente repartición de superficies para cada nudo:
Area Area
SCA1 0 19.0 40.0 %
SCA2 227 60.0 %
SCB1 354 28.5 20.0 %
SCB2 722 60.0 %
SCB3 991 20.0 %
SCC1 991 31.0 40.0 %
SCC2 1401 60.0 %
1
1
U c R i Sabemos que c Rh 6
n
Entonces
2 1
1
U Rh 3 i 2 Donde el Rh
n Pm
0.1
0.2
0.3
….
……
0.94 X XX
Qn Qn
Awp 8
Awp ' 8
Ø i
3
h3 i
8
0.3352 Ø 3 i
Luego Q 1
n
Caño de Hº n = 0.015
Pendiente i = 1‰ y para poner en numero un numero enter en la fórmula se realiza
el siguiente análisis en el que 1‰ = 1/1000 → i / 1000 =
0.0316 . i
8
Entonces de 1, tenemos: Q 0.706 Ø 3 i
3
Q 1 8
Depejando el Ø de la anterior, tendremos: Ø
i 0.706
0.375
Q
O lo que es lo mismo: Ø 1.14
i
Dimensionado
Areas t total I P tp Pe Q tr ØT ØC t. tr. V tr
has min mm/h mm min mm m3/seg m m min m/s
7.6 10 163
19.0 15 133 34 29 11.5 0.9 1.05 1.0 5 0.5
24.7 21 113 39 32 14.7 1.4 1.22 1.2 4 0.6
41.8 25 102 42 35 16.5 2.5 1.52 1.5 10 0.7
47.5 35 82 48 41 20.4 3.0 1.62 1.6 7 0.7
59.9 41 73 50 45 22.3 3.7 1.77 1.8 0 0.8
78.5 41 73 50 45 22.3 4.9 1.96 2.0 8 0.9
Vprom 0.71
Calcular la profundidad de Flujo de una canalización circular con los siguiente datos
So 0.001 m/m
Manning n : 0.013
En iteración 1, yj = 1.50 m h/Ø = 1.00
Función de forma del canal :
El método de Clark arrojo un valor de 6.40 m3/s contra los aproximadamente 5.54 m3/s
calculado manualmente según SCS. No obstante cabe aclarar que el método es
conservador
El valor del caudal máximo obtenido es de 6.00 m3/seg contra los aproximadamente 5.54
m3/s calculado manualmente según SCS.
Realizado por el Ing. DIEGO TAUS 2020 UT Nº7
REV. Colaborador Ing DAVID HARDOY AÑO
Archivo: [Link]
UNIVERSIDAD TECNOLÓGICA NACIONAL Pág.
ESCURRIMIENTO SUPERFICIAL 45
UT7 de
FLUJO EN CONDUCCIONES
Rev. A 45
MODELOS MATEMÁTICOS