Centro Azúcar 38(4): 42-46, octubre-dic.

, 2011

Estudio aerodinámico de los coeficientes

de sustentación y arrastre para pequeños
aerogeneradores.

Aerodynamic study of sustentation and

drag coefficients for small wind turbine

Ernesto Yoel Fariñas Wong* , Sergio Jauregui Rigo*.

* Universidad Central de Las Villas, Centro de Estudios Energéticos
y Tecnologías Ambientales CEETA. Carretera a Camajuaní km 5,5,
Santa Clara, Villa Clara, Cuba, CP: 54830. Telf: 5342281194.
Palabras claves: Pequeños aerogeneradores, coeficientes
aerodinámicos no lineales,
Key words: small wind turbine, nonlinear aerodynamic coefficients
*[email protected]

Resumen:
En la presente investigación se realizó una valoración sobre los métodos
para el análisis de los sistemas de orientación de en pequeñas
aeroturbinas de eje horizontal, con la finalidad de mejorar la captación de
la energía eólica actualizando los procedimiento metodológico para el
desarrollo de los sistemas de orientación de estas máquinas, mediante la
utilización de coeficientes de sustentación y arrastre de comportamiento
variable en función de la razón de aspecto.
Abstract:
In the present study was carried out an assessment of methods for the
analysis of guidance systems in small horizontal axis wind turbines, in
order to improve the capture of wind energy by updating the methodologi-
cal approach for the development of guidance systems of these ma-
chines, using lift and drag coefficients of variable behavior depending on
the aspect ratio.

42
 Centro Azúcar 38(4): 42-46, octubre-dic., 2011

Introducción deben mejorar el posicionamiento del rotor y la cola

Según Sagrillo (1997) el desarrollo de los pequeños e influir positivamente en la eficiencia de la máquina
aerogeneradores se ha visto limitado por la poca por la reducción de los problemas aeroelásticos.
profundidad en el análisis de los fenómenos Material y Método
aerodinámicos, por lo que puntualiza la necesidad La propuesta de Piggott (2001) como sistema de
de un análisis específico sobre los sistemas de orientación y protección en pequeños
orientación y protección, los que han sido en gen- aerogeneradores tropezó con la inexistencia de un
eral deficientes provocando problemas aeroelásticos aparato analítico para su desarrollo, condiciones que
en el rotor, la torre y el generador eléctrico, así como determinaron que el sistema fuera dimensionado
ruido aerodinámico, reflejo de las pérdidas en la mediante pruebas de campo por tanteo y error
captación de energía del viento. (Fariñas 2004) utilizando las dimensiones de otros
Gipe (2000) planteó que para revertir la situación sistemas, haciendo tedioso, impreciso y costoso
era necesario un estudio tanto analítico como ex- llegar a configurar el sistema. En base a esta situación
perimental de los sistemas de orientación y fue desarrollado el procedimiento de Mas y Batet
protección utilizados en estas máquinas. En este (2006), dicho procedimiento se basa en determinar
marco Piggott (2001) propone la utilización del los ángulos formados entre el eje del rotor (y) y el
“sistema de orientación y protección por momento de la cola (f) para diferentes velocidades del viento
de cola” (SOPMC) en aerogeneradores de eje hori- a partir de las características geométricas de la
zontal con rotor frente al viento, versiones de este aeroturbina, con el valor de estos ángulos se procede
sistemas se estudiaron en diferentes prototipos por a determinar la curva de potencia de la máquina.
Corbus (2006) y Bikdash (2004); estos reportes
concluyen que la utilización del sistema propuesto
por Piggott (2001) era viable, sin embargo, para
aumentar la eficiencia y atenuar los problemas
aeroelásticos presentes en los prototipos ensayados
era necesario el desarrollo de métodos que tuvieran
en consideración la aerodinámica para el cálculo de
la veleta de cola, debido a que esta se había
convertido en un componente activo en la regulación
y no solo como un elemento para la orientación o
decorativo.
En la Mas y Batet (2006) publican un procedimiento
para el cálculo y análisis de aerogeneradores con
SOPMC, este tiene en consideración las
recomendaciones antes expuestas, sin embargo, los
coeficientes aerodinámicos de sustentación y Figura 3 Esquema de fuerzas sobre la cola del
arrastre que influyen sobre la veleta del sistema de aerogenerador CEETA-01.
orientación y protección son asumidos constantes
para cualquier velocidad del viento y posición de la Hasta la presentación de está propuesta el
veleta, tratamiento que los autores consideran una procedimiento analítico aceptado para estos cálculos
limitación del procedimiento, debido a que el valor es el planteado por Mas y Batet (2006), e en opinión
de los coeficientes influyen en el posicionamiento del autor dicho procedimiento presenta limitaciones
del rotor y la cola respecto al viento, el cual si es relacionadas con el tratamiento aerodinámico de la
incorrecto provocará pérdidas en la captación de la veleta de cola, las cuales son corregibles mediante
energía del viento así como la ocurrencia de la utilización de métodos numéricos y coeficientes
problemas aeroelásticos. aerodinámicos variables. El procedimiento de Mas
Los autores consideran que es posible modificar el y Batet (2006) utiliza la teoría de Betz para el cálculo
procedimiento planteado por Mas y Batet (2006), del rotor omitiendo de esa forma las características
mediante la utilización de coeficientes de sustentación del rotor y su influencia sobre la corriente que toca
y arrastre de comportamiento variable los cuales la cola. En base a estas limitaciones se realiza un
43
 Centro Azúcar 38(4): 42-46, octubre-dic., 2011
estudio del comportamiento de los coeficientes Bartlett y Vidal (1955) refieren la ecuación 4 para
aerodinámicos variables para su inclusión en el lograr la curva de CL vs f mediante una aproximación
procedimiento. polinómica de segundo orden, la cual puede ser
Coeficiente aerodinámico de sustentación lineal empleada en geometrías tan disímiles como alas del-
Las veletas de cola de los sistemas de orientación y tas, placas rectangulares y alas rectangulares
protección de pequeñas aeroturbinas generalmente delgadas, presenta como restricción valores de AR
poseen figuras geométricas conocidas con < 3 y números de Reynolds por debajo de 104.
características aerodinámicas estudiadas en túneles
de viento. Debido a la similitud y los análisis de C L  C L  C DC 2 (4)
placas planas realizados en experimentos de
aerodinámica, en este trabajo se considera la veleta
de cola del SOPMC como una placa plana inmersa La restricción viene dada por la determinación de la
en un fluido, extendiendo así la teoría a la aplicación constante independiente CLf pendiente de la curva
práctica en el cálculo de aerogeneradores. de sustentación y se determina mediante la ecuación
A principios del siglo pasado Ludwing Prandtl 5, mientras que el coeficiente de flujo cruzado CDC
planteó (ecuación 1) como caracterizar el se determina por la ecuación 6, este último no
coeficiente de sustentación (CL) de una placa plana presenta restricciones en cuanto a la razón de
con AR<3 y flujo laminar. aspecto y número de Reynolds.
C L  1,97 AR 0, 72 (5)

CDC  2,82 AR 0,98 (6)

(1)
Polhamus (1966) utiliza para el cálculo del coeficiente
de sustentación en placas rectangulares, circulares y
d e ecuación 7, donde Kp es el coeficiente
l t a s ,

potencial de sustentación, el cual es función de AR y

L: longitud característica la forma del borde de ataque, Kv que es el coeficiente
Sc: área expuesta al flujo de vórtice de sustentación.
Sobre este mismo tema Katz y Plotkin (2001)
abordan el cálculo del coeficiente de sustentación CL  K p sen cos    v cos sen 2 (7)
estático, ecuación 3, para valores de AR<<1 y
ángulos de ataquesf<5º. Cuando se desprecia la acción del arrastre por
 fricción el coeficiente de arrastre (CD0) producto de
C L  AR (3) esta consideración puede ser calculado en la
2
El autor considera que no es viable la utilización de ecuación 8:
estas teorías en el análisis del SOPMC, debido a C D0  C L tan (8)
que el ángulof es por lo general mayor a cinco Para los cuales Polhamus (1966) definió los términos
grados, el flujo está perturbado por la presencia del Kp y Kv en función de la razón de aspecto, en las
rotor y la torre, mientras que la razón de aspecto ecuaciones 9 y 10.
para el aerogenerador CEETA-01 varía de 3 a 9.
Kp  0,0051 AR 2  0,2406 AR  1,31 (9)
Coeficiente de sustentación no lineal
En trabajos publicados en el siglo anterior por
Bartlett y Vidal (1955), Polharmus (1966) y Lamar Kv  3,41e 0 , 034 AR (10)
(1976) plantean distintos modelos matemáticos para
describir la no linealidad del coeficiente de Lamar (1976) extiende los trabajos de Polhamus
sustentación a través de la adición de términos (1966) a otras geometrías y a bordes laterales, los
relacionados con el ángulo de ataque.
44
 Centro Azúcar 38(4): 42-46, octubre-dic., 2011

vórtices de sustentación en los bordes laterales se CD0 es determinado por la ecuación 8 y el parámetro
estiman de la misma manera que los de succión en Di factor de arrastre inducido, ecuación 14, el cual
el borde de ataque, expresándose el coeficiente de depende del factor de eficiencia t .
sustentación mediante la ecuación 11. 1
Di  (14)
C L  K P sen cos 2   ( K v , le  K v , se ) cos sen 2  ( AR)
(11) Para cuerpos simétricos, el segundo término de la
ecuación 13 es el mismo que el de la ecuación 12,
El término Kv,le es igual a Kv de la ecuación 7, que pero la constante Di incluye la contribución del
representa la sustentación generada por el borde arrastre parásito y del arrastre inducido. La
de ataque y K v,se es la parte del vórtice de ecuación 13 es utilizada para aproximar el arrastre.
sustentación generado por el borde lateral. Mediante el análisis realizado es posible determinar
La utilización de la ecuación 11 tiene la dificultad los coeficientes de sustentación y de arrastre con un
que los coeficientes Kv,le y Kv,se aún se encuentran comportamiento no lineal para diferentes ángulos de
poco caracterizados y sus valores difieren ataques en base a una geometría predeterminada y
considerablemente de un reporte a otro, por lo que número de Reynolds calculado en función de la
los autores no recomienda aún la utilización de dicha velocidad y la magnitud característica de la veleta
ecuación, propone utilizar en el análisis de la cola de cola, siendo viable su introducción en el
de los pequeños aerogeneradores la ecuación 7 de procedimiento de cálculo desarrollada por Mas y
Polhamus (1966), la cual está mejor caracterizada Batet (2006) para el desarrollo de sistemas de
y ha sido empleada por varios autores (Wood, 1991; orientación en pequeños aerogeneradores, cálculo
Prado, 2007) en otras ramas de la ingeniería para el que debe redundar en el mejoramiento de las
cálculo de los coeficientes de sustentación con condiciones aeroelásticas de las pequeñas turbinas
buenos resultados. eólicas

Coeficiente de arrastre en función de la sustentación Conclusiones

La razón del coeficiente de arrastre y el de El análisis de los sistemas de orientación se ha
sustentación en la cola del aerogenerador es funda- basado fundamentalmente en ensayos de campo por
mental para que exista un buen control sobre la tanteo y error, mientras que los procedimientos
posición de la máquina. El arrastre está compuesto metodológicos desarrollados para su análisis
por el arrastre parásito y el arrastre inducido. El presentan limitaciones con la aerodinámica de la
arrastre parásito es generado por la fricción super- veleta de cola y con el análisis del rotor.
Es necesario reducir los problemas aerodinámicos
ficial, la cual es despreciable en superficies pulidas
y R e
6
<, mientras que el arrastre de forma o
1 0 de los sistemas de orientación expuestos en esta
arrastre de presión se convierte en una función del investigación por las simulaciones numéricas, debido
ángulo de ataque y del área expuesta al fluido. a que estos favorecen el incorrecto posicionamiento
El cálculo del arrastre inducido si es derivado de ladel rotor y la cola respecto a la dirección del viento
teoría del comportamiento lineal de la sustentación contribuyendo a la baja eficiencia de los pequeños
puede obtenerse mediante la ecuación 12. aerogeneradores.
El procedimiento metodológico planteado por Mas
CL
2 y Batet para el análisis de los pequeños
CD  (12) aerogeneradores con SOPMC presenta limitaciones
 ( AR)
debido al empleo de coeficientes aerodinámicos de
Siendo comúnmente expresada la relación entre los valor constante para el análisis de la cola, dichos
coeficientes de sustentación y el de arrastre por la coeficientes en opinión del autor deben ser variables
polar de arrastre en la forma de la ecuación 13 en función de la razón de aspecto de la veleta.

C D  C D 0  DiC L
2
(13) Bibliografía
1.Bartlett, G.E. y Vidal, R.J., (1955) “Experimental
investigation of influence of Edge Shape on The

45
 Centro Azúcar 38(4): 42-46, octubre-dic., 2011
Aerodynamic Characteristics of Low Aspect Ratio
W i n g s a t Journal of the Aeronauti-
L o w S p e e d s . ”

cal Sciences, 34, pp 34-41. ISSN: 0022-4650

2.Bikdash, M.; Chen, D. A. y Harb, M., (2004) “A
Hybrid Model of a Small Autofurling Wind Turbine”,
Journal of Vibration and Control, Vol. 7, pp.127-
148. ISSN: 1392-2114
3.Corbus, D.; Meadors, M., (2006) “Small Wind
Research Turbine”, October 2006, Final Report
Technical Report NREL/TP-500-38550, http://
www.osti.gov/bridge,
4.Fariñas, E., (2004) “Generador de imanes
permanentes para aplicaciones industriales”. Revista
Centroazucar, No 4, 2004, Universidad Central de
las Villas, Santa Clara, Cuba. 2004, pp 46-50.
ISBN: 02535757
5.Gipe, P., (2000) “Noise from Small Wind Tur-
bines: An Unaddressed Issue”, en Chelsea Green
6.Katz, J.; Plotkin, A., (2001.) Low-Speed Aero-
dynamics, Second Edition, Editorial Cambridge
7.University Press, Cambridge. Pp. 243 ISBN:
978-0-521-87144-0
8.Lamar, J.E., (1976) “Prediction of vortex flow
characteristics of wings at subsonic and supersomic
speeds”. Journal of Aircraft. 21, pp 56-79. ISSN:
0002-2667
9.Mas,V.E. y Batet, Ll.M., (2006.) Estudio
mecánico del aerogenerador T 100.
10.Piggott, H., (2001) Cómo construir un
generador eólico, traducción al castellano, Edito-
rial Scorading Wind, United Kingdom, pp. 73-108.
ISBN: 9780975197059
11.Polhamus, E.C., (1966.) “A concept of the vor-
tex lift of sharp-edge delta wings based on a lead-
ing-edge suction analogy”, NASA TN D-3767
12.Prado, R. A., (2007) Development of a code
for the viscous-inviscid interaction oriented to
turbomachinery, Tesis de Doctorado de la
Universidad de Buenos Aires, Área Ingeniería, Ar-
gentina.
13.Sagrillo, M., (1997) “Small Wind Generators and
Noise”, Windletter, February/March 1997. ISSN:
0036-8075
14.Wood D.H. (1991) A three-dimensional analy-
sis of stall-delay on a horizontal-axis wind turbine.
15.Journal of Wind Engineering and Industrial
Aerodynamics, 37, pp.1-14. ISBN: 89-97782-
58-9.

46