Construção, Desenvolvimento e Otimização de Concentrador e Difusor para Turbina Eólica - PDF
Construção, Desenvolvimento e Otimização de Concentrador e Difusor para Turbina Eólica
Download "Construção, Desenvolvimento e Otimização de Concentrador e Difusor para Turbina Eólica"
Nathan Amaral Salazar
1 Construção, Desenvolvimento e Otimização de Concentrador e Difusor para Turbina Eólica Jorge Alexandre Rente Paulo Relatório Final do Trabalho de Projeto apresentado à Escola Superior de Tecnologia e de Gestão Instituto Politécnico de Bragança para obtenção do grau de Mestre em Engenharia Industrial Ramo de Engenharia Eletrotécnica Julho de 2013
3 Construção, Desenvolvimento e Otimização de Concentrador e Difusor para Turbina Eólica Jorge Alexandre Rente Paulo Relatório Final do Trabalho de Projeto apresentado à Escola Superior de Tecnologia e de Gestão Instituto Politécnico de Bragança para obtenção do grau de Mestre em Engenharia Industrial Ramo de Engenharia Eletrotécnica Orientadores: Professor Doutor Luís Frölén Ribeiro Professor Doutor Valdemar Garcia Julho de 2013
5 Dedicatória À fundação Casa de Trabalho. i
7 Agradecimentos Gostaria de agradecer a todas as pessoas que estiveram presentes e que me apoiaram. Um agradecimento especial: ao Professor Doutor Luís Frölén Ribeiro pela orientação e ensinamentos concedidos ao longo do desenvolvimento deste projeto, assim como a dedicação e o incentivo; ao Professor Doutor Valdemar Garcia pela orientação, pelos conhecimentos transmitidos, pela dedicação e paciência na elaboração de todo o relatório; ao grande companheiro de trabalho Duarte Freitas, pelos conhecimentos, pela companhia e sobretudo pela amizade que construímos ao longo destes últimos anos; aos técnicos de laboratório da ESTiG, pela ajuda prestada e pela disponibilidade; ao Professor Doutor Miguel Paula pelos conhecimentos transmitidos, pela sua amizade e pela sua disponibilidade; ao Professor Doutor Vicente Leite pelos conhecimentos transmitidos e pela sua disponibilidade; à ESTiG pelo material nanciado e pelas facilidades concedidas para a realização deste projeto; aos meus pais e irmãos, pelo incentivo, paciência e pela ajuda concedida; à fundação Casa de Trabalho, pelo o apoio, pela motivação e pela educação nos 19 anos que estive acolhido, o que me permitiu acabar mais uma etapa do meu percurso. iii
8 iv _
9 v _ Resumo _ O objetivo principal deste trabalho consistiu na construção, desenvolvimento e otimização de um dispositivo de captação e aceleração do vento para aumentar a potência numa turbina eólica. Efetuou-se um estudo numérico e experimental com turbinas eólicas de eixo horizontal usando-se apenas turbinas e sistemas constituídos pelo conjunto turbina e dispositivo de captação e aceleração do vento. O dispositivo usado foi do tipo concentrador-difusor, C-D. O trabalho desenvolveu-se em 3 etapas: primeiro foram efetuadas simulações numéricas, seguidas por experiências com um modelo reduzido de turbina em túnel de vento e por m concluídas em experiências de campo. As turbinas usadas neste estudo eram constituídas por 6 pás, sendo o diâmetro de m na turbina usada nas simulações e no túnel de vento e de 0.91 m nas experiências de campo. Nas simulações numéricas e nos ensaios no túnel de vento a velocidade de aproximação do vento variou entre 5.7 e os 10 m/s. Com base nas simulações foram construídos dois modelos de C-D. Os resultados obtidos, quer nas simulações numéricas, quer nas experiências em túnel de vento, mostraram que a potência extraída foi sempre superior quando se utilizou o sistema turbina e C-D. Das simulações numéricas vericou-se um aumento máximo da velocidade do vento junto da turbina de 21% sendo o aumento médio de 15%. O aumento máximo na potência foi de 112% e o médio foi de 71%. Nos ensaios efetuados no túnel de vento vericaram-se aumentos na potência elétrica extraída para os dois dispositivos construídos: aumento máximo da potência foi de 153% para um aumento médio de 107%. No C-D 2 os aumentos foram de 187% e 124%, respetivamente. Nas experiências de campo foram apenas obtidos resultados preliminares da potência em função da velocidade do vento. Palavras-chave: Turbina eólica de eixo horizontal, concentrador-difusor
11 Abstract The main objective of this thesis was the construction, development and optimization of wind acceleration device to increase the output power of a wind turbine. Numerical and experimental studies of horizontal axis wind turbine were conducted without and with the acceleration device. The device was a concentrator-diuser (C-D). The work was developed in three stages: rstly numerical simulations were performed, followed by wind tunnel measurements of a scaled model turbine in the wind tunnel and nally eld experiments were carried out. A 6 bladed wind turbine was used in this study: m the diameter for the wind tunnel and 0.91 m in eld experiments. The ow eld velocity ranged from 5.7 and 10 m/s for the numerical simulations and for the experimental set-up. Numerical simulations yield information on a single turbine and with the turbine equipped with C-D and allowed to build 2 optimized C-Ds, latter installed in the wind tunnel. The power output results with the C-D were always larger with de C-D devise. For the numerical simulations the maximum, wind velocity increase was 21%, while the average increase was 15% when comparing to the no C-D case. The increase in maximum power was 112% and the average increase was 71%. The wind tunnel showed an increase in the electric power output for both devices built: for C-D 1 the increase of maximum power was 153% and the average increase of 107%; for C-D 2 the increase was 187% and 124%, respectively. Only preliminary results of the wind turbine power curve were obtained in the eld experiments. Keywords: Horizontal axis wind turbine, concentrator-diuser vii
13 Conteúdo Lista de Figuras xi Lista de Tabelas xv Lista de Abreviaturas xvii 1 Introdução Objetivos Enquadramento Implantação da tecnologia eólica no mundo Estrutura Teoria e Revisão Bibliográfica Energia eólica Teoria da energia eólica Dispositivos de captação e aceleração do vento Aspetos diferenciadores Materiais e Métodos Simulações numéricas Ensaios no túnel de vento Experiências de campo Resultados e Discussão Simulações numéricas Velocidade média do vento junto da turbina Resultados da potência Ensaios no túnel de vento Experiências de campo ix
14 x CONTEÚDO 5 Conclusões Conclusões Trabalhos futuros Bibliograa 49 Anexos 55 A Dimensionamento do poste 55 B Calibração do túnel de vento 59
15 Lista de Figuras Potência eólica instalada acumulada por ano no mundo [1] Potência eólica instalada acumulada na Europa [1] Representação esquemática de turbinas de eixo horizontal: (a) barlavento; (b) sotavento. Adaptado de [2] Representação esquemática de três turbinas de eixo vertical: (a) darrieus; (b) savonius; (c) giromill. Adaptado de [3] Esquema de análise do comportamento do vento numa turbina: (a) escoamento do vento; (b) distribuição de velocidade; (c) distribuição de pressão. Adaptado de [4] Cp em função TSR. Adaptado de [3] Exemplo de um perl de pá, modelo NACA [5] Representação 2D do perl NACA 8311 com as respetivas incógnitas de dimensionamento Representação esquemática de vários sistemas para captar e acelerar o vento: (a) concentrador; (b) tubo cilíndrico; (c) difusor. Adaptado de [6] Modelo 3D de dispositivo com turbina [7] Modelos de campo, turbina e dispositivo com turbina. Adaptado de [8] Modelos de campo, turbina e dispositivo com turbina em testes. Adaptado de [9] Modelo de campo, turbina e difusor [10] Modelo de estudo, turbina e C-D [11] Modelo de estudo, C-D [12] Sistema de captação e aceleração da torre Bahrain World Trade Center [13] xi
16 xii LISTA DE FIGURAS Representação esquemática de um sistema constituído por turbina, WT, e concentrador-difusor, C-D (a) Turbina Rutland 913; (b) gerador elétrico Distribuição da velocidade do vento em redor do perl NACA 8311, vista 2D Representação 3D do modelo de WT dimensionada Dispositivo C-D: (a) representação esquemática com a respetiva legenda; (b) vista de corte 3D Malha em 3D utilizada nas simulações Volume do cilindro junto à WT utilizado nas simulações para obter U m Túnel de vento usado nos ensaios Tubo de Pitot para medição da velocidade do vento Calibração do túnel de vento: velocidade do ar em função da frequência do ventilador Sequência de algumas fases da maquinagem e do acabamento das pás Turbina eólica utilizada nas experiências no túnel de vento Algumas etapas da construção dos C-D C-D 1 utilizado nas experiências no túnel de vento Representação 3D das partes do poste: (a) sapata; (b) união dos tubos; (c) xação da turbina Fases de construção inicial do poste: (a) furação das várias chapas; (b) fases de limpeza dos tubos e soldagem das chapas aos tubos; (c) pintura dos tubos e acessórios Fases da instalação da turbina e do poste no terreno da Quinta do Poulão: (a) xação da sapata; (b) xação da turbina e dos sensores; (c) xação de um tubo auxiliar para erguer o poste Turbina e poste xos no campo: (a) turbina e sensores; (b) sistema global instalado Sensores para análise do vento: (a) sensor de velocidade; (b) sensor de direção
17 LISTA DE FIGURAS xiii Representação esquemática do diagrama do circuito de controlo da potência fornecida à carga. Adaptado de [14] Visualização da carga xa ao poste, utilizada para dissipar a energia elétrica Sistema de aquisição de dados e sistema de controlo de potência: (a) representação esquemática do circuito; (b) placa impressa com os componentes integrados Distribuição de velocidades do vento em redor da turbina: (a) vista corte-trás; (b) vista corte lateral Distribuição de velocidades do vento em redor do sistema turbina e C-D 1: (a) vista corte-trás; (b) vista corte lateral Velocidade média adimensional do vento para várias simulações da primeira fase Velocidade média adimensional do vento em função da velocidade do vento para as várias simulações da segunda fase Potência adimensional para as várias simulações da primeira fase Potência adimensional em função da velocidade do vento para as várias simulações da segunda fase Potência adimensional em função da velocidade do vento para os ensaios da primeira fase no túnel de vento Potência adimensional em função da velocidade do vento para os ensaios da segunda fase no túnel de vento Potência em função da velocidade do vento para várias simulações numéricas Potência em função da velocidade do vento para vários ensaios no túnel de vento Potência média em função da velocidade do vento das experiências de campo da primeira fase Potência média em função da velocidade do vento das experiências de campo da segunda fase A.0.1Esquema estrutural do caso em estudo
18 xiv LISTA DE FIGURAS
19 Lista de Tabelas 1.1 Potência acumulada instalada nalguns países da Europa no nal de 2012 [1] Percentagem de energia elétrica total consumida na Europa em 2010, proveniente de energia eólica [15] Rugosidades típicas de vários tipos de superfícies do solo [16] Valores dos parâmetros resultantes do dimensionamento da turbina eólica utilizada Condições das simulações só com turbina e para sistema turbina e C-D 1 e C-D Código das simulações e condições geométricas dos C-D Código e condições experimentais para os ensaios no túnel de vento, quando se utilizou apenas a turbina e quando se utilizou o sistema turbina e C-D 1 e C-D xv
20 xvi LISTA DE TABELAS
21 Lista de Abreviaturas AWEA American wind energy association C-D CNC MDF MPP P-O SD TSR WT Concentrador-difusor Computer numerical control Medium density berboard Maximum power point Perturbação-observação Secure digital Tip speed ratio Wind turbine / Turbina eólica xvii
22 xviii LISTA DE ABREVIATURAS -
23 Capítulo 1 Introdução 1.1 Objetivos Neste trabalho de projeto de Mestrado em Engenharia Industrial pretendese projetar, desenvolver e otimizar um dispositivo de captação e aceleração do vento do tipo tubeira convergente-divergente de modo a que, aumentando a velocidade do vento, se possa melhorar a eciência de turbinas eólicas de pequenas dimensões. 1.2 Enquadramento O aumento do consumo de energia a nível nacional e mundial, ligado aos constantes aumentos na fatura energética e aos problemas ambientais associados à produção e utilização da energia de origem fóssil e nuclear, tem conduzido a um aumento exponencial da procura de energia a partir de fontes renováveis e pouco poluentes. Enquadra-se aqui, o aumento da procura da energia eólica, com o aumento da potência instalada em parques eólicos, em simultâneo com o aumento da investigação nesta área, de modo a melhorar a eciência na produção desta forma de energia. Este trabalho vem na sequência do projeto de licenciatura em Engenharia de Energias Renováveis desenvolvido pelos alunos Jorge Paulo e Duarte Freitas. Pretende-se cumprir os objetivos a seguir indicados: - Construção de modelos de turbinas eólicas de eixo horizontal, com 3 ou mais hélices e pers metálicos (ou de outro material) do tipo de tubo 1
24 2 1.2 Enquadramento convergente e divergente para que, quando colocados em turbinas, permitam aumentar a velocidade do vento e assim aumentar a eciência da turbina. - Realização de testes de simulação em computador para as turbinas e pers, recorrendo a programa comercial, com o objetivo de determinar variáveis como velocidades e pressões. - Teste experimental dos modelos de turbinas e pers, no túnel de vento existente no Laboratório de Mecânica de Fluidos e Hidráulica, obtendo-se valores de variáveis também obtidas na simulação em computador e outros parâmetros, como a tensão, corrente e a velocidade angular. - Construção de turbina eólica e perl, numa escala maior, semelhante à turbina testada em modelo reduzido e realização de testes experimentais. Neste relatório apresentam-se resultados preliminares obtidos apenas com a turbina. No futuro imediato serão realizadas experiências também com o dispositivo convergente-divergente de modo a obter resultados passíveis de ser comparados com os já obtidos nas simulações e nas experiências efetuadas em túnel de vento com o modelo reduzido de turbina e concentradordifusor. Este trabalho foi apresentado na conferência Windpower conference & exhibition 2013 realizado pela Associação de Energia Eólica da América (AWEA), que decorreu entre o dia 5 e 8 de Maio em Chicago, tendo sido apresentado em formato de poster e artigo [17] Implantação da tecnologia eólica no mundo Atendendo aos vários problemas associados à energia de origem fóssil, como o petróleo, como sejam o preço, a escassez e a poluição, os países voltam-se cada vez mais para fontes de energias renováveis e não poluentes. Destas, a energia proveniente do vento, energia eólica, é das que tem tido um aumento mais acentuado nos últimos anos. A potência eólica total instalada acumulada no mundo no nal de 2012 era de MW, (1135 MW em África e Médio Oriente, MW na Ásia, MW na Europa, 3505 MW na America Latina, MW na América do Norte e 3219 MW na Região do Pacíco) [1]. A potência eólica instalada acumulada a nível mundial tem vindo a subir bastante nos últimos anos. Nos últimos 12 anos passou de MW para
25 1.2 Enquadramento MW em 2012, como se pode observar na Figura Figura 1.2.1: Potência eólica instalada acumulada por ano no mundo [1]. Na Europa também se tem registado um aumento gradual da potência instalada acumulada, como se pode ver na Figura para os últimos 12 anos. Figura 1.2.2: Potência eólica instalada acumulada na Europa [1]. Na Tabela 1.1 constam os países da Europa que mais têm apostado na energia eólica, tendo Portugal em 2012 uma potência instalada acumulada de 4525 MW [1].
26 4 1.2 Enquadramento Tabela 1.1: Potência acumulada instalada nalguns países da Europa no nal de 2012 [1]. País Potência acumulada (MW) Alemanha Espanha Reino Unido 8445 Portugal 4525 Dinamarca 4162 A maior parte dos parques eólicos está instalada em terra (em inglês onshore). Contudo, dada a maior velocidade do vento, a instalação de parques eólicos no mar (em inglês oshore) tem vindo a aumentar a nível mundial. A título de exemplo, em Portugal a velocidade do vento varia entre 3.5 e 4.5 m/s em terra e entre 4.5 e 5.0 m/s no litoral, podendo chegar aos 6 a 7 m/s no mar [18]. A potência instalada acumulada no mar era de 5415 MW no nal de 2012 no mundo, tendo Portugal 2 MW instalados [1]. Do total da energia elétrica consumida em Portugal em 2010, 17% era produzida a partir de turbinas eólicas, o que garantia ao nosso país uma posição de destaque na Europa [15], ver Tabela 1.2. Tabela 1.2: Percentagem de energia elétrica total consumida na Europa em 2010, proveniente de energia eólica [15]. País Energia Produzida (%) Dinamarca 27 Portugal 17 Espanha 16 Alemanha 11 Reino Unido 6 França 3 A partir de 2003, Portugal apostou bastante na instalação de parques eólicos. No nal desse ano, a potência instalada de energia eólica era de 262 MW, em 2007 era de 2000 MW, em 2009 era de 3400 MW e em 2012 era de 4194 MW [19]. Em termos percentuais, a potência instalada vinda de parques eólicos
27 1.3 Estrutura 5 atingiu 21% em 2010, 22% em 2011 e 23% em 2012 da potência total instalada a nível nacional, representando 58% da potência total instalada proveniente de fontes renováveis [19]. A produção de eletricidade através deste recurso atingiu 17% em 2010, 18% em 2011 e 20% em 2012 da produção total [19]. 1.3 Estrutura Este relatório está estruturado em 5 capítulos. No capítulo 1 apresenta-se a introdução, referindo-se os objetivos deste trabalho, o enquadramento e apresenta-se a estrutura do trabalho. No capítulo 2 apresenta-se a revisão bibliográca. Neste capítulo referese a energia eólica, a teoria da energia eólica e apresenta-se um resumo dos dispositivos de captação e aceleração do vento. No capítulo 3 referem-se os materiais e métodos. Neste capítulo explicamse as simulações numéricas, os ensaios no túnel de vento e as experiências de campo. No capítulo 4 são apresentados os resultados e a sua discussão. Neste capítulo apresentam-se os resultados das simulações numéricas, dos ensaios no túnel de vento e das experiências de campo, velocidade média do vento e da potência. São comparados os resultados das várias fases. Por m, no capítulo 5 são apresentadas as conclusões e indicados os trabalhos futuros.
28 6 1.3 Estrutura
29 Capítulo 2 Teoria e Revisão Bibliográfica 2.1 Energia eólica Segundo os historiadores, já no ano 1700 a.c. tinham surgido as primeiras máquinas eólicas, no Oriente, que utilizavam a energia do vento para irrigação dos campos. Existem registos da utilização de energia eólica no Afeganistão no séc. VIII d.c. [20]. Na idade média eram usados moinhos de vento para moer o cereal e também no bombeamento de água para irrigação [21]. No início do século XX eram usados turbinas eólicas com 20 pás de baixa velocidade e com um elevado binário para bombear água [21]. Só nos anos 80 do século XX é que se passou a utilizar turbinas eólicas em massa para produção de energia elétrica distribuídas em parques eólicos [22]. O aproveitamento do vento para a produção de energia elétrica é feito por vários tipos de turbinas, dos quais se distinguem dois grupos principais: as de eixo horizontal, as mais comuns, e as de eixo vertical [2, 16, 23]. Nas turbinas de eixo horizontal podem distinguir-se dois tipos conforme a posição do rotor: a barlavento (em inglês upwind), quando o rotor está posicionado para o vento de frente e a sotavento (em inglês downwind), quando o rotor está posicionado com a parte de trás do rotor para o vento, conforme se ilustra na Figura [2, 23]. 7
30 8 2.1 Energia eólica (a) (b) Figura 2.1.1: Representação esquemática de turbinas de eixo horizontal: (a) barlavento; (b) sotavento. Adaptado de [2]. De entre os vários tipos de turbina de eixo vertical apresentam-se na Figura três das congurações mais comuns: darrieus, savonius e giromill [2, 23]. (a) (b) (c) Figura 2.1.2: Representação esquemática de três turbinas de eixo vertical: (a) darrieus; (b) savonius; (c) giromill. Adaptado de [3]. Qualquer tipo de turbina pode ter um número variado de pás acopladas no rotor, estando este número relacionado com vários aspetos, como por
31 2.1 Energia eólica 9 exemplo a sua potência e o local da instalação. Quanto à potência da turbina distinguem-se três grupos: o primeiro é designado de micro turbinas, até 1 kw de potência [24, 25]; o segundo é designado por pequenas turbinas (em Inglês, small turbine), de 1 a 100 kw [4, 24, 26, 27, 28] e por último temos as grandes turbinas (em Inglês, large turbine) com mais de 100 kw, normalmente na ordem dos MW [4, 26]. Quanto ao local, as turbinas podem estar instaladas no mar ou na terra. No mar podem estar xadas no fundo, normalmente perto da costa terrestre, ou podem ser utuantes, normalmente mais afastadas da costa e amarradas com cabos de aço ao fundo do mar [29]. Na terra podem estar instaladas no meio urbano ou afastadas deste. Como a produção de energia eólica aumenta com o aumento da velocidade do vento, a maioria das grandes turbinas estão instaladas nas montanhas altas, onde a velocidade é mais elevada. Normalmente, as micro e as pequenas turbinas estão instaladas em meios urbanos e também em locais remotos. Outro fator importante que inuencia a velocidade do vento é a rugosidade do solo. Na Tabela 2.1 apresentam-se alguns valores típicos consoante o tipo de superfície [16]. Tabela 2.1: Rugosidades típicas de vários tipos de superfícies do solo [16]. Descrição Rugosidade (m) Mar aberto Superfície com lama Pasto em terreno aberto 0.03 Terrenos agrícolas de produção baixa 0.1 Terrenos agrícolas de produção alta 0.25 Parques de estacionamento com árvores e arbustos 0.5 Obstáculos regulares (orestas, aldeias, bairros) 1 Centro das grandes cidades 2 Normalmente a velocidade do vento é medida a 10 m do solo [16]. Dado que a maioria das turbinas eólicas não se encontram a esta altura, o potencial eólico, isto é, a velocidade para uma determinada altura calcula-se através da seguinte equação [16]: u 2 = u 1. ln( h 2 d Z 0 ) ln( h 1 d Z 0 ) (2.1)
32 Teoria da energia eólica onde u 2 é a velocidade do vento para uma determinada altura h 2, u 1 é a velocidade do vento para a altura h 1 conhecida, Z 0 é o valor da rugosidade da superfície (ver Tabela 2.1) e d é o deslocamento da camada limite causado pelo obstáculo, que normalmente é zero [16]. Como nos meios urbanos existe uma maior rugosidade na superfície do solo, esta faz com que haja uma maior turbulência do vento fazendo baixar a sua velocidade [30]. O resultado é vericado numa perda na produção de energia elétrica. Com o objetivo de aumentar a energia elétrica produzida nos meios urbanos, vários autores [31, 32, 33, 34, 35, 36, 37, 38, 39, 40, 41, 42, 43, 44, 45] realizaram estudos de novos modelos de turbinas, estudos teóricos, estudos de viabilidade de turbinas nos meios urbanos e estudos de simulações numéricas de vários tipos de turbinas adaptadas para este m. No nal do ano de 2011 existiam cerca de pequenas turbinas no mundo, perfazendo um total de 576 MW de potência acumulada instalada [28]. Em 2009 foi registada uma produção 382 GWh com pequenas turbinas em todo o mundo [46]. No nal de Março de 2013, Portugal tinha 164 unidades de pequenas turbinas registadas, perfazendo um total de 564 kw de potência total acumulada [47]. O objetivo da nossa investigação coincide nesta temática, dado que se pretende obter um dispositivo que otimize a produção de energia elétrica nos meios urbanos. O dispositivo estudado reduz a turbulência do vento e faz aumentar a sua velocidade, aumentando assim a produção de energia elétrica. 2.2 Teoria da energia eólica A energia eólica tem como origem a energia solar. O movimento de massas de ar causadas pelas diferenças de temperatura na Terra faz com que haja este potencial [4]. De maneira a quanticar esse potencial, e portanto, a potência da turbina, é necessário conhecer o comportamento de algumas variáveis aerodinâmicas intervenientes no processo, como sejam a velocidade do vento e a pressão em torno da turbina. Na Figura representa-se esquematicamente a turbina e a corrente de escoamento (a), a distribuição de velocidade
33 2.2 Teoria da energia eólica 11 (b) e a distribuição de pressão (c). (a) (b) (c) Figura 2.2.1: Esquema de análise do comportamento do vento numa turbina: (a) escoamento do vento; (b) distribuição de velocidade; (c) distribuição de pressão. Adaptado de [4]. Tendo como referência a Figura e como apresentado em [4], antes do escoamento chegar à turbina, o vento tem velocidade U e uma pressão p 1. Junto da turbina, na parte da frente (secção 2), os valores correspondentes são U 2 e p 2. Quando passa da secção 2 para a 3, o ar sofre uma queda de pressão brusca, passando de p 2 para p 3, no entanto a velocidade mantémse, U 2 = U 3. A pressão aumenta depois gradualmente até à secção 4, onde p 4 = p 1. À medida que a pressão vai aumentando, de p 3 para p 4, a velocidade do vento vai diminuindo gradualmente, de U 3 a U 4. A energia que o turbina eólica (rotor) recebe do ar corresponde à queda de pressão nela registada, entre as secções 2 e 3, (Δp = p 3 p 2 ). Através da conservação de quantidade de movimento pode explicitar-se a força F v que o ar comunica ao rotor pode obter-se a partir da equação: F v = AUρ ar (U U 4 ) (2.2) em que A representa a área do rotor, ρ ar a massa volúmica do ar e U = U 2 = U 3 é a velocidade do ar ao nível do rotor. A força F v é igual e diretamente oposta à que o rotor exerce no ar, e que, sendo U 3 = U 2 tem por módulo A(p2 p3). Deste modo, obtém-se p 2 p 3 = ρ ar U(U U 4 ) (2.3) Por aplicação da equação de Bernouli entre 1 e 2 e 3 e 4, e atendendo a que U 2 = U 3 e que p 4 = p 1 obtém-se
34 Teoria da energia eólica p 2 p 3 = 1 2 ρ ar(u 2 U 2 4 ) (2.4) A comparação das equações 2.3 e 2.4 permite concluir que U na equação 2.2 é dado por U = U + U 4 2 (2.5) A equação 2.5 mostra que a velocidade no rotor é a média das velocidades da secção 1 e da secção 4. A potência útil da turbina, P u, corresponde à energia cinética comunicada ao rotor pelo ar que tem um caudal volúmico Q, é dador por P u = 1 2 ρ arq(u 2 U 2 4 ) (2.6) Se o rotor não existisse, a velocidade seria constante, U, atravessando a área A, resultando na potência disponível (P d ) P d = 1 2 ρ ar(au )U 2 (2.7) Deste modo, o rendimento teórico do rotor, Cp max, é denido como a razão entre P u e P d : Cp max = P u P d = 1 ρ 2 arq(u 2 U 2 4 ) 1 ρ 2 ar(au )U 2 = AU (U 2 U 2 4 ) AU 3 = = (U + U 4 )(U 2 U 2 4 ) 2U 3 = (1 + γ)(1 γ2 ) 2 (2.8) em que γ = U 4 U. A equação 2.8 tem um valor máximo para γ = 1 3, sendo o resultado conhecido como o limite de Betz: Cp max = ( )(1 1 9 ) 2 = 59.3% (2.9) O rendimento real, Cp, de uma turbina eólica depende de várias variáveis como por exemplo, o tipo de turbina, a geometria e o número de pás e a razão entre a velocidade tangencial da ponta da pá e a velocidade do vento, conhecida pela sigla TSR, tip speed ratio. Pode ser obtido a partir de um gráco como o da Figura
35 2.2 Teoria da energia eólica 13 Figura 2.2.2: Cp em função TSR. Adaptado de [3]. Para dimensionar a geometria pretendida das pás de uma turbina é necessário proceder-se ao cálculo de um conjunto de variáveis, conforme as equações de 2.10 a 2.17, a seguir apresentadas. As variáveis estão representadas na Figura e na Figura A velocidade tangencial na ponta da pá, v R, é dada por onde λ corresponde a TSR. v R = U λ (2.10) Através de λ obtém-se a velocidade angular ótima (Ω) da turbina pela equação 2.11 onde R é o raio do rotor. Ω = 30λU πr (2.11) No dimensionamento de uma turbina eólica procede-se à divisão do R em 10 ou 20 secções, designadas de r locais. Para cada r é necessário obter o TSR local (λ r ), obtendo-se pela manipulação da equação 2.11
36 Teoria da energia eólica λ r = Ωπr 30U (2.12) Através de λ r calcula-se a velocidade tangencial local, v r. A partir da equação 2.10 obtém-se v r = U λ r (2.13) Através de uma análise trigonométrica, obtém-se o ângulo local, ϕ r1, da velocidade resultante: ϕ r1 = arctan( U v r ) (2.14) O ângulo ótimo local, ϕ r, para extrair a potência ótima é dador por O ângulo de passo local, β r é dado por ϕ r = 2 3 ϕ r1 (2.15) β r = ϕ r α r (2.16) onde o α r é o ângulo de ataque local, o seu valor está relacionado com o coeciente de sustentação local, C Lr, e com o coeciente de arrasto local, C Dr, variando de acordo com o tipo de perl da pá, conforme se ilustra na Figura Figura 2.2.3: Exemplo de um perl de pá, modelo NACA [5]. O valor do α r deve ser escolhido de modo a obter-se um valor ótimo da razão C Lr C Dr.
37 2.3 Dispositivos de captação e aceleração do vento 15 Figura 2.2.4: Representação 2D do perl NACA 8311 com as respetivas incógnitas de dimensionamento. A corda local, c r, do perl local é obtida pela equação 2.17 c r = 16πr zc Lr sen 2 ( 1 3 ϕ 1) (2.17) A espessura local de cada parte da pá, e r, está relacionada com o modelo de perl escolhido, tendo uma relação xa com a corda local, c r. Existem vários tipos de pers que podem ser usados no desenho das pás [5]. Normalmente uma pá de eixo horizontal tem vários tipos de pers distribuídos ao longo do seu raio [23]. Para mais detalhes associados ao cálculo da potência e ao dimensionamento das pás, consultar os autores [16, 20, 23, 26, 27, 48, 49, 50]. 2.3 Dispositivos de captação e aceleração do vento Tem havido nas zonas urbanas de alguns países um aumento de instalações de pequenas turbinas [28, 46, 51]. Nesses locais, a velocidade do vento é normalmente inferior e a turbulência superior levando a uma menor produtividade. Existem estudos de turbinas com sistemas incorporados para aumentar a velocidade do vento. É sabido que a potência disponível (P d ) varia com a velocidade do vento (U ) elevada ao cubo, pelo que um pequeno
38 Dispositivos de captação e aceleração do vento aumento da U conduz a um aumento signicativo no valor da potência. Os dispositivos usados para captar e acelerar o vento que vai incidir numa turbina eólica (WT) são estruturas ocas normalmente de três tipos: concentrador, com o vento a ser captado por um tubo cuja secção reta diminui (ver Figura (a)); tubo cilíndrico (ver Figura (b)) e difusor em que a área da secção reta aumenta (ver Figura (c)) [6]. (a) (b) (c) Figura 2.3.1: Representação esquemática de vários sistemas para captar e acelerar o vento: (a) concentrador; (b) tubo cilíndrico; (c) difusor. Adaptado de [6]. _ Bet e Grassmann [7] realizaram simulações numéricas com um difusor e uma turbina de 0.56 m de diâmetro, conforme se mostra na Figura 2.3.2, tendo obtido um aumento de 100% na potência. Figura 2.3.2: Modelo 3D de dispositivo com turbina [7]. Mais tarde, Grassmann e colaboradores [8] realizaram ensaios de campo com duas turbinas de 0.56 m de diâmetro. Uma delas estava equipada com um difusor e a outra sem difusor (ver Figura 2.3.3). Obtiveram um aumento de 55% para velocidades elevadas e de 100% para baixas velocidades do vento.
39 2.3 Dispositivos de captação e aceleração do vento 17 Figura 2.3.3: Modelos de campo, turbina e dispositivo com turbina. Adaptado de [8]. Matsushima e colaboradores [9] realizaram simulações numéricas e construiram um protótipo com 1 m de diâmetro e 2 metros de comprimento (ver Figura 2.3.4). Figura 2.3.4: Modelos de campo, turbina e dispositivo com turbina em testes. Adaptado de [9]. Nas simulações numéricas a análise consistiu num difusor sem turbina no centro, registando um aumento máximo na velocidade do vento de 70%. Nos ensaios de campo obtiveram resultados da energia elétrica produzida com e sem difusor. O aumento máximo atingido foi de 140% e o aumento médio foi de 64%. Phillips e colaboradores [52] realizaram simulações numéricas de um
40 Dispositivos de captação e aceleração do vento difusor, já utilizado há alguns anos em experiências de campo. O rotor tem 7.3 m de diâmetro, conforme se mostra na Figura Figura 2.3.5: Modelo de campo, turbina e difusor [10]. Nos resultados numéricos, o aumento da velocidade do vento foi de 100% e o da potência foi de 300%. Contudo, referem que no campo o difusor não atingiu estes valores. _ Wang e colaboradores [11] realizaram um estudo aerodinâmico recorrendo a programas de simulações de um dispositivo do tipo C-D com uma turbina eólica. O tamanho do rotor foi de 0.91 m, conforme representado na Figura Figura 2.3.6: Modelo de estudo, turbina e C-D [11]. O aumento na velocidade do vento foi de 50% e o aumento da energia extraída foi de 125%. Setoguchi e colaboradores [12] zeram uma série de experiências com um dispositivo C-D, (ilustrado Figura 2.3.7), sem turbina com o objetivo de comprovar que existe uma variação na velocidade do vento devido à
41 2.3 Dispositivos de captação e aceleração do vento 19 alteração da parte envolvente do C-D. Obtiveram um aumento na velocidade do vento de 30%. Figura 2.3.7: Modelo de estudo, C-D [12]. Vários autores realizaram estudos de outros sistemas de captação e aceleração do vento indicados nas seguintes referências [53, 54, 55, 56, 57, 58, 59, 60]. Outros sistemas de captação podem ser incorporados nos próprios edifícios, onde o desenho do edifício vai captar e acelerar o vento. A título de exemplo, tem-se a torre Bahrain World Trade Center, como se pode ver na Figura com 3 turbinas de 225 kw cada. Figura 2.3.8: Sistema de captação e aceleração da torre Bahrain World Trade Center [13]. _
42 Dispositivos de captação e aceleração do vento Os dispositivos usados com as turbinas possibilitam a produção regular de energia elétrica em zonas de reduzida intensidade de vento durante um maior período de tempo. Também vão atenuar a turbulência do vento na turbina, levando a que a produção de eletricidade seja mais estável Aspetos diferenciadores Neste trabalho foi usado um dispositivo que combina um concentrador e um difusor (C-D) unidos na parte mais estreita, sendo a turbina colocada a meio, na parte mais estreita do dispositivo, Figura Pretende-se vericar que o valor da potência (P) obtido com o sistema turbina e C-D é superior ao obtido só com turbina. Figura 2.3.9: Representação esquemática de um sistema constituído por turbina, WT, e concentrador-difusor, C-D.
43 Capítulo 3 Materiais e Métodos Neste trabalho efetuou-se um estudo numérico e experimental de um modelo reduzido de turbina eólica de eixo horizontal equipada com um dispositivo do tipo C-D. Efetuaram-se também experiências de campo com uma turbina de eixo horizontal à escala real, de pequenas dimensões. Foram efetuadas simulações num programa de mecânica computacional, testes experimentais num túnel de vento utilizando-se um modelo de turbina e um sistema constituído por turbina e C-D e testes de campo utilizando apenas a turbina eólica. Como se pretende efetuar estudos num protótipo de C-D com uma turbina do tipo Rutland 913, com 6 pás, como se mostra na Figura (a), desenhou-se e construiu-se um modelo de turbina geometricamente semelhante. As dimensões do modelo da turbina e do C-D desenhados e posteriormente construídos, foram obtidas tendo em conta o tamanho do túnel de vento disponível para a realização dos ensaios experimentais e do tamanho do gerador elétrico, Figura (b). Como consequência o diâmetro (D) da turbina do modelo reduzido foi de m. 21
44 Simulações numéricas (a) (b) Figura 3.0.1: (a) Turbina Rutland 913; (b) gerador elétrico. 3.1 Simulações numéricas As simulações numéricas, quer da turbina quer do sistema turbina e C-D foram efetuadas no software comercial SolidWorks com o pack de mecânica de uidos computacional. Antes das simulações, dimensionou-se a turbina eólica recorrendo às equações descritas na secção 2.2 do capítulo 2. Tendo em conta a Figura usou-se T SR = 2. Para o parâmetro α r usou-se o valor de 5 obtido com base na análise dos parâmetros C Lr e o C Dr, com recurso a simulações numéricas efetuadas no programa SolidWorks. A título de exemplo, na Figura representa-se uma simulação numérica, onde v r = 70 m/s, U = 10 m/s e o tipo de perl é o NACA Figura 3.1.1: Distribuição da velocidade do vento em redor do perl NACA 8311, vista 2D.
45 3.1 Simulações numéricas 23 Por m, obteve-se a turbina eólica representada na Figura 3.1.2, semelhante ao modelo comercial Rutland 913 adquirido. Figura 3.1.2: Representação 3D do modelo de WT dimensionada. Os valores da geometria das pás apresentam-se na Tabela 3.1. Tabela 3.1: Valores dos parâmetros resultantes do dimensionamento da turbina eólica utilizada. Secção da pá r (mm) c r (mm) e r (mm) β r ( ) De seguida projetou-se o dispositivo C-D, conforme se ilustrada na Figura
46 Simulações numéricas (a) (b) Figura 3.1.3: Dispositivo C-D: (a) representação esquemática com a respetiva legenda; (b) vista de corte 3D. Para cada ensaio foi necessário denir vários parâmetros no programa SolidWorks, como por exemplo a velocidade do vento, o tamanho da malha, Figura 3.1.4, e a velocidade angular. Figura 3.1.4: Malha em 3D utilizada nas simulações. A malha tinha 0.5 x 0.4 x 0.4 m com mais de 1.8 milhões de células irregulares, onde 89% correspondiam ao uido e o restante às superfícies do sólido, WT e C-D. Inicialmente efetuaram-se simulações só com a turbina, sendo U constante com o valor de 10 m/s e a velocidade angular de 3517 rpm. De seguida efetuaram-se simulações com o sistema turbina e C-D, sendo U = 10 m/s e a velocidade angular de 6000 rpm. Como resultados, obtiveram-se os valores
47 3.2 Ensaios no túnel de vento 25 da potência e da velocidade média do vento, U m, em redor da turbina, num volume cilíndrico denido para o efeito conforme se representa na Figura Figura 3.1.5: Volume do cilindro junto à WT utilizado nas simulações para obter U m. Por m efetuaram-se simulações com a turbina e com o sistema turbina e C-D, sendo U e a velocidade angular variável. O valor de U variou entre 5.68 e 10 m/s e a velocidade angular entre 2013 e 5180 rpm. Como resultados, obtiveram-se os valores da potência e da U m em redor da turbina. 3.2 Ensaios no túnel de vento Os ensaios experimentais efetuados com a turbina e com o sistema turbina e C-D foram realizados num túnel de vento sub-sónico como se mostra na Figura O túnel de vento inclui um variador de frequência que permite realizar ensaios com várias velocidades do vento, U. Figura 3.2.1: Túnel de vento usado nos ensaios. O túnel foi previamente calibrado tendo-se registado o valor da velocidade do vento em função da frequência dada pelo variador de frequência. A
48 Ensaios no túnel de vento velocidade do vento foi medida com recurso a um tubo de Pitot conforme se mostra na Figura Figura 3.2.2: Tubo de Pitot para medição da velocidade do vento. As fórmulas utilizadas para calcular a velocidade do vento com o tubo de Pitot encontram-se no Anexo B. A curva de calibração resultante dos cálculos está representada no Figura Figura 3.2.3: Calibração do túnel de vento: velocidade do ar em função da frequência do ventilador. As pás da turbina com m de diâmetro (D) do rotor, foram construídas em MDF (medium density berboard ) e o spinner em alumínio. Para
49 3.2 Ensaios no túnel de vento 27 a sua construção foi utilizada uma máquina de controlo numérico sendo os acabamentos feitos à mão. Na Figura pode visualizar-se uma sequência de algumas fases da maquinagem das pás. Figura 3.2.4: Sequência de algumas fases da maquinagem e do acabamento das pás. Um suporte metálico foi utilizado para xar a turbina e o gerador no interior do túnel de vento. Na Figura (a) pode visualizar-se o conjunto turbina, gerador e suporte metálico e na Figura (b) a xação no túnel de vento. (a) (b) Figura 3.2.5: Turbina eólica utilizada nas experiências no túnel de vento. Com base nas simulações foram construídos dois C-D a usar nos ensaios experimentais: C-D 1 e C-D 2. Ambos os C-D têm L=D=0.108 m e F=10 mm. Outras características geométricas do C-D 1 são: α i = 8 e α o = 16 e do C-D 2 α i = 14 e α o = 18. Os dois C-D foram construídos em poliestireno extrudido (roofmate). Na sua construção foi utilizado a máquina CNC, sendo os acabamentos feitos à mão. Na Figura podem visualizar-se algumas etapas da sua construção.
50 Ensaios no túnel de vento Figura 3.2.6: Algumas etapas da construção dos C-D. Dois suportes metálicos foram construídos para os C-D de modo a xar os sistemas turbina e C-D no túnel de vento. Na Figura (a) pode visualizar-se o conjunto C-D com o suporte metálico e na Figura (b) pode visualizar-se o sistema turbina e C-D xados no interior do túnel de vento. (a) (b) Figura 3.2.7: C-D 1 utilizado nas experiências no túnel de vento. No túnel de vento foram, então, feitos ensaios experimentais só com a turbina, com o sistema turbina e C-D 1 e com o sistema turbina e C D 2. Os ensaios experimentais iniciavam-se com a escolha da frequência correspondente à velocidade do vento pretendida. Inicialmente efetuaram-se ensaios só com a turbina e por m efetuaramse ensaios com o sistema turbina e C-D. O valor de U variou entre 5.68 e 10 m/s. Para esta gama de velocidades registou-se a velocidade angular entre 1980 e 4601 rpm para os testes com a turbina, entre 2470 e 6334 rpm para o sistema turbina e C-D 1 e entre 2511 e 6022 rpm para o sistema turbina e C-D 2.
51 3.3 Experiências de campo 29 A medição das variáveis a medir foram com um voltímetro e um amperímetro. Foram registados valores da corrente (I) e tensão (V ) que permitiram depois calcular o valor da potência elétrica, P = V I. 3.3 Experiências de campo Para as experiências de campo foi adquirida a turbina do tipo Rutland 913, como referido e conforme se mostra na Figura (a). Nestas experiências foram obtidos resultados preliminares usando-se apenas a turbina e, portanto, sem dispositivo C-D. O equipamento para a realização destas experiências é constituído por quatro componentes: turbina; poste de sustentação; sistema de aquisição de dados e sistema de controlo de potência. A turbina foi adquirida e o poste foi dimensionado e construido para se adaptar à turbina. Também foi dimensionado e construido um sistema de aquisição de dados. Por último, foi dimensionado e construído o sistema de controlo de potência para este caso especíco. Este serve para otimizar a potência que a turbina pode fornecer. A turbina adquirida é composta pelo gerador elétrico, com uma potência máxima de 300 W, e pelo rotor, com 0.91 m de diâmetro e com seis pás. O poste foi previamente desenhado em 3D e é composto por três secções, conforme se mostra na Figura 3.3.1: zona da sapata; união dos tubos e xação da turbina.
52 Experiências de campo (a) (b) (c) Figura 3.3.1: Representação 3D das partes do poste: (a) sapata; (b) união dos tubos; (c) xação da turbina. O poste que sustenta a turbina é composto por dois tubos com 100 mm de diâmetro, 3 mm de espessura e com 6 m de altura cada. O sistema de xação é do tipo espiado e os cálculos associados à estabilidade do poste encontram-se no anexo A. Toda a estrutura foi construída no Laboratório de Tecnologia Mecânica da ESTiG e posteriormente foi colocada num terreno da Quinta do Poulão do IPB. Na Figura podem visualizar-se algumas fases da sua construção e na Figura a sua montagem no terreno. (a) (b) (c) Figura 3.3.2: Fases de construção inicial do poste: (a) furação das várias chapas; (b) fases de limpeza dos tubos e soldagem das chapas aos tubos; (c) pintura dos tubos e acessórios.
53 3.3 Experiências de campo 31 (a) (b) (c) Figura 3.3.3: Fases da instalação da turbina e do poste no terreno da Quinta do Poulão: (a) xação da sapata; (b) xação da turbina e dos sensores; (c) xação de um tubo auxiliar para erguer o poste. Na Figura pode visualizar-se a turbina pronta para as experiências de campo. (a) (b) Figura 3.3.4: Turbina e poste xos no campo: (a) turbina e sensores; (b) sistema global instalado. O sistema aquisição de dados é composto por 3 partes fundamentais: os sensores de corrente, de tensão, de velocidade e de direção do vento; o cartão SD e o microcontrolador, responsável por ler os sensores, processar e guardar os dados no cartão SD.
54 Experiências de campo Os sensores para análise do vento (ver Figura 3.3.5) são da Vector Instruments, sendo o modelo A100LK do sensor de velocidade e o modelo A200P do sensor de direção. Os sensores foram xados no poste a 11 m do solo. (a) (b) Figura 3.3.5: Sensores para análise do vento: (a) sensor de velocidade; (b) sensor de direção. O sensor de corrente usado foi o modelo ACS712 da marca Allegro e o sensor de tensão é dado por um divisor de tensão. A programação do microcontrolador foi feita em linguagem C e C++ usando como interface o programa Arduino. O sistema de aquisição de dados guarda-os a cada 2 segundos e cria um documento do tipo.txt por dia no cartão SD. O sistema de controlo de potência é composto por 3 partes principais: os sensores de corrente e da tensão; o microcontrolador e a carga. Os sensores de tensão e de corrente são os já referidos no sistema de aquisição de dados. O microcontrolador é o responsável por gerir a potência que a turbina eólica pode dar. A programação deste é feita em linguagem C usando como interface o programa MPLAB IDE. O algoritmo utilizado é o da perturbação-observação, P-O (para uma melhor compreensão consultar [61, 62]) associado a um conversor abaixador (em inglês buck converter [14, 63]). Na Figura pode-se observar o diagrama do conversor abaixador utilizado.
55 3.3 Experiências de campo 33 Figura 3.3.6: Representação esquemática do diagrama do circuito de controlo da potência fornecida à carga. Adaptado de [14]. A carga é um sistema que aproveita a energia produzida pela turbina eólica. Pode aproveitar a energia em diversas formas como por exemplo: iluminação, movimento de objetos, aquecimento ou armazenamento em baterias. Na Figura visualiza-se a carga utilizada instalada no poste. No caso em estudo a energia elétrica é dissipada em forma de calor. Figura 3.3.7: Visualização da carga xa ao poste, utilizada para dissipar a energia elétrica. O PCB (em inglês printed circuit board) do sistema de controlo e do sistema de aquisição de dados foram executados no programa Eagle (ver Figura (a)), estando ambos na mesma placa, conforme se visualiza na Figura (b). O processo de construção da placa de circuito impresso, conhecido como PCB é um processo químico, este foi realizado no Laboratório de Química da ESTiG.
56 Experiências de campo (a) (b) Figura 3.3.8: Sistema de aquisição de dados e sistema de controlo de potência: (a) representação esquemática do circuito; (b) placa impressa com os componentes integrados. O sistema de controlo da potência e aquisição de dados é alimentado por uma bateria de 12 volts. Depois da ligação e anação de todos os componentes procedeu-se ao registo dos parâmetros em análise. Num futuro próximo serão efetuadas as experiências com o sistema turbina e C-D. Como registos principais temos a velocidade do vento, U, corrente, I, e tensão, V, que permitiram calcular o valor da potência elétrica, P = V I.
57 Capítulo 4 Resultados e Discussão Apresentam-se as condições quer das simulações numéricas quer dos ensaios experimentais no túnel de vento, quer das experiências de campo e os resultados obtidos. Para melhor compreensão e análise dos resultados obtidos nas simulações numéricas, nos ensaios no túnel de vento e nas experiências de campo, a potência e a velocidade média do vento junto da turbina foram adimensionalizadas. A velocidade adimensional do vento é o quociente entre a velocidade média do vento do sistema turbina e C-D e a velocidade média do vento só da turbina: U mc D /U mw T. Do mesmo modo tem-se para a potência adimensional o quociente: P C D /P W T. 4.1 Simulações numéricas Para todas as simulações obtiveram-se valores da velocidade média do vento (U m ) junto da turbina e da potência (P). Os parâmetros U m e P foram obtidos em duas fases distintas. Na primeira fase o valor de U foi de 10 m/s, com 3517 rpm para as simulações só com turbina e 6000 rpm para a simulações efetuadas com o sistema turbina e C-D. Nesta fase foram realizadas 22 simulações. Na segunda fase o valor de U e da velocidade angular da turbina variaram de acordo com a Tabela 4.1. Nesta fase efetuaram-se 18 simulações. 35
58 Simulações numéricas Tabela 4.1: Condições das simulações só com turbina e para sistema turbina e C-D 1 e C-D 2. Velocidade angular (rpm) U (m/s) WT WT + C-D 1 WT + C-D Inicialmente obtiveram-se os resultados correspondentes só à turbina. Depois, nas simulações efetuadas com o sistema turbina e C-D usaram-se três valores do tamanho da ange do C-D: F=5 mm, F=10 mm e F= 15 mm. Para cada valor de F foram efetuadas 7 simulações, estando as condições geométricas de cada C-D representadas na Tabela 4.2. Tabela 4.2: Código das simulações e condições geométricas dos C-D. Simulações α i ( ) α o ( ) Para além da velocidade média junto da turbina, o software permitia obter também a distribuição de velocidades em redor da turbina e do C-D. A título de exemplo representa-se na Figura a distribuição de velocidades para uma simulação só com turbina e na Figura para uma simulação do sistema turbina e C-D 1, para U = 10 m/s.
59 4.1 Simulações numéricas 37 (a) (b) Figura 4.1.1: Distribuição de velocidades do vento em redor da turbina: (a) vista corte-trás; (b) vista corte lateral. (a) (b) Figura 4.1.2: Distribuição de velocidades do vento em redor do sistema turbina e C-D 1: (a) vista corte-trás; (b) vista corte lateral. Comparando a Figura e a Figura pode-se observar um aumento da velocidade do vento no C-D a partir do centro da turbina até à extremidade das pás. Como se verica um aumento da velocidade do vento na zona exterior-trás do C-D e uma diminuição da velocidade do vento na zona mais estreita do C-D (interior), existe uma variação de pressão entre estes 2 pontos. A variação de pressão leva a que o uxo de ar se mova da zona de maior pressão para a zona de menor pressão, levando a que haja um aumento da velocidade do vento na zona interior do C-D Velocidade média do vento junto da turbina Os resultados da velocidade média adimensional (U mc D /U mw T ) obtidos na primeira fase das simulações, de acordo com a Tabela 4.2, estão repre-
REDES DE SENSORES SEM FIO APLICADOS EM UMA TURBINA EÓLICA VERTICAL DE PEQUENO PORTE
REDES DE SENSORES SEM FIO APLICADOS EM UMA TURBINA EÓLICA VERTICAL DE PEQUENO PORTE Guilherme de Souza Araújo Faculdade de Engenharia Elétrica CEATEC/PUC-Campinas araujogs0@gmail.com Júlio César Penereiro