UNIVERSIDAD NACIONAL DEL ALTIPLANO
FACULTAD: CIENCIAS AGRARIAS
ESCUELA: TOPOGRAFIA Y AGRIMENSURA
RESUMEN
ANÁLISIS DEL FLUJO A TRAVÉS DE BIFURCACIONES EN CANALES
CURSO: HIDRAULICA II
DOCENTE: Ing. JASMANI YOVANI FLORES TEVES
ESTUDIANTE: ELVIS FRANS QUISPE MENDOZA
GRUPO: B
CODIGO: 170093
SEMESTRE: VI
PUNO_PERU
2020
INTRODUCCIÓN AL ANÁLISIS DEL FLUJO A TRAVÉS DE
BIFURCACIONES EN CANALES CON ESCURRIMIENTO A SUPERFICIE
LIBRE
RESUMEN Las estructuras hidráulicas de bifurcación, ampliamente utilizadas en
escurrimientos a superficie libre, permiten dividir el flujo que arriba a un sitio en dos o
más porciones. Este estudio comprende una revisión bibliográfica en la temática y la
aplicación correspondiente al prediseño de una estructura de división de flujo: una
bifurcación en el sistema hídrico Ludueña-Ibarlucea (Rosario, Argentina) cuyo objeto es
derivar caudales desde los cursos superficiales hacia un sistema de conductos cerrados a
construir, denominado Aliviador 3 Ramal Grandoli. Se llegó a describir el mecanismo
hidráulico que permite la derivación de la mitad del caudal transportado por el Canal
Ibarlucea a través de un canal de la mitad de su ancho y con una pendiente de fondo
mínima, de manera de cumplir una limitación de número de Froude implícita en el modelo
utilizado. Asimismo, se calcularon los coeficientes de pérdida de energía asociados. A
través de la metodología presentada se constata que es posible determinar un conjunto de
características geométricas de utilidad para el diseño preliminar, dejándose en claro que
para el diseño completo de tal estructura hidráulica resulta necesario recurrir al estudio
experimental mediante modelación física a fin de obtener una descripción más precisa del
flujo.
INTRODUCCIÓN
Una estructura de bifurcación o partición de flujo de un canal es una estructura que divide
el flujo en dos o más porciones, derivando parte del caudal que arriba a un sitio desde un
canal principal, pudiendo ser el caudal derivado incluso mayor al que continúa en la
alineación principal. La división de caudales es utilizada usualmente en sistemas de riego
y de distribución de plantas de tratamiento de agua. En estos últimos, canales
alimentadores dispuestos en grupos se enlazan a tanques de mezclado. Asimismo puede
encontrarse su aplicación en redes de desagües pluviales, tanto en grandes conductos
cerrados como en canales a cielo abierto. En obras de toma o en embocaduras de
conductos para entubamientos de cursos de agua se emplean estructuras de bifurcación
cuando por ejemplo se debe atravesar una zona densamente urbanizada. De la misma
manera, en obras de aliviadores de cursos de agua existentes se recurre a la extracción de
caudal del canal principal mediante conducciones cerradas para recuperar un determinado
nivel de protección ante inundaciones, recobrándose así capacidad de transporte de flujo.
Otra aplicación habitual de esta estructura hidráulica la constituye una boca de tormenta
en una red de drenaje pluvial ubicada en un punto intermedio. En principio, analizando
por consideraciones energéticas una estructura hidráulica de este tipo, se llega a que se
desarrolla, a lo largo del canal principal, una contracción en la bifurcación y una
expansión aguas abajo, lo que se traduce en un descenso del nivel en el canal de
derivación y un ascenso en el canal principal aguas abajo de la bifurcación. Esta
descripción general del comportamiento del flujo puede observarse en la Figura 1, que
además presenta los elementos generales de una bifurcación en un canal rectangular y en
donde pueden observarse las corrientes secundarias que aparecen en la misma
(particularmente una fuerte corriente de fondo hacia el canal de derivación), siendo δ el
ángulo de bifurcación e indicando con los subíndices 0, u y a a los parámetros hidráulicos
y geométricos (caudal Q y ancho b) en las secciones del canal principal antes y luego de
la bifurcación y del canal derivador, respectivamente. Un caso específico de este tipo de
estructura se presenta en el proyecto de saneamiento del sistema hídrico del Arroyo
Ludueña (Rosario, Santa Fe). En tal sistema, ante la necesidad de aliviar el caudal máximo
que actualmente llega al sistema de conductos existentes, se ha planteado una obra de
división del caudal transportado por el Canal Ibarlucea para derivarlo hacia un sistema de
conductos cerrados denominado Aliviador 3 Ramal Grandoli (a proyectar). En vistas de
aportar lineamientos teóricos para el diseño de esta obra hidráulica surge la necesidad de
analizar el flujo en bifurcaciones. El aporte se llevó a cabo a través de una recopilación
bibliográfica y luego mediante la aplicación de los enfoques considerados más acordes
con el tipo de estructura hidráulica planteada.
Figura 1. Bifurcación en un canal rectangular: (a) esquema general, (b) ( ) corrientes de superficie, ( ) corrientes de
fondo, (..….) línea de división de flujo, (- - - -) separación del flujo de las paredes (extraído de Hager, 2010).
BIFURCACIONES EN T.
Entendiendo bajo esta denominación a las bifurcaciones constituidas por un canal recto
el cual contiene un canal derivador en ángulo recto de ancho menor o igual que el primero,
para Fu < F0 < 0.75 (siguiendo en esta sección la misma nomenclatura de la anterior) y
sin presentarse condiciones de sumergencia en el canal derivador, Krishnappa y
Seetharamiath (1963) plantearon una razón de caudal derivado sin considerar relaciones
entre los tirantes del flujo aguas arriba y aguas abajo, según la ecuación (6), en la que βa
es la razón entre los anchos del canal de derivación y del canal principal, definida según
la ecuación (7) (Figura 5).
De modo de presentar una ecuación alternativa a (6) se recurrió a Lakshmana Rao et al.
(1968); citados en Hager (2010), que plantea una relación entre el caudal aguas abajo de
la bifurcación (Qu) y el caudal total que arriba a la misma (Q0):
Figura 5. Bifurcación en T: (a)
esquema en planta y
nomenclatura; (b) qa según F0
para Fa > 1 y βa (Krishnappa y
Seetharamiath, 1963; citado por
Hager, 2010); (c) efecto de
depresión de la superficie libre del
flujo en ingreso al canal lateral
(extraído de Hager, 2010).
Según Hager (2010), en flujos con Fa > 1/3 se desarrolla un resalto hidráulico y entonces
debe usarse Fa = 1/3 en la ecuación (8) para computar la relación de caudales Qu/Q0
teniendo en cuenta tal efecto. Por otra parte, el coeficiente de contracción Cc del canal
derivador depende de Fu (Figura 6) según:
donde Cc∞ se define como el coeficiente de contracción para Fu >> 0:
La sección de máxima contracción de flujo en el canal derivador se localiza
aproximadamente a ½ba a partir de la entrada al mismo aguas abajo. La longitud de la
separación del flujo de la pared lateral depende de Fu y Fa. Asimismo, la dirección media
del flujo en la bifurcación puede ser determinada mediante Cδ (coeficiente adimensional
que varía entre 0 y 1) en función de Fu, el cual puede obtenerse de la siguiente manera
(Hager, 2010):
Luego, el ángulo medio de bifurcación se calcula como el producto entre Cδ y δ, donde δ
= 90º para el caso de bifurcaciones en T.
Figura 6. Bifurcación en T: (a) esquema en planta y nomenclatura; (b) relación de repartición de caudales Qu/Q0; (c)
coeficiente de contracción Cc; (d) longitud relativa del despegue de la corriente en el ramal lateral La/ba; (b), (c) y (d)
como función de Fu y Fa; (….) resalto hidráulico en canal lateral (Lakshmana Rao et al., 1968; citado por Hager,
2010) (extraído de Hager, 2010)
División de flujo con derivador corto a 90º Ramamurthy y Satish (1988) propusieron un
modelo teórico simple para la división del flujo con un canal de derivación corto
emplazado perpendicularmente a un canal principal. Se consideraron varias relaciones de
ancho entre canales. El modelo fue desarrollado utilizando los principios de momento,
energía y continuidad y está basado en la existencia de condiciones de flujo libre en el
canal de derivación. El análisis hace uso de la similitud de la configuración del flujo entre
la bifurcación de un canal y un conducto con salida lateral con barrera (u orientador de
flujo) analizado bidimensionalmente para estimar el coeficiente de contracción del flujo
convergente al canal de derivación. Se relacionó la razón de descarga derivada con el
número de Froude (F) en el canal principal aguas abajo de la bifurcación. Este modelo
teórico fue validado con información de estudios experimentales efectuados
especialmente para el caso y de otros anteriores, resultando en la imposición de una
condición de flujo no sumergido en la zona de la derivación debido a condiciones de
control aguas abajo para satisfacer tal validez. La Figura 7 ilustra un esquema donde se
presentan las variables geométricas y de flujo relacionadas al problema: V, Q, h y g
designan la velocidad, el caudal, el tirante y la aceleración de la gravedad,
respectivamente; se denomina con B al ancho del canal principal y con L al ancho del
canal de derivación (debe ser L ≤ B). Los subíndices 1, 2, 3* y 3 refieren a las secciones
indicadas en el esquema. La sección 3 está localizada a una distancia cercana a 4L desde
la entrada a la derivación.
Figura 7. Flujo en un canal con bifurcación (adaptado de Ramamurthy y Satish, 1988).
Se plantea la determinación de la dependencia de Q3/Q1 con F2 (o F1) y L/B, para la cual
Ramamurthy y Satish (1988) proponen relacionar el coeficiente de contracción del flujo
Cc con el parámetro de velocidad:
Sin embargo, existe una dificultad para encontrar tal dependencia debido a las
complejidades asociadas al flujo a superficie libre, y entonces se hace uso de la similitud
de la configuración del flujo con un conducto con salida lateral con orientador de flujo
(Figura 8) que permite estimar Cc para el canal de derivación a partir de tal analogía.
Debido a la compleja naturaleza del flujo en el canal de derivación sólo pueden hallarse
valores globales de Cc, representativos de la sección entera contraída para el modelo del
canal de derivación. Luego, Q3/Q1 se relaciona con F2 (o F1) para valores dados de L/B
(Ramamurthy y Satish, 1988). Tanto resultados provenientes de estudios anteriores
(Krishnappa y Seetharamiah, 1963; citados por Hager, 2010; Law, 1965; Sridharan, 1966)
como resultados de experiencias realizadas en Ramamurthy y Satish (1988) fueron
utilizados para verificar el modelo propuesto. Las experiencias provenientes de Law
(1965) estuvieron restringidas a los casos de sección de máxima contracción en la
derivación no sumergida debido al control del flujo por una compuerta al final del canal
de derivación. Law (1965) y Sridharan (1966) indican que el flujo a la entrada de la
derivación está generalmente no sumergido cuando F3 > 0.3. En este caso se aumentó tal
limitación a un valor de F3 ≥ 0.35 (Ramamurthy y Satish, 1988).
Figura 8. Variación de Cc con η1 2 en bifurcaciones
laterales con orientadores de flujo (adaptado de
Ramamurthy y Satish, 1988).
La configuración general del flujo entrante al canal de derivación se aproxima al patrón
de flujo lateral que pasa por un conducto con salida lateral con barrera analizado
bidimensionalmente, en el cual L representa el ancho de la salida y B el ancho del
conducto (Figura 8). En ambas configuraciones el flujo de aproximación a la derivación
se separa del canal principal y comienza a contraerse hasta llegar a un ancho CcL en la
sección de la máxima contracción.
Esta similitud entre las dos configuraciones de flujo en cuanto a la contracción del flujo
bifurcado permite estimar el coeficiente de contracción Cc para el canal de derivación
utilizando el modelo del conducto bidimensional con salida lateral con orientador de flujo
(Ramamurthy y Satish, 1988). Con el objeto de desarrollar el modelo teórico del flujo en
el canal de derivación se formula un conjunto de hipótesis: los canales principal y de
derivación son rectangulares y horizontales, L ≤ B y el flujo en el canal principal es
subcrítico; la pérdida de energía entre las secciones 1 y 3* es despreciable; la distribución
de presiones es hidrostática en las secciones 1 y 2 y la presión media actuando a lo largo
de la superficie divisoria MN es el promedio de las alturas piezométricas en M y N; la
pérdida por fricción es despreciable entre 1 y 2; cuando el flujo no está sumergido a la
entrada de la derivación, las condiciones de flujo son aproximadamente críticas en la
sección 3* donde ocurre la máxima contracción; Cc para el modelo de conducto de salida
lateral con barrera analizado bidimensionalmente provee una buena estimación para el
coeficiente de contracción global en el canal de derivación cuando los valores de L/B y
η1 son los mismos en los dos casos (Ramamurthy y Satish, 1988).
Aplicando conservación de la energía entre las secciones 1 y 3*, definiendo el parámetro
de velocidad η1 = V1/V3*, aplicando continuidad entre las secciones y conservación de
momento en la dirección del flujo principal y operando matemáticamente, Ramamurthy
y Satish (1988) llegan a un conjunto de ecuaciones que permiten resolver el modelo. La
solución es descrita a continuación:
1. dado F1, se obtiene η1 de la siguiente ecuación:
2. de la solución existente del problema del conducto analizado bidimensionalmente
se obtiene Cc entrando con η1 y L/B; con L/B, η1 y Cc se calcula Q3/Q1 para el
canal de derivación empleando la siguiente ecuación:
3. de la ecuación mostrada a continuación se puede obtener la relación de tirantes
antes y después de la derivación en el canal principal (h2/h1):
4. el número de Froude aguas abajo F2 se obtiene de la siguiente ecuación:
5. por último, dada una relación h-Q del canal de derivación se debe comprobar que
F3 > 0.35.
Si la desigualdad F3 > 0.35 no se cumple puede observarse de los resultados
experimentales que las condiciones del flujo en la sección de máxima contracción
se tornan claramente subcríticas y por lo tanto cualquier cambio en la condición
del flujo en el canal de derivación viaja aguas arriba y afecta drásticamente el flujo
en el canal principal. Esto puede observarse en la Figura 9, donde se ilustra la
relación de tirantes h2/hc (siendo hc es el tirante crítico) para distintos F3. La
gráfica se construyó con el modelo teórico proporcionado anteriormente y los
resultados experimentales se graficaron como puntos, observándose una adecuada
aproximación entre valores computados y observados. La condición de flujo no
sumergido en la derivación en la zona de la bifurcación se encuentra normalmente
en canales con derivaciones cortas o con derivaciones en pendiente. Más
generalmente, las condiciones del flujo en la derivación pueden depender de la
naturaleza del control en la sección aguas abajo de la bifurcación. Cuando el flujo
en la sección de máxima contracción (sección 3*) es crítico, el flujo
inmediatamente aguas abajo de esta sección se torna supercrítico y expansivo.
Dependiendo del control aguas abajo el flujo puede permanecer supercrítico o
formarse un resalto hidráulico. Es por esto entonces que cuando el flujo sea crítico
en la sección de la máxima contracción, los cambios en el canal de derivación no
afectarán significativamente el flujo en el canal principal. Asimismo resulta
oportuno explicar que según el modelo teórico, para valores fijos de L/B, Cc
aumenta a medida que η1 disminuye (Ramamurthy y Satish, 1988).
En la Figura 10 se contrastaron los resultados teóricos de F2 en función de Q3/Q1
para distintas L/B con los resultados experimentales y se halló que coinciden para
valores de F2< 0.7 (Ramamurthy y Satish, 1988).
Por último, se expone un diagrama de flujo explicando un método iterativo de
cálculo del caudal derivado basado en las gráficas de F2 en función de Q3/Q1 para
distintas L/B de la Figura 10
Figura 9. Variación de h2/hc con F3 (extraído de Ramamurthy y Satish, 1988).
Figura 10. Variación de Q3/Q1 con F2 para diferentes valores de L/B; resultados de estudios previos
citados por Ramamurthy y Satish, 1988; ○ resultados de las experiencias presentadas por Ramamurthy y
Satish, 1988 (adaptado de Ramamurthy y Satish, 1988).
