Se llevó a cabo la simulación de la ruptura del borde de un tanque de enfriamiento
usando dos métodos: Uno de volúmenes finitos denominado CARPA y un algoritmo en diferencias finitas centradas. El algoritmo en volúmenes finitos CARPA utilizado con el pre y post procesador GiD mostró ser una herramienta muy poderosa en la animación de planicies de inundación, útiles en la interpretación de resultados tanto en el espacio como en el tiempo, para la protección civil, así como para definir las
posibles zonas afectadas debido a fenómenos como el de la ruptura de bordos.Postprint (published version

Aragón Hernández, José Luis

Arganis Juarez, Maritza Liliana

Bladé i Castellet, Ernest

Cruz Gerón, Juan Ansberto

De Luna Cruz, Faustino

Fuentes Mariles, Oscar Arturo

Mendoza Facundo, José Enedino

Sánchez Juny, Martí

Spanish

UPCommons. Portal del coneixement obert de la UPC

 Simulación de la Onda de Aavenida  por Ruptura de Tanque de Enfriamiento  Fuentes1, O.A., Arganis, M.L.1, Bladé, E. 2, Sánchez, M.2,   De Luna, F.1, Cruz, J.A.1, Mendoza, J. E.1, Aragón, J.L.1  1Instituto de Ingeniería, Universidad Nacional Autónoma de México, 2 E.T.S. d’Eginyers de Camins Canals I Ports, Campus Nord- Barcelona, España.    Resumen   Se llevó a cabo la simulación de la ruptura del borde de un tanque de enfriamiento usando dos métodos: Uno de volúmenes finitos denominado CARPA y un algoritmo en diferencias finitas centradas. El algoritmo en volúmenes finitos CARPA utilizado con el pre y post procesador GiD mostró ser una herramienta muy poderosa en la animación de planicies de inundación, útiles en la interpretación de resultados tanto en el espacio como en el tiempo, para la protección civil, así como para definir las posibles zonas afectadas debido a fenómenos como el de la ruptura de bordos.   Frases Clave  Ruptura de un dique, CARPA, Gid, volúmenes finitos, diferencias finitas.     Introducción  Cuando el bordo de un embalse se rompe o falla, se libera un gran volumen de agua en muy poco tiempo y viaja aguas abajo a alta velocidad lo cual puede causar graves pérdidas humanas y daños materiales. La estimación de  las zonas inundadas y sus profundidades, así como las velocidades del flujo tiene importantes aplicaciones porque permite conocer los efectos sobre las poblaciones existentes aguas abajo, edificaciones así como en las áreas agrícolas e industriales; o incluso permite formular medidas de mitigación de daños. Para obtener el hidrograma de salida debido a la falla es importante llevar a cabo pruebas de laboratorio; para simular los efectos aguas abajo de dicho hidrograma se requieren modelos matemáticos que deben ser capaces de representar lo mejor posible el fenómeno físico.  Aureli y Mignosa, (2001) hicieron un modelo de laboratorio de un dique para analizar la forma de la brecha y el comportamiento del flujo bidimensional en la superficie del agua aguas abajo del dique; posteriormente ellos simularon las condiciones experimentales, usando el método de captura de choque de Mc Cormack, de tipo explícito y predictor-corrector, con una precisión de segundo orden en tiempo y espacio, el cual resuelve las ecuaciones de flujo bidimensional para canales y que incluye un término suavizador de las oscilaciones numéricas que presentan comúnmente las diferencias finitas de segundo orden para los choques y para números de Fraude cercanos a uno. Finalmente los autores hicieron comparaciones entre las áreas medidas y calculadas, concluyendo la viabilidad del esquema de Mc Cormack en la ruptura de diques o en el análisis de ruptura de presas. Hay varios estudios relacionados con modelos matemáticos que intentan representar el fenómeno físico de una onda de avenida. (Fennema et al., (1986), Brufra y García-Navarro, (2000), Michaud, et al.,( 2005)). Nguyen et al, 2006, validaron con éxito un modelo bidimensional de flujo superficial usando un UFVM (método de volúmenes finitos no estructurado) en varios casos entre ellos la idealización de una ruptura parcial y total de una presa.  En este estudio se analizó la eventual ruptura de un tanque de enfriamiento de la Planta termoeléctrica José López Portillo, ubicada en Río Escondido, Coah., México. Dicha planta tiene cuatro unidades de 300 MW que utilizan carbón como combustible; el sistema de enfriamiento es de tipo abierto y requirió la construcción de un tanque con una superficie de cerca de 300 ha (3x106 m2). El tanque fue construido sobre el terreno natural utilizando diques de una longitud de 13 km (13x103 m), con una altura hasta de 10 m, el agua utilizada es bombeada utilizando un acueducto con una longitud de 30 km (30x103 m), desde el río Bravo.                      Figura 1. Ubicación del  sitio de estudio, Río Escondido, Coah.            METODOLOGÍA Modelo bidimensional del flujo a superficie libre   Las ecuaciones de Saint Venant, en su forma conservativa, aplicables al análisis del flujo bidimensional a superficie libre, se pueden escribir  de la siguiente manera:                                                                                                   [1] En la ecuación anterior:      ;      ;                           [2] Donde:  es el vector de variables de flujo,  es el tensor de flujo y  es el término independiente o término fuente.   El primer término de la ecuación 1 representa la variación temporal local de la masa y la cantidad de movimiento; el segundo término representa la variación espacial de los flujos de dichas cantidades; el término independiente representa la ganancia o pérdida de masa y de cantidad de movimiento por unidad de tiempo en un volumen diferencial. Debido a que la masa es constante en el problema, la primera componente del vector de variables independientes es cero. Estas ecuaciones son válidas para flujos con velocidades verticales pequeñas, pendientes del fondo del cauce suaves y cuando las dimensiones horizontales predominan sobre la vertical.  Las ecuaciones de Saint Venant constituyen un sistema de ecuaciones diferenciales en derivadas parciales hiperbólico, cuasi-lineal y con término independiente.   Método de volúmenes finitos  El modelo CARPA utiliza el método de volúmenes finitos que ha probado ser muy útil para desarrollar esquemas de alta resolución en dos dimensiones para dichas ecuaciones. Los volúmenes finitos permiten discretizar el dominio con regiones irregulares y por lo tanto adaptadas a los contornos y las ecuaciones se pueden usar en su forma integral.  Integrando las ecuaciones de Saint Venant sobre un volumen finito se llega a la expresión siguiente para la evolución temporal de las variables dependientes:                                                                         [3] Donde:  es el tensor numérico de flujo, o sea  evaluado en la pared  del polígono i,  es el número de lados;  es el vector normal exterior a la pared wk y  es su longitud.  representa la integral del término independiente sobre el volumen finito. En la ecuación anterior, la expresión del flujo numérico , en función de los valores que toman las variables en los elementos próximos, es lo que marca la diferencia entre esquemas numéricos.      Modelos de simulación Modelo CARPA en ambiente GiD  El algoritmo numérico en volúmenes finitos utilizado por el programa CARPA  está basado en el esquema WAF TDV que se puede interpretar como una extensión de los sistemas de ecuaciones asociados al esquema desarrollado por Lax-Wendroff,  o dicho de otra manera, es extensión con aproximación de segundo orden del esquema de Roe; se basa  en el método de Godunov junto con los también llamados aproximación a la solución de Riemann de Roe (Bladé y Gómez, 2006). Dentro del ambiente Gid el problemtype conocido como CARPA  se selecciona; con dicho problemtype se pueden importar o crear modelos digitales del terreno; posteriormente se le establecen condiciones de frontera al modelo, considerando una análisis en una o dos dimensiones o una combinación de ellas.  Se asignan condiciones iniciales además de que se especifica si el problema es en una o en dos dimensiones; adicionalmente se pueden proporcionar el coeficiente de rugosidad de Manning, según el tipo de material, o incluso según el tipo de uso de suelo. Los datos del cálculo también se proporcionan y con ello se puede ejecutar el cálculo.    Modelo en diferencias finitas  Para el flujo en el cauce se considera el método de Sánchez y Fuentes para flujo no permanente (Fuentes y Sánchez, 1986); este método considera como flujo unidimensional al movimiento del agua en el cauce natural; se limita al flujo en cauces y a volúmenes de control constantes en el espacio en los que solo se modifica el tirante hidráulico. El modelo matemático resuelve las ecuaciones de flujo no permanente a superficie libre en un cauce natural.  Las ecuaciones de flujo no permanente a superficie libre para cauces no prismáticos se consideran de la siguiente manera:                                                                                                                   [4]                                                                                                [5]  Donde: Y  es la elevación de la superficie del agua respecto a un plano horizontal de referencia (suma del tirante más la cota de plantilla), en m; A, área hidráulica, en m²;   Velocidad media, en m/s; B Ancho de la superficie libre,en m.  g,  aceleración de la gravedad, en m/s2, Sf , pendiente de la línea de energía, adimensional, X,  distancia longitudinal, en m, T , tiempo en s.  Para representar en diferencias finitas a las ecuaciones 4 y 5 se toman en cuenta los volúmenes de control mostrados en la siguiente figura.                Figura 2. Secciones a lo largo del cauce   Con ayuda del modelo en diferencias finitas para el modelo de flujo no permanente en el cauce se determina  Yi en la etapa de cálculo j+1 así como las velocidades de llegada (Vi,j+i) y de partida (Ui,j+1).  En las secciones inicial y final deben fijarse las condiciones de frontera, las cuales se determinan de acuerdo con el problema en estudio.  El carácter dinámico de las inundaciones en las planicies hace necesario emplear modelos matemáticos que, por lo menos, incluyan ecuaciones de flujo en dos dimensiones (Fuentes, et al 1997) a través de las ecuaciones de conservación de cantidad de movimiento  y ecuación de continuidad:                                                                                                 [6]                                                                                                [7]                                                                                                                       [8]   donde: , , pendiente de fricción en las direcciones  y  adimensional; , aceleración de la gravedad, en m/s2; , , componentes de la velocidad en las direcciones  y , adimensional, , nivel de la superficie libre del agua con respecto al nivel del terreno, en m, , coeficiente de rugosidad según la fórmula de Manning, en s/m1/3 , , , , direcciones del sistemas de ejes cartesiano derecho; , tiempo, en s.  Para calcular el flujo en una planicie de inundación se debe resolver el sistema de ecuaciones anterior considerando condiciones iniciales y de frontera; dado que no existe un método analítico para encontrar la solución, se propone un método numérico de diferencias finitas.  Sea el área de inundación en proyección horizontal dividida en celdas con longitud a lo largo del eje x  y del eje y  (Figura 3).  Figura 3. Arreglo de celdas considerado en el método numérico para el área de inundación     Ecuaciones de conservación de cantidad de movimiento  Expresando las ecuaciones 6 y 7 como:                                                                                           [9]                                                                                          [10]  donde: . Así, para el componente de la velocidad en dirección del eje X, se puede expresar en diferencias finitas del modo siguiente:                                                                                  [11]  donde: ,  y   y   Para resolver la Ecuación  11 se consideran dos casos  y .  a)  La velocidad  tiene que ser positiva con lo cual la ecuación 11 resulta una ecuación de segundo grado, cuya solución es:                                                                           [12] b) La velocidad  tiene que ser negativa, con lo cual la ecuación 11 queda como una ecuación de segundo grado, cuya solución es:                                                                              [13] Siguiendo un razonamiento similar para el componente de la velocidad en dirección del eje  (ecuación 10), tenemos:                                                                                 [14]  siendo: ,  y   y   De esta forma para:                  ;                                                        [15]                  ;                                                           [16]  Ecuación de continuidad La ecuación 8 se puede expresar en diferencias finitas  y al ordenar términos se obtiene:                                            [17]  Con las ecuaciones 12 a 17 se obtienen los valores de las componentes de la velocidad en dos direcciones (  y  ) así como la elevación de la superficie del agua , , en el tiempo  para las celdas ubicadas en el interior de la zona donde ocurre la inundación, considerando un gasto de entrada a cualquier celda en la malla expresado como un gasto unitario.   Determinación del hidrograma experimental de salida por la brecha   El tanque de enfriamiento analizado está compuesto por varios diques: el más externo, el central y un deflector, construido en su mayoría por arcilla y con un filtro central de arena como protección contra posibles filtraciones o grietas y una pared de enrocamiento. En la Figura 4 se ve una vista en planta y  así como la curva elevaciones capacidades. Un cuerpo de agua se localiza al norte del tanque y se considera útil para regular la avenida producida ante la posible ruptura de un dique. Para estimar la forma del hidrograma de salida se hicieron pruebas de laboratorio obteniendo resultados experimentales y posteriormente con ayuda de similitud dinámica (Le=64) se obtuvo la forma del hidrograma del prototipo (Figura 5), dicho hidrograma se utilizó para alimentar los modelos de simulación.  Figura 4. Planta general del tanque de enfriamiento y curva elevaciones capacidades   Figura 5. Determinación del hidrograma experimental       APLICACIÓN Datos de entrada a los modelos de simulación  En la Tabla 1 se presentan los datos de entrada más representativas aplicados en los modelos de simulación. Tabla 1. Datos de entrada Concepto CARPA-GiD 2D II UNAM Malla Tin 8253 elementos (triángulos) Grid 100x100 Método Volúmenes finitos Diferencias finitas Dt (s) 1.0 1.5 Tiempo de cálculo (s) 1,360 (11.5 h) 6,240 (11.5 h) n manning 0.032 0.040 Lenguaje Fortran, TCL/TK Visual Basic    RESULTADOS Mapas de profundidades  En la Figura 6 se presenta una comparación de resultados obtenidos por los dos métodos, para el instante 11.5 h (hora y media después de que se presenta el pico en el hidrograma de ruptura), tiempo en el cual inicia la atenuación de la avenida.  Con el Gid se obtiene un mapa directamente en el cual se presenta la escla gráfica de profundidades, además de que se pueden efectuar animaciones de la simulación; mientras que con el programa de flujo bidimensional  se requiere de trabajo adicional con empleo de otros softwares para poder conseguir el mapa presentado, que contiene menores detalles.                              Figura 6. Comparación de mapas de profundidades t=11.5 horas    CONCLUSIONES  Se comparó la aplicación del programa CARPA en el ambiente GiD con un modelo en volúmenes finitos y de un modelo de flujo bidimensional en diferencias finitas desarrollado en visual Basic. Con el programa CARPA  se lograron obtener mapas detallados  de las profundidades en distintos instantes de las posibles zonas afectadas, además de tener la posibilidad de efectuar simulaciones con el mismo software. Las herramientas analizadas permiten obtener información sumamente útil en materia de protección civil.  REFERENCIAS   Aureli, F. and P. Mignosa, (2001). “Comparison Between Experimental and Numerical Results of 2D Flows due to Levee-Breaking”. XXIX IAHR Congress Proceedings, Beijing, China, pp 16-21.  Bladé, C. E. and V.M. Gómez (2006). Modelación del flujo en lámina libre sobre cauces naturales. Análisis integrado en una y dos dimensiones. Monografía CIMNE No. 97, U.P.C., Barcelona, Spain, 227 pp.        Brufra, P. and P. García-Navarro (2000). ”Two dimensional dam break flow simulation”. Int. J. Numer. Meth. Fluids , 33, pp 35-57.   Fennema, R.J. and M.H. Chaudhry (1986).  “Explicit Numerical Schemes for Unsteady Free-Surface Flows with Shocks”. Wat. Resour. Res., 22 No.13, pp 1923-1939.   Michaud, J., Johnson, C., Iokepa, J. and Marohnic, J. (2005). “Methods for Estimating the Impact of Hypothetical Dam Flood Analysis”. Environmental Sciences Research, Vol. 59, Chemistry for the Protection of the Environment, 4, pp 195-199.         Nguyen, D.K, Shi, Y.E., Wang, S.S.Y. and Nguyen, T.H. (2006). “2DShallow-Water  model using unstructured finite-volumes methods”. Journal of Hydraulic Engineering. ASC. March, 2006 pp 258-269.       Fuentes M., O.A. y Sánchez B., J.L. (1986) "Cálculo de avenidas provocadas por la ruptura de una  presa", XII Congreso Latinoamericano de Hidráulica, Sao Paulo, Brasil, octubre.12 pp.   

Simulación de la onda de avenida por ruptura de tanque de enfriamiento

Se llevó a cabo la simulación de la ruptura del borde de un tanque de enfriamiento

usando dos métodos: Uno de volúmenes finitos denominado CARPA y un algoritmo en diferencias finitas centradas. El algoritmo en volúmenes finitos CARPA utilizado con el pre y post procesador GiD mostró ser una herramienta muy poderosa en la animación de planicies de inundación, útiles en la interpretación de resultados tanto en el espacio como en el tiempo, para la protección civil, así como para definir las

posibles zonas afectadas debido a fenómenos como el de la ruptura de bordos

UPCommons

https://upcommons.upc.edu/bitstream/2117/14748/1/5436509.pdf

Simulación de la onda de avenida por ruptura de tanque de enfriamiento

Abstract

Similar works

Full text

Available Versions

UPCommons. Portal del coneixement obert de la UPC

UPCommons