Determinación de las características aerodinámicas de perfiles para aerogeneradores de baja potencia utilizando Dinámica de Fluidos Computacional

LÓPEZ, EZEQUIEL; ZAPPA, ANDRÉS; WILD CAÑÓN, CARLOS

La Plata

Congreso; XXIII Congreso sobre Métodos Numéricos y sus Aplicaciones ENIEF 2017; 2017

Facultad de Ingeniería, Universidad Nacional de La Plata - Asociación Argentina de Mecánica Computacional

La teoría de elemento de pala y conservación de la cantidad de movimiento (BEM, por Blade Element-Momentum Theory) constituye la base usualmente empleada para el diseño de Aerogeneradores de Baja Potencia. Esta teoría requiere del conocimiento a priori de los coeficientes de sustentación y arrastre de los perfiles aerodinámicos que conforman la sección de las palas, los cuales suelen obtenerse de forma experimental. Cuando no se dispone de datos experimentales, una opción es utilizar técnicas CFD (Computational Fluid Dynamics) a fin de obtener las características aerodinámicas del perfil mediante la simulación del flujo en el rango de operación. Las palas de los Aerogeneradores de Baja Potencia trabajan en un régimen de flujo considerado de bajo número de Reynolds, típicamente por debajo de 5e5. En este régimen es posible que la capa límite laminar se separe y, bajo ciertas circunstancias, vuelva a adherirse corriente abajo formando lo que se conoce como burbuja de separación laminar (LSB, por Laminar Separation Bubble). La formación de las LSBs puede provocar una disminución del coeficiente de sustentación del perfil y un aumento del coeficiente de arrastre, como así también generar un bucle de histéresis de estos coeficientes en la región del stall. La simulación por CFD de este tipo de flujos requiere por lo tanto del empleo de modelos de turbulencia que sean capaces de representar dichos fenómenos. En este trabajo se comparan dos modelos de turbulencia de tipo EVM (Eddy Viscosity Model) para las ecuaciones RANS (Reynolds-Average Navier-Stokes): una versión modificada del modelo SST (Shear Stress Transport) k-omega, en la cual se emplea la versión para bajos números de Reynolds del modelo k-omega, y el modelo transicional kT-kL-omega actualizado con las últimas correcciones introducidas por sus autores. Las modificaciones a ambos modelos fueron implementadas en la suite OpenFOAM(R), software que se utilizó para simular el flujo alrededor de un perfil Wortmann FX63-137 para ángulos de ataque comprendidos entre -6 y 12 grados, y números de Reynolds entre 1e5 y 5e5. Los resultados obtenidos son contrastados con datos experimentales y numéricos disponibles en la literatura.