FLUJO GRAFUALMENTE VARIADO
Introducción
El flujo variado puede ser clasificado como rápidamente variado o gradualmente variado.
En el primer caso (rápidamente variado) la profundidad de flujo cambia abruptamente en
una distancia comparativamente corta, por ejemplo, en un resalto hidráulico. En el otro
caso, se requieren distancias mayores para que alcancen a desarrollarse los perfiles de
flujo gradualmente variado, en la ingeniería es importante su estudio, pues determinar
cuándo se tendrá una profundidad hidráulica mayor será de gran beneficio para reducción
de costos, diseños de canales, recubrimiento, si aplicamos las ecuaciones conocidas, solo
se obtendrán resultados que se aproximen a la verdad, pero será suficiente para describir
condiciones como: tipo de perfil, profundidad, profundidad normal, velocidad, entre otras
Objetivos
Entender el comportamiento del flujo gradualmente variado y la influencia de los
controles que lo generan.
Aplicar modelos matemáticos desarrollados para el cálculo de perfiles de flujo
gradualmente variado, contrastándolos con las mediciones realizadas en el
laboratorio
Analizar perfiles de flujo experimentalmente y compararlos con los resultados
teóricos
Caracterizar el tipo de perfil producido dependiendo del tipo de flujo obtenido antes
y después del mismo
Metodología y datos
En un principio, se dispone un flujo libre por un canal, lo que nos otorga lo que en el
laboratorio llamaremos “Profundidad normal” , es decir, la profundidad que se obtiene sin
ningún tipo de control, además de diligenciar las distintas cualidades geométricas del canal
tales como: ancho (B) , largo (L), elevación aguas arriba (Z1), elevación aguas abajo (Z2),
pendiente (So), y por ultimo la altura de aforo de caudales (H); dichos datos serán de vital
importancia para caracterizar nuestro tipo de flujo y propiedades geométricas de la siguiente
manera
Caudal y caudal unitario
El caudal se realiza por el método convencional de aforo de un tubo patronado, junto con
una curva de patronamiento de formula conocida 𝑄 = 17.896 ∗ 𝐻 2.416 se tiene según el dato
obtenido (H) del tubo patronado
𝑯 = 7.5𝑐𝑚
𝑸 = 17.896 ∗ 𝐻 2.416
𝑸 = 17.896 ∗ (7.5𝑐𝑚)2.416
𝑐𝑚3
𝑸 = 2327.577 𝑆
Además, como el canal es de tipo rectangular se dispone a encontrar el caudal unitario
(q), esto para facilitar cálculos posteriores de la siguiente manera
𝑄
𝒒= 𝑏
𝑐𝑚3
2327.577
𝒒= 14𝑐𝑚
𝑆
𝑐𝑚2
𝒒 = 166.255 𝑠
Profundidad critica
Este dato es fundamental, pues con el caracterizaremos el tipo de flujo inicial, es decir
correspondiendo la profundidad normal o antes de producirse el perfil gracias a el control
impuesto
3 𝑞2
𝒀𝒄 = √981 𝑐𝑚/𝑠2
3 𝑐𝑚2 2
(166.255 )
𝒀𝒄 = √ 981 𝑐𝑚/𝑠𝑠 2
𝒀𝒄 = 3.04cm
Esta profundidad critica será nuestro primer indicativo para caracterizar el tipo de flujo
presente antes de formarse el perfil pertinente, a su ves clasificarlo dependiendo de la altura
que se presenta después de producirse
Así, como Yn<Yc se presenta inicialmente un Flujo supercrítico, esta caracterización
también puede realizarse por el numero Froude, o la pendiente critica, forma que usaremos
posteriormente
Área y Área critica
Es importante conocer el área critica de la sección pues con ella y el perímetro critico se
calcularán a groso modo las pendientes críticas, las cuales determinarán el tipo de flujo
después del perfil
𝑨𝒓𝒆𝒂 = 𝐵 ∗ 𝑌 ó 𝑨𝒓𝒆𝒂𝒄𝒓𝒊𝒕𝒊𝒄𝒂 = 𝐵 ∗ 𝑌𝑐
𝑨𝒓𝒆𝒂 = 14 𝑐𝑚 ∗ 1.7 𝑐𝑚 = 23.8 𝑐𝑚2
𝑨𝒓𝒆𝒂𝒄𝒓𝒊𝒕𝒊𝒄𝒂 = 14 𝑐𝑚 ∗ 3.04 𝑐𝑚 = 42.6 𝑐𝑚2
Perímetro y Perímetro critico
Se obtiene de la suma de los lados mojados, ya sea en un comportamiento normal o critico
𝑷𝒆𝒓𝒊𝒎𝒆𝒕𝒓𝒐 = 𝐵 + 2𝑌 ó 𝑷𝒆𝒓𝒊𝒎𝒆𝒕𝒓𝒐𝒄𝒓𝒊𝒕𝒊𝒄𝒂 = 𝐵 + 2𝑌𝑐
𝑷𝒆𝒓𝒊𝒎𝒆𝒕𝒓𝒐 = 14 𝑐𝑚 + 2 (1.7 𝑐𝑚) = 17.4 𝑐𝑚
𝑷𝒆𝒓𝒊𝒎𝒆𝒕𝒓𝒐𝒄𝒓𝒊𝒕𝒊𝒄𝒂 = 14 𝑐𝑚 + 2 (3.04𝑐𝑚) = 20.086 𝑐𝑚
Radio hidráulico y radio hidráulico critico
Es la razón entre el área mojada y el perímetro mojado, tal y como se muestra a
continuación
𝑹 𝑜 ′ 𝑹𝒄
𝐴𝑟𝑒𝑎 23.8𝑐𝑚2
𝑹 = 𝑃𝑒𝑟𝑖𝑚𝑒𝑡𝑟𝑜 = 17.4𝑐𝑚
= 1.368𝑐𝑚
𝐴𝑟𝑒𝑎𝑐𝑟𝑖𝑡𝑖𝑐𝑎 42.60𝑐𝑚2
𝑹𝒄 = 𝑃𝑒𝑟𝑖𝑚𝑒𝑡𝑟𝑜𝑐𝑟𝑖𝑡𝑖𝑐𝑜
= 20.086𝑐𝑚 = 2.12𝑐𝑚
Velocidad y Velocidad critica
se determina la velocidad media teórica por medio de la división del caudal entre el área de
la siguiente manera
𝐶𝑎𝑢𝑑𝑎𝑙 𝐶𝑎𝑢𝑑𝑎𝑙
𝑽= ó 𝑽𝒄 =
𝐴𝑟𝑒𝑎 𝐴𝑟𝑒𝑎𝑐𝑟𝑖𝑡𝑖𝑐𝑎
𝑐𝑚3
2327.577 𝑐𝑚
𝑽= 23.8 𝑐𝑚2
𝑠
= 97.797 𝑠
𝑐𝑚3
2327.577 𝑐𝑚
𝑽𝒄 = 𝑠
= 54.64
42.6 𝑐𝑚2 𝑠
Coeficiente de Manning
La suposición básica en el análisis del flujo gradualmente variado es que la rapidez de
cambio de la perdida de carga con respecto a la distancia a lo largo del canal en una cierta
sección transversal, dicha rapidez de cambio se rige principalmente por la fórmula de
Manning para el caso tenemos
2 1
𝒏 = (𝑅 3 ∗ 𝑆𝑜 2 )/𝑉
2 1
1.368𝑐𝑚3 ∗0.020412
𝒏= 97.9797𝑐𝑚/𝑠
𝒏 = 0.002
Con dicho valor encontraremos la pendiente critica, la cual nos servirá para caracterizar
nuestro flujo aguas abajo
Pendiente critica
𝑉𝑐∗𝑛 2
𝑺𝒇 = (𝑅𝑐 2/3 )
𝑐𝑚 2
54.64 ∗0.002
𝑺𝒇 = ( 𝑠
2 )
2.123
𝑺𝒇 = 0.00355
cómo experimentalmente se obtuvo que la profundidad después de producido el perfil fue
𝑌 ≥ 4.8se tiene un remanso subcrítico
Pues Y>Yc>Yn correspondiente a 4.8cm>3.04cm>1.7cm y además So>Sc>0 igual a
0.02041>0.00355>0 lo que nos da el tipo de perfil S1, para caracterizar el perfil se prosigue
a realizar el método tramo a tramo y el método de integración grafica
Anexo de datos y cálculos para el desarrollo del laboratorio
DATOS INICIALES
Ancho del canal b 14 cm
Elevación Aguas Arriba Z1 84 cm
Elevación Aguas Abajo Z2 79 cm
Longitud del tramo L 245 cm
Profundidad normal del agua Yn 1.7 cm
Altura de patronamiento H 7.5 cm
Pendiente So 0.02041
Caudal Q 2327.577 cm3/s
Caudal unitario q 166.255 cm2/s
Profundidad critica Yc 3.04 cm
Área A 23.8 cm2
Perímetro P 17.4 cm
Radio hidráulico Rh 1.368 cm
Velocidad inicial Vi 97.797 cm/s
Coeficiente de Manning n 0.002
Pendiente critica Sc 0.00355
Tipo de flujo (inicial) Flujo supercrítico
TIPO DE PERFIL PRODUCIDO
S1
Métodos de calculo
Para el cálculo de perfiles de flujo gradualmente variado se utiliza la ecuación dinámica. Se
deben hacer algunas suposiciones, entre ellas:
Se consideran tramos de análisis relativamente pequeños, de tal forma que se
pueda considerar flujo uniforme y así determinar la pendiente de fricción utilizando
una ecuación de resistencia al flujo, usualmente Manning.
La pendiente del canal es pequeña, por ende, la profundidad de flujo vertical es
aproximadamente la misma profundidad perpendicular al fondo, es decir que no se
requiere corregir la profundidad de flujo por la pendiente.
El coeficiente de rugosidad es independiente del tirante hidráulico y constante en
todo el tramo en consideración.
Para conocer la variación de la profundidad del flujo gradualmente variado en relación con
la longitud del canal ya sea hacia aguas arriba o aguas abajo de la sección de control, se
emplean métodos teóricos aproximados entre los cuales los más usados son: el método
tramo a tramo y el de integración gráfica.
Método tramo a tramo
Método que se fundamenta en dividir el canal en tramos, y se fundamenta principalmente
en la conservación de la energía y las perdidas de la misma a lo largos de los distintos
tramos de la siguiente manera
𝑉12 𝑉22
𝑍1 + 𝑌1 + = 𝑍2 + 𝑌2 + + 𝑺𝒇 ∗ ∆𝒙
2𝑔 2𝑔
𝑉𝑖 2
donde se tiene Z1 − Z2 = So ∗ ∆x ; y energía 𝐸𝑖 = 𝑌𝑖 + , reemplazando por tanto tenemos
2𝑔
𝐸2−𝐸1
∆𝑥 = 𝑆𝑜−𝑆𝑓
, por tanto, tendremos que calcular la velocidad para cada corte transversal, y
por ello se necesitara el área del mismo, dato que se calcula con las alturas encontradas
por el milímetro y la base del canal de la siguiente manera
Al dividirte los tramos en 20 cm se tiene la siguiente información
DATOS
L (cm) Y (cm)
0 10
20 9.5
40 9.1
60 8.9
80 8.5
100 8.1
120 7.8
140 7.3
160 7
180 6.5
200 6.1
220 5.5
240 4.8
De esta manera, encontramos el área para cada tramo transversal, teniendo que la base
del canal es de 14 cm se tiene
𝐴𝑟𝑒𝑎 𝑇𝑅𝐴𝑀𝑂 = 𝑌𝑖 ∗ 𝐵
𝐴𝑟𝑒𝑎 𝑇𝑅𝐴𝑀𝑂 = 10𝑐𝑚 ∗ 14𝑐𝑚
𝐴𝑟𝑒𝑎 𝑇𝑅𝐴𝑀𝑂 = 140 𝑐𝑚2
Como se conocen los datos de caudal, para encontrar la velocidad media del tramo bastara
con dividir dicho caudal entre el área de la siguiente manera
𝑄
𝑉𝑖 = 𝐴𝑖
𝑐𝑚3
2327.577
𝑉𝑖 = 140 𝑐𝑚2
𝑠
𝑐𝑚
𝑉𝑖 = 16.6 𝑠
Una ves se tiene la velocidad por tramo se procede a calcular la energía del mismo con la
ecuación anteriormente mencionada
𝑉𝑖 2
𝐸𝑖 = 𝑌𝑖 + 2𝑔
𝑐𝑚 2
(16.6 )
𝐸𝑖 = 10𝑐𝑚 + 𝑠
𝑐𝑚
2∗980 2
𝑠
𝐸𝑖 = 10.14 𝑐𝑚
Ahora se desea conocer el gradiente hidráulico, para esto usaremos la siguiente formula
2
𝑄∗𝑛
𝑆𝑓 = ( 2 )
4.64∗𝐴𝑚∗𝑅𝑚3
Así que será necesario calcular el área media, el perímetro medio y con ello el radio
hidráulico medio, la profundidad media estada rada por
𝑌𝑖+𝑌𝑖+1
𝑌𝑚 = 2
10𝑐𝑚+9.5𝑐𝑚
𝑌𝑚 = 2
𝑌𝑚 = 9.75𝑐𝑚
Por tanto, el área media será
𝐴𝑚 = 𝑏 ∗ 𝑌𝑚
𝐴𝑚 = 14𝑐𝑚 ∗ 9.75𝑐𝑚
𝐴𝑚 = 136.5 𝑐𝑚2
El perímetro medio será
𝑃𝑚 = 𝑏 + 2 ∗ 𝑌𝑚
𝑃𝑚 = 14𝑐𝑚 + 2 ∗ 9.75 𝑐𝑚
𝑃𝑚 = 33.5 𝑐𝑚
Y consigo se prosigue a calcular el radio hidráulico medio de la siguiente manera
𝐴𝑚
𝑅𝑚 =
𝑃𝑚
136.5 𝑐𝑚2
𝑅𝑚 = 33.5 𝑐𝑚
𝑅𝑚 = 4.07 𝑐𝑚
Finalmente sustituyendo en la formula anterior tendremos
2
𝑄∗𝑛
𝑆𝑓 = ( 2)
4.64∗𝐴𝑚∗𝑅𝑚3
𝑐𝑚3
2
2327.577 ∗0.002
𝑆𝑓 = ( 𝑠
2 )
4.64∗136.5𝑐𝑚2 ∗(4.07𝑐𝑚)3
𝑆𝑓 = 0.00001
Una ves que se conoce el gradiente hidráulico se prosigue a calcular el valor del
desplazamiento en X, esto para describir el perfil de forma teórica
𝐸2−𝐸1
∆𝑥 = 𝑆𝑜−𝑆𝑓
10.14−9.66
∆𝑥 = 0.02041−0.00001
∆𝑥 = 23.7620556 𝑐𝑚
Como comprobación final, el acumulado de todos los ∆𝑥, es decir 𝐿 = ∑ ∆𝑥, para nuestro
caso se obtuvo una longitud teórica 240 ≈ 237.25
método de integración grafica
El método tiene como base la expresión diferencial presentada en la ecuación siguiente
𝒅𝒀 𝑆𝑜−𝑆𝑓
𝒅𝑿
= 1−𝐹𝑟2
al momento de considerar tramos se obtiene la siguiente expresión
∆𝑥 1−𝐹𝑟 2
∆𝑦
= 𝑆𝑜−𝑆𝑓
Debido a que se necesita el número de Froude y el gradiente hidráulico para calcular la relación
tendremos que realizarlo de igual forma que en el caso anterior (método tramo a tramo)
2
𝑄∗𝑛
𝑆𝑓 = ( 2 )
4.64∗𝐴∗𝑅3
𝑐𝑚3
2
2327.577 ∗0.002
𝑆𝑓 = ( 𝑠
2 )
4.64∗140𝑐𝑚2 ∗(4.12𝑐𝑚)3
𝑆𝑓 = 0.00001
Y para el numero de Froude obtendremos la siguiente expresión
𝑉𝑖
𝐹𝑟 = , expresando esto en la forma que la formula lo requiere tenemos
√2∗981
𝑉𝑖 2
𝐹𝑟 2 = 981𝑐𝑚
2∗ 2
𝑠
16.6𝑐𝑚 2
( )
𝐹𝑟 2 = 𝑠
981𝑐𝑚
2∗ 2
𝑠
𝐹𝑟 2 = 0.02817623
O expresando en términos del caudal se tiene
𝑄 2 ∗𝑏
𝐹𝑟 2 = (𝑔∗𝐴3 )
2
𝑐𝑚3
(2327.577 ) ∗14𝑐𝑚
𝑠
𝐹𝑟 2 = 𝑐𝑚
(981 2 ∗(140𝑐𝑚2 )3 )
𝑠
𝐹𝑟 2 = 0.02817623
Como las variables A y Sf son en función de la profundidad se puede expresar todo de la
siguiente manera
𝑑𝑥
𝑑𝑦
= 𝐹(𝑌)
𝑑𝑥 = 𝐹(𝑌)𝑑𝑦
E integrando a ambos lados se tiene
𝑌𝑛
𝑥 = ∫𝑦𝑜 𝐹(𝑌)𝑑𝑦
Expresándolo de otra forma obtenemos
𝒅𝒀 𝑆𝑜−𝑆𝑓
𝒅𝑿
= 1−𝐹𝑟2
Lo que a su vez es equivalente a
∆𝑥 1−𝐹𝑟 2
=
∆𝑦 𝑆𝑜−𝑆𝑓
Para nuestro caso será
∆𝑥 1−0.02817623
=
∆𝑦 0.02041−0.00001
∆𝑥
∆𝑦
= 𝐹(𝑌) = 47.63407949
Debido a que el área de cada trapecio formado para cumplir la integración grafica es
∆𝐴 = ∆𝑋 ∗ ∆𝑌
Despejando tenemos
∆𝑥
∆𝑦
= 𝐹(𝑌)
∆𝑥 = 𝐹(𝑌)∆𝑦
En caso de ser el área media del trapecio se tendrá
∆𝑥𝑚 = (𝐹(𝑌𝑖 ) + 𝐹(𝑌𝑖+1 ))∆𝑦
∆𝑥𝑚 = 𝐹(𝑌𝑚 )∆𝑦
Así tendremos últimamente
∆𝐴𝑚 = ∆𝑥𝑚 ∗ ∆𝑌
Con ∆𝑌 = 𝑌𝑖 − 𝑌𝑖+1
Ahora bien, en el caso de nuestro ejercicio se tiene
47.634+47.4063
𝐹(𝑌𝑚 ) =
2
𝐹(𝑌𝑚 ) = 47.5023
Con ∆𝑌 = 10 − 9.5 = 0.5𝑐𝑚
Por tanto
∆𝑥 = 𝐹(𝑌)∆𝑦
∆𝑥 = 𝐹(𝑌)∆𝑦
∆𝑥 = 47.5023 ∗ 0.5𝑐𝑚 = 23.76011 𝑐𝑚
Realizando la misma comprobación 𝐿 = ∑ ∆𝑥 se obtiene 240 ≈ 231.71dando como satisfactorio el
resultado
Como paso final, se anexan las tablas de los cálculos pertinentes por los dos métodos a continuación
DATOS CALCULOS PATA METODO TRAMO A TRAMO
L (cm) Y (cm) A (cm2) V (cm/s) E (cm) Ym (cm) Am (cm2) Pm(cm) Rm (cm) Sf △x (cm) L (cm)
0 10 140.0 16.6 10.14 9.75 136.5 33.50 4.07 0.00001 23.7620556 23.7620556
20 9.5 133.0 17.5 9.66 9.3 130.2 102.60 1.27 0.00004 18.9452811 42.7073367
40 9.1 127.4 18.3 9.27 9 126.0 97.00 1.30 0.00004 9.43788682 52.1452235
60 8.9 124.6 18.7 9.08 8.7 121.8 262.40 0.46 0.00015 18.9020823 71.0473058
80 8.5 119.0 19.6 8.69 8.3 116.2 18.30 6.35 0.00001 18.6377281 89.6850338
100 8.1 113.4 20.5 8.31 7.95 111.3 23.40 4.76 0.00001 13.880683 103.565717
120 7.8 109.2 21.3 8.03 7.55 105.7 15.12 6.99 0.00001 22.8985883 126.464305
140 7.3 102.2 22.8 7.56 7.15 100.1 2341.88 0.04 0.00544 18.5013546 144.96566
160 7 98.0 23.8 7.29 6.75 94.5 179.76 0.53 0.00022 22.486306 167.451966
180 6.5 91.0 25.6 6.83 6.3 88.2 15.64 5.64 0.00001 17.3958418 184.847808
200 6.1 85.4 27.3 6.48 5.8 81.2 35.40 2.29 0.00004 25.1819555 210.029763
220 5.5 77.0 30.2 5.97 5.15 72.1 27.70 2.60 0.00004 27.217129 237.246892
240 4.8 67.2 34.6 5.41 COMPARACION: 240 ≈ 237.25
DATOS CALCULOS POR EL METODO DE INTEGRACION GRAFICA
L (cm) Y (cm) A (cm2) P (cm) R (cm) Sf Fr2 F (Y) F (Y)m △Y=Yi-Yi-1 △x (cm) L (cm)
0 10 140.00 34.00 4.12 0.00001 0.02817623 47.6340749 47.5202309 0.5 23.7601154 23.7601154
20 9.5 133.00 33.00 4.03 0.00001 0.03286337 47.4063868 47.2964109 0.4 18.9185644 42.6786798
40 9.1 127.40 32.20 3.96 0.00001 0.03739028 47.186435 47.1238021 0.2 9.42476042 52.1034402
60 8.9 124.60 31.80 3.92 0.00001 0.03996804 47.0611692 46.9174938 0.4 18.7669975 70.8704378
80 8.5 119.00 31.00 3.84 0.00001 0.04588029 46.7738185 46.6003103 0.4 18.6401241 89.5105619
100 8.1 113.40 30.20 3.75 0.00001 0.05301855 46.4268021 46.2722798 0.3 13.8816839 103.392246
120 7.8 109.20 29.60 3.69 0.00001 0.05937438 46.1177575 45.8002814 0.5 22.9001407 126.292386
140 7.3 102.20 28.60 3.57 0.00001 0.0724293 45.4828052 45.2464301 0.3 13.573929 139.866316
160 7 98.00 28.00 3.50 0.00002 0.08214644 45.010055 44.5123626 0.5 22.2561813 162.122497
180 6.5 91.00 27.00 3.37 0.00002 0.10259893 44.0146703 43.4903815 0.4 17.3961526 179.518649
200 6.1 85.40 26.20 3.26 0.00002 0.12413475 42.9660927 41.8645232 0.6 25.1187139 204.637363
220 5.5 77.00 25.00 3.08 0.00003 0.16935375 40.7629537 38.6795616 0.7 27.0756931 231.713056
240 4.8 67.20 23.60 2.85 0.00004 0.25477639 36.5961695 COMPARACION: 240 ≈ 231.71
Para seguir con el desarrollo del laboratorio se prosigue a graficar los perfiles, tanto por el método
experimental como por el teórico, obteniendo
Perfiles de flujo
10
Profundidad hidraulica Y (cm)
9
4
250.000 200.000 150.000 100.000 50.000 0.000
L (cm)
Metodo Tramo a Tramo Metodo Integracion grafica Experimental
Análisis de resultados
Al comparar las curvas obtenidas por el método experimental y teórico, ve aprecia una clara
variación en el método experimental, por eso es de vital importancia realizar estos ajustes por
métodos teóricos para así normalizar el comportamiento
Conclusiones
El método de integración gráfica fue para el caso la que más cambia del método
experimental, esto se debe a que método de integración gráfica difiere de los demás
métodos en situaciones en que los tirantes de inicio y llegada varían mucho, este
fue el caso y se puede observar que no hay mucha similitud de los resultados en
comparación al método tramo a tramo, quien dio mejores resultados respecto a lo
experimental
La curva de remanso apreciada fue la de tipo S1, la cual por observación concordó
con el tipo dado en el desarrollo del laboratorio
