CARACTERIZACIÓN DE UN EXTRACTOR DE AIRE

ACCIONADO POR UN ROTOR EÓLICO DE EJE VERTICAL

Proyecto desarrollado con el objetivo de obtener el título de:

INGENIERA MECÁNICA

POR:

HANNY JULIET CALLEJAS REYES

ÁLVARO E. PINILLA, Ph.D

PROFESOR ASESOR

UNIVERSIDAD DE LOS ANDES

FACULTAD DE INGENIERÍA
DEPARTAMENTO DE INGENIERÍA MECÁNICA
BOGOTÁ, COLOMBIA
DICIEMBRE 16 DE 2011
 TABLA DE CONTENIDO

ÍNDICE DE FIGURAS ..............................................................................................................................iv

ÍNDICE DE TABLAS ................................................................................................................................ v
ÍNDICE DE GRÁFICAS ............................................................................................................................vi
AGRADECIMIENTOS.............................................................................................................................. vii
RESUMEN ............................................................................................................................................. 1
1. INTRODUCCIÓN ............................................................................................................................ 3
2. RELACIÓN CON TRABAJOS PREVIOS EN LA UNIVERSIDAD DE LOS ANDES .................................. 7
3. OBJETIVOS ................................................................................................................................... 8
4. CARACTERIZACIÓN DE LA TURBINA EÓLICA .............................................................................. 9
4.1. PROCEDIMIENTO EXPERIMENTAL ..................................................................................... 10
4.1.1. CURVAS CARACTERÍSTICAS DE ROTORES EÓLICOS ..................................................... 10
4.1.2. MONTAJE EXPERIMENTAL ............................................................................................. 12
4.1.3. DISEÑO DEL EXPERIMENTO PARA DETERMINAR LA TURBINA DE MEJOR RENDIMIENTO:
APLICACIÓN DEL MÉTODO DE TAGUCHI ....................................................................................... 18
4.1.4. CURVAS ADIMENSIONALES PARA OTRAS CONFIGURACIONES DE LA TURBINA ............. 21
4.1.5. CURVAS ADIMENSIONALES PARA UNA TURBINA COMERCIAL ...................................... 22
4.2. RESULTADOS Y ANÁLISIS ................................................................................................. 25
4.2.1. TURBINA DE MEJOR RENDIMIENTO: RESULTADOS DEL MÉTODO DE TAGUCHI ............ 25

4.2.2. CURVAS ADIMENSIONALES PARA OTRAS CONFIGURACIONES DE LA TURBINA ........... 31

4.2.3. CURVAS ADIMENSIONALES PARA UNA TURBINA COMERCIAL ..................................... 34
5. DETERMINACIÓN DE LA CAPACIDAD DE EXTRACCIÓN ............................................................. 37
5.1. PROCEDIMIENTO EXPERIMENTAL ..................................................................................... 38
5.1.1. CAUDAL DE EXTRACCIÓN ............................................................................................. 38
5.1.2. MONTAJE EXPERIMENTAL ............................................................................................. 38
5.1.3. VENTILADOR EXTRACTOR ............................................................................................ 40
5.2. RESULTADOS Y ANÁLISIS DE RESULTADOS ...................................................................... 42
5.2.1. CURVA DE VELOCIDAD ANGULAR VS. VELOCIDAD DE VIENTO .................................. 42
5.2.2. CURVA DE CAUDAL DE EXTRACCIÓN DE AIRE VS. VELOCIDAD DE VIENTO ................. 43

ii
 5.2.3. CURVAS DE VELOCIDAD DE EXTRACCIÓN DE AIRE VS. VELOCIDAD ANGULAR CON Y
SIN VENTILADOR ............................................................................................................................ 43

6. CONCLUSIONES ......................................................................................................................... 45
7. RECOMENDACIONES Y TRABAJO FUTURO ................................................................................ 48
8. BIBLIOGRAFÍA ........................................................................................................................... 50
9. ANEXOS ..................................................................................................................................... 51

iii
ÍNDICE DE FIGURAS

Figura 1 Potencia instalada en energía eólica en el mundo Fuente: (BRITISH PETROLEUM, 2011)
---------------------------------------------------------------------------------------------------------------------------------3
Figura 2 Instalación de Extractores Eólicos en ambiente urbano. Fuente: (INDUSTRIAS GM, s.f.)
Medellín, Colombia. ------------------------------------------------------------------------------------------------------5
Figura 3 Extractores Eólicos Empresa EXOVENT. Barranquilla, Colombia. Fuente: (EXOVENT,
s.f.) ----------------------------------------------------------------------------------------------------------------------------5
Figura 4 Creación de corrientes de aire y funcionamiento del extractor eólico. Fuente: (EXOVENT,
s.f.) ----------------------------------------------------------------------------------------------------------------------------6
Figura 5 a) Turbina instalada en el soporte estructural, se observa la brida y el rodamiento
superior. b) Soporte en aluminio. c) Rodamiento inferior ------------------------------------------------------9
Figura 6 Sección transversal turbina deje vertical y diagrama de velocidades. Tomado y Adaptado
de (MURCIA, 2007, pág. 11) y (PINILLA, Notas de Aerodinámica, 2010, págs. 8-3) --------------------- 10
Figura 7 Montaje experimental de medición de presión, humedad relativa y temperatura. -------- 14
Figura 8 Motor eléctrico acoplado al eje del rotor ------------------------------------------------------------- 15
Figura 9 Montaje experimental Péndulo Trifilar ----------------------------------------------------------------- 17
Figura 10 Configuraciones seleccionadas de la turbina evaluadas experimentalmente --------------- 20
Figura 11 Configuraciones de prueba para turbinas de pocos álabes (baja robustez) ----------------- 22
Figura 12 Características geométricas del turbo extractor comercial. Fuente: (MORALES A., 2003,
pág. 18) -------------------------------------------------------------------------------------------------------------------- 23
Figura 13 Montaje experimental en el túnel de viento para la turbina eólica comercial -------------- 24
Figura 14 Péndulo trifilar para medición de momento de inercia de la turbina comercial ----------- 24
Figura 15 Esquema vista superior de disposición de álabes para configuraciones de turbina de
pocos álabes y situaciones de bloqueo con respecto a la dirección del viento incidente. ------------ 34
Figura 16 Esquema de montaje experimental en el túnel de viento para medición de la velocidad
de extracción ------------------------------------------------------------------------------------------------------------- 38
Figura 17 Arriba: Montaje experimental para medición de caudal de aire extraído en el túnel de
viento. Abajo: Vista inferior al interior del ducto, se observa el sensor de velocidad de aire. ------ 40
Figura 18 Izquierda: Ventilador acoplado al eje de la turbina. Derecha: Ventilador en movimiento
------------------------------------------------------------------------------------------------------------------------------- 41

iv
ÍNDICE DE TABLAS

Tabla 1 Potencia instalada de energía eólica por lugares del mundo. Fuente: (BRITISH PETROLEUM,
2011) ................................................................................................................................................... 4
Tabla 2 Rangos y Resolución de instrumento de medición .............................................................. 12
Tabla 3 Factores de diseño y de ruido y sus respectivos niveles del Método de Taguchi ............... 20
Tabla 4 Arreglo ortogonal del diseño del experimento con los parámetros de diseño.................... 21
Tabla 5 Resultados de los experimentos para coeficiente de momento-par CT .............................. 25
Tabla 6 Tabla de Respuesta para Coeficiente de Momento-Par promedio para cada nivel de cada
variable de diseño del Método de Taguchi ....................................................................................... 26
Tabla 7 Resultados del experimento: Turbina con mayor coeficiente de momento-par ................. 29
Tabla 8 Constante de tiempo de respuesta a primer orden y aceleración angular para
configuraciones de 8 álabes. Ver en anexos las curvas de respuesta. .............................................. 33
Tabla 9 Comparación de variables cuantitativa entre la turbina diseñada por Murcia y la turbina
comercial ........................................................................................................................................... 35
Tabla 10 Anemómetro de hilo caliente. Nota: El anemómetro arroja datos de velocidad de aire a
nivel del mar, es decir con una densidad de aire de 1.225/ . Los valores arrojados por éste
fueron ajustados a la densidad del aire promedio de Bogotá 0.89/ . ..................................... 39
Tabla 11 Condiciones de diseño del ventilador axial extractor........................................................ 41
Tabla 12 Cálculo del momento de inercia para la configuración 8 álabes, con ángulo de ataque 50°
........................................................................................................................................................... 52
Tabla 13 Momentos de Inercia para configuraciones de turbina Uniandes probadas
experimentalmente y de turbina Comercial ..................................................................................... 53

v
ÍNDICE DE GRÁFICAS

Gráfica 1 Curva de calibración del motor: Velocidad angular en función del voltaje ....................... 14
Gráfica 2 Señal de voltaje del motor eléctrico en función del tiempo ............................................. 15
Gráfica 3 Velocidad angular vs. Tiempo ............................................................................................ 16
Gráfica 4 Caracterización de motor eléctrico para la turbina comercial ......................................... 23
Gráfica 5 Resultados según media de CT promedio por cada nivel y cada variable de diseño ....... 26
Gráfica 6 Potencia Mecánica promedio para cada nivel de cada variable de diseño según método
de Taguchi ......................................................................................................................................... 27
Gráfica 7 Coeficiente de Rendimiento Promedio para cada nivel de cada variable de diseño según
método de Taguchi............................................................................................................................ 28
Gráfica 8 Coeficiente de Rendimiento vs. Velocidad Específica de la turbina de mejor rendimiento.
 = 1.5 × 10 ................................................................................................................................ 29
Gráfica 9 Coeficiente de Momento-Par vs. Velocidad Específica de la turbina de mejor
rendimiento. Re = 1.5 × 10 ........................................................................................................... 30
Gráfica 10 Potencia Mecánica Rotacional vs. Velocidad Específica de la turbina de mejor
rendimiento.  = 1.5 × 10 ........................................................................................................... 30
Gráfica 11 Torque vs. Velocidad Específica de la turbina de mejor rendimiento.  = 1.5 × 10 30
Gráfica 12 Curvas de coeficiente de rendimiento Cp para configuraciones con bajo número de
álabes con ángulo de ataque de 70°.  = 1.5 × 10 ..................................................................... 31
Gráfica 13 Curvas de coeficiente de momento-par CT para configuraciones con bajo número de
álabes con ángulo de ataque 70°.  = 1.5 × 10 .......................................................................... 32
Gráfica 14 Comparación entre turbinas Uniandes y Comercial. Curvas Características a velocidad
de viento 5 m/s. ................................................................................................................................ 36
Gráfica 15 Comparación entre turbinas Uniandes y Comercial. Curvas Características a velocidad
de viento 10 m/s. .............................................................................................................................. 37
Gráfica 16 Curva de velocidad angular vs. Velocidad de viento ....................................................... 42
Gráfica 17 Curva de Caudal de extracción vs. Velocidad de viento .................................................. 43
Gráfica 18 Velocidad de aire extraído en el ducto de succión en función de la velocidad angular de
la turbina con y sin ventilador extractor. .......................................................................................... 44

vi
AGRADECIMIENTOS

Quiero agradecer a mi familia por apoyarme de manera incondicional a lo largo de

toda mi vida y por darme tantos momentos inolvidables llenos de cariño. Gracias
por todo su empeño y gran esfuerzo para darme la oportunidad de estudiar en
ésta universidad. A mi papá por darme siempre un ejemplo de responsabilidad y
honestidad en cada cosa que hago y por todo su interés en cada proyecto que
hice a lo largo de estos años. A mi mamá por darme todo su cariño siempre
mediante consejos y palabras que sólo una mamá sabe cómo y cuándo decir, que
me llenaron de fuerzas y muchas ganas de seguir adelante y por exigirme y
enseñarme a hacer siempre todas las cosas con cariño e impecables. A mi
hermana por cuidarme, por escucharme cuando lo necesité, por estar pendiente
de cómo me iba y por enseñarme cosas de la vida no ingenieril.

Quiero agradecer a mi asesor Álvaro Pinilla, por confiar en mí, por su paciencia y
por enseñarme a tener paciencia, por sus consejos, por su interés constante en
éste proyecto, por motivarme aun más y ayudarme a seguir estudiando y
aprendiendo en éste campo del aprovechamiento de recursos renovables de
energía dentro de mi proyecto de vida.

También, un especial agradecimiento a Omar Rodríguez, la persona que hace

que hasta lo que parece imposible)sea posible en el laboratorio de fluidos, por
toda su ayuda, por estar siempre disponible con todo su talento, por sus consejos
y buen humor, que hacen ameno el trabajo allí.

A Pao y Diana por ser unas excelentes compañeras de tesis, laboratorio y carrera.

Finalmente a todas las personas que me preguntaron en algún momento cómo iba
la tesis, que estuvieron pendientes y que me dieron una sonrisa. A todos los que
conocí a lo largo de la carrera, que me hicieron crecer como persona y como
ingeniera.

vii
RESUMEN

En el presente proyecto se desarrollaron pruebas experimentales a extractores de

aire accionados por rotores eólicos de eje vertical. Dichos turbo-extractores son
utilizados en la actualidad comúnmente para ventilar recintos industriales o
domésticos ahorrando un gran consumo de energía eléctrica. Sin embargo, poco
se conoce en cuanto a su diseño y rendimiento según diferentes variables como:
forma de los álabes, diámetro, condiciones de operación, comportamiento del flujo
a través de los álabes entre otros.

Las pruebas experimentales se realizaron en el túnel de viento cerrado TVIM-49-

60-1x1 de la Universidad de los Andes, Bogotá, Colombia.

En primer lugar, se evaluaron dos turbinas eólicas tipo Darrieus desde el punto de
vista de su capacidad de conversión de energía eólica en energía mecánica
rotacional. Una de ellas fue diseñada y manufacturada por Juan Pablo Murcia
(MURCIA, 2007) en su proyecto de grado llamada en éste proyecto como Turbina-
Uniandes y la otra es una turbina comercial.

Para la turbina Uniandes se diseñó un experimento con aplicación del método de

Taguchi para caracterizarla según tres variables de diseño: número de álabes,
ángulos de ataque y velocidades de viento. El diseño experimental permitió
conocer la influencia de cada una de éstas variables en el desempeño general de
la turbina en términos de coeficientes adimensionales de rendimiento y de
momento-par, potencia mecánica desarrollada y torque. A su vez se encontró la
turbina de mejor desempeño dada la mejor combinación de variables
mencionadas para maximizar los coeficientes de rendimiento y de momento-par.
Ésta se comparó con la turbina comercial mencionada de mayor diámetro.

En segundo lugar, se evaluaron las turbinas en cuanto a su capacidad de

extracción de aire en función de la velocidad de viento. Adicionalmente, se acopló
1
un ventilador extractor al eje de la turbina Uniandes para ver su influencia en el
caudal extraído. Se encontró que la turbina comercial no tiene un diseño
aerodinámico de sus álabes con respecto a las condiciones de operación. Dado lo
anterior, la turbina Uniandes, presenta un mejor rendimiento que la comercial, y su
caudal de extracción es competitivo en comparación con la comercial pese a tener
un menor tamaño.

Finalmente, el ventilador acoplado al eje de la turbina Uniandes, no tiene un efecto

positivo sobre el caudal de extracción disminuyéndolo. Es recomendable que de
ser instalado un ventilador en el ducto de la turbina, sea independiente de ésta. De
preferencia, el ventilador no debe depender de la velocidad rotacional de la turbina
ni de su potencia mecánica de entrada, pues las condiciones de operación
cambian según la velocidad de viento a la que opere la turbina, lo que hace
incluso difícil su diseño.

2
 1. INTRODUCCIÓN

Durante los últimos 15 años, la obtención de energía a partir de fuentes

renovables ha crecido significativamente y de manera exponencial. A finales de
2010 la potencia instalada de energía eólica en el mundo era cerca de 26 veces la
del año 1997 y en cuanto a energía solar en 2010 es cerca de 90 veces de lo que
se estaba generando 13 años atrás (BRITISH PETROLEUM, 2011).

Figura 1 Potencia instalada en energía eólica en el mundo Fuente: (BRITISH PETROLEUM, 2011)

Las fuentes renovables de energía aportan hoy un porcentaje importante en la

generación eléctrica de la mayoría de países europeos, Estados Unidos y Japón
principalmente. La mayoría de países en América Latina se están moviendo
también creando políticas e investigación en torno a las renovables, algunos
países como México, Brasil, Argentina y Costa Rica tienen varios proyectos
eólicos en funcionamiento.

En la siguiente tabla se puede observar la potencia instalada en energía eólica por

lugares del mundo para diciembre de 2010.

3
 Potencia Instalada
Lugar (Países que más aportan)
[MWe]

Norteamérica (U.S.A , Canadá) 44 285

Latinoamérica (Argentina, Brasil, México, Costa
Rica)
2 705
Europa (Alemania, España, Italia, Francia, Reino
Unido, Dinamarca)
87 674
África (Egipto, Morocco, Tunisia) 1 113
Asia (China, India, Japón) 63 645
Tabla 1 Potencia instalada de energía eólica por lugares del mundo. Fuente: (BRITISH PETROLEUM, 2011)

El viento como fuente de energía renovable, puede ser aprovechado gracias a

dos grandes grupos de rotores eólicos: Rotores de eje horizontal y rotores de eje
vertical. Los rotores de eje horizontal son comúnmente asociados con
aerogeneradores tripala ubicados cerca a las costas y/o en lugares con poca
edificación. También, las aerobombas ubicadas en sitios rurales para el bombeo
de agua pozo a animales o personas son accionadas por rotores de eje horizontal.
En contraposición, los rotores de eje vertical son menos reconocidos y populares.

Los rotores eólicos de eje vertical tienen aplicaciones, al igual que los de eje
horizontal, en generación eléctrica, bombeo, desalinización de agua, entre otros
(ISLAM, 2006), pero pueden llegar a ser más competitivos en aplicaciones como
extracción de aire de recintos, generación de electricidad urbana o mezclado de
sólidos o fluidos en sitios rurales. Lo anterior, debido a que estos rotores
presentan una ventaja, con respecto a los de eje horizontal, al no necesitar que el
viento tenga una dirección predominante a lo largo de un periodo de tiempo en
cierto lugar, esto abre grandes posibilidades de tener aplicaciones urbanas en
donde dada la existencia de edificaciones de alturas irregulares el comportamiento
del viento es arbitrario y no puede alcanzar altas velocidades.

Adicionalmente, los rotores de eje vertical presentan ventajas en cuanto a

simplicidad de operación pues no requieren caja de reducción y se facilita su

4
mantenimiento
o por su fácil acceso
acceso, pues son ubicados a bajas alturas o sobre
edificaciones.

Como ya se mencionó, una de las aplicaciones de los rotores de eje vertical es la

extracción de aire. Ésta puede significar un gran ahorro en consumo energético
tanto a nivel doméstico como a nivel industrial para ventilación de recintos,
recintos y son
ubicados en los techos de estas edificaciones
edificaciones. En las Figuras 2 y 3 se puede
observar una instalación urbana de turbo
turbo-extractores
extractores eólicos con rotor de eje
vertical.

Figura 2 Instalación de Extractores Eólicos en ambiente urbano. Fuente: (INDUSTRIAS GM, s.f.) Medellín,
Colombia.

Figura 3 Extractores Eólicos Empresa EXOVENT. Barranquilla, Colombia. Fuente: (EXOVENT, s.f.)

El viento incidente crea fuerzas de arrastre sobre los álabes de la turbina,

turbina creando
un momento par al eje. La potencia disponible de viento se convierte en potencia
mecánica rotacional y por
or el movimiento de los álabes se crea una presión de
succión que extrae el aire acumulado en el techo y lo expulsa al exterior a través

5
de los álabes. Al extraer el aire acumulado en los techos, se crean corrientes
naturales de aire que ventilan el recinto.

Figura 4 Creación de corrientes de aire y funcionamiento del extractor eólico. Fuente: (EXOVENT, s.f.)

Varias características determinan el rendimiento de éste tipo de turbinas eólicas, a

saber: diámetro (LAI, 2003), perfiles aerodinámicos de los álabes, velocidades de
viento, potencia mecánica al eje, el efecto de techos inclinados en el desempeño
de la turbina (LIEN & AHMED, 2010) y se presenta adicionalmente un fenómeno
de doble flujo o flujo cruzado que aún no ha sido completamente entendido, pues
dos fluidos se encuentran al tiempo en contacto con los álabes: el aire extraído y
el viento incidente.

En el presente trabajo, se realizaron diferentes pruebas experimentales a

extractores eólicos para entender su funcionamiento y la manera como afectan
diferentes variables de diseño y operación su desempeño general.

6
 2 . R E L A C I Ó N C O N T R A B A JO S PR E V I O S E N L A U N I V E R SI D A D D E
LOS AN DES

En el año 2003, Felipe Morales realizó una caracterización de un turboextractor

eólico comercial (MORALES A., 2003). Morales, diseñó una serie de
experimentos para caracterizar el sistema como: Su capacidad de transformación
de energía eólica en energía mecánica y su capacidad como sistema de
extracción de aire. Después, en su proyecto de grado Juan Pablo Murcia
(MURCIA, 2007) , modificó e implementó un método de diseño para turbinas
eólicas de eje vertical a bajos números de Reynolds basado en el Método de
Múltiples Tubos de Corriente. Una vez diseñada y manufacturada, Murcia realizó
más de 100 pruebas experimentales a la turbina para 11 diferentes
configuraciones posibles variando el número de álabes y 10 ángulos de ataque
posibles a una velocidad de viento promedio de 7.65 m/s.

Para las pruebas de las turbinas que caracterizaron Murcia (MURCIA, 2007) y
Morales (MORALES A., 2003), el suministro de viento se realizó con ayuda de un
ducto de 300mm de diámetro con un ventilador con una frecuencia máxima de giro
de 60Hz. Este montaje no fue completamente adecuado pues no podía brindar las
condiciones óptimas y suficientes para la caracterización de dichas turbinas en
cuanto a la calidad del flujo de viento, rango de velocidades de viento y tamaño del
ducto en comparación con el tamaño de las turbinas.

7
 3 . O B JE T I V O S

OBJETIVO GENERAL:

Caracterizar un extractor de aire eólico en el túnel de viento TVIM-49-60-1x1

OBJETIVOS ESPECÍFICOS:

- Mejorar estructuralmente la turbina para el montaje experimental en la

sección de pruebas del túnel de viento.
- Caracterizar la turbina eólica de eje vertical para diferentes configuraciones
según número de álabes, ángulos de ataque y velocidades de viento.
- Caracterizar y evaluar la capacidad de extracción de aire del turbo-extractor
- Visualizar con ayuda de humo las líneas de corriente de los fluidos a través
de la turbina.

8
 4. CARACTERIZACIÓN DE LA TURBINA EÓLICA

La turbina eólica diseñada y construida por (MURCIA, 2007) fue mejorada

estructuralmente y fue construido
construido, adicionalmente, un soporte para el montaje
experimental en la sección de pruebas del túnel de viento. El eje de diámetro de ¼
(in) de material acero plata
lata (AISI O1) se encuentra soportado por dos rodamientos
radiales en la parte inferior de la turbina. Los rodamientos se ubican en el soporte
de aluminio y una brida anclada al piso de la sección de pruebas soporta el
conjunto.

Figura 5 a) Turbina instalada en el soporte

oporte estructural
estructural,, se observa la brida y el rodamiento superior. b) Soporte en
aluminio. c) Rodamiento inferior

La turbina consta de 16 álabes con ángulos de ataque variables de 0° a 100° en

intervalos de 10°. Un total de 11 con
configuraciones
figuraciones son posibles según número de
álabes y su ubicación (MURCIA, 2007, pág. 15) y a su vez 10 ángulos de ataque
diferentes son posibles, para un total de 110 turbinas diferentes.

Un diagrama de velocidade
velocidades sobre el perfil del álabe de la turbina se muestra en la
el ángulo de ataque α de los álabes que se
Figura 6. Así mismo, se muestra e

9
quiere variar bajo diferentes condiciones de operación según el diseño
experimental.

Figura 6 Sección transversal turbina deje vertical y diagrama de velocidades. Tomado y Adaptado de (MURCIA, 2007,
pág. 11) y (PINILLA, Notas de Aerodinámica, 2010, págs. 8-3)

Adicionalmente, el rendimiento de cada turbina depende de la velocidad de viento

a la cual opera. Se desea encontrar la turbina de mejor configuración No. de
álabes y Ángulo de ataque que presente el mayor coeficiente de Momento-Par y
su velocidad de viento óptima de operación.

4.1. PROCEDIMIENTO EXPERIMENTAL

4.1.1. CURVAS CARACTERÍSTICAS DE ROTORES EÓLICOS

Existen tres coeficientes adimensionales que permiten cuantificar el desempeño

de un rotor eólico: El coeficiente de rendimiento CP que relaciona la potencia
mecánica rotacional al eje y la potencia disponible de viento, el coeficiente de
Momento-Par CT y la velocidad específica λ que relaciona la velocidad tangencial
de los álabes de la turbina y la velocidad no perturbada de viento. Con éstos

10
coeficientes es posible obtener las dos curvas adimensionales características de
los rotores eólicos: CT vs. λ y CP vs. λ.
Los coeficientes se pueden expresar como sigue (PINILLA, Notas del Curso
Electivo en Energía Eólica. Segundo semestre de 2011., 2011, págs. 3-21),

  
 =  = =
1 1  
2   2  
Donde P es la potencia mecánica al eje, ρ es la densidad del aire, V es la
velocidad no perturbada del viento, A es el área transversal proyectada de la
turbina, T es el momento-par, ω la velocidad angular del rotor y R el radio máximo
del extremo de los álabes al eje de rotación de la turbina.
Adicionalmente, con base en la relación,

 =∙

es posible encontrar una relación similar con los coeficientes adimensionales,

 =  ∙ 

Las curvas características se pueden obtener mediante un método indirecto a

partir del análisis de la respuesta dinámica de un rotor eólico. Con la respuesta en
velocidad angular del rotor en función del tiempo a una velocidad de viento de
entrada, se puede obtener la aceleración angular al derivar dicha velocidad
angular. Con el cálculo de la inercia y la aceleración angular, se obtiene el torque.
Con estos datos, y mediciones de radio, área transversal y velocidad de viento ya
es posible encontrar las curvas adimensionales.
Es importante mencionar que con el método indirecto, sólo se tienen mediciones
de torque de arranque y de velocidad angular de desboque. En contraste, el
método directo implica mediciones con un sensor de torque en diferentes puntos
de velocidad angular variando la carga.

11
 4.1.2. MONTAJE EXPERIMENTAL

Se realizó un montaje experimental que permitiese medir adecuadamente las

variables involucradas en los coeficientes adimensionales.

EQUIPOS E INSTRUMENTOS DE MEDICIÓN

- Túnel de viento TVIM – 49 - 60 - 1x1. Disponible en el Laboratorio de Dinámica de Fluidos

ML 032. Universidad de los Andes. Colombia.

- Manómetro digital:

VAISALA BAROMETRIC PRESSURE TRANSFER STANDARD PTB330TS

Izquierda: Sensor de humedad relativa y de temperatura. Derecha: Manómetro digital

Fabricante: VAISALA. Imagen tomada de 1

INSTRUMENTO RANGO RESOLUCIÓN

Manómetro digital 5000 hPa 1 Pa
Termopar -10°C a 40°C 0.01 °C
Sensor de humedad relativa Hasta antes de punto de condensación 0.01 %
Tabla 2 Rangos y Resolución de instrumento de medición

DENSIDAD DEL AIRE

La densidad del aire se determinó a partir de tres variables de entrada: presión

atmosférica, temperatura y humedad relativa. Los cálculos se basaron en el
procedimiento expuesto por P. Giacomo (GIACOMO, 1982). Así, la densidad se
puede expresar como

1
Manual del Usuario VAISALA. Disponible en Línea:
http://www.vaisala.com/Vaisala%20Documents/User%20Guides%20and%20Quick%20Ref%20Guides/PTB33
0TS_User_Guide_in_English.pdf Recuperado en Noviembre de 2011.

12
  !" !'
= $1 − &' (1 − )*
#  !"

Donde,
ρ : Densidad del aire (kg/m3)
!" : Masa molar del aire húmedo (kg mol-1)
!' : Masa molar del agua (kg mol-1)
&' : Fracción molar del vapor de agua
R: Constante de los gases ideales (R=8.1314 (J mol-1 K-1))
Z: Factor de compresibilidad
T: Temperatura del aire (K)

VELOCIDAD NO PERTURBADA DEL VIENTO

La velocidad del viento se obtuvo con mediciones de presión estática y absoluta

en un tubo de Pitot conectado al manómetro digital.

La velocidad se halló dada la siguiente relación,

2(, − -. // , : Presión absoluta [Pa]

=( ) -. : Presión estática [Pa]

: Densidad del aire instantánea [kg/m3]
V: Velocidad del viento [m/s]

A continuación, se observa el montaje experimental en la sección de pruebas del

túnel de viento para medición de velocidad de viento.

13
 Figura 7 Montaje experimental de medición de presión, humedad relativa y temperatura.

VELOCIDAD ANGULAR

Mediante un motor eléctrico acoplado al eje como generador, se adquirió una

señal de voltaje relacionada directamente con la velocidad angular del eje del
rotor. La curva de calibración del motor se muestra en la Gráfica 1. Así mismo el
motor acoplado al eje se muestra en la Figura 8.

Gráfica 1 Curva de calibración del motor: Velocidad angular en función del voltaje

14
 Figura 8 Motor eléctrico acoplado al eje del rotor

El motor fue conectado a un sistema de adquisición de datos con tiempos de

muestreo de 0.1s. Se adquirió la respuesta del rotor a partir del reposo hasta su
velocidad angular en estado estable. La respuesta del rotor se modeló como un
sistema de primer orden.

Gráfica 2 Señal de voltaje del motor eléctrico en función del tiempo

15
Así, con la conversión de voltaje a velocidad angular, y el modelo de respuesta a
sistema de primer orden es posible obtener una gráfica de velocidad angular en
función del tiempo,

Gráfica 3 Velocidad angular vs. Tiempo

ÁREA TRANSVERSAL Y RADIO

Se utilizó el software SolidEdge® para calcular el área proyectada de la turbina

dando como resultado un área de 0.03819 [m2]. El radio máximo de los álabes de
la turbina es 0,15 [m].

MOMENTO-PAR

Como se mencionó anteriormente, la aceleración angular se obtuvo derivando la

función del modelo de primer orden de la velocidad angular en función del tiempo.
El momento de inercia se determinó con el método del péndulo trifilar (Figura 9)
(MABIE & OCVIRK, págs. 321-322). Los datos obtenidos de inercias para cada
configuración estudiada se encuentran en el anexo A.

El momento-par en función del tiempo se obtiene gracias a la siguiente relación,

(01 = 2 ∙ 3(01

16
 Figura 9 Montaje experimental Péndulo Trifilar

NÚMERO DE REYNOLDS

El número de Reynolds Re es un parámetro adimensional que relaciona las

fuerzas
as inerciales sobre las fuerzas viscosas de un fluido, y está definido para éste
proyecto con base en la cuerda del perfil aerodinámico, para la turbina eólica de
eje vertical así,
4
=
5

Donde,
: Velocidad de viento (m/s)

4: Cuerda del perfil aerodinámico del álabe

: Densidad del aire (kg/m3)
5: Viscosidad dinámica (6 ∙ 7// )

La viscosidad dinámica del aire en Bogotá se estableció como

5 = 0,000017034 6 ∙ 7/ (RUIZ & VANEGAS, 2002, pág. 4)
4).

17
 4.1.3. DISEÑO DEL EXPERIMENTO PARA DETERMINAR LA TURBINA DE MEJOR
RENDIMIENTO: APLICACIÓN DEL MÉTODO DE TAGUCHI

El número de álabes, el ángulo de ataque y la velocidad de viento son las tres

variables que influyen en el desempeño de la turbina. Para encontrar la turbina de
mejor desempeño, se debe encontrar la influencia de cada una de dichas variables
en el coeficiente de momento-par. Como se mencionó anteriormente, según las 11
configuraciones posibles y los 10 posibles ángulos de ataque se deberían realizar
un total de 110 experimentos. Un total de 110 curvas de rendimiento serían
calculadas para una velocidad de viento específica. Esto quiere decir, que si se
desea saber la influencia de la velocidad del viento probando las turbinas a 2 o 3
velocidades diferentes, se deberían realizar 220 o 330 experimentos
respectivamente. La realización de tal número de experimentos resulta
inconveniente, exhaustivo y tomaría bastante tiempo llevarlos a cabo. El objetivo
de la caracterización es obtener la mejor turbina y la velocidad de viento óptima,
luego es necesario utilizar un método de diseño robusto de experimentos con
aplicación del método de Taguchi.
El método de Taguchi fue propuesto en Japón por el Dr. Genichi Taguchi en
1993. El método busca optimizar productos y procesos minimizando la variación
camino a mejorar la calidad. (DIETER, 2009). Se trata de un método de
optimización con base en la manipulación de variables de diseño y variables de
ruido, que en la aplicación de experimentos estadísticamente diseñados es
conocido como diseño robusto.

El diseño robusto permite encontrar los valores óptimos de cada una de las
variables que influyen en la optimización de cierto parámetro. En primer lugar, se
debe realizar un diseño de parámetros, experimentos estadísticamente diseñados
con base en un diseño factorial fraccionado (Fractional Factorial Design). Los
diseños tradicionales, no fraccionados, varían un parámetro a la vez para
encontrar su influencia. En comparación, los diseños fraccionados disminuyen el
número de experimentos que se deben llevar a cabo al realizar una combinación
ordenada de los parámetros de diseño y de ruido entre sí.
18
Se diseñaron entonces dos arreglos ortogonales para los parámetros de diseño y
para los parámetros de ruido.
Los parámetros de diseño son los factores de control que afectan la calidad del
desempeño de la turbina, a saber: No. de Álabes, Ángulo de Ataque y Velocidad
de viento.
Los factores de ruido se refieren a aquellas variables que influyen en el resultado
del desempeño de la turbina pero que son incontrolables ó no es práctico intentar
controlarlos. Para la turbina se ha identificado un factor de ruido de tipo externo: la
temperatura del aire de la sección de pruebas del túnel de viento. Se pudieron
medir dos niveles para este factor de ruido: Baja Temperatura (T<19°C, de 7 a
10am en la sección de pruebas) y Alta Temperatura (T>20°C, de 2 a 5pm en la
sección de pruebas). En el anexo 9.3 se muestra una adquisición de temperaturas
a lo largo de 24 horas en la sección de pruebas del túnel de viento para determinar
dichos horarios de medición.

Cada parámetro de diseño y de ruido debe tener ciertos niveles de operación. Se

han seleccionado 3 niveles para cada uno de los parámetros de diseño y dos
niveles para el parámetro de ruido.

Con base en los resultados obtenidos por (MURCIA, 2007), se seleccionaron 3

niveles con los que se obtuvieron buenos desempeños de la turbina en sus
pruebas experimentales: Las configuraciones de 8, 12 y 16 álabes fueron
seleccionadas. Así mismo, los ángulos de calaje de 50°, 70° y 90° fueron
seleccionados. Y los niveles de velocidad de viento a probar son 5m/s, 7.5m/s y
10 m/s. A continuación, se muestran imágenes de cada una de las turbinas
probadas.

19
 Figura 10 Configuraciones seleccionadas de la turbina evaluadas experimentalmente

En la Tabla 3 se muestran las variables de diseño y de ruido y sus respectivos

niveles a analizar. También, en la Tabla 4 se muestran los arreglos ortogonales
para los parámetros de diseño y de ruido.

Tabla 3 Factores de diseño y de ruido y sus respectivos niveles del Método de Taguchi

20
 Tabla 4 Arreglo ortogonal del diseño del experimento con los parámetros de diseño.

El arreglo ortogonal de los factores de diseño consta de 9 filas x 3 columnas, un

total de 9 experimentos deben ser llevados a cabo sólo variando dichos
parámetros. Adicionalmente, cada experimento debe ser realizado en cada nivel
del factor de ruido. En total, se desarrollaron 18 experimentos. A manera de
comparación, si se hubiera variado un parámetro a la vez, se habrían tenido que
ejecutar 3 = 27exp. × 2 niveles de ruido = 54 experimentos. El número de
experimentos se redujo a 1/3 con ayuda del método de Taguchi.

4.1.4. CURVAS ADIMENSIONALES PARA OTRAS CONFIGURACIONES DE LA TURBINA

Algunas configuraciones de la turbina ofrecen un tipo de rotores distintos de pocos

álabes y baja robustez. Se realizaron pruebas experimentales a 8 configuraciones
adicionales de la turbina, con un ángulo de ataque de 70° y velocidad de viento de
5m/s.
Las configuraciones de prueba se muestran en la Figura 11. El cálculo de los
momentos de inercia asociados a estas configuraciones mediante el método del
péndulo trifilar se encuentra en los anexos.

21
 Figura 11 Configuraciones de prueba para turbinas de pocos álabes (baja robustez)

4.1.5. CURVAS ADIMENSIONALES PARA UNA TURBINA COMERCIAL

Existen en el mercado una gran cantidad de turbo-extractores de diferentes

tamaños y materiales que son utilizados tanto a nivel doméstico como a nivel
industrial a pesar de no tener conocimiento suficiente acerca de cómo diseñarlos
de manera óptima para su función. Se desea comparar la turbina eólica diseñada
por (MURCIA, 2007) con una turbina comercial adquirida años atrás en Colombia
por la Universidad de los Andes y caracterizada por Morales (MORALES A., 2003)
en su proyecto de grado. Dada la forma de los álabes y su proceso de
manufactura, a partir de láminas de acero dobladas, se puede intuir que no existe
un diseño aerodinámico previo de los álabes ni de las condiciones de operación de
la turbina.
Se hallarán entonces las curvas características de la turbina comercial en el túnel
de viento TVIM-49-60-1x1, se compararán con las curvas reportadas por Morales

22
(MORALES A., 2003) en su proyecto y se compararán con las de la turbina
diseñada por Murcia.

Las características geométricas de la turbina comercial se muestran a

continuación,

Figura 12 Características geométricas del turbo extractor comercial. Fuente: (MORALES A., 2003, pág. 18)

Utilizando la misma metodología de caracterización con la otra turbina, se adquirió

una señal de voltaje mediante un motor eléctrico para obtener la respuesta en el
tiempo de la velocidad angular de la turbina. La caracterización del motor se
muestra en la siguiente gráfica,

Gráfica 4 Caracterización de motor eléctrico para la turbina comercial

23
El montaje experimental en el túnel de viento se observa en la Figura 12.

Figura 13 Montaje experimental en el túnel de viento para la turbina eólica comercial

Se realizaron mediciones del momento de inercia de la turbina, datos que se

encuentran registrados en los anexos. Con el momento de inercia y la aceleración
angular, fue posible encontrar el torque en función del tiempo y las curvas
adimensionales características.

Figura 14 Péndulo trifilar para medición de momento de inercia de la turbina comercial

24
 4.2. RESULTADOS Y ANÁLISIS

4.2.1. TURBINA DE MEJOR RENDIMIENTO: RESULTADOS DEL MÉTODO DE TAGUCHI

Se ejecutaron los 18 experimentos de la matriz de diseño de Taguchi y se

obtuvieron las curvas adimensionales características para cada uno de éstos.
Para encontrar cómo afectan las variables de diseño al coe
coeficiente
ficiente de momento-
momento
par y encontrar la combinación que lo maximiza se tabularon los datos de CT
máximo de cada experimento en la Tabla 5.

Tabla 5 Resultados de los experimentos para coeficiente de momento

momento-par CT

Se calculó una tabla de respuesta promedio por cada nivel para cada variable de
diseño.
Así, para el primer nivel de la variable Ángulo tenemos,

0.2896 X 0.3925 X 0.4064

6IJK 1: 50° NO,PÁRSTUV = = 0.363
0
3

Para el tercer nivel de la

a variable Velocidad de Viento tenemos,

0.4064 X 0.3913 X 0.2844

6IJK 3: 10/7 NO,PYZ[.Y\Z]^_ = = 0.361
3

25
El procedimiento se repite para cada uno de los niveles por cada variable de
diseño obteniendo los resultados presenta
presentados en la Tabla 6,

Tabla 6 Tabla de Respuesta para Coeficiente de Momento

Momento-Par
Par promedio para cada nivel de cada variable de diseño
del Método de Taguchi

Estos valores son graficados obteniendo la Gráfica 5.

Gráfica 5 Resultados según media de CT promedio por cada nivel y cada variable de diseño

De la Gráfica 5 se puede concluir lo siguiente,

Para la variable ‘Ángulo’,, el comportamiento es no lineal. El nivel de menor

variabilidad, es decir,
r, el que tiene una menor pendiente es el nivel 2. Éste nivel a
su vez maximiza el coeficiente de momento
momento-par.
par. En nivel óptimo de ángulo es el
nivel 2, álabes con ángulo de ataque a 70°.

La variable ‘No. De Álabes’ tiene un comportamiento lineal. Y el nivel

vel que máxima
el coeficiente momento-par
par es el nivel 3, 16 álabes.

26
Finalmente, la variable ‘Velocidad de viento’, presenta un comportamiento no
lineal. Aunque el coeficiente de momento-par permanece invariante para los
niveles 2 y 3, en el nivel 1 tiene un ligero aumento de un 5%.
Un análisis similar se puede realizar para otras variables de interés como potencia
mecánica (Gráfica 6) y coeficiente de rendimiento CP (Gráfica 7), obteniendo los
siguientes resultados:

Gráfica 6 Potencia Mecánica promedio para cada nivel de cada variable de diseño según método de Taguchi

De la Gráfica 6 se puede concluir que,

Para la variable ‘Ángulo’, el comportamiento tiende a ser invariante. Se observa

que el nivel 1-50° maximiza la potencia mecánica al eje pero de manera no muy
significativa. Se puede intuir que el ángulo no afecta la potencia mecánica al eje
alcanzada por la turbina.

La variable ‘No. De Álabes’ tiene un comportamiento lineal. Y el nivel que máxima

el coeficiente momento-par es el nivel 3, 16 álabes.

Finalmente, la variable ‘Velocidad de viento’, presenta un comportamiento lineal. A

mayor velocidad de viento mayor potencia mecánica al eje alcanza la turbina, esto

27
en un intervalo de 0 a 10 m/s. El nivel que maximiza la potencia mecánica
rotacional es el nivel 3, una velocidad de viento de 10 m/s.

Gráfica 7 Coeficiente de Rendimiento Promedio para cada nivel de cada variable de diseño según método de Taguchi

De la Gráfica 7 se puede concluir lo siguiente,

Para la variable ‘Ángulo’, el comportamiento no lineal aunque su variación entre

niveles es menor. Se observa que el nivel 2-70° maximiza el coeficiente de
rendimiento y a su vez garantiza el hecho de ser el nivel con más invariabilidad. El
nivel óptimo para maximizar el coeficiente de rendimiento es el nivel 2, 70°.

La variable ‘No. De Álabes’ tiene un comportamiento lineal, tal y como sucede con
el coeficiente de momento-par y la potencia mecánica. Y el nivel que máxima el
coeficiente de rendimiento es el nivel 3, 16 álabes.

Finalmente, la variable ‘Velocidad de viento’, presenta un comportamiento

invariante en general. Al igual que lo que sucede con el coeficiente de momento-
par, el coeficiente de rendimiento permanece invariante para los niveles 2 y 3. Sin
embargo, en la velocidad del nivel 1, 5 m/s, se observa un leve aumento.
Con base en el análisis anterior, es evidente que para las variables Ángulo y No.
De Álabes los niveles 2 y 3, respectivamente, son los que maximizan tanto el

28
coeficiente de momento-par como la potencia mecánica y el coeficiente de
rendimiento. La variable Velocidad de viento, permanece invariante para los 3
niveles en los dos coeficientes adimensionales, pero en ambos, se observa un
ligero aumento en el nivel 1, 5 m/s, que corresponde a su velocidad de viento de
diseño (MURCIA, 2007).
En resumen, los niveles óptimos para maximizar el coeficiente de momento-par y
de rendimiento de la turbina son:

Variable de Diseño Nivel Óptimo

Ángulo de ataque 70°
Número de Álabes 16
Velocidad de Viento 5 (m/s)
Tabla 7 Resultados del experimento: Turbina con mayor coeficiente de momento-par

La configuración óptima encontrada de ángulo y número de álabes coincide con la

reportada por Murcia en su caracterización (MURCIA, 2007).

Las curvas características para esta turbina a 5m/s de velocidad de viento se

muestran a continuación,

Gráfica 8 Coeficiente de Rendimiento vs. Velocidad Específica de la turbina de mejor rendimiento.

b (c⁄d1∗g,gbh (c1∗g,h ijk/cl m
`a = = n. bt × ngu
n,ohl p ngqb (rs∙d1

29
Gráfica 9 Coeficiente de Momento-Par vs. Velocidad Específica de la turbina de mejor rendimiento. vw = n. bt × ngu

Gráfica 10 Potencia Mecánica Rotacional vs. Velocidad Específica de la turbina de mejor rendimiento. `a = n. bt ×
ngu

Gráfica 11 Torque vs. Velocidad Específica de la turbina de mejor rendimiento. `a = n. bt × ngu

30
 4.2.2. CURVAS ADIMENSIONALES PARA OTRAS CONFIGURACIONES DE LA TURBINA

Las curvas adimensionales para las configuraciones de pocos álabes permiten

estudiar el rendimiento de las turbinas según la disposición de los álabes. Para
algunas configuraciones, los álabes se encuentran equidistantes mientras que
para otras no. Las curvas se muestran en las Gráficas 12 y 13.

Gráfica 12 Curvas de coeficiente de rendimiento Cp para configuraciones con bajo número de álabes con ángulo de
ataque de 70°. `a = n. bt × ngu

Las configuraciones con mejor rendimiento son la Conf. 10, la Conf. 8-1 y la Conf.
8-2. En la Gráfica 7, se concluyó que el coeficiente de rendimiento aumenta
cuando el número de álabes aumentan y esto se puede evidenciar en la Gráfica
12, en la cual el rendimiento aumenta conforme aumenta el número de álabes.

31
La turbina de 8 álabes ofrece 3 configuraciones distintas (Conf. 8-1, Conf. 8-2 y
Conf. 8-4) variando la manera en la que están dispuestos los álabes (Figura 11).
Cuando los álabes se encuentran equidistantes como en la configuración 8-1 la
turbina alcanza un mayor rendimiento, seguida de la turbina de la configuración 8
álabes agrupados de a dos (8-2), en donde a pesar de estar en parejas, la parejas
son equidistantes. Finalmente, cuando la turbina tiene 8 álabes agrupados de a
cuatro, tiene el menor rendimiento. Lo mismo sucede para éstas configuraciones
al obtener de mayor a menor coeficiente de momento-par respectivamente en la
Gráfica 13.

La misma tendencia se observa en las dos configuraciones de la turbina de 4

álabes. Cuando los álabes se encuentran equidistantes la turbina alcanza un
mayor rendimiento y un mayor coeficiente de momento-par que cuando están
agrupados de a dos.

Gráfica 13 Curvas de coeficiente de momento-par CT para configuraciones con bajo número de álabes con ángulo de
ataque 70°. `a = n. bt × ngu

32
Como se concluyó de la Gráfica 5, se obtiene un mayor coeficiente de momento-
par conforme aumenta el número de álabes de la turbina en general. Sin embargo,
en la Gráfica 13 se observa que la configuración 4-1 sobrepasa el CT de la
configuración 6-3 para velocidades específicas menores a 0.15.

La velocidad de viento, el área proyectada de la turbina, el radio y la densidad del

viento son constantes para las gráficas 12 y 13. Así, la única variable que afecta el
coeficiente de momento-par es el torque. El torque depende de la inercia y de la
aceleración angular del rotor. Al comparar únicamente las 3 configuraciones de 8
álabes, la inercia no varía significativamente pues la masa se conserva. Así, la
aceleración angular juega el papel más importante. Ésta es la derivada de la
respuesta en el tiempo de la velocidad angular que fue modelada como una
respuesta de primer orden, por lo cual la aceleración depende de la constante de
tiempo τ.

Constante de Aceleración
Turbina
tiempo τ (s) Angular (rad/s2)
Conf. 8-1 3,10 5,93
Conf. 8-2 3,30 5,00
Conf. 8-4 3,55 4,78
Tabla 8 Constante de tiempo de respuesta a primer orden y aceleración angular para configuraciones de 8 álabes.

En la tabla 8 se puede observar que a mayor τ, menor aceleración angular, y

como consecuencia menor torque, lo cual explica el menor coeficiente de
momento-par y de rendimiento de las configuraciones 8-4 y 8-2.

La disposición de los álabes influye directamente en el tiempo de respuesta del

rotor hasta llegar a la velocidad en estado estable. Dos situaciones de bloqueo en
el arranque de la turbina pueden presentarse cuando el viento incide sobre el rotor
como se muestra en la Figura 15. Al agrupar los álabes se dejan espacios por los
cuales el viento tiene la posibilidad de pasar y dificultar el movimiento rotacional.

33
 Figura 15 Esquema vista superior de d
disposición de álabes para configuraciones de turbina de pocos álabes y
situaciones de bloqueo con respecto a la dirección del viento incidente.

En la situación de Bloqueo 1, una corriente de viento atraviesa internamente la

turbina, dificultando el movimiento de los álabes y como consecuencia el arranque.
En la situación de Bloqueo 2, corrientes de aire rodean la turbina pasando a través
de gran parte de su circunferencia bloqueando el movimiento. A medida que la
turbina empiece a rotar cada vez a mayor velocidad
velocidad,, menos cantidad de aire podrá
atravesar la circunferencia de la turbina y ésta logrará mantenerse en movimiento.
m
Así mismo, según la disposición de los álabes, entre menos espacio exista por el
cual el aire pueda atravesar la turbina, mayor velocidad angular estable tendrá la
turbina a cierta velocidad de viento. Estas situaciones de bloqueo son las que
extienden
tienden el tiempo de respuesta del rotor aumentando su constante de tiempo y
presentando una aceleración angular menor hasta llegar a su velocidad angular
estable.

4.2.3. CURVAS ADIMENSIONALES PARA U

UNA TURBINA COMERCIAL

En la Tabla 9 se muestra una comparación entre las dos turbinas para algunas
variables de interés,

34
 TABLA DE COMPARACIÓN

TURBINA

Uniandes Comercial
Conf. 16-70°
2
Inercia (kg m ) 0,00613 0,06476
Masa (kg) 0,83 2,38
No. Álabes 16 30
Material de los álabes ABS Aluminio
Diámetro (m) 0,3 0,43
Altura de los álabes (m) 0,15 0,259
2
Área Frontal (m ) 0,03819 0,08838
Tabla 9 Comparación de variables cuantitativa entre la turbina diseñada por Murcia y la turbina comercial

Se realizaron mediciones a una velocidad de viento de 5 m/s. La comparación

entre las curvas características de las dos turbinas se muestra en la Gráfica 14.
Se observa en general que aunque dimensionalmente la turbina Uniandes es más
pequeña, el coeficiente de rendimiento supera el alcanzado por la turbina
comercial bajo la misma velocidad de viento. En la Gráfica 14 se muestra
adicionalmente el coeficiente de momento
momento-par
par para las dos turbinas en función de
la velocidad específica, teniendo un may
mayor coeficiente la turbina Uniandes.
Uniandes

En general que la velocidad angular en estado estable de la turbina Comercial es

mucho menor que la alcanzada por la turbina Uniandes.. Una característica que
afecta directamente la velocidad angular es el material de los álabes y la masa
que le agregan al sistema, aumentando la inercia y por consecuencia el torque de
arranque. Así mismo el valor de la constante de tiempo de respuesta del sistema
de primer orden a una entrada escalón es mayor, y la respuesta en el tiempo
tiem es
más lenta. La turbina Uniandes tiene un menor torque de arranque que la turbina
Comercial.

35
 Gráfica 14 Comparación entre turbinas Uniandes y Comercial. Curvas Características a velocidad de viento 5 m/s.

A una velocidad de viento de 5 m/s las potencias máximas alcanzadas por las
turbinas son iguales cercanas a 0,17W. La potencia mecánica aumenta si la
velocidad de viento aumenta, así en las curvas para una velocidad de 10 m/s

36
(Gráfica 15) se puede observar que la turbina Comercial alcanza una potencia de
1.2W y supera a la de la turbina Uniandes (0,9W).

Gráfica 15 Comparación entre turbinas Uniandes y Comercial. Curvas Características a velocidad de viento 10 m/s.

El momento-par generado en la turbina Comercial es mayor pues su área

transversal y su masa es mayor, así como su inercia. Dado que este torque es
mayor, aunque la velocidad angular sea menor, esto puede explicar que su
potencia mecánica iguale y hasta sobrepase a la desarrollada por la turbina
Uniandes aunque su rendimiento siga siendo menor.

5. DETERMINACIÓN DE LA CAPACIDAD DE EXTRACCIÓN

En este capítulo se muestra un procedimiento experimental para encontrar la

capacidad de extracción de aire de las dos turbinas estudiadas en este proyecto.
Para esto se construyó un Plenum como volumen de control y se midió el caudal
de aire extraído a través del ducto de las turbinas para diferentes velocidades de
viento. Adicionalmente, se diseño y acopló al eje de la turbina Uniandes un
ventilador extractor y se evaluó el caudal extraído con éste.

37
 5.1. PROCEDIMIENTO EXPERIMENTAL

5.1.1. CAUDAL DE EXTRACCIÓN

A través de la construcción de un volumen de control se puede medir la velocidad

del aire extraído a través de un ducto de las turbinas y el caudal de aire extraído
se puede calcular como,
x = Vz{|} ∙ ~z{€z
Donde,

x: Caudal de aire extraído ( ⁄71

-"‚ƒ" : Área transversal del ducto de salida ( 1
"‚O„ : Velocidad del aire a través del ducto de salida (⁄71

5.1.2. MONTAJE EXPERIMENTAL

Se construyó un volumen de control en forma de prisma rectangular con un área

inferior de entrada de aire y otra en la parte superior conectada al ducto de las
turbinas como se muestra en el esquema de la Figura 16,

Figura 16 Esquema de montaje experimental en el túnel de viento para medición de la velocidad de extracción

38
El plenum se compone de un marco en madera y de paredes cubiertas con
plástico. Tiene de sección cuadrada de dimensiones 600 x 600 x 800 (mm). El
área de entrada es de 142 x 142 (mm) igual al área de salida al ducto en la parte
superior es ‡ ƒˆ‰., = (0.16 1 = 0.02 .
† †

La medición de la velocidad se realizó mediante un anemómetro de hilo caliente

para medición de velocidades bajas (menores a 1 m/s) cuyos datos técnicos se
encuentran en la Tabla 10. La sonda del anemómetro fue ubicada en la sección
transversal del ducto como se muestra en la Figura 17, y se midió la velocidad a
diferentes distancias desde la pared del ducto hasta el centro para obtener una
distribución de velocidades a través del mismo.

Anemómetro de Hilo Caliente

(Hot Wire Thermo-Anemometer)

Imagen obtenida de 2
Fabricante EXTECH Instruments
Temperatura de Operación 0 a 50°C
Rango 0.2 a 20.0 m/s
Resolución 0.1 m/s
Tabla 10 Anemómetro de hilo caliente. Nota: El anemómetro arroja datos de velocidad de aire a nivel del mar, es
decir con una densidad de aire de n. ŠŠb jk⁄cl . Los valores arrojados por éste fueron ajustados a la densidad del
aire promedio de Bogotá g. ‹h jk⁄cl.

En montaje experimental en la sección de pruebas del túnel de viento se muestra

en la siguiente figura,

2
Heavy Duty Hot Wire Thermo-Anemometer: Airflow meter with telescoping probe designed to fit into
small openings. Disponible en: http://www.extech.com/instruments/product.asp?catid=1&prodid=39

39
 Figura 17 Arriba: Montaje experimental para medición de caudal de aire extraído en el túnel de viento. Abajo: Vista
inferior al interior del ducto, se observa el sensor de velocidad de aire.

5.1.3. VENTILADOR EXTRACTOR

Se diseñó un ventilador axial con base en el Análisis simplificado de diseño de

ventilador/bomba axial (PINILLA, Notas de Aerodinámica, 2010, págs. 8-1 a 8-8).

Las características de diseño del ventilador se muestran en la siguiente tabla,

40
 CONDICIONES DE DISEÑO DEL VENTILADOR
DIÁMETRO (mm) 140
PRESIÓN (Pa) 0,05
VEL. ROT.(RPM) 191
No. DE ASPAS 3
VELOCIDAD ESPECÍFICA NS 10,15
DIÁMETRO ESPECÍFICO ∆ 0,4
PERFIL PLACA PLANA
DIÁMETRO DEL CUBO (mm) 25,2
Tabla 11 Condiciones de diseño del ventilador axial extractor

El diámetro del ventilador está limitado dimensionalmente por el diámetro del

ducto de succión de la turbina el cual es de 160mm. Se seleccionó un diámetro de
140mm tomando en cuenta defectos dimensi
dimensionales
onales de manufactura tanto del
ducto como del ventilador.

La velocidad angular de operación corresponde a la velocidad en el punto de

mejor rendimiento de la turbina para una velocidad de viento de 10 m/s (Ver anexo
9.2).
Las aspas fueron moldeadas en lámina metálica de ferrosilicio y el cubo fue
manufacturado en aluminio y fue acoplado al eje de la turbina (Figura 18).
18)

Figura 18 Izquierda: Ventilador acoplado al eje de la turbina. Derecha: Ventilador en movimiento

41
 5.2. RESULTADOS Y ANÁLISIS DE RESULTADOS

Las curvas presentadas a continuación corresponden a mediciones realizadas a la

turbina Comercial y a la mejor configuración de la turbina Uniandes (Conf. 16
álabes – 70°).

5.2.1. CURVA DE VELOCIDAD ANGULAR VS. VELOCIDAD DE VIENTO

La turbina Uniandes tiene una velocidad de arranque menor (0.8 m/s) que la
turbina Comercial (3.4 m/s). Así mismo llega a velocidades angulares mayores
para la misma velocidad de viento como se observa en la Gráfica 16.

Gráfica 16 Curva de velocidad angular vs. Velocidad de viento

La masa de la turbina Uniandes es casi tres veces menor a la masa de la turbina

Comercial. La masa influye en la velocidad angular alcanzada por la turbina al
modificar su inercia.

42
 5.2.2. CURVA DE CAUDAL DE EXTRACCIÓN DE AIRE VS. VELOCIDAD DE VIENTO

Gráfica 17 Curva de Caudal de extracción vs. Velocidad de viento

La Gráfica 17 muestra el caudal de aire extraído en función de la velocidad de

viento. Pese a su tamaño, la turbina Uniandes es competitiva con respecto a la
Comercial. La reducción en tamaño y volumen de la turbina Uniandes se
compensa al alcanzar velocidades angulares mayores. Adicionalmente, cabe
anotar que la turbina Comercial tiene una velocidad de arranque de 3.4 m/s
(Gráfica 16), y en la Grafica 17 se observa que extrae aire incluso cuando su
velocidad angular es cero a 2.1 m/s de viento. Esto se debe a que la velocidad de
viento crea un diferencial de presión suficiente entre el exterior y el plenum para
inducir el movimiento del aire desde el interior del plenum hacia el exterior.

5.2.3. CURVAS DE VELOCIDAD DE EXTRACCIÓN DE AIRE VS. VELOCIDAD ANGULAR

CON Y SIN VENTILADOR

Los resultados de las pruebas se muestran en la Gráfica 18. El ventilador influye

negativamente en la capacidad de extracción de la turbina disminuyendo la
velocidad del aire en el ducto de succión.

43
Gráfica 18 Velocidad de aire extraído en el ducto de succión en función de la velocidad angular de la turbina con y sin
ventilador extractor.

Dado lo anterior, se puede concluir que la turbina extrae un buen caudal de aire
sólo con la presión de succión que crea en el ducto sin necesidad de acoplar un
ventilador directamente al eje de la turbina. La velocidad angular de la turbina
cambia en función de la velocidad de viento que depende del entorno en el que se
encuentre la turbina en un rango considerable. Este cambio en el punto de
operación de velocidad angular del ventilador dificulta su diseño y para ciertas
velocidades de viento el ventilador presenta un efecto negativo sobre el caudal de
extracción. (Lai, 2005) presenta un trabajo acerca de la mejora en la capacidad de
un turbo-extractor accionado por la turbina de eje vertical y el accionamiento de un
ventilador en el ducto mediante un panel solar.

44
 6 . C O N C L U SI O N E S

Respecto a la caracterización de la turbina eólica en cuanto a su capacidad de

conversión de energía eólica potencial en energía mecánica rotacional:

- La aplicación del método de Taguchi al diseño de experimentos permitió

diseñar un experimento factorial fraccionado que redujo el número de
experimentos a llevar a cabo en dos tercios. Así mismo, se pudo estudiar la
influencia de tres variables de diseño como son: No. de álabes, ángulo de
ataque y velocidad de viento bajo diferentes niveles de prueba, sobre el
desempeño de la turbina eólica.

- Se encontró la configuración de número de álabes y ángulo de ataque que

maximiza los coeficientes adimensionales de la turbina Uniandes:
Configuración 16 álabes con ángulo de ataque a 70°. Los resultados
coinciden cualitativamente con lo reportado por Murcia (MURCIA, 2007).

- Cuantitativamente los resultados de las curvas características obtenidas en

el presente trabajo son diferentes a las reportadas por Murcia (MURCIA,
2007). La velocidad angular estable para cierta velocidad de viento
disminuyó a la mitad. Así mismo la turbina Comercial disminuyó a la mitad
su velocidad angular con respecto a lo reportado por Morales (MORALES
A., 2003). Un factor que pudo haber influido en estos resultados es el
montaje experimental para la caracterización. Pues en el túnel de viento
TVIM-49-60-1x1 de la Universidad de los Andes, se cuenta con un flujo de
viento controlado, ordenado y con variaciones de presión mínimas.
Adicionalmente, el viento rodea a la turbina completamente en toda su
periferia incluyendo el ducto. En el montaje de Murcia (MURCIA, 2007) y
Morales (MORALES A., 2003), el ducto disponible era abierto y del tamaño
de las turbinas, de tal manera que éstas eran impactadas por un chorro de
viento de su tamaño. Así mismo, para el presente trabajo se realizó

45
 medición de densidad de aire instantánea para cada prueba y se midieron
todas las inercias asociadas a las configuraciones de prueba de las
turbinas.
- Las turbinas con mayor número de álabes presentan un mejor desempeño.

- A mayor velocidad de viento, la potencia mecánica desarrollada por las

turbinas es mayor. El límite entonces está relacionado con la resistencia
mecánica y estructural de la turbina. Sin embargo, la probabilidad de que
estos extractores eólicos operen a velocidades de viento de 10 m/s o mayor
es mínima. Realizar mediciones a velocidades de viento mayores sin
innecesarias.

- Las turbinas con álabes equidistantes presentan un mejor rendimiento que

las turbinas con álabes agrupados. Dicha disposición de los álabes afecta
directamente el tiempo de respuesta del rotor a cierta velocidad de viento y
como consecuencia su aceleración angular disminuye y consecuentemente
el torque y la potencia mecánica desarrollada.

- El material de los álabes modifica la inercia de la turbina, el torque y la

velocidad angular estable de la turbina modificando sus curvas de
rendimiento.

- Los coeficientes adimensionales, momento-par y potencia mecánica

desarrollados por la turbina Uniandes son mejores que los desarrollados
por la Comercial. Lo que confirma que la turbina Comercial no cuenta con
un diseño aerodinámico de fondo a diferencia de la Turbina Uniandes. A su
vez, la turbina Comercial presenta un mayor momento-par de arranque, lo
que no es deseable para operar a velocidades bajas de viento.

46
Respecto a la capacidad de extracción:

- La capacidad de extracción de aire de la turbina Uniandes es similar a la

de la turbina Comercial a diferentes velocidades de viento pese a ser de
menor tamaño lo que la hace competitiva.

- La turbina Uniandes no requiere de ventilador acoplado directamente al eje

para extraer mayor cantidad de aire. Incluso éste puede actuar
negativamente sobre su desempeño. Si, según la aplicación, el usuario
requiere de mayor caudal de extracción de aire cuando hay bajas o nulas
velocidades de viento puede accionar un ventilador extractor conectado a
cualquier fuente de electricidad, completamente independiente de la
turbina. Lo importante, es que éste ventilador pueda ser prendido o
apagado cuando el usuario lo desee, y no dependa de la potencia mecánica
de la turbina ni de su velocidad angular. Pues dependiendo de la velocidad
de viento y velocidad angular de la turbina éste puede actuar tanto de
manera neutral como negativamente sobre el caudal de extracción.

- A una velocidad de viento promedio de 5 m/s que se puede encontrar en las

ciudades, los turboextractores pueden extraer del orden de 40 l/s.

47
 7 . R E CO M E N D A CIO N E S Y T RA BA J O F U T U RO

VISUALIZACIÓN DE LAS LÍNEAS DE CORRIENTE A TRAVÉS DE LA TURBINA

Durante el desarrollo del proyecto se logró realizar una visualización con ayuda de
una máquina generadora de humo para observar el comportamiento del fluido de
aire extraído para ciertas
rtas velocidades de viento.

En la siguiente ilustración se observan algunos fotogramas de un video para una

velocidad de viento de 1 m/s. La cámara tomó las fotografías con una frecuencia
de 30fps. Se puede observar que la turbina extrae el aire y lo expulsa
e en
paquetes o bocanadas de aire.

Ilustración 1 Visualización de flujo con ayuda de humo con una velocidad de viento de 1 m/s. La velocidad angular de
la turbina es 38 RPM valor promedio.

El caudal arrojado por la cámara humo y la presión no fueron fáciles de controlar.

Por lo cual, se observó que éste caudal modificó la velocidad angular de la turbina
aumentándola en comparación con la velocidad angular de la turbina esperada y

48
caracterizada anteriormente para dicha velocidad de viento. Sin embargo, esta
velocidad angular modificada por el caudal de humo se mantuvo constante.

El flujo de viento en el túnel de viento a velocidades bajas de operación como la

mostrada en la ilustración no es completamente uniforme (RAMIREZ, 2010), por
esto la visualización pudo ser afectada por la variabilidad de flujo en la sección de
pruebas del túnel.

Se recomienda, para visualizar de una manera más definida las líneas de corriente
a diferentes velocidades de viento (incluyendo 5m/s, que corresponde a su
velocidad de diseño), utilizar la técnica del cable de humo (Smoke-wire) de Batill
Mueller 1983.

CAPACIDAD DE EXTRACCIÓN

En cuanto a la capacidad de extracción sería conveniente introducir una variable

más al experimento controlando la temperatura dentro del plenum. Así es posible
cuantificar el aporte de extracción de aire mediante flujo libre de aire hacia el
exterior debido al diferencial de temperaturas.

49
 8 . B I B L I O G R A FÍ A

BRITISH PETROLEUM. (2011). Statistical Review of World Energy 2011. Recuperado el 16 de 10 de

2011, de www.bp.com

DIETER, G. (2009). Engineering design (4ta ed.). Boston: McGraw-Hill.

EXOVENT. (s.f.). Sistema de Ventilación Eólico. Recuperado el Mayo de 2011, de EXOVENT:

http://www.exovent.com/sistema.htm

GIACOMO, P. (1982). Equation for the Determination of the Density of Moist Air. Metrología 18 ,
33-40.

INDUSTRIAS GM. (s.f.). Extractor de Turbina E-2100. Recuperado el 04 de 2011, de INDUSTRIAS

GM: http://www.industriasgm.com/

ISLAM, M. (2006). Aerodynamic models for Darrieus-type straight-bladed vertical axis wind
turbines. ELSEVIER , Renewable and Sustainable Energy Reviews 12 (2008)1087-1109, 23.

LAI, C.-m. (2003). Experiments on the ventilation efficiency of turbine ventilators used for buildings
and factory ventilation. ELSEVIER , ENERGY AND BUILDINGS , 6.

LIEN, S., & AHMED, N. (2010). Effect of Inclined roof on the airflow associated with a wind driven
turbine ventilator. ELSEVIER. , Energy and Buildings 43 (2011) 358-365, 8.

MABIE, H., & OCVIRK, F. (1973). Mechanical and Dynamics of Machinery. New York: Wiley.

MORALES A., F. O. (2003). CARACTERIZACIÓN DE UN TURBOEXTRACTOR EÓLICO. Bogotá,

Colombia: Proyecto de Grado Ingeniería Mecánia - Universidad de los Andes.

MURCIA, J. P. (2007). DISEÑO, CONSTRUCCIÓN Y CARACTERIZACIÓN DE UN EXTRACTOR

ACCIONADO POR UN ROTOR EÓLICO DE EJE VERTICAL. Bogotá, Colombia.: Proyecto de Grado
Ingeniería Mecánica - Universidad de los Andes.

MURCIA, J. P., & PINILLA, Á. (2007). DESIGN AND LABORATORY EXPERIMENTAL TESTS OF A
TURBINE VENTILATOR. ASME International Mechanical Engineering Congress and Exposition , 5.

PINILLA, Á. (2010). Notas de Aerodinámica. Bogotá: Departamento de Ingeniería Mecánica -

Universidad de los Andes.

PINILLA, Á. (2011). Notas del Curso Electivo en Energía Eólica. Segundo semestre de 2011. Bogotá:
Universidad de los Andes.

RAMIREZ, M. (2010). Comportamiento del flujo de aire en el túnel de viento TVIM-49-60-1x1.

Bogotá: Proyecto de Grado Ingeniero Mecánico.
50
ROSS, I. (2011). Wind Tunnel blockage corrections: Review and application to Savonius Vertical-
axis Wind Turbines. . Journal of Wind energy and Industrial Aerodynamics 99(2011) 523-538 , 16.

RUIZ, M. A., & VANEGAS, S. (Mayo de 2002). Modelamiento mediante la dinámica de fluidos
computacional de un chorro plano libre. Obtenido de MEF. Simulación y Análisis Ltda.:
http://www.mefsimulacion.com/pdf/General/CFDJet.PDF

TAYLOR, J. R. (1997). An introduction to error analysis : the study of uncertainties in physical

measurements (Segunda ed.). Sausalito, Calif. : University Science Books, c1997. .

WHITE, F. (2006). Mecánica de Fluidos (6ta ed.). Madrid, España: McGrawGill.

9 . ANEXOS

51
ANEXO 9.1: MEDICIONES DE INERCIA PARA DISTINTAS TURBINAS PROBADAS

EXPERIMENTALMENTE

Para cada configuración se midió con un cronómetro en el péndulo trifilar el tiempo en dar 10
oscilaciones. Se halló un promedio del periodo de oscilación y se propagó la incertidumbre
asociada a cada variable en la siguiente ecuación para el cálculo de la inercia (MABIE & OCVIRK,
1973):

Se siguió el procedimiento mostrado a continuación con todas las demás configuraciones cuyos
datos se encuentran tabulados en la tabla 13:

CONFIGURACIÓN 8- 50
No Alabes 8
Ángulo 50
masa (kg) 0.6394
L (m) 0.49
Porcentaje de
R (m) 0.1 Error
Inercia J 0.00383 kg m2 ± 0.00017 4%

N° Osc t (s) T (s) Propagación de Incert.

10 10.89 1.09 Ut 0.02 s
10 10.87 1.09 Um 0.0001 kg
10 10.82 1.08 UT 0.02 s
10 10.89 1.09 UR 1.00E-03 m
10 10.83 1.08 UL 5.00E-03 m
10 10.87 1.09 dT/dJ 0.0070 s
10 10.9 1.09 dm/dJ 0.0060 kg
10 10.87 1.09 dR/dJ 0.0766 m
10 10.85 1.09 dL/dJ -0.0078 m
10 10.84 1.08 UJ 0.00017 kg m2
10 10.96 1.10
10 10.84 1.08
10 10.8 1.08
10 10.88 1.09
PROMEDIO 1.09
Desviación 0.0040
t0.05,n-1 2.1604
IC 95% Error ± (s) 0.0023
Tabla 12 Cálculo del momento de inercia para la configuración 8 álabes, con ángulo de ataque 50°

52
Tabla 13 Momentos de Inercia para configuraciones de turbina Uniandes probadas experimentalmente y de turbina
Comercial

ANEXO 9.2: CURVAS CARACTERÍSTICAS TURBINA UNIANDES PARA VELOCIDAD DE VIENTO

VI 10 M/S.

53
54
ANEXO 9.3: VARIACIÓN DE LA TEMPERATURA EN LA SECCIÓN DE PRUEBAS DEL TÚNEL DE VIENTO
A LO LARGO DE 24 HORAS.

55