Numerical simulations of transition and long-term response of a wind turbine airfoil

Este artículo presenta simulaciones numéricas utilizando Nek5000 y EllipSys para analizar la transición y la respuesta del flujo a largo plazo de un perfil alar FFA-W3 a números de Reynolds bajos, validando la capacidad del código EllipSys para capturar inestabilidades de Kelvin-Helmholtz e identificando una modulación lenta y periódica del coeficiente de fuerza normal vinculada a oscilaciones de baja frecuencia y al posible estallido de burbujas.

Lo esencial

Imagina una pala de aerogenerador como un ala gigante que gira. Al igual que el ala de un avión, necesita un flujo de aire suave para funcionar de manera eficiente. Pero cuando el viento la golpea desde ciertos ángulos, el aire puede "atascarse" y separarse de la superficie, creando un desorden caótico y arremolinado. Este artículo es como un experimento de túnel de viento de alta tecnología, pero en lugar de usar un modelo físico, los investigadores construyeron uno virtual dentro de una supercomputadora para observar exactamente cómo se comporta este aire.

Aquí está la historia de su descubrimiento, desglosada en conceptos simples:

1. El túnel de viento virtual

Los investigadores quisieron estudiar una sección específica de una enorme pala de aerogenerador (de una turbina de 10 MW). Utilizaron dos programas informáticos diferentes, NEK5000 y ELLIPSYS, para simular el flujo de aire sobre esta pala.

Piensa en NEK5000 como una cámara de alta gama, ultra precisa, que captura cada pequeño detalle pero es muy lenta y costosa de ejecutar. ELLIPSYS es como una cámara ligeramente más rápida y eficiente. El equipo primero tuvo que demostrar que la cámara "más rápida" (ELLIPSYS) podía ver lo mismo que la de "alta gama". Encontraron que, si bien ELLIPSYS perdía algunas ondulaciones diminutas y tenues en el aire suave (porque suavizaba las cosas un poco de más), era excelente para capturar los grandes remolinos caóticos que realmente importan para el rendimiento de la pala.

2. ¿Qué tan ancho debe ser el túnel?

Antes de ejecutar las simulaciones largas, tuvieron que determinar qué tan ancho debía ser su "túnel de viento" virtual. Si el túnel es demasiado estrecho, podría apretar el aire y crear resultados falsos. Si es demasiado ancho, desperdicia potencia de cómputo.

Probaron un túnel "estrecho" (10% del ancho del ala) contra un túnel "ancho" (20% del ancho).

• La analogía: Imagina que observas el flujo de un río. Si solo miras una franja estrecha del río, ¿te pierdes las grandes olas?
• El resultado: Encontraron que el túnel estrecho era en realidad suficiente. Las grandes olas y remolinos se formaron perfectamente en el espacio estrecho. Esto significó que podían ahorrar mucho tiempo de computación usando la caja de simulación más pequeña y estrecha.

3. La "burbuja" y la "aleta"

La parte más interesante del estudio ocurrió en la parte superior del ala (el lado de succión).

• La burbuja de separación: A medida que el aire fluye sobre el ala, se desprende por un momento, creando una pequeña bolsa de aire recirculante llamada "Burbuja de Separación Laminar" (LSB, por sus siglas en inglés). Piensa en esto como un pequeño remolino temporal en la superficie del ala.
• La inestabilidad: Dentro de esta burbuja, el aire no se queda quieto; vibra y se enrolla en ondas (como las ondas en un estanque). Los investigadores observaron cómo crecían estas ondas. Descubrieron que el principal "rodillo" en esta burbuja es un tipo de inestabilidad llamada modo Kelvin-Helmholtz.
• El descubrimiento: Confirmaron que el programa informático "más rápido" (ELLIPSYS) podía predecir con precisión cómo crecían estas ondas y cómo se comportaba la burbuja, coincidiendo con los resultados del programa ultra preciso.

4. El pulso lento y rítmico (La gran sorpresa)

Después de validar sus herramientas, dejaron que la simulación corriera durante muchísimo tiempo (equivalente a 50 veces el paso del aire frente al ala). Aquí es donde encontraron algo especial.

Mientras el aire se agitaba con movimientos rápidos y caóticos, notaron un pulso rítmico muy lento en la fuerza que empuja el ala.

• La analogía: Imagina un redoble de tambor. El movimiento rápido y caótico del aire es como un redoble de tambor rápido y agudo. El pulso lento que encontraron es como un latido profundo y lento que ocurre una vez cada 48 segundos (en el tiempo de la simulación).
• El efecto: Este latido lento hizo que la fuerza sobre el ala oscilara hacia arriba y hacia abajo aproximadamente un 10.5%.
• La conexión con los aerogeneradores reales: Cuando tradujeron esto a un aerogenerador real que gira, se dieron cuenta de que este pulso lento ocurre una vez cada 7.7 rotaciones completas de la pala.

5. ¿Por qué sucede esto?

Los investigadores creen que este pulso lento es causado por un ciclo de "entrada en pérdida" (el aire se queda atrapado o stalling) y "salida de pérdida" (el aire se libera o unstalling) en el ala.

• El ciclo: El aire se queda atrapado, creando una gran burbuja. Luego, algo provoca que la burbuja estalle y el aire se adhiera suavemente de nuevo. Luego, la presión aumenta de nuevo, la burbuja se forma y el ciclo se repite.
• El detonante: Sospechan que esto sucede porque el aire gira hacia atrás con tanta fuerza en el ala que crea un estado de "inestabilidad absoluta", una forma elegante de decir que el aire es tan turbulento que no puede evitar oscilar por sí mismo.

6. La conclusión fundamental

Este estudio es una historia de éxito para el modelado por computadora. Demostraron que un programa informático más rápido y eficiente (ELLIPSYS) puede ser confiable para estudiar la compleja física de los aerogeneradores, siempre que se verifique primero contra el "estándar de oro".

Descubrieron que, incluso en una pala de aerogenerador gruesa, existe un "respiro" de flujo de aire muy lento y rítmico que ocurre aproximadamente cada 8 rotaciones de la pala. Este respiro hace que la sustentación (la fuerza que hace girar la turbina) suba y baje significamente. Comprender este ritmo lento es crucial porque, aunque no romperá la turbina inmediatamente, estos cambios lentos y grandes en la fuerza podrían fatigar los materiales tras muchos años de operación.

En resumen: Construyeron un túnel de viento virtual, demostraron que una computadora rápida podía hacer el trabajo y descubrieron que las palas de los aerogeneradores tienen un "latido" lento y rítmico causado por la formación y explosión de burbujas de aire en su superficie.

Resumen técnico

Resumen Técnico: Simulaciones Numéricas de la Transición y Respuesta a Largo Plazo de un Perfil Aerodinámico de Aerogenerador

Planteamiento del Problema
El desprendimiento del flujo en las palas de los aerogeneradores, impulsado por ángulos de incidencia elevados y geometrías de perfiles gruesos, impacta significativamente el rendimiento y la carga estructural. Estudios previos han identificado múltiples frecuencias naturales en flujos desprendidos, incluyendo el desprendimiento de estelas, modos Kelvin-Helmholtz (KH) y oscilaciones de baja frecuencia (LFOs). Las LFOs se caracterizan por fuertes variaciones de sustentación y a menudo se atribuyen al estancamiento y destancamiento periódico del flujo, potencialmente vinculados a la inestabilidad absoluta dentro de burbujas de separación laminar (LSBs). Si bien las LFOs han sido observadas en diversos estudios de perfiles aerodinámicos, su ocurrencia a bajos ángulos de ataque (AoA) en perfiles de aerogeneradores gruesos, y los mecanismos específicos que las impulsan en estas condiciones, requieren mayor investigación. Este estudio se centra en la transición y la respuesta a largo plazo de una sección del perfil aerodinámico de la turbina de referencia DTU 10-MW (serie FFA-W3) a un número de Reynolds de cuerda ($Re_c$) de $1 \times 10^5$.

Metodología
La investigación emplea simulaciones de grandes remolinos (LES) con resolución de pared sin modelos de subrejilla para resolver las ecuaciones de Navier-Stokes incompresibles. Se utilizan dos solvers:

1. NEK5000: Un código de método de elementos espectrales (SEM) conocido por su mínima disipación numérica y alta precisión.
2. ELLIPSYS: Un código de volúmenes finitos desarrollado en DTU/Risø, tradicionalmente utilizado para RANS, pero aplicado aquí como LES implícito.

El estudio implica una validación sistemática de ELLIPSYS frente a NEK5000. Los pasos metodológicos clave incluyen:

• Sensibilidad del Dominio: Evaluación del ancho del dominio en la dirección de la envergadura ($L_z$) comparando simulaciones con $L_z = 0.1c$ y $L_z = 0.2c$ para asegurar la resolución de las estructuras de flujo relevantes.
• Validación del Solver: Comparación de los campos de flujo medio, la evolución de las perturbaciones y las características espectrales entre ambos solvers.
• Análisis de Estabilidad: Utilización de Ecuaciones de Estabilidad Parabolizadas (PSE) y Descomposición Ortogonal Propia Espectral (SPOD) para analizar el crecimiento de las inestabilidades (modos de Tollmien-Schlichting y KH) y la estructura espacial de los modos coherentes.
• Integración a Largo Plazo: Uso del solver ELLIPSYS, computacionalmente más eficiente, para simular la evolución del flujo durante 50 pasos de paso ($T = 50c/U_\infty$) para capturar fenómenos de baja frecuencia.

Resultados Clave

• Ancho del Dominio: Se determinó que un ancho de envergadura del 10% de la cuerda ($L_z = 0.1c$) es suficiente para capturar la evolución de las principales perturbaciones y reproducir las estadísticas de flujo promedio, coincidiendo con los resultados de un dominio más ancho ($L_z = 0.2c$).
• Validación del Solver: ELLIPSYS demostró una estrecha concordancia con NEK5000 en los campos de flujo medio y la evolución de las perturbaciones más amplificadas. Sin embargo, ELLIPSYS subestimó la amplitud de las ondas de Tollmien-Schlichting (TS) en la capa límite adherida debido a una mayor disipación numérica. Por el contrario, ELLIPSYS predijo con precisión la evolución del modo KH dentro de la burbuja de separación laminar (LSB), alineándose bien con las predicciones de PSE.
• Mecanismo de Transición: La transición en el lado de succión es impulsada por la inestabilidad KH que se desarrolla en la capa de cizalladura separada. El modo KH más amplificado tiene una frecuencia de $f \approx 15$, y su evolución aguas abajo es capturada adecuadamente por ambos solvers y por PSE.
• Oscilaciones de Baja Frecuencia (LFOs): Las simulaciones a largo plazo revelaron una modulación lenta del coeficiente de fuerza normal ($C_N$) con una amplitud del 10.5% y un periodo de 48 pasos de paso. Esto corresponde a una frecuencia $f = 0.021$ y un número de Strouhal $St = 0.0012$.
  • Esta frecuencia es aproximadamente 720 veces más baja que la frecuencia de la inestabilidad KH primaria.
  • En el contexto de la turbina DTU 10-MW, este periodo equivale a aproximadamente 7.7 rotaciones de la pala.
  • El $St$ observado es inferior a los valores reportados típicamente en la literatura para otros perfiles, ocurriendo en un ángulo de ataque efectivo relativamente pequeño de $3.1^\circ - 3.3^\circ$.

Significancia y Reivindicaciones
El artículo establece que ELLIPSYS es una herramienta viable para LES de pared de flujos transicionales, ofreciendo una alternativa computacionalmente eficiente a los códigos de elementos espectrales para integraciones de largo tiempo. El estudio confirma que el estancamiento y destancamiento periódico del flujo, potencialmente desencadenado por la inestabilidad absoluta debido al alto flujo inverso (hasta -17% en el lado de succión), es responsable de las LFOs observadas.

Un hallazgo significativo es la identificación de LFOs en un perfil grueso (24% de espesor) a un bajo ángulo de ataque, una condición distinta de los estancamientos de perfiles delgados típicamente asociados con tales fenómenos en la literatura. Los autores señalan que el flujo inverso en ambos lados del perfil es suficiente para permitir la inestabilidad absoluta, lo que puede impulsar el estallido periódico de la LSB hacia un flujo totalmente estancado. El estudio contribuye a la comprensión de cómo la geometría de un perfil grueso influye en el inicio y las características de las oscilaciones de baja frecuencia, destacando que estos fenómenos pueden ocurrir a ángulos de incidencia y números de Strouhal más bajos de lo documentado previamente para perfiles estándar.