Introducción

En la actualidad, la creciente demanda de energía ha impulsado diversas e innovadoras maneras de aprovechar el recurso eólico para la producción de electricidad. Fei, Mai y Li (2012) describen un sistema de recolección de energía eólica basado en las vibraciones mecánicas que se inducen en una cinta flexible debido a la interacción del fluido con la estructura. Perkovic, Silva, Ban, Kranjcˇevic´ y Duic´ (2013) presentan un mecanismo que consiste en un volantín diseñado para transformar la energía cinética del viento, a gran altitud, en energía mecánica.

Por otra parte, las turbinas eólicas convencionales (artefactos que extraen energía cinética del viento) están en constante rediseño e investigación para obtener cada vez mejores productos en cuanto a eficiencia, confiabilidad y armonía con el entorno (Burton, Sharpe, Jenkins & Bossanyi, 2011).

Las turbinas eólicas de eje horizontal son las más comúnmente diseñadas hoy en día (Manwell, McGowan & Rogers, 2009). Poseen una torre que soporta el eje de rotación a una determinada altura sobre el suelo y paralelo a este, de modo que se aprovecha la energía del viento que viaja a una velocidad mayor en tanto más alta sea la torre que soporta a la turbina.

Existe una amplia gama de formas de turbinas de eje vertical; entre ellas destaca la Darrieus, que ampliamente fue estudiada en Estados Unidos y Canadá en las décadas de 1970 y 1980 (Manwel et al., 2009). Aún hoy en día se presentan modelos matemáticos y computacionales para el análisis de este tipo de turbinas (Islam, Ting & Fartaj, 2008).

En el campo de las turbinas eólicas de eje horizontal, recientemente se ha explorado la alternativa de utilizar, en lugar de aspas, cilindros girando que aprovechan el efecto Magnus para generar la fuerza de sustentación (Bychkov, Dovgal & Kozlov, 2007). La Figura 1 muestra el aerogenerador conocido como turbinas Magnus, mostrando sus parámetros geométricos y cinemáticos (diámetro D y longitud L del cilindro, velocidad angular de la turbina y del cilindro ω) así como el viento incidente U _∞.

El efecto Magnus es una fuerza que actúa sobre un cuerpo en rotación sumergido en un fluido en movimiento; la dirección de esta fuerza es perpendicular a la dirección del flujo, según se muestra en la Figura 2. Esto se debe a una diferencia de presión causada por la diferencia de velocidad que existe en el fluido entre la cara superior e inferior; la modificación en el campo de velocidades sucede al girar el cilindro arrastrando consigo el fluido a su alrededor por la condición de no deslizamiento.

Bychkov et al. (2007) presentan una investigación teórica y experimental que propone la operación óptima de una turbina eólica que utiliza cilindros girando en lugar de aspas. Sus resultados para un cilindro en rotación muestran un crecimiento lineal para el coeficiente de sustentación C _L , que es muy superior al correspondiente de un perfil alar.

Ito, Kawashima, Nishizawa, Ushiyama & Komatinovic (2007) fabrican un prototipo de una turbina Magnus-Savonius, en la que los cilindros de la turbina Magnus son reemplazados por rotores Savonius. De este modo, se evita la necesidad de hacer girar los cilindros con energía externa, pues los rotores Savonius generarían su propia rotación por arrastre. Una de sus conclusiones es que este aparato tendrá un coeficiente de potencia mucho menor al de turbinas similares con aspas, pero a cambio el costo de fabricación es muy bajo.

Bychkov (2008) presenta una alternativa para la turbina Magnus, en la cual la base del cilindro no gira. En este caso, emplea tapas laterales en los cilindros, para suprimir la vorticidad en las puntas de estos.

Gono, Rusek & Hrabcˇík (2009) describen la turbina Magnus manufacturada por la empresa Mecaro Co. Destacan que en ese modelo se emplean cilindros con aletas espirales, para potenciar el efecto Magnus generado con la rotación de los cilindros.

Luo, Huang & Wu (2011) presentan una solución analítica para el desempeño de una turbina Magnus. Desprecian el efecto del arrastre para el cálculo del torque desarrollado por los cilindros, por lo que su modelo corresponde a una turbina en un flujo no viscoso.

Sun, Zhuang, Cao, Huang & Wu (2012) llevan a cabo una simulación numérica para estudiar el efecto de modificar la forma de los cilindros en una turbina Magnus, empleando ANSYS CFX. La simulación se lleva a cabo con solo un cilindro, aprovechando la periodicidad rotacional del problema e imponiendo adecuadamente las condiciones de borde. Utilizan el modelo de turbulencia SST k - ω por su precisión en el tratamiento de la capa límite ante gradientes de presión adversos. Obtienen que un cono truncado, con radio creciente alejándose del eje de rotación de la turbina, produce un coeficiente de potencia mayor que un cilindro circular. Pero para esto requiere mayores valores de ω.

Sedaghat (2014) reporta que para el caso de la turbina Magnus resulta clave lograr la reducción de la relación entre arrastre y sustentación, si se busca mejorar la eficiencia.

Massaguer, Massaguer, Pujol, Comamala y Velayos (2014) proponen generar la rotación de los cilindros de la turbina Magnus por efectos centrífugos y autorotación, evitando la necesidad de invertir energía externa para generar el efecto Magnus. Entre otros aspectos, sus pruebas concluyen que este tipo de turbina genera menos ruido que las convencionales.

La revisión bibliográfica evidencia que la turbina Magnus ha sido objeto de estudio durante los últimos años. Uno de los resultados es que la turbina tiene eventuales ventajas respecto de las tradicionales que usan aspas. No obstante, ningún autor presenta una curva de potencia, teórica o experimental, para una turbina Magnus, ni su comportamiento ante modificaciones geométricas o cinemáticas; estos aspectos forman parte del presente estudio y en eso radica su originalidad e innovación.

Metodología

Modelamiento y método propuesto

La teoría del Disco Actuador, descrito en detalle por Burton et al. (2011), explica que el viento disminuye su velocidad al atravesar un generador eólico. Para modelar el intercambio de energía cinética entre el aire y una turbina eólica se define el factor de inducción axial a, el cual es una medida de qué tanto disminuyó su velocidad el viento al salir de la turbina y, por ende, cuantifica cuánta energía entregó el aire al generador.

Mediante argumentos de conservación de momentum lineal y la ecuación de Bernoulli, Burton et al. (2011) demuestran que la fuerza horizontal (denotada como T, del inglés Thrust) ejercida por la turbina al aire y viceversa es:

(1)

Donde U _∞ representa la velocidad del viento libre, lejos de la turbina, ρ es la densidad del fluido y A _d es el área del disco que barren las aspas.

Este disco puede ser analizado en pequeños anillos, tal como lo propone la teoría BEM (Blade Element Momentum) descrita en detalle por Burton et al. (2011). Al realizar este tipo de análisis, la ecuación 1 puede ser escrita para un anillo diferencial a una distancia r como:

(2)

La ecuación 2 es válida para a< 0.4 y dentro de ella se tiene el coeficiente de empuje (C _T ).

(3)

Por consideraciones de conservación de momentum angular, Burton et al. (2011) obtienen que el diferencial de torque Q generado por el viento en cada disco es:

(4)

Donde a´ es el factor de inducción tangencial, que modela el cambio en la velocidad tangencial del viento al atravesar la turbina eólica. Aguas arriba el viento viaja horizontalmente, pero dado que interacciona con la turbina y la hace girar, se genera una fuerza de reacción que dota al viento de velocidad tangencial aguas abajo de la turbina.

La ecuación 2 pierde validez cuando la turbina opera en un estado de estela turbulenta, que ocurre cuando a> 0.4. Para compensar este efecto, se emplean dos correcciones anidadas (Moriarty y Hansen, 2005). La de Glauert reemplaza el por:

(5)

Donde F está definido según el modelo de corrección de puntas de la siguiente manera:

(6)

Aquí B representa la cantidad de aspas que componen la turbina y, de acuerdo con la figura 3, α se escribe como:

Esta corrección de puntas también se debe aplicar de forma directa en las ecuaciones 2 y 4, lo que resulta de la siguiente forma:

Por otra parte, las velocidades y fuerzas aerodinámicas que actúan sobre un elemento de aspa tal como se muestra en la figura 3, permiten escribir expresiones para T y Q en función de las magnitudes de las fuerzas de sustentación FL y arrastre FD. No obstante, estas dos fuerzas pueden ser escritas en función de los coeficientes adimensionales CL y CD de forma diferencial (Burton et al., 2011) como:

Donde c es la cuerda del aspa y Un se obtiene de la Figura 3.

Para determinar los factores de inducción axial y tangencial se combinan las expresiones para la fuerza horizontal según la conservación de momentum lineal (ecuación 8), el torque según la conservación de momentum angular (ecuación 9) y las fuerzas aerodinámicas (Figura 3 y ecuación 10) para obtener:

Las ecuaciones 11 y 12 se reescriben de manera que todos los términos estén a un lado del igual, del otro lado se colocan los residuos R1 y R2, que idealmente deben ser igual a cero. Pero al igualar dos modelos que toman en cuenta simplificaciones y desprecian efectos tridimensionales y turbulencia entre otros, se debe buscar el punto en el que ambos residuos alcancen un mínimo. De modo que se discretiza el disco actuador en anillos, y para cada posición r se buscan los valores de a y a´, que se determinan cuando se cumple que:

El método que se propone e implementa en un código para resolver la ecuación 13 es no iterativo, por tanto se le ha llamado método directo y consiste en evaluar todas las posibles combinaciones de a y a´ en la ecuación 13, considerando que ambas pueden tomar valores entre 0 y 1. El producto de esta evaluación es una superficie que representa el valor de la suma de los residuos, que están asociados a fuerza y momento que predice de forma no compensada alguno de los dos modelos empleados.

La solución válida para el sistema es el punto donde la superficie de residuo total alcanza el valor mínimo posible. Por lo que una vez evaluadas todas las combinaciones, un algoritmo extrae de forma directa los factores de inducción axial y tangencial que finalmente gobiernan el rendimiento de la turbina.

Luego de resolver la ecuación 13 para cada valor de r, se determina la potencia P:

Una vez discretizado el disco en fracciones de radio, se calcula la potencia como:

El proceso descrito se debe repetir para todo el rango de velocidad de viento en estudio.

Validación del método propuesto

Para validar el método directo propuesto se han utilizado como datos de entrada las características de la turbina Fase II (Ceyhan, 2008) del NREL (Laboratorio Nacional de Energía Renovable de EE.UU.), que se compone de tres aspas de 0.4572 m de cuerda, girando a 71.3 revoluciones por minuto (RPM), con un diámetro de 10.06 m y una potencia de 19.8 kW. Los resultados comparativos de los valores predichos por el método directo y los valores reales de la turbina se muestran en la Figura 4.

Se encuentra una adecuada concordancia entre los resultados obtenidos por el método directo propuesto y las mediciones experimentales para la turbina Fase II. Se aprecia una ligera sobreestimación de la potencia para velocidades de viento bajas y altas, menores a 9 m/s y mayores a 19 m/s. En el rango restante, el método directo presenta una subestimación en la predicción de potencia. No obstante, se captura adecuadamente la tendencia de la curva de potencia y sus valores puntuales, considerando que el modelo utilizado no toma en cuenta efectos como rugosidad, pandeo y vibración de las aspas o turbulencia y variación del ángulo de incidencia lateral del viento, entre otros.

Adecuación del método directo para turbinas Magnus

El análisis de la turbina Magnus se lleva a cabo con el método directo validado anteriormente. Las dos diferencias fundamentales que implican una modificación en el planteamiento del modelo son la geometría, al reemplazar el perfil alar por un cilindro, y la rotación del cilindro.

Es por esto que en las ecuaciones 11 y 12, el concepto de cuerda debe ser sustituido por el diámetro del cilindro. Por otra parte, los coeficientes CL y CD son independientes del ángulo de ataque, dada la simetría cilíndrica de sus aspas; estos coeficientes, en el caso de la turbina Magnus, son función de la velocidad relativa de rotación del cilindro ω=ωrc/U∞.

Para estudiar el desempeño de la turbina Magnus se emplea el método directo propuesto, prestando especial atención a la selección de las curvas de CL y CD, pues son radicalmente distintas para un cilindro en comparación con cualquier perfil aerodinámico y, por ende, contienen la clave para lograr posibles mejoras respecto de las turbinas tradicionales en algún campo de aplicación particular.

Resultados y discusión

Se programa el método propuesto en el software comercial Matlab y se toma como dato de entrada las curvas de CL y CD según las mediciones experimentales de Bychkov et al. (2007).

Los resultados de las simulaciones realizadas mediante el método directo, para una turbina de 2 cilindros, se muestran en la Figura 5, donde se han combinado cilindros de 15 y 25 cm de diámetro con velocidades de rotación en el rango de 10 a 4000 RPM y con turbinas de 1 a 3 m de radio con velocidades angulares de entre 10 y 30 RPM.

El principal hallazgo de esta combinación de variables es que, en general, se requiere alta velocidad de rotación del cilindro, del orden de 103 RPM, para lograr extraer potencia del viento, dado que las líneas celeste y rosado (ω = 10 y 100 RPM, respectivamente) están sobre el eje de cero potencia en todos los arreglos de variables para turbinas pequeñas.

Por esta razón, en la Figura 6 se presenta el resultado de una nueva combinación de parámetros donde ω = 1, 4, 8 y 12 x 103 RPM, con el fin de determinar el comportamiento de la curva de potencia en relación con el cambio en todas las otras variables.

Con este nuevo escenario es posible definir la tendencia de la curva de potencia de la turbina Magnus. Se destaca que al aumentar la velocidad angular de los cilindros, el parámetro al que es más sensible el sistema, aumenta la potencia promedio y su punto máximo en relación con la velocidad del viento incidente.

La necesidad de una alta velocidad de rotación del cilindro se debe a que la curva del coeficiente de sustentación del cilindro en rotación no se inicia en el origen de coordenadas, sino que para valores menores a ω= 0.5 se obtiene CL = 0 (Bychkov et al., 2007).

Además, para que exista el efecto Magnus se requiere un adecuado balance entre la velocidad del viento incidente y la rotación del cilindro. Por una parte, una corriente de viento muy tenue no genera ninguna sustentación, por otro lado, tampoco habrá sustentación si la velocidad del viento es excesiva, en ese caso resulta despreciable la diferencia de velocidad del viento en las caras del cilindro y el sistema se puede aproximar como si el cuerpo cilíndrico no estuviese girando.

Por otra parte, si el cilindro gira a muy baja velocidad angular, no modifica las líneas de corriente del viento y, por ende, no se genera sustentación, pero también se obtiene cero sustentación si el cilindro gira con una velocidad angular excesiva. En este caso, el fenómeno se debe a que en la vecindad del cilindro las líneas de corriente corresponden prácticamente a las de un cilindro girando dentro de un fluido en reposo y no es posible percibir el efecto del flujo de viento externo.

Considerando que ω combina la velocidad del viento con la velocidad angular del cilindro y que según estudios previos con ω= 2 se obtiene la potencia máxima de la turbina, es razonable que al aumentar la velocidad angular de los cilindros se incremente la potencia capturada, pues se genera más sustentación y a la vez que el punto de máxima potencia ocurre para una velocidad mayor del viento incidente, ya que se debe mantener el balance para mantener constante el adimensional ω= 2.

Se observa también en las figuras 5 y 6 que al aumentar la longitud de los cilindros se obtiene una mayor potencia de la turbina, del mismo modo que en las turbinas eólicas convencionales. Esto es fácilmente pronosticable, tomando en cuenta que la potencia disponible en el viento es proporcional al área barrida por los cilindros.

Cabe destacar que, a diferencia del aumento en la velocidad de rotación de los cilindros, en este caso no ocurre una modificación significativa en el valor de velocidad del viento donde da inicio la extracción de energía para cada combinación particular de las otras variables, por lo que hacer la turbina más grande no la hace ineficiente para bajas velocidades del viento, como sí ocurre cuando se aumenta ω.

En cuanto al diámetro de los cilindros, se encuentra que al aumentar su valor también crece la potencia máxima entregada por la turbina, pero la extracción de energía se inicia a velocidades del viento mayores y el punto de máxima potencia también se corre a la derecha en el gráfico de potencia en función de la velocidad del viento.

Este es un comportamiento similar al obtenido cuando se modifica ω. La mejora en la captura de potencia se debe, por una parte, a que las fuerzas aerodinámicas sobre las aspas de turbinas eólicas son directamente proporcionales a la cuerda (ver ecuaciones 10 y 11) y en el caso de las turbinas Magnus, el diámetro corresponde con el concepto de cuerda de un aspa convencional.

Por otra parte, considerando la definición de ω, que en su numerador contiene tanto a ω como al radio del cilindro, resulta razonable que un aumento en el radio cause un efecto similar al de aumentar la velocidad angular del cilindro, pues en ambos casos se requiere un aumento en el denominador (que es la velocidad del viento) para mantener ω= 2 constante y operar en el punto óptimo.

La variación en la velocidad angular de la turbina modifica la curva de potencia de modo que, conforme se aumenta se requiere una mayor velocidad del viento para que la turbina extraiga energía, pero a la vez se obtiene como beneficio que la cota máxima de potencia se incremente.

Sin embargo, para valores de similares a los que se emplean en turbinas eólicas convencionales, la potencia máxima se alcanza a valores de U _∞ muy altos (20-25 m/s), lo cual no es práctico para la mayoría de las aplicaciones. Por tanto, de manera preliminar se puede indicar que la turbina Magnus ofrece ventajas para moderados o bajos, tal como proponen autores como Bychkov et. al (2007), quienes señalan que la turbina Magnus opera de forma adecuada a un tercio de la velocidad angular de las eólicas convencionales, resultado que se obtuvo también en este trabajo.

En cuanto a la cantidad de cilindros de la turbina y su efecto en el comportamiento, el análisis demuestra que la tendencia es la misma que al aumentar el diámetro del cilindro o la velocidad angular de la turbina. No obstante, cada cilindro representa más material, mayor complejidad de construcción y un nuevo consumo de energía para generar la rotación, por lo que en este parámetro se debe incluir una nueva variable, que es el gasto energético necesario para la operación de la turbina en función de la cantidad de cilindros.

Conclusiones

Se requiere una velocidad de rotación del orden de 10³ RPM de las aspas cilíndricas de una turbina tipo Magnus para lograr extraer potencia del viento, lo que coincide con lo obtenido en investigaciones anteriores.

La variable que gobierna el comportamiento de esta turbina es la velocidad relativa de rotación del cilindro ω=ωr _c /U _∞ . En concordancia con otros autores, se encuentra que para un rendimiento óptimo, la turbina debe operar en valores cercanos a ω= 2.

La turbina Magnus resulta eficiente con un bajo número de cilindros, dado que al aumentar la cantidad de aspas cilíndricas, disminuye el máximo coeficiente de potencia alcanzable.

La eficiencia de este tipo de turbinas es directamente proporcional al diámetro del cilindro para un radio de turbina fijo, lo que no había sido reportado anteriormente en la bibliografía.

Acorde con lo señalado por otros autores, la turbina Magnus opera de forma adecuada a velocidades angulares menores que las eólicas convencionales de eje horizontal.

Recomendaciones para trabajo futuro

Agradecimientos

Se agradece al Instituto Tecnológico de Costa Rica y al Ministerio de Ciencia, Tecnología y Telecomunicaciones de Costa Rica, por el financiamiento parcial para esta investigación. Al Arq. Diego Murillo Calderón, por su colaboración con los diagramas.

Bibliografía

Fechas de Publicación

Histórico