Spanish (pdf)
Article in xml format
Article references
Automatic translation
Send this article by e-mail
Cited by SciELO 
Cited by Google 
Similars in
SciELO 
Similars in Google 
Permalink
INTRODUCCIÓN
El Caribe colombiano ha experimentado fenómenos erosivos desde finales del siglo XVIII, y esta situación ha empeorado desde la década de 1970 debido a causas naturales y humanas (Correa, Alcántara-Carrió y González, 2005; Posada y Henao, 2008). Esta situación tiene un efecto negativo en la gente que vive en las áreas costeras, que es aproximadamente 12,5 % de la población del país, puesto que el retroceso de la línea de costa y de los arrecifes ocasiona pérdida de manglares, de tierra, y daños a las construcciones humanas (Alosnso, Sierra-Correa, Arias-Isaza y Fontalvo, 2003, Posada y Henao, 2008).
Entre los fenómenos naturales que ocasionan la erosión se encuentran: viento; mareas (que en el Caribe colombiano se caracterizan por ser mixtas tipo semidiurnas, con predominancia de los componentes diurnos, y micromareal con amplitudes menores a 0,5 m); aumento del nivel medio del mar y el oleaje, el cual tiene dirección principal oeste, periodos de entre 6 a 9 s, y alturas entre 1 m y hasta 3 m en condiciones de tormenta, y pueden producir escapamientos de 10 cm a 50 cm, en el sentido de las corrientes (Correa, Alcántara-Carrió y González, 2005; Posada y Henao, 2008; Invemar y Ministerio de Ambiente, Vivienda y Desarrollo Territorial, 2003; Lonin y Escobar, 2009).
Para proteger las zonas costeras de la erosión, alrededor del mundo se han realizado diferentes propuestas de estructuras tanto sumergidas como flotantes, las cuales se han analizado con modelos numéricos. Un ejemplo de esto se tiene en el estudio realizado por Dorrel (EP 0 628 113 B1, 2004), de un arrecife artificial formado por estructuras huecas, con el cual predijo la elevación de la superficie libre antes y después de la estructura. Por su parte, Jiang, Liang, Zhu y Liu (2016); Liu, Zhang, Li y Xie (2013); Mendonça et al. (2011); Srisuwan y Rattanamanee (2015), realizaron simulaciones numéricas de arrecifes artificiales y evaluaron su capacidad para atenuar el oleaje, para la protección costera. Chen, Wang, Dong y Zheng (2012) adelantaron una modelación hidrodinámica en 2D de un rompeolas flotante; Diamantoulaki y Angelides (2010) analizaron el efecto de la configuración sobre la efectividad de rompeolas flotantes; Ji et al. (2006) estudiaron diferentes estructuras flotantes en condiciones de oleaje regular; Michailides y Angelides (2012) usaron modelos tridimensionales para evaluar el comportamiento de estructuras flotantes; Peng, Lee, Shin y Mizutani (2013) compararon parámetros entre modelos físicos y numéricos.
El objetivo de este artículo es investigar el comportamiento hidrodinámico aproximado de dos tipos de estructuras de disipación de oleaje, sumergida y flotante. Para esto se seleccionaron cinco estructuras básicas sumergidas (prisma triangular, cubo, esfera, cono, y cilindro), configuraciones grupales de estas y una placa flotante como estructuras de protección costera. Su desempeño fue evaluado para la reducción de la altura del oleaje bajo condiciones típicas de oleaje del Caribe colombiano, a diferentes profundidades, para establecer cuál estructura presenta el mejor funcionamiento para las diferentes condiciones. Inicialmente se explica la metodología; es decir, se detallan las geometrías, configuración del modelo numérico y escenarios de simulación; posteriormente, se realiza el análisis de resultados considerando el efecto de la profundidad, geometría y los campos de velocidad en la reducción del oleaje transmitido por cada geometría; por último, se mencionan las conclusiones del trabajo.
METODOLOGÍA
Cinco estructuras básicas sumergidas, prisma triangular, cubo, esfera, cono y cilindro, configuraciones grupales de estas geometrías y una estructura flotante fueron analizadas en un modelo CFD, con el software Flow3D con el fin de lograr una primera aproximación a condiciones asociadas hacia estas diferentes geometrías mediante modelación CFD. En la figura 1 se muestran las estructuras a evaluar.
Fuente: elaboración propia.
Figura 1 Configuraciones de estructuras a evaluar a) Geometrías básicas; b) Configuraciones grupales; c) Estructura flotante
El dominio computacional fue el siguiente:
Geometrías básicas: 40 m de largo, 3 m de ancho y alto de entre 5 m y 6 m.
Configuración grupal: 60 m de largo, 6,4 m de ancho y alto de entre 5 m y 6 m.
Estructura flotante: 60 m de largo, 6,4 m de ancho y alto de 8 m.
Debido a la extensión del dominio computacional y con la finalidad de reducir costos computacionales durante las simulaciones se usaron dos bloques de mallas: una malla gruesa de celdas de 0,10 m que abarca todo el dominio computacional, y una malla fina que se concentra en la geometría, para así obtener resultados con mayor precisión alrededor de esta, con celdas de 0,05 m. Para la estructura flotante, se usó un solo bloque de malla con celdas de 0,15 m teniendo en cuenta que por lo menos existan 4 celdas en la longitud más pequeña de la estructura para lograr un óptimo renderizado durante las simulaciones.
Se simuló un tren de olas lineal (basado en la teoría de Airy, la cual asume que el tren de olas proviene de un reservorio de fondo plano) como condición de frontera de entrada que se desplaza a lo largo del eje X, con una condición de frontera outflow en la salida del modelo, junto con condiciones de frontera laterales tipo simetría y una condición de frontera superior de tipo presión. Los módulos físicos usados para esta simulación fueron el renormalization group model (RGM) como modelo de turbulencia, y air entrainment.
Con respecto a las variables de simulación, para cada una de las estructuras se evaluaron dos profundidades (2 m y 3 m), de tal forma que estas se mantuvieran en condiciones sumergidas, y dos alturas de ola con su respectivo periodo asociado (altura de 0,88 m con un periodo de 13,58 s y altura de 2,71 m con un periodo de 7,91 s), las cuales son representativas de las condiciones del Caribe colombiano, específicamente en las playas de los departamentos de La Guajira y Magdalena, donde se ha evidenciado la problemática de erosión costera en los últimos años. Por lo que se obtienen 20 casos de simulación para las geometrías básicas (tabla 1).
Tabla 1 Casos a evaluar para las geometrías básicas y configuraciones grupales, y flotantes
Fuente: elaboración propia.
Para las configuraciones grupales se evaluaron los mismos 20 casos que se usaron para las geometrías básicas. Para la estructura flotante se evaluaron cuatro casos, con las mismas profundidades y alturas de ola analizadas para las geometrías básicas y las configuraciones grupales (tabla 2).
Se tomaron dos puntos de medición, antes y después de la estructura, en donde se midió la elevación de la superficie libre. A partir de esto, se calcularon las alturas de ola antes y después de la estructura; con estos datos se obtuvo el coeficiente de transmisión por medio de la ecuación (1).
En donde Hi y Ht corresponden a las alturas de olas incidentes y transmitidas, respetivamente, por la estructura evaluada (CERC, 1984). Debido a que esta investigación es una primera aproximación, no se realizó una calibración con modelos físicos. Si bien una modelación física es necesaria para poder calibrar el modelo numérico, estas simulaciones se hacen como una primera aproximación con la cual es posible obtener tendencias, comportamientos generales y órdenes de magnitud para las estructuras evaluadas, y así facilitar una mejor planeación para pruebas físicas que se puedan llevar a cabo en el futuro.,
RESULTADOS
En el presente estudio se realizaron simulaciones numéricas con el fin de evaluar el efecto sobre el coeficiente de transmisión de las olas, producido por diversas geometrías de estructuras sumergidas y flotantes.
Geometrías básicas individuales
Se analizaron los resultados para establecer el efecto de la profundidad y geometría en el flujo.
Efecto de la profundidad
En general, para el prisma triangular, cubo, esfera y cilindro, los menores coeficientes de trasmisión se presentan en el caso 1, en donde se tienen la menor profundidad, y la menor altura de ola; los mayores coeficientes de trasmisión están en los casos 2 y 4, en los que se tiene la mayor altura de ola. Para una misma profundidad se encontró que a una mayor altura de ola, los valores para el coeficiente de trasmisión aumentan.
El cono presenta un comportamiento irregular con respecto a las demás estructuras. La gran mayoría de los puntos presentan un coeficiente de transmisión por arriba de 0,6, y los mínimos valores del coeficiente de transmisión se dan cuando la profundidad del agua es la mayor evaluada (figura 3). Se utilizaron gráficos de cajas y bigotes para ilustrar los resultados, en los cuales la línea central representa la mediana, los extremos de la caja los percentiles 25 y 75, y los bigotes llegan hasta los puntos que no son atípicos; es decir, aquellos que son mayores que el percentil 75 más 1,5 por la diferencia entre el percentil 75 y 25, o menores que el percentil 25 menos 1,5 por la diferencia entre el percentil 75 y 25.
Efecto de la geometría
Para analizar el efecto de la geometría se realizó una comparación entre las cinco estructuras básicas, evaluándolas en las profundidades seleccionadas.
Los resultados para todas las estructuras en una profundidad de 2 m se muestran en la figura 4. En estos casos, las estructuras presentan comportamientos diferentes. El coeficiente de transmisión para el cono es el mayor en comparación con las demás estructuras. Esto es consistente con el hecho que la relación entre el área proyectada ocupada y el área total para el cono es la menor de todas las estructuras, con un valor de 0,5. Esto da lugar a que el agua se desplace entre los espacios laterales vacíos del cono, posiblemente generando estos resultados.
Fuente: elaboración propia.
Figura 4 Coeficiente de transmisión para geometrías básicas para h=2,00 m
Las estructuras para las cuales el área proyectada ocupada es mayor presentan menores valores de coeficiente de transmisión, como el prisma triangular, cubo o cilindro, las cuales tienen un valor de 1 de área proyectada, mientras que la esfera tiene un valor de 0,78, con menores resultados.
La geometría que presenta los menores valores de coeficiente de transmisión es el prisma triangular, en la condición de menor altura de ola. En la figura 5 se pueden apreciar los resultados de todas las estructuras para una profundidad de 3 m.
Fuente: elaboración propia.
Figura 5 Coeficiente de transmisión para geometrías básicas para h=3,00 m
En este caso, todas las estructuras se comparan a la misma profundidad de 3 m. El prisma triangular presenta los menores valores del coeficiente de transmisión. Adicionalmente, a medida que el área proyectada ocupada de las estructuras crece, los valores para el coeficiente de transmisión disminuyen.
En general, las estructuras básicas con las cuales hay una mayor disminución de la altura del oleaje transmitido, en orden descendente son: prisma triangular, cubo, cilindro, esfera y cono.
Campos de velocidad y esfuerzo cortante
Para la interacción fluido/estructura, se sabe que el comportamiento del fluido varía de acuerdo con la estructura en el medio. En la figura 6, se ilustran los campos de velocidad alrededor de la estructura, después del paso de la primera ola.
En dicha figura también se muestra el efecto que cada estructura tiene en el fluido; en el prisma triangular, el cubo y el cilindro, se crean vórtices al lado derecho, en la esfera se crea un pequeño vórtice en el lado de la estructura y en el cono no se crearon vórtices significativos. Estos vórtices representan los patrones de velocidad generados debido a la interacción del fluido con la estructura, y pueden indicar zonas en donde es posible que puedan presentarse en un futuro problemas de socavación. La diferencia entre los vórtices que se dan en cada estructura puede deberse a las características geométricas de las estructuras que producen cambios en el comportamiento del flujo alrededor de la estructura.
Relacionado con los vórtices está la velocidad del fluido alrededor de la estructura, la cual incrementa los esfuerzos cortantes en el lecho; de esta manera se ocasionarían problemas de erosión y socavación. En la figura 7 se ilustran los esfuerzos cortantes alrededor de la estructura después del paso de una ola, los cuales pueden asociarse con la velocidad. Las geometrías con un área proyectada de 1 (prisma triangular, cubo y cilindro) muestran los esfuerzos cortantes máximos al lado derecho de las estructuras, al nivel de la base, mientras que la esfera y el cono señalan los esfuerzos cortantes máximos a ambos lados de las estructuras.
En la figura 8 se observan las líneas de flujo para el prisma triangular, cono y esfera, respectivamente. Para las geometrías que ocupan todo el espacio en la sección transversal, las líneas de flujo son rectas; el área que recibe el impacto del flujo es mayor, y puede asociarse a un menor coeficiente de transmisión. Para las geometrías que no ocupan todo el espacio en la sección transversal, las líneas de flujo rodean la estructura y los valores del coeficiente de transmisión tienden a ser mayores.
Configuraciones grupales
Se analizaron los resultados de los escenarios evaluados para las geometrías básicas en grupo. Se estudió el efecto del cambio en la profundidad y de la geometría de cada grupo.
Efecto de la profundidad
Entre los casos evaluados se tiene que la profundidad afecta el coeficiente de transmisión en los casos 1 y 3, es decir, en los que se evaluó la menor altura de ola y presentan menores valores en la menor profundidad. El cono, de igual manera que para las geometrías básicas individuales, muestra un comportamiento irregular en comparación con el resto de las geometrías (figura 9).
Efecto de la geometría
Los valores obtenidos para los coeficientes de transmisión son similares para las diferentes configuraciones, tanto para la condición de profundidad de 2 m, como para la profundidad de 3 m. En promedio el grupo de estructuras cónicas muestra el menor desempeño, con valores promedio de coeficiente de transmisión de 0,83 en el caso de la profundidad de 2 m (figura 10), y valores promedio de 0,80 para el coeficiente de transmisión en la profundidad de 3 m.
Fuente: elaboración propia
Figura 10 Coeficiente de transmisión para configuraciones grupales para h=2,00 m
En general, de las configuraciones evaluadas, las que presentan los mayores porcentajes de reducción de flujo de energía medio local en orden descendente son: configuración de cubos, cilindros, prismas triangulares, esferas y conos, aunque los resultados obtenidos para el prisma triangular, el cilindro y la esfera son muy similares.
Campos de velocidad y líneas de flujo
El campo de flujo alrededor del grupo de estructuras se ve afectado debido a la interacción de las líneas de flujo con las caras de las estructuras. En la figura 11 se observan las líneas de flujo alrededor del grupo de cubos (lado izquierdo) y conos (lado derecho). El color de las líneas de flujo corresponde a la variable de velocidad.
En el caso del grupo de cubos existe una mayor interacción y, por ende, una mayor afectación sobre las trayectorias de las líneas de flujo, en comparación con las trayectorias en el caso del grupo de conos. El grupo de cubos obtuvo los menores valores de coeficiente de transmisión, lo cual puede obedecer a la alta cantidad de colisiones y cambios de trayectorias que se dan entre los espacios del grupo de cubos. Por el contrario, el grupo de estructuras cónicas presentó altos valores del coeficiente de transmisión, lo cual se puede asociar a las pocas colisiones frontales de las líneas de flujo producto del cambio suave en la dirección de las líneas de flujo debido a la superficie redondeada de las caras de los conos.
Estructura flotante
La figura 12 muestra, a través de un diagrama de cajas, los resultados obtenidos en cuanto al coeficiente de transmisión para la estructura flotante.
Los valores más bajos de reducción de energía del oleaje para la estructura flotante ocurrieron con una altura de ola de 0,88 m con una profundidad de 3,0 m, mientras que para la misma condición de profundidad, pero con una altura de ola mayor (2,71 m), se observó una gran variación en los valores del coeficiente de transmisión en cada una de las olas que se trasmitieron sobre la estructura flotante.
CONCLUSIONES
Se tienen las siguientes conclusiones, de acuerdo con los resultados obtenidos:
Entre las estructuras básicas, la que presenta un menor coeficiente de transmisión es el prisma triangular.
La estructura cónica es la que posee el mayor coeficiente de transmisión tanto para la estructura básica como para la configuración grupal.
Todas las estructuras para la configuración grupal, excepto el cono, tienen un comportamiento similar, aunque el menor coeficiente de transmisión se da con la configuración grupal de cubos.
La estructura flotante tiene menores coeficientes de transmisión cuando la estructura se encuentra más cerca de la superficie.
Para una misma profundidad, la configuración grupal cúbica tiene un mejor desempeño que la estructura flotante.
