Análisis aerodinámico del perfil FX 60-126 en una turbina eólica y

comprobación de la potencia mediante modelos analíticos y numéricos

Autor: Luis Manuel Mendoza Gomero

Grupo Revisor: Ríos Roldán Y.1; San Martín Alfaro J.2; Sánchez Salazar S.3

1. Estudiante de ingeniería mecánica del X ciclo UNT. ([email protected])

2. Estudiante de ingeniería mecánica del X ciclo UNT. ([email protected])
3. Estudiante de ingeniería mecánica del X ciclo UNT. ([email protected])

Resumen

En esta investigación se describe en primera instancia el proceso de determinación de valores más reales de los
coeficientes aerodinámicos del perfil Wortmann FX 60-126 para mejorar la confiabilidad de los resultados al
obtener el coeficiente de potencia de dicha turbina. La utilidad de esta investigación radica en la profundización
de los coeficientes aerodinámicos del perfil seleccionado, esto con la finalidad de obtener resultados
suficientemente confiables para validar el coeficiente de potencia obtenido de la turbina eólica mediante el uso de
diferentes modelos analíticos y métodos numéricos. Para ello hacemos uso de modelos analíticos -numéricos
(XFoil) y herramientas computacionales (Ansys-CFX) para contrastar los coeficientes aerodinámicos obtenidos
con los experimentales proporcionados por el Instituto de Aerodinámica y Gasdinámica de la Universidad de
Stuttgart. Luego se genera y valida el código propio (Python) elaborado en base a la teoría de BEM (Blade Element
Momentum) para el análisis de la turbina eólica. Posteriormente, haciendo uso de las fuentes aerodinámicas del
perfil y de las condiciones de operación y geometría de la pala obtenemos la potencia y el empuje que genera el
Waira 12. Una vez obtenidos los valores de potencia y empuje proveniente del modelo analítico -numérico, se
procede a realizar el modelado del rotor de la turbina eólica y posterior a ello su simulación fluido dinámica
computacional (Ansys-CFX). Los valores resultantes de potencia y empuje obtenidos a través modelo de BEM en
base a las diferentes fuentes aerodinámicas (experimental, Xfoil y CFD), se encuentran próximos entre sí, sin
afectar drásticamente los resultados finales del modelo.Palabras Clave: análisis de fallas, ejes, desgaste, ensayos
mecánicos, fatiga.

Keywords: Perfil Wortmann, coeficiente de potencia, métodos analíticos-numéricos

Abstract.

In this investigation, describe in the first instance the process of determining more real values of the aerodynamic
coefficients of the Wortmann FX 60-126 profile to improve the reliability of the results by obtaining the power
coefficient of said turbine. The usefulness of this radical investigation in the deepening of the aerodynamic
coefficients of the selected profile, this with the modification of obtaining sufficiently reliable results to validate
the coefficient of power obtained from the wind turbine through the use of different analytical models and
numerical methods. For this we use analytical-numerical models (XFoil) and computational tools (Ansys-CFX) to
compare the aerodynamic coefficients obtained with the experiments provided by the Institute of Aerodynamics
and Gasodynamics of the University of Stuttgart. Then the own code (Python) developed based on the BEM theory
(Blade Element Momentum) for the wind turbine analysis is generated and validated. Subsequently, using the
aerodynamic sources of the profile and the operating conditions and geometry of the blade, we obtain the power
and thrust generated by the Waira 12. Once the power and thrust values from the analytical-numerical model are
used , we will proceed to perform the wind turbine rotor modeling and after its computational dynamic fluid
simulation (Ansys-CFX). The resulting power and thrust values obtained through the BEM model based on the
different aerodynamic sources (experimental, Xfoil and CFD), are close to each other, without drastically affecting
the final results of the model. Key Words: failure analysis , axles, wear, mechanical tests, fatigue.

Keywords: Wortmann profile, power coefficient, analytical-numerical methods

I. INTRODUCCIÓN

El sistema energético se encuentra en un punto

crítico ya que el consumo masivo de combustibles
fósiles, que son la principal fuente de energía, es El presente trabajo de investigación pretende
mayor responsable de las emisiones que están determinar un método aumento de confiabilidad de
comenzando a provocar cambios climáticos severos. resultados de la potencia realizados por método
Para lo cual las energías renovables son una analítico y numérico. Esto con el fin de precisar los
alternativa real para cubrir la demanda energética resultados y conocer el máximo rendimiento de la
que va en aumento y en reducir las emisiones de turbina a trabajar. En primer lugar se utilizará los
gases contaminantes que agravan el problema del softwares ANSYS-CFX y XFOIL para encontrar los
cambio climático. Estas energías, frente a las fuentes coeficientes aerodinámicos y corroborarlos con los
tradicionales que producen emisiones de del Instituto de Stuttgart. Obteniendo estos valores
contaminantes, como el carbono gas, tienen se procederá a analizar la turbina Waira mediante el
características que actualmente las hacen viables y método numérico BEM y el software ANSYS-CFX.
competitivas debido al decremento de sus costos de A causa de las bondades que representa el uso de las
producción durante los últimos años. turbinas eólicas como fuente de producción de
Actualmente la sociedad se encuentra más energía, en el estudio realizado busca profundizar en
consciente y preocupada por los impactos el estudio del rotor de la turbina eólica en vistas a
ambientales que ocurren a nivel global [1]. En el incentivar futuras investigaciones y promover la
caso del gas natural, se ha observado que la inversión en el campo de la energía eólica.
electricidad generada por las centrales de ciclo
combinado, presenta costos totales mayores respecto REFERENCIAS
a la salud humana y externalidades ambientales, en [1] Muro y Fernández de la Oz, 2013
comparación con la energía eólica. [2]. Para revertir [2] McCubbin y Sovacool, 2013
el impacto ambiental producido por la generación de [3] Kumar et al., 2016
electricidad se han implementado centrales [4] Ministerio de Energía y Minería del Perú.
eléctricas que emplean fuentes renovables de (2016). “Atlas Eólico del Per ‘u”.
energía. La energía eólica es una fuente de energía [5] Comité de operación económica del sistema
renovable, limpia, abundante y con un alto potencial interconectado nacional. (2016). “Informe de la
a nivel global. Para transformar la energía eólica en operación anual del SEIN”.
electricidad se construyen parques eólicos. Un
parque eólico (central o granja eólica) es el sitio II. METODOLOGÍA
donde se instala un grupo de turbinas eólicas para
generar electricidad. En los últimos años, la Para el desarrollo de la investigación se realizó en
demanda mundial de energía renovable, primera instancia el análisis y simulación
principalmente la energía eólica, ha aumentado [3]. aerodinámica del perfil Wortmann FX 60-126, así
(Kumar et al., 2016) como también el análisis y simulación aerodinámica
En nuestro país, desde mediados de los años ochenta, de la turbina eólica Waira 12.
se realizaron las primeras instalaciones de turbinas
eólicas. En Villa El Salvador se instaló una turbina 1. ANÁLISIS Y SIMULACIÓN
prototipo de 500 W a cargo de ITINTEC, posterior a AERODINÁMICA DEL PERFIL
ello, ELECTROPERU gestionó las centrales eólicas WORTMANN FX 60-126
de Yasila y Toran con turbinas eólicas con capacidad
El perfil aerodinámico seleccionado en el diseño de
de 3.5, 4 y 12 kW. En los años noventa se instalaron
la turbina eólica es el FX 60 − 126, diseñado en el
dos turbinas eólicas, una en Malabrigo de 250 kW y
año 1960 y pertenece a la familia de perfiles
otra en Marcona de 450kW, siendo las de mayor
elaborados por el Dr. Franz Xavier Wortmann.
potencia hasta la fecha [4]. En el 2016 la generación
de energía eléctrica a nivel nacional fue de 48 326.4 Los coeficientes aerodinámicos experimentales
GWh, de la cual 2 373.4 GWh provino de recursos fueron obtenidos del Sttutgarter Profilkataloge I [25]
energéticos renovables (RER), representando el desarrollado por el Instituto de Aerodinámica y Gas-
4.91% del total. Dentro de ella, la energía eólica dinámica de la Universidad de Sttutgart.
representó el 2.18% del total generado, que equivale
a 1054.1 GWh. [5] También se obtuvo los coeficientes aerodinámicos
A causa de las bondades que representa el uso de las de sustentación y arrastre para diversos números de
turbinas eólicas como fuente de producción de Reynolds mediante la aplicación del Software XFoil,
energía, el presente trabajo busca profundizar en el para lo cual se consideró un flujo incompresible,
estudio del rotor de la turbina eólica en vistas a también la viscosidad del aire que está a una
incentivar futuras investigaciones y promover la temperatura de 25° C y una intensidad turbulenta de
inversión en el campo de la energía eólica. 1, la cual representa una alta intensidad de
turbulencia.
Finalmente, se realizó la simulación fluido dinámica Para obtener los coeficientes de sustentación y
del perfil. Para lo cual, dado el perfil definido, se arrastre se simularon las 3 mallas bajo las mismas
procedió a definir un volumen de control usando la condiciones de operación de numero de Reynolds
herramienta SpaceClaim. 1.5 x 105 y para el ángulo de ataque α = 4°.

Figura 1: Dimensión del volumen de control.

Luego para generar las diversas mallas se utilizó la
herramienta ICEM ANSYS, las cuales fueron Figura 4: Condiciones de frontera del perfil
orientadas a mantener un skewness menor a 0.9, se aerodinámico.
generaron 3 mallas, gruesa, media y fina. 2. ANÁLISIS Y SIMULACIÓN
AERODINÁMICA DE LA TURBINA
WAIRA 12
La turbina eólica Waira 12 diseñada y gestionada por
la empresa Waira Energía SAC y financiada por
Innóvate Perú y con la colaboración de la Pontificia
Universidad Católica del Perú- PUCP, Redondos
S.A. y Fundición Ferrosa S.A.C. pertenece al
proyecto 217 FINCyT - FIDECOM PIPEA 201. Los
parámetros de diseño proporcionados por la empresa
Figura 2: Estructuración de las mallas. Waira Energía SAC se presentan en la tabla 1.

Tabla 1: Parámetros de diseño turbina Waira 12.

Figura 3: Malla alrededor del perfil aerodinámico.
Para el análisis de la turbina se seguirá el modelo
En las condiciones de frontera se indicó como fluido BEM (Blade Element Momentum), este modelo
de trabajo el aire a condiciones ambientales, Air at calcula la potencia y el empuje generado por una
25°C en el software Ansys CFX, la densidad y la turbina eólica de eje horizontal ingresando los
viscosidad también están establecidas a la parámetros de diseño.
temperatura de 25°C. El modelo de turbulencia
seleccionado para la simulación del perfil es el κ−ω Validación del modelo BEM
SST, el cual es una versión mejorada del modelo κ − La validación del código se realizó comparando los
ω de Wilcox. [21]. El número de Reynolds es un resultados del modelo de BEM con los obtenidos en
factor que varía en las diferentes simulaciones, los las siguientes publicaciones las cuales
cuales serán: 7 x 105, 1 x 106, 1.5 x 106 y 2 x 106. denominaremos: OPT HAWT NACA 4418 - The
Y la longitud característica del perfil aerodinámico Performance Test of Three Different Horizontal
c=1. Otro parámetro que variará en las diferentes Axis Wind Turbine [7] y OPT HAWT NACA 4412
simulaciones es el ángulo de ataque, los cuales serán: - Small Wind Turbines [8]
0, 2, 4, 6, 8, 10, 12, 14, 16, 18. La presión de
referencia en esta investigación será la atmosférica,  La turbina eólica descrita por Hsiao et al [7], la
así se define la presión relativa como 0 Pa. cual denominaremos OPT HAWT NACA 4418,
Finalmente, la rugosidad del perfil aerodinámico la emplea el perfil aerodinámico NACA 4418 a lo
cuál es 0 debido a que tiene la condición de no largo de toda la pala. Para el cálculo de potencia
deslizamiento y su superficie se le considera suave. y empuje mediante el modelo de BEM se usaron
las condiciones de operación experimentales
 descritas por Hsiao [7], las cuales son V = 10m/s Debido al recurso computacional y a
y λ = 4.92. recomendaciones bibliográficas [16] la simulación
se realizó tomando un tercio del rotor. Luego se
Los coeficientes aerodinámicos requeridos para procedió a generar los volúmenes de control
el código se obtuvieron del XFoil para una (dominio rotatorio y dominio estacionario)
intensidad de turbulencia de 1, el cual cuenta
con el perfil NACA 4418 en su base de datos Después de diversos análisis se decidió usar el
para generar la información aerodinámica dominio estacionario tipo difusor, para ahorrar
necesaria. Finalmente se procesan los resultados recurso computacional. Finalmente, a través de las
y se obtuvieron los valores de la potencia y del herramientas de revolución y operaciones booleanas
empuje de la turbina en estudio haciendo uso de obtenemos los volúmenes de control dominio
las teorías de Hansen [2] y Wood [8]. rotatorio (Figura 7a) y dominio estacionario (Figura
7b).
 La turbina eólica descrita en Wood [8], la cual
denominaremos OPT HAWT NACA 4412,
emplea a lo largo de toda su pala el perfil
aerodinámico NACA4412. Para el cálculo de la
potencia y del empuje mediante el modelo de
BEM se usaron las condiciones de operación
descritas por Wood [8], donde V = 10 m/s y
TSR = 10 y los coeficientes aerodinámicos
obtenidos por el Xfoil para una intensidad de
turbulencia de 1. Figura 7: Volúmenes de Control.

Simulación fluido dinámica de la turbina eólica De manera similar a la generación de mallas para el
Waira 12 perfil Wortman Fx 60-126 se crearon 4 mallas
diferentes: mallas basta, gruesa, media y fina. Estas
De acuerdo a la descripción de la turbina eólica se mallas son no estructuradas debido a la complejidad
procedió a la modelación. Se inició con la creación de la pala.
de las estaciones de la pala de la turbina eólica y en
cada estación se insertó el perfil Wortman FX 60- En las condiciones de frontera se indicó que el
126. Posteriormente estos perfiles fueron escalados modelo de turbulencia seleccionado es el κ − ω SST
y rotados para las magnitudes de cuerda y torsión y las propiedades del aire a 25◦ C son ρ = 1.185
descritas en la Tabla 1. kg/m3 y v = 1.54 x 10-5 m2/s. En esta sección no
utilizaremos el número de Reynolds, las variables
serán la velocidad de flujo libre (V) y la velocidad
de rotación de la turbina eólica (ω). Los valores de ω
son: 40, 60, 80, 100, 120, 140, 160, 180 RPM. La
presión de referencia y rugosidad será análoga a las
definidas para el perfil.
Debido a la selección de simular un tercio del rotor,
las caras que se encuentran a 120° se definieron con
la condición de periodicidad, así como se aprecia en
Figura 5: Estaciones de la pala de la turbina eólica. la figura 9.

Posteriormente, se optó por tomar la sección del

perfil aerodinámico escalado para que mantenga la
máxima proximidad a la pala generada por la
empresa Waira Energía. Se finalizó el modelado con
la creación del Hub, donde obtuvimos finalmente el
rotor de la turbina eólica en la Figura 6.

Figura 8: Condiciones de frontera dominio rotatorio.

Figura 6: Vista trimétrica de la turbina eólica.

Figura 9: Condiciones de frontera dominio
estacionario.
Para obtener los resultados de potencia y empuje se
simularon las 4 mallas bajo las condiciones de
diseño presentadas anteriormente donde la velocidad
de flujo del aire es V = 10.5 m/s y la velocidad
rotacional es ω = 120 RPM. Y luego se procede a
encontrar los errores relativos en las diferentes
mallas.
Finalmente se elige la malla que tenga las
distribuciones de presión y esfuerzo cortante más
próximas a los resultados del perfil bidimensional.

III. RESULTADOS
IV. CONCLUSIONES
 Los coeficientes aerodinámicos obtenidos de la
simulación mediante el software Ansys CFX
son próximos a los valores obtenidos por el
XFoil y también a los datos de la fuente
experimental en el rango de ángulos de ataque
de 0◦ a 12◦.
 Se concluye que el método de BEM es una
buena herramienta para la fase de diseño
debido a que el máximo punto de extracción de
potencia se encuentra aproximadamente en su
punto de diseño en la simulación
computacional.
 Los resultados obtenidos por la simulación
computacional son suficientemente confiables
debido a que no están sujetos a
simplificaciones del modelo de BEM, el cual
no toma en cuenta la influencia radial y los
efectos tridimensionales en su proceso de
iteración.
 Considerando una eficiencia mecánica de 0.99
y una eficiencia eléctrica de 0.9, las potencias
obtenidas mediante el modelo de BEM y la
simulación computacional se vuelven
respectivamente Pe BEM = 31.5 kW y Pe CFX
= 31.7 kW, donde ambas son mayores a la
potencia nominal de 25 kW del Waira 12.

V. REFERENCIAS BIBLIOGRÁFICAS