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

Este artigo apresenta simulações numéricas utilizando o Nek5000 e o EllipSys para analisar a transição e a resposta de fluxo de longo prazo de um aerofólio FFA-W3 em baixos números de Reynolds, validando a capacidade do código EllipSys de capturar instabilidades de Kelvin-Helmholtz e identificando uma modulação lenta e periódica do coeficiente de força normal ligada a oscilações de baixa frequência e ao potencial rompimento de bolhas.

O essencial

Imagine uma pá de turbina eólica como uma asa gigante e giratória. Assim como a asa de um avião, ela precisa de um fluxo de ar suave para funcionar com eficiência. Mas, quando o vento atinge a pá em certos ângulos, o ar pode ficar "preso" e se separar da superfície, criando uma confusão caótica e turbulenta. Este artigo é como um experimento de túnel de vento de alta tecnologia, mas em vez de usar um modelo físico, os pesquisadores construíram um modelo virtual dentro de um supercomputador para observar exatamente como esse ar se comporta.

Aqui está a história da descoberta deles, dividida em conceitos simples:

1. O Túnel de Vento Virtual

Os pesquisadores queriam estudar uma fatia específica de uma enorme pá de turbina eólica (de uma turbina de 10 MW). Eles usaram dois programas de computador diferentes, NEK5000 e ELLIPSYS, para simular o fluxo de ar sobre essa pá.

Pense no NEK5000 como uma câmera de alto nível, ultraprecisa, que captura cada pequeno detalhe, mas é muito lenta e cara para operar. O ELLIPSYS é como uma câmera ligeiramente mais rápida e eficiente. A equipe primeiro teve que provar que a "câmera mais rápida" (ELLIPSYS) conseguia ver as mesmas coisas que a "de alto nível". Eles descobriram que, embora o ELLIPSYS perdesse algumas ondulações minúsculas e tênues no ar suave (porque suavizava demais as coisas), ele era excelente em capturar os grandes redemoinhos caóticos que realmente importam para o desempenho da pá.

2. Quão Largo o Túnel Precisa Ser?

Antes de rodar as simulações longas, eles tiveram que descobrir o quão largo o seu "túnel de vento" virtual precisaria ser. Se o túnel for muito estreito, ele pode espremer o ar e criar resultados falsos. Se for muito largo, desperdiça poder computacional.

Eles testaram um túnel "estreito" (10% da largura da asa) contra um túnel "largo" (20% da largura).

• A Analogia: Imagine observar um rio fluindo. Se você olhar apenas para uma faixa estreita do rio, você perde as grandes ondas?
• O Resultado: Eles descobriram que o túnel estreito era, na verdade, suficiente. As grandes ondas e redemoinhos se formaram perfeitamente no espaço estreito. Isso significou que eles puderam economizar muito tempo de computador usando a caixa de simulação menor e mais estreita.

3. A "Bolha" e a "Aba"

A parte mais interessante do estudo aconteceu no topo da asa (o lado de sucção).

• A Bolha de Separação: À medida que o ar flui sobre a asa, ele se desprende por um momento, criando uma pequena bolsa de ar recirculante chamada "Bolha de Separação Laminar" (LSB). Pense nisso como um pequeno redemoinho temporário na superfície da asa.
• A Instabilidade: Dentro dessa bolha, o ar não fica apenas parado; ele vibra e se enrola em ondas (como ondulações em um lago). Os pesquisadores observaram essas ondas crescerem. Eles descobriram que o principal "rolamento" nesta bolha é um tipo de instabilidade chamada modo Kelvin-Helmholtz.
• A Descoberta: Eles confirmaram que o programa de computador "mais rápido" (ELLIPSYS) podia prever com precisção como essas ondas cresciam e como a bolha se comportava, coincidindo com os resultados do programa ultrapreciso.

4. O Pulso Lento e Rítmico (A Grande Surpresa)

Após validar suas ferramentas, eles deixaram a simulação rodar por um tempo muito longo (equivalente a 50 vezes o fluxo de ar passando pela asa). Foi aqui que encontraram algo especial.

Enquanto o ar estava agitado com movimentos rápidos e caóticos, eles notaram um pulso rítmico muito lento na força que empurra a asa.

• A Analogia: Imagine uma batida de tambor. O turbilhonamento rápido do ar é como um rufar de tambor rápido e agudo. O pulso lento que eles encontraram é como um batimento cardíaco profundo e lento que acontece uma vez a cada 48 segundos (no tempo de simulação).
• O Efeito: Esse batimento cardíaco lento fez com que a força na asa oscilasse para cima e para baixo em cerca de 10,5%.
• A Conexão com Turbinas Reais: Quando traduziram isso para uma turbina eólica real, perceberam que esse pulso lento acontece uma vez a cada 7,7 rotações completas da pá.

5. Por Que Isso Acontece?

Os pesquisadores acreditam que este pulso lento é causado por um ciclo de o ar "estolar" (ficar preso) e "desestolar" (soltar-se) na asa.

• O Ciclo: O ar fica preso, criando uma grande bolha. Então, algo desencadeia o estouro da bolha, e o ar se reatacha suavemente. Em seguida, a pressão aumenta novamente, a bolha se forma e o ciclo se repete.
• O Gatilho: Eles suspeitam que isso acontece porque o ar está girando para trás tão fortemente na asa que cria um estado de "instabilidade absoluta" — uma forma sofisticada de dizer que o ar está tão turbulento que não consegue evitar oscilar por conta própria.

6. A Conclusão

Este estudo é uma história de sucesso para a modelagem computacional. Eles provaram que um programa de computador mais rápido e eficiente (ELLIPSYS) pode ser confiável para estudar a física complexa de turbinas eólicas, desde que você o verifique contra o "padrão ouro" primeiro.

Eles descobriram que, mesmo em uma pá de turbina eólica espessa, existe um "respiro" lento e rítmico do fluxo de ar que ocorre aproximadamente a cada 8 rotações da pá. Esse respiro faz com que a sustentação (a força que gira a turbina) suba e desça significamente. Compreender esse ritmo lento é crucial porque, embora possa não quebrar a turbina imediatamente, essas oscilações lentas e grandes de força podem, eventualmente, causar fadiga nos materiais ao longo de muitos anos de operação.

Em resumo: Eles construíram um túnel de vento virtual, provaram que um computador rápido poderia fazer o trabalho e descobriram que as pás das turbinas eólicas têm um "batimento cardíaco" lento e rítmico causado por bolhas de ar que se formam e estouram em sua superfície.

Resumo técnico

Resumo Técnico: Simulações Numéricas de Transição e Resposta de Longo Prazo de um Aerofólio de Turbina Eólica

Definição do Problema
A separação do escoamento em pás de turbinas eólicas, impulsionada por ângulos de incidência elevados e geometrias de aerofólios espessos, impacta significativamente o desempenho e a carga estrutural. Estudos anteriores identificaram múltiplas frequências naturais em escoamentos separados, incluindo desprendimento de esteira (wake shedding), modos de Kelvin-Helmholtz (KH) e oscilações de baixa frequência (LFOs). As LFOs são caracterizadas por fortes variações de sustentação e são frequentemente atribuídas ao estol e destol periódico do escoamento, potencialmente ligados à instabilidade absoluta dentro de bolhas de separação laminar (LSBs). Embora as LFOs tenham sido observadas em diversos estudos de aerofólios, sua ocorrência em baixos ângulos de ataque (AoA) em aerofólios de turbinas eólicas espessos, e os mecanismos específicos que as impulsionam nessas condições, requerem investigação adicional. Este estudo foca na transição e na resposta de longo prazo de uma seção de um aerofólio da turbina de referência DTU 10-MW (série FFA-W3) a um número de Reynolds de corda ($Re_c$) de $1 \times 10^5$.

Metodologia
A pesquisa utiliza simulações de grandes escalas (LES) com resolução de parede, sem modelos de subgrade-scale, para resolver as equações de Navier-Stokes incomprimíveis. Dois solvers são utilizados:

1. NEK5000: Um código de método de elemento espectral (SEM) conhecido pela dissipação numérica mínima e alta precisão.
2. ELLIPSYS: Um código de volumes finitos desenvolvido na DTU/Risø, tradicionalmente usado para RANS, mas aqui aplicado como LES implícito.

O estudo envolve uma validação sistemática do ELLIPSYS contra o NEK5000. As principais etapas metodológicas incluem:

• Sensibilidade do Domínio: Avaliação da largura do domínio de envergadura ($L_z$) comparando simulações com $L_z = 0.1c$ e $L_z = 0.2c$ para garantir a resolução das estruturas de escoamento relevantes.
• Validação do Solver: Comparação dos campos de escoamento médio, evolução de perturbações e características espectrais entre os dois solvers.
• Análise de Estabilidade: Utilização de Equações de Estabilidade Parabolizadas (PSE) e Decomposição Ortogonal Proporcional Espectral (SPOD) para analisar o crescimento de instabilidades (modos de Tollmien-Schlichting e KH) e a estrutura espacial de modos coerentes.
• Integração de Longo Prazo: Utilização do solver ELLIPSYS, computacionalmente mais eficiente, para simular a evolução do escoamento ao longo de 50 passagens de corda ($T = 50c/U_\infty$) para capturar fenômenos de baixa frequência.

Principais Resultados

• Largura do Domínio: Uma largura de envergadura de 10% da corda ($L_z = 0.1c$) foi considerada suficiente para capturar a evolução das principais perturbações e reproduzir as estatísticas de escoamento médio, coincidindo com os resultados de um domínio mais largo ($L_z = 0.2c$).
• Validação do Solver: O ELLIPSYS demonstrou concordância próxima com o NEK5000 para campos de escoamento médio e a evolução das perturbações mais amplificadas. No entanto, o ELLIPSYS subestimou a amplitude das ondas de Tollmien-Schlichting (TS) na camada limite aderente devido à maior dissipação numérica. Por outro lado, o ELLIPS_YS previu com precisão a evolução do modo KH dentro da bolha de separação laminar (LSB), alinhando-se bem com as previsões de PSE.
• Mecanismo de Transição: A transição no lado de sucção é impulsionada pela instabilidade KH que se desenvolve na camada de cisalhamento separada. O modo KH mais amplificado possui uma frequência de $f \approx 15$, e sua evolução a jusante é bem capturada por ambos os solvers e pelo PSE.
• Oscilações de Baixa Frequência (LFOs): Simulações de longo prazo revelaram uma modulação lenta do coeficiente de força normal ($C_N$) com uma amplitude de 10,5% e um período de 48 passagens de corda. Isso corresponde a uma frequência $f = 0.021$ e um número de Strouhal $St = 0.0012$.
  • Esta frequência é aproximadamente 720 vezes menor que a frequência primária da instabilidade KH.
  • No contexto da turbina DTU 10-MW, este período equivale a aproximadamente 7,7 rotações da pá.
  • O $St$ observado é inferior aos valores tipicamente relatados na literatura para outros aerofólios, ocorrendo em um ângulo de ataque efetivo relativamente pequeno de $3,1^\circ - 3,3^\circ$.

Significância e Alegações
O artigo estabelece que o ELLIPSYS é uma ferramenta viável para LES com resolução de parede de fluxos transicionais, oferecendo uma alternativa computacionalmente eficiente aos códigos de elementos espectrais para integrações de longo prazo. O estudo confirma que o estol e destol periódico do escoamento, potencialmente desencadeado pela instabilidade absoluta devido ao alto fluxo reverso (até -17% no lado de sucção), é responsável pelas LFOs observadas.

Uma descoberta significativa é a identificação de LFOs em um aerofólio espesso (24% de espessura) em um baixo ângulo de ataque, uma condição distinta dos estóis de aerofólios finos tipicamente associados a esses fenômenos na literatura. Os autores observam que o fluxo reverso em ambos os lados do aerofólio é suficiente para permitir a instabilidade absoluta, que pode impulsionar o estouro periódico da LSB em um escoamento totalmente estolado. O estudo contribui para a compreensão de como a geometria de aerofólios espessos influencia o surgimento e as características das oscilações de baixa frequência, destacando que esses fenômenos podem ocorrer em ângulos de incidência e números de Strouhal mais baixos do que o documentado anteriormente para aerofólios padrão.