Introducción

Las turbinas eólicas son usualmente instaladas en ubicaciones de baja turbulencia, sin embargo, el desarrollo urbano junto con la creciente ocupación de campos eólicos en tierra ha creado un mercado para la instalación de turbinas en lugares de mayor irregularidad superficial, como zonas boscosas o cerca de comunidades desarrolladas ⁽¹⁾. En estas regiones la velocidad del viento es variable, registrándose moderada (0.15>TI>0.10) a alta (TI≥0.15) intensidad de turbulencia (TI, por su nombre en inglés) ⁽²⁾. TI se define como la razón entre la desviación estándar muestral y la velocidad promedio del viento en un mismo periodo ⁽³⁾.

Es por esta razón que numerosos estudios se enfocan en el efecto de la turbulencia sobre el rendimiento de las turbinas eólicas. Se ha realizado una comparación entre las turbinas de eje horizontal con las de eje vertical; encontrándose que las turbinas de eje horizontal en pequeña escala superan su contraparte vertical, sin embargo, estas requieren mayor mantenimiento debido a ser más afectadas por las condiciones del viento ⁽⁴⁾. El efecto del viraje en condiciones turbulentas también es investigado, donde con la presencia de un concentrador de flujo se consiguió mejor rendimiento de la turbina cuando el viraje se encontraba en un ángulo entre 0° y 20° con vientos en velocidades moderadas a altas ⁽⁵⁾, ⁽⁶⁾, ⁽⁷⁾. Hay estudios enfocados en los efectos de la turbulencia sobre microturbinas aumentadas con difusores donde se obtiene in incremento en el coeficiente de rendimiento comparado con su contraparte regular ⁽⁸⁾, ⁽⁹⁾. Finalmente, la TI se relaciona con la recuperación de estela; a mayor TI se consigue una recuperación más veloz debido a la energía cinética turbulenta de estela ⁽¹⁰⁾, ⁽¹¹⁾.

Otros estudios que lidian con turbinas en flujos turbulentos reportan que bajo estas condiciones la turbina puede producir ligeramente más potencia que bajo flujos más laminares ⁽¹²⁾. También, a bajas velocidades la alta turbulencia puede resultar en un aumento en la cantidad de energía producida ⁽¹³⁾. Otros investigadores concluyen que la TI no tiene un efecto significativo sobre la potencia producida por las turbinas ⁽¹⁴⁾.

El fenómeno de alta turbulencia en perfiles aerodinámicos usualmente implica la comparación entre simulaciones numéricas y mediciones experimentales ⁽¹⁵⁾ y necesita el uso de modelos muy específicos de turbulencia para considerar las particularidades de la zona límite, así como los bordes del perfil ⁽¹⁶⁾. La optimización de la geometría del perfil es un área donde las simulaciones computacionales, junto con resultados experimentales, han ayudado a determinar que una geometría optimizada atrasa la separación del flujo, lo cual incrementa la potencia obtenida por 7.5 % (11). En ⁽¹⁷⁾ tres diferentes geometrías de generadores de vórtices se pusieron a prueba, una rectangular, triangular y una Clark-Y; esto debido a que los generadores de vórtices pueden ser utilizados para optimizar el flujo en las aspas de las turbinas, esto puede mejorar la energía anual producida de 1% a 3%. No obstante, otros estudios concluyen que las simulaciones numéricas no son suficientes para describir adecuadamente las condiciones del viento alrededor del perfil; en una de estas se descubrió una relación inversa entre la velocidad del flujo del jet y el coeficiente de arrastre ⁽¹⁸⁾.

La literatura a menudo considera la TI solamente en el eje de giro de la turbina dentro del túnel de viento o en pruebas en campo, dejando por fuera la variación de la TI en otras partes del área barrida por el rotor, obviando que esta puede variar en este espacio. ⁽¹⁹⁾ muestra una colección de puntos de medición en el área barrida por el rotor, pero no detalla como es que varía la intensidad de turbulencia en dicha área.

Como se indica por ⁽²⁰⁾, los estudios que lidian con las fluctuaciones con la velocidad del viento deben ser muy detallados para describir correctamente las condiciones reales del viento y sus efectos sobre el rendimiento de la turbina; también se da a saber que las turbinas eólicas a pequeña escala no han sido ampliamente investigadas. Más aún, ⁽⁶⁾ reconoce que hay una ausencia de literatura que se enfoque sobre los efectos de bajos números de Reynolds sobre la TI. En la poca literatura disponible donde se considera el efecto de la TI sobre la producción de potencia los rangos varían de bajo a moderados, como lo es el caso de ⁽¹⁹⁾ en donde se usan valores de TI de 1.4 %, 8 % y 13.5 %.

Es de mucha importancia saber en detalle las características del flujo a través de una turbina, especialmente en condiciones de alta turbulencia. Así, este estudio presenta un análisis detallado de la TI dentro de un túnel de viento sobre el área barrida por un rotor a distintas distancias de la cuadricula de madera.

Materiales y métodos

Equipo

Las pruebas fueron llevadas a cabo en un túnel de viento de capa límite en el Instituto de investigación para mecánica aplicada (RIAM, por su nombre en inglés) en la universidad de Kyushu en Japón. Una ilustración esquemática del túnel se muestra en la Figura 1.

Figura 1 Túnel de viento de capa límite Universidad de Kyushu Fuente 21 

La sección de prueba es de 15 m de largo, 3.6 m de ancho y 2 m de alto. Las características técnicas del túnel se pueden consultar en ⁽²¹⁾. Hay 9 paneles removidos en la sección de prueba para reducir el efecto de bloqueo ⁽²²⁾, ⁽²³⁾, así como los paneles laterales en las secciones 3, 4 y 5 del túnel de la figura 1. El panel 3 en el lado derecho (visto desde la dirección del flujo) está parcialmente removido. Esta configuración del túnel ha sido utilizada con éxito en investigaciones como ⁽²⁴⁾ y ⁽²⁵⁾.

La turbulencia fue generada con la cuadricula de madera de la figura 2, con las siguientes dimensiones: M = 156.2 mm, b = 43.5 mm y t = 29.9 mm.

Figura 2 Cuadricula de turbulencia 

La velocidad del viento fue medida con un anemómetro de hilo caliente manufacturado por Kanomax Japan Inc., Modelo 0251R-T5. Consiste en un filamento de Tungsteno de 5 µm de diámetro. Cada punto fue medido por 60 s, con una frecuencia de muestreo de 1 kHz con un filtro pasa bajas de 200 Hz. Más detalles de estos parámetros, calibración de equipo y especificaciones están disponibles en ⁽⁵⁾.

Montaje experimental

La cuadrícula de turbulencia se colocó en el marco entre los paneles 2 y 3 del túnel. Cuatro secciones transversales fueron determinadas, localizadas a 0.5 m, 0.93 m, 1.32 m y 3.1 m desde la cuadricula. Las secciones de medición no están equidistantes debido a que la turbulencia no varía linealmente dentro del túnel ⁽²⁶⁾.

Cinco puntos de medición del viento fueron simultáneamente montados en cada sección transversal, como se puede observar en la figura 3. Cuatro de los cinco puntos están separados 0.5 m del eje del túnel (su centro), localizados arriba, abajo y a ambos lados de él, mientras que el quinto punto se encuentra en lo que sería el eje de la turbina. Estos puntos fueron identificados como “Center, Up, Left, Down y Right” (sus nombres en inglés), donde su nombre corresponde a su ubicación relativa al eje del túnel. Es importante mencionar que la Figura 3 muestra una turbina en el túnel, pero esta está para fines ilustrativos solamente, las pruebas se realizaron en un túnel sin turbina, solo contaba con la cuadricula para generar la turbulencia. También, en la figura 3 el túnel se observa como cuadrado, pero el utilizado es un túnel rectangular.

Figura 3 Puntos de medición en cada sección transversal 

La TI es calculada como la desviación estándar muestral dividida por el promedio de las mediciones tomadas en un periodo de prueba ⁽²⁷⁾, donde un periodo corresponde a 60 s. Luego, este número es multiplicado por 100 para expresarlo como un porcentaje. Esta fórmula se observa en la ecuación 1.

(1)

Aquí, representa la velocidad en cada medición, es el promedio de las mediciones y n es la cantidad total de datos recolectados cada 60 s. Finalmente, en el caso en que alguna de las mediciones relacionadas con los puntos no se adhiera a la tendencia que muestren los otros, un análisis más profundo será realizado de la velocidad y la TI del punto atípico.

Resultados

Velocidad del viento

Midiendo solamente la velocidad del viento en cada uno de los puntos, se genera la figura 4.

Figura 4 Velocidad del viento promedio en cada uno de los puntos de medición 

La tendencia de la velocidad del viento es estabilizarse a un valor alrededor de 9 m/s. Los puntos más cercanos a la cuadricula muestran la diferencia más grande en velocidades, lo cual sucede debido a la cuadricula alterando el flujo (en las especificaciones del túnel, su TI regular es menor a 0.5 %, de acuerdo con la figura 1) causando anisotropía en la velocidad que es evidente en el área transversal. Conforme aumenta la distancia desde el eje, disipación viscosa y transferencia de momento permiten resultados más similares en el área barrida del rotor. Esta disminución en velocidades no sucede para el punto Up (arriba del eje), lo cual es investigado más a fondo luego.

Intensidad de turbulencia

En la figura 5 se presentan los resultados sobre intensidad de turbulencia en cada punto luego de haberse aplicado la fórmula (1).

Figura 5 Intensidad de turbulencia en cada uno de los puntos de medición 

Justo como sucedió con la velocidad del viento, la TI tiende a decrecer a mayores distancias de la cuadrícula de madera. No obstante, a diferencia de las velocidades de la figura 4 que mostraban valores altos a 1.32 m de la cuadrícula y luego empezaban a converger a valores de alrededor de 9.25 m/s; el comportamiento de la TI tiende a decrecer claramente. Esto sucede debido a la manera en la que se calcula la intensidad de turbulencia. Como la TI se calcula usando la desviación estándar de la velocidad alrededor de un punto, alta TI se puede obtener cuando el viento es, en promedio, más lento pero su desviación estándar en ese punto es mayor. Como se puede esperar, hay una variación alta de la velocidad instantánea cerca de la cuadrícula y, como se mencionó antes, esta mayor variabilidad da camino a que haya un flujo más estable conforme aumenta la distancia desde la cuadrícula.

No obstante, en la figura 5 el punto “Up” presenta un comportamiento completamente distinto al que presentan los otros cinco en la sección transversal de 3.1 m. A continuación, un análisis más profundo se realizará sobre el punto “Up” en específico.

El caso del punto “Up”

Ambas las figuras 4 y 5 muestran como el punto Up difiere en el comportamiento de su velocidad y TI respectivamente al resto de los puntos. En la figura 4, su velocidad no converge de la misma manera a como lo hacen los otros. En la figura 5, hay un claro aumento en la TI del punto “Up” comparado a la disminución en los otros.

Para comprender ambos fenómenos, es importante primero recordar las dimensiones del túnel discutidas en la sección de Equipo. Adicionalmente, hay que tomar en cuenta que el túnel no está completamente cerrado para evitar el efecto de bloqueo. En específico, es necesario considerar que el punto “Up” está tan solo medio metro debajo del exterior del túnel. Esto es algo que no sucede en los otros puntos, el punto a la izquierda y derecha del eje se encuentran a 1.3 m del exterior mientras que el punto inferior está a medio metro del piso. Esto significa que el punto “Up” es el que se encuentra más cerca de las secciones abiertas del túnel.

Si bien este es un túnel de circuito cerrado, al tener una sección parcialmente abierta algo del flujo puede salir del túnel temporalmente a través del panel removido en la sección 3 y luego ser succionado de nuevo en la sección 5 (ver figura 1), uniéndose al flujo dentro del túnel. Esta particularidad sucede principalmente debido a que el aire tiende a evitar obstáculos en su camino, por ejemplo, los instrumentos de medición.

Como el punto “Up” es el más cercano al exterior, el flujo dirigido hacia afuera es registrado por el anemómetro, lo cual explica por qué la velocidad es mayor que en el resto de los puntos en la figura 4.

El comportamiento de la TI es más complejo y necesita un mayor análisis, por esto se realizó una nueva serie de mediciones alrededor del punto “Up” con una variación vertical. Los resultados de estas mediciones se presentan en la figura 6.

Figura 6 TI y velocidad del viento en el punto Up en la sección transversal localizada a 31 m de la cuadrícula 

Como se puede observar, los nuevos puntos de medición fueron colocados a 300 mm, 400 mm, 500 mm y 600 mm desde el eje del túnel. En la figura 6 la velocidad tiende a disminuir conforme la distancia vertical aumenta. Eventualmente, la distancia fuera del túnel sería de 0 m/s.

La TI aumenta considerablemente a lo largo de la componente vertical, lo cual sugiere que la combinación de las vigas en el techo del túnel (ver figura 2), con el aire entrante ocasionado por la sección semi abierta, afectan la variabilidad en la velocidad del viento. Al suceder esto, se pierde la consistencia con los otros 4 puntos principales.

Si el punto “Up” se colocara a 0.4 m por encima del eje de la turbina en vez de 0.5 m en la sección transversal de 3.1 m, las mediciones realizadas coincidirían más con las tomadas en los otros puntos. Esta pequeña disminución de 10 cm entre puntos de medición se traduce a un significativo cambio en la TI. Esto apoya la idea de que estudios que lidian con fluctuaciones en la velocidad deben ser realizados bajo situaciones muy controladas para poder conocer las verdaderas condiciones del viento ⁽²⁰⁾.

Es importante considerar que la dirección del viento cambia hacia el exterior del túnel. Como los anemómetros solo pueden medir la velocidad a lo largo de la dirección del túnel, este cambio de dirección genera mediciones erráticas. Esto conlleva a la ilusión de una TI aún mayor.

Con todo lo anterior en mente, esta nueva sección analiza la desviación entre el promedio de las mediciones y los datos tomados en solamente el eje del túnel.

Conclusiones y recomendaciones

En este estudio, una cuadrícula de turbulencia fue usada para analizar la velocidad del viento y la intensidad de turbulencia a lo largo de diferentes secciones transversales de un túnel, tomando en cuenta el área barrida por el rotor de una turbina eólica.

Al examinar la velocidad del viento, diferencias importantes en los valores medidos en los cinco puntos en el área barrida fueron encontradas en regiones de alta turbulencia cerca de la cuadrícula, pero existe una convergencia de velocidades en las secciones transversales más alejadas de la cuadrícula lo cual llevó a una reducida TI.

Investigando la TI obtenida, se encontró que esta sigue lo establecido por la literatura, la TI no tiene un comportamiento lineal decreciente conforme se aleja de la cuadrícula.

Adicionalmente, se demostró que el porcentaje de error puede alcanzar hasta 15.2 % si la TI se considera solamente en el eje de la turbina, en vez de obtener el promedio de esta en toda el área barrida.

También, el valor de TI puede cambiar rápidamente en cortas distancias, especialmente en regiones cercanas al exterior del túnel de viento.

La principal conclusión de esta investigación es que para describir adecuadamente el flujo que incide sobre las turbinas eólicas, es necesario tomar en cuenta el promedio de las mediciones dentro del área barrida por el rotor de la turbina y no limitarlas a medir exclusivamente sobre el eje de la turbina.

Como recomendación para trabajo futuro, es necesario medir el efecto sobre el coeficiente de potencia cuando la velocidad del viento y la TI se registran en el eje de la turbina en comparación con cuando se toma en cuenta el promedio sobre el área barrida.

Contribución de los autores

Los autores declaran un porcentaje de contribución al presente artículo de la siguiente forma: Gustavo Richmond Navarro 80%; Carlos Arias Arguedas 10% y Pedro Casanova Treto 10%.