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

Este estudo emprega Simulações Numéricas Diretas geometricamente resolvidas para revelar que o aumento do número de pás em turbinas eólicas de eixo vertical acelera a ruptura dos vórtices de estol dinâmico por meio de interações pá-vórtice, fazendo com que a esteira próxima transite mais rapidamente para a dinâmica de corpos rombudos e demonstrando que o número de pás, e não a razão de velocidade na ponta, é o fator primário que governa essa transição e as características de entrada a jusante.

O essencial

Imagine uma turbina eólica de eixo vertical (VAWT) não como uma máquina gigante, mas como um cata-vento girando em pé, de pé, em um rio de ar. Este artigo é como um filme em alta definição e câmera lenta que dá zoom tão de perto nas pás que podemos ver as "correntes de ar" invisíveis girando ao seu redor. Os pesquisadores usaram um supercomputador para simular exatamente como o ar se comporta logo atrás das pás giratórias, um local chamado de "esteira próxima".

Aqui está a história do que eles descobriram, explicada de forma simples:

A Visão Geral: Por que olhar para a "Esteira Próxima"?

A maioria das turbinas eólicas gira como uma hélice de avião (eixo horizontal). Mas estas de eixo vertical giram como um centrifugador de salada. Os pesquisadores queriam saber: O que acontece com o ar imediatamente após passar pelas pás?

Eles descobriram que o ar não flui suavemente atrás da turbina. Em vez disso, torna-se caótico. As pás cortam o ar, criando redemoinhos gigantes e giratórios chamados Vórtices de Estol Dinâmico (DSVs). Pense neles como redemoinhos gigantes e invisíveis que são levantados pelas pás e arrastam-se atrás da turbina.

A Principal Descoberta: O "Número de Pás" Importa Mais que a Velocidade

A equipe testou turbinas com uma pá, duas pás e três pás. Eles também alteraram a velocidade de rotação das turbinas.

Aqui está a reviravolta surpreendente: Quantas pás a turbina tem importa muito mais do que a velocidade com que ela gira.

• A Turbina de Uma Pá (O Solista):
  Imagine um dançarino solitário girando em uma sala. Quando ele gira, ele cria um redemoinho gigante e poderoso de ar atrás dele. Este redemoinho gigante (o DSV) permanece forte e intacto por muito tempo. É como uma nuvem pesada e de movimento lento que leva muito tempo para se dissipar. Como essa nuvem gigante permanece, o ar atrás da turbina permanece "desordenado" e caótico por uma longa distância.

• A Turbina de Três Pás (O Trio):
  Agora imagine três dançarinos girando próximos uns dos outros. Enquanto um dançarino gira, o ar que ele acabou de perturbar (o redemoinho) é imediatamente atingido pelo próximo dançarino.
  Os pesquisadores descobriram um novo mecanismo que chamam de "Impacto de Vórtice".

  • A Analogia: Pense no redemoinho gigante como uma bolha de sabão. Na turbina de uma pá, a bolha flutua inteira para longe. Na turbina de três pás, a próxima pá age como um alfinete, estourando a bolha prematuramente.
  • O Resultado: O redemoinho gigante e desordenado é esmagado em pequenas bolhas inofensivas (vórtices menores) antes que possa viajar longe. O ar atrás da turbina de três pás torna-se "calmo" e ordenado muito mais rápido do que a versão de uma pá.

O Efeito "Engarrafamento"

O artigo também explica que ter mais pás cria um pouco de um "engarrafamento" para o vento.

• Com mais pás, há menos espaço aberto para o vento fluir através do centro da turbina.
• Isso força o vento a passar ao redor da turbina, como a água fluindo ao redor de uma pedra em um riacho.
• Isso altera o comportamento da esteira. Em vez de ser dominada por redemoinhos giratórios gigantes (que ocorrem com baixos números de pás), a esteira começa a se comportar como a esteira atrás de uma pedra sólida simples (um "corpo rombudo").
• Por que isso é bom: A esteira "semelhante a uma pedra" se recupera (torna-se ar liso novamente) muito mais rápido do que a esteira "redemoinho".

O Teste de "Auto-Similaridade"

Os pesquisadores queriam saber exatamente quando o ar desordenado atrás da turbina se transforma novamente em ar liso e previsível. Eles usaram um truque matemático chamado "análise de auto-similaridade".

• A Analogia: Imagine tentar prever a forma de uma pluma de fumaça. No início, a fumaça é uma bagunça caótica de redemoinhos. Mas, eventualmente, ela se estabiliza em uma forma previsível, de curva em sino.
• Eles descobriram que as turbinas de três pás se estabilizam nessa forma previsível muito mais cedo do que as turbinas de uma pá. A "desordem" desaparece mais rápido.

O Que Isso Significa para o Futuro (De Acordo com o Artigo)

O artigo menciona especificamente que essas descobertas são importantes para parques eólicos onde as turbinas são colocadas muito próximas umas das outras.

• Se você colocar uma segunda turbina logo atrás da primeira, o ar que ela "respira" depende fortemente de quantas pás a primeira turbina tem.
• Se a primeira turbina tiver poucas pás, a segunda turbina será atingida por redemoinhos gigantes e caóticos.
• Se a primeira turbina tiver muitas pás, o ar já terá "curado" e se tornado mais suave quando atingir a segunda turbina.

Em resumo: Este artigo usou uma simulação por computador superpoderosa para mostrar que adicionar mais pás a uma turbina eólica de eixo vertical age como um "quebrador de redemoinhos". Ele esmaga os vórtices gigantes e desordenados de ar em pedaços minúsculos, permitindo que o vento recupere seu fluxo suave muito mais rápido. Isso é crucial para o projeto de parques eólicos onde as turbinas estão agrupadas muito próximas.

Resumo técnico

Resumo Técnico: Simulação Numérica Direta da Dinâmica da Esteira Próxima de Turbinas Eólicas de Eixo Vertical

Declaração do Problema
Enquanto as Turbinas Eólicas de Eixo Horizontal (HAWT) dominam a pesquisa atual, as Turbinas Eólicas de Eixo Vertical (VAWT) oferecem vantagens distintas para aplicações offshore, urbanas e remotas, incluindo centros de gravidade mais baixos e insensibilidade direcional. Um benefício teórico chave das VAWTs é o seu potencial para maiores densidades de potência em arranjos de espaçamento reduzido devido a regiões de esteira mais curtas e recuperação de esteira aprimorada. No entanto, a esteira próxima de uma VAWT é dominada por características de fluxo complexas associadas ao estol dinâmico, especificamente a formação e o desprendimento de Grandes Vórtices de Estol Dinâmico (DSVs). Diferentemente da esteira distante, que se comporta como uma esteira de corpo rombudo, a esteira próxima é altamente dependente de parâmetros geométricos e operacionais, como a razão de velocidade na ponta ($\lambda$) e a solidez estática ($\sigma$). Dados experimentais existentes frequentemente carecem de resolução próxima à geometria ou em três dimensões, enquanto Simulações de Grandes Vórtices (LES) de alta fidelidade acopladas ao Método da Linha de Atuador (ALM) falham em resolver geometricamente a camada limite da lâmina. Consequentemente, o ALM não consegue capturar a dinâmica da bolha de separação laminar (LSB), a formação precisa do DSV ou as interações críticas lâmina-vórtice (BVI) que modulam essas características.

Metodologia
Para abordar essas limitações, os autores apresentam Simulações Numéricas Diretas (DNS) geometricamente resolvidas de VAWTs Darrieus de lâminas retas utilizando o framework do método espectral/hp, Nektar++. O estudo investiga três configurações de turbina com uma, duas e três lâminas (perfis aerodinâmicos NACA 0018) em uma faixa de razões de velocidade na ponta ($\lambda = 1.5, 2.5, 4.5$). Este espaço de parâmetros corresponde a solididades estáticas ($\sigma$) de 0,2, 0,4 e 0,6, e uma solidez dinâmica ($\sigma_D$) variando de 0,47 a 0,94.

As simulações empregam uma formulação de Quadro de Referência Móvel (MRF) para lidar com a geometria rotativa sem remalhamento dinâmico. O fluxo é modelado como equações de Navier-Stokes incompressíveis em um número de Reynolds baseado na corda de $Re_c = 10^4$, um valor escolhido para representar ambientes urbanos/extraterrestres e garantir a viabilidade numérica, mantendo a relevância para a física fundamental. A resolução da malha satisfaz os requisitos de DNS ($\Delta y^+ \lesssim 1$) com 32 planos de Fourier na direção da envergadura. Viscosidade Espectral de Extinção (SVV) é aplicada para estabilidade numérica em regiões sub-resolvidas. A metodologia é validada contra dados experimentais de Battisti et al., mostrando forte concordância nos perfis de velocidade média no tempo na direção do fluxo em locais de esteira próxima.

Principais Contribuições e Resultados

1. Resolução do Estol Dinâmico e BVI:
   A DNS captura com sucesso o ciclo completo de estol dinâmico, incluindo o enrolamento da camada limite, formação do DSV, separação e convecção. Uma descoberta nova é a identificação de um mecanismo de impacto de vórtice impulsionado por BVI.

   • Configuração de Uma Lâmina: O DSV forma-se, separa-se e convecta-se para a esteira próxima como uma estrutura coerente de grande escala, causando aceleração significativa do fluxo no lado de sotavento.
   • Configuração de Duas Lâminas: O aumento da solidez reduz a convecção do DSV, fazendo com que ele difunda e se quebre em estruturas menores antes de atingir a esteira próxima.
   • Configuração de Três Lâminas: O espaçamento reduzido entre lâminas introduz BVI no início da passagem a montante. Essas interações ocorrem durante a fase inicial do estol dinâmico (desenvolvimento da bolha de separação laminar), perturbando a LSB e reduzindo o tamanho eventual do DSV. Além disso, à medida que a lâmina se aproxima da zona a jusante, o impacto direto de vórtices de esteiras a montante causa a quebra prematura do DSV em múltiplos núcleos distintos. Este mecanismo acelera a transição da dinâmica de vórtice coerente para a dinâmica de corpo rombudo.

2. Recuperação da Esteira e Transição para Dinâmica de Corpo Rombudo:
   O estudo quantifica a evolução a jusante da esteira. Maior solidez (mais lâminas) resulta em déficits de velocidade iniciais mais severos, mas inicia a recuperação mais cedo e mais rapidamente.

   • Baixa Solidez ($N_{lâminas}=1$): A recuperação é governada pela mistura transversal induzida por grandes DSVs, levando a uma recuperação mais lenta e assimetria persistente da esteira.
   • Alta Solidez ($N_{lâminas}=3$): A quebra de DSVs impulsionada por BVI remove as estruturas coerentes que distinguem a esteira da VAWT de um corpo rombudo. Consequentemente, a esteira transita mais rapidamente para mecanismos de recuperação associados à camada de cisalhamento, que são inerentemente mais rápidos.

3. Análise de Auto-Similaridade:
   Estendendo a análise de auto-similaridade para a esteira próxima (anteriormente limitada à esteira distante em estudos com ALM), os autores quantificam a taxa na qual a esteira perde a dependência de estruturas coerentes geradas pelas lâminas e colapsa em uma solução Gaussiana auto-similar.

   • A análise revela que o número de lâminas é mais influente do que a razão de velocidade na ponta na definição da taxa dessa transição.
   • Maiores contagens de lâminas levam a um início mais rápido da auto-similaridade. Curiosamente, a taxa de transição atinge o pico em $\lambda=2.5$ em vez do $\lambda$ mais alto, sugerindo um efeito de platô além de certa solidez dinâmica.

Significância e Alegações
O artigo alega que a abordagem de DNS geometricamente resolvida é essencial para identificar os mecanismos específicos pelos quais o número de lâminas influencia a aerodinâmica das VAWTs. A significância primária reside na descoberta da quebra prematura de DSV induzida por BVI, um mecanismo que acelera a transição para a dinâmica de esteira de corpo rombudo independentemente de efeitos de bloqueio.

Os autores afirmam que essas descobertas têm implicações diretas para o projeto de arranjos de turbinas de espaçamento reduzido e configurações de pares acoplados. Como o fluxo de entrada experimentado por um rotor a jusante depende fundamentalmente do número de lâminas da turbina a montante, pares com baixo número de lâminas encontrarão grandes DSVs coerentes, enquanto pares com alto número de lâminas encontrarão um campo mais difuso e perturbado. O artigo conclui que DNS de configurações em pares é justificada para compreender plenamente como essas condições de entrada distintas afetam o desempenho a jusante, notando que a recuperação mais rápida da esteira e o início mais precoce da auto-similaridade em solididades mais altas podem permitir uma modelagem analítica mais precisa de tais arranjos em distâncias a jusante mais curtas.