See discussions, stats, and author profiles for this publication at: https://www.researchgate.

net/publication/292293548

Diseño del Rotor de una Turbina Eólica con Alto Rendimiento Aerodinámico
(Numerical design of a wind turbine rotor blades with a high aerodynamic
performance)

Conference Paper · October 2006

CITATIONS READS

0 3,608

5 authors, including:

Nelson Garcia-Polanco Pedro Guillermo Pieretti

Whirlpool EMEA Simon Bolívar University
63 PUBLICATIONS 64 CITATIONS 6 PUBLICATIONS 7 CITATIONS

SEE PROFILE SEE PROFILE

All content following this page was uploaded by Nelson Garcia-Polanco on 30 January 2016.

The user has requested enhancement of the downloaded file.

 1

Diseño del Rotor de una Turbina Eólica con Alto Rendimiento

Aerodinámico

A. Graterol1, G. Mayoral1, N. García1, H. Díaz2, P. Pieretti2

Universidad Simón Bolívar
Caracas Apdo. 89000, Venezuela
1
Laboratorio de Mecánica de Fluidos
2
Instituto de Energía-INDENE
e-mail: anibal_graterol@yahoo.es, paragerald@gmail.com, ngarcia@usb.ve

Resumen

El presente trabajo de investigación tiene como objetivo el diseño aerodinámico del rotor de una turbina eólica de 1000 W,
que permita generar electricidad de forma rentable en las zonas rurales del país que no cuenten con un servicio eléctrico
adecuado. Para ello el diseño se basará en zonas con vientos clase 3 (velocidad promedio entre 5.1 y 5.6 m/s), la cual
coincide con las velocidades de viento presentes en muchas zonas de Venezuela. Para realizar el diseño se estudiará en
primer lugar cuáles son las variables que intervienen en el comportamiento aerodinámico de una turbina eólica y se
determinará de qué manera puede realizarse un diseño que permita una alta eficiencia aerodinámica con vientos clase 3; esto
se hará en parte mediante el uso de un algoritmo de cálculo (basado en la teoría de momento y la teoría de elemento de pala)
que fue adaptado a las condiciones estudiadas y que estima el rendimiento de la turbina en estado estable. A partir del
diseño aerodinámico teórico se construirá un modelo a escala 1:2 que permita estudiar el comportamiento real de la turbina
y poder determinar los requerimientos y dificultades asociadas a la construcción, a partir de los cuales se pueda recomendar
los materiales a utilizar y un proceso de fabricación de los álabes que permitan una construcción sencilla y económica sin
detrimento de la calidad del producto final.

Palabras claves: Diseño Aerodinámico de Hélices, Turbina Eólica, Teoría del Elemento de Pala.

Introducción en aerogeneradores eólicos en el país, dado que ya

prácticamente se ha llegado al tope de energía hidroeléctrica
Últimamente la sociedad mundial ha ido cobrando aprovechable y el ahorro de combustibles fósiles cobra
conciencia sobre el gran impacto que tienen las actividades mucha importancia ante los altos precios actuales del
humanas en la naturaleza y se están dando los primeros petróleo en el mercado internacional.
pasos hacia el desarrollo sustentable. En este ámbito la El aprovechamiento del recurso eólico puede realizarse
energía eólica ha presentado un desarrollo muy importante: desde dos perspectivas: 1) desarrollando grandes plantas
el año 2005 concluyó con una capacidad instalada mundial eólicas en las zonas con mayor potencial eólico del país,
de 52,322 MW, lo cual representa un crecimiento del 25% como lo son la Península de la Guajira y la Península de
respecto al 2004 [1] y se espera que para los próximos años Paraguaná, lo cual implicaría el uso de turbinas de gran
el crecimiento del sector se mantenga por encima del 15% tamaño diseñadas para velocidades de viento altas, las
anual. Esto ha sido posible gracias a las investigaciones que cuales tienen altas eficiencias y 2) desarrollando
se han venido realizando a nivel mundial en el área de aerogeneradores destinados a cubrir las necesidades de
generación eólica, las cuales han permitido disminuir pequeñas comunidades e incluso viviendas individuales en
considerablemente los costos de producción de energía zonas con servicio eléctrico nulo o deficiente, esto
eléctrica a través de aerogeneradores. En la actualidad implicaría el uso de turbinas de pequeño tamaño y
continúan las investigaciones dirigidas a hacer más rentable adaptadas a velocidades de viento bajas, de modo que
la generación eólica y se ha iniciado el desarrollo de permitan ser usadas en un mayor número de lugares. Este
turbinas eólicas capaces de generar energía en zonas con desarrollo tendría gran aplicación en las zonas costeras del
vientos promedios bajos (entre 5.1 y 5.6 m/s) de forma país y los estados andinos.
rentable, aunque la mayor parte de los desarrollos se Esta investigación esta orientada al diseño de
realizan para zonas con vientos clase 4 (entre 5.6 y 6.0 m/s). turbinas eólicas (como la mostrada en la figura 1) que
En Venezuela sin embargo, el aprovechamiento puedan ser empleadas en desarrollos basados en la segunda
eólico es prácticamente nulo, situación ante la cual es muy perspectiva descrita, por lo cual el diseño debe satisfacer los
importante que se inicien las investigaciones y desarrollos siguientes requerimientos:
 2

- Tener la capacidad de generar en forma eficiente en zonas turbina

con vientos clase 4. N Velocidad de rotación en RPM
- Se debe evitar que el número de Reynolds en el álabe sea Pot Potencia en el eje de la turbina
menor a 500,000, dado que para tales condiciones pequeñas r Posición radial local
perturbaciones como cambios en la turbulencia del viento R Radio de la turbina
incidente, vibraciones en el álabe o rugosidad en la Re Número de Reynolds
superficie del mismo pueden producir grandes variaciones ro Radio de la raíz del álabe
en el comportamiento aerodinámico del perfil [2]. U Velocidad tangencial del segmento de álabe
- El aerogenerador debe ser de bajo y fácil mantenimiento. V Velocidad del viento en el borde de ataque del
- Poder funcionar sin el requerimiento de supervisión álabe
especializada. Vdiseño Velocidad de viento de diseño
- Poder ser construido con la tecnología presente en el país. Vtang Velocidad tangencial inducida en el viento
- Ser económico para poder competir con sistemas de Vwind Velocidad del viento libre
generación diesel. W Velocidad relativa
z Número de álabes
α Ángulo de incidencia
εmax Error máximo admisible
θ Ángulo de deflexión de los álabes
λ Radio de velocidades local (local speed ratio)
λo Radio de velocidades en la punta (tip speed
ratio)
ρ Densidad del aire
υ Viscosidad cinemática del aire
φ Angulo de la velocidad relativa respecto al
plano de rotación de la turbina
Figura 1. Diseño del Aerogenerador de tres palas. ω Velocidad de rotación

Nomenclatura Formulación del Problema

a Coeficiente de inducción axial
A Área barrida por los álabes El problema en estudio puede dividirse en dos
ap Coeficiente de inducción tangencial partes:1) Estudio de las variables que intervienen en el
c Cuerda comportamiento de una turbina eólica y determinación de
Cd Coeficiente de arrastre una geometría del álabe adaptada a las condiciones de
Cl Coeficiente de sustentación operación (Potencia generada = 1000 W, Velocidad de
Cp Coeficiente de potencia viento promedio = 5.8 m/s). 2) Diseño mecánico completo
Ct Coeficiente de empuje de un aerogenerador de 1000 W y determinación de un
dA Diferencial fuerza de arrastre proceso de fabricación del los álabes que sea sencillo y
dF Diferencial fuerza resultante de sustentación y permita un acabado de alta calidad.
arrastre Para ello el estudio se realizará en las siguientes
dFa Diferencial fuerza resultante en la dirección del etapas: 1) Modelaje teórico de los álabes mediante un
eje de la turbina algoritmo de cálculo que permita el estudio en condiciones
dFu Diferencial fuerza resultante en la dirección del estacionarias y obtención de la geometría óptima de una
movimiento del álabe turbina eólica, a la vez que se estudia la idoneidad de
dM Diferencial de momento en la dirección de distintos perfiles según las condiciones de flujo presentes en
rotación el álabe. 2) Construcción de un modelo que permita
dMa Diferencial de momento en la dirección axial determinar las características reales de la turbina y conocer
dR Ancho del segmento de álabe las dificultades asociadas con la construcción de los álabes.
dS Diferencial fuerza de sustentación 3) Propuesta de una metodología experimental que permita
F Factor de pérdida en la punta medir la potencia en la turbina para: distintas velocidades de
Fa Fuerza total en la dirección del eje de la viento, ángulos de calado y radios de velocidad en la punta,
turbina y determinar la localización de la separación del flujo (en
Fu Fuerza total en la dirección de la velocidad caso de presentarse). 4) En base a los resultados obtenidos
tangencial en (1) y (2), realizar el diseño mecánico completo de la
M Momento total en la dirección de rotación de la turbina y determinación de un proceso de fabricación de los
turbina álabes.
Ma Momento total en la dirección axial de la
 3

Modelaje teórico por lo cual la potencia mecánica necesaria en el eje de la

turbina es de 1372 W.
El modelaje teórico fue basado en la teoría de momento y la
teoría de elemento de álabe, con la cual se discretiza el Variables de entrada:
álabe en segmentos dR y para cada uno de ellos, aplicando ρ, ν, Potnecesaria, R
conservación de momento para volúmenes de control
Vdiseño, λo, Cl, Cd, z, dR, εmax, ro
anulares que abarquen el segmento de álabe en estudio, es
posible determinar las velocidades del flujo de aire a través
del mismo. Con las condiciones de flujo conocidas es r = ro, Fu = 0, Fa = 0
posible determinar las fuerzas actuantes sobre el álabe M = 0, Ma = 0
utilizando las ecuaciones de fuerza de sustentación y
arrastre para perfiles aerodinámicos en cada uno de los
segmentos. Para la determinación del flujo del álabe fue
tomada en cuanta la relación empírica desarrollada por Cálculo de U y λ
Glauert, incluyendo las perdidas en la punta del álabe. El
modelaje teórico aplicado tiene las siguientes limitaciones a = 0.2 (valor semilla)
[3]: 1) Se considera que la velocidad de viento libre es igual
para toda el área barrida por la turbina. 2) Se asume que las
fuerzas actuantes en el álabe son bidimensionales, es decir
se desprecia el flujo en la dirección radial. 3) Se desprecia el Cálculo de ap, V, φ, c, F, Ct y a*
efecto de un álabe sobre el otro. 4) Para el cálculo del a = a*
coeficiente de inducción axial se desprecia el efecto del
arrastre.
|a - a*| ≤ εmax NO
En la Figura 2 se muestran las velocidades del flujo en el
álabe y las fuerzas actuantes según el modelaje realizado. SI

Cálculo de dFu, dFa, dM

dMa y Re

Fu = Fu + dFu, Fa = Fa + dFa
M = M + dM, Ma = Ma + dMa

r = r + dR

r>R? NO

SI

Cálculo de Pot y Cp

FIN
Figura 2. Triángulo de velocidades y fuerzas
actuántes sobre un segmento de álabe.
Figura 3. Diagrama de flujo
El algoritmo calcula el comportamiento de la turbina para del algoritmo de cálculo.
un radio dado, para obtener el radio que permita producir la
potencia deseada es necesario tantear hasta conseguirlo.
Los valores de dR y εmax fueron variados hasta que su
Hay que tener en cuenta que la potencia en la turbina debe
valor no intervenía en los resultados obtenidos por el
ser mayor a la potencia eléctrica requerida debido a las
algoritmo. El valor de ro fue fijado de modo que permitiese
pérdidas en el generador y sistemas de control. El sistema
tener un espacio suficiente para acoplar los álabes a el eje
de control eléctrico estaría conformado por un convertidor y
de giro sin comprometer la resistencia mecánica del acople
un inversor. Se consideró que la eficiencia del generador,
y que además permitiese reducir de forma suave la cuerda
convertidor e inversor es de 90% para cada uno de modo
del álabe desde c(ro) hasta el diámetro de la barra sujetadora
que la eficiencia global del sistema eléctrico sería de 72.9%,
 4

del álabe. Sin embargo, dado que es necesario tener el tomaron estos rangos de valores para cada una de estas
diseño aerodinámico completo para conocer las fuerzas variables porque abarcan los que en forma general se
aerodinámicas actuantes sobre el álabe que son las que van presenta en la práctica. Para esto se puede emplear la
a determinar las dimensiones mínimas del acople, se tomó siguiente ecuación:
un valor base de ro para realizar los estudios aerodinámicos
 0.67
λo ⋅ z 1.92⋅ λo ⋅ z Cd
2 (2)
y luego este fue modificado al realizar el diseño mecánico. Cp , max 0.593⋅  − ⋅
Se tomó como base ro = 0.075 m.  (
 1.48 + z
0.67 )
− 0.04 ⋅ λo + 0.0025λo
⋅
2 1 + 2⋅ λo⋅ z Cl 

Determinación de los parámetros de diseño: Esta es una ecuación empírica que permite estimar el
La determinación de las demás variables de entrada hasta máximo Cp que se puede obtener en función de λo, z y
ahora no conocidas como son λo, z, Cl y Cd se hará Cl/Cd. En la Figura 4 se puede ver la gráfica de Cp,max en
tomando en cuenta los requerimientos planteados función de λo para z = 3 Cl/Cd = 120, 95 y 60; en la Figura
anteriormente para nuestro aerogenerador. Las 5 para z = 2. En primer lugar se puede observar que el
implicaciones de estos requerimientos desde un punto de coeficiente de potencia aumenta a medida que aumenta
vista técnico serían: Cl/Cd, es mayor para tres álabes que para dos y para la
configuración de tres álabes alcanza su máximo entre λo
- Alta eficiencia para vientos clase 4: se debe procurar un entre 6.5 y 10 mientras para la configuración de dos álabes
alto coeficiente de potencia, se debe evitar que se produzca alcanza su máximo entre 7.5 y 12.
separación laminar en el álabe y se debe procurar que el
perfil seleccionado tenga un rango amplio de ángulo de Ya que se sabe como influye el coeficiente Cl/Cd en el
incidencia antes de entrar en pérdida. coeficiente de potencia, es necesario estudiar su influencia
en el número de Reynolds, lo que se puede hacer utilizando
- Buen funcionamiento con números de Reynolds bajos: se el algoritmo de cálculo (Ver Tablas 1 y 2).
debe seleccionar un perfil aerodinámico que tenga buenas
características a bajos Re y que minimice los problemas de El número de Reynolds se calculó con la siguiente
separación laminar. ecuación:
- Poco y fácil mantenimiento: se debe evitar que los álabes W⋅c (3)
Re
sean muy esbeltos o que posean secciones muy delgadas, así ν
mismo se debe evitar que la cuerda no sea muy pequeña en
las puntas del álabe de modo que los álabes sean Los valores mostrados en las tablas 1 y 2 fueron calculados
mecánicamente resistentes, y puedan resistir ráfagas de usando ρ = 1.204 Kg/m3 υ = 1.15 x 10^-5 m2/s, Pot necesaria =
vientos fuertes y que resistan el fenómeno de fatiga debido a 1372 W, V diseño = 5.8 m/s, Cl = 0.8, λo = 7.
las cargas cíclicas presentes en la turbina. Además es
necesario que la turbina posea la mayor estabilidad
mecánica posible. En las tablas 1 y 2 se puede observar que el número de
Reynolds aumenta a medida que disminuye el cociente
- Construcción fácil: se deben considerar los mismos
Cl/Cd debido al aumento de las dimensiones de los álabes
aspectos que en el punto anterior, además si es posible se
por la disminución del coeficiente de potencia,
debe optar por perfiles que tengan su lado de alta presión
presentándose en promedio una diferencia de 6.2% entre el
recta. También hay que tomar en cuenta el tamaño de la
caso de Cl/Cd = 120 y el de Cl/Cd = 40. También se puede
turbina no sería económicamente adecuado utilizar algún
observar que existe una diferencia considerable en los
mecanismo de variación del ángulo de inclinación de los
números de Reynolds presentes en turbinas de tres y dos
álabes.
álabes, dado que las dimensiones de una turbina de dos
Cabe destacar que de las variables por determinar la que álabes tienen que ser mayores para poder extraer la misma
más caracteriza el comportamiento de la turbina es el radio cantidad de energía del viento que una de tres; para λo = 7
de velocidad en la punta (λo), que viene definido por: en la turbina de dos álabes se presenta en promedio un
aumento en el número de Reynolds de 50% respecto a la
ω⋅ R (1) turbina de tres álabes.
λo
V
En las figuras 6 y 7 se presenta con detalle la distribución de
Reynolds y cuerda a lo largo del álabe para dos y tres
álabes, λo = 7, Cl/Cd = 70 y Cl = 0.8.
Tomando en cuenta los requerimientos planteados, es
necesario determinar que valores de las variables de entrada
desconocidas permiten cumplir con los requerimientos
planteados. Esto se hará calculando Cp, Re y c para: λo є
(3,12), z є (2, 3), Cl є (0.4, 1.2) y Cl/Cd (60, 120). Se
 5

Figura 6. Número de Reynolds vs radio(m)

Figura 4. Coeficiente de potencia en función para 3 y 2 álabes, λo = 9, Cl/Cd = 80 y Cl = 0.8.
de λo para turbina de 3 álabes.

Figura 5. Coeficiente de potencia en función

de λo para turbina de 2 álabes.
Figura 7. Cuerda vs radio(m) para 3 y 2 álabes,
λo = 9, Cl/Cd = 80 y Cl = 0.8.
Tabla 1. Número de Reynolds para turbina de tres álabes.
De las figuras 4.3 y 4.4 se puede ver que la geometría y el
Reynolds / 105 comportamiento del flujo en turbinas de dos y tres álabes
Cl/Cd | r/R 0.2 0.4 0.6 0.8 1 tienen diferencias significativas. La configuración de dos
40 3.36 3.65 3.72 3.75 3.75
álabes tiene la ventaja que los álabes tienen mayores
dimensiones por lo cual, como ya se había mencionado, el
70 3.24 3.52 3.57 3.61 3.62 número de Reynolds es 50% mayor a la configuración de
90 3.21 3.48 3.54 3.57 3.58 tres álabes, y por lo visto en las tablas 4.1 y 4.2, en el
proceso de diseño va a ser de importancia procurar
120 3.17 3.44 3.50 3.53 3.53
aumentar el número de Reynolds en el álabe, ya que en
varios casos este cayó por debajo de 500.000. Por otra parte
Tabla 2. Número de Reynolds para turbina de tres álabes. las mayores dimensiones del álabe proporcionan la ventaja
adicional de que le brinda mayor resistencia mecánica a los
Reynolds / 105 álabes. Otra ventaja de la configuración de dos álabes es
Cl/Cd | r/R 0.2 0.4 0.6 0.8 1 que se facilita considerablemente la construcción de la
turbina. En lo que respecta a la configuración de tres álabes,
40 5.05 5.49 5.59 5.63 5.65 esta tiene la ventaja adicional de que brinda una mayor
70 4.86 5.28 5.37 5.41 5.42 estabilidad mecánica, lo cual facilita evitar los problemas de
90 4.81 5.22 5.31 5.35 5.36 vibraciones. Del análisis hecho hasta ahora ya es posible
fijar el primer parámetro de diseño: para las condiciones de
120 4.76 5.17 5.26 5.29 5.31 funcionamiento en que operará la turbina la configuración
que brinda más ventajas es la de dos álabes, aunque hay que
mencionar que a la hora de la construcción del prototipo (en
 6

un proyecto futuro) será necesario un balanceo minucioso disminuyen demasiado. Aun falta estudiar el valor o rango
para evitar problemas de vibraciones. óptimo de Cl. En la Tabla 3 se muestra los valores de Cp
para distintos Cl. En las Figuras 11 y 12 se muestra la
Ya con el número de álabes seleccionado, se estudiará la
influencia de Cl en el Re y c.
velocidad específica. En las figuras 8 y 9 se muestran la
influencia de λo en el número de Reynolds y la cuerda en el
álabe, respectivamente para z = 2, Cl/Cd = 70, Cl = 0.8. En
estas se puede ver que a medida que disminuye λo aumenta 250
el número de Reynolds y la cuerda, de manera que valores
bajos de λo ayudan a evitar los problemas de separación 200
laminar y mejoran la resistencia mecánica y facilidad de
construcción. También valores bajos de λo disminuyen los 150

RPM
problemas de vibraciones y los problemas de incrustaciones
de partículas sólidas en el álabe, al funcionar este a menores 100
velocidades tangenciales. Pero por otro lado valores altos de
λo facilitan conseguir un generador que funcione a la 50
velocidad de giro de la turbina, simplifica la caja de
velocidades en caso de ser empleada y también disminuyen 0
los momentos torsores en el eje de la turbina. 4 6 8 10 12
λo

Figura 10. N (RPM) vs λo.

Tabla 3. Coeficiente de potencia vs. Coeficiente de

sustentación.
Cl 0.4 0.8 1.2
Cp 0.524 0.524 0.524

Figura 8. Número de Reynolds vs radio(m)

para λo = 6, 9 y 12, tres álabes, Cl/Cd = 80 Cl = 0.8.

Figura 11. Número de Reynolds vs radio (m) para

Cl = 0.4, 0.8 y 1.2, λo = 8, z = 2 y Cl/Cd = 70.

Figura 9. Número de Reynolds vs radio(m)

para λo = 6, 9 y 12, tres álabes, Cl/Cd = 80 y Cl = 0.8.

En la Figura 10 se muestra el gráfico de N vs λo. A partir de

esto consideramos adecuado tomar λo = 8, de modo que N ≈
150 RPM y el número de Reynolds y la cuerda no
 7

Cos estos criterios se estudiaron distintos perfiles usados

normalmente para aplicaciones a bajas velocidades. Dos de
los que cumplían mejor con nuestros requerimientos son los
FX 63-137 y E387. Las curvas Cl vs α de estos perfiles son
mostradas la Figuras 12 y 13, y las curvas Cl/Cd vs α en las
Figuras 14 y 15.

FX 63-137 (C)
2
1,8
1,6
1,4 Re = 100,019
1,2
1
Re = 200.082
0,8

Cl
0,6
Figura 12. cuerda vs radio (m) para 0,4
Re = 349,958
Cl = 0.4, 0.8 y 1.2, λo = 8, z = 2 y Cl/Cd = 70. 0,2
0 Re = 499,818
-0,2
De la tabla 3 se puede ver que el valor de Cl no influye en el -10 -5 -0,4
0 5 10 15 20
coeficiente de potencia y de los gráficos en las figuras 10 y -0,6
11 se puede ver que al disminuir Cl aumenta la cuerda y el alfa
número de Reynolds, por lo cual se puede concluir que
conviene un perfil aerodinámico con un coeficiente de Figura 13. Curva Cl vs α para el perfil FX 63-137
sustentación bajo.
De esta manera para el diseño de los álabes se requiere:
1,4
- Alto coeficiente Cl/Cd.
1,2
E387 (E )
- Número de álabes a determinar según los resultados
1
obtenidos con las pruebas con modelos.
0,8
- λo = 8 Re = 100,029
0,6 Re = 199,879
- Bajo coeficiente Cl.
Cl

Re = 300,289
0,4
Re = 349,865
0,2
Selección del perfil Re = 459,961
0 Re = 500,094
Del análisis anterior se sabe que necesitamos un perfil que -10 -5 0 5 10 15
tenga un alto valor de Cl/Cd y un bajo valor de Cl, además -0,2
no debe ser propenso al fenómeno de separación laminar.
-0,4
Por otro lado, tomando en cuenta que el ángulo de alfa
inclinación del álabe no podrá ser variado durante el
Figura 14. Curva Cl vs α para el perfil E387
funcionamiento, se necesita que el perfil soporte el mayor
rango de variaciones posible en el ángulo de calado sin que
el álabe entre en pérdida.
En primer lugar se buscarán perfiles que tengan un buen
comportamiento para las condiciones de diseño, para
Reynolds entre 1 x 10^5 y 5 x 10^5. Para ello el orden de
prioridades propuesto será:
1) Amplio rango de α para el cual el perfil no entre en
pérdida y Cl/Cd se mantenga alto.
2) Pocos problemas de separación laminar.
3) Bajo coeficiente de sustentación.
4) Alto coeficiente Cl/Cd.
 8

FX 63-137 (C) Desventajas: (1) El rango en el cual Cl/Cd se mantiene alto

120,00 no es muy grande, para Re = 350,000 este se mantiene por
encima de 60 sólo entre 1° ≤ α ≤ 8.5°.
100,00
Re = 100,019 Ventajas: (1) Se puede ver que no hay mucha diferencia
80,00
Re = 200,082 entre las curvas de Cl vs α para distintos Re, lo cual nos
60,00 Re = 349,958 muestra que para los números de Re presentados en la curva
Cl/Cd

Re = 499,818 no se produce separación laminar. (2) No posee coeficientes

40,00 de sustentación muy altos, se mantiene menor a 1 para α ≤
5° y un valor máximo de 1.25 para α = 10°.
20,00
En base a lo anterior se considera que el perfil más
0,00
adecuado es el FX 63-137, dado que es el que mejor
-10 -5 0 5 10 15 20
-20,00 satisface el requerimiento prioritario. Para este perfil el
cociente Cl/Cd se mantiene por encima de 80 para un rango
alfa
de 13° mientras que para el perfil E387 este rango es de 8°,
lo cual representa una diferencia de 40%.
Figura 15. Curva Cl/Cd vs α para el perfil FX 63-137
Estimación del ángulo de calado óptimo
E387 (E ) Para la determinación del ángulo de calado óptimo de los
120,00 álabes usando el perfil FX 63-137 es necesario tomar en
cuenta: (1) para radios de velocidades en la punta distintas a
100,00 Re = 100,029 la nominal el ángulo de incidencia se incrementa
Re = 199,879 considerablemente (Ver Figura 17), (2) el coeficiente de
80,00
Re = 300,289 sustentación del perfil seleccionado es relativamente alto,
60,00 Re = 349,865 (3) el perfil puede presentar problemas de separación
Cl/Cd

Re = 459,961 laminar, especialmente si los números de Reynolds son

40,00 Re = 500,094 bajos por debajo de 200,000.
20,00

0,00
-10 -5 0 5 10 15
-20,00
alfa
Figura 16. Curva Cl/Cd vs α para el perfil E387

Estudiando las gráficas se puede ver:

Para el perfil FX 63-137:
Desventajas: (1) La diferencia entre la curva Cl vs α para Re
= 100,019 (Ver Figura 12) con las curvas para los demás
Reynolds nos indica que se produce separación laminar; Figura 17. Diferencia del ángulo de incidencia para
esta separación puede producirse a Reynolds mayores distintos λ con la el presente en condición nominal
debido a irregularidades o vibraciones en el álabe, o
también si el viento al llegar al alabe es muy turbulento. (2) En base a lo descrito anteriormente y con el fin de evitar
Se puede observar que (Ver Figura 12) que el coeficiente de separación laminar se fijará el número de Reynolds para la
sustentación es bastante alto, para las curvas de Re ≥ condición nominal en 500,000 para r = ro, que es la sección
200,000 es mayor a 1 para α ≥ 2, con un valor máximo de en donde el número de Reynolds es menor. Como se dijo
1.7 para α = 11°. anteriormente ro va a depender del diseño mecánico del
Ventajas: (1) El cociente Cl/Cd se mantiene alto para un prototipo el cual no es descrito en este paper. El valor de ro
amplio rango de α, para Re = 350,000 este se mantiene por final es de 0.35 m. Conociendo el valor de ro, z, λo y ya
encima de 60 para -2° ≤ α ≤ 11°. seleccionado el perfil ya se puede buscar un ángulo de
calado que permita cumplir con la condición del número de
Para el perfil E387: Reynolds y que tenga un cociente de Cl/Cd lo más alta
 9

posible. Haciendo esto obtenemos que los requerimientos geometría teórica y luego fue modificada levemente para
planteados se satisfacen para α = -3.06, para el cual Cl = garantizar que el área obtenida por el diseño teórico y la real
0.486 y Cl/Cd = 44.59. de los modelos fuese igual.
Estos álabes se probarán en las configuraciones de dos y
Pruebas con modelos tres álabes, aunque tomando en cuenta que la configuración
de dos álabes va a estar funcionando en condiciones fuera
Las pruebas con modelos serán llevadas a cabo con el
de diseño dado que este se realizó para tres álabes.
objetivo de determinar:
La construcción de los álabes se realizó siguiendo el
1) Estudio del comportamiento real de la turbina. Se
siguiente procedimiento:
realizarán las curvas de potencia generada vs velocidad del
viento para distintos radios de velocidad en la punta y - Se construyó la forma del álabe en anime, lo cual se hizo
distintos ángulos de inclinación de los álabes. en tres secciones de forma que fuese más sencilla su
construcción. Estas secciones abarcarían los siguientes
2) Dificultades asociadas a la construcción de los álabes.
radios: de ro = 0.125 m a r = 0.485 m, de r = 0.485 m a r =
Para esto es necesario: 0.725 m y de r = 0.725 m a r = 0.825 m, Para la
construcción se hizo la forma del álabe en negativo en
Para prueba 1: a) Construcción y montaje de modelos. forma para cada uno de los extremos de las secciones
Para la prueba 2: a) Igual que (1.a). descritas, luego estas fueron colocadas en piezas de anime
de longitud igual a la de cada una de las secciones, con
b) Generador eléctrico. diferencias de ángulos entre ellos iguales a las diferencias
c) Dispositivo electrónico de medición y registro de de ángulo de incidencia para las condiciones nominales.
velocidad del viento, dirección del viento, velocidad de Luego las piezas de anime fueron cortadas por medio de un
rotación de la turbina, potencia eléctrica generada y cable caliente que era apoyado y desplazado en cada uno de
orientación de la turbina. los perfiles de fórmica.
De estos aspectos analizaremos en detalle la construcción de - Las secciones de anime ya cortadas correspondientes a
los modelos. cada álabe fueron pegadas juntas.
1.a) Construcción de modelos: Se construyeron tres álabes - Las formas de los álabes en anime fueron forradas con
diseñados para las siguientes condiciones nominales: láminas de madera balsa de 1 mm de espesor.
Vdiseño = 8 m/s, z = 3, λo = 7, Cl = 0.8, Cl/Cd = 100. Las - Los álabes fueron reforzados en la raíz, la unión de las
características del modelo fueron: ro = 0.125 m, R = 0.825 secciones, el borde de ataque y el borde de fuga con fibra de
m, la distribución de cuerda y ángulo de inclinación para las vidrio tejida.
condiciones nominales son mostrados en la Figura 21.
- Se agregó una capa de de mastique la cual fue lijada
30 cuidadosamente con lo cual se logró obtener una superficie
c (cm) muy lisa y regular.
25
fi (grados)
20 - Por último se forraron los álabes con MONOKOTE, con el
objetivo de garantizar que la superficie fuese
15
completamente liza y que soportase la humedad y la lluvia.
10
Los álabes van sujetos a un buje de madera que permite
5 fijarlo pero que también permite variar su ángulo de
0 inclinación. El buje va fijado a un eje metálico el cual se
-5
0 0.2 0.4 0.6 0.8 1 conecta por medio de un mecanismo de cadena y poleas al
generador. Este mecanismo a su vez va sujeto a la torre a
r (m)
través de un eje que permite que la turbina se mantenga
Figura 19. Distribución de cuerda y ángulos en el álabe. alineada en la dirección del viento, para lo cual se va a
contar con una veleta colocada en la parte trasera del
aerogenerador.
Cómo se puede ver las distribuciones de cuerda y ángulo de
inclinación fueron linealizadas. Esto se hizo debido a que el
método de construcción seleccionado (que es explicado más RESULTADOS
adelante), no permite la fabricación de los álabes con la Diseño teórico: Para las condiciones finales expuestas en la
geometría teórica. Esta linealización de la geometría fue sección de modelaje teórico se obtuvo que el diseño final
realizada procurando la menor desviación posible a la fue:
 10

Tabla 4. Resultados del diseño teórico del álabe

λo z ro (m) R (m) α Cl Cl/Cd
8 2 0.35 2.913 -3.06° 0.486 44.59 CONCLUSIONES
Para el diseño de turbinas eólicas de potencias menores o
iguales a 1 KW, es muy importante el estudio del número de
Reynolds en los álabes, dado que las dimensiones de estos
pueden ser lo suficientemente reducidas como para que, con
un diseño no apropiado, el número de Reynolds sea menor a
500,000 (cabe recordar que para Re ≤ 500,000 pequeñas
perturbaciones como cambios en la turbulencia del viento
incidente, vibraciones en el álabe o rugosidad en la
superficie del mismo pueden producir grandes variaciones
en el comportamiento en el comportamiento de un perfil, lo
cual no es deseado en la turbina). Ante esto, se debe
procurar incrementar la cuerda de los álabes, lo cual
favorece el uso de turbinas de dos álabes en vez de tres, con
un radio de velocidad en la punta de de diseño medias
(alrededor de 8).
Adicionalmente las dimensiones reducidas del álabe
dificultan la construcción de los mismos y su resistencia
Figura 20. Distribución de ángulo de incidencia. mecánica no es muy buena. Con la construcción de los
modelos se pudo observar que para evitar estos problemas
se debe procurar que la cuerda no sea menor a 5 cm.
En lo que respecta al proceso de fabricación, la construcción
de los modelos mostró que la construcción de los álabes
debe ser realizada con moldes, dado que si se realizan por
separado de forma manual, se presentan algunas diferencias
entre los álabes lo cual produciría vibraciones durante el
funcionamiento de la turbina.

REFERENCIAS
[1] “Record Year For Wind Energy: Global Wind Power
Market Increased By 43% In 2005”. Awea News. Global
Wind Energy Council. 17 / 02 / 2006 22:34
Figura 21. Distribución de cuerda.
[2] Spera, David A. 1994, “Wind Turbine Technology”.
ASME Press. EEUU.
Construcción de modelos: Con la construcción de los álabes
se pudo observar: [3] Manwell J. F., McGowan J. G., Rogers A. L. 2002.
- El método de fabricación utilizado no es el más adecuado, “Wind Energy Explained”. John Wiley & Sons.
dado que implica la fabricación desde cero y por separado
de cada uno de los álabes por lo cual: (1) se presentan [4] Moriarty P. J., Hansen A.C. 2005. “AeroDyn Theory
algunas diferencias entre los álabes lo cual producirá que Manual” Technicall Report. National Renewable
tengan comportamientos aerodinámicos distintos entre ellos, Laboratory. EEUU.
que a su vez producirá vibraciones en la turbina (2) y la
fabricación de los álabes requiere mucho tiempo.
- Para dimensiones de cuerda menores a 5 cm se hace muy
difícil lograr la forma del perfil aerodinámico deseada y
disminuye mucho la resistencia mecánica del álabe.

View publication stats