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

Edgar Orlando Ladino Moreno, Germán Ricardo Santos Granados, Cesar Augusto García Ubaque

• español

  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.

• English

  Context: Describing mathematically the hydraulic behavior in control structures leads to nonlinear coupled differential equations, which do not have analytical solution in the majority of engineering problems. However, it is possible to obtain approximate solutions from the finite volume method (FVM). This method converts a continuous medium with infinite variables into a discrete medium with established geometries and certain boundary conditions.

  Methods: Using the Ansys Fluent 17.0 software, a 1:70 scale 2D numerical model was developed without considering the effect of stacks and abutments on the weir under permanent flow conditions of the control structure and the radial gate system of the El Quimbo dam. The hydrodynamic analysis of the structure was carried out under free surface conditions and with different hydraulic loads, for openings of 704 meters above sea level, 706.9 meters above sea level, 709.4 meters above sea level, 712 meters above sea level and 724.6 meters above sea level. Subsequently, the numerical model was validated based on the comparison between the numerical results and the theoretical data made by both Ingetec and the tests in the physical model developed by the National University of Colombia in Manizales.

  Results: For a return period of 100 years, the numerical model estimated a sheet of 0.0306 m corresponding to 2,143 m in the El Quimbo structure. With this arrangement, 80.7996 Liters per second entered the control volume and 80.7324 Liters per second left, giving rise to a variation of 0.08316%. The largest water sheet variation of the numerical model with respect to the observed sheet is 3.1% (0.258 m), which corresponds to the increasing probable maximum for an opening of gates with elevation (724.6 m) and with a discharge of 290 L/s. Similarly, the higher speeds developed by the flow in the landfill originate at the exit of the deflector, corresponding to the abscissa K0 + 5,19 m. According to the numerical model, it was determined that the impact of the jet in the channel is made with a maximum speed of 6.15 m/s.

  Conclusions: The relation between the hydraulic load with respect to the design load must not exceed 1.33, which agrees with that recommended by the USBR. Also, the cavitation condition found in the fast can be considered incipient cavitation. In landfill design, the H/HD ratio directly affects the appearance of cavitation phenomena on the structure, the development of high velocities in the fast, the deflector, and the discharge capacity of the landfill.