1. INTRODUCCIÓN
En el informe presentado por la ex primera ministra de noruega Gro Harlem Brundtland ante la ONU se nombra por primera vez el término desarrollo sostenible, el cual se define como: "la satisfacción de las necesidades de la generación presente sin comprometer la capacidad de las generaciones futuras para satisfacer sus propias necesidades" (ONU, 1987). De igual forma, según las Naciones Unidas, el cambio climático supone efectos adversos, referidos a cambios en el medio ambiente físico o en la biota con consecuencias nocivas significativas en la composición, la capacidad de recuperación o la productividad de los ecosistemas naturales o sujetos a ordenación, o en el funcionamiento de los sistemas socioeconómicos, o en la salud y el bienestar humano (Pinto y Salazar, 2014). Por consiguiente, uno de los objetivos mundiales en pro del desarrollo sostenible es "Garantizar el acceso a una energía asequible, fiable, sostenible y moderna para todos" (PNUD, 2015), por lo tanto, los países en vía de desarrollo tienen un mayor compromiso para lograr esta meta. Paralelamente, el Congreso Colombiano decretó mediante la ley 1715 del 2014 "promover el desarrollo y la utilización de las fuentes no convencionales de energía, principalmente aquellas de carácter renovable" (COLOMBIA, 2014). Por todo lo anterior, es importante desarrollar tecnología que mejore los sistemas de generación y transformación de energía renovable.
Colombia está en un punto intermedio entre los trópicos de cáncer y capricornio, lo que indica que es una zona de alta radiación solar e incidencia de vientos alisios del noreste, vientos constantes que recorren estas latitudes a una velocidad media de 5.55 m/s. Sin embargo, la orografía y otros factores modifican las corrientes atmosféricas cambiando la condición natural de los vientos a velocidades inferiores de la media (GÓMEZ, 2012).
En consecuencia, se estudia la forma de aumentar la potencia de un sistema eólico, utilizando para ello el Principio Venturi en los vientos que impulsan el aerogenerador. Asimismo, se utilizan las velocidades del viento presentes en el territorio Colombiano iguales e inferiores a 3.5 m/s (Ver figura 1), siendo esta la velocidad mínima promedio a la que arrancan los aerogeneradores comerciales.
El concentrador de flujo eólico fundamenta su diseño en un túnel cónico con dos secciones circulares abiertas y transversales a la dirección del viento, en el cual se aplica el efecto Venturi; de tal forma, que el aerogenerador se ubica en la sección menor del túnel donde se ve aumentada la velocidad del viento, por el principio de continuidad. Al respecto concierne decir, que el potencial eólico por unidad de área es proporcional a la densidad y a la velocidad del viento al cubo, lo que indica, que si se aumenta la velocidad del viento en el paso por el aerogenerador, se incrementa en un orden cúbico el potencial eólico disponible.
2. METODOLOGÍA
Para realizar el estudio del concentrador eólico se asumieron tres etapas:
2.1 Análisis experimental
El análisis se realiza bajo un modelo a escala del sistema eólico, teniendo presente el análisis de antecedentes e investigativo (Bell & Mehta). En el experimento se evalúan las velocidades del viento tanto en la entrada como en la salida del concentrador con el objetivo de identificar el perfil de velocidades circundantes que se generan en el concentrador para una posterior validación. El ensayo está compuesto por los siguientes elementos:
Un anemómetro CK120197.
Un concentrador eólico con un ángulo de ataque de 20.5° y con una relación de áreas de 4.
Un ventilador de piso industrial de marca Home Elements HEHAF-18.
En las figuras 3 y 4 se ilustran los perfiles de velocidades obtenidos del experimento.
De las figuras 3 y 4 se observa que el perfil P6 es alterado por la disposición del concentrador eólico, dando como resultado el perfil P6". Esto indica, que las velocidades en la entrada del concentrador son inferiores a las velocidades circundantes en el ambiente (Ver tabla 1). Sin embargo, las velocidades a la salida del concentrador son superiores en orden de 0.9 m/s a las velocidades circundantes del ambiente.
2.2 Modelado
Para determinar la longitud de operación se siguen los parámetros de Bell & Mehta. Donde:
Por lo tanto, para no estar al mínimo de generar números de Reynolds altos se promedian los límites de la longitud de operación (Longitud larga: genera el crecimiento de la capa límite. Longitud corta: genera disturbios en la uniformidad del flujo a medida que la sección transversal cambia (MOLINA, 2012) adquiriendo una longitud óptima para la geometría del concentrador. Permitiendo obtener la tabla 2.
Para obtener el perfil de contracción se acude a la ecuación polinomial de quinto orden de Bell & Mehta, la cual tiene como objetivo principal minimizar la longitud de contracción de manera que se obtenga un mínimo de Reynolds a la salida.
Para adquirir el perfil de la contracción se utiliza el software de Matlab, el cual permite procesar los datos. Una vez desarrollado el código en Matlab se obtienen las coordenadas y se exportan a Solid Edge (Ver figura 5 y 6, se ilustran los perfiles de contracción para cada relación de áreas).
Los datos exportados son 1000 coordenadas en los ejes X y Y. Se adicionan 1000 coordenadas de Z cero para ubicar el perfil en el boceto de solid edge.
Una vez almacenados los datos se procede a realizar el modelo del concentrador de flujo eólico mediante las herramientas de Solid Edge (Ver figura 7).
El flujo másico es el mismo en todas las secciones del concentrador de flujo eólico. Por ende, se debe cumplir la ecuación de conservación de masa para el sistema. (Ver figura 8).
Asumiendo que las propiedades del fluido permanecen constantes en el tiempo y partiendo de que el concentrador eólico tiene dos superficies de control; entonces, se puede despreciar la parte temporal y simplificar la ecuación a:
Una vez conocidas las velocidades y las perdidas por presión según la tabla de Kisieliev, se realiza un balance de energía para determinar la presión de salida mediante la siguiente ecuación.
2.3 SIMULACIÓN
El análisis de fluidos en el concentrador de flujo eólico se realiza mediante simulaciones en CFD, ANSYS Workbench.
El dominio computacional tiene una forma cilíndrica con medidas de 6 veces la longitud de la geometría aguas abajo, 3 veces aguas arriba y con un diámetro de 5 veces de la sección trasversal del sólido (Ver figura 9).
El mallado híbrido está conformado por elementos tetraédricos para el volumen del fluido en general y prismáticos para la captura del comportamiento de la capa límite adyacente a la pared del concentrador. Los elementos de la malla global están conformados por: tamaño mínimo de 4 mm, tamaño máximo de 120 mm, tasa de crecimiento para las celdas en el dominio de 1.2 y una transición suave con 5 capas de inflación con una tasa de crecimiento 0.172 (Ver figura10). No se presentó la necesidad de realizar mallas locales.
La cantidad de los nodos y elementos están sujetos a la disposición de la geometría y el dominio computacional. De los cuales, 3348420 son tetraédricos y 662920 prismáticos, conectados por 952537 nodos.
La calidad ortogonal alcanzó un rango muy bueno con un promedio de 0.86. Donde, 0 es inaceptable y 1 es excelente. De igual forma, la oblicuidad obtuvo un rango excelente con un promedio de 0.22. Donde, 1 es inaceptable y 0 es excelente.
Análogamente, se establece el Shear Stress Transporte SST K-ω como modelo de turbulencia con un fluido de densidad 1.2 kg / m 3 y viscosidad dinámica de 1.8x10-5 N s/m 2 . Simultáneamente, la velocidad de entrada se establece constante con un valor de 3 m/s y es perpendicular a la superficie de la entrada del dominio computacional con una intensidad de turbulencia porcentual de 2.6. También, se establece una presión a la salida de 0 pa.
Se selecciona el algoritmo de solución SIM-PLE-C para el acoplamiento presión - velocidad, el cual se programa para 1000 iteraciones, en el cual la solución presenta convergencia a las 516 iteraciones con un criterio de convergencia de 1e- 6 para velocidades en X.
Como resultado de la simulación se muestran seis líneas internas (series) de las velocidades del fluido en el concentrador (Ver figura 11), con una velocidad inicial de 3 m/s. las velocidades en la entrada (0 m) del concentrador están en el orden de 1.7 m/s y a la salida (1.04 m) del mismo alcanzan valores de 3.9 m/s. Sin embargo, 20 cm después de la salida del concentrador, estas velocidades están por el orden de 4 m/s. Las velocidades a la salida presentan una tolerancia aproximada de 0.05 m/s, lo que indica, que estas velocidades muestran una homogeneidad en toda el área de salida.
En la figura 12 se ilustran las posibles trayectorias que pueden seguir las partículas del fluido para una velocidad inicial de 3 m/s, estas trayectorias asignan magnitudes superiores a 4 m/s. El flujo se establece como estado estacionario.
Los mapas de contorno (Ver figura 13 y 14) ilustran que las velocidades en la entrada del concentrador son inferiores a las velocidades circundantes del ambiente (3 m/s). Lo que indica, que el concentrador en la entrada actúa como una resistencia al paso del fluido. Sin embargo, las velocidades máximas están por el orden de 4 m/s cerca a la salida del concentrador, debido al efecto Venturi. Estas velocidades son homogéneas a cierta distancia en la cual se podría ubicar el rotor del aerogenerador.
Las figuras 15 y 16 ilustran las presiones presentes en las líneas internas del concentrador, estas presiones toman valores positivos antes de 0.8 m en el interior del ducto, debido a que las velocidades del fluido han disminuido en la entrada. Después de esta distancia toman valores negativos dado que las velocidades han aumentado por encima de los 3 m/s. Lo que indica, que a medida que incrementa la velocidad disminuye la presión.
RESULTADOS
Para la validación de los resultados se toman los datos obtenidos en el desarrollo experimental y los obtenidos en la simulación de prueba. En el cual, el promedio de las velocidades circundantes del ambiente de la parte experimental está en 2.7 m/s y el promedio de velocidades a la salida del concentrador "modelo" está en un valor de 3.6 m/s, el cual se coteja con el promedio de velocidades a la salidad de la simulación de prueba de 3.4 m/s. determinando un error relativo porcentual de 5.5.
Una vez realizada la simulación de prueba y la validación de los datos, se efectúan las simulaciones para el diseño del concentrador. Donde, inicialmente se realizan 6 simulaciones con diferentes relaciones de áreas.
Análogamente, se comparan los resultados de las seis simulaciones entre sí determinado que el Perfil del Concentrador 4 (PC4) presenta velocidades considerables a la salida con menores dimensiones de diseño (volumen y material). Sin embargo el PC5, PC6 y PC7 presenta un ligero aumento a la salida, pero no justifica tal aumento respecto al volumen del concentrador (Ver figura 17).
Posteriormente, se realizan dos simulaciones más al PC4, con la diferencia que se lleva al límite la longitud de operación del concentrador con base en los parámetros de Bell & Mehta. En el cual, se logra determinar que el Perfil del Concentrador 4 Corto (PC4C) genera mayores velocidades a la salida del orden de 4 m/s, optimizando así la longitud del mismo (Ver figura 18).
DISCUSIÓN
La validación de las simulaciones en los diseños finales se asume sobre el comportamiento semejante en la solución del algoritmo de acoplamiento presión-velocidad a la simulación de prueba. No obstante, se debe tener presente que esta última se validó con el análisis experimental, donde dicho análisis se valora de manera empírica debido a las limitaciones de los recursos.
CONCLUSIONES
Se realizaron un total de 9 simulaciones, la primera de ellas está orientada a la validación del análisis experimental, las siguientes 6 están enfocadas a identificar la relación de área que aporte mayores velocidades a la salida del concentrador eólico con los mínimos recursos de diseño y las 2 últimas están encaminadas a calcular el aumento de la velocidad respecto a la variación de la longitud de operación de Bell & Mehta. Por último, se determina que el concentrador PC4C presenta un aumento de 1 m/s las velocidades del fluido (Aire), llevando la potencia eólica disponible de a, lo que representa un incremento 2.37 veces.
El concentrador eólico PC4C tiene como dimensiones, una longitud axial de 1.04 m con una relación de áreas de 4, siendo su diámetro de salida de 1.17 m. Lo que permite menores volúmenes de diseño respecto al PC4, PC4L, PC5, PC6 y PC7.
La mayoría de los aerogeneradores comerciales de baja potencia requieren una velocidad mínima de arranque de 3.5 m/s. No obstante, las velocidades del viento en el territorio de Colombia están por el orden de los 3 m/s según el mapa de vientos de la UPME. Esto implicaría que se requieren mayores velocidades para poner en marcha a los aerogeneradores. Sin embargo, el concentrador eólico permitiría poner en funcionamiento a los aerogeneradores con estas velocidades, debido a que este aumenta los vientos a 4m/s. Por lo tanto, se generarían mayores tiempos en el funcionamiento.