Modelación numérica integrada de los procesos hidráulicos en el drenaje urbano

Resumen

El crecimiento continuo de las ciudades, o la ocurrencia de un evento extraordinario de lluvia, en muchas ocasiones sobrepasa las medidas de diseño de algunos de los elementos del drenaje urbano (principalmente elementos de captación y colectores) provocando que el sistema de drenaje de una población sea incapaz de desalojar el escurrimiento por dicho evento. El resultado final, es una inundación. Las consecuencias de ello son problemas de circulación para vehículos y personas, daños materiales, e incluso la pérdida de vidas humanas, en función principalmente de los valores de las variables hidráulicas como el calado y la velocidad del flujo. En este contexto, es clara la necesidad de contar con herramientas numéricas para reproducir dichos eventos y proponer actuaciones. Por este motivo, en este trabajo se presenta el desarrollo de un modelo numérico integrado uni/bidimensional, para el cálculo de los procesos hidráulicos que ocurren en el drenaje urbano de forma totalmente acoplada. Por un lado el drenaje urbano se compone de diferentes elementos y por el otro el modelo numérico integrado se forma de diferentes modelos numéricos, el resultado de ello y para efectos prácticos y de modelación numérica son modelos de tejados, calles, elementos de captación, conductos y unión de conductos. Cada modelo dispone de al menos una opción de modelación numérica, que permite calcular el proceso hidráulico en el elemento correspondiente. De esta forma, el modelo numérico integrado se puede utilizar de forma conjunta, o cada uno de los modelos de forma separada. El modelo de tejados emplea las ecuaciones de la onda cinemática en una dimensión y las ecuaciones de aguas poco profundas en dos dimensiones para modelar el proceso de transformación lluviaescurrimiento y la propagación del flujo en los tejados. Por su parte el modelo de calles únicamente emplea las ecuaciones de aguas poco profundas en dos dimensiones para calcular el proceso de transformación lluvia?escurrimiento y la propagación del flujo en las calles. Es importante destacar que, para que la modelación sea más realista, el flujo proveniente de los tejados puede descargar ya sea de forma puntual o de forma uniforme tanto en las calles como en los colectores. En los colectores y en las uniones de una red de alcantarillado ocurren fenómenos hidráulicos complejos, de ahí que se cuente con diferentes alternativas para su cálculo. De esta forma, en el modelo de conductos, se emplean las ecuaciones de Saint Venant en una dimensión para el cálculo del flujo en lámina libre. Para el cálculo de flujo en presión se dispone de tres opciones: el método de la ranura de Preissmann, dos modelos de dos ecuaciones (ecuaciones de flujo en lámina libre y ecuaciones de flujo en presión): uno considera flujo incompresible en presión y otro flujo compresible en presión. Además, se dispone de una formulación para considerar el efecto que el aire atrapado (ya sea delante de un frente de onda de presión o entre dos frentes de onda de presión) tiene sobre la hidrodinámica del flujo. En el modelo de la unión de conductos se cuenta con dos aproximaciones: la formulación en una dimensión, empleada por la mayoría de modelos y una formulación en dos dimensiones basada en las ecuaciones de aguas poco profundas en dos dimensiones. La interacción entre el sistema de drenaje mayor (tejados, calles, áreas verdes, áreas grises, canales naturales y canales artificiales, etc.) y el sistema de drenaje menor (colectores, uniones, compuertas, etc.), conocida como drenaje urbano dual, es un proceso complejo pero fundamental, se realiza en los puntos físicos correspondientes a la localización espacial de los elementos de captación (fuentes y sumideros). De esta forma en el modelo de los elementos de captación, los caudales de intercambio de este proceso, se pueden estimar con la metodología Flumen y/o empleando las ecuaciones tipo orificio y vertedero. Desde el punto de vista numérico, para la solución de las ecuaciones de flujo en lámina libre en una dimensión, las ecuaciones de flujo en presión en una dimensión y las ecuaciones de aguas poco profundas en dos dimensiones, se emplean esquemas numéricos en volúmenes finitos, capaces de calcular discontinuidades y frentes de onda en la solución, sin necesidad de técnicas adicionales. En cambio, para la solución de las ecuaciones de la onda cinemática en una dimensión, se utiliza un esquema en diferencias finitas. Además, se usa el método de las características para la imposición de condiciones de contorno y el método de Newton?Raphson para la solución de sistemas de ecuaciones no lineales. Para cada uno de los modelos numéricos implementados en primer lugar se evalúa su funcionamiento y desempeño, posteriormente se valida y en algunos casos se realiza alguna aplicación, primero de forma separada y posteriormente empleando más de un modelo numérico Para ello, se emplean casos de referencia, resultados de otros modelos numéricos, soluciones analíticas y ensayos de laboratorio existentes en la literatura técnica. Los resultados de este proceso indican que los modelos numéricos son capaces de reproducir de forma aceptable y con un nivel de precisión adecuado los tipos de flujo modelados: flujo en lámina libre en régimen lento y en régimen rápido, flujo en presión positiva y negativa, además de las transiciones que se pueden presentar entre estos flujos, como frente seco mojado y cambio de flujo en lámina libre a flujo en presión (flujo mixto), incluso con la presencia de aire atrapado. Por otra parte, el modelo numérico integrado se aplica a dos cuencas urbanas totalmente distintas entre sí, por su geometría, las características de los elementos del drenaje urbano, y las condiciones y variedad de flujos que ocurren. Los resultados obtenidos permiten decir que el modelo es capaz de resolver esta variedad de flujos de forma razonable y funcionar de forma integrada entre los elementos del drenaje urbano. Ello permite pensar que, el modelo numérico integrado se puede utilizar en el cálculo de los procesos hidráulicos del drenaje urbano, con un desempeño satisfactorio.