1. Introducción
La capacidad de producir energía a partir de recursos renovables es una línea de investigación en creciente desarrollo, la cual busca dejar la dependencia de recursos fósiles y tecnologías que contaminan al medioambiente. Enmarcados en el Objetivo de Desarrollo Sostenible (ODS) 7 para energía asequible y no contaminante, la energía proveniente de mares y océanos posee un enorme potencial que podría ser explotado de una manera más eficiente y económica. Además, los patrones existentes en el océano son más predecibles que otras fuentes de energía renovables de naturaleza estocástica (Bellec, 2017). Sin embargo, es necesario que en los próximos años exista un crecimiento acelerado en cuanto a tecnología para que este tipo de energía contribuya de forma importante en el cuidado ambiental y el desarrollo sostenible (Chowdhury et al., 2021). La implementación de modelos matemáticos y simulación numérica para diseñar y predecir el funcionamiento de mecanismos de transformación de energía es una técnica muy utilizada en la actualidad, ayudando a mejorar la eficiencia y rendimiento de estos.
Con el objetivo de modelar y predecir el comportamiento de las olas en el mar se crean los llamados Numerical Wave Tanks (NWT), herramienta que con una creciente accesibilidad a dispositivos con alto poder computacional se vuelve muy útil en la ingeniería marina para diversos campos de estudio (Marques et al., 2018). Estos proporcionan una valiosa herramienta de pruebas en la experimentación y optimización de dispositivos generadores de energía, además de ser de alta fiabilidad cuentan con una buena rentabilidad por no requerir de costos de instalación física. De acuerdo con Windt et al. (2018) con la tecnología en constante crecimiento, los modelos numéricos pueden ser mejorados y optimizados para reducir tiempo de simulación y costo computacional.
Varios estudios han utilizado la dinámica de fluidos computacional para desarrollar modelos numéricos para dichos tanques. Alamian et al. (2015) presentaron un modelo de tanque de olas para estudiar olas regulares generadas por un pistón en un tanque 2D, por medio del método BEM, el cual ayudó mucho para calibrar tanques experimentales. Kim et al. (2016) utilizando CFD realizaron estudios de cargas de olas sobre estructuras submarinas fijas, implementando un NWT con gran éxito. Bruinsma et al. (2018) estudiaron la interacción fluido-estructura para estructuras flotantes con amarre. Para esto utilizaron un NTW donde la propagación se realizó mediante teoría de olas Lineal y 2nd Order Stokes. Otros estudios se centraron en optimizar tanques numéricos para mejorar la predicción del comportamiento de las olas al implementar diferentes técnicas de simulación y de mallado (Sun et al, 2016; Penalba et al, 2018; Chen, et al., 2019). En los últimos años se han utilizado modelos CFD para caracterizar la interacción fluido-sólido en dispositivos de generación energía que utilizan aire, con resultados satisfactorios y que han ayudado en el diseño de dichos mecanismos (Andrade-Terán, 2022; Rodríguez et al, 2021).
En el estudio de convertidores de energía a partir de las olas llamados ‘dispositivos WEC’, por sus siglas en inglés (Wave Energy Converter), en español denominados ‘tecnología undimotriz’; los NWT se introducen como una técnica útil y fiable para simular el funcionamiento de los convertidores. Zhi et al. (2018) aplicaron la técnica de NWT en el estudio del funcionamiento de un Overtopping Wave Energy Converter (OWEC) utilizando CFD con olas regulares mostrando resultados con alta fiabilidad respecto al experimental. E. J. Ransley ha participado en múltiples estudios (Ransley, 2015; Ransley, Greaves, Raby, Simmonds y Hann, 2017; Ransley, et al., 2017) de modelos numéricos para predecir el funcionamiento de WEC en diferentes condiciones. En dichos estudios se utilizaron modelos con malla dinámica aplicando 6DOF y modelos multifásicos VOF para la interacción aire-agua en el tanque numérico, tomando como principal modelo de WEC el WaveStar Machine uno de los dispositivos más prometedores para la transformación de energía a partir de las olas, que fue ideado en el año 2000 por los entusiastas de la navegación Niels and Keld Hansen y que ha sido desarrollado constantemente en las últimas décadas (Piera, 2019). Sin embargo, en Latinoamérica han sido escasos los estudios referentes a este tipo de tecnología y esfuerzos por asimilar propuestas existentes. Por lo tanto, el propósito de este estudio es realizar simulaciones de la dinámica existente cuando interactúa el WEC del tipo WaveStar Machine con un flujo de olas controlado, para predecir el funcionamiento del convertidor en condiciones típicas de la región.
2. Metodología
Para analizar la interacción fluido-sólido en la transformación de energía proveniente de las olas, se reproduce un modelo numérico CFD basado en el modelo físico instalado en Coast Ocean Wave Basin (Uneversity of Plymouth, 2021) escalado 1:5. Del modelo operacional INOCAR-SWAN (Instituto Oceanográfico de la Armada del Ecuador, 2022) se obtienen datos típicos de estados de mar en las costas ecuatorianas, para ser estudiados. Varios períodos de olas son simulados partiendo del reposo, para analizar la evolución en el tiempo del flujo y el desplazamiento que se provoca en el WEC por los diferentes estados de olas.
2.1 Geometría dispositivo y tanque de olas
Las pruebas se realizan tomando como dispositivo generador un WEC del tipo Wavestar Machine, el cual consiste en un absorbedor flotante hemisférico, un pistón hidráulico, un rotor y una estructura fija al fondo del mar o a una plataforma flotante (Jakobsen et al., 2016). El movimiento del flotador mientras sube y baja por las olas se transfiere a través de la hidráulica al rotor para producir electricidad (Kramer, Laurent y Peter, 2011). La geometría del dispositivo utilizada se basa en el modelo instalado en la Universidad Plymouth con las dimensiones descritas en los trabajos de Ransley et al. (2017) y de Windt et al. (2020) de un brazo simple del Wavestar Machine, como muestra la figura 1.
El dispositivo mecánicamente opera bajo las siguientes condiciones físicas y mecánicas (se toma como origen de coordenadas el centro geométrico de la boya):
- Centro de gravedad: (-0.17; 0.35; 0) m
- Centro de rotación: (-1.57; 1.68;0) m
- Masa total: 220 kg
- Inercia aproximada:124 kg m.
El tanque de olas, figura 2, donde se colocará el WEC para las simulaciones geométricamente es modelado con las siguientes dimensiones: 16 m de largo, 3 m de profundidad, 3.75 m de altura libre y 6 m de ancho. El dominio computacional del tanque es simplificado implementado un eje de simetría en la dirección de no incidencia de las olas. El sistema de conversión de energía será dispuesto en el tanque numérico, aplicando la simplificación propuesta por Windt et al. (2020) donde se coloca geométricamente solo la boya flotante, ya que es la única parte del WEC que entra en contacto directamente con el fluido; esta se encuentra instalada estratégicamente 4.97 m de la entrada del flujo y 2.6 m del fondo del tanque en la posición de flotabilidad estable, 0.4 m de la boya sumergida.
2.3 Simulación numérica
Numéricamente se modela la propagación de olas sobre el dispositivo generador de energía simulando la interacción fluido-solido. Inicialmente el WEC es fijado en la posición de flotabilidad estable para simular el flujo de olas con la boya como obstáculo, para así conseguir datos necesarios para validar la dinámica del fluido y optimizar el cálculo del NWT. Luego se habilita un único grado de libertad rotacional en el eje z al WEC, para simular el movimiento que causan las olas al sólido.
2.3.1 Ecuaciones de gobierno
El flujo transitorio en el NWT es resuelto por simulación numérica directa (DNS) mediante ecuaciones de conservación de masa (Ecuación 1) y de cantidad movimiento (Ecuación 2) para flujo laminar.
Además, un modelo multifásico de volumen de fluido (VOF) es utilizado para caracterizar la interacción de las fases (aire-agua) en la superficie libre en donde se tiene α = 1 cuando la celda está llena de agua y α = 0 llena de aire. La ecuación 3 es la que rige la variación de fracción de volumen.
Donde:
: son las fuerzas gravitacionales y fuerzas externas ejercidas [N]
t: tiempo [s].
: viscosidad del fluido [Pa s].
La teoría Second Order Stokes es utilizada para la generación de olas en el modelo multifásico VOF, esta se rige por la ecuación diferencial parcial de Laplace, ecuación 4 (Maâtoug y Ayadi, 2016).
Donde:
X: coordenada horizontal de incidencia directa para la ola
: número de ola angular [rad/m], = ( donde es la longitud e ola
:corresponde a la inclinación de la ola [rad].
Dinámica del sólido
La teoría 6DOF calcula movimientos tridimensionales a partir de las fuerzas y momentos aplicados a la superficie (ANSYS Inc., 2013). El movimiento del WEC en el tanque se modela limitando a un único grado de libertad rotacional con respecto al eje Z, la expresión que rige esta teoría está descrita por la ecuación 5.
Donde:
L: es el tensor de inercia.
: el vector de momento del cuerpo [Nm].
es la velocidad angular del cuerpo [rad/s].
2.3.2 Malla
La malla propuesta presenta refinamiento en el área de influencia de las olas y las zonas cercanas al sólido para lograr mejor precisión del flujo y disminuir la cantidad de celdas necesarias para caracterizar el fenómeno. Esta se construye mediante celdas hexaédricas estructuradas mediante la técnica de mallado por bloques utilizando ANSYS ICEM CFD, la malla obtenida es mostrada en la figura 3.
Virtualmente se establece una región de malla deformable suficiente para un libre movimiento de la boya y otra región con las celdas cercanas al sólido que siguen el movimiento pasivamente sin deformarse. Aplicando técnicas de Smoothing y Layering se permite la deformación de las celdas con el movimiento de la boya, sin perder la calidad de la malla. La discretización estructurada por bloques permite calidad alta en la malla, donde únicamente el 0.006 % son celdas con calidad mínima de entre 0.35 a 0.5, el resto de la malla presenta calidad superior a 0.55; estas celdas poseen un índice de ortogonalidad superior al 0.49. La independencia de malla se realizó monitoreando un punto en el dominio del NTW aguas arriba de la boya, logrando independencia en cuanto a presión en el punto monitoreado con una malla con 239 832 celdas y 256 779 nodos.
2.3.3 Configuración de solución
Las simulaciones CFD han sido resueltas utilizando ANSYS Fluent con doble precisión para control de decimales. El método utilizado para el acople Presión-Velocidad en el esquema PISO (Pressure-Implicit with Splitting of Operators) activando skewnnes-neighbor con correcciones, para mejorar la convergencia cuando se deforme la malla. En la discretización espacial se empleó: el esquema PRESTO! para la presión, Second Order Upwind para el momento y un esquema Compresive para la discretización de la fracción de volumen. La pared superior del tanque se establece como salida libre a la presión atmosférica, simulando un tanque abierto. Para la condición de entrada se activa la teoría de OpenWave, estableciendo un estado inicial tipo Flat para un tanque sin olas al iniciar la simulación. La superficie libre del líquido en 0 m y el fondo del tanque a -3 m permitiendo 3.75 m para el flujo e interacción del aire. La condición de salida del flujo se establece con Pressure-outlet activando Open Chanel manteniendo el fondo del tanque a -3 m y nivel de agua en 0 m.
Las condiciones para E.O.2 y E.O.3 son escaladas 1:5 de los datos reales obtenidos tanto para la altura significativa como para el período promedio de los estados de mar típicos durante el mes de abril. Datos recolectados de las boyas virtuales en las zonas cercanas a las costas de Esmeraldas y al golfo de Guayaquil, respectivamente, figura 4.
Los estados de oleaje utilizados, enlistados en la tabla 2, se establecen en las condiciones de contorno mediante Second Orden Stokes parametrizando las condiciones de mar necesarias.
2.4 Transmisión del movimiento
Con el WEC libre se controla un grado de libertad rotacional al registras la posición angular (θ) del centro de gravedad de toda la estructura cada time step respecto de la posición inicial del WEC. Dicho moviente de la boya provoca el desplazamiento del cilindro hidráulico (Xc) el cual se comprime y se extiende con el flujo de las olas transmitiendo aceite a presión al rotor. La ecuación 5, es la utilizada para trasladar geométricamente el movimiento de la boya hacia el pistón hidráulico. Las dimensiones y ángulos de interés considerados en el cálculo están presentados en la figura 1.
Donde:
a= 1.385 [m]
b= 0.89728 [m]
3. Resultados y discusión
El tanque inicialmente con aguas tranquilas con la boya en su estado estable de flotabilidad es irradiado con olas durante 10s, logrando extender la propagación de las olas a lo largo de todo el tanque, necesario para analizar la dinámica del conjunto fluido-solido.
3.1. Validación
La elevación de la superficie libre (FSE), a través del tiempo del WEC fijado en la posición estable, simulando el estado E.O. 1 obtenida de la simulación numérica es comparada con las pruebas de flujo realizadas experimentalmente, como se muestra en la figura 5. Se toma como puntos de control dos puntos del dominio computacional, los cuales se encuentran a 4.77 m aguas arriba (figura 5a) y 2.04 m aguas abajo (figura 5b) respecto al centro geométrico de la boya. En rojo entrecortado se encuentran los datos del estudio experimental y en línea negra continua, datos del numérico realizado en el presente estudio.
Se puede observar que las curvas obtenidas son muy similares a las presentadas a partir de las pruebas experimentales, demostrando que todas las simplificaciones realizadas al modelo numérico no disminuyen la precisión y confiabilidad del cálculo. Los resultados presentan un error de predicción respecto al modelo experimental de entre 3.6 y 4.7 % en el flujo de olas E.O. 1, comparando los máximos y mínimos de los valores obtenidos. La propagación de las olas hasta encontrarse con la boya muestra una congruencia temporal de décimas de segundo, variación de altura de olas de menos de 1 cm y luego de la boya la dispersión es mínima siendo el error de flujo aceptable para el estudio.
La ecuación 5 es aplicada para calcular el desplazamiento del cilindro hidráulico Xc, transformando el desplazamiento angular libre de la estructura, a un solo punto de absorción de energía. La figura 6, presenta una comparación de los datos obtenidos para el desplazamiento del cilindro en la simulación E.O. 1 con los datos del experimental (Ransley et al., 2017).
El desplazamiento del cilindro Xc en el numérico muestra un ligero desfase temporal respecto al experimental, en cuanto a la amplitud existe un error máximo del 6 %. En un lapso aproximado de 3s, los resultados oscilan de una forma irregular, esto se debe a que el dispositivo no alcanza la estabilidad al iniciar la simulación (Ransley et al., 2017).
3.2 Respuesta del dispositivo a diferentes estados de olas
El libre movimiento del brazo del WEC con centro de giro en el punto pivote por la interacción con el fluido es ilustrado en la figura 7, donde se establece una isosurface bajo el criterio VOF= 0,5. Esta representa la elevación de la superficie libre del agua (FSE), mostrando contornos de velocidad del fluido alrededor de la boya.
Como se observa en la figura 7, existe mayor velocidad de flujo en el estado de olas 3 (E.O. 3) donde existe mayor altura en las olas con períodos más cortos. Para los estados 1 y 2 (E.O. 2) se observa flujo con menor velocidad, donde la altura de las olas es menor en períodos más largos. Sin embargo, el desplazamiento que provocan los estados E.O.1 y E.O.2 son superiores a los provocados por las olas del E.O.3, como lo demuestra la gráfica en la figura 8. Esto podría deberse a que un flujo muy rápido no permite desarrollar completamente la propagación de la ola. Períodos muy bajos provocan desorden en las olas, mientras que períodos largos permiten a la ola mantener mayor energía mientras se desplazan (Kayac Autovaciable, 2019).
La respuesta del dispositivo, al impacto de las Olas sobre la boya, se ve reflejado en el desplazamiento angular (θ) del brazo del WEC con respecto al punto pivote. Durante las simulaciones se ha registrado a través del tiempo este desplazamiento (figura 8), partiendo desde el punto de flotabilidad estable se presenta la oscilación del brazo para cada prueba durante los 10 s simulados.
Estos resultados sugieren mejor desempeño del convertidor con los estados 1 y 2, los cuales presentan olas con períodos más largos. Asimismo, este desplazamiento angular es trasladado al movimiento lineal del cilindro hidráulico (Ecuación 5), donde se denota mejor transformación de energía desde las olas en E.O.1 y E.O.2 que en E.O.3. La figura 9 muestra el desplazamiento que se genera en el cilindro hidráulico Xc, en respuesta al movimiento del WEC provocado por el paso de los diferentes tipos de olas estudiadas, denotando las condiciones más favorables para la transformación de la energía a partir de las olas.
El flujo de olas E.O.2 provoca un desplazamiento del pistón Xc en promedio 5.5 % superior al dato del estado de control E.O.1, significando energéticamente una potencia superior en dicho porcentaje por cada flotador que se instale en el sistema del WEC. Por otro lado, el estado E.O.3 muestra un decaimiento en promedio del 67 % en la dinámica del pistón absorbedor por lo que se reduciría significativamente la potencia hidráulica del mismo.
3.3 Transferencia de energía
La potencia (Po), por metro de columna de agua, existente en los estados de mar estudiados se calculan utilizando la ecuación 6 (Perez, 2018). La potencia mecánica en el brazo del WEC en respuesta, es calculado aplicando la ecuación 7, considerando el valor del momento (M) ejercido respecto al punto pivote y la velocidad angular del brazo analizada a partir del desplazamiento (Fernández et al., 2018).
Donde:
H: Altura de Ola.
T: Período de Ola.
La tabla 2 presenta la potencia disponible de las olas en los estados E.O.1, E.O.2 y E.O.3 y la potencia mecánica del brazo del WEC provocada por dichos estados.
El dispositivo con el estado de mar E.O.3 presenta menor potencia mecánica a pesar de existir propagación de olas con altura significativa superior al resto. El estado de mar E.O.2 aporta una de potencia de oleaje inferior, sin embargo, la potencia mecánica generada en el brazo del WEC para este estado es 50 % superior a la existente con el estado E.O.3.
4. Conclusiones y recomendaciones
La metodología de simulación numérica utilizada permitió obtener alta precisión en los resultados del desplazamiento del sólido y mínima distorsión en la propagación de las olas. Obteniendo así resultados con alta fiabilidad en el estudio del funcionamiento del WEC, presentando un error de predicción de milímetros en la propagación de las olas y de 1 a 2 cm en la dinámica del pistón actuador. Esta precisión en los resultados permitió analizar la generación de energía que se conseguiría con el WEC funcionando bajo las condiciones en las costas de Ecuador.
Las simulaciones denotan que el flujo de olas con altura significativa alta y períodos más largos favorecen la transformación de energía undimotriz. El mayor desplazamiento del pistón absorbedor de energía en los estados E.O.1 y E.O.2 lo demuestran. En estos, el brazo del WEC presenta un desarrollo completo de la dinámica con desplazamiento angular entre 3 y 4 ° y hasta 4 cm de desplazamiento en el cilindro hidráulico a escala 1:5. También se observa que, a pesar de presentar una alta disponibilidad energética proveniente de las olas, un estado de mar con flujo muy rápido (olas de períodos cortos y altura pronunciada) no asegura la generación de energía elevada. Estos estados ocurren en puntos con poca profundidad de mar, como lo es el punto estudiado cerca al golfo de Guayaquil, presentando aproximadamente un 50 % de potencia mecánica inferior al punto en la Costa de Esmeraldas.
Es recomendable realizar más estudios de funcionamiento, replicando el presente estudio, en diferentes zonas de la región en Sudamérica para poder empezar a explotar el potencial existente en los océanos y mares. A la par es necesario realizar estudios del impacto que causaría el dispositivo al ecosistema donde se lo instale, ya que le principal objetivo es desarrollar dispositivos de generación de energía con mínimo impacto ambiental.