Machine Translated by Google

energias
Artigo

Projeto, fabricação e teste de desempenho de um 100-W

Turbina Eólica de Eixo Vertical de Lâmina Helicoidal em Baixo

Taxa de velocidade de ponta

Dowon Han 1 1 2
, Young Gun Heo 1,2, Nak Joon Choi 3 , Sang Hyun Nam e , Kyoung Ho Choi
Kyung Chun Kim 1,* ID
1 Escola de Engenharia Mecânica, Universidade Nacional de Pusan, Busan 46241, Coréia;
[email protected] (DH); [email protected] (YGH); [email protected] (SHN)
2 DNDE Co., Busan 48059, Coréia; [email protected] ICT Co.,
3
Busan 48059, Coréia; [email protected] * Correspondência:
[email protected]; Tel.: +82-51-510-2324

Recebido: 6 de maio de 2018; Aceito: 7 de junho de 2018; Publicado: 11 de junho de 2018

Resumo: Uma turbina eólica de pás helicoidais de eixo vertical de 100 W foi projetada, fabricada e testada em um
túnel de vento. Uma relação de velocidade de ponta relativamente baixa de 1,1 foi direcionada para uso em um
ambiente urbano com uma velocidade de vento nominal de 9 m/s e uma velocidade de rotação de 170 rpm. As
dimensões básicas foram determinadas através de um método de projeto baseado no momento de acordo com o protocolo IEC 61400-2.
A potência foi estimada por um modelo matemático que leva em consideração o desempenho aerodinâmico da forma
da pá NACA0018. A sustentação e o arrasto da pá em relação ao ângulo de ataque durante a rotação foram calculados
usando simulação de dinâmica de fluidos computacional (CFD) 2D para levar em consideração a região de estol. A
potência média calculada pelo modelo foi de 108,34 W, o que atende à meta de saída de 100 W. A turbina eólica
fabricada foi testada em um grande túnel de vento de circuito fechado, e as saídas de potência foram medidas para
determinadas velocidades do vento. Na condição de projeto, a potência medida foi de 114,7 W, 5,9% maior que a do
modelo matemático. Este resultado valida o método de dimensionamento proposto e a estimativa de potência pelo
modelo matemático.

Palavras-chave: turbina eólica de eixo vertical; lâmina helicoidal; baixa relação ponta-velocidade; projeto e fabricação; teste de
túnel de vento

1. Introdução

O consumo de combustíveis fósseis aumentou, resultando em altas emissões de CO2 e sérias mudanças
climáticas. A pesquisa sobre energia renovável está ativamente em andamento para resolver esses problemas
ambientais e em antecipação ao esgotamento dos combustíveis fósseis [1]. A energia eólica é uma fonte de energia
renovável ambientalmente amigável que não causa poluição ambiental, e seu uso está se espalhando rapidamente em
todo o mundo [2].
Desde o Protocolo de Kyoto em 1997 até a Convenção do Clima de Paris em 2015, mais países estão
se preocupando com as mudanças climáticas devido à industrialização. Para resolver este problema, metas
específicas foram definidas e a geração de energia eólica é considerada um sistema de geração de energia ideal.
A energia eólica suprime as emissões de CO2 em 828 g/kWh em comparação com a geração de energia a carvão. Na
geração de energia eólica, a emissão de poluentes como dióxido de carbono e metano é de 1/50 a 1/100 das outras
fontes de energia, assumindo que a velocidade do vento seja superior a 8 m/s [3].
A pesquisa sobre geração de energia eólica tem, portanto, sido ativamente perseguida. No início, a pesquisa
em geradores de turbinas eólicas horizontais de médio e grande porte foi o foco principal. No entanto, devido a
fatores que afetam o meio ambiente, como o ruído, tais aerogeradores são difíceis de instalar perto de residências e

Energias 2018, 11, 1517; doi:10.3390/en11061517 www.mdpi.com/journal/energies

Machine Translated by Google

Energias 2018, 11, 1517 2 de 17

têm efeitos negativos no ecossistema, como o movimento de pássaros [4]. Em contraste, pequenas turbinas eólicas
podem ser instaladas perto de áreas residenciais adjacentes a cargas de energia, e pesquisas relevantes foram
realizadas [5-7].
Um gerador de turbina eólica pode ter um eixo de rotação vertical ou horizontal. Um gerador de energia eólica
de eixo vertical é vantajoso para instalação em centros de cidades porque não é afetado pela direção do vento tanto
quanto um gerador de energia eólica de eixo horizontal. É de fácil manutenção, pois não necessita de estrutura
complicada, como dispositivos de guinada [8,9]. No caso de aerogeradores de eixo horizontal , o ângulo de ataque
devido à rotação do aerogerador é constante. Muitos estudos foram conduzidos sobre a previsão das características
aerodinâmicas das pás, e muitas tecnologias proprietárias foram estabelecidas [10-13]. No entanto, no caso de
turbinas eólicas de eixo vertical, o ângulo de ataque devido à rotação da turbina eólica muda continuamente.

Portanto, é essencial desenvolver um processo de verificação de saída para uma turbina eólica de eixo vertical.
O método mais óbvio de verificação é a experimentação. No entanto, devido a limitações espaciais, este método é
limitado a turbinas eólicas muito pequenas. Modelos matemáticos têm sido estudados para várias previsões de saída
usadas no projeto de turbinas eólicas de eixo horizontal, mas eles não levam em consideração o ângulo de inclinação
das pás verticais ou a localização do suporte ou apenas estudos parciais foram feitos.

Darrieus propôs o conceito de uma turbina eólica de eixo vertical em 1931 [14]. A primeira abordagem
simplificada é o modelo numérico de tubo de fluxo único proposto por Templin [15]. As características do aerofólio
foram calculadas usando a teoria dos elementos de pá [15]. A saída de todo o rotor é a mesma que o desempenho
de uma única pá com o comprimento da corda de todas as pás do rotor. Essa abordagem nos permite prever o
desempenho do rotor em termos de torque médio por revolução do rotor.
Strickland estendeu vários modelos de tubo de fluxo único para refletir pequenos tubos de fluxo que preservam a teoria do
momento e do elemento de pá [16]. O desempenho foi determinado considerando o número de Reynolds local em vez da
média, conforme proposto por Templin. Paraschivoiu propôs um modelo duplo de tubo de fluxo múltiplo (DMST) [17]. O
modelo é dividido na região a montante do rotor entre 0 e 180ÿ e na região a jusante após 180ÿ para cada tubo de fluxo. O
DMST pode prever razoavelmente o desempenho de uma turbina eólica de eixo vertical com uma pequena solidez e um
pequeno comprimento de corda em comparação com o raio do rotor. Islã e outros. comparou e analisou três modelos
mecânicos para projetar uma turbina eólica Darrieus com pás retas e prever o desempenho. O modelo DMST, o modelo de
vórtice livre e o modelo em cascata foram comparados [18].

Sutherland e outros. propôs um modelo de tubo de corrente e um modelo de vórtice que pode analisar a
resposta aerodinâmica de uma turbina eólica de eixo vertical usando um modelo matemático [19]. Wang et ai.
aplicaram um modelo de painel de vórtice 2D (VPM2D) a uma turbina eólica de eixo vertical de pá reta, mas este
modelo não conseguiu refletir a estol dinâmico [20]. Brusca et ai. analisaram a relação entre a relação de aspecto de
uma turbina eólica de eixo vertical com pás retas usando um código de cálculo baseado em um modelo de tubo de
fluxo múltiplo [21]. Eles concluíram que uma razão de aspecto baixa tem algumas vantagens sobre uma razão de
aspecto alta e enfatizaram que o fator de potência foi afetado pela solidez e pelo número de Reynolds do rotor.
Testes de campo ou testes de túnel de vento também foram realizados para verificar o desempenho das turbinas eólicas.
Sheldahl conduziu um teste de campo em uma turbina eólica de eixo vertical do tipo Darrieus com um aerofólio NACA 0012
e comparou com os resultados de um teste de túnel de vento ideal [22]. Bedon et ai. relataram resultados de testes de campo
para uma pá helicoidal de 1 m de diâmetro com um aerofólio NACA 0018 [23]. Recentemente, Cheng et al. realizou uma
simulação de campo de fluxo 2D de um VAWT helicoidal com quatro pás por meio de uma grande simulação de turbilhão (LES).
Eles mostraram que a variação do ângulo de ataque (AOA) e a interação lâmina-vigília sob diferentes condições de razão de
velocidade de ponta são as duas principais razões para a saída de potência do VAWT helicoidal [24].
Neste estudo, uma turbina eólica foi projetada com base em uma turbina eólica de eixo vertical do tipo elevador.
A saída inicial do projeto é de 100 W e a taxa de velocidade de ponta alvo é 1,1, que é menor do que a taxa de 4-6
de turbinas eólicas convencionais de eixo vertical. Para a relação de velocidade de ponta convencional, os ângulos
de ataque máximo e mínimo são reduzidos. Portanto, o ângulo de ataque não atinge a região de estol, de modo que
os coeficientes de sustentação e arrasto podem ser facilmente previstos e o modelo DMST convencional
Machine Translated by Google

Energias 2018, 11, 1517 3 de 17

pode ser usado. No entanto, nosso VAWT de destino foi projetado para uso urbano e possui uma relação de
velocidade de ponta muito baixa para reduzir o ruído. Na baixa taxa de velocidade de ponta, o ângulo de3 de 17
ataque
está sujeito à região de estol. Para levar em consideração a região de estol, os coeficientes de sustentação e
arrasto no ângulo de ataque correspondente foram medidos para reduzir o ruído. Na baixa taxa de
velocidade de ponta, o ângulo de ataque está sujeito à região de estol. Para calculado usando o método 2D CFD.
estol, o
usando
região Como
experimentais.
ede
coeficiente
pode
o método
estol.
resultado,
serum
aplicado
CFD
de
bom
Sem
sustentação
a
2D.
acordo
saída
para
qualquer
Como de
projetar
com
resultado,
potência
eabordagem
arrasto
os VAWT
dados
prevista
anosaída
em
experimentais.
correspondente
de uma
alta
foi
deconfirmada
fidelidade,
relação
potência
Sem
de
ângulo
prevista
como
velocidade
com
qualquer
de
uma
projeto
foiataque
confirmada
boa
abordagem
deCFD
consideração
ponta
de 3D,
acordo
baixa,
com
de
foram
alta
com
nosso
incluindo
da
calculados
fidelidade,
os
região
modelo
dados
a decomo
DMST
lâminas
de helicoidais
velocidade
projeto de 3Dhorizontal
deeixo
ponta CFD,
baixa,foram selecionadas
incluindo
, nosso modelo para
a regiãoDMST
de baixo
estol.
pode ruído
turbina e baixaépara
ser aplicado
eólica, flutuação
difícil de VAWT
projetar saída.
projetar Ao
em contrário
o tamanho uma umadeturbina
de relação um
eólica de eixo vertical com uma saída alvo As pás helicoidais foram selecionadas para baixo ruído e baixa
flutuação de saída. Ao contrário de um eixo horizontal por causa da mudança constante do ângulo de ataque
uma do devido
aerofólio
saída
com alvo Aà fórmula
rotação de
do projeto
rotor. turbina
não reflete
eólica,
asécaracterísticas
difícil projetar o
aerodinâmicas
tamanho de uma
do aerofólio
turbina eólica
da pá,de
tão
eixo
separadas
vertical
verificação por causaéda
de saída mudança constante do ângulo de ataque do aerofólio devido à rotação do rotor. A
necessária.

A saída do projeto de pás de turbinas eólicas de eixo vertical foi analisada, a fórmula de projeto não
projetoUm
reflete as
decaracterísticas
rotor
pás de
realturbinas
foi fabricado
aerodinâmicas
eólicase de
umeixo
teste
dovertical
aerofólio
de túnel
foide
da
analisada
vento
pá, portanto,
foipor
realizado
meio
saída
dapara
separada
obtenção
verificação.
por
da curva
umAmodelo
saída
de saída
do
matemático.
em
diferentes velocidades de vento. A potência de saída na velocidade do vento projetada e o modelo
matemático. Um rotor real foi fabricado e um teste de túnel de vento foi realizado para obter a velocidade de
rotação do projeto confirmada por experimentos. a curva de saída em diferentes velocidades de vento.
A potência de saída na velocidade do vento projetada e na velocidade de rotação projetada foi confirmada por
experimentos.

2. Projeto de turbina eólica de eixo vertical de pá helicoidal de 100

W 2. Projeto de turbina eólica de eixo vertical de pá helicoidal
uma
de
necessário
alta,
como
inclinação.
relativamente
100
o
desvantagem
o
dispositivo
WInNo
No
um
Embora
caso
caso
alta,
grande
deem
de
ode
aguinada
formato
eficiência
uma
número
que
umaa turbina
turbina
direção
dade
de
lâmina
dispositivos,
eólica
geração
eólica
do e
évento
de
um
complicado.
deseja
eixo
eixo
dispositivo
écomo
limitada.
relativamente
horizontal,
horizontal,
um
Háde
dispositivo
No
também
inclinação.
um
caso
é grande
dade
forma
Embora
guinada
número
é complicado.
ae
deeficiência
um
dispositivos
dispositivo
Hádetambém
geração
é de
necessário,
uma
seja
a estrutura
direção desvantagem
do simples,
vento
é simples, éem
[8].ecentro que apara
ée vantajoso
vantajoso
porque direção
não
para do vento
instalação
háinstalação
restrição écentro
noem
da limitada.
direção
umada No
cidade caso
cidade
do vento de
[8].umhá
umporque aerogerador
aerogerador
nenhuma de eixo vertical,
de eixorestrição
vertical,da
a estrutura é

Os tipos típicos de pás de turbinas eólicas de eixo vertical são Darrieus, gyro-mill, Savonius e helicoidais Os
tipos típicos de pás de turbinas eólicas de eixo vertical são Darrieus, gyro-mill, Savonius e pás
helicoidais. O tipo helicoidal é vantajoso porque a faixa de flutuação da saída é menor do que as lâminas. O tipo
Darrieus
menos helicoidal
Darrieus
menos
ruído
ou giro-mill,
do
ruído
ou
éque
vantajoso
gyro-mill,
doum
eque
orotor
desempenho
umeporque
Savonius,
orotor
desempenho
Savonius,
a faixa
deque
partida
deéde
que
flutuação
umpartida
automática
rotor
é umdo
da
automática
rotor
tipo
saída
édo
arrasto.
melhor.
tipo
é menor
éarrasto.
melhor.
Também
também
do que
É menos
tem
possui
as menos
páscarga
asconvencionais
lâminas
carga
mecânica
mecânica
convencionais
e e

A Figura 1 mostra um esquema básico para o projeto aerodinâmico de uma turbina eólica de eixo vertical. A Figura 1 mostra um esquema básico para o projeto aerodinâmico de uma turbina eólica de eixo vertical.

A equação
eólica (1)
[25].
A equação representa
(1) representaaapotência dovento
potência do ventoqueque flui para
flui para o da
o rotor rotor da turbina
turbina eólica [25].

Figura 1. Parâmetros básicos de projeto do aerogerador vertical.

Figura 1. Parâmetros básicos de projeto do aerogerador vertical.
Machine Translated by Google

Energias 2018, 11, 1517 4 de 17

1 3
P = ÿAÿV 2 (1)

A equação (2) é a potência mecânica gerada pela rotação do rotor da turbina eólica.

P = Tÿ (2)

A proporção da potência convertida da potência do rotor para o fluxo de energia eólica é chamada de coeficiente de
potência, que é um conceito de eficiência de conversão de energia aerodinâmica. Teoricamente, o valor máximo do coeficiente de
potência é 0,593 em uma turbina eólica de eixo horizontal, conhecido como limite de Betz. O limite de Betz é derivado da teoria
do momento do disco do atuador e é o máximo teórico assumindo que o fluxo é estacionário, invíscido e irrotacional [26]. A turbina
Darrieus é uma típica turbina eólica de eixo vertical do tipo elevador e tem um coeficiente de potência máxima de cerca de 0,4 a
uma taxa de velocidade de ponta de 5 [9].

A equação (3) mostra a potência (P) do aerogerador considerando o coeficiente de potência

(CP) e eficiência de transmissão de potência (ÿ). Aÿ é a área de varredura do rotor (ver Equação (5)).

1 3
P = ÿAÿV 2 ÿCP (3)

A relação de velocidade de ponta (ÿ) está intimamente relacionada com o coeficiente de potência. A razão da velocidade da
ponta é definida como a razão entre a velocidade da ponta da pá e a velocidade do vento na qual a ponta da pá se move com rotação,
conforme mostrado na Equação (4).
Rÿÿ
= (4)
EM

Todos os rotores de turbinas eólicas têm uma relação de velocidade de ponta ideal com potência máxima. A
relação ótima está relacionada com a mudança da velocidade do vento de entrada. A área de varredura do rotor (Aÿ)
é determinada pelo raio e altura da turbina eólica.

Aÿ = 2RH (5)

A área varrida pelo vento deve considerar a altura do rotor (H) e a relação de aspecto em relação ao raio (R). Quanto maior
o raio do rotor, maior o torque gerado, mas quanto maior o comprimento da escora, menor a estabilidade estrutural. No entanto,
quando a altura do rotor é maior, o torque gerado é menor e a velocidade de rotação do rotor deve ser aumentada para obter a
mesma potência de saída. A relação de aspecto (AR) pode ser expressa como

H
COM = ÿ H = AR × 2R (6)
2R

A solidez (ÿ) é uma importante variável que determina o desempenho dos aerogeradores. A solidez é definida como a razão
entre a área projetada total (NC) da pá do rotor e a área rotacional do rotor da turbina eólica. A área projetada é a projeção na
seção vertical do eixo rotativo e pode ser expressa como

NC
ÿ= (7)
2ÿR

O comprimento da corda da pá (C) pode ser calculado usando a solidez. O comprimento da corda é o comprimento do
aerofólio e é uma importante variável de projeto porque o torque gerado muda de acordo com o comprimento da corda.

Geradores de energia eólica urbana devem operar em baixa velocidade com baixo ruído. Uma turbina eólica
Savonius pode girar com uma relação de velocidade de ponta inferior a 1, mas a vibração e o ruído são severos
devido às características de um rotor do tipo arrasto. Entre as turbinas eólicas de eixo vertical tipo elevador, pás helicoidais
Machine Translated by Google

Energias 2018, 11, 1517 5 de 17

as turbinas eólicas têm uma faixa mais estreita de flutuação de saída em comparação com as turbinas eólicas
5 de 17
Darrieus e gyro-mill, e sua eficiência é maior devido à maior área varrida pelo vento.
O objetivo deste estudo é projetar uma pá de turbina eólica de eixo vertical de baixa velocidade com uma velocidade de ponta velocidade de 9m/s.

Projetamos a turbina eólica com a taxa de velocidade definida pelo projetista de 1,1 a uma velocidade nominal do vento de 9 m/s. Projetamos a turbina

eólica com o designer Classe S após modificar a Classe 1 identificada na IEC 61400-2: “Requisitos de projeto para pequenos ventos definidos Classe S

após modificar
parâmetros a Classe
básicos 1 identificada
do pequeno sistema na
de IEC 61400-2:
turbinas “Requisitos
eólicas de2projeto
.” A Figura mostrapara
uma pequenas turbinas
imagem dos eólicas”.
parâmetros A Figura
básicos 2 mostra sistema
do pequeno uma imagem dos eólica
de turbina
da IEC IEC
velocidade 61400-2
nominal [27]. Aselecionada
do vento velocidade nominal do vento
foi de 9 m/s, que éselecionada foi de 9 m/s,
inferior à velocidade que édo
nominal menor
ventoque a velocidade ànominal
correspondente Classedo
1. vento 61400-2
A altura [27].
do cubo do Arotor é

e ode 8 m.de
perfil Considera-se cisalhamento
velocidade (V(Z)) do vento,
é mostrado correspondente
na Equação à Classe
(8). perfil (V(Z)) 1. A altura
é mostrado nado cubo do(8).
Equação rotor é de 8 m. O cisalhamento do vento é considerado

Figura 2. Força e velocidades atuando na turbina eólica vertical para vários ângulos de azimute. Figura 2. Força e velocidades atuando na turbina eólica vertical para vários ângulos de azimute.

a
V(Z) = Vre f Z/Zre f () = (8)
(/) (8)

O índice de força de cisalhamento do vento (a) pode variar de acordo com a rugosidade da superfície do terreno [28].
O índice de força de cisalhamento do vento (a) pode variar de acordo com a rugosidade da superfície do terreno
Neste estudo, 0,23 foi selecionado para a condição de uma floresta ou um pequeno número de edifícios.
[28]. Neste estudo, 0,23 foi selecionado para a condição de uma floresta ou um pequeno número de edifícios.
Amáxima
Tabela 1
doresume os parâmetros
vento, Tabela deos
1 resume projeto aplicados
parâmetros neste estudo.
de projeto A velocidade
A velocidade máxima do vento, a intensidade da turbulência e a inclinação adimensional são as mesmas da Classe 1. A saída alvo é de 100 W, a
aplicados neste estudo.
intensidade da turbulência e a inclinação adimensional são as mesmas da Classe 1. A saída alvo é de 100 W e a densidade do ar é 1,225 kg/
m3 , conforme especificado na IEC 61400-2. A eficiência da transmissão de potência e a densidade do ar é de 1,225 kg/m3 , conforme especificado na IEC
do61400-2.
projeto, assumida
A eficiência
como
de transmissão
95%, e o coeficiente
de energiade foi
potência
assumida
médio
como
da 95%
turbina
e ofoicoeficiente
definido como
de potência
0,15. A médio
saída do
da projeto,
turbina foi
a densidade
definido como
do ar,
0,15.
a velocidade
A saída do
vento projetada e a eficiência são definidas na Equação (3). densidade, velocidade do vento de projeto e eficiência são definidos na Equação (3).

Tabela 1. Parâmetros básicos do aerogerador vertical. Tabela 1. Parâmetros básicos do aerogerador vertical.

Classe SWT S
Classe SWT

[m/s] S
50
Vre f [m/s] [m/
s]
Velocidade [m/ 50

s] ÿ [-] 99
I15 [-] [-] 0,18 0,18

a 2
[-] 2
Machine Translated by Google

Energias 2018, 11, 1517 6 de 17

A equação (9) é usada para calcular o raio do rotor da turbina eólica através da relação entre a área varrida
pelo vento e a razão de aspecto. Assim, foi escolhido um raio de 0,55 m e uma altura de 1,43 m .

Aÿ
R= (9)
4 × AR

Foi escolhida uma taxa de velocidade de ponta de 1,1. No caso de uma turbina eólica do tipo Darrieus, o
coeficiente de potência máximo pode ser alcançado em taxas de velocidade de ponta entre 4 e 6. Na velocidade do
vento projetada de 9 m/s com um raio de rotor de 0,55 m, a velocidade de rotação está entre 630 e 950 rpm, o que
é inadequado para uso no centro da cidade. Uma velocidade rotacional de 170 rpm foi derivada aplicando o raio, a
velocidade do vento projetada e a razão da velocidade de ponta na Equação (4). O número de pás foi escolhido
como 4 e a solidez foi definida como 0,3, que foi substituída na Equação (7) para determinar o comprimento da corda de 0,25 m.
A Tabela 2 mostra as especificações da pá do rotor helicoidal. Embora o raio e a altura do rotor da
turbina eólica de eixo vertical possam ser derivados, a relação de saída do aerofólio da pá, o número
de pás e a forma da pá tridimensional não são refletidos na fórmula do projeto. Mais pesquisas sobre
as características de fluxo são necessárias.

Tabela 2. Especificações da lâmina.

Item Descrição

tipo de rotor Helicoidal

Saída de potência nominal 100W

Velocidade nominal do vento 9 m/s

Coeficiente de potência 0,15

Área varrida 1,57 m2

Proporção da tela 1,3

Raio do rotor 0,55 m

Altura do Rotor 1,43 m

Velocidade de rotação 170 rpm

Solidez 0,3

Comprimento do acorde 0,25 m

Número de lâminas
aerofólio 4 NACA0018

A equação de projeto para a turbina eólica de eixo vertical não reflete os fatores para o aerofólio da pá,
portanto, é necessária uma previsão adicional da potência de saída. Neste estudo, a potência aerodinâmica do rotor
do aerogerador foi investigada aplicando um aerofólio NACA 0018 e um modelo matemático usando as forças de
sustentação e arrasto do aerofólio de acordo com o ângulo de ataque. Ao contrário da pá de uma turbina eólica de
eixo horizontal, que tem um ângulo de ataque fixo, o ângulo de ataque varia para uma turbina eólica de eixo vertical,
dependendo do ângulo de rotação do rotor [29] . A Figura 2 apresenta o vetor velocidade de ponta e os vetores
sustentação e arrasto gerados pela rotação da pá da turbina.
O ângulo de ataque muda com o vetor de velocidade da ponta da pá e o vetor de velocidade do vento influxo.
A soma vetorial (W) do vetor de velocidade de ponta e do vetor de velocidade do vento entrante (V) é calculada pela Equação
(10). O valor máximo ocorre em ÿ = 0ÿ , e o valor mínimo ocorre em ÿ = 180ÿ .

2
W = V2 ÿ ÿ sen2 ÿ + cos2 ÿ = V 1 + 2ÿ cos ÿ + ÿ2 (10)

O ângulo de ataque (ÿ) é o ângulo entre a soma vetorial e a direção do comprimento da corda.
À medida que a soma vetorial muda, o ângulo de ataque tem um valor positivo na região a montante do rotor e um
valor negativo na região a jusante. O ângulo de ataque pode ser expresso pela Equação (11).

pecado eu
ÿ = tanÿ1 (11)
cos ÿ + ÿ
Machine Translated by Google
7 de 17

ÿpecado
7 de 17
Energias 2018, 11, 1517 ÿ=ÿ (11)
ÿ + ÿcos

àConforme
velocidade
ângulo
medida mostrado
deque
ataque
de aponta. ana
páégira equação,
Conforme
razão
de acordo acom
variável
de velocidade
mostrado a taxa
na
dedominante
equação,
de
ponta. que
velocidade
A aFiguraafeta
variável o ângulo
de3 ponta.
mostra
dominante deque
ataque
Aoproporção.
ângulo afetaéoa
de ataque
Ado
ângulo
for
velocidade
direção
velocidade
alcance
ocorre
Figura
ovetor
vetor
de
porque
do
do
3
deataque
de
mostra
vetor
da
de
ângulo
velocidade
velocidade
ponta.
ponta.
quanto
de
que
ode
velocidade
ângulo
Isso
Quanto
o
ataque
maior
de
aerofólio
da
ocorre
avanço
de
ponta
for
maior
que
ataque
de
oporque
recebe
vetor
o
avanço
daaaerofólio
pá,
razão
àde
medida
adurante
maior
eda
direção
velocidade
ade
direção
pá,
recebe
velocidade
será
que
ea da
orazão
aângulo
durante
da
soma
da
pá
taxa
soma
ponta
de
gira
de
de
do
de
velocidade
aponta,
de
de
velocidade
vetor
rotação.
ataque
daacordo
vetores
pá,se
menor
maior
torna-se
aproxima
rotação.
de
com
se
da
ponta,
oserá
aproxima
ponta.
a
alcance
taxa
menor.
Quanto
da
amenor
taxa
Isso
de
direção
do
da
Ode
maior
o
alcance
Arotor.
faixaododo ângulo
rotor. de ataque do aerofólio é de cerca de ±66,5° em um ataque
ângulo de ataque do aerofólio é de cerca de ±66,5ÿ em uma rotação de rotação do torna-se menor.

Figura 3. Variação do ângulo de ataque em uma revolução da pá para diferentes taxas de velocidade de ponta. Figura 3. Variação do ângulo de ataque em uma revolução da pá para diferentes taxas de velocidade de ponta.

A Tabela 3 mostra o ângulo de ataque em ângulos de rotação do rotor de 0° a 180° quando a taxa de velocidade de ponta A
Tabela 3 mostra o ângulo de ataque em ângulos de rotação de rotor de 0° a 180° quando a taxa de velocidade de ponta é
1,1. é 1.1.

Tabela 3. Ângulo de ataque na rotação da pá (ÿ = 1,1).

Tabela 3. Ângulo de ataque em uma revolução da pá (ÿ = 1,1(.

Azimute
Ângulo de azimute [ÿ ]
Ângulo [°] SIM [°] AOA [ÿ ]

00 0,000
0,000
10 4,789
10 4,789
20 9,575
20 9,575
30 14,354
4030 14,354 19,122
5040 23.875
60 19.122
50 28.607
70 23.875
60 33.311
80 28.607
9070 37.976

10080 33.311 42.589 37.976 47.127

110 51.558
90 42.589
120 55.828
100 47.127
130 59.840
110 51.558
140 63.400
120
150 55.828 66.069
160
130 66,568
170 59,840
140 59,298
180 63,400 0,000
Machine Translated by Google
8 de 17

Os coeficientes de sustentação e arrasto (CL, CD) do NACA 0018 na Tabela 3 são definidos nas Equações (12)
Energias 2018, 11, 1517 8 de 17
e (13), respectivamente [30].

ÿÿ
Os coeficientes
NACA 0018 na Tabela
de sustentação
3 são definidos
e=arrasto
ÿÿ nas
(CL,
Equações
CD) do (12) (12) 1 ÿ 2
e (13), respectivamente [30].
CL = (12)
12
ÿÿFL
ÿcHW2
= ÿÿ (13)
1 FD
ÿ2
CD = (13)
ÿcHW2
12
Nessas equações, a sustentação (FL) e o arrasto (FD) são obtidos por meio da análise CFD 2D, que foi Nessas
equações, a sustentação (FL) e o arrasto (FD) são obtidos pela análise CFD 2D, que foi realizado para investigar
o desempenho aerodinâmico das pás do VAWT. A Figura 4 mostra o realizado para investigar o desempenho
aerodinâmico das pás do VAWT. A Figura 4 mostra o domínio de análise de fluxo bidimensional e a posição da pá. O
campo de fluxo foi expandido para domínio de análise de fluxo bidimensional e a posição da lâmina. O campo de fluxo
foi expandido para cerca de 15D pela largura de 7D com base no diâmetro do rotor (D), 3D da entrada do campo de
fluxo para cerca de 15D pela largura de 7D com base no diâmetro do rotor (D), 3D do entrada do campo de fluxo para
o centro do rotor e ao redor do centro do rotor para a saída. O centro de análise de sensibilidade da grade do rotor e
ao redor do centro do rotor até a saída. A análise de sensibilidade da grade foi realizada antes da análise de fluxo
bidimensional. O sistema de grade selecionado tem 1.162.500 executados antes da análise de fluxo bidimensional. O
sistema de grade selecionado possui 1.162.500 nós, nós, 575.142 elementos e o valor máximo de Y+ é 2,86. Um
sistema de grade mais denso é construído em torno dos 575.142 elementos e o valor máximo de Y+ é 2,86. Um sistema
de grade mais denso é construído em torno do rotor do rotor , incluindo uma parte da região a montante do rotor, de
modo que as características de fluxo geradas, incluindo uma parte da região a montante do rotor, de modo que as
características de fluxo geradas ao redor da pá possam ser bem descrito. Na superfície da lâmina, uma grade de
prisma foi densamente a lâmina pode ser bem descrita. Na superfície da pá, uma grade de prisma foi densamente
arranjada para prever o gradiente de velocidade e a separação do fluxo devido ao efeito viscoso da ordem para prever
o gradiente de velocidade e a separação do fluxo devido ao efeito viscoso do fluido. fluido.

(a)

(b)

Figura 4. Cont.
Machine Translated by Google

Energias 2018, 11, 1517 9 de

17 9 de 17

(c)

(d)

Figura 4. Geometria 2D da turbina eólica vertical e sistema de malha; (a) Geometria do campo de fluxo 2D,
incluindo a Figura 4. Geometria 2D da turbina eólica vertical e sistema de malha; (a) Geometria do campo de fluxo
2D incluindo quatro pás, (b) posição das pás turbina eólica vertical 2D, (c) sistema de malha do campo de fluxo 2D,
e (d) malha de quatro pás, (b) posição das pás turbina eólica vertical 2D, ( c) sistema de malha do campo de fluxo
2D, e (d) sistema próximo à
Lâmina 1. sistema de malha próximo à Lâmina 1.

A Tabela 4 resume as condições de contorno aplicadas à análise CFD 2D. A entrada foi definida como Tabela 4 resume as condições de contorno aplicadas à análise CFD 2D. A entrada foi ajustada para 9 9

m/s, que é a velocidade do vento projetada, e as condições de saída foram condições de contorno de pressão. m/s, que é a velocidade do vento projetada, e as condições de saída eram condições de contorno de pressão.

A área ao redor do rotor é dividida em áreas separadas, e a área do rotor recebe uma condição de rotação A área ao
redor do rotor é dividida em áreas separadas, e a área do rotor recebe uma rotação e a área de descanso recebe uma
condição
uma
adesiva. Adeanálise
condição parada. foiArealizada
de parada. velocidade
porde
A velocidade rotação
análise foi de
de rotação foi170
transiente de arpm
e 170 eae
lâmina
rpm superfície
lâminada
decondição
foia submetida superfície
à condiçãoefoi
aadesiva.
área de Adescanso
submetida análise recebem
à condição
foi realizada
por transiente
entrada, análise
o intervalo
e intervalodede
tempo
tempofoifoi
dado
dado
umumtempo
tempo
correspondente
correspondente aoao
ângulo
ângulo
dede
rotação
rotaçãododo .rotor
Para
rotor dedeo fluido
1ÿ1°. Parade ar
com densidade
densidade
análise URANS de
de 1,225 1,225
foi kg/m3 kg/m3
realizada ae
a 1 como
atm 1uma
atm°C
25 e especificado
25 ÿCtransiente.
análise especificado
na IECna IECaplicado
O tempo 61400-2
61400-2 efoi aplicadodae turbulência
foi a intensidade
a intensidade oda
fluido de entrada,
turbulência ar 18%.
foi defoi18%.
de com
A
A análise URANS foi realizada como um intervalo transiente de cerca de 9,8 × 10ÿ4 s, correspondendo a 1ÿ de rotação
do rotor. O modelo de turbulência adotado análise. O intervalo de tempo foi de cerca de 9,8 × 10ÿ4 s, correspondendo a
1° de rotação do rotor. A turbulência é
o mesmo
mesmo modelo
modelo de
de turbulência
turbulência SST.
SST. modelo adotou o
Machine Translated by Google

Energias 2018, 11, 1517 10 de 17

Tabela 4. Condições de contorno para análise CFD 2D.

Item Descrição
Entrada 9m/s
Tomada Pressão atmosférica (condição de abertura)
Lado (d) Condição de deslizamento livre
Superfície da lâmina Condição antiderrapante
Velocidade de rotação 170 rpm
Figura 4. Geometria 2D da turbina eólica vertical e sistema de malha; (a) Geometria do campo de fluxo 2D incluindo
Direção de rotação Sentido anti-horário
quatro pás, (b) posição das pás turbina eólica vertical 2D, (c) sistema de malha do campo de fluxo 2D e (d) sistema de
malha perto da pá 1.

A Figura
de5do
e arrasto mostra
ataque os em
bidimensional
aerofólio coeficientes
obtido
emaestado
cada de
partirângulo
da sustentação
estacionário.
Figura
de ataque e os
5 mostra arrasto
obtido do aerofólio
coeficientes
a partir emCFD
da análise cada ângulo
de sustentação
a análise CFD bidimensional em estado estacionário.

pá Variação
de 0ÿ
Figura
(NACA0018)
5. a 180ÿ do coeficiente
Figura 5. Variação
calculado por simulação desustentação
dos coeficientes de sustentaçãoe arrastoeem
CFD 2D. calculado arrasto
um ânguloem um ângulo
de azimute
por simulação
da pá de de
0° a azimute
CFD 2D. da
180° (NACA0018)

O coeficiente normal (CN) e o coeficiente tangencial (CT) são gerados a partir da lâmina por
O coeficiente
por normal
usando o coeficiente (CN) e eoocoeficiente
de sustentação tangencial
coeficiente de arrasto (CT)
e calculado são
usando gerados
as Equações (14)aepartir
(15). da lâmina
ÿsinÿ + ÿcosÿ
usando o coeficiente de sustentação e o coeficiente = ÿ e calculado usando as Equações (14) e
de arrasto (15).
(14)

CN = CL cos ÿ + CD sen ÿ (14)

CT = CL sen ÿ ÿ CD cos ÿ (15)

A força normal (FN) e força tangencial (FT) da pá podem ser calculadas através dos
coeficientes normal e tangencial usando as Equações (16) e (17).
1
FN(ÿ) = 2 ÿcHW2CN (16)

1
FT(ÿ) = 2 ÿcHW2CT (17)

A saída de potência pode finalmente ser calculada usando o torque da lâmina (Equação (18)) e o ângulo
velocidade usando a força tangencial:

1
T(ÿ) = ÿcHW2CTR (18)
2
Machine Translated by Google

Energias 2018, 11, 1517 11 de 17

As saídas de potência instantânea e média do rotor projetado são dadas pelas Equações
(19) e (20).
P(ÿ) = T(ÿ) × ÿ (19)
N pi

Pavimentar = p P(ÿ)dÿ (20)

2 0

3. Teste de túnel de vento

Um rotor de turbina eólica de eixo vertical helicoidal de 100 W foi fabricado com base nas dimensões do
projeto na Tabela 2. Os rotores, cubos e suportes foram projetados e são estruturalmente estáveis de acordo
com a IEC 61400-2. O eixo de rotação conecta os cubos superior e inferior e é projetado para suportar a flexão
causada pelo vento. Um tubo de aço carbono (50A × Sch. 40) com um diâmetro externo de 60,5 mm foi usado
para tubulação de pressão.
O suporte é uma parte importante para conectar uma pá ao cubo. É um dos componentes que recebe a maior carga. O
suporte suporta o peso do rotor e a força centrífuga da rotação das pás. O projeto foi feito considerando a posição da pá e a
posição do cubo. O cubo é uma peça que fixa o eixo rotativo e o suporte. Quando o rotor é girado, o suporte é projetado para
não se mover na direção de rotação. Por fim, foram fabricadas lâminas helicoidais com material FRP, que possui excelente
conformabilidade. Todas as peças foram projetadas para fabricação. A Figura 6 mostra a turbina fabricada.

12 de 17

Figura 6. Turbina eólica de eixo vertical de pás helicoidais de 100 W.

Figura 6. Turbina eólica de eixo vertical de pás helicoidais de 100 W.
Machine Translated by Google

Energias 2018, 11, 1517 12 de 17

O teste de túnel de vento foi realizado no grande centro de túnel de vento de circuito fechado na
Universidade Nacional de Chonbuk, na Coréia. As dimensões da seção de teste são 5 m (L) × 2,5 m (A) × 20 m
(L), e a velocidade do vento pode ser controlada até 30 m/s. Uma estrutura de favo de mel e uma malha quadrada
são instaladas dentro do túnel de modo que a intensidade de turbulência do fluido de entrada pode ser inferior a
1,5% e a distribuição média da velocidade de fluxo pode ser inferior a ± 2%. A Figura 7 mostra uma foto da seção
de teste com a turbina eólicaFigura
fabricada instalada.
6. Turbina eólica de eixo vertical de pás helicoidais de 100 W.

Figura 7. Seção de teste do túnel de vento.

túnel
Figurade7.vento.
Seção de teste do

A velocidade do vento é controlada através de um painel de controle do ventilador e confirmada através

de um anemômetro instalado no teto da seção de teste do túnel de vento. O teste foi realizado após a
estabilização da velocidade do vento de acordo com o anemômetro. A energia elétrica do gerador de turbina
eólica foi medida por um dissipador de energia. A tensão e a corrente de eletricidade do gerador foram medidas
enquanto o dissipador de energia absorve a energia gerada. Além disso, a frequência trifásica, tensão e corrente
do gerador foram medidas usando um osciloscópio. A velocidade de rotação do aerogerador pode ser estimada
através da frequência trifásica (fe) do gerador,

NmNP 120 fe fe = ÿ Nm
= 120 (21)
POR EXEMPLO

onde Np é o número de pólos do gerador.

O teste do túnel de vento é dividido na medição inicial da velocidade do vento e na medição da
geração de energia. A velocidade inicial do vento é a velocidade do vento na qual a turbina eólica de eixo
vertical está se movendo desde a paralisação e é medida sem resistência elétrica. A velocidade inicial do
vento é medida com um aumento gradual da velocidade do vento com intervalos de 0,5 m/s e mantida por
5 min para estabilização em cada velocidade do vento.
A cada velocidade do vento, foram medidas a tensão, a corrente e a frequência trifásica gerada pelo
gerador do aerogerador. Usando um multímetro e um osciloscópio, as curvas de tensão e corrente ao longo
do tempo foram todas medidas e a saída calculada. A frequência medida foi usada para obter a velocidade
de rotação da turbina eólica da Equação (21). O teste foi realizado com aumentos na velocidade do vento
de 1 m/s a partir da velocidade inicial do vento. Em cada seção de velocidade do vento, o teste foi realizado
alterando a taxa de serviço do controlador. A faixa de saída, incluindo o ponto de potência máxima em uma
determinada velocidade do vento, foi examinada. Quando todas as condições foram alteradas, os resultados
foram registrados após a obtenção de um estado estabilizado. A Tabela 5 fornece uma breve descrição do
condições de teste.
Machine Translated by Google

Energias 2018, 11, 1517 13 de 17

Tabela 5. Cenário para ensaio em túnel de vento.

Condição de teste Velocidade do vento [m/s] Condição do medidor de controle inteligente

Iniciando o teste de velocidade do vento 0~ Taxa de serviço = 0

5
6
7

máx. teste de saída de energia 8 Variação da taxa de serviço

9
10
11

4. Resultados e discussão
Durante a rotação das pás, o torque gerado foi obtido pela Equação (18), conforme Figura 8. O torque máximo
14 de 17
ocorreu entre ângulos de rotação de 50 e 70ÿ . Após ângulo de 70°, condições de operação acima de 260 rpm. O VAWT
, o ângulo
construído de ataque
foi projetado paraatinge a região
ter uma dede
potência estol seguido
saída de 100por
Wa uma
umarápida diminuição
velocidade do de
do vento torque.
9 m/sAcima
a uma de 90° de
velocidade
ângulo,
do o 170
rotor de que rpm
significa
(TSRfluxo reverso,
= 1,1). Como oo experimento
valor do torque é quase para
foi realizado zero validação
e a lâminadanão cria nenhuma
potência força
no ponto de de levantamento.
projeto,
O valorgráficos
geramos de saídacomfinal calculado
potências a partir
até 260 rpm. das Equações (19) e (20) é 108,34 W e o coeficiente de potência
é 0,154, o que mostra uma discrepância de 8,34% do valor de projeto, mas é maior que a saída desejada.

Figura 8. Variação do torque de uma pá por revolução. Figura 8. Variação do torque de uma pá por revolução.

A Figura 9 mostra as curvas de potência derivadas de valores experimentais em cada velocidade do vento.
Todos os rotores de turbinas eólicas têm pontos de saída máxima em cada velocidade do vento, e o método de
rastreamento do ponto de potência máxima (MPPT) é usado para controlar a potência máxima do vento em cada velocidade do vento.
Este método minimiza a perda de desempenho da turbina eólica [31]. O ponto máximo de saída aumenta acentuadamente
à medida que a velocidade do vento aumenta porque a saída aumenta proporcionalmente ao cubo da velocidade do vento,
conforme mostrado na Equação (1). Devido à estabilidade estrutural, o teste não foi realizado em condições de operação
acima de 260 rpm. O VAWT construído foi projetado para ter uma potência de saída de 100 W a uma velocidade do vento
de 9 m/s a uma velocidade do rotor de 170 rpm (TSR = 1,1). Desde que o experimento foi
Machine Translated by Google

Energias 2018, 11, 1517 14 de 17

realizada para a validação da potência no ponto de projeto, geramos gráficos com a Figura 8.
saídas de potência até 260 Variação
rpm. do torque de uma pá por revolução.

Figura 9.
resultados
de
obtida
ventono Gráfico
teste dede
(símbolos:
do teste; potência
túnel
Figura
linhas: dede
decurva
9.
vento acordo
Gráfico comresultados
(símbolos:
ajuste). a velocidade
de potência de acordo do vento
do teste;
com obtida
linhas: no de
a velocidade
curva teste
do de túnel
ajuste).
vento

A Tabela 6 resume a potência máxima e as revoluções por minuto em cada velocidade do vento. A
potência de saída ocorre quando a velocidade de rotação é controlada para manter as revoluções nominais
de 170 rpm. O controle de saída e o controle de velocidade são funções que devem ser aplicadas de
acordo com IEC 61400-12-1: “Medições de desempenho de energia de turbinas eólicas produtoras de
eletricidade”. Eles são aplicados de várias maneiras para refletir a operação segura de geradores de
turbinas eólicas [32-35].

Tabela 6. Potência de saída e coeficiente de potência de acordo com a velocidade do vento obtida no teste de túnel de vento.

Velocidade do Velocidade de Potência da saída Observação

Cp
Vento [m/s] rotação [rpm] [EM]

~4 25.6 0,05 0,001 velocidade inicial 3,5 m/s

5 111 10.44 0,087 -
6 131 26h30 0,127 -
7 161 55,74 0,169 -
8 182 97,92 0,199 94,7 W (Cp = 0,192) a 170 rpm 114,7 W
9 215 160,2 0,228 (Cp = 0,16) a 170 rpm 135,8 W (Cp =
10 255 251,9 0,262 0,14) a 171 rpm 156,2 W (Cp = 0,12) a
11 258 304.4 0,238 170 rpm

A turbina eólica de eixo vertical começa a 3,5 m/s, e a potência aumenta à medida que a velocidade
do vento aumenta. A partir de 8 m/s, a saída máxima é gerada ao exceder a velocidade de rotação
projetada de 170 rpm. A saída máxima na velocidade do vento projetada de 9 m/s é de 160,2 W e a
velocidade de rotação é de 215 rpm. Quando a carga é controlada em 170 rpm, a saída é de 114,7 W, que
é maior que a saída desejada. O coeficiente de potência de projeto da turbina eólica é de 0,15 e a saída é maior do que
Machine Translated by Google

Energias 2018, 11, 1517 15 de 17

o coeficiente de potência projetado começando em 7 m/s. Quando a turbina eólica é operada a 9 m/s e
170 rpm, o coeficiente de potência é 0,163, que é maior que o valor de projeto de 0,15.
A Tabela 7 compara a potência de saída da turbina eólica prevista pelo modelo matemático
com a potência medida no experimento do túnel de vento. Na condição de projeto, a potência
medida foi de 114,7 W, 5,9% maior que a do modelo matemático. Este resultado valida o método de
projeto proposto e a estimativa de potência pelo modelo matemático pode ser útil para projetar um
VAWT de baixa velocidade com uma precisão razoável. A maior potência do teste de túnel de vento
pode ser resultado do efeito de confinamento devido às paredes do túnel de vento.

Tabela 7. Potência e coeficiente de potência obtidos a partir do modelo matemático e teste de túnel de vento.

Método Potência [W] Cp

Modelo matemático 108,34 0,154
Teste de túnel de vento 114,7 0,163

O tamanho da pá foi determinado pelas Equações (1)–(9), mas é impossível projetar a mudança do
aerofólio e o ângulo de torção da pá. No processo de estimativa de potência por meio do modelo
matemático, a saída pode ser prevista refletindo a sustentação e o arrasto do aerofólio [36].
No entanto, é impossível considerar a mudança de saída devido à esteira ocorrendo na faixa após 180
graus de azimute. Mais estudos devem ser feitos para investigar as estruturas de fluxo associadas à
lâmina helicoidal rotativa.

5. Conclusões

A fórmula básica do projeto rendeu as dimensões do rotor da turbina eólica com uma potência
aerodinâmica de 100 W a uma velocidade nominal do vento de 9 m/s e uma taxa de velocidade de ponta de
1,1. O torque devido à rotação do rotor pode ser calculado aplicando as forças de sustentação e arrasto
derivadas dos resultados CFD 2D. A saída média foi calculada como 108,34 W e a saída desejada de 100 W
foi atendida. A turbina projetada foi fabricada e um teste de túnel de vento foi realizado. A variação de saída
de acordo com a velocidade do rotor foi medida a cada velocidade do vento. Quando a velocidade do vento
incidente é de 9 m/s na velocidade de rotação de 170 rpm, a potência medida foi de 114,7 W e o método de projeto foi validado.
No entanto, o método de projeto não pode prever a variação de potência devido ao número de pás, ao
ângulo de torção da pá helicoidal, ao ângulo de passo e à posição do suporte. Mais pesquisas devem ser
realizadas para diferentes detalhes de geometria do rotor helicoidal.

Contribuições do autor: DH projetou, fabricou e testou o VAWT e escreveu o artigo. A YGH realizou simulação CFD 2D e analisou
os dados. A NJC forneceu a ideia do projeto e a SHN realizou experimentos em túneis de vento. KHC supervisionou a fabricação de
VAWT e KCK supervisionou o projeto e a papelaria.

Agradecimentos: Este trabalho foi apoiado pela concessão da National Research Foundation of Korea (NRF) financiada pelo
governo coreano (MSIT) através do Global Core Research Center for Ships and Offshore Plants (GCRC-SOP, No. 2011-0030013).

Conflitos de Interesse: Os autores declaram não haver conflito de interesse.

Referências

1. Tong, KC Aspectos técnicos e econômicos de um parque eólico offshore flutuante. J. Wind Eng Ind. Aerodyn. 1998,
74–76, 399–410. [CruzRef]
2. Ahmed Shata, AS; Hanitsch, R. Avaliação do potencial de energia eólica e geração de eletricidade na costa do Mar Mediterrâneo
no Egito. Renovar. Energy 2006, 31, 1183–1202. [CruzRef]
3. Associação Europeia de Energia Eólica. Energia Eólica e Política Climática da UE; EWEA: Bruxelas, Bélgica, 2011.
4. Greening, B.; Azapagic, A. Impactos ambientais de micro-turbinas eólicas e seu potencial para contribuir
à meta de mudanças climáticas do Reino Unido. Energia 2013, 59, 454–466. [CruzRef]
Machine Translated by Google

Energias 2018, 11, 1517 16 de 17

5. Syngellakis, K.; Clemente, P.; Cace, J. Questões Administrativas e de Planejamento para Pequenas Turbinas Eólicas em Áreas
Urbanas; Comissão Europeia: Bruxelas, Bélgica, 2006.
6. Ji, HS; Qiang, L.; Bico, JH; Mieremet, R.; Kim, KC Efeito da direção do vento nas estruturas próximas da esteira de uma pá de
turbina eólica em espiral de Arquimedes. J. Vis. 2016, 19, 653–665. [CruzRef]
7. Kim, KC; Ji, HS; Kim, YK; Lu, Q.; Bico, JH; Mieremet, R. Estudo experimental e numérico das características aerodinâmicas de uma
turbina eólica espiral de Arquimedes. Energies 2014, 7, 7893–7914. [CruzRef]
8. Scheurish, F.; Fletcher, TM; Brown, RE A influência da curvatura da pá e da torção da pá helicoidal no desempenho de uma turbina
eólica de eixo vertical. In Proceedings of the 48th AIAA Aerospace Sciences Meeting Incluindo o New Horizons Forum and
Aerospace Exposition, Orlando, FL, EUA, 4–7 de janeiro de 2010.
9. Hau, E. Turbinas eólicas, Fundamentos, Tecnologias, Aplicação, Economia, 2ª ed.; Springer: Berlim, Alemanha, 2006.
10. Manwell, JF; McGowan, JG; Rogers, AL Wind Energy Explained — Teoria, Projeto e Aplicação; John Wiley & Sons Ltd.: Hoboken,
NJ, EUA, 2002.
11. Choi, NJ; Nam, SH; Jeong, JH; Kim, KC Estudo numérico sobre o arranjo de turbinas de eixo horizontal em
um parque eólico: Efeito da distância de separação. J. Vento Eng. Ind. Aerodyn. 2013, 117, 11–17. [CruzRef]
12. Choi, NJ; Nam, SH; Jeong, JH; Kim, KC Estudo CFD de mudanças de potência aerodinâmica com variação de espaçamento entre
turbinas para um parque eólico offshore de 6 MW. Energias 2014, 17, 7483–7498. [CruzRef]
13. Heo, YG; Choi, NJ; Choi, KH; Ji, HS; Kim, KC Estudo CFD sobre a potência aerodinâmica de uma turbina eólica aumentada de
110 kW. Construção de energia. 2016, 129, 162–173. [CruzRef]
14. Darrieus, G. Turbine tendo seu eixo rotativo transversal ao fluxo da corrente. Patente dos EUA 1835018,
8 de dezembro de 1931.

15. Templin, RJ Teoria do Desempenho Aerodinâmico para Turbina Eólica de Eixo Vertical NRC, Laboratório Técnico
Relatório; LTR-LA-160; Conselho Nacional de Pesquisa, Canadá: Ottawa, ON, Canadá, 1974.
16. Strickland, JH A Turbina Darrieus: Um Modelo de Previsão de Desempenho Usando Múltiplos Streamtubes; Laudo Técnico
Laboratorial ; AREIA74-0431; Sandia National Laboratory: Livermore, CA, EUA, 1975.
17. Paraschivoiu, I. Modelo Double-Multiple Streamtube para estudar turbinas eólicas de eixo vertical.
J. Propuls. Power 1988, 4, 370-377. [CruzRef]
18. Islam, M.; Ting, DS-K.; Fartaj, A. Modelos aerodinâmicos para vento asix vertical de lâmina reta tipo Darrieus
turbinas. Renovar. Sus. Energy Rev. 2008, 12, 1087–1109. [CruzRef]
19. Sutherland, HJ; Berg, DE; Ashwill, TD Uma Retrospectiva da Tecnologia VAWT, Sandia Report; SAND2012-0304;
Sandia National Laboratories: Livermore, CA, EUA, 2012.
20. Wang, LB; Zhang, L.; Zeng, ND Um modelo de painel de vórtice 2-D de fluxo potencial: Aplicações à turbina de maré de lâmina
reta de eixo vertical. Conv. de energia Gerenciar 2007, 48, 454–461. [CruzRef]
21. Brusca, S.; Lanzafame, R.; Messina, M. Projeto de uma turbina eólica de eixo vertical: como a relação de aspecto afeta
desempenho da turbina. Int. J. Energy Environ. Eng. 2014, 5, 330–340. [CruzRef]
22. Sheldahl, RE Comparação de campo e túnel de vento Darrieus Wind Turbine Data. J. Energy 1981, 5, 254–256.
[CruzRef]
23. Bedon, G.; Castelli, MR; Benini, R. Ensaios experimentais de uma turbina eólica de eixo vertical com pás torcidas.
ICMIME 2013, 2013, 384-387.
24. Cheng, Q.; Liu, X.; Ji, HS; Kim, KC; Yang, B. Análise aerodinâmica de uma turbina eólica de eixo vertical helicoidal.
Energias 2017, 10, 575. [CrossRef]
25. Burton, T.; Sharpe, D.; Jenkins, N.; Bossanyi, E. Manual de Energia Eólica; John Wiley & Sons, Ltd.: West Sussex,
Reino Unido, 2001.

26. Tong, W. Geração de Energia Eólica e Projeto de Turbinas Eólicas; WIT Press: Ashurst, Reino Unido, 2010.
27. IEC 61400-2. Requisitos de Projeto para Pequenas Turbinas Eólicas, 2ª ed.; Padrão Internacional: Genebra,
Suíça, 2006.
28. Ray, ML; Rogers, AL; McGowan, JG Análise de modelos de cisalhamento do vento e tendências em diferentes terrenos.
Em Proceedings of the American Wind Energy Association Windpower, Amherst, MA, EUA, 29 de janeiro de 2014.
29. El Kasmi, A.; Masson, C. Um modelo k-ÿ estendido para fluxo turbulento através de turbinas eólicas de eixo horizontal.
J. Vento Eng. Ind. Aerodyn. 2008, 96, 103–122. [CruzRef]
30. Eriksson, S.; Bernhoff, H.; Leijon, M. Avaliação de diferentes conceitos de turbinas para energia eólica. Renovar. Sus.
Energy Rev. 2008, 12, 1419–1434. [CruzRef]
31. Abo-Khalil, AG; Lee, D.-C. Controle MPPT de Sistemas de Geração Eólica Baseado na Velocidade Estimada do Vento
Usando SVR. Ind. Elétron. IEEE 2008, 55, 1489-1490. [CruzRef]
Machine Translated by Google

Energias 2018, 11, 1517 17 de 17

32. IEC 61400-12-1. Turbinas eólicas: Medições de desempenho de energia de turbinas eólicas produtoras de eletricidade, 1ª ed.;
Padrão Internacional: Genebra, Suíça, 2005.
33. Chinchilla, M.; Arnaltes, S.; Burgos, JC Controle de geradores de ímãs permanentes aplicados a sistemas eólicos de velocidade
variável conectados à rede. Conv. de energia IEEE 2006, 21, 130–135. [CruzRef]
34. Muljadi, E. Butterfield, Geração de turbina eólica de velocidade variável controlada por Pitch. Ind. Appl. IEEE 2001, 37,
240–246. [CruzRef]
35. Wang, Q.; Chang, L. Um algoritmo inteligente de extração de potência máxima para velocidade variável baseada em inversor
sistemas de turbinas eólicas. Elétron de Potência. IEEE 2004, 19, 1242–1249. [CruzRef]
36. Karbasian, RH; Kim, KC Investigações numéricas sobre a estrutura do fluxo e o comportamento dos vórtices no estol dinâmico
de um hidrofólio de inclinação oscilante. Oceano Eng. 2016, 127, 200–211. [CruzRef]

© 2018 pelos autores. Licenciado MDPI, Basel, Suíça. Este artigo é um artigo de acesso aberto
distribuído sob os termos e condições da Creative Commons Attribution

(CC BY) (http://creativecommons.org/licenses/by/4.0/).