Simulación numérica CFD de la estructura de control y del sistema de compuertas radiales - represa El Quimbo

Abstract

Contexto: Describir matemáticamente el comportamiento hidráulico en estructuras de control conlleva a ecuaciones diferenciales acopladas no lineales, las cuales no cuentan con solución analítica en la mayoría de los problemas de ingeniería. No obstante, es posible obtener soluciones aproximadas a partir del método de volúmenes finitos (FVM). Este método convierte un medio continuo con variables infinitas en un medio discreto con geometrías establecidas y condiciones de contorno determinadas.
 Método: Se desarrolló un modelo numérico en 2D sin considerar el efecto de las pilas y los estribos en el azud, a escala 1:70, bajo condiciones de flujo permanente de la estructura de control y del sistema de compuertas radiales de la represa El Quimbo, a partir de la implementación el software Ansys Fluent 17.0. El análisis hidrodinámico de la estructura se realizó bajo condiciones de superficie libre y con diferentes cargas hidráulicas, para aperturas de 704 msnm, 706,9 msnm, 709,4 msnm, 712 msnm y 724,6 msnm Posteriormente, se validó el modelo numérico a partir de la comparación de los resultados numéricos con los datos teóricos realizados por Ingetec y los datos de las pruebas realizadas en el modelo físico desarrollado por la Universidad Nacional de Colombia (sede Manizales).
 Resultados: Para un periodo de retorno de 100 años, el modelo numérico estimó una lámina de 0,0306 m que corresponden a 2,143 m en la estructura El Quimbo. Con esta disposición, ingresó al volumen de control un total de 80,7996 L/s y salieron 80,7324 L/s, originando una variación de caudal de 0,08316 %. La mayor variación de lámina de agua del modelo numérico con respecto a la lámina observada es de 3,1 % (0,258 m), la cual corresponde a la creciente máxima probable para una apertura de compuertas con cota (724,6 msnm) y con una descarga de 290 L/s. De igual forma, las mayores velocidades desarrolladas por el flujo en el vertedero se originan en la salida del deflector, correspondiente a la abscisa K0+ 5,19 m. Según el modelo numérico, se determinó que el impacto del chorro en el canal se hace con una velocidad máxima de 6,15 m/s.
 Conclusiones: La relación entre la carga hidráulica con respecto a la carga de diseño no debe superar 1,33, lo cual a su vez concuerda con lo recomendado por la USBR. Asimismo, la condición de cavitación encontrada en la rápida puede ser considerada cavitación incipiente. En el diseño de vertederos, la relación H/HD incide de manera directa en la aparición de fenómenos de cavitación sobre la estructura, en el desarrollo de altas velocidades en la rápida, en el deflector y en la capacidad de descarga del vertedero.