Direct Numerical Simulation of Vertical-Axis Wind Turbine Near-Wake Dynamics

Este estudio emplea simulaciones numéricas directas con resolución geométrica para revelar que aumentar el número de álabes en turbinas eólicas de eje vertical acelera la descomposición de los vórtices de estallido dinámico mediante interacciones álabe-vórtice, provocando que la estela cercana transicione más rápidamente hacia la dinámica de cuerpos romos y demostrando que el número de álabes, en lugar de la relación de velocidad en la punta, es el factor principal que gobierna esta transición y las características del flujo de entrada aguas abajo.

Lo esencial

Imagina una turbina eólica de eje vertical (VAWT) no como una máquina gigante, sino como un molinillo de viento girando en posición vertical dentro de un río de aire. Este artículo es como una película de alta definición en cámara lenta que hace zoom tan de cerca sobre las palas que podemos ver las «corrientes de aire» invisibles que giran a su alrededor. Los investigadores utilizaron una supercomputadora para simular exactamente cómo se comporta el aire justo detrás de las palas giratorias, un lugar llamado «estela cercana».

Aquí está la historia de lo que descubrieron, explicada de forma sencilla:

La Gran Imagen: ¿Por qué observar la «estela cercana»?

La mayoría de las turbinas eólicas giran como una hélice de avión (eje horizontal). Pero estas de eje vertical giran como un centrifugador de ensalada. Los investigadores querían saber: ¿Qué le sucede al aire inmediatamente después de pasar por las palas?

Descubrieron que el aire no fluye suavemente detrás de la turbina. En cambio, se vuelve caótico. Las palas cortan el aire, creando remolinos gigantes y giratorios llamados Vórtices de Pérdida Dinámica (DSV). Imagina estos como remolinos gigantes e invisibles que son levantados por las palas y se arrastran detrás de la turbina.

El Descubrimiento Principal: La «Cantidad de Palas» Importa Más que la Velocidad

El equipo probó turbinas con una pala, dos palas y tres palas. También cambiaron la velocidad a la que giraban las turbinas.

Aquí está el giro sorprendente: Cuántas palas tiene la turbina importa mucho más que qué tan rápido gira.

• La Turbina de Una Pala (El Solista):
  Imagina a un bailarín solitario girando en una habitación. Cuando gira, crea un remolino gigante y poderoso de aire detrás de él. Este remolino gigante (el DSV) se mantiene fuerte e intacto durante mucho tiempo. Es como una nube pesada y lenta que tarda mucho en disiparse. Debido a que esta nube gigante permanece, el aire detrás de la turbina se mantiene «desordenado» y caótico durante una gran distancia.

• La Turbina de Tres Palas (El Trío):
  Ahora imagina a tres bailarines girando muy cerca entre sí. Mientras un bailarín gira, el aire que acaba de perturbar (el remolino) es inmediatamente golpeado por el siguiente bailarín.
  Los investigadores descubrieron un nuevo mecanismo al que llaman «Impacto de Vórtice».

  • La Analogía: Imagina el remolino gigante como una burbuja de jabón. En la turbina de una pala, la burbuja se aleja flotando intacta. En la turbina de tres palas, la siguiente pala actúa como un alfiler, estallando la burbuja prematuramente.
  • El Resultado: El remolino gigante y desordenado se rompe en burbujas diminutas e inofensivas (vórtices más pequeños) antes de poder viajar lejos. El aire detrás de la turbina de tres palas se vuelve «tranquilo» y ordenado mucho más rápido que la versión de una pala.

El Efecto de «Atasco de Tráfico»

El artículo también explica que tener más palas crea un poco de un «atasco de tráfico» para el viento.

• Con más palas, hay menos espacio abierto para que el viento fluya a través del centro de la turbina.
• Esto obliga al viento a ir alrededor de la turbina, como el agua que fluye alrededor de una roca en un arroyo.
• Esto cambia el comportamiento de la estela. En lugar de estar dominada por los remolinos giratorios gigantes (que ocurren con pocas palas), la estela comienza a comportarse como la estela detrás de una roca sólida simple (un «cuerpo romo»).
• Por qué esto es bueno: La estela «tipo roca» se recupera (vuelve a ser aire suave) mucho más rápido que la estela «tipo remolino».

La Prueba de «Auto-similitud»

Los investigadores querían saber exactamente cuándo el aire desordenado detrás de la turbina vuelve a convertirse en aire suave y predecible. Utilizaron un truco matemático llamado «análisis de auto-similitud».

• La Analogía: Imagina intentar predecir la forma de una columna de humo. Al principio, el humo es un caos desordenado de remolinos. Pero eventualmente, se asienta en una forma predecible, con forma de curva de campana.
• Descubrieron que las turbinas de tres palas se asientan en esta forma predecible mucho antes que las turbinas de una pala. El «desorden» desaparece más rápido.

Qué Significa Esto para el Futuro (Según el Artículo)

El artículo menciona específicamente que estos hallazgos son importantes para parques eólicos donde las turbinas están colocadas muy cerca entre sí.

• Si colocas una segunda turbina justo detrás de una primera, el aire que «respira» depende en gran medida de cuántas palas tenga la primera turbina.
• Si la primera turbina tiene pocas palas, la segunda turbina recibe el impacto de remolinos gigantes y caóticos.
• Si la primera turbina tiene muchas palas, el aire ya se ha «curado» y se ha vuelto más suave para cuando llega a la segunda turbina.

En resumen: Este artículo utilizó una simulación por computadora superpotente para mostrar que añadir más palas a una turbina eólica de eje vertical actúa como un «rompe-remolinos». Destroza los vórtices de aire gigantes y desordenados en pedazos diminutos, permitiendo que el viento recupere su flujo suave mucho más rápido. Esto es crucial para diseñar parques eólicos donde las turbinas están apretadas muy juntas.

Resumen técnico

Resumen Técnico: Simulación Numérica Directa de la Dinámica de la Estela Cercana de Turbinas Eólicas de Eje Vertical

Planteamiento del Problema
Mientras que las Turbinas Eólicas de Eje Horizontal (HAWT) dominan la investigación actual, las Turbinas Eólicas de Eje Vertical (VAWT) ofrecen ventajas distintivas para aplicaciones offshore, urbanas y remotas, incluyendo centros de gravedad más bajos e insensibilidad direccional omnidireccional. Un beneficio teórico clave de las VAWT es su potencial para lograr mayores densidades de potencia en arreglos de espaciamiento cercano debido a regiones de estela más cortas y una recuperación de estela mejorada. Sin embargo, la estela cercana de una VAWT está dominada por características de flujo complejas asociadas con la pérdida aerodinámica dinámica, específicamente la formación y el desprendimiento de Grandes Vórtices de Pérdida Aerodinámica Dinámica (DSV). A diferencia de la estela lejana, que se comporta como una estela de cuerpo romo, la estela cercana depende altamente de parámetros geométricos y operativos, como la relación de velocidad de punta ($\lambda$) y la solidez estática ($\sigma$). Los datos experimentales existentes a menudo carecen de resolución cerca de la geometría o en tres dimensiones, mientras que las Simulaciones de Grandes Remolinos (LES) de alta fidelidad acopladas con el Método de Línea de Actuador (ALM) no resuelven geométricamente la capa límite de la pala. En consecuencia, el ALM no puede capturar la dinámica de la burbuja de separación laminar (LSB), la formación precisa del DSV ni las interacciones críticas pala-vórtice (BVI) que modulan estas características.

Metodología
Para abordar estas limitaciones, los autores presentan Simulaciones Numéricas Directas (DNS) geométricamente resueltas de VAWT tipo Darrieus de palas rectas utilizando el marco del método espectral/hp, Nektar++. El estudio investiga tres configuraciones de turbina con una, dos y tres palas (perfiles alares NACA 0018) en un rango de relaciones de velocidad de punta ($\lambda = 1.5, 2.5, 4.5$). Este espacio de parámetros corresponde a solididades estáticas ($\sigma$) de 0.2, 0.4 y 0.6, y una solidez dinámica ($\sigma_D$) que varía de 0.47 a 0.94.

Las simulaciones emplean una formulación de Marco de Referencia Móvil (MRF) para manejar la geometría rotatoria sin remallado dinámico. El flujo se modela mediante las ecuaciones de Navier-Stokes incompresibles a un número de Reynolds basado en la cuerda de $Re_c = 10^4$, un valor elegido para representar entornos urbanos/extraterrestres y asegurar la viabilidad numérica manteniendo la relevancia para la física fundamental. La resolución de la malla satisface los requisitos de DNS ($\Delta y^+ \lesssim 1$) con 32 planos de Fourier en la dirección de la envergadura. Se aplica Viscosidad de Disipación Espectral (SVV) para la estabilidad numérica en regiones subresueltas. La metodología se valida contra datos experimentales de Battisti et al., mostrando un fuerte acuerdo en los perfiles de velocidad media en el tiempo en la dirección del flujo en ubicaciones de estela cercana.

Contribuciones y Resultados Clave

1. Resolución de la Pérdida Aerodinámica Dinámica y las BVI:
   La DNS captura con éxito el ciclo completo de pérdida aerodinámica dinámica, incluido el enrollamiento de la capa límite, la formación del DSV, la separación y la convección. Un hallazgo novedoso es la identificación de un mecanismo de impacto de vórtice impulsado por las BVI.

   • Configuración de una Pala: El DSV se forma, separa y convecta hacia la estela cercana como una estructura coherente a gran escala, causando una aceleración significativa del flujo en el lado de sotavento.
   • Configuración de dos Palas: La solidez aumentada reduce la convección del DSV, provocando que se difunda y se descomponga en estructuras más pequeñas antes de alcanzar la estela cercana.
   • Configuración de tres Palas: La reducción del espaciamiento entre palas introduce BVI al inicio del paso de barlovento. Estas interacciones ocurren durante la fase temprana de la pérdida aerodinámica dinámica (desarrollo de la burbuja de separación laminar), perturbando la LSB y reduciendo el tamaño eventual del DSV. Además, a medida que la pala se acerca a la zona de sotavento, el impacto directo de vórtices procedentes de estelas aguas arriba causa la ruptura prematura del DSV en múltiples núcleos distintos. Este mecanismo acelera la transición de la dinámica de vórtices coherentes a la dinámica de cuerpo romo.

2. Recuperación de la Estela y Transición a la Dinámica de Cuerpo Romo:
   El estudio cuantifica la evolución aguas abajo de la estela. Una mayor solidez (más palas) resulta en déficits de velocidad iniciales más severos, pero inicia la recuperación antes y más rápidamente.

   • Baja Solidez ($N_{palas}=1$): La recuperación está gobernada por la mezcla transversal inducida por grandes DSV, lo que lleva a una recuperación más lenta y a una asimetría persistente de la estela.
   • Alta Solidez ($N_{palas}=3$): La descomposición de los DSV impulsada por las BVI elimina las estructuras coherentes que distinguen la estela de la VAWT de un cuerpo romo. En consecuencia, la estela transita más rápidamente hacia mecanismos de recuperación asociados a la capa de corte, que son intrínsecamente más rápidos.

3. Análisis de Auto-similitud:
   Extendiendo el análisis de auto-similitud a la estela cercana (anteriormente limitado a la estela lejana en estudios con ALM), los autores cuantifican la tasa a la cual la estela pierde la dependencia de las estructuras coherentes generadas por las palas y colapsa sobre una solución gaussiana auto-similar.

   • El análisis revela que el número de palas es más influyente que la relación de velocidad de punta para establecer la tasa de esta transición.
   • Mayores conteos de palas conducen a un inicio más rápido de la auto-similitud. Curiosamente, la tasa de transición alcanza su máximo en $\lambda=2.5$ en lugar del $\lambda$ más alto, lo que sugiere un efecto de meseta más allá de cierta solidez dinámica.

Significado y Afirmaciones
El artículo afirma que el enfoque de DNS geométricamente resuelto es esencial para identificar los mecanismos específicos mediante los cuales el número de palas influye en la aerodinámica de las VAWT. El significado principal radica en el descubrimiento de la ruptura prematura del DSV inducida por BVI, un mecanismo que acelera la transición a la dinámica de estela de cuerpo romo independientemente de los efectos de bloqueo.

Los autores afirman que estos hallazgos tienen implicaciones directas para el diseño de arreglos de turbinas de espaciamiento cercano y configuraciones de pares acoplados. Dado que el flujo de entrada experimentado por un rotor aguas abajo depende fundamentalmente del número de palas de la turbina aguas arriba, los pares con bajo número de palas encontrarán grandes DSV coherentes, mientras que los pares con alto número de palas encontrarán un campo más difuso y perturbado. El artículo concluye que las DNS de configuraciones emparejadas están justificadas para comprender completamente cómo estas condiciones de entrada distintas afectan el rendimiento aguas abajo, señalando que la recuperación de estela más rápida y el inicio más temprano de la auto-similitud a mayores solididades pueden permitir una modelización analítica más precisa de dichos arreglos a distancias aguas abajo más cortas.