Desenvolvimento de Uma Turbina Eólica de Eixo Vertical
Desenvolvimento de Uma Turbina Eólica de Eixo Vertical
Desenvolvimento de Uma Turbina Eólica de Eixo Vertical
=5m/s ................. 18
Figura 7 ngulo de ataque ao longo da posio da p , em funo do coeficiente de velocidade
perifrica ............................................................................................................................................. 19
Figura 8 Representao de foras numa VAWT giromill para um perfil simtrico ......................... 22
Figura 9 - Coeficiente de sustentao C
L
em funo do ngulo de ataque para um perfil com
oscilao sinusoidal para diferentes frequncias
(27)
............................................................................. 23
Figura 10 Perfil NACA0015 em condies de perda dinmica
(30)
................................................... 24
Figura 11 Visualizao experimental do escoamento numa VAWT. Re=3x10
3
, =2, imagens com
80x80mm
(25)
.......................................................................................................................................... 24
Figura 12 Descrio da evoluo dos vrtices iniciais formados no bordo de ataque (1-a) e no
bordo de fuga (1-b) e dos vrtices secundrios formados no bordo de ataque (2-a) e no bordo de
fuga (2-b). Re=3x10
3 (25)
. ........................................................................................................................ 25
Figura 13 Contornos de vorticidade a 415 RPM, V
=4m/s, =4 .......... 55
Figura 28 Magnitude da velocidade do escoamento para a configurao 5, V
=4m/s, =4 (detalhe
na 4 rotao) ........................................................................................................................................ 56
Figura 29 Obteno das seces em carto canelado (esquerda) e colocao das seces
alinhadas numa base para a produo do molde das ps (direita) ...................................................... 65
Figura 30 Meia asa aps cura, retirada do molde (esquerda) e as duas meias asas prontas para
serem unidas para formar uma p (direita) .......................................................................................... 66
Figura 31 P finalizada ..................................................................................................................... 66
Figura 32 Ps feitas no Dep. de Aeroespacial finalizadas ................................................................ 67
Figura 33 Turbina finalizada (sem instrumentos) ............................................................................ 68
Figura 34 Detalhe do apoio do taqumetro, travo e brao das ps ............................................... 68
Figura 35 Detalhe do acoplamento das ps aos braos .................................................................. 68
Figura 36 Seco de p para ensaio de testes no-destrutivos. Apresentam-se 3 sensores
acsticos e o local de destruio por fadiga
(58)
..................................................................................... 70
iii
Figura 37 Resultados do ensaio no destrutivo de emisses acsticas. Os pontos negros indicam a
localizao calculada da fonte de emisso acstica
(58)
......................................................................... 70
Figura 38 Resultados do ensaio no destrutivo de leitura de infra-vermelhos. O aumento relativo
de temperatura de uma zona da p fornece importantes informaes sobre a concentrao de
esforos na mesma
(58)
........................................................................................................................... 70
Figura 39 Teste de mquina elctrica ............................................................................................. 77
Grficos
Grfico 1 Distribuio probabilstica da intensidade da velocidade do vento (esquerda) e da
direco do vento (direita) para Ponta Delgada
(8)
.................................................................................. 7
Grfico 2 Potncia elctrica elica instalada
(9)
................................................................................. 7
Grfico 3 - Valores tpicos para Cp
(6)
.................................................................................................. 13
Grfico 4 - Nmero de Reynolds mdio para vrias configuraes ................................................... 21
Grfico 5 Relao entre o alongamento (aqui apresentado em funo do raio) e o Cpmax para
uma turbina com 1 p NACA0015, =0,0833 e Re=360000
(36)
............................................................. 28
Grfico 6 Teste de Cp com diferentes end plates para uma VAWT subaqutica de 3 ps com
r=45cm, NACA 63(4)-021, =0,435, h=0,5m e V
=1,5 m/s
(36)
.............................................................. 28
Grfico 7 Potncia mdia a 5m/s vs potncia nominal em VAWTs de pequena dimenso ........... 42
Grfico 8 Preo em funo da potncia real em VAWTs de pequena dimenso............................ 42
Grfico 9 Preo em funo da rea para VAWTs de pequena dimenso ....................................... 43
Grfico 10 Potncia mdia a 5m/s em funo da rea A para VAWTs de pequena dimenso ...... 43
Grfico 11 Alongamento AR de VAWTs de pequena dimenso ...................................................... 44
Grfico 12 Comparao de resultados para o modelo k RNG
(35)
.................................................. 48
Grfico 13 - C
L
em funo de Re e (at 20) para vrios perfis
(42)
.................................................. 49
Grfico 14 C
D
em funo de Re e (at 20
o
) para vrios perfis
(42)
................................................. 50
Grfico 15 Evoluo dos resultados ao longo das rotaes da turbina .......................................... 54
Grfico 16 Resultados CFD para configurao 5 ............................................................................. 57
Grfico 17 Cp max em funo da velocidade do vento para todas as configuraes ..................... 58
Grfico 18 - |Fn| em funo de para configurao 5 com V
=8 m/s ............................................. 59
Grfico 19 - |Fn| em funo de para vrias configuraes a V
=8 m/s, =5 ................................. 59
Grfico 20 |Fn| em funo de para configurao 5 com V
=8 m/s ............................................. 60
Grfico 21 Medio 1 do binrio aplicado pelo travo ................................................................... 71
Grfico 22 Medio 2 do binrio aplicado pelo travo ................................................................... 72
Grfico 23 Medio 3 do binrio aplicado pelo travo ................................................................... 72
iv
Grfico 24 Funcionamento de VAWTs ............................................................................................ 73
Grfico 25 Resultados experimentais obtidos para a configurao 5 ............................................. 75
Grfico 26 Resultados experimentais obtidos para a Configurao 4 ............................................ 75
Quadros
Quadro 1 Caractersticas das turbinas da Figura 2
(2) (4)
................................................................... 10
Quadro 2 Vantagens e desvantagens de HAWTs e VAWTs
(2) (4) (10)
................................................. 12
Quadro 3 Evoluo das centrais elicas dos Aores
(20) (23)
.............................................................. 15
Quadro 4 - Potncia elica a instalar nos Aores at 2014
(20)
........................................................... 16
Quadro 5 Valores para
mx
e
max
em funo do coeficiente de velocidade perifrica ............. 19
Quadro 6 Valores de C
pmax
em funo de de acordo com a teoria de Glauert-Schmitz ............... 27
Quadro 7 Variao de preos na mesma turbina ............................................................................ 42
Quadro 8 Descrio das configuraes simuladas .......................................................................... 57
Quadro 9 Resultados experimentais e numricos .......................................................................... 76
v
Lista de abreviaes
A rea;
AR Alongamento;
c
Corda da p;
D
C
Coeficiente de resistncia;
L
C
Coeficiente de sustentao;
P
C
Coeficiente de potncia;
COE Custo da energia (do ingls Cost Of Energy);
D
Resistencia;
EWEA
European Wind Energy Association;
f Frequncia;
F
Fora;
n
F Fora normal;
t
F Fora tangencial;
HAWT
Turbina elica de eixo horizontal (do ingls Horizontal Axis Wind Turbine);
I
Momento de inrcia;
IST
Instituto Superior Tcnico;
L
Sustentao;
LC Custos totais ao longo da vida (do ingls Lifetime Costs);
LEO Energia produzida ao longo da vida (do ingls Lifetime Energy Output);
M Binrio;
n Nmero de ps;
P Potncia;
r Raio;
Re Nmero de Reynolds;
UE Unio Europeia;
V Velocidade;
VAWT Turbina elica de eixo vertical (do ingls Vertical Axis Wind Turbine);
in cut
V
_
Velocidade de arranque;
out cut
V
_
Velocidade mxima;
vi
i
V Velocidade induzida;
p
V Velocidade do escoamento na p;
rated
V Velocidade nominal;
V Velocidade do vento;
o ngulo de ataque;
u Posio angular;
Coeficiente de velocidade perifrica (em ingls tip speed ratio);
Densidade;
o Factor de bloqueamento (em ingls solidity);
e Velocidade angular;
u
Viscosidade cinemtica.
1
1 Apresentao
1.1 Objectivos e enquadramento
O desenvolvimento desta tese baseou-se inicialmente na proposta de uma empresa para a criao
de uma unidade hbrida de fornecimento de energia elctrica, composta por um sistema fotovoltaico
e elico. Pretendia-se estudar a viabilidade da integrao de um sistema destes num poste de
iluminao, sem ligao rede elctrica. Do estudo preliminar efectuado concluiu-se que tal sistema
no economicamente vivel face aos postes de iluminao convencionais em locais com acesso
rede elctrica. Esta concluso levou ao abandono da ideia de desenvolvimento deste sistema.
No entanto, deste estudo, resultou a ideia de que as turbinas elicas de pequena dimenso teriam
preos relativamente elevados. Esta ideia levou proposta, da mesma empresa, para o
desenvolvimento de uma turbina elica de eixo vertical (VAWT) de pequenas dimenses e baixo
custo.
O processo de construo escolhido teria de ser flexvel devido incerteza no ritmo de produo
inicial e minorando aumentos de custos na empresa como a compra de maquinaria dedicada. Outros
factores da maior importncia seriam uma baixa manuteno, fcil instalao, montagem e
transporte. Baixo rudo e uma esttica apelativa seriam tambm factores importantes a considerar.
Entretanto, a empresa resolveu abandonar o projecto. Da em diante este viria a ser desenvolvido
como tese de mestrado, por parte do IST. A escolha do Departamento de aeroespacial justificada
pela similaridade de reas abordadas entre a indstria aeronutica e as turbinas elicas,
nomeadamente o estudo de asas. Em outras universidades como TuDelft ou Griffith, a relao entre
a indstria elica e a indstria aeronutica tem-se aprofundado pelo mesmo motivo
(1)
.
No incio do desenvolvimento do trabalho havia alguma falta de conhecimento em relao a este
tema. Inicialmente, o objectivo seria o desenvolvimento completo e tendencialmente ptimo de uma
VAWT. Com a reviso bibliogrfica efectuada e tendo em conta o estado da arte concluiu-se que esta
uma rea num estado prematuro de desenvolvimento e investigao
(2)
. Assim sendo, no seria
2
possvel, no contexto de uma tese de mestrado, executar os objectivos inicialmente propostos. Desta
forma, os objectivos foram adaptados. Aps esta adaptao inicial, fixou-se como objectivo desta
tese efectuar o estudo preliminar para o desenvolvimento de uma VAWT e fornecer uma base de
trabalho para o desenvolvimento de outros trabalhos nesta rea.
Para cumprir o objectivo proposto, procedeu-se ao estudo de vrios parmetros considerados
fundamentais no desenvolvimento de VAWTs como os tipos de perfis usados nas ps, o nmero de
ps ou o alongamento das mesmas
(2)
. Iniciou-se o estudo atravs de uma reviso bibliogrfica e uma
sntese do estado da arte. De seguida procedeu-se ao desenvolvimento de modelos tericos e
obteno de resultados. Recolheram-se dados de modelos comerciais de VAWTs actualmente
disponveis para estudar as tendncias da indstria e fazer a comparao com os resultados obtidos
e o estado da arte. Foram obtidos resultados numricos para vrias configuraes de VAWTs usando
o modelo kw SST transitional com o software Fluent. Por fim construiu-se um prottipo
propositadamente para este trabalho a partir do qual se obtiveram resultados experimentais que
foram usados para a validao dos resultados numricos e tericos obtidos.
3
1.2 Estrutura e organizao do texto
Este texto est organizado em 8 captulos que reflectem aproximadamente a evoluo cronolgica
do trabalho desenvolvido.
O captulo 1 apresenta os objectivos do trabalho, enquadra o seu desenvolvimento e caracteriza a
estrutura e a organizao do texto.
O captulo 2 fornece uma introduo energia elica, que se pretendeu tecnicamente simples.
feita uma breve reviso histrica e uma apreciao de alguns conceitos simples mas fundamentais
como a caracterizao do vento e a sua potncia. De seguida faz-se uma breve reviso dos sistemas
usados para extraco e aproveitamento da energia elica. Termina-se com uma apresentao sobre
a energia elica nos Aores, cobrindo a sua evoluo histrica, tecnologias utilizadas e perspectivas
futuras.
No captulo 3 prossegue-se a sntese do estado da arte j iniciada no captulo 2 mas com maior
nfase tcnico e com total destaque para o desenvolvimento de turbinas elicas de eixo vertical.
Dada a extenso da informao apresentada, dividiu-se o captulo em 11 subcaptulos. No primeiro
subcaptulo faz-se uma introduo e nos seguintes apresentado o estado da arte para cada
parmetro de funcionamento considerado.
No captulo 4 faz-se uma reviso de modelos comerciais disponveis de VAWTs. A partir desta
reviso constri-se uma base de dados que usada para a obteno e comparao de vrios
parmetros.
No captulo 5 apresenta-se a anlise aerodinmica de VAWTs. O captulo encontra-se divido em 2
subcaptulos. No primeiro subcaptulo faz-se uma sntese do estado da arte enquanto no segundo
subcaptulo apresenta-se o trabalho efectuado nesta rea.
O captulo 6 dedica-se ao tema da construo de turbinas elicas. No primeiro subcaptulo faz-se
uma sntese do estado da arte e no segundo subcaptulo descreve-se a construo do prottipo
desenvolvido propositadamente para este trabalho.
O captulo 7 dedica-se ao teste de VAWTs. Seguindo a estrutura dos captulos anteriores, o captulo
encontra-se dividido em dois subcaptulos. No primeiro subcaptulo feito um resumo do estado da
arte e no segundo subcaptulo descrito o teste e os resultados obtidos com o prottipo construdo.
O texto termina no captulo 8 que se encontra dividido em 3 subcaptulos. No primeiro subcaptulo
apresentam-se as concluses do trabalho efectuado, no segundo sugerem-se alguns trabalhos a
desenvolver no futuro e no terceiro faz-se uma breve apreciao da reaco das pessoas ao
prottipo em funcionamento.
4
2 Introduo
2.1 O vento e a energia elica
Energia elica a energia do vento que tem sido usada h sculos como uma fonte de energia pela
humanidade. Esta forma de energia pode ser convertida em formas mais teis de energia atravs,
por exemplo, de turbinas elicas para a produo de energia elctrica. Outras aplicaes referem-se
transmisso directa de energia mecnica da turbina para o bombeamento de gua ou arejamento
de tanques de aquicultura
(3)
.
J no sc. XVII a.C., os babilnios usaram moinhos de vento em projectos de irrigao. Mais tarde,
inventores Persas trouxeram inovaes e novos conceitos e, no sc. VII, gegrafos rabes relataram a
existncia de moinhos de vento semelhantes s modernas portas giratrias. Estes moinhos de eixo
vertical ainda existiam em pases como o Afeganisto em 1980, estimando-se que teriam uma
potncia de 75 hp e conseguiriam moer uma tonelada de trigo em 24h
(3)
.
Os moinhos de vento europeus datam do sc. XII. Os sistemas originais tinham de ser alinhados
manualmente com o vento mas rapidamente surgiram inovaes e o moinho de vento Europeu
sofreu um grande desenvolvimento e prosperou: no fim do sc. XIX existiam mais de 30000 moinhos
de vento na Europa, usados principalmente para a moagem de gro e o bombeamento de gua
(3) (4)
.
Uma das primeiras tentativas para gerar electricidade a partir do vento, com uma turbina elica de
eixo vertical (VAWT), foi feita nos EUA em 1888 por Charles Brush. Em 1922, S.J. Savonius apresentou
a VAWT Savonius e em 1931 Georges Darrieus patenteou a VAWT Darrieus e giromill
(4)
.
O aparecimento do motor a vapor seguido do desenvolvimento do motor trmico parecia ter
ditado o fim da energia elica. Porm, na dcada de 60, este conceito ressurgiu e foi finalmente
aceite no incio da dcada de 90. Vrios factores foram responsveis por esta transformao, como a
crescente necessidade de energia, a compreenso do limite de utilizao de combustveis fsseis e os
seus efeitos adversos no ambiente. Havia o potencial representado pelo vento, disponvel em todo o
lado e j comprovadamente utilizado para o transporte e potncia mecnica. Uma revoluo
tecnolgica permitia a utilizao de novas tecnologias, como materiais compsitos, aumentando a
5
capacidade de produo elica e a reduo de custos. Por fim, a vontade poltica de implementao
destas tecnologias permitiu apoios para o seu desenvolvimento e investigao, alteraes legais para
permitir a sua interconexo na rede de distribuio e apoios para a sua utilizao
(2)
.
O vento, na Terra, consiste no movimento do ar devido a diferenciais de presso atmosfrica. Esses
diferenciais so causados principalmente pelo irregular efeito da radiao solar e pelo movimento de
rotao do prprio planeta
(2)
. Perto da superfcie terrestre, a velocidade do vento diminui por atrito
ocorrendo um diferencial vertical de velocidade do vento. Existem vrios modelos para a sua
descrio, sendo comum a utilizao de um perfil logartmico
(2)
.
As caractersticas locais do vento dependem essencialmente das caractersticas geomorfolgicas do
local. Nos locais terrestres, prximos de gua, o vento tende a soprar, de dia, da gua para terra e
noite em sentido inverso. Tal deve-se diferena entre a capacidade trmica mssica dos meios,
sendo superior na gua, fazendo com que esta varie a sua temperatura mais lentamente que o solo.
Nos locais com orografia acentuada, as caractersticas do vento so essencialmente locais e
principalmente ditadas pela topografia do terreno. Por exemplo, no caso de vales, durante o dia, a
incidncia solar nas montanhas circundantes rapidamente aquece o ar frio prximo destas causando
ventos locais do vale para a montanha. Durante a noite, as montanhas circundantes arrefecem mais
rapidamente que o vale, fazendo com que o vento predominante seja da montanha para o vale
(5)
.
O potencial elico a potncia dada pela energia cintica do vento numa dada rea por unidade de
tempo:
3
2
1
= AV P
elica
,
(1)
sendo
elica
P a potncia (potencial elico) em W,
a densidade do ar em kg/m
3
,
V a velocidade
do vento em m/s e A
a rea perpendicular ao escoamento considerada em m
2 (2)
. Note-se a
importncia de
V
com uma dependncia cbica.
diminui com a altitude. De acordo com a atmosfera ISA, varia entre 1,2250kg/m
3
ao nvel do
mar para 1,0065kg/m
3
a 2 km de altitude
(6)
. Note-se tambm a importncia da humidade relativa e
da temperatura do ar cujos aumentos fazem diminuir .
Existem vrias cartas que apresentam o potencial elico entre as quais a apresentada na Figura 1.
A maior parte destas cartas foram desenvolvidas por organizaes governamentais e apresentam
condies muito generalistas resultando de extrapolaes de medies locais e considerando
modelos muito simples de topografia. Outras foram desenvolvidas por rastreamento por satlite com
uma resoluo, geralmente, muito deficiente.
Existem vrios programas que se dedicam simulao das condies de vento com maior detalhe.
Um exemplo de um destes programas uma das componentes do programa Green Islands onde
6
usando o software dinamarqus WAsP (Wind Atlas analysis and application Program), uma base de
dados da topografia das ilhas dos Aores com uma elevada resoluo e os dados recolhidos nas
estaes do Instituto de Meteorologia, pretende-se obter uma carta que permita avaliar o potencial
deste recurso na regio. Tal permitir uma melhor orientao na poltica energtica regional.
Figura 1 Potencial elico por unidade de rea na Europa
(7)
A medio das caractersticas do vento num local feita com recurso a equipamentos
especializados que medem normalmente 4 parmetros: A velocidade do vento (intensidade e
direco), a temperatura e a presso atmosfrica. A frequncia f de amostragem, tambm outra
componente muito importante e s com frequncias elevadas se pode avaliar em detalhe a
turbulncia do vento de um local. Os sistemas de medio usados pelas agncias meteorolgicas
muitas vezes utilizam frequncias demasiado baixas para o estudo especfico do vento para produo
de electricidade. No entanto, a sua elevada disponibilidade e longos registos histricos (note-se que
o perodo mnimo de amostragem em climatologia define-se como 5 anos) tornam-nos um elemento
base em muitos estudos. Nesse caso da maior importncia a informao sobre as caractersticas
dos dados recolhidos como a altitude, escala de tempo ou calibrao.
7
0%
5%
10%
15%
20%
25%
N
E
S
W
Os dados recolhidos so normalmente apresentados atravs da distribuio probabilstica da
intensidade da velocidade do vento e a distribuio probabilstica da direco do vento como o
Grfico 1.
Grfico 1 Distribuio probabilstica da intensidade da velocidade do vento (esquerda) e da direco do vento
(direita) para Ponta Delgada
(8)
A nvel mundial, o aproveitamento de energia elica para produo de electricidade tem
aumentado cerca de 25% anualmente, desde o incio da dcada de 90 (Grfico 2). Esse aumento tem
sido mais significativo na Europa e essencialmente em Espanha, Dinamarca e Alemanha. No futuro
espera-se que esta tendncia se mantenha principalmente com o aproveitamento de rea offshore
(rea no mar)
(9)
.
Grfico 2 Potncia elctrica elica instalada
(9)
O desenvolvimento desta indstria resultou na criao de mais de 150000 postos de trabalho
(directos e indirectos) na Europa em 2007, a maioria deles com elevada formao. O rpido
desenvolvimento da indstria elica tem levado falta de profissionais qualificados principalmente
em Espanha, Alemanha e Dinamarca (onde se concentram 75% dos postos de trabalho). Espera-se
0%
5%
10%
15%
20%
25%
30%
35%
40%
45%
0 1,5 4,1 6,7 9,3 11,8 14,4 17,0
P
r
o
b
a
b
i
l
i
d
a
d
e
Intensidade da Velocidade do vento m/s
Jan
Jul
0
20
40
60
80
100
120
1
9
9
0
1
9
9
1
1
9
9
2
1
9
9
3
1
9
9
4
1
9
9
5
1
9
9
6
1
9
9
7
1
9
9
8
1
9
9
9
2
0
0
0
2
0
0
1
2
0
0
2
2
0
0
3
2
0
0
4
2
0
0
5
2
0
0
6
2
0
0
7
2
0
0
8
P
o
t
n
c
i
a
e
l
c
t
r
i
c
a
e
l
i
c
a
i
n
s
t
a
l
a
d
a
G
W
Ano
UE Resto do mundo
8
que em Portugal, o nmero de empregos nesta rea venha a disparar, dada a atractividade do
mercado, a existncia de profissionais qualificados e os baixos nveis salariais
(10)
.
A eficincia das turbinas elicas normalmente caracterizada pelo coeficiente de potncia que
representa a fraco de potncia do vento que extrada pela turbina:
3
2
1
= =
AV
P
P
P
C
turbina
elica
turbina
P
,
(2)
sendo
P
C
o coeficiente de potncia e
turbina
P
a potncia no rotor da turbina
(6)
.
O custo da energia de acordo com a sua origem num dado tempo COE (do ingls Cost Of Energy)
normalmente obtido em cntimos de dlares ou cntimos de euros por kWh de energia. O clculo
feito pela diviso dos custos totais, ajustados para o mesmo perodo de tempo LC (do ingls
Lifetime Costs) pela energia total obtida ao longo da vida de um equipamento LEO (do ingls
Lifetime Expected Output)
1 (9) (11)
:
LEO
LC
COE =
(3)
Os LC de um equipamento compreendem as seguintes componentes:
- Equipamento;
- Instalao;
- Transporte;
- Manuteno;
- Financiamento;
- Combustvel;
- Desmantelamento;
- Impostos ou apoios;
- Licenciamento
(11)
.
O clculo e anlise destas componentes so muito difceis alm de que existem muitos parmetros,
como o preo do petrleo ou taxas de carbono, que permitem um ajustamento conforme cada
estudo fazendo com que os resultados divulgados sejam muito variveis, tornando difcil uma
correcta anlise
(11)
. Walford
(12)
apresenta uma descrio de mtodos para a compreenso e reduo
de custos com a manuteno de turbinas elicas. Note-se que, por vezes, apresentado o custo da
energia anual, no tendo em conta a variao dos custos ao longo da vida de um equipamento
(normalmente aumento da manuteno e diminuio da eficincia) ou algumas componentes como
o desmantelamento. A insensibilidade para esta questo leva muitas vezes a escolhas incorrectas
sobre o tipo de equipamento a instalar.
1
A falta de termos universais para a definio destes parmetros levou s denominaes indicadas pelo autor mas na
literatura outras podem ser encontradas como, por exemplo, LEC (do ingls Lifetime Energy Capture) ao invs de LEO
(24)
.
9
Existe uma tendncia clara para a diminuio do COE das energias renovveis com a energia
elica, hdrica, solar trmica (considerada como energia primria e no para a produo de
electricidade), geotrmica e biomassa mostrando alguma competitividade actualmente. Com o
aumento dos custos dos combustveis fsseis e a aplicao de taxas sobre a libertao de gases de
efeito de estufa, muitas destas tecnologias tornar-se-iam prtica corrente
(2) (11)
. Estudos apontam
tambm para uma forte influncia da componente elica de produo de energia no preo da
electricidade. Este efeito, denominado merit order effect, ocorre pelo facto da produo de energia
elctrica de origem elica ter elevados custos iniciais (aproximadamente 75%) mas baixos custos
durante o seu funcionamento e um baixo custo marginal (mudana no custo de produo com a
variao da quantidade produzida). Tal leva a que, dentro dos limites de capacidade de produo de
energia elica, o preo se mantenha relativamente constante
(13) (11)
. Com uma maior coeso das
diferentes redes de distribuio de energia na Europa seria possvel uma maior diminuio do custo
da electricidade com a transmisso de electricidade de origem elica entre os vrios pases
(14)
. O
COE da energia elica tem vindo a diminuir de 11c/kWh em 1987 para 5c/kWh em 2004 em
zonas terrestres (custos ajustados a 2006)
(11)
. Para comparao, de acordo com o mesmo estudo, o
COE para centrais trmicas a carvo de 3,5c/kWh e 4,5c/kWh para centrais trmicas a gs
natural sem ciclo combinado (custos ajustados a 2006)
(11)
. Para mais dados, consultar os estudos
Mott Macdonald
(15)
ou California energy comission
(16)
.
Existem tambm muitos factores sociais ou econmicos indirectos associados aos diferentes tipos
de energia sendo o mais bvio a poluio mas sendo tambm importante referir a produo
descentralizada (no caso, por exemplo, da energia elica ou de mini-hdricas) e o consequente
aumento do equilbrio da distribuio de riqueza e diminuio do preo da energia produzida (por
competitividade entre concorrentes). A armazenagem de gua no caso da energia hdrica ou o
aproveitamento de materiais anteriormente considerados como lixo (como na biomassa),
incrementa o rendimento de certas actividades, principalmente no sector primrio
(2)
. Dados estes
resultados, espera-se que as energias renovveis, e em especial a energia elica, venham a ter uma
cada vez maior contribuio para a produo energtica mundial
(9) (5) (2)
.
10
2.2 Tecnologia e desenvolvimento
As turbinas elicas podem ser de dois tipos: tipo drag (resistncia) e tipo lift (sustentao). Nas
turbinas tipo drag, a fora motriz a resistncia aerodinmica enquanto nas de tipo lift a
sustentao aerodinmica. A utilizao de cada conceito implica diferenas no seu desenvolvimento
(2)
. Um dos parmetros mais importantes a relao entre a intensidade da velocidade das ps da
turbina e a velocidade do vento
=
V
r e
,
(4)
sendo e
a velocidade angular da turbina em rad/s e
r o raio em m
(2)
.
As turbinas tipo drag, teoricamente, nunca circulam mais rpido que o escoamento ( 1 < ), ao
contrrio do que acontece nas de tipo lift. As turbinas tipo drag tm normalmente um arranque a
V
mais baixo mas tm grande resistncia aerodinmica a
V e regime de
funcionamento constantes.
A turbina colocada a uma altura considervel,
exigindo uma maior resistncia estrutural do poste e
dificultando a manuteno.
Requer sistema de alinhamento com o vento.
O movimento das ps e do rotor funcionam como
um giroscpio gigante exercendo grandes foras no
caso de mudana de direco do vento.
Sistema normalmente mais ruidoso por implicar um
maior e pela origem do rudo estar mais alta.
VAWT
Vantagens Desvantagens
No necessita de alinhamento com o vento.
Componentes, como, por exemplo, o gerador, mais
prximas do cho, facilitando a montagem e
manuteno e exigindo uma menor resistncia
estrutural do poste.
Normalmente menos ruidoso por ter um menor
e pela origem do rudo estar prxima do solo.
Sistema ainda pouco desenvolvido permitido
maiores retornos para investigao.
A colocao a baixa altitude no permite aproveitar
ventos muito intensos.
Cargas aerodinmicas cclicas, induzindo fadiga ou
mesmo destruio (embora tal tambm acontea nas
HAWT mas de forma mais limitada e controlvel).
Maior rea de ps (uma componente normalmente
cara) que nas HAWT para a mesma potncia.
Requer, normalmente, sistema de arranque (para
VAWTs tipo lift)
(1)
.
Quadro 2 Vantagens e desvantagens de HAWTs e VAWTs
(2) (4) (10)
A indstria de energia elica tem tido um grande desenvolvimento recentemente, fruto da procura
e apoios prestados. Apesar de todo o desenvolvimento, errado considerar esta como uma indstria
madura sendo que ainda h muitas incertezas ao nvel das vrias reas intervenientes como, por
exemplo, a meteorologia, aerodinmica ou materiais. Outras reas mais especficas tais como
estratgias de manuteno, desenho de parques elicos ou integrao de redes inteligentes chegam
a ser contemporneas e, em parte, consequncia deste grande desenvolvimento. Subsistem ainda
tantas reas a investigar que actualmente no bvia a escolha ptima das dimenses de uma
turbina elica para um dado local
(19)
.
A necessidade de maior desenvolvimento nesta rea chegou a ser decretada pela Unio Europeia
UE com vista a facilitar a introduo de energias renovveis na Europa.
Trs tendncias tm dominado o mercado:
- Aumento das dimenses das turbinas elicas;
- Aumento da eficincia das turbinas elicas;
- Diminuio do COE (embora esta tendncia se tenha alterado muito recentemente, tal
como previsto em muitos estudos, que previam que as cargas cclicas nas ps de uma HAWT
levariam a um aumento do COE com o grande aumento das dimenses das mesmas)
(7)
.
13
Ao nvel das turbinas de pequenas dimenses, tem havido uma grande expresso por parte das
pessoas em geral na construo e utilizao destes sistemas, particularmente do tipo VAWT
Savonius, facilmente explicvel pelo conceito simples de funcionamento e pela sua facilidade de
construo (colocando as duas partes de um barril cortado em sentidos opostos fixas a um veio
vertical). Ao nvel comercial, j h muito tempo existem solues para mercados muito especficos
(como, por exemplo, zonas isoladas e barcos) mas comeam agora a aparecer sistemas
economicamente mais viveis face electricidade da rede. Estas turbinas so actualmente adaptadas
para outras situaes como para o aproveitamento da energia das mars ou de energia hdrica como
o exemplo da turbina GHT (do ingls Gorlov helical turbine).
O limite de Betz o valor mximo terico de
P
C
(6)
:
% 59
27
16
max
~ =
P
C .
(5)
Foi calculado considerando uma turbina ideal, com infinito nmero de ps, sem perdas. Na
realidade, vrios factores levam a que as turbinas actuais no atinjam este valor como a existncia de
atrito aerodinmico, a utilizao de um nmero finito de ps, as perdas nas pontas das mesmas
(end/tip losses), o momento transmitido ao escoamento que passa pela turbina, entre outros
(5).
Analisando o momento transmitido ao escoamento que passa pela turbina, pode-se obter um valor
de
max P
C mais detalhado como mostra o Grfico 3 (apresentado como limite de Betz corrigido). As
demonstraes destes resultados podem ser encontradas em detalhe nas referncias apresentadas,
sendo aqui apenas expostos os resultados para melhor compreenso de algumas escolhas e
avaliaes que sero feitas mais frente. Note-se que os resultados apresentados pelo Grfico 3
esto algo desfasados dos resultados actuais, tendo as VAWTs Savonius actuais uma curva
semelhante american multiblade
(17)
.
Grfico 3 - Valores tpicos para Cp
(6)
Nas redes elctricas actuais, o consumo dirio de electricidade tem fortes oscilaes ocorrendo
mnimos durante a noite, aproximadamente entre as 24 e as 5h, e mximos entre as 19 e as 22h.
14
Estas oscilaes so indesejveis pois implicam a utilizao parcial de geradores conforme a
necessidade, diminuindo a sua rentabilidade. Uma soluo usual a utilizao de tarifas multi-
horrias pelas empresas de distribuio de electricidade. Aplicando-se maiores taxas no consumo de
electricidade durante os perodos de maior procura e fazendo o inverso nos perodos de menor
procura, altera-se o comportamento do consumidor diminuindo as oscilaes no consumo
(20)
.
Com a introduo de energia renovvel na rede elctrica torna-se necessrio controlar as
oscilaes na produo de electricidade. Para a diminuio destas oscilaes, aposta-se na
diversificao das fontes energticas e na previso de condies de produo. De resto, necessrio
a implementao de sistemas de armazenagem de energia desde que economicamente viveis
(20)
.
Em pequenas aplicaes usam-se normalmente sistemas de baterias mas a sua aplicao em
grandes dimenses no efectuada devido ao seu elevado custo. Existe muita investigao a ser
efectuada nesta rea actualmente. Prev-se tambm, no futuro, a utilizao indirecta de baterias no
armazenamento de energia atravs da utilizao de veculos elctricos inteligentes. Uma aplicao
destas poderia levar ao equilbrio da procura de electricidade diria e a um novo paradigma para o
automvel passando a ser uma fonte de receita, comprando electricidade durante a noite, quando
esta se encontra a preos reduzidos, e vendendo parte durante o dia a preos mais elevados.
Outra tecnologia usada a compresso de ar. Embora actualmente no seja uma tecnologia muito
aplicada, tem uma elevada investigao associada. A utilizao de sistemas de inrcia, preservando a
energia como energia cintica tem sido aplicada com grande sucesso em sistemas com uma
frequncia de oscilao elevada como no caso das turbinas elicas, sendo que actualmente, parte
das turbinas elicas de grandes dimenses trazem incorporados sistemas desta natureza. A
possibilidade de produzir e armazenar hidrognio para reconverso futura em electricidade vista
como uma tecnologia de grande potencial principalmente se a adopo do carro a hidrognio se
realizar. uma tecnologia que no se encontra disseminada devido s perdas de energia durante o
processo, s dificuldades no armazenamento e transporte do hidrognio e aos custos associados,
principalmente, s clulas de combustvel.
Por fim, existe o armazenamento de energia potencial nas barragens seja pelo controlo de gua
escoada ou at pelo bombeamento de gua de jusante para montante, armazenando-a para
utilizao futura. Estes sistemas de bombeamento so utilizados em conjunto com energia elica ou
utilizando o excesso de produo elctrica de origem nuclear. actualmente o sistema com maior
disseminao para grandes armazenamentos, estando a tecnologia disponvel e testada
(20)
.
Embora aqui se tenha tentado fornecer uma introduo sucinta energia elica e s tecnologias
envolvidas, para mais informao, uma excelente introduo energia elica dada por Sahin
(21)
.
Descries mais exaustivas so dadas por Manwell
(5)
ou EWEA
(7)
.
15
2.3 A energia elica nos Aores
O aproveitamento da energia elica para produo de electricidade teve incio, nos Aores, em
1988 com a instalao do parque elico do Figueiral, na ilha de Santa Maria, constitudo por 9
aerogeradores de 30 kW. Os resultados positivos obtidos e o empenhamento do Governo Regional
nas energias renovveis levaram instalao de numerosos parques como indicado no Quadro 3.
As instalaes esto sob a superviso do grupo EDA Electricidade dos Aores, empresa
responsvel pela produo e distribuio de electricidade no arquiplago
(20) (22)
.
Nome Localizao
Potncia
total MW
N
turbinas
Potncia
nominal kW
Fornecedor Modelo
Ligao
rede
Serra do Cume Terceira 4,5 5 900 Enercon E-44 2008
Pico Pico 1,8 6 300 Enercon E30 2005
Lomba dos frades Faial 1,8 6 300 Enercon E30 2002
Figueiral Santa
Maria
0,9 3 300 Enercon E30 2002
Boca da vereda Flores 0,6 2 300 Enercon E30 2002
Serra branca
(ampliao)
Graciosa 0,6 2 300 Enercon E30 2002
Pico da urze (2
ampliao)
So Jorge 0,6 2 300 Enercon E30 2002
Santa Maria
(desactivao)
Santa
Maria
-0,27 -9 30 AEROMAN 14/30 2001
Pico da urze
(ampliao)
So Jorge 0,15 1 150 NORDTANK NTK150 1994
Serra branca Graciosa 0,2 2 100 NORDTANK NTK150 1992
Pico da urze So Jorge 0,40 4 100 NORDTANK NTK150 1991
Santa Maria Santa
Maria
0,27 9 30 AEROMAN 14/30 1988
Quadro 3 Evoluo das centrais elicas dos Aores
(20) (23)
A poltica energtica da regio tem sofrido alteraes tendo-se criado metas ambiciosas de 50% de
energia elctrica renovvel em 2015 e de 75% em 2018. Para alcanar estes objectivos, o governo
regional tem apoiado trs grandes vertentes:
- A actualizao dos parques energticos do grupo EDA para produo renovvel;
- A liberalizao do mercado de micro-gerao permitindo e, essencialmente,
desburocratizando a venda de energia rede por particulares (PROENERGIA);
- O apoio investigao no domnio das energias renovveis e do potencial dos Aores nesta
rea, com o apoio do programa MITPortugal.
Todas estas iniciativas resultaram na formao do programa Green Islands e de outros organismos
reguladores como o ARENA. Importa tambm salientar outros mecanismos como o caso da
certificao energtica de edifcios
(20)
.
Tendo em conta esta poltica energtica, o grupo EDA tem procedido ampliao e criao de
centrais hdricas, elicas (Quadro 4) e geotrmicas.
16
Nome Localizao Potncia instalada kW Potncia a instalar kW Potncia final kW
Figueiral Santa Maria 900 660 1560
Graminhais So Miguel 0 9000 9000
Serra do cume Terceira 4500 4500 9000
Serra branca Graciosa 800 660 1260
Pico da urze So Jorge 1150 990 2140
Pico Pico 1800 660 2460
Faial Faial 0 4500 4500
Boca da Vereda Flores 600 0 600
Total 9750 20970 30520
Quadro 4 - Potncia elica a instalar nos Aores at 2014
(20)
Estes investimentos esto j contabilizados na Figura 3.
Figura 3 - Origem da energia elctrica produzida nos Aores
(20) (22)
A elevada penetrao de energias renovveis nos Aores e as oscilaes de produo que estas
apresentam, levou necessidade de implementao de sistemas de regulao. Tais sistemas so,
como se viu anteriormente, normalmente necessrios, quando existe uma forte componente de
energias renovveis. No caso dos Aores, recorreu-se a sistemas de inrcia nas ilhas com maior
penetrao (Flores e Graciosa). Note-se que outras ilhas como So Miguel, que tambm tm uma
elevada taxa de penetrao, utilizam fontes de energia renovvel mais estveis e previsveis, como a
energia hdrica ou geotrmica, no necessitando de sistemas de estabilizao da rede elctrica. Esta
utilizao tem demonstrado a possibilidade de utilizao de elevadas penetraes de energias
renovveis em redes elctricas de vrias dimenses, tornando os Aores num laboratrio vivo,
exemplo para todo o mundo
(20)
.
Com o desenvolvimento arquitectado at 2014, espera-se que a penetrao das energias
renovveis atinja valores ainda maiores (50%), em redes elctricas maiores (como na ilha de S.
Miguel com uma populao de 130000 habitantes) e com elevada diversificao. Considerando
novamente S. Miguel, ter-se- energia elica, hdrica, geotrmica e fotovoltaica por particulares. Este
ser um grande ponto de desenvolvimento e esto a ser preparados novos mecanismos de regulao
como hdricas reversveis e instalao de redes inteligentes
(20)
.
0%
20%
40%
60%
80%
100%
2
0
0
0
2
0
0
1
2
0
0
2
2
0
0
3
2
0
0
4
2
0
0
5
2
0
0
6
2
0
0
7
2
0
0
8
2
0
0
9
2
0
1
4
O
r
i
g
e
m
d
e
e
n
e
r
g
i
a
e
l
c
t
r
i
c
a
p
r
o
d
u
z
i
d
a
Ano
Geotrmica Hdrica Elica No-renovvel
17
3 Parmetros de funcionamento
3.1 Introduo aos sistemas VAWT
As turbinas elicas tm uma curva de potncia caracterstica como mostra a Figura 4.
Figura 4 Curva de potncia caracterstica de turbinas elicas
(24)
Na zona I, a velocidade do vento demasiado baixa para que ocorra a produo de energia. A
partir da velocidade de arranque,
in cut
V
_
, na zona II, a turbina inicia a produo de energia mas com
valores inferiores produo nominal. Nesta zona (II), teoricamente, a produo aumenta com o
cubo da velocidade do vento at atingir a velocidade nominal,
rated
V . A partir desta velocidade, a
turbina funciona na zona III, de produo nominal, onde a produo de electricidade
aproximadamente constante. A transio da zona III para a zona IV d-se quando se atinge a
velocidade mxima,
out cut
V
_
, a partir da qual no ocorre produo de electricidade pois o
funcionamento neste regime seria prejudicial para a turbina
(24)
.
Apesar das VAWT tipo giromill serem o tipo mais simples de turbina elica, a sua anlise
aerodinmica bastante complexa
(3)
. Considere-se uma turbina de raio R , girando com uma
velocidade angular e , submersa num escoamento com velocidade
=5m/s
A intensidade da componente normal do escoamento nas ps
n
V no depende da velocidade de
rotao da turbina e tendo uma curva igual para uma igual velocidade do vento
V . A intensidade
da componente tangencial do escoamento
t
V altera-se com a mudana de
V e do coeficiente de
19
velocidade perifrica . O aumento de , leva ao aumento de
t
V . Como
n
V no se altera, conclui-
se que com o aumento de ,
t
V torna-se a componente dominante no escoamento, diminuindo a
amplitude do ngulo de ataque o , como se ver de seguida.
O ngulo de ataque pode ser obtido como
(3)
|
|
.
|
\
|
=
t
n
V
V
1
tan o . (8)
Substituindo os resultados anteriores, obtm-se
|
|
.
|
\
|
+
=
|
|
.
|
\
|
+
=
u
u
u
u
o
) cos(
) (
tan
) cos(
) (
tan
1 1
sen
V V
sen V
. (9)
Conclui-se que o ngulo de ataque o depende apenas da posio da p u e do coeficiente de
velocidade perifrica , sendo independente da velocidade do vento
V .
A Figura 7 mostra o ao longo de u para diferentes .
Figura 7 ngulo de ataque ao longo da posio da p , em funo do coeficiente de velocidade perifrica
1 2 3 4 5 6 7 8 9
max
o
90 30 20 14,5 11,55 9,6 8,22 7,18 6,38
max o
u 180 120 110 105 101 100 98,5 97 96,5
Quadro 5 Valores para
mx
e
max
em funo do coeficiente de velocidade perifrica
Como se viu anteriormente, a variao de o tanto maior quanto menor
(25)
. Verifica-se
tambm que, para 3 > , < 20
max
o tornando os efeitos de perda nas ps desprezveis.
A intensidade da velocidade do escoamento no perturbado relativo p dada pela expresso
2 2
n t p
V V V + =
(3)
. (10)
-90
-75
-60
-45
-30
-15
0
15
30
45
60
75
90
0 90 180 270 360
n
g
u
l
o
d
e
a
t
a
q
u
e
Posio
1
2
3
4
5
6
7
20
Substituindo as expresses anteriores obtm-se
+ + =
2 2
) ( ) cos ( u u e sen V V r V
p
+ + =
2 2
) ( ) cos ( u u sen V V V V
p
+ + + =
u u u
2 2 2 2 2 2 2
cos cos 2 sen V V V V V
p
. 1 cos 2
2
+ + =
u V V
p
(11)
A mdia da velocidade numa revoluo dada por
. 1
2
+ =
V V
p
(12)
Os valores mximos e mnimos so dados por
( ) 1 1 2
2
max
+ = + + =
V V V
p
, (13)
( ). 1 1 2
2
min
= + =
V V V
p
(14)
A variao relativa da intensidade da velocidade dada por
( ) ( )
+
+
= A
1
1 1
2
V
V V
V
p
1
2
2
+
= A
p
V
.
(15)
Portanto
p
V A
no depende de
V .
Note-se que
= =
V V
p
lim lim
0
2
lim lim = = A
p
V
.
O conhecimento dos parmetros que afectam o escoamento na p da maior importncia. Como
se viu anteriormente, um maior coeficiente de velocidade perifrica diminui a variao do ngulo
de ataque o e a variao relativa da velocidade do escoamento no perturbado relativo p
p
V A .
Este conhecimento ser da maior importncia na escolha do regime de funcionamento da turbina
pois a variabilidade de o e
p
V est relacionada com uma maior fadiga dos equipamentos
(2)
.
A partir da intensidade mdia da velocidade
p
V , possvel obter o nmero de Reynolds mdio Re
para a configurao em estudo considerando
c
V
u
1
Re
2
+
=
, (16)
sendo u a viscosidade cinemtica do escoamento (considerada
5
10 5 , 1
m
2
/s) e c a corda das
ps da turbina.
21
O Grfico 4 ilustra Re para vrias configuraes.
Grfico 4 - Nmero de Reynolds mdio para vrias configuraes
As foras aerodinmicas resultantes do escoamento nas ps
p
V podem ser obtidas pela
interpolao de grficos do coeficiente de sustentao
L
C e de resistncia
D
C de acordo com o
ngulo de ataque o e o nmero de Reynolds Re em cada posio u . No entanto, no fcil
encontrar este tipo de informao para os regimes de funcionamento tpicos de uma VAWT (baixos
Re
e elevados o ). Um dos trabalhos pioneiros nesta rea foi de Sheldahl (26). No entanto, como se
viu anteriormente, para 3 > , os valores de o so baixos o suficiente para se poder usar a maioria
das medies disponveis. Com a obteno das foras de sustentao L e de resistncia D, obtm-
se as componentes de fora normal
n
F
e tangencial
t
F
nas ps dadas por
(3)
( ) ( ) o o cos D Lsen F
t
+ = , (17)
( ) ( ) o o Dsen L F
n
+ = cos . (18)
A Figura 8 ilustra esta situao para um perfil simtrico.
A partir das relaes apresentadas conclui-se que o ngulo de ataque o preponderante na fora
tangencial
t
F produzida. A relao apresentada anteriormente de o com o coeficiente de
velocidade perifrica que define o regime de funcionamento da turbina. A baixos , o
elevado, fazendo com que a p entre em perda e diminua a relao sustentao/resistncia D L . A
elevados , o baixo e a componente tangencial de L reduzida. Existe um ptimo
intermdio onde o suficientemente elevado para que a componente tangencial de L seja
considervel e suficientemente baixo para que a relao D L seja elevada.
De acordo com as concluses anteriores, um perfil adequado para uma VAWT giromill deveria
entrar em perda a elevados o , o que permitiria a sua utilizao em ngulos elevados (baixo ),
22
maximizando a componente tangencial de L , aumentando
t
F . Por outro lado, a relao D L
deveria ser elevada para diminuir a resistncia de cada p. Existem perfis que apresentam uma
relao D L bastante elevada num dado regime de funcionamento. Esta caracterstica usada em
aeronutica no desenvolvimento de perfis de asas e o mesmo se pode aplicar no desenvolvimento de
VAWTs. Esta caracterstica denominada como drag bucket dada a curva de D em funo de L
apresentar um rebaixamento numa dada regio de funcionamento fazendo lembrar um balde (em
ingls bucket).
Figura 8 Representao de foras numa VAWT giromill para um perfil simtrico
A fora tangencial mdia dada por
}
=
t
u
t
2
0
2
1
d F F
t t
(3)
. (19)
O binrio mdio no veio da turbina dado por
r F n M
t
=
(20)
sendo n o nmero de ps e r o raio da turbina
(3)
.
A potncia mdia da turbina dada por
e M P = (21)
p
p
p
23
sendo e a velocidade angular da turbina
(3)
.
Nas ps da turbina, a rpida variao do ngulo de ataque o e da intensidade do escoamento
p
V
faz com que a resposta da p seja dinmica, e portanto diferente dos resultados normalmente
publicados em regimes estticos. Essa diferena assinalvel e depende de um grande nmero de
factores. A prpria perda (stall), passa a ser denominada perda dinmica (dynamic stall). O fenmeno
de perda dinmica caracterizado por um atraso na entrada em perda produzindo valores de
sustentao L , resistncia D
ou binrio M diferentes dos obtidos com perda esttica,
nomeadamente com histerese. Para baixos valores do nmero de Mach, caracterizado pela
formao de um vrtice no bordo de ataque e sucessiva passagem sobre a superfcie de baixa
presso da p. afectado por parmetros relacionados com o movimento da p como o tipo de
movimento, frequncia ou ngulo de ataque mximo. tambm afectado por parmetros
relacionados com a camada limite como o perfil usado ou o nmero de Reynolds Re
(2) (27) (28) (29)
. Em
geral, quanto maior o ngulo de ataque mximo, maior o crescimento do vrtice de perda dinmica
e mais acentuados os seus efeitos. Relativamente frequncia do movimento, acima de um certo
valor, ocorre um perodo duplo de variao dos parmetros, por exemplo, de L , com o ngulo de
ataque o como mostra a Figura 9.
Figura 9 - Coeficiente de sustentao C
L
em funo do ngulo de ataque para um perfil com oscilao sinusoidal para
diferentes frequncias
(27)
No caso de uma VAWT, o ngulo de ataque o varia, como se viu anteriormente, e varia tanto mais
quanto menor o coeficiente de velocidade perifrica . Desta forma, tem uma grande influncia
no regime de perda dinmica da VAWT e, embora a perda dinmica leve a uma maior eficincia da
turbina, provoca tambm maior fadiga e rudo
(25)
.
24
Figura 10 Perfil NACA0015 em condies de perda dinmica
(30)
Nobuyuki
(25)
, apresenta uma descrio do escoamento atravs de uma VAWT da qual se
apresentam alguns resultados atravs da Figura 11 e da Figura 12.
Como se pode constatar, ocorrem alteraes importantes no escoamento a montante do poste.
Estas alteraes no escoamento vo trazer importantes consequncias para o desempenho das ps a
jusante do poste. Paraschivoiu
(2)
indica que a dificuldade em simular tais alteraes no escoamento,
um grande entrave no desenvolvimento de modelos para VAWTs.
Figura 11 Visualizao experimental do escoamento numa VAWT. Re=3x10
3
, =2, imagens com 80x80mm
(25)
25
Figura 12 Descrio da evoluo dos vrtices iniciais formados no bordo de ataque (1-a) e no bordo de fuga (1-b) e
dos vrtices secundrios formados no bordo de ataque (2-a) e no bordo de fuga (2-b). Re=3x10
3 (25)
.
Larsen
(30)
, sintetiza e apresenta modelos para o clculo numrico de perda dinmica.
Tem havido uma grande procura de mtodos para o controlo de perda dinmica. Bons resultados
(e possivelmente aplicveis numa VAWT) tm sido obtidos com sistemas de sopro
(31)
ou controlo do
ngulo de ataque usando actuadores piezocermicos
(32)
. Outros dispositivos e tcnicas de controlo
so apresentados mais frente.
26
3.2 rea de varrimento
A rea A de varrimento de uma VAWT determinada pela sua altura, dimetro e forma das ps. A
obteno de dimenses ptimas tem sido uma preocupao constante em engenharia de turbinas
elicas, cuja resposta reside normalmente no estudo de factores econmicos
(2)
. Vrios estudos
foram feitos nesta rea mas essencialmente para HAWTs. Para VAWTs Darrieus, os estudos
apresentados referem-se normalmente a grandes dimenses (reas de varrimento superiores a 1000
m
2
). Estes estudos concluram que, assumindo uma estrutura e materiais similares, a massa (e
consequentemente o custo) era proporcional ao cubo do dimetro da VAWT (Darrieus) sendo que o
mesmo acontecia para o custo dos sistemas de transmisso. Tal explicvel pela menor
disponibilidade de componentes de grandes dimenses, dificuldade associada sua produo,
transporte e instalao e menor capacidade de produo em massa de componentes de grande
dimenso. No entanto, como se viu anteriormente, a energia total capturada proporcional rea
da turbina (para um alongamento 1 = AR , ser proporcional ao dimetro ao quadrado). Assim,
pode-se concluir que no existe benefcio num dimetro maior
(2)
.
s dimenses e peso da turbina esto tambm associados os custos de transporte e instalao.
Relativamente ao transporte, foram contactados transitrios que informaram que, dentro de um
tecto mximo de dimenses, os custos so proporcionais ao volume da encomenda, respectivo peso
e custo do carregamento do respectivo material a transportar. Os pesos considerados so
relativamente elevados. O volume contabilizado, no mnimo, ao m
3
. No caso do carregamento, o
ideal que o material seja transportvel por uma pessoa sem ajuda, e que caiba num veculo de
transporte convencional. Assim, para evitar custos elevados no transporte, o ideal seria uma turbina
composta por componentes que no ultrapassem individualmente os 40 kg. As dimenses mximas
deveriam ser inferiores a 1,20 m de forma a caber numa palete EUR (ISO1), cujas dimenses so 1,20
por 0,80 m.
Desta forma, para reduzir custos no transporte, conclui-se que a turbina dever ser composta por
vrias componentes, cada uma com um comprimento mximo de cerca de 1,10 m e um peso mximo
de 40 kg. Estas condies so tambm muito favorveis para a instalao da turbina.
27
3.3 Limite de Betz e perdas
Como indicado anteriormente, o limite de Betz estipula que
% 59
max
~
P
C .
No entanto, na prtica, o valor mximo obtido inferior a 45%
(33)
. A teoria de Betz assume vrias
condies: o fluido considerado homogneo e incompressvel, no existe resistncia, a turbina tem
um nmero infinito de ps e o escoamento no ganha qualquer rotao. Na prtica, a rotao
transmitida ao escoamento aps passar pela turbina, a resistncia aerodinmica e a existncia de um
nmero finito de ps e respectivas perdas levam a uma diminuio do valor de
max P
C obtido
(5)
.
Escoamento com rotao
Considerando a rotao transmitida ao escoamento pela turbina obtm-se uma expresso que
relaciona
max P
C com o coeficiente de velocidade perifrica . O Quadro 6 apresenta alguns valores
obtidos e no Grfico 3 possvel comparar a importncia do funcionamento a elevados
(33) (34)
.
0,5 1,0 1,5 2,0 2,5 5,0 7,5 10,0
max P
C 0,289 0,416 0,477 0,511 0,533 0,570 0,581 0,585
Quadro 6 Valores de C
pmax
em funo de de acordo com a teoria de Glauert-Schmitz
As end/tip losses ocorrem pela passagem de fluido da zona de maior presso para a de menor
presso no extremo de uma asa. Est tambm associado formao de vrtices nestas zonas. Este
efeito de perda atenuado pela utilizao de ps com alongamento AR superior e, para HAWTs,
considera-se
|
.
|
\
|
= q
n
tip
84 , 1
1 , (22)
sendo n o nmero de ps e o coeficiente de velocidade perifrica
(33)
. Este o principal factor
apontado para as diferenas entre as simulaes CFD em 2D e 3D
(35)
. Robert
(35)
conclui que os
vrtices nas pontas das ps tm maior intensidade nas posies em que a p tem maior lift (embora
com um ligeiro atraso devido ao tempo necessrio para o escoamento responder a essas condies).
Esta concluso pode ser verificada na Figura 13 onde a intensidade dos vrtices na ponta das asas
superior na posio a
(35)
.
Ye
(36)
conclui (embora para uma VAWT para uso hdrico) que para um alongamento 3 > AR as
end/tip losses so desprezveis, como mostra o Grfico 5.
Comparou tambm as velocidades induzidas nas ps, em vrias zonas, para uma simulao 3D com
os resultados da simulao 2D, sendo que a diferena entre estas aumenta com a proximidade
ponta da asa. Concluiu assim que as end/tip losses so o efeito tridimensional mais importante.
Ye
(36)
testou ainda a utilizao de end-plates obtendo os resultados apresentados no Grfico 6.
28
Figura 13 Contornos de vorticidade a 415 RPM, V
=1,5 m/s
(36)
29
Este resultado justifica a utilizao de end plates nas pontas das ps para 25 , 2 > de forma a
diminuir a formao de vrtices, tcnica j utilizada na aviao pela adio de winglets nas asas.
Alm deste aspecto, este resultado permite avaliar a importncia das end/tip losses. Ao usar end-
plates, Ye
(36)
considera que a boa aproximao entre os modelos numricos 2D e os resultados
experimentais permitem concluir que este o factor predominante de efeitos tridimensionais.
Existem ainda perdas relacionadas com a interaco entre as ps que so muito difceis de avaliar.
Robert
(35)
, tambm conclui que, sob certas condies, as ps podem vir a sofrer influncia da sua
prpria esteira como mostra a Figura 14 onde a seco W representa a esteira da p c. Os mesmos
resultados so descritos por Nobuyuki
(25)
como se viu anteriormente na Figura 12.
Figura 14 Em certas condies de rotao, as ps esto sujeitas sua prpria esteira (W da p c)
(35)
As perdas em outras componentes da turbina so maioritariamente atribudas aos rolamentos e
aos braos de suporte
(35)
.
30
3.4 Alongamento
O alongamento (em ingls aspect ratio) AR definido como
D
h
AR = , (23)
sendo h a altura da turbina e D o dimetro da mesma
(2)
. Por vezes, AR definido em funo do
raio ao invs do dimetro.
Inicialmente, as VAWTs tinham um baixo AR de forma a minimizar o comprimento das ps e, em
geral, a quantidade de material usado para uma dada rea
(2)
. No entanto, existem importantes
perdas nas extremidades das ps e aumentando AR , diminui-se o valor relativo destas perdas.
Como se viu anteriormente (Grfico 5), para um 3 > AR , estas perdas so desprezveis.
Outro factor importante a variao da intensidade do vento com a altitude. Para uma turbina
suficientemente grande este um efeito importante. Paraschivoiu
(2)
apresenta um modelo terico
para o clculo da influncia do alongamento no funcionamento da turbina (devido variao da
intensidade do vento com a altitude) e conclui que a diferena na potncia produzida pode ser
superior a 10%, mas sem comparar com resultados experimentais.
O aumento do AR leva a um aumento da velocidade angular e (para manter a mesma
velocidade de escoamento nas ps
p
V e consequentemente o mesmo coeficiente de velocidade
perifrica ) e a uma diminuio do binrio M para a mesma potncia
(2)
. Como se ver mais
frente, a escolha do gerador determina a velocidade de funcionamento. Como os geradores
necessitam, em geral, de velocidades elevadas, pode-se optar por um sistema de baixa velocidade e
elevado binrio (baixo AR ), usando uma caixa de velocidades ou um sistema de elevada velocidade
(elevado AR ) sem caixa de velocidades.
A escolha deste parmetro vai depender essencialmente do regime de funcionamento desejado
(dependente essencialmente do tipo de gerador escolhido) e de um estudo do custo associado ao
mesmo.
31
3.5 Coeficiente de velocidade perifrica
O coeficiente de velocidade perifrica , considerado um factor crucial no desenvolvimento de
uma turbina elica
(33) (34)
. necessrio saber o ptimo para se avaliar as condies de
funcionamento, potncia obtida e o retorno do investimento efectuado
(33) (34)
. Para um sistema
HAWT,
optimo
normalmente entre 7 e 8 para uma turbina com 3 ps
(33) (34)
. Para um sistema VAWT,
optimo
normalmente entre 3 e 5
(35)
.
Paraschivoiu
(2)
refere que a curva caracterstica do coeficiente de potncia
p
C em funo de
pode ser explicada considerando 3 regimes de funcionamento. Para valores de abaixo do valor de
optimo
(onde
p
C mximo), ocorre perda nas ps. Para valores de acima do valor de
optimo
,
predominam efeitos secundrios (embora no indique quais). Para valores perto de
optimo
, existe
uma zona de transio onde todos os fenmenos esto presentes mas em menor escala. A baixos ,
as ps no interagem suficientemente com o escoamento enquanto a altos , as ps comeam a
sofrer do escoamento de outras ps, como mostra a Figura 12
(35)
.
Como os geradores tm melhores eficincias a velocidades mais elevadas, um mais elevado
ajuda na eficincia obtida no gerador
(2)
. Associado a mais elevado, esto maiores valores de rudo
(5)
.
32
3.6 Factor de bloqueamento e nmero de ps
O factor de bloqueamento (do ingls solidity), uma caracterstica crucial no desenvolvimento de
uma turbina e considera-se
r nc = o , (24)
sendo n o nmero de ps, c a corda de cada p e r o raio da turbina
(36)
.
Alguns resultados indicam que um menor factor de bloqueamento o resulta numa maior
eficincia
P
C e num coeficiente de velocidade perifrica ptimo
optimo
superior
(35)
. Paraschivoiu
(2)
conclui que
optimo
diminui com o aumento de o e que o coeficiente de potncia mximo
max P
C ,
para uma dada velocidade do vento
associado
(2)
. Tradicionalmente, usaram-se perfis simtricos em VAWTs pois como o ngulo de ataque
o muda de sinal a cada 180, pensava-se que usando perfis com curvatura, quaisquer ganhos
obtidos com o positivo seriam muito inferiores eficincia perdida com o negativo. No entanto,
resultados recentes como os obtidos por Kirke
(1)
, mostram o contrrio.
37
3.9 Gerador
As VAWTs de grandes dimenses usam, normalmente, uma caixa de velocidades e um gerador de
induo. A caixa de velocidades e o sistema de transmisso representam normalmente 50% do custo
de uma VAWT de grandes dimenses. Este elevado custo tem levado tentativa de desenvolvimento
de geradores para funcionamento a baixas rotaes, que permitam uma acoplagem directa
turbina, sistemas estes conhecidos como direct drive
(2)
. Note-se que existem perdas importantes
com a utilizao de uma caixa de velocidades que podem atingir os 6%
(39)
. Normalmente, os
geradores necessitam de elevadas velocidades de rotao e para um funcionamento eficiente.
Como se viu anteriormente, possvel controlar os parmetros de funcionamento da turbina, por
exemplo, atravs do alongamento AR . A soluo para os geradores convencionais pode ser uma
turbina de elevada velocidade e (baixo factor de bloqueamento o , elevado AR ) ou uma turbina
de baixa velocidade e acoplada a uma caixa de velocidades. No caso da turbina Turby o gerador
usado do tipo manes permanentes direct-drive devido sua simplicidade e robustez
(40)
.
No caso das VAWTs, a colocao do gerador no solo ou a uma altura mais reduzida permite uma
maior liberdade no peso e dimenses do gerador
(2)
.
A escolha de um gerador est profundamente ligada ao tipo de ligao que se pretende. VAWTs
podem ser ligadas a bancos de baterias e nesse caso pretende-se um gerador DC ou AC com
posterior rectificao. No caso de sistemas para ligao rede, pretende-se condies mais
especficas, nomeadamente de frequncia e tenso. Recentemente, o desenvolvimento da
electrnica de potncia tem permitido uma grande flexibilidade dos sistemas a usar
(2) (21)
.
38
3.10 Arranque
Uma das grandes desvantagens das VAWTs tipo lift a dificuldade para o arranque automtico das
mesmas. No caso de VAWTs ligadas rede elctrica, possvel efectuar o seu arranque funcionando
o gerador como motor ou usando um motor prprio para o efeito. No caso de turbinas de pequenas
dimenses, a incluso de tal caracterstica dispendiosa. No caso de VAWTs que no estejam ligadas
rede elctrica esta caracterstica torna-se uma grande desvantagem
(1)
. Tm sido propostas vrias
solues para facilitar o arranque de VAWTs. Usualmente, as turbinas de menores dimenses
incorporam elementos de turbinas do tipo drag como o caso da VAWT Ropatec ilustrada em
anexo. No entanto, Kirke
(1)
, demonstra que a incorporao destes elementos com reas inferiores a
10% da rea frontal da VAWT, tm pouco efeito. As VAWT Windsmile usam um perfil patenteado, J-
blade (Figura 21), que de acordo com a companhia, permite um compromisso entre VAWTs tipo lift
e drag.
Figura 21 Perfil J-blade
Kirke
(1)
indica que as solues actualmente com maior viabilidade so a utilizao de ps com
curvatura, ps flexveis ou ps que permitam a variao automtica do ngulo de ataque de forma
passiva (passive pitch control). A utilizao de mecanismos para o controlo activo do ngulo de
ataque das ps (active pitch control) revela-se demasiado difcil e dispendiosa para VAWTs de
pequenas dimenses
(1)
.
39
3.11 Outros dispositivos
A presena de uma torre central afecta o escoamento ao longo da turbina. Considerando que a
esteira da turbina tem uma largura L ,
L
D V C
V
p poste Dp
esteira
4
= A (25)
sendo
esteira
V A a variao da velocidade do escoamento na zona considerada da esteira,
Dp
C o
coeficiente de atrito do poste,
poste
V a velocidade do escoamento no poste e
p
D o dimetro do
poste
(2)
. White
(6)
indica que o valor de
Dp
C aproximadamente 1,2 em regime laminar e 0,3 em
regime turbulento.
Embora as ps sejam as principais responsveis pelas foras de atrito na turbina, a torre e os
suportes ajudam a reduzir o coeficiente de potncia mximo
max P
C . Os suportes so necessrios
para fornecer rigidez s ps e diminuir as tenses na turbina.
A variao de
P
C pode ser calculada como
b
t
D p
S
S
C C
3
0
4
1
= A , (26)
sendo
0 D
C o coeficiente de atrito do brao usado, o coeficiente de velocidade perifrica,
t
S a
rea varrida pela turbina e
b
S a rea total dos braos
(2)
.
Tem havido uma activa investigao no desenvolvimento de mtodos para controlo da potncia
produzida pela turbina. No caso de HAWTs, esse controlo feito pela actuao em quatro
parmetros: o ngulo de ataque das ps o (normalmente denominado pitch control), a velocidade
do rotor, a dimenso da turbina e as caractersticas aerodinmicas dos perfis usados. Actualmente,
as tcnicas mais utilizadas referem-se ao pitch control, usando um controlo activo do ngulo de
incidncia das ps e pelo desenvolvimento de perfis
(5)
. Kirke
(1)
, apresenta um estudo exaustivo de
sistemas de controlo de ngulo de ataque activos (active pitch control) e passivos (passive pitch
control) e conclui que para VAWTs de pequenas dimenses, os sistemas passivos trazem grandes
vantagens. Um destes resultados apresentado na Figura 22.
40
Figura 22 C
P
em funo de numa VAWT com variao sinusoidal de para vrias amplitudes
(1)
Outra tcnica a utilizao de velocidade de funcionamento varivel, da qual est dependente o
funcionamento do gerador. A alterao da dimenso do rotor seria uma tcnica de grande utilidade
permitindo o aproveitamento de ventos a baixas velocidades e a sobrevivncia com ventos de
grande velocidade, reduzindo as dimenses das ps.
Por fim existem tcnicas de controlo activo de escoamento (active flow control) que j so h
muito estudadas na indstria aeronutica e que tm grande aplicabilidade em turbinas elicas
(24)
.
A entrada em perda das ps pode ser alterada pela aplicao de geradores de vrtices que so
pequenas alhetas que se comportam como asas finitas de cujos bordos marginais emanam vrtices.
Esses vrtices promovem uma transio da camada limite de laminar a turbulento, aumentam a
quantidade de movimento da mesma e atrasam, ou evitam, a entrada em perda da p
(2) (37)
. Vrios
estudos foram efectuados nesta rea mas a grande variedade de parmetros em jogo ainda no
permite uma concluso final sobre a efectividade da utilizao destes sistemas
(2)
. No entanto,
estudos apontam para um possvel aumento da produo anual da ordem dos 10% e para a
diminuio da sensibilidade da turbina acumulao de sujidade ou gelo nas ps. A colocao dos
geradores de vrtices deve ser feita a 30% da corda e de ambos os lados da p
(2)
.
Uma descrio exaustiva dos mtodos de controlo em desenvolvimento ou j correntemente
usados em turbinas (essencialmente HAWTs) feita por Johnson
(24)
.
41
4 Sistemas similares
Com base na referncia (41), obtiveram-se dados de vrias VAWTs comercialmente disponveis.
Estes dados foram compilados numa base de dados disponvel no Anexo 1. A obteno desta base de
dados permite a avaliao de importantes parmetros no desenvolvimento de VAWTs.
As turbinas elicas tm normalmente presente uma indicao da produo nominal, que ocorre
velocidade nominal,
rated
V . Esta velocidade normalmente elevada (10-16 m/s) levando a que os
valores de potncia apresentados no sejam fiveis para a estimao de produo energtica de uma
turbina. O mtodo mais correcto para a estimao de produo energtica passa pela apreciao da
curva de potncia e da distribuio da probabilidade da velocidade do vento no local da instalao.
Alm deste aspecto, as curvas de potncia fornecidas pelos diferentes construtores no tm em
conta as mesmas condies. Alguns modelos trazem medidas da energia produzida em funo da
velocidade mdia do vento. Estes resultados diferem na distribuio usada (note-se que mesmo com
a mesma mdia, os resultados so diferentes para distribuies diferentes, especialmente tendo em
conta a dependncia cbica da potncia do vento
elica
P com a velocidade do mesmo
V ). Assim,
para uma melhor apreciao dos modelos existentes, obteve-se para cada um a produo energtica
estimada num local com
V mdio de 5 m/s.
A produo energtica mdia foi convertida para potncia, sendo apresentada nos resultados
seguintes como potncia mdia a 5m/s. Comeou-se por comparar a potncia nominal com a
potncia mdia a 5m/s, como apresentado no Grfico 7. Como se verifica, existe uma grande
diferena entre a potncia nominal e a energia realmente produzida pela turbina para s m V / 5 =
.
Como se viu anteriormente, importante que a turbina desenvolvida tenha um preo competitivo.
Assim, analisou-se o preo das diferentes turbinas e concluiu-se que tal comparao difcil pois os
preos apresentados dependem das circunstncias. Estas circunstncias so maioritariamente a
incluso, ou no, de impostos, do custo de transporte, montagem, garantias/extras e o local.
42
Grfico 7 Potncia mdia a 5m/s vs potncia nominal em VAWTs de pequena dimenso
Para se ter uma ideia da dificuldade associada a uma justa comparao de preos, veja-se o
exemplo do Quadro 7, onde se comparam os custos da mesma turbina (Windspire standard 1,2kW)
indicados por diferentes fornecedores.
Fornecedor Custo
Stony Plain, Canada (centro-oeste)
6523 Preo no local com impostos;
767 - Base de beto;
<11500 - Sistema completo com montagem;
Dugald, Canada (centro sul)
385 - Autorizao de engenharia;
14000 - Instalao completa sem transporte ou
impostos;
Costa Rica (Amrica Central) 8950 - Instalao completa com impostos.
Dublin, Irlanda 14541 - Instalao completa com impostos.
Rennes, Frana (Nordoeste) 13500 - Instalao completa sem impostos.
Quadro 7 Variao de preos na mesma turbina
Como se pode ver a variao de preos muito grande. As turbinas so produzidas em Michigan,
nos EUA e as diferenas de preos entre pases seria de esperar, tal como acontece com a maioria
dos produtos. ainda mais interessante verificar as diferenas no mesmo pas.
Mesmo considerando as dificuldades em obter um custo nas mesmas condies, tentou-se obter o
preo da turbina completa no local de produo (sem montagem ou transporte) em funo da
potncia mdia a 5m/s como mostra o Grfico 8.
Grfico 8 Preo em funo da potncia real em VAWTs de pequena dimenso
y = 0,07x + 24
0
100
200
300
400
500
600
0 1000 2000 3000 4000 5000 6000 7000
P
o
t
n
c
i
a
m
d
i
a
a
5
m
/
s
W
Potncia nominal W
y = 16x + 5292
0
5000
10000
15000
20000
25000
0 100 200 300 400 500 600
P
r
e
Potncia mdia a 5m/s W
43
Note-se a disparidade de preos existente. Certamente que o mtodo de avaliao usado e a
dificuldade em obter preos nas mesmas condies para todas as VAWT pesam nas diferenas
encontradas mas, tal como todos os produtos, existem vrias gamas de qualidade disponveis. O
facto de este mercado estar ainda em crescimento, facilita tambm a existncia de modelos e
fornecedores que dominam pequenos nichos onde a concorrncia ainda no existe. De qualquer
forma, obtm-se uma base de comparao, que ser um guia no desenvolvimento desta turbina.
Para a situao presente, voltou-se a utilizar a base de dados de VAWTs de pequenas dimenses
para se obter uma relao entre o custo da turbina e a respectiva rea, apresentado no Grfico 9.
Grfico 9 Preo em funo da rea para VAWTs de pequena dimenso
A baixa abrangncia de reas consideradas no permite uma avaliao completa da evoluo do
preo com a rea. No entanto, para VAWTs de pequenas dimenses, os custos parecem ser
proporcionais rea. Note-se que os preos da base de dados no consideram os custos de
transporte e instalao.
Outra comparao importante a da potncia mdia a 5m/s em funo da rea A da turbina,
apresentada no Grfico 10.
Grfico 10 Potncia mdia a 5m/s em funo da rea A para VAWTs de pequena dimenso
y = 1295x + 1918
0
5000
10000
15000
20000
25000
30000
35000
40000
0 5 10 15 20
P
r
e
Area m
y = 25x + 42
0
100
200
300
400
500
600
0 5 10 15 20
P
o
t
n
c
i
a
m
d
i
a
a
5
m
/
s
W
rea m
Mximo
para 5 m/s
44
Novamente, existe uma grande disparidade entre os valores recolhidos, pelas mesmas razes
apresentadas anteriormente.
O Grfico 11 apresenta a altura em funo do dimetro para as turbinas analisadas.
Grfico 11 Alongamento AR de VAWTs de pequena dimenso
O grfico apresenta resultados para 42 turbinas, curvas de rea (a cheio) e de alongamento AR (a
tracejado) de acordo com o dimetro e altura. Embora as VAWTs com 1 > AR sejam cerca de dois
teros do total, existem muitos exemplares com baixo AR . Esta disparidade tem a ver,
essencialmente, com o tipo de gerador, como se viu anteriormente.
Das VAWTs de pequenas dimenses analisadas, a maioria apresenta um gerador direct drive de
manes permanentes. A maioria produzida pelos prprios fabricantes. H, no entanto, geradores
deste gnero, especialmente desenvolvidos para turbinas elicas j disponveis no mercado.
O nmero de ps varia entre 3 a 5, sendo 5 o valor mais comum. No esto disponveis
informaes sobre parmetros como o factor de bloqueamento ou os tipos de perfis usados.
Os dados recolhidos permitem uma rpida avaliao de vrios parmetros. Por exemplo, uma
turbina com 1 a 5 m
2
dever ter uma potncia mdia de 70 a 170 W para s m V / 5 =
. O preo
dever ser entre 6000 a 8000.
0
1
2
3
4
5
6
7
8
0 1 2 3 4 5 6 7 8
A
l
t
u
r
a
m
Dimetro m
AR=2
AR=1
AR=0,5
AR=3
AR=0,33
A=2m
2
A=6m
2
A=16m
2
45
5 Anlise aerodinmica
5.1 Estado da arte
Os modelos aerodinmicos so essenciais para a deduo de parmetros ptimos de
funcionamento e para a previso da eficincia da turbina antes da sua produo. Actualmente, no
existe um procedimento terico completo para o desenvolvimento de uma turbina elica, residindo
grande parte do desenvolvimento na obteno de resultados de vrias configuraes por clculo
numrico
(17)
. Vrios modelos foram desenvolvidos para VAWTs, sendo os mais usados actualmente o
double/multiple streamtube model, o vortex model e o cascade model
(3)
. As principais componentes
destes modelos, por habitual ordem de desenvolvimento, so
- Clculo da velocidade
p
V em funo da posio u de acordo com o coeficiente de
velocidade perifrica , raio da turbina r e velocidade de escoamento livre
V ;
- Clculo das foras aerodinmicas na p;
- Clculo da velocidade induzida
i
V
(velocidade do escoamento no interior da turbina,
diferente da velocidade do escoamento livre
=
AV
F n r
C
t
P
e
(28)
}
=
t
u
t
2
0
2
1
d F F
t t
(29)
( ) ( ) ( )
}
=
t
u o o
t
e
2
0
3
cos
2
1
2
1
d D Lsen
AV
n r
C
P
(30)
Trabalhando em termos de coeficientes,
2
2
1
p
L
cV
L
C
= , (31)
2
2
1
p
D
cV
D
C
= , (32)
o
o D D D
C C C + =
0
, (33)
o
o L L
C C = . (34)
Obtm-se
( )
3 0 2 1
2
1
+ + + =
d d l P
C C C
AV
nc r
C
o o
t
e
, (35)
sendo
( )
( )
( )
( )
( ) ( ) u u
u
u
u
u
t
d
sen
tg sen
sen
tg 1 cos 2
cos cos
2 1
2
0
1
1
+ +
|
|
.
|
\
|
|
|
.
|
\
|
+
|
|
.
|
\
|
+
= +
}
, (36)
( )
( )
( )
( )
( ) ( ) u u
u
u
u
u
t
d
sen
tg
sen
tg 1 cos 2
cos
cos
cos
2 1
2
0
1
2
+ +
|
|
.
|
\
|
|
|
.
|
\
|
+
|
|
.
|
\
|
+
= +
}
, (37)
( )
( )
( ) ( )
}
+ +
|
|
.
|
\
|
|
|
.
|
\
|
+
= +
t
u u
u
u
2
0
2 1
3
1 cos 2
cos
cos d
sen
tg . (38)
Os integrais foram calculados usando uma calculadora grfica Ti-89 para vrios valores de . Com
os resultados dos integrais, o coeficiente de potncia dado pela frmula seguinte,
52
( ) ( ) ( )
0
2
1 2 44 , 1
D D L P
C C C
A
nc
C + + + = t t
o o
(39)
Em concluso, este modelo relativamente simples de implementar, bastando saber para a
condio de clculo que se pretende ( Re ) quais as caractersticas de
o L
C ,
o D
C e
0 D
C dos perfis
usados. O modelo no tem em conta a interferncia dos diferentes elementos entre si nem
caractersticas dinmicas pois os parmetros de aerodinmica das ps so obtidos atravs de curvas
estticas. Outra caracterstica da maior importncia a utilizao da relao linear de
L
C e
D
C com
o . Tomando como exemplo o perfil NACA0012, a partir do Grfico 13 e do Grfico 14 compreende-
se que esta relao s permanea vlida at 6 ~ o
para
4
10 4 Re = ou 10 ~ o para
5
10 6 , 3 Re = . De acordo com o Quadro 5, isto corresponde a 5 , 9 >
para
4
10 4 Re =
e
5 , 5 > para
5
10 6 , 3 Re = . Estas limitaes pem em causa a utilizao deste modelo
principalmente para VAWTs com ps com corda c reduzida e com baixos
de funcionamento. No
entanto, dada a simplicidade do mesmo, fez-se aqui a sua apresentao. Incluindo elementos que
permitam uma avaliao do efeito da turbina sobre o escoamento, por exemplo, atravs de clculos
de momento como acontece com o streamtube model, pode haver a possibilidade de
desenvolvimento de um modelo expedito e com bons resultados.
5.2.3 CFD Computational Fluid Dynamics
Aps a compreenso inicial dos mecanismos de funcionamento bsicos de VAWTs, iniciou-se a
simulao 2D para vrias configuraes em CFD. Pretendia-se fazer a simulao de vrias
configuraes que permitissem estudar vrios parmetros como os perfis usados, o nmero de ps e
o factor de bloqueamento. Pretendeu-se tambm simular configuraes que seriam idnticas s que
seriam testadas no prottipo construdo para o efeito. Como se ver mais frente, o prottipo
construdo de pequenas dimenses para reduzir os custos e facilitar a sua montagem e teste. Assim
definiu-se um raio r da turbina de 0,5m. As dimenses das ps estavam limitadas pelo processo de
construo usado e pela necessidade de reduo de custos e tempo de execuo pelo que optou-se
por definir a corda das ps c de 5 ou 10 cm. O nmero de ps n varia entre 2 e 5 dados ser esta a
gama de valores mais comum na indstria. Usaram-se os perfis NACA0012 e NACA0018 pela
disponibilidade de informao existente e pela recomendao de investigadores do IST. A lista de
configuraes simulada apresentada no Quadro 8. Para cada configurao, obteve-se o binrio M
e as foras no eixo paralelo ao escoamento
x
F e no eixo perpendicular ao escoamento
y
F .
Optou-se pela utilizao do modelo kw SST, transitional devido aos bons resultados obtidos
anteriormente por investigadores do IST e por Wang
(29)
.
53
O modelo foi implementado usando o programa Fluent. As malhas foram obtidas com o programa
Gambit. As malhas foram desenvolvidas tendo em conta vrios parmetros para assegurar a sua
qualidade usando um baixo equiangle skew, baixas variaes de tamanho de clulas adjacentes
(smoothness) e baixo alongamento das clulas. Junto s superfcies das ps e do poste central,
implementou-se uma camada fina de clulas de forma a reproduzir eficientemente as condies de
parede. Nas outras regies usou-se uma malha mais aberta de forma a diminuir os tempos de clculo
necessrios. A Figura 25 e a Figura 26 apresentam as caractersticas das malhas usadas. Todos os
clculos foram efectuados usando um computador porttil pessoal. O tempo dispendido para o
clculo e preparao de todas as condies testadas foi superior a 6 meses. No Anexo 3 apresentam-
se as condies de malha utilizadas e os tempos de clculo para cada configurao simulada.
Figura 25 Malha usada para CFD (esquerda) e detalhe da zona central para a configurao 5 (direita)
Figura 26 Detalhe da malha em torno dos perfis para a configurao 5 (esquerda) e em torno do poste central
(direita)
Iniciou-se os clculos para uma velocidade do vento s m V / 5 =
a um coeficiente de velocidade
perifrica 5 = . A partir desse ponto, em geral, foi-se alterando em 1 at o coeficiente de
potncia
P
C atingir aproximadamente 0. De seguida refinaram-se as condies para os valores
54
mximos de
P
C . Como se ver adiante, aps os resultados experimentais, foram ainda obtidas
novas condies, principalmente a mais baixos. De seguida obtiveram-se os mesmos resultados
para valores de
V , diferentes de 5
m/s, foram de 8m/s e 11 m/s. Em alguns casos apenas se obtiveram resultados para
correspondente ao
max P
C para s m V / 5 =
, 4 = .
Para cada condio foi calculado um tempo suficiente para que a turbina realizasse 4 rotaes. Os
resultados eram obtidos a partir da quarta rotao. Tal necessrio pois os clculos iniciam-se com
uma condio inicial diferente do sistema estabilizado. Assim necessrio deixar correr um tempo
de clculo suficiente para que os transientes iniciais ocorram e os resultados estabilizem.
Os resultados obtidos apontam para uma correcta simulao dos fenmenos fsicos. Ocorre a
formao alternada de vrtices a jusante do poste central. Existe tambm, como esperado, um
dfice da velocidade do escoamento
=4m/s, =4
=2
=4 =6
=8 =10
=12 =14
=16
x
y
=18
56
Figura 28 Magnitude da velocidade do escoamento para a configurao 5, V
V m/s e 8 =
V diferentes, ou seja, a
elevados as curvas traduzem uma forma sinusoidal com intensidade varivel mas f
aproximadamente constante. O Grfico 19 apresenta a evoluo de
n
F para as restantes
configuraes.
Grfico 18 - |Fn| em funo de para configurao 5 com V
=8 m/s,
optimo
Como se pode verificar, o comportamento sinusoidal verifica-se para as configuraes com 3
(configuraes 5 e 7)e 4 ps (configurao 3). As restantes configuraes apresentam uma curva
aproximadamente sinusoidal mas com distores. Nestes casos,
n
F tem mais que uma frequncia,
sendo mais difcil a anlise do comportamento dinmico da turbina.
0
10
20
30
40
50
60
70
0,00 1,05 2,09 3,14 4,19 5,23 6,28
F
n
N
rad
Conf.2, =4,5
Conf.3, =4
Conf.4, =3,5
Conf.5, =3
Conf.6, =3,5
Conf.7, =5
60
O Grfico 20 apresenta a intensidade da fora normal
n
F , em funo de de acordo com os
valores apresentados no Grfico 18.
Embora estes resultados pudessem ser avaliados a partir do Grfico 18, desta forma compreende-
se mais facilmente a evoluo de
n
F com . Note-se a amplitude de valores obtidos: Para 5 = ,
a turbina com 1 = A m
2
sujeita a esforos de intensidade entre cerca de 30N e 90N. A obteno
destes resultados da maior importncia para o dimensionamento de turbinas.
Grfico 20 |Fn| em funo de para configurao 5 com V
=8 m/s
Em resumo, obtiveram-se as curvas de funcionamento de vrios parmetros, para diversas
configuraes. Comparando os resultados obtidos entre si foi possvel verificar algumas
caractersticas importantes das VAWT, algumas j descritas anteriormente e aqui verificadas.
Nomeadamente, a existncia de um
optimo
o e a relao de o com . Um dos resultados mais
interessantes ser a comparao dos perfis NACA0018 e NACA0012. De acordo com a literatura
existente, o perfil NACA0018 deveria apresentar melhores resultados mas verificou-se o contrrio.
da maior importncia a validao dos resultados aqui obtidos atravs de testes prticos para que
sirvam de base de trabalho para outros estudos.
-0,5
-0,4
-0,3
-0,2
-0,1
0
0,1
0,2
0,3
0,4
0,5
0
10
20
30
40
50
60
70
80
90
100
0 1 2 3 4 5
C
p
|
F
n
|
N
Mximo
aximo
Mdia
aximo
Mnimo
aximo
61
6 Construo do prottipo VAWT
6.1 Reviso de materiais e processos de construo
Embora o objectivo seja apenas a produo de um prottipo para o teste de diversas configuraes
de funcionamento, fez-se uma ligeira reviso sobre o desenvolvimento destes sistemas.
Paraschivoiu
(2)
refere que uma grande parte dos problemas de fadiga apresentados em VAWTs
devem-se a conexes mal desenvolvidas e que a utilizao conjunta de sistemas de aperto de maior
rea e de adesivos de enchimento resultaram em grandes melhorias. Indica tambm que os
desenvolvimentos mais promissores na rea das VAWTs so a utilizao de cabos de apoio (em ingls
guy cables), colunas treliadas (em ingls truss columns) e materiais compsitos.
Relativamente s ps, actualmente, so usados essencialmente materiais compsitos com fibras de
vidro ou carbono em resinas de polyester ou epoxy. O interior das ps preenchido com materiais
sandwich como PVC, PET ou madeira de balsa e o exterior com materiais poliuretnicos. Nestas
condies, as ps tm uma vida til de 20-25 anos
(43)
resistindo a mais de
8
10 ciclos de
extenso/compresso no caso das HAWTs e a condies de temperaturas extremas, humidade,
chuva, neve, gelo, radiao solar, relmpagos e salinidade
(44)
. As ps em material compsito tm
uma estrutura exterior que define a sua forma, uma estrutura central que a principal responsvel
pela resistncia da p e o restante espao preenchido por compsito sandwich cujo propsito
ajudar na manuteno da forma desejada e transferir os esforos para a estrutura central
(45)
.
Alguns factores so de importante considerao no desenvolvimento das ps de uma turbina
elica. Os relmpagos so uma preocupao, principalmente no caso de turbinas que se destaquem
em altura no ambiente circundante. Normalmente este problema resolvido colocando condutores
na turbina desde o topo at base, onde se espera que a energia do relmpago se dissipe. No
entanto, continua a existir a possibilidade de destruio da turbina, principalmente pelo facto de a
sbita passagem de corrente por estes condutores transformar qualquer humidade em vapor
(46)
.
O tratamento da turbina no fim de vida tambm da maior importncia porque pode afectar os
custos da turbina e um parmetro tido cada vez mais em conta pelos consumidores,
62
principalmente quando se trata de equipamentos associados crescente conscincia ambiental,
como so as energias alternativas. Os melhores processos, por ordem decrescente, so a reduo de
produo de resduos, o que se consegue aumentando a vida til das componentes e reduzindo a
quantidade de material usado. De seguida o melhor processo a reutilizao. Em terceiro lugar tem-
se a reciclagem, que ainda um processo difcil mas em grande desenvolvimento para materiais
compsitos (espera-se que em 2030, perto de 225000 toneladas de ps de turbinas necessitem de
ser tratadas anualmente). Outras solues consistem na deposio em aterro sanitrio, uma
alternativa cada vez com menor possibilidade de uso dada a tendncia da maioria dos pases
industrializados para a reduo desta alternativa. Outra alternativa a incinerao, que a opo
mais usada actualmente, mas que apresenta problemas como a libertao de poluentes, a produo
de cinzas que necessitam de ser recicladas ou incorporadas em aterro sanitrio. Antes da
incinerao, existe ainda a necessidade de reduzir as dimenses das ps, um processo difcil e de
elevado custo, levando ao consumo de energia e libertao de poluentes
(43)
.
Assim, a possibilidade de reciclagem importante, sendo que actualmente apenas cerca de 30% de
um compsito plstico reforado com fibra reciclvel como novo material compsito. O restante
pode ser usado na indstria de construo como material de enchimento
(43)
. A melhor opo, no
caso de turbinas de pequena dimenso, evitar materiais plsticos reforados com fibras usando, no
seu lugar, plsticos facilmente reciclveis (como PET) e uniformizando os materiais usados na turbina
e facilitando a sua separao
(43)
. Outra possibilidade a utilizao de materiais alternativos como
bambu que actualmente desfeito e combinado com resinas sintticas
(43)
. O bambu um material
de grande resistncia e durabilidade. O seu rpido desenvolvimento e baixo impacto ambiental
tornam-no uma alternativa economicamente atraente e ecologicamente responsvel.
No desenvolvimento de uma turbina importa tentar reduzir os custos mesmo que tal implique uma
diminuio da eficincia de funcionamento, desde que se obtenha um produto final equilibrado. De
acordo com os custos associados a equipamentos de produo de energia, podemos adaptar para o
caso das VAWT de pequenas dimenses, obtendo-se os seguintes custos:
- Equipamento
o Turbina;
o Gerador;
o Electrnica e ligaes;
o Suporte;
- Instalao;
- Transporte;
- Manuteno;
- Desmantelamento;
- Financiamento;
63
- Licenciamento;
- Impostos ou apoios
(2)
.
Ao nvel do equipamento, importa escolher uma configurao que permita uma fcil construo.
Deve-se tambm manter ao mnimo o nmero de peas diferentes para permitir uma economia de
escala. Uma forma de reduo de custos, no caso da utilizao de materiais compsitos, poder ser
feita pela combinao de fibras de vidro e de carbono
(47)
.
Voltando aos custos totais do equipamento a utilizar, temos o custo do gerador, controlvel
essencialmente pelo tipo de gerador a utilizar. A tendncia actual, principalmente no caso de VAWTs
de pequenas dimenses, a utilizao de geradores do tipo direct drive e de manes permanentes
pelo seu menor custo, menor complexidade de electrnica associada e menor manuteno. A nica
forma de diminuir mais o custo desta componente pela compra em grandes quantidades ou pelo
desenvolvimento de um sistema prprio para o efeito. O custo da electrnica est associado
converso das caractersticas da energia, que tambm afectado, tal como o gerador, pela economia
de escala ou pelo desenvolvimento de um sistema prprio para o efeito.
Os custos de instalao devero ser reduzidos evitando a necessidade de equipamentos prprios
como gruas. O equipamento dever ser adequado para que a instalao implique apenas a
preparao do terreno com a construo de uma base em beto e a colocao de forma manual ou
utilizando um automvel puxando um cabo ligado ao equipamento. A facilidade de montagem dos
equipamentos elctricos tambm da maior importncia reduzindo-se a necessidade de
ferramentas especializadas e diminuindo-se o tempo de trabalho de pessoal especializado.
Inicialmente, o equipamento dever estar de acordo com as normas de instalao para
equipamentos da mesma categoria. Ainda no existe um licenciamento uniformizado para turbinas
de pequenas dimenses
(48)
.
A capacidade de transporte do equipamento importante. Para isso necessrio analisar todas as
dimenses de transporte, desde a possibilidade de transporte martimo em paletes estandardizadas
ao transporte manual de componentes no local de instalao. Estas questes levam-nos para um
desenvolvimento modular atravs de peas pequenas e leves cujo encaixe leve ao equipamento final.
A montagem no local seria uma grande vantagem para o transporte mas um inconveniente caso a
montagem seja complicada e/ou implique pessoal e/ou equipamento especializado.
A manuteno dever ser mnima. Embora no seja fcil encontrar referncias s principais falhas
nestes equipamentos, os prottipos montados pelos laboratrios Sandia eram normalmente
afectados por problemas nos rolamentos
(2)
.
O desmantelamento de uma turbina destas dimenses em princpio no traz grandes problemas
mas um sistema modular e a utilizao de materiais reciclveis de baixo impacto ecolgico seriam
sempre preferveis.
64
A nvel de financiamento, o baixo valor destes equipamentos no permite grandes alteraes. A
nica alterao possvel ser a alterao do custo inicial pela produo de um equipamento com
menor vida til ou menor eficincia. Outra forma de diminuio do custo pelo apoio do estado
atravs, por exemplo, da diminuio de impostos. Isso implica normalmente a certificao do
equipamento de acordo com vrias normas o que diminuiria tambm a necessidade de
licenciamento do equipamento tornando mais barato a sua instalao.
O factor combustvel no importante nestes sistemas excepto pela utilizao de electricidade
para provocar o arranque das turbinas, quando aplicvel.
Os parmetros de funcionamento como a velocidade de arranque
in cut
V
_
e a velocidade mxima
out cut
V
_
vo determinar os regimes de funcionamento da turbina e consequentemente as condies
extremas que as componentes da turbina vo ter que suportar. Como essas condies extremas vo
ser usadas no desenvolvimento das componentes, conclui-se que a escolha das condies de
funcionamento essencial para o desenvolvimento da turbina. Veers
(49)
refere que esses parmetros
de funcionamento devem ser seleccionados de forma a optimizar a energia capturada e a vida til da
turbina e fornece um mtodo para o fazer. Uma primeira considerao, comprovada por resultados
experimentais, que os danos acumulados enquanto a turbina est parada so desprezveis. De
seguida, usando funes de densidade de probabilidade da energia e dos danos da turbina em
funo da velocidade do vento, possvel uma anlise dos regimes de velocidade pretendidos. A
maior dificuldade est na obteno da funo de densidade de probabilidade para os estragos na
turbina em funo da velocidade do vento. Tal deve-se s condies extremas de funcionamento das
turbinas elicas, sujeitas a uma grande quantidade de esforos muito diferentes entre si.
Existem, no entanto, trabalhos nesta rea como Kelley
(50)
que apresenta mtodos computacionais
com resultados prximos dos valores experimentais. Um estudo de grande interesse e abrangente
sobre esta temtica dado por Nijssen
(51)
.
O programa de simulao de estruturas Cosmos normalmente integrado no pacote Solidworks
usado em muitas empresas, inclusive pelo Sandia laboratories com muito bons resultados
(52)
.
No Anexo 1 possvel verificar alguns dos materiais utilizados e as configuraes escolhidas por
diversos fabricantes.
65
6.2 Construo do prottipo
Pretendia-se testar com o prottipo uma ampla gama de configuraes semelhantes s simuladas
em CFD. Para tal seria necessrio que o prottipo permitisse a alterao do nmero de ps e os tipos
de ps utilizadas (com diferentes perfis e corda). Pretendia-se tambm um sistema de baixo custo e
simples para que a sua instalao e manuseamento fosse acessvel e o seu tempo de execuo curto.
Na construo do prottipo, as componentes mais difceis de construir so as ps. Estas requerem
um molde onde a fibra assente e curada. Normalmente so construdos dois moldes sendo a p
fabricada em duas peas que so unidas. Tal como referido anteriormente, no caso das VAWT
giromill o facto de a p ter seco constante sem toro facilita a sua construo. Os moldes podem
ser construdos unindo vrias seces de material cortadas. Usando, por exemplo, pranchas de
madeira com as seces cortadas e unidas, obtm-se moldes para ps de pequenas dimenses de
grande qualidade com um baixo custo, usando tcnicas, ferramentas e materiais amplamente
disponveis
(53).
Neste caso, comeou-se por usar seces de carto canelado recortadas
manualmente. As seces foram alinhadas numa prancha de madeira como mostra a Figura 29.
Figura 29 Obteno das seces em carto canelado (esquerda) e colocao das seces alinhadas numa base para a
produo do molde das ps (direita)
De seguida forrou-se as seces de carto canelado com pedaos de carto para fortalecer a unio
das diferentes seces e por fim forrou-se com uma folha de carto canelado para obter a superfcie
onde aplicar a fibra, ou seja, o molde final.
Como se tinha visto que a fibra aderia fortemente ao carto e que o carto absorvia alguma da
resina perdendo as suas propriedades, aplicou-se uma folha de papel de alumnio por cima do carto
canelado e fibrou-se. Fizeram-se duas meias asas pelo mesmo processo, cortaram-se os excessos de
material, descolou-se as folhas de alumnio e obtiveram-se as peas apresentadas na Figura 30.
Por fim uniram-se as duas meias ps para formar uma p usando a mesma resina usada na
construo das meias ps, obtendo-se o resultado apresentado na Figura 31.
Como se verifica, o produto final tem elevadas tolerncias que podem ser explicadas por vrios
factores. Em primeiro lugar, a utilizao de carto canelado impe algumas dificuldades durante o
recorte dos perfis pois a estrutura interna do carto anisotrpica e, conforme a regio de corte, os
66
perfis utilizados para o molde ficam com maior ou menor resistncia estrutural. Nos locais onde foi
aplicada cola para a construo do molde, o carto canelado alterou as suas caractersticas,
tornando-se mais malevel, como se tivesse humidade. A utilizao de perfis de suporte apenas a
cada 10 cm, deixou uma grande rea do molde sem suporte, na qual a superfcie de carto canelado
abateu ligeiramente. Ainda se tentou preencher as concavidades com material de enchimento como,
por exemplo, espuma de poliuretano, mas sem sucesso. Por fim, como este trabalho foi desenvolvido
nos Aores, onde os materiais compsitos disponveis so maioritariamente usados para a
manuteno de embarcaes navais, usou-se a nica fibra disponvel. Esta, tinha caractersticas
desadequadas para o trabalho em curso. Foi usada uma fibra de vidro de grande espessura, tornando
muito difcil a sua aplicao em zonas com baixo raio de curvatura como o bordo de ataque.
Problemas semelhantes e algumas solues so apresentados por Derek
(52)
entre os quais destacam-
se a utilizao de materiais de enchimento (como, por exemplo, madeira de balsa) para ajudar a
manter a forma desejada e a utilizao de moldes de elevada qualidade.
Figura 30 Meia asa aps cura, retirada do molde (esquerda) e as duas meias asas prontas para serem unidas para
formar uma p (direita)
Figura 31 P finalizada
Para ps mais pequenas, torna-se mais difcil a manuteno de parmetros de tolerncia aceitveis
(54)
. Tentou-se tambm obter um molde em silicone a partir do molde anterior de carto. Tal
facilitaria a produo de um elevado nmero de ps dada a flexibilidade do material usado e o facto
de o silicone no aderir ao material compsito. O silicone lquido de baixa viscosidade disponvel era
demasiado caro. Assim tentou-se obter silicone lquido de baixa viscosidade a partir de silicone
comum dissolvendo-o. Aps vrias tentativas conseguiu-se dissolver o silicone em diluente sinttico
67
e usou-se o material obtido para fazer um molde. Ao solidificar, o molde de silicone encolheu pelo
que a sua aplicao no foi possvel.
As ps acabaram por ser produzidas no departamento de Aeroespacial por um mtodo j utilizado
anteriormente para a produo de asas de pequenas aeronaves. Cortou-se com fio quente, numa
mquina prpria, o interior das ps a partir de uma placa de espuma de isolamento trmico usada
habitualmente em construo civil. O interior de espuma foi usado como molde para a aplicao de
compsito de fibra de vidro. No fim cortaram-se os excessos de material obtendo-se as ps como
mostra a Figura 32.
O comprimento das ps estava limitado pela mquina de corte em cerca de 1m. Optou-se por fazer
as ps com esta dimenso mxima para aumentar o alongamento AR da turbina, diminuindo a
importncia de efeitos tridimensionais como end/tip losses. A avaliao destes efeitos seria depois
contabilizada pelo teste de ps com comprimento inferior (cortando as ps originais).
Figura 32 Ps feitas no Dep. de Aeroespacial finalizadas
As restantes componentes da turbina exigiam uma complexidade de construo que s seria
possvel numa oficina especializada. Fez-se o projecto da turbina com vista sua construo em liga
de alumnio para facilitar o seu manejamento durante a fase experimental deste trabalho. O trabalho
foi entregue a uma oficina escolhida de acordo com o oramento apresentado. A construo total da
turbina demorou cerca de 3 meses, muito mais tempo que o inicialmente indicado de 2 semanas. A
inexistncia dos materiais necessrios, a demora para os mesmos chegarem aos Aores e sucessivos
erros de construo levaram ao atraso. O desenho inicial foi fortemente alterado pelas razes
anteriores e, muitas vezes, por livre arbtrio da prpria oficina. Apesar de todos os inconvenientes, a
turbina finalizada permite, tal como se pretendia, o teste de vrias configuraes. No entanto, o
dimetro do poste central de 5 cm, diferente do testado em CFD de 10 cm. A turbina apresentada
na Figura 33, Figura 34 e Figura 35.
Como se pode ver pelas figuras, as diferentes configuraes eram montadas em dois discos que
permitiam configuraes at 6 ps. O acoplamento realizado por 2 conjuntos porca-parafuso por
cada brao. A turbina fixa ao poste por 3 conjuntos porca-parafuso. O acoplamento dos braos com
68
as ps feito usando 2 conjuntos porca-parafuso por brao, fixos em braadeiras metlicas de
alumnio (que facilmente se adaptam forma do perfil).
Figura 33 Turbina finalizada (sem instrumentos)
Figura 34 Detalhe do apoio do taqumetro, travo e brao das ps
Figura 35 Detalhe do acoplamento das ps aos braos
O poste fixo parede por braadeiras metlicas. Na base acoplado a uma dobradia permitindo
baixar toda a turbina para mudar a configurao pretendida ou outras tarefas.
69
7 Teste do prottipo VAWT
7.1 Resumo de testes existentes
Existem vrios resultados experimentais de VAWTs, principalmente do tipo Darrieus. Alguns destes
resultados j foram apresentados anteriormente mas existem outros que so agora apresentados
dada a relevncia na aplicao deste trabalho.
Para testar a eficincia de uma VAWT, necessrio controlar o coeficiente de velocidade perifrica
, aplicando um binrio M de intensidade conhecida, contrrio ao produzido pela turbina.
necessrio medir a velocidade do vento
r
i
o
N
m
Peso kg
Resultados
Mdia
72
O binrio M
medido varia consoante o peso P utilizado. Tal deve-se existncia de esforos nos
rolamentos que levam a um aumento de M . Desta forma, sabendo que diferentes configuraes da
turbina apresentam diferentes foras normais exercidas sobre o poste, conclui-se que M variar
conforme a configurao usada. Note-se tambm a maior variao nos resultados obtidos com pesos
maiores, devido diminuio do tempo t . Procedeu-se a um novo ensaio. Colocou-se a turbina na
posio vertical e, usando um berbequim, acelerou-se a mesma at se atingir uma velocidade
angular
0
e em rad/s constante. De seguida desacoplou-se o berbequim e, sabendo o tempo t em s
que a turbina demora a parar, obteve-se o binrio M em N.m a partir da frmula
t
I
M
0
e
= , (41)
sendo I o momento de inrcia da turbina sem ps em kg.m
2
. Obteve-se o Grfico 22.
Grfico 22 Medio 2 do binrio aplicado pelo travo
O binrio medido aproximadamente constante, ao contrrio das medies feitas inicialmente.
Desta forma no parece haver uma variao importante de M com
e
.
Usou-se o mesmo
procedimento para calibrar o travo. Os resultados obtidos so apresentados no Grfico 23.
Grfico 23 Medio 3 do binrio aplicado pelo travo
medida que se aumenta o binrio aplicado pelo travo, torna-se cada vez mais difcil a medio
do mesmo. A turbina pra mais rapidamente, dificultando a medio de t . Para evitar esta situao,
-0,11
-0,1
-0,09
-0,08
-0,07
-0,06
-0,05
-0,04
0,00 10,00 20,00 30,00 40,00 50,00
B
i
n
r
i
o
N
m
Velocidade angular inicial rad/s
Resultados
Mdia
-0,35
-0,3
-0,25
-0,2
-0,15
-0,1
-0,05
0 0,2 0,4 0,6 0,8 1 1,2
B
i
n
r
i
o
N
m
Peso kg
Resultados
Mdia
Regresso linear de 2
ordem
y=-0,101x
2
-0,075x-0,083
73
tentou-se usar uma velocidade angular inicial
0
e superior mas para atingir valores elevados de
rotao com um grande binrio aplicado exigiam-se tenses demasiado altas no sistema de
acoplamento berbequim-turbina. Este sistema, que inicialmente se pretendia que funcionasse com
rodas dentadas, acabou por ser feito pelo acoplamento de uma superfcie de borracha. Este material,
quando sujeito s grandes tenses exigidas desgasta-se muito rapidamente tornando a sua operao
muito difcil. Este factor levou aos valores mximos medidos. No entanto, a partir dos resultados
obtidos, foi possvel obter uma medida do binrio aplicado em funo da fora aplicada ao travo.
Note-se ainda que, mesmo com as elevadas velocidades aplicadas turbina (as medies foram
feitas sem as ps), nunca se notou qualquer tipo de movimento que indicasse um desequilbrio na
sua construo. Assim assumiu-se que o eixo principal de inrcia estava alinhado com o eixo de
rotao da turbina.
Com a medio destes parmetros, obtiveram-se os dados necessrios para o teste das diferentes
configuraes da turbina. As medies foram feitas ao ar livre, numa zona elevada, livre de
obstculos. Foi difcil conciliar as medies com o clima dada a altura do ano (incio dos testes em
Dezembro) e as condies climatricas altamente variveis dos Aores. Para as medies era
necessrio que estivesse de dia, sem chuva, com vento aproximadamente constante. Estas condies
raramente se apresentaram.
Para os testes, aplicava-se um peso conhecido no travo e aguardava-se que as condies de vento
e rotao da turbina permanecessem constantes durante um curto intervalo de tempo e obtinham-
se os valores de velocidade do vento
V , menor o declive da
regresso linear respectiva, o que seria de esperar pois
2
3
2
1
=
V
C
AV
M
C
P P
e
.
(42)
O tratamento dos dados recolhidos, nomeadamente o clculo de propagao de erros, feito no
Anexo 4. Existe uma condio de equilbrio para 8 =
V m/s e 083 , 0 = M Nm em
073 , 0 440 , 0 exp = erimental . Para comparao com os resultados numricos, procedeu-se
interpolao dos resultados obtidos em CFD (Anexo 3) de 4 , 0 = e 55 , 0 = para 083 , 0 = M Nm
tendo-se obtido 386 , 0 =
CFD
. Tal condio encontrada para 083 , 0 = M Nm mas, tal como se viu
anteriormente, esta configurao apresenta uma fora normal
n
F que exercida sobre os
rolamentos afectando M . De acordo com o Grfico 20, a intensidade de
n
F
estaria
aproximadamente entre 10 e 35 N o que de acordo com o Grfico 21 corresponde a um M superior
a 0,2. Durante os testes, a medio de
n
F poderia ser realizada, por exemplo, pela aplicao de
extensmetros ou acelermetros no poste principal. No entanto, sem uma medio de
n
F a
comparao de resultados feita para 083 , 0 = M Nm.
Obtiveram-se tambm resultados experimentais para a configurao 4: NACA0018, 2 ps com
corda de 10cm. Tal como na situao anterior, procedeu-se a uma reviso dos resultados obtidos em
CFD. Os resultados experimentais so apresentados no Grfico 26.
Grfico 26 Resultados experimentais obtidos para a Configurao 4
5,0E-04
1,0E-03
1,5E-03
2,0E-03
0,15 0,20 0,25 0,30 0,35 0,40 0,45
C
p
6
7
8
9
10
V
m/s
76
Tal como para a configurao 5, procedeu-se ao tratamento dos dados recolhidos (Anexo 4).
Calculou-se uma condio de equilbrio para 8 =
V m/s e 083 , 0 = M Nm em
043 , 0 346 , 0 exp = erimental . Para comparao com os resultados numricos, procedeu-se
interpolao dos resultados obtidos em CFD (Anexo 3) de 3 , 0 = e 4 , 0 = para 083 , 0 = M Nm
tendo-se obtido 300 , 0 =
CFD
.
Os resultados obtidos para ambas as configuraes so apresentados no Quadro 9.
Teste
Condio de equilbrio para V
=8m/s, M=0,083Nm
CFD
experimental
Configurao 4 0,300 0,3460,043
Configurao 5 0,386 0,4400,073
Quadro 9 Resultados experimentais e numricos
Para validar os resultados obtidos importa fazer uma avaliao dos erros obtidos durante as
experincias. Em primeiro lugar, a turbina testada um sistema tridimensional, sujeito a efeitos que
no so medidos durante as simulaes bidimensionais em CFD como, por exemplo, tip vortices, ou a
existncia de braos de suporte. As ps construdas tm tolerncias, rugosidade, no so rgidas, e
no esto perfeitamente tangenciais trajectria descrita. A prpria turbina no rgida, afectando
as condies de escoamento. Esta situao, de forte vibrao da turbina, de tal forma grave que
inviabiliza o uso da mesma para velocidades angulares superiores s utilizadas durante os ensaios. O
poste central tem menos 5 cm de raio que os modelos obtidos em CFD e possui estruturas acopladas,
como os discos de suporte dos braos, que afectam o escoamento. Relativamente s condies de
teste, o local utilizado, afecta o escoamento mesmo que se tenha tido o cuidado de escolher uma
localizao privilegiada, elevada e com poucos obstculos circundantes. O ar tem propriedades
variveis como a temperatura ou a densidade que afectam as suas caractersticas. Alm disso, o
vento, varivel em intensidade e direco, o que se tentou evitar recolhendo valores apenas
quando as condies de vento se mantivessem relativamente estveis durante o intervalo de
medio. Relativamente aos equipamentos utilizados, o anemmetro o nico equipamento
calibrado. no entanto, um equipamento com alguma inrcia funcionando como filtro passa-baixo. A
velocidade angular foi medida com recurso a um taqumetro de bicicleta, sem qualquer certificado
de calibrao. O cronmetro utilizado para determinar o tempo do ensaio tambm no tinha
qualquer certificado de calibrao. De forma a eliminar erros aleatrios fizeram-se bastantes ensaios
para cada condio testada de forma a avaliar a qualidade da medida e a sua preciso
(59)
. Nos
grficos apresentados, muitos dos resultados obtidos encontram-se sobrepostos devido resoluo
dos instrumentos utilizados.
77
Importa referir que no foi feita uma anlise de erros aos valores experimentais de M dado que
parte das medies (por exemplo, o peso de certas componentes) foi efectuada na oficina onde se
construiu a turbina pelos funcionrios da mesma. No houve acompanhamento das medies
efectuadas nem foi possvel ter a certeza de quais os instrumentos utilizados nas medies.
Tendo em conta estas consideraes, o facto dos resultados numricos coincidirem ou serem
muito prximos do intervalo dos resultados experimentais, parece apontar para a validao dos
resultados numricos. No entanto h que ter em ateno que tais resultados foram obtidos apenas
para duas configuraes, numa gama muito restrita de valores de velocidade do vento
V
e
coeficiente de velocidade perifrica . Tal deveu-se, como indicado anteriormente, dificuldade em
obter condies de teste adequadas.
Durante a realizao deste trabalho, procedeu-se ainda a outros tipos de testes e experincias
relacionados com o mesmo mas que, devido a restries de espao so apenas apresentados muito
sumariamente. A ideia inicial de utilizar uma mquina elctrica para testar a turbina foi posta em
prtica pela aquisio de mltiplos tipos de motores e geradores em ferro-velho/empresas de
reciclagem. A maioria encontrava-se em ptimas condies, funcionando sem problemas e foram
obtidos gratuitamente excepo do trabalho necessrio para os encontrar e desmontar dos
respectivos equipamentos originais. Fizeram-se testes extensivos com o objectivo de obter as curvas
de funcionamento de alguns com vista sua utilizao, como mostra a Figura 39.
Figura 39 Teste de mquina elctrica
Foi utilizado um berbequim de velocidade varivel acoplado ao motor. Os terminais do motor eram
ligados a um potencimetro e a um multmetro. A medio da velocidade de rotao era realizada
com recurso ao mesmo taqumetro de bicicleta que foi usado na turbina. Foram ainda realizados
outros ensaios a diversos materiais e tcnicas para a construo, por exemplo, de ps.
78
8 Concluses e trabalho futuro
8.1 Concluses
O trabalho desenvolvido permitiu uma profunda compreenso do tema estudado, tendo-se
retirado vrias concluses:
- As VAWTs, especialmente do tipo lift, so sistemas inerentemente complexos, de difcil
anlise. So sistemas relativamente pouco estudados, com pouca informao disponvel e
com falta de dados e modelos para o seu estudo.
- A simulao de VAWTs em CFD, usando o modelo kw SST, oferece uma descrio de
fenmenos fsicos esperada e com resultados semelhantes aos obtidos
experimentalmente. No entanto, a gama restrita de resultados prticos obtidos no
permite analisar o modelo completamente. De qualquer forma, os resultados em CFD
obtidos so uma importante fonte de informao para trabalhos nesta rea.
- Das configuraes testadas, a configurao 4 (2 ps NACA0018 com mm c 100 = e
m r 5 , 0 = ) foi a que apresentou melhores resultados. No entanto, a gama de condies
testadas foi limitada e no foi possvel comprovar todos os resultados numricos com
ensaios prticos. A escolha de uma configurao ptima para o desenvolvimento de uma
turbina depende de muitos outros factores como resistncia estrutural ou factores
socioeconmicos que no foram totalmente avaliados.
- A obteno de resultados experimentais de VAWTs tem de ser feita usando equipamentos
criteriosamente desenvolvidos. As foras aplicadas na turbina durante o seu
funcionamento so intensas e difceis de quantificar e a existncia de fenmenos, como a
vibrao observada neste trabalho, afectam as medies obtidas. Os instrumentos usados
devem ter resoluo e taxa de amostragem suficiente para a avaliao dos parmetros
pretendidos.
- O teste em condies reais difcil, principalmente devido natureza varivel do vento.
Assim, os sistemas desenvolvidos para o teste de VAWTs em condies reais devem ser
79
munidos de instrumentos automatizados que permitam a recolha de dados durante longos
perodos de tempo sem a presena de um utilizador. A prpria turbina deve ser preparada
para resistir aos elementos, podendo ficar no exterior durante longos perodos de tempo.
Embora no fizesse parte dos objectivos desta tese, a concluso mais pertinente talvez se prenda
com o futuro das VAWT. Na opinio do autor, estes sistemas oferecem grandes vantagens em
relao s HAWT, principalmente em locais com ventos fortemente variveis como o caso de
ambientes urbanos. Tal comprovado pelos diversos resultados obtidos e recolhidos e pelo forte
desenvolvimento recente de VAWTs comerciais. As assinalveis diferenas obtidas no coeficiente de
potncia
P
C com a mudana de pequenos parmetros, parecem promissoras para que haja a
possibilidade de serem descobertas configuraes com caractersticas desejveis que tornem as
VAWTs uma soluo comprovada e de grande valor.
80
8.2 Trabalho futuro
Este trabalho permitiu uma maior sensibilidade para as questes inerentes ao desenvolvimento de
VAWTs. Sendo um conceito relativamente pouco estudado, existem vrias reas onde ser do maior
proveito um estudo mais aprofundado, facilmente integrvel em trabalhos de mestrado ou
doutoramento. O forte crescimento do sector da energia elica e, em particular, o esperado
desenvolvimento da micro-gerao, tornar estas reas de estudo de grande interesse, promovendo
uma grande visibilidade dos estudos efectuados e respectivos centros de investigao. Importa
salientar a capacidade do Dep. de Aeroespacial para o desenvolvimento destas competncias, dadas
as similaridades entre as reas de estudo, a competncia do corpo docente e a existncia de um
laboratrio com condies adequadas. Tal aposta seria certamente recompensada dada a grande
procura de recursos humanos e conhecimento nesta rea. Em especial, o desenvolvimento de
sistemas menos convencionais como as VAWTs, dever trazer grandes retornos iniciais
(2)
.
Alguns dos trabalhos possveis para desenvolvimento so:
- Estudo dos modelos existentes para a previso do comportamento de VAWTs. Actualmente,
apesar de vrios modelos terem sido desenvolvidos, a sua aplicabilidade e certeza ainda
difcil de quantificar. Sugere-se o estudo de vrios modelos, inclusive computacionais e
comparao com resultados experimentais. No caso de modelos computacionais, importa
tentar quantificar diferentes parmetros como, por exemplo, as caractersticas da malha.
Um grande objectivo seria facilitar a escolha de modelos a usar em trabalhos posteriores.
- Desenvolvimento de uma bancada de testes para VAWTs, facilmente adaptvel a diferentes
configuraes. Sugere-se o desenvolvimento de um sistema semelhante ao que foi
desenvolvido nesta dissertao para a obteno de resultados experimentais. Espera-se, no
entanto, o desenvolvimento e implementao de sistemas de instrumentao que permitam
uma recolha de dados automtica e de elevada qualidade. Este sistema poderia ser
desenvolvido para utilizao exterior, em tnel de vento ou em outras condies
possibilitando a implementao de muitos estudos relativos a vrios parmetros como os
tipos de perfis usados, solidity o , alongamento AR , inrcia da VAWT, rudo, etc.
- Estudo de modelos de fadiga e ensaios no-destrutivos para aplicao em turbinas elicas. O
desenvolvimento de sistemas que permitam uma diminuio dos custos de manuteno e o
aumento do conhecimento sobre as condies reais da turbina so de grande importncia
na indstria de energia elica. Alguns trabalhos do Dep. de aeroespacial so actualmente
desenvolvidos nesta rea, como a implementao de sistemas de instrumentao em
compsito.
81
- Comparao de VAWTs e HAWTs no s do ponto de vista terico de comparao de
vantagens e desvantagens como tem largamente sido feito at hoje, mas tambm do ponto
de vista, energtico, por exemplo, atravs da simulao de resultados de diferentes turbinas
em diferentes ambientes.
- Desenvolvimento de perfis especficos para VAWTs. Embora tenha havido algum
desenvolvimento nesta rea, as condies especficas de funcionamento de VAWTs,
atingindo-se elevados ngulos de ataque o e o funcionamento a nmeros de Reynolds Re
relativamente baixos, exigem o estudo, e possvel desenvolvimento, de perfis para estas
condies. A inexistncia de curvas, comprovadas por ensaios prticos, para a maior parte
dos perfis em regimes de funcionamento de VAWTs, um grande entrave actualmente.
82
8.3 Notas finais
Durante a execuo deste trabalho, houve sempre uma grande curiosidade em torno do
desenvolvimento deste tipo de equipamento das pessoas em geral. Principalmente aps a execuo
da turbina e a sua fixao no local de testes, eram frequentes as questes e os pedidos de
informao acerca do seu funcionamento e da possibilidade de adquirir um equipamento do gnero.
Em geral, as pessoas mostravam-se muito optimistas em relao a este tipo de tecnologia e em
relao energia elica em geral. Estas observaes, embora sem medida cientfica, esto em
consonncia com os resultados obtidos por Dalton
(60)
, que indica que existe, em geral, uma forte
aceitao das energias renovveis e da energia elica em particular. Indica ainda que os turistas
podem ser um forte incentivo para a implementao de energias renovveis e que, apesar das
pessoas com mais idade apresentarem melhores conhecimentos relativamente a este assunto e
serem os maiores utilizadores do turismo ecologicamente responsvel, so os mais jovens que esto
mais dispostos a instalar este tipo de sistemas.
83
Bibliografia
Existem poucas referncias relativamente ao desenvolvimento de VAWTs tipo giromill. Os estudos
mais prximos deste tipo foram realizados na dcada de 70/80 pelo laboratrio Sandia, cujos
resultados podem ser encontrados em www.sandia.gov. Os restantes resultados foram obtidos a
partir do motor de busca www.sciencedirect.com. Foram ainda consultados alguns trabalhos nas
bibliotecas do Instituto Superior Tcnico ou atravs do sistema Fnix, principalmente teses de
mestrado. Foram usados vrios livros, a maioria generalistas sobre as vrias reas em estudo. Livros
especficos usados foram apenas as referncias (5) e (2), sendo o ltimo o nico especfico sobre o
desenvolvimento de VAWTs e ambos difceis de encontrar no mercado nacional.
Importante informao foi tambm obtida atravs de vrias agncias como a EWEA, atravs do site
www.ewea.org, seminrios atendidos e do contacto com profissionais do sector.
Bibliografia
1. Kirke, Brian Kinloch. Evaluation of self-starting vertical
axis wind turbines for stand-alone applications. PhD thesys,
Griffith University. 1998.
2. Paraschivoiu, Ion. Wind turbine design (with emphasis
on Darrieus concept). cole Polytechnique de Montral.
ISBN 2-553-00931-3, 2002.
3. Islam, Mazharul, Ting, David e Fartaj, Amir.
Aerodynamic models for Darrieus-type straight-bladed
vertical axis wind turbines. Renewable & sustainable energy
reviews. 2008, Vols. 12 p.1087-1109.
4. Eriksson, Sandra, Bernhoff, Hans e Leijon, Mats.
Evaluation of different turbine concepts for wind power.
Renewable & sustainable energy reviews. 2006, Vols. 12
p.1419-1434.
5. Manwell, J. F., McGowan, J. G. e Rogers, A.L. Wind
energy explained (theory,design and application). John
Wiley & Sons Ltd. ISBN 0-470-84612-7, 2002.
6. White, Frank M. Fluid Mechanics. McGraw Hill. ISBN
978-0-07-128645-9, 6 edio.
7. EWEA. Wind energy: the facts. EWEA. 2009.
8. Instituto de Meteorologia. Resumos climatolgicos de
aerdromo para os aeroportos do arquiplago dos Aores.
Instituto de Meteorologia. 2009.
9. EWEA. Oceans of opportunity - Harnessing Europe's
largest domestic energy resource. EWEA. 2009.
10. . Wind at work: Wind energy and job creation in
the EU. EWEA.
11. . The economics of wind energy. EWEA. Maro
2009.
12. Walford, Christopher. Wind turbine reliability:
Understanding and minimizing wind turbine operation and
maintenance costs. Sandia Laboratories report SAND2006-
1100. 2006.
13. EWEA. Wind energy and electricity prices - Exploring
the "merit order effect". EWEA. 2010.
84
14. . Integrating wind - Developing Europe's power
market for the large-scale integration of wind power.
EWEA. 2009.
15. Mott Macdonald. UK electricity generation costs
update. s.l. : Department of Energy and Climate Change
(UK), Junho 2010.
16. California Energy Comission. Comparative costs of
California central station electricity generation
technologies. California Energy Comission report. 2007.
17. Saha, U.K., Thotla, S. e Maity, D. Optimum design
configurantion of Savonius rotor through wind tunnel
experiments. Journal of wind engineering and industrial
aerodynamics. 2008, Vols. 96 p.1359-1375.
18. Riegler, Hannes. HAWT versus VAWT: Small VAWTs
find a clear niche. REFOCUS. 2003, Vol. Julho/Agosto 2003.
19. Castro, Rui M. G. Introduo energia elica. IST.
1998.
20. Furtado, Antnio Aguiar. Energia e desenvolvimento
sustentvel - Plano de investimentos da EDA em energias
renovveis. Apresentao Green Islands em Ponta Delgada.
2010.
21. Sahin, Ahmet. Progress and recent trends in wind
energy. Progress in energy and combustion science. 2004,
Vols. 30 p.501-543.
22. EEG - Empresa de electricidade e gaz. EEG. [Online]
www.eeg.eda.pt.
23. Rodrigues, lvaro. Wind farms in Portugal. INEGI.
2007.
24. Johnson, Scott, van Dam, C.P. e Berg, Dale. Active
load control techniques for wind turbines. Sandia
Laboratories report SAND2008-4809. 2008.
25. Fujisawa, Nobuyuki e Shibuya, Satoshi. Observations
of dynamic stall on Darrieus wind turbine blades. Journal of
wind engineering and industrial aerodynamics. 2001, Vols.
89 p.201-214.
26. Sheldahl, Robert e Klimas, Paul. Aerodynamic
characteristics of seven symmetrical airfoil sections through
180-degree angle of attack for use in aerodynamic analysis
of vertical axis wind turbines. Sandia laboratories report
SAND80-2114. 1981.
27. Sarkar, Sunetra e Venkatraman, Kartik. Influence of
pitching angle of incidence on the dynamic stall behavior of
a symmetric airfoil. European journal of mechanics. 2008,
Vols. 27 p.219-238.
28. Geissler, W. e Haselmeyer, H. Investigation of
dynamic stall onset. Aerospace science and technology.
2006, Vols. 10 p590-600.
29. Wang, Shengyi, et al., et al. Numerical investigations
on dynamic stall of low Reynolds number flow around
oscillating airfoils. Computers & fluids. 2010, Vols. 39
p.1529-1541.
30. Larsen, J.W., Nielsen, S.R. e Krenk, S. Dynamic stall
model for wind turbine airfoils. Journal of fluids and
structures. 2007, Vols. 23 p.959-982.
31. Sun, M. e Sheikh, S.R. Dynamic stall suppression on
an oscillating airfoil by steady and unsteady tangential
blowing. Aerospace science and technology. 1999, Vols. 6
p.355-366.
32. Thakkar, D. e Ganguli, R. Single-crystal piezoceramic
actuation for dynamic stall suppression. Sensors and
actuators A. 2006, Vols. 128 p.151-157.
33. Yurdusev, M.A., Ata, R. e etin, N.S. Assessment of
optimum tip speed ratio in wind turbines using artificial
neural networks. Energy. 2006, Vols. 31 p.2153-2162.
34. etin, N.S., et al., et al. Assessment of optimum tip
speed ratio for wind turbines. Mathematical and
computational applications. 2005, Vols. 10 N.1 p.147-154.
35. Howell, Robert, et al., et al. Wind tunnel and
numerical study of a small vertical axis wind turbine.
Renewable Energy. 2010, Vols. 35 p.412-422.
36. Li, Ye e Calisal, Sander. Three-dimensional effects
and arm effects on modellong a vertical axis tidal current
turbine. Renewable energy. 2010, Vols. 35 p.2325-2334.
37. Brederode, Vasco de. Fundamentos de aerodinmica
incompressvel. Edio do autor. ISBN 972-97402-0-8, 1997.
38. Klimas, Paul. Tailored airfoils for vertical axis wind
turbines. Sandia laboratories report SAND84-1062. 1992.
39. Ashwill, Thomas. Measured data for the Sandia 34-
meter vertical axis wind turbine. Sandia laboratories report
SAND91-2228. 1992.
40. van Bussel, Gerard. The development of Turby, a
small VAWT for the built environment. Global wind energy
conference. 2004.
41. AllSmallWindTurbines. [Online] 2010.
http://www.allsmallwindturbines.com/.
42. Sheldahl, Robert. Comparison of field and wind
tunnel Darrieus wind turbine data. Sandia laboratories
report SAND80-2469. 1981.
43. Larsen, Kari. Recycling wind turbine blades.
Renewable energy focus. 2009, Vols. Jan/Fev 2009 p.70-73.
44. Kensche, Cristoph W. Fatigue of composites for wind
turbines. International journal of fatigue. 2006, Vols. 28
p.1363-1374.
45. Reuterlov, Stefan. The role of sandwich composites
in turbine blades. Reinforced plastics. 2002, Vols. Maro
2002 p.32-34.
46. Wallace, Jack e Dawson, Mark. O&M strategies: wind
turbine blades. Renewable energy focus. 2009, Vols.
Maio/Junho 2009 p.36-41.
47. Ong, Cheng-Huat e Tsai, Stephen. The use of carbon
fibers in wind turbine blade design: a seri-8 blade example.
Sandia laboratories report SAND2000-0478. 2000.
48. Sharman, David. Getting real about small wind. Part3.
Renewable energy focus. 2010, Vols. Maro/Abril 2010
p.36-39.
49. Veers, Paul S. A general method for fatigue analysis
of vertical axis wind turbine blades. Sandia laboratories
report SAND82-2543. 1983.
50. Neil D. Kelley, Herbert J. Sutherland. Damage
estimates from long-term structural analysis of a wind
turbine in a U.S. wind farm environment. s.l. : Wind energy,
1997.
51. Nijssen, R. P. L. Fatigue life prediction and strenght
degradation of wind turbine rotor blade composites. Sandia
laboratories report SAND2006-7810P. 2007.
52. Berry, Derek. TX 100 manufacturing - final project
report. Sandia laboratories report SAND2007-6066. 2007.
53. Habali, S.M. e Saleh, I.A. Local design, testing and
manufacturing of small mixed airfoil wind turbine blades of
glass fiber reinforced plastics. Part II: Manufacturing of the
blade and rotor. Energy conversion & management. 2000,
Vols. 41 p.281-298.
85
54. Clausen, P.D. e Wood, D.H. Research and
development issues for small wind turbines. Renewable
energy. 1999, Vols. 16 p.922-927.
55. Zayas, Jose R., Jones, Perry L. e Holman, Adam. CX-
100 and TX-100 blade field tests. Sandia laboratories report
SAND2005-7454. 2005.
56. Beattie, Alan G. e Rumsey, Mark. Non-destructive
evaluation of wind turbine blades using an infrared camera.
Sandia laboratories report SAND98-2824C. 1998.
57. Sutherland, Herbert, et al., et al. The application of
non-destructive techniques to the testing of a wind turbine
blade. Sandia laboratories report SAND93-1380. 1994.
58. Rumsey, Mark e Paquette, Joshua. Structural health
monitoring of wind turbine blades. s.l. : Sandia laboratories,
2007.
59. Azinheira, J. e Fonseca, A. Instrumentao. IST.
60. Dalton, G. J., Lockington, D. A. e Baldock, T. E. A
survey of tourist attitudes to renewable energy supply in
Australian hotel accomodation. Renewable energy. 33
p.2174-2185, 2008.
61. Veers, Paul. Modeling stochastic wind loads on
vertical axis wind turbines. Sandia Laboratories report
SAND83-1909. 1984.
62. EWEA. Pure Power - Wind energy targets for 2020
and 2030. EWEA. 2009.
63. Patrikar, Rajendra M. Modeling and simulation of
surface roughness. Journal of applied surface science. 2004,
Vols. 228 p.213-220.
64. Jureczko, M., Pawlak, M. e Mezyk, A. Optimisation of
wind turbine blades. Jornal of materials processing
technology. 2005, Vols. 167 p.463-471.
65. Marsh, George. Bigger blades - the carbon option.
Reinforced plastics. 2002, Vols. Maro 2002 p.20-30.
66. Swiatecki, Sebastian. Turbine makers make their own
blades. Reinforced plastics. 2002, Vols. Maro 2002 p.36-40.
67. Habali, S.M. e Saleh, I.A. Local design, testing and
manufacturing of small mixed airfoil wind turbine blades of
glass fiber reinforced plastics. Part I: Design of the blade
and root. Energy conversion and management. 2000, Vols.
41 p.249-280.
68. Underwood, C.P., et al., et al. Renewable-energy
clusters for remote communities. Applied energy. Janeiro
2007, Vols. 84 p.579-598.
69. Moreno, Ferran Gual. Design of a small wind
generator. Dissertao para obteno do grau de mestre
em engenharia mecnica - IST. 2008.
70. Neto, Joo A. Yield optimization based on wind
resource. Dissertao para a obteno do grau de mestre
em engenharia mecnica - IST. 2010.
71. Penedo, Ricardo Jorge Marques. Optimizao
aeroelstica de uma turbina elica urbana. Dissertao
para obteno do grau de mestre em engenharia
aeroespacial - IST. 2009.
72. Prowell, Ian e Veers, Paul. Assessment of wind
turbine seismic risk: Existing literature and simple study of
tower moment demand. Sandia Laboratories report
SAND2009-1100. 2009.
73. Kadlec, Emil G. Characteristics of future vertical axis
wind turbines. Sandia laboratories report SAND79-1068.
1982.
74. Brumioul, Nicolas. Evaluation of aerodynamic criteria
in the design of a small wind turbine with the lifting line
model. Dissertao de mestrado para a obteno do grau
de mestre em engenharia aeroespacial - IST. 2010.
75. Westerholm, Gary. The wind energy frontier.
Refocus. 2004, Vols. Julho/Agosto 2004 p.42-45.
76. Marsh, George. Tilting at windmills. Refocus. 2005,
Vols. Setembro/Outubro 2005 p.37-42.
77. Islam, Mazharul, Fartaj, Amir e Ting, David. Current
utilization and future prospects of emerging renewable
energy applications in Canada. Renewable & sustainable
energy reviews. 2004, Vols. 8 p.493-519.
78. Fortunato, B., Mummolo, G. e Cavallera, G.
Economic optimisation of a wind power plant for isolated
locations. Solar energy. Fevereiro 1997, Vols. 60 N.6 p.347-
358.
79. Bishop, Justin e Amaratunga, Gehan. Evaluation of
small wind turbines in distributed arrangement as
sustainable wind energy option for Barbados. Energy
conversion & management. 2008, Vols. 49 p.1652-1661.
80. Habali, S. e Saleh, I. Design and testing of small mixed
airfoil wind turbine blades. Renewable energy. 1995, Vols. 6
N.2 p.161-169.
81. Bechly, M. e Clausen, P. Structural design of a
composite wind turbine blade using finite element analysis.
Computers & structures. 1997, Vols. 63 N.3 p.639-646.
82. Lanzafame, R. e Messina, M. Power curve control in
micro wind turbine design. Energy. 2010, Vols. 35 p.556-
561.
83. Ata, R. e Kocygit, Y. An adaptive neuro-fuzzy
inference system approach for prediction of tip speed ratio
in wind turbines. Expert systems with applications. 2010,
Vols. 37 p.5054-5460.
87
Anexos
89
Anexo 1 Sistemas VAWT similares
Criou-se uma base de dados de VAWTs de pequenas dimenses, disponveis no mercado que foi utilizada ao longo deste trabalho.
Marca Modelo
Potncia
nominal
Produo
Potncia
real
Preo Altura Peso
Vel. de
arranq
ue
Vel.
nomin
al
Material Gerador Observaes rea
W kWh/ano W m m kg m/s m/s
m
Aeolos
V 300w 300
1,4 1,2 55 2
PM.
1,7
V 600w 600
1,6 1,3 76 2
2,1
V 1000w 1000
2,4 2 140 2,5
4,8
V 3000w 3000
3,6 3 295 3
10,8
V 5000w 5000
3,8 4 370 3
15,2
V 10kw 10000
6,6 6 620 3
39,6
Aerocatcher
Vertikalrotor 0.3 300
998 1,4 1,2 50
Sistema hbrido
Darrieus + Savonius.
1,7
Vertikalrotor 0.5 500
1790
CleanField Energy CleanField 3,5kW 3000 2000 228
3 2,75 245
Ps de
FV.
DD PM
trifsico.
Darrieus. 8,3
FlexiEnergy
Flex500vawt 500 265 30
35 3 17 Al. B DD. Descontinuado.
Black300 vawt 300
780 1,6 1 20 2 20 FV. PM. Darrieus 3 ps. 1,6
Goodridge
engineering
Forgen 500 500 29 3 490 0,31 0,2
2 20
PM DD.
Savonius 3 ps.
Elevada resistncia.
0,1
Forgen 1000 1000 79 9 742 0,465 0,3
0,1
HelixWind
S322
1962 224 6345 2,65 1,21 135 3,6
Al e ao.
PM.
Savonius helicoidal
3,2
S 594
3362 384 11762 4,87 1,21 635 3,6
5,9
D100
600 68 3714 1 1 35 3,5
FV e ao. Darrieus helicoidal
1,0
D361
3168 362 8500 1,9 1,9 130 1,5
3,6
Hi-vawt
DS 300 300 91 10 1010 1,209 1,245 30 3 15,5
Al e ao. PM trifsico.
Hibrido Savonius
Darrieus 3 ps.
1,5
DS 1500 1500 876 100 5950 3,9 2,8 350 2,5 15 10,9
DS 3000 3000 2200 251 8960 5,6 4 680 2,5 15 22,4
SEaB Energy Windbuster 5000 5000 571 22000 6 3,25
3,5 16
FV, Al e
ao.
DD PM.
Darrieus helicoidal
com 5 ps.
19,5
Windspire
Standard 1200 2000 228 5000 6,1 1,2 283 3,8 12,5
Al e ao. B PM DC.
7,3
Extreme 1200 2000 228 5000 4 1,2 257 3,8 17 4,8
Quietrevolution
qr12
50000 5708
12 7
Lanamento em 2012.
84,0
qr2.5
3800 434 11500 2,5 2,5
6,3
90
qr5
10000 1142 35000 5 3,1
4,5 16 FC. DD PM.
Preo apenas da
turbina e controlador.
J comercializado.
15,5
Ropatec
Easy vertical 1000 306 35
1,15 1,8 130 3
Darrieus com 3 ps.
2,1
Simply vertical 3000 1314 150
2 3,3 450 3
6,6
Maxi vertical 6000 2628 300
2,5 4,7 760 3
11,8
Big Star vertical 20000 4380 500
4,3 8 3600 3
Darrieus com 5 ps. 34,4
Turby Turby 2500 1250 143 18189 2,7 2 136
Darrieus helicoidal
com 3 ps.
5,4
UrbanGreenEnergy
UGE-1K 1000 1400 160 5770 2,7 1,8 175 3 30
FV e FC. PM.
4,9
UGE-4K 4000 5000 571 15400 4,4 3 444 3,5 30 13,2
UGE-600 600 800 91 4420 1,6 1,38 81,6 3,5 32 2,2
Venco Power
Twister 300T 300 105 12 4400 1 1 35 3,5 25
FV e ao. DD PM.
Darrieus helicodal
com 3 ps.
1,0
Twister 1000T 1000 876 100 6750 1,9 1,9 150 3,5 20 3,6
Vertikon H50 50000 17520 2000 175000 12,5 12
2,5
FV e ao.
Assncrono
com caixa
de
velocidades.
Darrieus com 3 ps. 150,0
Vertical Wind
Energy
3kw 3000 3000 342 15000 2,5 2,5
Darrieus com 3 ps. 6,3
WePOWER
Falcon 600W 600 700 80 4300 1 1,3 88 2,7
DD PM. Darrieus com 5 ps.
1,3
Falcon 1.2kW 1200 1576 180 6190 2 1,78 160 2,7
3,6
Falcon 3.4kW 3400 1752 200 15477 3,6 3 633 2,7
10,8
Falcon 5.5kW 5500 3933 449 23216 4,6 4 985 2,7
18,4
Falcon 12kW 12000 8760 1000 41000 6,2 6 1905 2,7
37,2
Windside
WS - 0.3 A 1800 260 30 5570 1,5 0,3 87 3
PM DD.
VAWT do tipo drag
para condies
extremas de
temperatura e
intensidade do vento.
0,5
WS - 0,3 B
1,5 0,3 43 2,8
0,5
WS - 0,15 B
38 3,8
0,0
WS-2B
1720 196 24060 2 0,5 550 2
1,0
WS - 4A
3440 393 46000 4,7 1 1000 1,9
PM DD.
4,7
WS - 4B
4,7 1 700 1,5
4,7
WindSmile
200W vawt 200 131 15 2200 1 1 50 1,3
Al e ao.
PM DD DC.
Darrieus com 4 ps.
1,0
1kW vawt 1000 438 50 6679 1,6 1,5 60 2,5 20
PM DD.
2,4
5kW vawt 5000 1752 200 22554 2,5 3,2 700 3 20 Darrieus com 3 ps. 8,0
1 - Preo apenas da turbina. 3 - PM-Permanent magnets, DD-Direct drive, DC-Direct current, B-Sem escovas (brushless).
2 - Preo da turbina em Portugal. 4 - FV- Compsito de fibra de vidro, FC-Compsito de fibra de carbono, Al-Alumnio
91
92
Anexo 2 Clculos
t
dt
d
o e e
e
o + = =
0
dt
du
e =
2
0 0 0
2
1
t t t
dt
d
o e u u o e
u
+ + = + =
o I M =
I
Mt
2
2
= u
I
rMt
l r l
2
2
= = u
t
M M = Pr
) (Pr
2 2
Pr
2
) (Pr
2
2
t
t
t
M r
lI
t
rt
lI
M
I
t M r
l
= =
=
disco tubo a
I I I 2 + =
2
1
mr I
tubo
=
( )
| |. .
2
2
2
1
2
2
m kg
r r m
I
disco
=
sendo m a massa,
1
r o raio interior e
2
r o raio exterior. Para kg m
tubo
320 , 2 = , m r 048 , 0
1
= ,
m r 15 , 0
2
= e kg m
disco
275 , 0 = obtm-se
2
. 00534528 , 0 m kg I
tubo
=
2
. 00277695 , 0 m kg I
disco
= .
Com estes valores obtm-se,
2
. 01089918 , 0 m kg I
a
=
.
t
I
M t
I
M
0
0
e
e e = + =
93
Anexo 3 Curvas de potncia obtidas por CFD
Quadro resumo
Configurao Perfil
Corda
[mm]
Nmero
de ps
Factor de
bloqueamento
Tempo de
clculo
para 0,1 s
[h]
Nmero
de
condies
obtidas
Tempo de
clculo
total [dia]
Nmero
de clulas
da malha
1
NACA0018
50
3 0,3 4 24 28,4 75000
2 2 0,2 2,6 16 9,8 55920
3 4 0,4 5 12 13,4 98484
4
100
2 0,4 2,2 21
2
20,5 46252
5 3 0,6 3 20 37,4 57138
6
50
5 0,5 5,8 15 18,8 90600
7 NACA0012 3 0,3 4,4 21 17,6 75049
Para qualquer uma destas configuraes, o dimetro da turbina de 1 m e o poste central tem um
dimetro de 10 cm.
Configurao 1
Perfil NACA0018, corda de 50mm, 3 ps e factor de bloqueamento de 0,3.
P
C
2 3 4 4,5 5 6 7
V
2
-0,239 -0,365
-0,556 -0,258 -0,400
3
-0,035
4
0,358
5 -0,069 -0,010 0,251 0,309 0,219 -0,024
8 -0,021 0,043 0,241 0,257 0,168 -0,106
11 -0,014 0,071 0,275 0,269 0,217
M
2 3 4 4,5 5 6 7
V
2
-0,097 -0,112
-0,136 -0,053 -0,070
3
-0,021
4
0,389
5 -0,264 -0,025 0,481 0,525 0,335 -0,030
8 -0,211 0,282 1,181 1,12 0,659 -0,345
11 -0,254 0,873 2,55 2,213 1,607
2
Das quais 3 so a mesma condio mas com medies diferentes. Ou seja, o nmero de pontos no grfico
ser 18.
94
Configurao 2
Perfil NACA0018, corda de 50mm, 2 ps e factor de bloqueamento de 0,2.
P
C
3 4 4,5 5 6
V
2
-0,385
3
0,147
4
0,101
0,317 0,194
5 -0,010 0,179 0,245 0,251 0,015
8 0,014 0,198 0,256 0,238 0,14
11
0,269
M
3 4 4,5 5 6
V
2
-0,094
3
0,081
4
0,124
0,311 0,159
5 -0,037 0,343 0,418 0,384 0,019
8 0,094 0,969 1,116 0,932 0,458
11
2,217
Configurao 3
Perfil NACA0018, corda de 50mm, 4 ps e factor de bloqueamento de 0,4.
P
C
3 4 4,5 5 6
V
3
0,021
4
0,425
5 0,010 0,296 0,331 0,226 -0,143
8 0,083 0,261 0,259 0,214
11
0,314
M
3 4 4,5 5 6
V
3
0,013
4
0,462
5 0,027 0,567 0,563 0,346 -0,189
8 0,543 1,278 1,129 0,838
11
2,905
95
Configurao 4
Perfil NACA0018, corda de 100mm, 2 ps e factor de bloqueamento de 0,4.
P
C
0,3 0,4 2 3 3,5 4 5 6
V
2
0,068
3
0,244
4
0,353
5 0,027 0,278 0,366 0,343 0,168 -0,197
8 0,001 -0,001 0,056 0,293 0,358 0,342 0,178 -0,142
11
0,383
M
0,3 0,4 2 3 3,5 4 5 6
V
2
0,024
3
0,193
4
0,494
5 0,105 0,709 0,801 0,657 0,256 -0,251
8 0,083 -0,051 0,548 1,913 2,006 1,676 0,699 -0,466
11
4,059
Configurao 5
Perfil NACA0018, corda de 100mm, 3 ps e factor de bloqueamento de 0,6.
P
C
0,4 0,55 1 1,5 2 3 3,5 4 5
V
2
0,176
3
0,257
4
0,308
5
0,069 0,373 0,339 0,242 -0,120
6 0,001
7 0,004
8 0,001 -0,005 -0,018 0,064 0,093 0,334 0,328 0,225 -0,125
11
0,361
M
0,4 0,55 1 1,5 2 3 3,5 4 5
V
2
0,062
3
0,203
4
0,504
5
0,265 0,951 0,742 0,463 -0,184
6 0,026
7 0,151
8 0,060 -0,191 -0,358 0,832 0,91 2,185 1,835 1,104 -0,49
11
3,823
96
Configurao 6
Perfil NACA0018, corda de 50mm, 5 ps e factor de bloqueamento de 0,5.
P
C
2 3 3,5 4 4,5 5
V
3
-0,144
4
5
0,049 0,199 0,245 0,152 0,000
8
0,100 0,228 0,217
-0,032
11 -0,019 0,129 0,258 0,245
0,011
M
2 3 3,5 4 4,5 5
V
3
-0,099
4
5
0,124 0,436 0,469 0,258 -0,005
8
0,655 1,276 1,062
-0,126
11 -0,357 1,589 2,732 2,267
0,080
Configurao 7
Perfil NACA0012, corda de 50mm, 3 ps e factor de bloqueamento de 0,3.
P
C
3 4 4,5 5 6 7
V
2
-0,560
3
0,254
4
0,313
5 -0,002 0,249 0,235 0,32 0,135 -0,185
8 0,020 0,189 0,294 0,317 0,146 -0,137
11 0,031 0,22 0,333 0,340 0,168 -0,079
M
3 4 4,5 5 6 7
V
2
-0,137
3
0,140
4
0,306
5 -0,006 0,476 0,401 0,491 0,172 -0,203
8 0,129 0,927 1,280 1,242 0,477 -0,384
11 0,381 2,034 2,743 2,523 1,040 -0,420
97
Anexo 4 Anlise de resultados experimentais
Cada medida efectuada no corresponde a um nmero mas sim a um intervalo. Tal ocorre devido
ao processo de medio, ao erro do medidor e incerteza dos instrumentos utilizados. Quando se
efectuam clculos com essas medies, a sua incerteza afecta o resultado das operaes. Chama-se a
este processo propagao de erro. Considere-se uma varivel A, funo f das variveis X e Y .
Ambos X e Y so medies que correspondem ao intervalo
x
x X c + = e
y
y Y c + = ,
respectivamente, sendo c o erro associado a cada varivel. Neste caso, ) , ( Y X f A = sendo
a
a A c + = e
y x a
y
f
x
f
c c c
c
c
+
c
c
= . De acordo com a frmula anterior, obtm-se as seguintes
frmulas para as principais operaes:
( )
b a b a
b a b a c c c c + + = + ) ( ) (
( )
b a b a
b a b a c c c c + = ) ( ) (
( )
b a b a
a b b a b a c c c c + = . ) ).( (
|
.
|
\
|
+
|
.
|
\
|
=
b
a
b
a
b
a
b b
a
b
a
c
c
c
c
2
) (
) (
( )
a a
a a a c c 2 ) (
2 2
=
( )
a a
a a a c c
2 3 3
3 ) ( =
Procedeu-se tambm aos clculos tendo em ateno os algarismos significativos. No caso da adio
e subtrao conserva-se um nmero de casas decimais igual ao da parcela com menor nmero de
casas decimais. Nas multiplicaes e divises, o resultado ter o mesmo nmero de algarismos
significativos do factor ou divisor com menor nmero de algarismos significativos. Para exemplificar
o tratamento de dados efectuado, apresenta-se em detalhe os clculos de um dos ensaios
experimentais para a configurao 5:
Dados recolhidos
Tempo de durao do ensaio
1 , 0 0 , 10 = t s (cronmetro com mostrador digital: erro=resoluo)
Leitura inicial e final do anemmetro
5 870 =
i
Pos m, (anemmetro com mostrador analgico: erro=metade da resoluo)
5 945 =
f
Pos m,
Velocidade linear do sensor do taqumetro
1 , 0 5 , 3 = Vel km/h, (taqumetro com mostrador digital: erro=resoluo)
Raio do sensor do taqumetro
0005 , 0 10 6 , 1
2
1
= d m, (rgua graduada: erro=metade da resoluo)
Raio das ps da turbina
0005 , 0 500 , 0 = r
98
Clculos
Obteno da velocidade do vento
1 , 1 5 , 7
1 , 0 0 , 10
10 75
=
t
Pos Pos
V
i f
m/s
Converso da velocidade linear do sensor do taqumetro para m/s
s m h km Vel / 03 , 0 97 , 0 / 1 , 0 5 , 3 = =
Obteno da velocidade angular da turbina
2 , 0 1 , 6
0005 , 0 160 , 0
03 , 0 97 , 0
1
=
= =
d
Vel
e rad/s
Clculo da velocidade linear das ps da turbina
1 , 0 1 , 3 ) 0005 , 0 500 , 0 ( ) 2 , 0 1 , 6 ( = = = r V
p
e m/s
Clculo do coeficiente de velocidade perifrica
07 , 0 41 , 0
1 , 1 5 , 7
1 , 0 1 , 3
=
= =
V
V
p
Resumo dos resultados obtidos para a configurao 5
Ensaio 36 39 49 54 55 58 60
0,410,07 0,370,07 0,350,06 0,480,08 0,440,07 0,440,07 0,590,09
073 , 0 440 , 0 =
Resumo dos resultados obtidos para a configurao 4
Ensaio 1 2 4 6 7 8 10
0,310,04 0,330,04 0,310,04 0,400,05 0,400,05 0,340,04 0,330,04
043 , 0 346 , 0 =