O presente trabalho utiliza o método Constructal Design para desenvolver o estudo numérico da configuração de materiais de alta condutividade térmica em forma de “Y” que minimiza a resistência ao fluxo de calor, quando áreas ocupadas pelos materiais de alta e baixa condutividades são mantidas constantes. Para a solução numérica da equação diferencial da difusão do calor e suas respectivas condições de contorno, foi utilizado o software MATLAB ®, mais especificamente a ferramenta PDETOOL, Partial Differential Equations Tool. O objetivo deste trabalho é a minimização da resistência térmica do sistema gerador de calor com baixa condutividade térmica com a utilização de vias em formato de Y com material de alta condutividade térmica e volume constante, sendo variáveis os comprimentos e espessuras do material dos ramos simples e bifurcados. Todas as possibilidades geométricas foram avaliadas e a geometria ótima foi aquela que conduziu a menor resistência térmica. Duas condições são apresentadas, a primeira tem os ramos e a base da geometria “Y” com igual condutividade térmica. Os resultados para esta configuração mostram que existem valores específicos para os graus de liberdade que minimizam a resistência térmica. Nesse caso, os ramos se degeneraram e a configuração ótima tem a forma de um “V”. A segunda configuração apresenta combinações de condutividade térmica diferentes, para os ramos e a bases. Para estes casos obteve-se um valor otimizado próximo de 1 para a razão entre os comprimentos dos ramos simples e bifurcados, indicando que a configuração otimizada tem realmente a forma de um “Y” o que demonstra a dependência entre a geometria e as condições impostas pelo meio. Embora o design inicial do Y possa assumir diversas configurações, quando comparado o primeiro design com o design final, no caso do Y com iguais condutividades térmicas se conseguiu uma melhora superior a 28% e no caso do Y com condutividades diferentes mais de 30 %. Finalmente, este trabalho mostra que a geometria otimizada é aquela que melhor distribui as imperfeições, isto é, os pontos quentes (pontos de temperatura máxima), o que está de acordo com o princípio da ótima distribuição das imperfeições. / The present work used the method Constructal Design to develop numerical analyses of pathways of high thermal conductivity in "Y" shape which minimizes the thermal resistance when areas occupied by the materials of high and low conductivities are kept constant. For the numerical solution of the differential equations of heat diffusion and their boundary conditions, we used the MATLAB ® software, specifically the PDETOOL tool. The aim was to minimize the thermal resistance of the heat generator system with low thermal conductivity with the use of Y-shaped pathways with high thermal conductivity and constant volume, with variable lengths and thicknesses of material from stem and forked branches. All geometric possibilities were evaluated and the optimal geometry was that which resulted in lower thermal resistance. Two conditions were studied. In the first one the stem and branches of the "Y" have equal thermal conductivity. The results for this configuration show that there are specific values for the degrees of freedom to minimize the thermal resistance. In this case, the branches have degenerated and the optimum configuration has the shape of a "V". The second configuration offers different combinations of thermal conductivity, for branches and bases. For these cases we obtained a optimized value close to 1 for the ratio between the lengths of stem and bifurcated branches, indicating that the optimized configuration actually has the shape of a "Y" which shows the dependency of geometry and condition imposed by the environment. Although the initial design of Y can take various configurations, when compared the first design to the final design, in the case of Y with equal thermal conductivity it this improvement was achieved an improvement greater than 28% and in the case of Y with different conductivities over 30%. Finally, this study showed that the optimized geometry is the one that better distributes imperfections, this is, hot spots (points of maximum temperature), which is in accordance with the principle of the optimal distribution of imperfections.

Submitted by Raquel Vergara Gondran (raquelvergara38@yahoo.com.br) on 2016-04-11T17:52:57Z No. of bitstreams: 1 tese_daniele.pdf: 1693586 bytes, checksum: a89420a4bb70aea354ae5524f20ac563 (MD5) / Approved for entry into archive by cleuza maria medina dos santos (cleuzamai@yahoo.com.br) on 2016-05-04T20:39:58Z (GMT) No. of bitstreams: 1 tese_daniele.pdf: 1693586 bytes, checksum: a89420a4bb70aea354ae5524f20ac563 (MD5) / Made available in DSpace on 2016-05-04T20:39:58Z (GMT). No. of bitstreams: 1 tese_daniele.pdf: 1693586 bytes, checksum: a89420a4bb70aea354ae5524f20ac563 (MD5) Previous issue date: 2009 / Ao contrário dos bioprocessos submersos, que são amplamente utilizados e estudados, bioprocessos em estado sólido (BES) ainda são carentes de estudos de modelagem e simulação, o que aponta para o grande potencial de otimização. A dificuldade no aperfeiçoamento de BES está associada a problemas com a dissipação do calor gerado pelas atividades metabólicas do microrganismo durante o crescimento. Esta dificuldade na transferência de calor dentro do biorreator pode levar a zonas de altas temperaturas, que afetam adversamente a produtividade. A modelagem matemática é uma ferramenta essencial para otimizar bioprocessos. Através de modelos matemáticos é possível otimizar as variáveis operacionais para controle do bioprocesso e também analisar o design do biorreator. A otimização geométrica, de acordo com a Teoria Constructal, visa melhorar o desempenho do biorreator através, por exemplo, de minimizar a temperatura no interior do leito a níveis ótimos para o cultivo. O presente trabalho apresenta projetos de biorreatores para BES, todos com geometria otimizada, obtidos a partir de experimentação numérica, através de um software de computational fluid dynamics (CFD). O modelo matemático utilizado era preditivo e significativo ao nível de confiança de 95%. A otimização geométrica foi apresentada em função das condições operacionais do cultivo. Para o biorreator de coluna e leito fixo com paredes isoladas, foram apresentadas as geometrias ótimas em função da velocidade, da vazão e da temperatura do ar de admissão. Para uma temperatura do ar de admissão de 29,5 ºC, as configurações ótimas ((D/L)opt) variaram entre 1,0 e 2,4 para uma faixa de velocidade de admissão do ar entre 0,003 e 0,006 m s-1 . Relacionando com vazão, as razões mostraram-se ótimas entre 2,2 ≤ (D/L) ≤ 2,6 quando operando sob 3,3 a 3,5 10-5 m3 s-1 . Outro biorreator estudado foi o biorreator modular, composto de módulos elementares com geometria otimizada, sendo adaptável a diferentes escalas de produção e de fácil montagem. As configurações ótimas dos módulos de geometria retangular e seção quadrada foram apresentadas para diferentes volumes de módulos, em função da temperatura e da velocidade do ar de admissão. Foi observado que o volume máximo do módulo sem resfriamento externo é 5 L, para uma velocidade do ar de admissão acima de 0,0045 m s-1 e temperatura inferior ou igual 29,0 ºC O último biorreator proposto foi o biorreator hollow, semelhante a um biorreator de coluna e leito fixo, porém com um duto oco inserido nele. O duto interno tem inúmeros furos perpendiculares às suas paredes, mas sua saída é isolada, permitindo que o ar penetre no meio poroso. A geometria otimizada do biorreator hollow foi apresentada em função da fração de volume do duto interno, da razão entre os diâmetros de entrada e saída do duto interno, da vazão e da temperatura do ar de admissão. Em comparação com o biorreator de coluna convencional de mesmas dimensões e sob mesmas condições operacionais, o biorreator hollow apresentou temperatura máxima mais baixa, demonstrando que o projeto é eficiente para resfriar o meio poroso. Concluiu-se, enfim, no presente trabalho, que a geometria é um parâmetro importante e a sua otimização pode beneficiar o desempenho do biorreator. / Unlike the submerged bioprocesses, that was wildly used and studied, solid state bioprocess (SSB) are still poorly studied with respect to modeling and simulation, what indicates a big potential for optimization. The difficulty in the BES improvement is associated to problems with the dissipation of heat generated by metabolic activities of microorganisms during growth. This difficulty in transferring heat from the bioreactor could lead to areas with high temperature, which usually affect the productivity adversely. The mathematical modeling is an essential tool for optimizing bioprocesses. Using mathematical models it is possible to optimize operational variables to control the bioprocess and also explore the design of the bioreactor. The geometric optimization, according Constructal Theory, aims to improve the performance of the bioreactor through, for example, minimizing the temperature inside the bed to optimum levels for the bioprocess. The present work presents designs of bioreactors to SSB, all with optimized geometry, obtained from numerical experiments, by computational fluid dynamics (CFD) software. The mathematical model used has been predictive and significant at .95 level of confidence. The geometric optimization was presented as function of operational conditions of the cultivation. For the column fixed bed bioreactor with isolated wall, the optimal configurations are shown as function of flow, velocity and temperature of inlet air. For a inlet air temperature of 29.5 ºC, the optimal configurations ((D/L)opt) varied between 1.0 e 2.4 to a range of inlet velocity between 0.003 e 0.006 m s-1 . Relating with the volumetric flow, the optimal ratios presented between 2.2 ≤ (D/L) ≤ 2.6 when operating under 3.3 a 3.5 10-5 m3 s-1 . Other studied bioreactor was the modular bioreactor, consisting of elementary modules with optimized geometry, being adaptable to different scales of production and easy assembly. The optimal configurations of the modules with rectangular geometry and square section were shown depending on the volume of modules and the temperature and velocity of inlet air. It was observed that the maximum volume of the module without external cooling was 5 L, for a inlet velocity upper 0.0045 m s-1 and temperature smaller or equal to 29.0 ºC. The last proposed bioreactor was the hollow bioreactor, similar to a column fixed bed bioreactor, but with an empty duct inserted on it. The internal duct has innumerable holes perpendicular to its wall (the inlet port), but its end is insulated, allowing the air penetrates into the porous medium. The optimized geometry of hollow bioreactor was presented in function of the volume fraction of internal duct, the ratio between the diameters of inlet and outlet of the internal duct, the flow rate and temperature of the inlet air. Comparing with the conventional column bioreactor with the same configuration and same operational conditions, the hollow bioreactor showed a lower maximum temperature. This demonstrates that the project is efficient at cooling the porous medium. Finally, it was concluded that the geometry is an important parameter and its optimization can benefit the performance of the bioreactor.

Trocadores de Calor Solo-Ar (TCSA) são dispositivos usados para melhorar o condicionamento térmico de ambientes construídos. Eles funcionam através da ventilação do ar por um ou mais dutos enterrados, utilizando o solo como fonte ou sumidouro de calor. Em virtude da defasagem entre as temperaturas do ar e das camadas superficiais da Terra, é possível resfriar o ar no verão e aquecê-lo no inverno. Seus princípios de operação baseiam-se na mecânica dos fluidos e transferência de calor, áreas onde a teoria Construtal tem sido usada para melhorar o desempenho, ou reduzir imperfeições, de vários sistemas térmicos, após uma estratégica análise de suas geometrias (ou de seu design), também chamada de método Design Construtal que se fundamenta na lei construtal. Dessa lei, projetos em engenharia devem começar, por exemplo, descobrindo as arquiteturas que facilitam o escoamento entre uma fonte pontual e um volume, ou vice-versa. Explorando esses conceitos para TCSA, onde a corrente é o calor que flui entre os dutos e o solo, esta tese objetiva centralmente avaliar o desempenho térmico desses dispositivos em função de possíveis desenhos de suas estruturas. Particularmente, isso foi feito: (1) usando um duto e variando seu diâmetro e vazão de ar; (2) inserindo novos dutos, até um total de cinco, mantendo a vazão de ar constante, e estudando diversas disposições geométricas. Numericamente, foram empregados dois modelos testados e validados a partir de dados experimentais, que foram simulados no código comercial de dinâmica dos fluidos computacional, FLUENT, o qual baseia-se em métodos de volumes finitos. As geometrias e malhas foram geradas no software GAMBIT. Dos resultados numéricos, esta tese também avançou desenvolvendo novos modelos para o conceito de potencial térmico instantâneo, que é uma medida das diferenças de temperatura entre a saída e a entrada dos dutos de TCSA. Descobriu-se que estes modelos podiam ser simplificados em termos de funções senoidais, facilitando análises e definições analíticas subsequentes, criando um quadro resumo para grandes volumes de dados simulados. Das avaliações dos resultados nos estudos com um duto, foram encontradas faixas para a obtenção de rendimentos térmicos anuais estimados em pelo menos 80%, com mudanças no diâmetro e/ou vazão do ar. Analisando múltiplos dutos, foram reveladas melhorias significativas de rendimento com: (1) o aumento no número de dutos; (2) a redução, até certos limites, da razão entre os espaçamentos verticais e horizontais entre eles; (3) a ampliação da razão entre o volume da instalação e o volume do domínio computacional. Destaca-se, ainda, que em todos os casos foram calculados balanços energéticos, com a estimativa das taxas (e quantidade) de calor trocado, bem como da energia elétrica consumida pelos ventiladores, mostrando sua viabilidade no condicionamento de ar com redução do consumo de energia elétrica. / The Earth-Air Heat Exchangers (EAHE) are devices used to improve the thermal conditions of built environments. They work by blowing the air inside buried ducts to use the soil as heat source or sink. Due to a phase difference between the air and ground temperatures, it is possible to cool the air in the summer and heat it in the winter. These operating principles are based on the areas of fluid mechanics and heat transfer, where the Constructal theory has been used to improve the performance, or reduce imperfections, of several thermal systems, after a strategic evaluation of their geometries (or their design), also called Constructal Design method, which is anchored by the constructal law. From this law, engineering projects should begin, for example, finding the architectures which ease the flow between a source point and a volume, or vice versa. Bringing these concepts to EAHE, where the current is the heat which flows between the ducts and the soil, this thesis mainly aims to assess the thermal performance of EAHE due to possible designs of their layouts. Particularly, this was done by: (1) adopting one duct and varying its diameter and/or the air flow; (2) using up to five ducts, keeping the air flow constant, and studying various geometric shapes. Numerically, two models were used, which have been tested and validated from experimental data, and the simulations were done in the CFD software FLUENT, which is based in the finite volume methods. The geometries and meshes were generated by the software GAMBIT. From the numerical results, this thesis also pursued developing new models for the concept of the instantaneous thermal potential, which is a measure of the temperatures differences between the ducts outlet and inlet. It was found that these models can be simplified in terms of sinusoidal functions, helping analysis, definitions, and creating a framework to summarize large volumes of simulated data. From the evaluation of the data in the studies with one duct, ranges for changes in their diameter and/or air flow were found to obtain at least 80% of estimated annual efficiency. Analyzing multiple ducts, significant improvements in efficiency were also obtained by: (1) increasing the number of ducts; (2) reducing, to some extent, the ratio between their vertical and horizontal spacings; (3) increasing the ratio between the installation volume and the computational domain. It should also be highlighted that in all cases the energy balances were computed, to find the rates (and amounts) of exchanged heat, as well as the electric energy consumed by the fans, showing the economic viability of using EAHE as devices for air conditioning.

A conversão da energia das ondas dos oceanos em energia elétrica é uma alternativa para o problema da falta de combustíveis fósseis. Uma das possibilidades de aproveitamento é através de dispositivos cujo princípio de funcionamento é o de Coluna de Água Oscilante (CAO). No presente trabalho o objetivo é, através da modelagem computacional e do emprego do Design Construtal, maximizar a potência hidropneumática de um dispositivo conversor de energia das ondas do mar do tipo CAO. São analisados diferentes eixos da restrição física, no formato elíptico, que representa a turbina, e duas formas geométricas na região de transição entre a câmara hidropneumática e a chaminé do dispositivo CAO: trapezoidal e semicircular. Considerando um domínio bidimensional, as restrições para estes problemas são: Área da restrição elíptica (AR), Área total do dispositivo (AT) e razão entre a área da restrição elíptica e a área total (ϕn). Os graus de liberdade analisados são: a razão entre os comprimentos dos eixos da restrição elíptica (d1/d2) para o caso da restrição física da turbina, o ângulo de inclinação da parede (α) para o caso com região de transição trapezoidal, o raio (r) e H2/l (razão entre altura e comprimento da chaminé de saída da câmara CAO) para o caso com região de transição semicircular. Para a solução numérica é empregado um código de dinâmica dos fluidos computacional, FLUENT®, baseado no Método de Volumes Finitos (MVF). O modelo multifásico Volume of Fluid (VOF) é aplicado no tratamento da interação água-ar. O domínio computacional é representado por um tanque de ondas com um dispositivo CAO acoplado. Os resultados obtidos indicam que, para o estudo da região de transição trapezoidal o desempenho do conversor tem aproximadamente o mesmo desempenho para todas as geometrias estudadas. A região de transição semicircular, apresenta resultados para os quais foi possível otimizar a potência hidropneumática. O estudo da turbina indica que foi possível determinar uma geometria capaz de converter a energia da onda incidente ao dispositivo, sem que ocorresse a obstrução do escoamento de ar na chaminé do dispositivo CAO. Assim, mostra-se a relação entre o método Design Construtal e o clima de ondas na definição das dimensões que maximizam a potência hidropneumática. / The conversion of ocean’s wave energy into electrical energy is an alternative for the scarcity of fossil fuels. One of the possibilities of energy use is through devices, whose operating principle is the Oscillating Water Column (OWC). In this work the aim is, through computer modeling and the Constructal Design, to maximize hydropneumatic power of a power converter device type OWC. Different axes of physical constraint with elliptical shape, representing the effect of the turbine , are analyzed. Two geometric shapes in the transition region between the hydropneumatic chamber and the chimney OWC device, trapezoidal and semicircular, are also analyzed. Considering a two-dimensional domain the restrictions for this problem are: Elliptical restriction area (AR), Total area device (AT) and the ratio between the area of the elliptical restraint and the total area (ϕn). The considered degrees of freedom are: the ratio between the lengths of the axes (d1/d2) of the elliptical restraint, for the turbine’s physical constraint case, the inclination angle (α) of the wall for the trapezoidal transition case, and the radius (r) and H2/l (ratio between height and length of output chimney CAO) for the semicircular transition region case. For the numerical solution, a commercial code of computational fluid dynamics, FLUENT®, which is based on the Finite Volume Method (FVM), is employed. The multiphase model Volume of Fluid (VOF) is applied in the treatment of water-air interaction. The computational domain is represented by a wave tank with a fixed OWC device. The obtained results indicate that, for the study of the trapezoidal transition region, the performance of converter don’t seems to be compensatory only by changing the geometry of the trapezoidal area. However, for the semicircular transition region, it was possible to optimize a hydropneumatic power. The study of turbine effect indicates a geometry capable of converting the energy of the incident wave to the device, without obstructing the air flow in the chimney of de OWC, showing the relationship between the Constructal Design method and the wave climate in the definition of the dimensions that maximize the hydropneumatic power.

Elementos estruturais como as placas finas fazem parte de um grande número de aplicações nas mais diversas áreas da engenharia e são de grande importância para a engenharia naval e aeronáutica, na construção de cascos de embarcações e estruturas offshore, e na construção de fuselagens de aviões. Por constituírem-se em um elemento estrutural esbelto, estão sujeitas a um comportamento mecânico diferenciado denominado de flambagem, proveniente de um carregamento de compressão uniaxial. O fenômeno da flambagem pode ser dividido em flambagem elástica e elasto-plástica, sendo dependente de aspectos dimensionais, construtivos e/ou operacionais. A inclusão de perfurações em placas provoca uma redistribuição de suas tensões internas, afetando não apenas a sua resistência, mas também as suas características de flambagem. Neste trabalho, desenvolveu-se a análise do comportamento mecânico de placas finas perfuradas de aço, simplesmente apoiadas em suas bordas, e submetidas à compressão. Serão analisados dois graus de liberdade: H/L e H0/L0. Para H/L, serão analisadas placas com H/L = 1,00 e H/L = 0,50, sendo que H e L representam, respectivamente, a largura e o comprimento da placa. Para H0/L0, serão analisadas infinitas possibilidades, sendo que H0 e L0 representam, respectivamente, a largura e o comprimento da perfuração. As placas utilizadas possuem espessura (h) de 10,00 mm e perfuração centralizada. Quanto às perfurações, estas serão dos tipos: oblonga longitudinal, oblonga transversal, elíptica, retangular, losangular, hexagonal longitudinal e hexagonal transversal. Ainda em relação às perfurações, serão consideradas as seguintes frações ϕ = 0,08; 0,10; 0,15; 0,20 e 0,25, sendo que (ϕ) corresponde ao volume da perfuração. Para a determinação das cargas crítica e última de flambagem, foi utilizada a simulação numérica com o auxílio do software Ansys®, que é baseado no método dos elementos finitos. A aplicação do método Design Construtal, possibilitou a determinação das geometrias ótimas para todos os tipos de perfurações, todos os valores de (ϕ) e para todas as relações de H/L. Os resultados obtidos mostram que há influência do tipo, da forma e do tamanho da perfuração na definição das curvas limites à flambagem e das curvas à flambagem elasto-plástica. Foi possível definir, para cada tipo de perfuração e para todos os valores de (ϕ), os pontos de transição entre a flambagem elástica e à elasto-plástica, assim como os pontos que definem os valores máximos para o fator TLNMáx (tensão limite normalizadora). / Structural elements such as thin plates are part of a large number of applications in various areas of engineering and are of great importance for marine and aerospace engineering, construction and offshore structures hulls, and the construction of airplane fuselages. Being a slender structural element, they are subject to a different mechanical behavior known as buckling, caused by a compressive loading. The phenomenon of buckling can be divided in elastic and elasto-plastic buckling, being dependent dimensional, construction and / or operational aspects. The inclusion of perforations in plates causes a redistribution of its internal stress, affecting not only their resistance but also their buckling characteristics. In this work it was performed the analysis of the mechanical behavior of thin perforated steel plates, simply supported on its edges, and subjected to compression. In the analysis it was considered two degrees of freedom: H/L and H0/L0. For H/L will be analyzed plates with H/L = 1.00 and H/L = 0.50, wherein H and L represent respectively the width and length of the plate. There are endless possibilities for the relation H0/L0. The studied plates have a thickness (h) of 10.00 mm and centralized perforation. The following types of perforation will be used: longitudinal oblong, transverse oblong, elliptical, rectangular, diamond, longitudinal hexagonal and transverse hexagonal. Also in relation to perforations, it will be considered the following fractions (ϕ = 0.08; 0.10; 0.15; 0.20 and 0.25), wherein (ϕ) corresponds to the volume ratio of the perforation. For determining the critical and ultimate buckling loads it was utilized numerical simulation with the assistance of Ansys® software, which is based on the finite element method. The application of the Constructal Design method of this study made it possible to determine the optimal geometries for all types of perforations, for all values of (ϕ) and all the relations H/L. The results show that there is an influence of the perforation type, shape and size, in defining the limit curves of the buckling and the curves of the elasto-plastic buckling. It was also possible to define, for each type of perforation and for all (ϕ) values, the transition points between elastic and elasto-plastic buckling; as well as the points that define the maximum values for the TLNMáx factor (normalized limit stress).

O presente trabalho destina-se ao estudo numérico da geometria de tubos de seção elíptica que facilite a transferência de calor adimensional e diminua a queda de pressão adimensional (Δ̃) sofrida pelo escoamento. O método aplicado é o Design Construtal, que visa determinar a geometria que apresentará a menor resistência ao escoamento, ou seja, busca-se determinar a razão de aspecto da elipse (=⁄) que favorece a transferência de calor e diminui a queda de pressão do escoamento. O fluido empregado neste estudo apresenta características de viscoplasticidade. A relação entre a tensão cisalhante e a taxa de deformação obedece ao modelo de Herschel-Bulkley modificado. Considera-se que o escoamento é incompressível, laminar, bidimensional, externo e ocorre em regime permanente. A solução numérica do problema proposto é realizada com um código comercial baseado no método dos volumes finitos. É investigada a influência do índice de potência, , sobre a seção elíptica que facilita o escoamento e, para isso, este índice é variado de 0,4 a 1. A influência dos números de Reynolds (√), Herschel-Bulkley (√) e Prandtl (√) sobre o comportamento do escoamento também é avaliada. √ é variado de 1 a 40, √ é variado de 1 a 100 e √ é variado de 0,1 a 100 Os resultados mostram que, para um escoamento com √=1, √=1 e √=1, o aumento do índice de potência influencia negativamente na transferência de calor adimensional e a seção elíptica, que maximiza esta transferência de calor adimensional, tende a ser mais alongada na direção do escoamento. Já e influenciam positivamente na transferência de calor adimensional. Para um escoamento com √=1, √=1, =0,4 conclui-se que com o aumento de a razão de aspecto ótima (q,opt), do ponto de vista térmico, diminui. Quando é considerado um escoamento com √=1, √=1, =0,4 conclui-se que q,opt diminui com o aumento de , ou seja, a elipse torna-se mais alongada no sentido do escoamento. A variação de √ em um escoamento com √=1, √=1, =0,4 mostra que o aumento deste parâmetro acarreta em aumento da taxa de transferência adimensional e de Δ̃. / The present work is aimed at the numerical study of the geometry of elliptic section tubes that facilitates the dimensionless heat transfer and decreases the dimensionless pressure drop (Δ̃) suffered by the flow. The applied method is the Construtal Design, which aims to determine the geometry that will present the least resistance to the flow, that is, to determine the aspect ratio of the ellipse (=⁄) that favors heat transfer and decreases the flow pressure drop. The fluid used in this study has viscoplasticity characteristics. The relationship between shear stress and strain rate follows the modified Herschel-Bulkley model. It is considered that the flow is incompressible, laminar, two-dimensional, external and occurs in steady state. The numerical solution of the proposed problem is carried out with a commercial code based on the finite volume method. The influence of the power index, n, on the elliptical section facilitating the flow is investigated, and for this, the index is varied from 0.4 to 1. The influence of the Reynolds number (√), Herschel-Bulkley number (√) and Prandtl number (√) on the flow behavior is also evaluated √ is varied from 1 to 40, √ is varied from 1 to 100 and √ is varied from 0.1 to 100. The results show that for a flow with √=1, √=1 and √=1, the increase of the power index negatively influences the dimensionless heat transfer and the elliptic section, which maximizes this dimensionless heat transfer, tends to be more elongated in the direction of flow. Already √ and √ influence positively the dimensionless heat transfer. For a flow with √=1, √=1, =0.4 it is concluded that with the increase of √ the optimum aspect ratio (q,opt), from the thermal point of view, decreases. When a flow is considered with √=1, √=1, =0.4 it is concluded that q,opt decreases with the increase of √, that is, ellipse becomes more elongated in the flow direction. The variation of √ in a flow with √=1, √=1, =0.4 shows that the increase of this parameter causes an increase of the dimensionless transfer rate and Δ̃.

O princípio operacional do dispositivo de galgamento consiste de uma estrutura que utiliza uma rampa para direcionar as ondas incidentes para o reservatório. A água armazenada retorna para o oceano após a passagem por uma turbina que está acoplada a um gerador de energia elétrica. O presente trabalho propõe dois estudos numéricos a respeito de um conversor de energia das ondas do mar do tipo galgamento. Para ambos os casos, o objetivo do estudo é a aplicação do método Design Construtal na definição da melhor forma para a rampa de modo a maximizar a massa de água que entra no reservatório, conduzindo a uma maior geração de energia elétrica. O grau de liberdade b/B, isto é, a razão entre a base superior e a base inferior da rampa trapezoidal, foi otimizado, mantendo-se fixos a área total do tanque de ondas, a área da rampa e as características da onda. Para a análise numérica do princípio de funcionamento deste dispositivo foi empregado um domínio computacional tridimensional (3D), gerado através do software GAMBIT, onde o conversor é acoplado a um tanque de ondas regulares. A solução das equações de conservação e a equação do transporte da fração volumétrica foi realizada com o código comercial de Dinâmica dos Fluidos Computacional FLUENT, que é baseado no Método de Volumes Finitos (MVF). Aplica-se o modelo multifásico Volume of Fluid (VOF) no tratamento da interação água-ar. Para o primeiro estudo, as características da onda regular empregada estavam em escala de laboratório. Os resultados mostraram que houve uma razão ótima (b/B)o = 0.43, que maximiza a quantidade de água que entra no reservatório para o caso estudado. Para ambos os casos, a razão ótima foi encontrada para o extremo inferior do grau de liberdade, além dos resultados apontarem um aumento significativo na massa admitida no reservatório e, por consequente, um maior aproveitamento das ondas incidentes. / The operational principle of an overtopping device consists of a structure which utilizes a ramp to direct incident waves to the reservoir. The stored water returns to the ocean after passing through a turbine that is coupled to an electric generator. This work proposes two numerical studies of a WEC of sea waves of the type overtopping. In both cases, the objective of the study is the application of Constructal Design method to define the best geometry of the ramp which maximizes the mass of water entering the reservoir, leading to increase the generation of electricity. The degree of freedom b/B, that is, the ratio between the upper base and the lower base of the trapezoidal ramp, has been optimized, keeping fixed the total area of the wave tank, the area of the ramp and the wave characteristics. For the numerical analysis of the working principle of this device it was used a three-dimensional computational domain (3D) generated by GAMBIT software where the device is inserted to a tank of regular waves. The solution of conservation equations and equation of transport of the volumetric fraction was carried out with the Commercial Code of Computational Fluid Dynamics FLUENT, which is based on Finite Volume Method (FVM). It was applied the multiphase model Volume of Fluid (VOF) in the treatment of the interaction water-air. For the first study, the characteristics of the employed regular wave were on a laboratory scale. The results showed that there were an optimal ratio (b/B)o = 0.43, which maximizes the amount of water entering the reservoir for the case study. For the second study, the characteristics of the regular wave were employed at actual scale and the results showed that there was an optimum ratio (b/B)o = 0.38, which maximizes the amount of water entering the reservoir for the case study. In both cases, the optimum ratio is found for the extreme lower of freedom of degree and the results showed a significant increase in the mass allowed in the reservoir and, consequently, larger use of the incident waves.

O presente trabalho apresenta um estudo numérico bidimensional sobre a otimização da geometria de um dispositivo conversor de energia das ondas do mar em energia elétrica. O objetivo principal é, através da modelagem computacional de um dispositivo cujo principio de funcionamento é o de Coluna de Água Oscilante (CAO) e do emprego de Constructal Design, maximizar a conversão da energia das ondas do mar em energia elétrica. Essa técnica é baseada na Teoria Constructal. O aspecto inédito deste trabalho, em relação aos estudos disponíveis na literatura, é o fato de levar em conta o clima de ondas de uma dada região e, a partir disso, dimensionar o dispositivo de modo que ele tenha um desempenho otimizado. Para tanto, foi empregado o método Constructal Design, os graus de liberdade empregados são: H1/L (razão entre a altura e o comprimento da câmara CAO) e H3 (profundidade de submersão do dispositivo CAO). A relação H2/l (razão entre altura e comprimento da chaminé de saída da câmara CAO) é considerada um parâmetro fixo. Foram realizados estudos levando em conta uma onda em escala de laboratório e um espectro de ondas real. Foi também realizado um estudo sobre a influência da perda de carga da turbina através de uma restrição física. Para a solução numérica foi empregado um código comercial de dinâmica dos fluidos computacional, FLUENT®, baseado no Método de Volumes Finitos (MVF). A geometria e a geração a malha foi realizada no software GAMBIT®. O modelo multifásico Volume of Fluid (VOF) é aplicado no tratamento da interação água-ar. O domínio computacional é representado por um tanque de ondas com o dispositivo CAO acoplado. Os resultados obtidos mostram que é possível estabelecer uma razão de H1/L ótimo, conhecendo-se o clima de ondas, ou seja, o recomendável é que esta razão seja igual a quatro vezes a altura da onda dividido pelo comprimento da onda incidente. / The present work presents a two-dimensional numerical study about the geometric optimization of an ocean Wave Energy Converter (WEC) into electrical energy. The main goal is, through computational modeling of a device whose operating principle is the Oscillating Water Column (OWC) and from employment Constructal Design, to maximize the conversion of energy of ocean waves into electricity. This technique is based on Constructal Theory. The inedited aspect of this work comparing to the available studies is that it takes into account the wave climate of a given region to design the device so that it achieves optimum performance. Constructal Design is employed varying the degrees of freedom H1/L (ratio between the height and length of OWC chamber) and H3 (lip submergence). While the relation H2/l (ratio between height and length of chimney) is kept fixed. Studies were performed considering a wave on a laboratory scale and a spectrum of real waves. Yet a study of the influence of the turbine pressure losses was performed using a physical constraint. For the numerical solution it is used the Computational Fluid Dynamic (CFD) commercial code FLUENT®, based on the Finite Volume Method (FVM). The geometry and mesh generation was performed in GAMBIT ® software. The multiphasic Volume of Fluid (VOF) model is applied to tackle with the water-air interaction. The computational domain is represented by an OWC device coupled with the wave tank. The results show that it is possible to establish a relationship of H1 / L optimum, if the wave climate is know. It is recommended that this ratio be equal to four times the height of the wave divided by the length of the incident wave.

O presente trabalho tem como objetivo maximizar a potência hidropneumática convertida num dispositivo do tipo Coluna d'Água Oscilante (CAO). Para fazê-lo, o método Constructal Design é aplicado para aprimorar a geometria e a profundidade de submersão do dispositivo. No desenvolvimento do método Constructal são propostos e analisados três graus de liberdade: H1/L (razão entre a altura e comprimento da câmara do dispositivo CAO), H2/l (razão entre a altura da câmara e o comprimento da chaminé) e H3 (profundidade de submersão do dispositivo CAO). As restrições do problema (parâmetros constantes) são a área da câmara A1 e a área total do dispositivo CAO A2. O domínio computacional consiste de um dispositivo CAO inserido num tanque que é submetido a ondas na escala real. A malha é desenvolvida no software Ansys Icem®. O código de Dinâmica dos Fluidos Computacional Ansys Fluent® é empregado para encontrar a solução numérica a qual é baseada no método dos Volumes Finitos. O modelo multifásico Volume of Fluid (VOF) é usado na interação das fases água-ar. Os resultados indicam que a potência hidropneumática máxima obtida é de 190 W para razões de H1/L, H2/l e H3 iguais a 0,135, 6,0 e 9,5 m respectivamente. Por outro lado, o menor valor obtido da potência hidropneumática é de quase 11 W, o que mostra a utilidade do método Constructal, para fornecer uma relação entre o clima de ondas de um lugar determinado e as dimensões ótimas do dispositivo CAO. / The present work aims to maximize the hydropneumatic power converted in an Oscillating Water Column (OWC) device. To do this, Constructal Design is applied to optimize its geometry and submergence. For the development of Constructal method, it has been proposed and analyzed three degrees of freedom: H1/ L (ratio between the height and length of OWC chamber), H2/l (ratio between height and length of chimney), and H3 (submergence). The problem constraints (fixed parameters) are total area of the OWC chamber A1 and total area of OWC device A2. The computational domain consists of an OWC inserted in a tank where waves in a real scale are generated. The mesh is developed in ANSYS ICEM®. The Computational Fluid Dynamics code FLUENT® is used to find the numerical solution which is based on Finite Volume Method (FVM). The multiphasic Volume of Fluid (VOF) model is applied to tackle with the water-air interaction. The results show that the maximum hydropneumatic power obtained was 190 W for H1/L, H2/l e H3 ratios equal to 0.135, 6.0 and 9.5 m respectively. In contrast, the smaller value obtained for the hydropneumatic power is almost 11 W. So, it shows the utility of Constructal Method which provides a relationship between the wave climate of a particular place and the optimal dimensions for the OWC device.

Este trabalho trata da aplicação do método Design Construtal para investigar a transferência de calor através de caminhos de alta condutividade térmica com geometrias definidas. O objetivo é obter a configuração que reduz a temperatura máxima em excesso do sistema considerando que as áreas ocupadas pelos materiais de alta e baixa condutividade são tratadas como constantes. Assim, o objeto de estudo é um volume de área finita onde ocorre a geração de calor. O escoamento da energia térmica para fora do volume é feito através de um caminho condutor de alta condutividade térmica. O trabalho considerou a resistência térmica de contato entre o elemento condutivo e o corpo gerador de calor, onde um terceiro material com resistência térmica equivalente à resistência de contato é interposto entre os dois primeiros. Na solução da equação da difusão do calor, foi realizado um tratamento numérico através de um código baseado em elementos finitos e utilizando o toolbox PDETool, Partial Differential Equations Tool, que pertence ao aplicativo comercial MatLab®. O tratamento numérico foi realizado considerando-se caminhos condutivos com geometrias em forma de "I" e em forma de "T", mantendo-se as frações de área constantes e variando-se os comprimentos dos materiais de alta condutividade e os da resistência térmica de contato. A otimização geométrica foi feita considerando-se os graus de liberdade existentes para cada geometria, onde os valores otimizados para a situação ideal, ou de acoplamento térmico perfeito, foram comparados para os resultados envolvendo a resistência térmica de contato (RTC). Os resultados indicam que a RTC pode aumentar a temperatura máxima em excesso, assim como tem efeito significativo sobre as ótimas configurações calculadas quando a resistência de contato é levada em consideração para ambas as configurações "I" e "T" estudadas. / This work applies Constructal Design to investigate the heat transfer through high conductive pathways with defined geometries. The objective is to find the configuration which reduces the maximal excess of temperature considering the areas with high and low thermal conductivity are constants. Thus, the object studied here is a volume with a finite area and heat generation. The outside heat flux is conducted through a high thermal conductive pathway. Here, special attention is given to the thermal contact resistance between the high conductive pathway and the solid body, where a third material with a thermal resistance equivalent to the thermal contact resistance is inserted between them. A numerical treatment was given in order to solve the heat diffusive equation. It was used a numerical code based on finite elements and the toolbox – PDETool, Partial Differential Equations Tool, which is part of the MatLab® applicative. The numerical treatment was achieved considering "I" and "T" geometries for the high conductive pathways keeping the areas fraction constants and varying the lengths of both high conductive and the equivalent thermal contact layer materials. The optimization was performed considering the degrees of freedom of each geometry, where the optimized values for the ideal situation, i.e., perfect thermal contact were compared with the results considering the thermal contact resistance. The results indicate that the thermal contact resistance can increase the excess of temperature, as well as it has a significant effect on the optimal configurations when using perfect thermal contact or taking into account the thermal contact resistance for "I" and "T" shaped geometries.