Turbinas eólicas usualmente são agrupadas em grandes parques, reduzindo o custo de instalação, transmissão da energia e manutenção periódica. A superposição das esteiras sobre turbinas adjacentes normalmente reduz consideravelmente a capacidade total, objeto de estudo de Micrositing. Porém, por vezes o “efeito Venturi” ocasionado pelas turbinas à montante induz maior velocidade às turbinas adjacentes aumentando o potencial eólico disponível nas linhas consecutivas. De forma inovadora empregar o Design Constructal de Bejan, o modelo do disco atuador genérico e a Dinâmica dos Fluidos Computacional (CFD) para obter a melhor disposição geométrica das turbinas em uma área plana e não rugosa, com foco à maior potência extraída por área de turbinas instaladas. Para tal, modelar e predizer o comportamento da esteira é fundamental, assim como conhecer os modelos de esteira e a aplicabilidade dos métodos empregados. O Design Constructal é a fonte dos parâmetros geométricos base das simulações: o espaçamento entre as turbinas e as razões de diâmetros. Após 64 simulações semi-iterativas e mais de 60 iterativas verifica-se que o maior ganho em potência disponível por área é de 7,37% para a configuração V = 7m/s, S = 3D, d/D = 0.5, L = 3D e 8,48% para a configuração V = 11m/s, S = 3D; d/D = 0.25 & 0.5, L= 0.75D, valor relativo à execução de somente um diâmetro de 100 metros. / Usually wind turbines are grouped in large parks, reducing the cost of installation, energy transmission and periodic maintenance. But the overlapping of the aerodynamical wakes on adjacent turbines reduces the total capacity, Micrositing study. However, the "Venturi effect" caused by the turbines upstream sometimes increases the speed to the adjacent turbines increasing the wind potential available in straight lines. Innovatively employing the Design Constructal Bejan, the model of the actuator disc and Computational Fluid Dynamics (CFD) to search the best geometrical layout of the turbines on a roughless and flat area, focus on higher power extracted by area. To do this, model and predict the wake of behavior is fundamental, as well as know the aerodynamical wakes models and the applicability of the methods employed. The Design Constructal is the source of the simulation’s parameters: spacing between the turbines and the diameter’s ratio. After concluded 64 semi-iterative and iterative simulations, and more than 60 verifies, the best gain in available power per area is 7.37% for the configuration V = 7 m/s; S = 3d; d/D = 0.5; L = 3D. And the gain of 8.48% for the configuration V = 11m/s, s = 3D; d/D = 0.25 & 0.50; L = 0.75D, comparing to the implementation of just 100 meters diameter.

A maioria dos líquidos encontrados na natureza são não newtonianos e o estudo do seu comportamento reológico tem uma importância significante em diferentes áreas da engenharia. Entre eles, existe uma classe de fluidos que exibem pequena deformação aparente quando sujeitos a um nível de tensões inferior a uma tensão limite de escoamento, referido como comportamento viscoplástico. Nesta classe de materiais, alguns apresentam também comportamento elástico quando submetidos a baixas taxas de cisalhamento. A presente Tese tem como objetivo o estudo numérico de escoamentos bidimensionais em regime permanente de fluidos elasto-viscoplásticos através de uma expansão-contração planar. O modelo mecânico é definido pelas equações de conservação de massa e de balanço de momentum acopladas ao modelo elasto-viscoplástico proposto nesta Tese. Esta modelagem é aproximada por um método de elementos finitos multi-campos estabilizado baseado na metodologia de Galerkin mínimos-quadrados que possui como variáveis primais os campos de tensão extra polimérica, velocidade e pressão. As condições de compatibilidade entre os sub-espaços de elementos finitos para tensão extra-velocidade e velocidade-pressão são violadas, permitindo assim a utilização de interpolações de igual ordem. O método estabilizado foi implementado no código de elementos finitos para fluidos não newtonianos em desenvolvimento no Laboratório de Mecânica dos Fluidos Aplicada e Computacional (LAMAC) da Universidade Federal do Rio Grande do Sul. Nesta Tese é adotada uma metodologia alternativa para a definição das zonas rígidas do escoamento como sendo a posição onde a taxa de cisalhamento é igual a um valor dado pela relação de parâmetros reológicos do fluido, especificamente a tensão limite de escoamento e a viscosidade newtoniana para baixas taxas de cisalhamento. Nas simulações numéricas realizadas, o tempo de relaxação adimensional, o número de salto, o coeficiente power-law, a vazão adimensional e a massa específica adimensional são variados de forma a avaliar de que modo influenciam na dinâmica de escoamentos elastoviscoplásticos. Os resultados obtidos estão qualitativamente de acordo com a literatura, atestando a estabilidade da formulação empregada. / Non-Newtonian fluids are the majority of liquids found in nature and the study of their rheological behavior has a significant importance on different areas of engineering. Among them, there is a class of materials that exhibits little apparent deformation when subjected to a stress level behind an yield stress, referenced as viscoplastic material. In this class of materials, some fluids also exhibit elastic behavior at low shear rates. The present work aimed to a numerical study of two-dimensional steady state laminar flows of elasto-viscoplastic fluids through a planar expansion-contraction cavity. The mechanical model was defined by the mass conservation and momentum balance equations coupled to the elasto-viscoplastic model porposed in this work. This modeling has been approximated by a stabilized multi-field finite element method based on the Galerkin least-squares methodology, having as primal variables the elastic extra-stress component, velocity and pressure fields. In this way, the compatibility conditions between the extra-stress-velocity and pressure-velocity (Babuška- Brezzi condition) finite element subspaces are violated, allowing to use equal-order finite element interpolations. The stabilized method has been implemented in the finite element code for non-Newtonian fluids under development at the Laboratory of Applied and Computational Fluid Mechanics (LAMAC) of the Federal University of Rio Grande do Sul. An alternative methodology is adopted to define the yield surface as the position where the strain rate is equal to a value given by the relation of the rheological parameters of the fluid, namely the yield stress and the viscosity at low shear rates. In the performed numerical simulations, the non-dimensional relaxation time, the jump number, the power-law coefficient, the non-dimensional flow rate, and the non-dimensional density are varied in order to evaluate their influence on the elasto-viscoplastic fluid dynamics. All results found are in qualitatively accordance with the affine literature, and attesting the good stability features of the formulation.

Este trabalho tem como propósito investigar através do método Design Construtal a influência do espaçamento (S0) entre cilindros de seção transversal elíptica na maximização da densidade de transferência de calor em um escoamento externo sob efeito de convecção forçada. A razão de aspecto (r) entre os eixos vertical e horizontal dos cilindros elípticos também é um parâmetro avaliado. O estudo proposto é assumido bidimensional, incompressível, laminar e permanente. O regime de escoamento é dirigido por uma diferença de pressão ΔP, que se mantém através do domínio e é governada pelo número de Bejan (Be). Foram avaliados escoamentos com quatro diferentes números de Bejan, Be = 102, 103, 5.103, 104. O fluido que escoa através do domínio possui as propriedades termofisicas definidas pelo número de Prandtl, Pr = 0,72. O método Design Construtal, associado à busca exaustiva, foi empregado para determinar as restrições, graus de liberdade e o objetivo na avaliação geométrica do sistema A solução numérica das equações de conservação de massa, quantidade de movimento e energia foram resolvidas baseadas no método de volumes finitos, através do código comercial de dinâmica dos fluidos computacional FLUENT®. As geometrias e malha do domínio computacional foram desenvolvidas no pacote GAMBIT®. Como resultados, obteve-se que os casos ótimos apresentaram resultados consideravelmente melhores do que as demais configurações, onde se obteve ganhos de desempenho na densidade de transferência de calor de 50% a 97% em relação às configurações de menor desempenho estudadas. Além disso, foi possível verificar que o sistema adapta sua geometria ótima para cada condição de escoamento, a fim de proporcionar a melhor arquitetura de fluxo para atender ao objetivo térmico de maximizar a transferência de calor em um menor espaço físico, atendendo os princípios da Teoria Construtal. / This work investigates, through the Construtal Design method, the influence of the spacing (S0) between cylinders with elliptic cross in the maximization of the heat transfer density in an external flow with forced convection. The aspect ratio (r) between the vertical and horizontal axes of the elliptical cylinders is also evaluated. The proposed study is assumed twodimensional, incompressible, laminar and permanent. The flow regime is directed by a pressure difference ΔP, which is governed by the Bejan number (Be). The flows were evaluated for different values of the Bejan number, Be = 102, 103, 5.103, 104. The fluid flowing through the domain has its thermophysical properties defined by Prandtl number, Pr = 0.72. The Construtal Design method, associated with the exhaustive search, was used to determine the restrictions, degrees of freedom and objective in the geometric evaluation of the system. The numerical solution of the mass conservation, momentum and energy equations is solved based on the finite volume method, using the commercial fluid dynamics software FLUENT ®. The geometries and mesh of the computational domain were developed in the GAMBIT® package. The results show that the optimal cases performs considerably better than the other configurations, with an increase in the heat transfer density of 50% to 97% in comparison to the performance of lower level cases studied. In addition, it was possible to verify that the system adapts its optimal geometry to every flow condition in order to provide a better flow architecture that meets the thermal objective of maximizing a heat transfer in a smaller physical space in agreement with the principles of the Constructal Theory.

O coletor solar de tubos evacuados une uma alta absortividade de radiação solar a um ótimo grau de isolamento térmico. Estas características, aliadas a um custo relativamente baixo, fazem deste tipo de coletor o mais utilizado no mundo. Por isso, diversos tipos de abordagens estão sendo adotadas para descrever seu comportamento térmico. Nesta linha, este trabalho foi desenvolvido através de abordagem experimental e teórica, sendo a última subdividida em numérica por CFD e analítica. A abordagem experimental contou com a construção de uma bancada para medida de temperatura e radiação, em um coletor de 24 tubos evacuados acoplados em um reservatório de 178 L. A abordagem por CFD utilizou um modelo tridimensional transiente. Com o modelo numérico validado, utilizando resultados experimentais, simularam-se diferentes condições de operação, em termos de ângulo de inclinação, fluxo de calor sobre o coletor, tamanho do reservatório e temperatura de entrada da água. Em relação à abordagem analítica, esta é dividida em: modelo de irradiância, modelo do tubo e modelo do reservatório. O modelo de irradiância determina a distribuição da radiação solar ao longo da circunferência do tubo. Parte desta radiação é absorvida pelo coletor e transferida para água. Esta quantidade é determinada com o modelo do tubo, que é baseado no método de resistências térmicas. O modelo do reservatório descreve o comportamento térmico da água em seu interior, tanto em aquecimento quanto em resfriamento, analisando a interação energética com o coletor e com o meio externo. O desenvolvimento do modelo do tubo passa pela avaliação da vazão mássica entre o tubo e o reservatório, além disso, o coeficiente de transferência de calor por convecção no interior do tubo deve ser determinado. Tais variáveis são determinadas a partir de uma correlação para o número de Reynolds, a qual foi obtida com resultados da avaliação por CFD e é função dos números de Nusselt, Grashof e Prandtl. A correlação proposta apresentou bom ajuste com os resultados numéricos. Com a bancada de ensaio experimental foram feitas medidas de temperatura da água no reservatório ao longo de alguns dias. Para as mesmas condições do experimento, a temperatura média da água no reservatório foi estimada com resultados da integração dos modelos de irradiância, do tubo e do reservatório. A diferença entre os resultados experimental e teórico foi de 4,8% para a energia acumulada. / The evacuated tube solar collector combines high solar radiation absorptivity to a great thermal insulation degree. These characteristics, combined with a relatively low cost, make this type of collector the most used in the world. Therefore, various types of approaches are being adopted to describe its thermal behavior. In this way, this work was developed through experimental and theoretical approaches, the latter being subdivided into numeric, by CFD, and analytical analysis. For the experimental approach a test bench was built. The tests was carried on a solar collector with 24 evacuated tubes coupled to a 178 L tank, measuring temperature and solar radiation. The CFD approach used a transient three-dimensional model. After the numerical model validation using experimental data, simulations was carried over different operating conditions in terms of angle, heat flux on the collector, tank size and water inlet temperature. The analytical approach is divided into: irradiance model, tube model and tank model. The irradiance model determines the irradiance distribution of solar radiation along the circumference of the tube. Part of this radiation is absorbed by the collector and transferred to water, this amount is determined with the tube model, using the thermal resistance method. The tank model describes the thermal behavior of inside water, both in heating and in cooling, analyzing energy interaction with the collector and the external environment. The development of the tube model involves the assessment of the mass flow rate between the tube and the tank, furthermore the convection heat transfer coefficient inside the tube must be determined. These variables are determined from a correlation for the Reynolds number, which was obtained with evaluation results by CFD. Proposed Reynolds number is a function of the Nusselt, Prandtl and Grashof numbers. The correlation presented a good agreement with the numerical results. Using the experimental test bench the water temperature was measured into the tank over a few days. For the same experiment conditions, the average temperature of the water into the tank was estimated by results of integration of irradiance, tube and tank models. The difference between the experimental and theoretical results was 4.8% for the accumulated energy.

Publicación a texto completo no autorizada por el autor / El documento digital no refiere aseso / Calcula y proyecta un sistema de protección contra incendio a base de agua y espuma en la zona de dique y tanques de almacenamiento de hidrocarburos del terminal de almacenamiento. Determina el caudal y presión de la bomba contra incendio principal. Determina los diámetros de la red de tuberías del sistema de contra incendio por medio de los cálculos hidráulicos. Determina el volumen de reserva de agua en la cisterna para abastecer el sistema contra incendio. / Trabajo de suficiencia profesional

A simulação numérica de escoamento tem se desenvolvido em diferentes aspectos para melhorar sua exatidão e se aproximar cada vez mais da representação de fenômenos reais. O desenvolvimento científico sobre escoamentos turbulentos multifásicos dependem da simulação numérica, onde o custo computacional destas ainda é um fator determinante, principalmente quando se visa resultados mais precisos. Uma das grandes dificuldades é que códigos numéricos de alta ordem, utilizados para simular escoamentos bifásicos, produzem erros de dispersão numérica presentes principalmente na interface, já que nesta posição se apresenta uma descontinuidade nas propriedades físicas dos fluidos envolvidos. O presente trabalho tem o objetivo de propor uma formulação de alta precisão para escoamentos bifásicos. Para isso, realiza-se modificações em um código já existente, denominado Incompact3d que possui a possibilidade de paralelização para até 100:000 núcleos computacionais, utilizando uma decomposição 2D. Entre estas modificações destacam-se a implementação de uma nova formulação do método Level Set, além de pequenas modificações nas equações de Navier-Stokes, como a adição do termo de gravidade, a segregação do termo de difusividade e a adição de um termo referente a força de tensão superficial. Ainda propõem-se uma nova consideração do termo de pressão, que é separado em pressão não hidrostática e pressão hidrostática. Diversos testes de verificação e validação foram realizados, apresentando a capacidade de se realizar simulações numéricas de alta ordem para escoamentos bifásicos. Com esta metodologia, foi identificar características de fenômenos turbulentos bifásicos. / Flow numerical simulation has been developed in di erent ways to improve its results accuracy and have a better approximation to real phenomenon. To scientific development turbulent two-phase flow depends of numerical simulation, where computational cost is an important factor, mainly when accurate results are required. High order two-phase flow codes issue is that they produce numerical dispersion errors that are amplified with the presence of interface, since in its position exist physical proprieties discontinuity. This research main objective is to propose a high precision formulation to two-phase flow. For that, we modified Incompact3d code that has a 2D decomposition making parallelized calculation up to 100:000 computational cores. Some modifications are the implementation of a new Level Set Method formulation; some minor Navier-Stokes equation modifications as the addition of gravity term; the segregation of di usive term; and the addition of surface tension forces. We also proposed a new pressure term consideration where we separate it in hydrostatic and non-hydrostatic pressure. Many verification and validation tests where developed presenting the ability to perform high order numerical simulations for two-phase flows. Also it was possible to identify multiphase turbulent structures in the phenomena simulated.

Correntes de densidade são fluxos gravitacionais gerados pela diferença de densidade entre dois fluidos. Correntes de turbidez fazem parte de uma sub-classificação das correntes de densidade, na qual o fluido mais denso tem, na sua composição, partículas em suspensão. Muitos trabalhos numéricos já estudaram a dinâmica das correntes de turbidez, mas, nenhum dos encontrados aplicou mudanças de relevo concomitantes com a simulação, causadas pela sedimentação das próprias partículas da corrente e nem alterou o relevo após a passagem de cada evento em um domínio tridimensional. O presente trabalho pretende analisar a alteração no relevo de fundo causada por uma corrente de turbidez. No código Incompact3d, as equações de Navier-Stokes, Continuidade e Transporte e Difusão são resolvidas em uma malha cartesiana tridimensional. A condição inicial adotada é a de Lock-Exchange. As simulações realizadas utilizaram Simulação Numérica Direta (DNS). O código utiliza um esquema compacto centrado de sexta ordem, em diferenças finitas, para o esquema espacial, e Adams-Bashfort de terceira ordem para o esquema temporal. A validação do código foi realizada comparando-se com trabalhos experimental e numéricos. A análise das diferentes proporções granulométricas mostrou que quanto maior é a quantidade de material grosso na condição inicial, maior será seu depósito para um dado tempo. Em consequência, mais relevante se torna a consideração da alteração do relevo de fundo. Além disso, quanto maior o fator de compactação do sedimento, maior será o erro de não considerar a atualização de fundo. Os resultados também apontaram que os erros médios ao não considerar a atualização do fundo são da ordem de 4% da massa de depósito em 20 tempos adimensionais, para os parâmetros utilizados. Ao se propagar uma corrente de turbidez sobre o depósito de outra, os erros se mostram menores. / Gravity currents are gravitational fluxes triggered by density di erence between two fluids. A sub-classification of those are turbidity currents, in which the denser fluid is composed by the lighter fluid plus suspended particles. Many papers had shown turbidity currents dynamics, although none of the papers found had applied changes in the simulated topography due to deposit during the own simulation, neither they had altered a 3D domain topography after each flux, applying the changes caused by the previous current. The present dissertation aims to analyse the turbidity current dynamics alteration caused by the influence of its own deposit, altering the topography during the very simulation. The analysis is conducted in a polidispersed turbidity current. The Incompact3d code solves Navier-Stokes, continuity and transport-di usion equation, in a tridimensional cartesian mesh. Lock-exchange was chosen to be the initial condition. Direct Numerical Simulations (DNS) are performed. Sixth order compact finite-di erence schemes are used on the spatial domain, while third order Adams-Bashfort is applied for the temporal evaluation. Comparisons with numerical and experimental papers were performed for code verification. Results showed the coarser the particles on the starting lock-exchange, the higher its deposit is, and the more the terrain will be altered. Nevertheless, the bigger the compacting factor, the bigger the error of not considering bathymetry alteration. Results also point that the average errors of not considering the update are in order of 4% on the mass deposit, after 20 dimensionless times, for the used parameters. When a current propagates over the deposit of a previous one, these errors are smaller.

Simulações numéricas do escoamento atmosférico em microescala constituem o foco principal deste estudo. Estas simulações são abordadas tendo em vista aplicações para o setor eólico, em especial para avaliações de produção de energia em parques eólicos. Existem diversas categorias de incertezas associadas às estimativas de produção de energia para um projeto eólico, mas na maioria dos casos, a incerteza associada ao modelo de escoamento é a maior e mais relevante de todas. Dentro do setor eólico, o termo “modelo de escoamento” refere-se à ferramenta numérica utilizada para extrapolar o recurso eólico medido na posição das torres anemométricas (e sensores remotos) até as posições projetadas para os aerogeradores. Diversos autores sugerem através de estudos comparativos que os modelos tipo “CFD RANS k-ε” atualmente representam o “estado da arte” para aplicações em parques eólicos e são os mais utilizados comercialmente no setor. Contudo, o escoamento atmosférico livre é intrinsicamente turbulento, e a dinâmica dos escoamentos turbulentos é um campo científico que ainda não foi totalmente dominado pelo conhecimento humano. O presente estudo demonstra que a maioria dos “modelos de escoamento” atualmente disponíveis possuem pontos fracos, em especial quando aplicados em simulações do escoamento atmosférico livre sobre áreas com topografia e rugosidade complexas Uma das fraquezas presentes na maioria dos modelos de microescala para escoamento atmosférico é a “incapacidade” de simular com precisão o escoamento que ocorre durante períodos de “estabilidade atmosférica”. Diversos locais com elevado potencial eólico apresentam ciclos durante os quais as características do escoamento são afetadas pela ocorrência de estratificação térmica dentro da Camada Limite Atmosférica. Tendo como objetivo principal melhorar as simulações do escoamento nestas condições, propõe-se através deste estudo algumas modificações na modelagem “CFD RANS k-ε” tradicionalmente empregada. Dentre estas, destacam-se a inclusão de um perfil estratificado de temperatura potencial como condição de contorno, a inclusão dos efeitos das forças de empuxo no equacionamento “k-ε” e a solução simultânea das equações para balanço de energia e para o fluxo de temperatura potencial. Este modelo foi chamado de “RANS estável”. Para validação deste modelo foram utilizadas cinco torres anemométricas instaladas em um local com topografia complexa. Estas torres foram montadas e instrumentadas conforme as melhores práticas internacionais Os dados anemométricos registrados por essas torres demonstram a presença de ciclos diários de estabilidade atmosférica. Os erros de previsão cruzada foram calculados comparando-se as previsões de cada modelo com as medições reais registradas na posição das torres. O erro global médio de previsão cruzada entre torres anemométricas obtido com a composição dos modelos RANS “estável + neutro” foi de 3,8% enquanto o erro obtido apenas com o modelo RANS k-ε tradicional foi de 5,2%. Para o modelo linear WAsP, amplamente utilizado no setor eólico, o erro foi de 7,1%. Além dos erros de previsão cruzada entre torres, os perfis verticais de velocidade e os fatores de aceleração direcionais obtidos com a composição dos modelos RANS “estável + neutro” também sugerem que esta é uma alternativa versátil e promissora para capturar os ciclos de estabilidade atmosférica utilizando simulações numéricas em regime permanente. / Microscale numerical simulations of the atmospheric wind flow are the central focus of this study. These simulations are analysed from the wind energy perspective. Special attention is given to the usage and application of these simulations in energy production assessments for proposed wind farms. There are multiple uncertainty categories associated with energy production forecasts for future wind farms. However, in most cases the uncertainty factors related with wind flow modelling are the largest and most relevant of them all. The wording “flow model” refers to the numerical simulations (or “models”) that are used to extrapolate the anemometric data recorded at meteorological masts positions to the proposed wind turbine positions. Several authors have demonstrated through comparative studies that the “CFD RANS k-ε” models currently represent the “state of the art” when it comes to microscale wind flow simulations targeted at wind farms. Nonetheless, the atmospheric wind flow is turbulent by nature, and the dynamics of turbulent flows represent one of the scientific fields that have not yet been fully dominated by the human knowledge. The present study demonstrates that the majority of flow models currently available to mankind still lack in precision, even more so when it comes to modelling free atmospheric wind flow over complex terrain. One of the major weak spots of most microscale wind flow models is their inability to precisely simulate the wind flow that occurs during periods of atmospheric stability Numerous locations with large potential for wind energy production present cyclic periods of thermal stratification inside the atmospheric boundary layer. These cycles alter the dynamics and characteristics of the wind stream. With the purpose of improving wind flow simulations under stable atmospheric conditions, some modifications to the standard “RANS k-ε” model implementation are proposed. The most significant of these modifications are the usage of a potential temperature profile among the boundary conditions, the inclusion of the buoyancy forces in the “k-ε” equations and the simultaneous solution of the equations for energy balance and for potential temperature transport. This “modified” model was named “stable RANS”. It was validated using five well mounted meteorological masts installed in a location with complex topography. The anemometric data measured by these site masts suggest the existence of strong daily cycles of atmospheric stability. Cross prediction errors were calculated by comparing the forecasts (outputs) from each flow model against real wind data measured at each mast position The global average cross prediction error yielded by the RANS “stable +neutral” model was around 3,8%, whereas the error yielded by the traditional “RANS k-ε” implementation was near 5,2%. For the linear model WAsP the error was calculated to be 7,1%. In addition to cross prediction errors, the vertical wind speed profiles and speed-up factors calculated with the RANS “stable +neutral” model composition also suggest that it is a promising and versatile alternative for capturing the effects from atmospheric stability on wind flow using steady state numerical simulations.

A furação é um processo caracterizado por apresentar dificuldades em relação à formação e remoção do cavaco da região de corte, e à geração de altas temperaturas devido ao cisalhamento e encruamento do material no fundo do furo, principalmente, na usinagem do aço inoxidável austenítico. Com isso a presença do fluido de corte com boas propriedades de refrigeração e lubrificação faz-se necessária para preservar a vida da broca e garantir um acabamento satisfatório do furo. Contudo, com a tendência mundial de diminuição da quantidade de fluido de corte dispendido, tem-se buscado técnicas diferentes de aplicação que utilize o lubrirrefrigerante em quantidades reduzidas. Associado a essas novas técnicas, vem sendo empregado cada vez mais o uso de partículas sólidas adicionadas ao fluido com o intuito de aumentar sua eficiência. Um tipo de partícula, que vem sendo aplicada em diversos setores produtivos devido a suas ótimas propriedades mecânicas, químicas e elétricas, é o grafeno. Sendo assim, pretende-se analisar o comportamento da furação do aço inoxidável austenítico AISI 304 utilizando flocos de multicamadas de grafeno dispersos em fluido de corte empregando o projeto de experimentos Box-Behnken A aplicação do fluido é feita externamente utilizando a técnica de quantidade reduzida com três vazões: 1,5 l/h; 2,0 l/h e 2,5 l/h. O comportamento do processo é avaliado pelos esforços de corte (força de avanço e momento torsor), pela rugosidade (média e total), e pelos desvios (dimensional e de circularidade) na entrada e na saída do furo. O desempenho do grafeno é avaliado pela comparação dos resultados da sua utilização com os de outras duas condições de aplicação do fluido de corte sem adições: em quantidade reduzida e em abundância. Teve-se como resultado que todas as variáveis respostas analisadas sofreram alguma influência da velocidade de corte e/ou da vazão do fluido aplicado em quantidades reduzidas. Os benefícios da lubrificação e refrigeração do grafeno são sentidos em alguns resultados, como nos menores valores médios de rugosidade e de desvio dimensional. Mas, na análise do processo como um todo, não se afirma que o fluido de corte com flocos de multicamadas de grafeno aplicado externamente em quantidade reduzida traz mais benefícios que o sem adições. / Drilling is a machining process characterized by difficulties in the formation and removal of the chip from the cutting region and the generation of high temperatures due to shearing and hardening of the material at the bottom of the hole, especially in the machining of austenitic stainless steel. Thus the presence of cutting fluid with good cooling and lubrication properties is necessary to preserve the drill life and ensure a satisfactory hole finish. However, with the worldwide trend of decreasing the amount of cutting fluid expended, different application systems have been sought which use the coolant in reduced amounts. Associated with these new techniques, the presence of solid particles added to the fluid has been increasingly used in order to enhance its efficiency. One of these particles, which have been applied in several productive sectors due to its excellent mechanical, chemical and electrical properties, is graphene. Therefore, it is intended to analyze the drilling behavior of AISI 304 austenitic stainless steel using multilayer graphene flakes dispersed in the cutting fluid using the Box- Behnken Design The application of the fluid is done externally by means of reduced quantity lubricant with three flows: 1.5 l/h, 2.0 l/h and 2.5 l/h. The behavior of the process is evaluated via thrust force and torque, average and total roughness, and dimensional and circularity deviations at the input and output holes. The performance of graphene is evaluated by comparing the results of its use with two other cutting fluid conditions: quantity reduced and abundance. It was found that all the analyzed responses variables had some influence of the cutting speed and/or the flow of the applied fluid in reduced quantities. The benefits of graphene lubrication and cooling are felt in some results, such as the lower values of average roughness and dimensional deviation. However, in the analysis of the process as a whole, it is not possible to state that the multilayer graphene flakes cutting fluid externally applied brings more benefits than without additions.

A crescente demanda de energia dos dias de hoje, aliada aos graves problemas ambientais advindos da sua produção tem estimulado a procura de formas alternativas de energia. Dentre as diversas fontes de energia, as renováveis ganham expressão particularmente significativa. Neste contexto, a energia eólica revela-se uma aposta correta, a nível econômico, energético e ambiental. Por outro lado, a microprodução é apontada como uma das soluções para os problemas supracitados. Assim, a combinação de ambas as áreas tem sido alvo de especial interesse. Este trabalho visa caracterizar o vento em ambiente urbano, com particular interesse na possibilidade de seu aproveitamento energético. Um dos principais aspectos para determinar a viabilidade técnica e econômica de uma instalação eólica é a avaliação precisa da distribuição das velocidades de vento na área de aproveitamento. A instalação de turbinas eólicas em áreas urbanas tem determinado a necessidade de aprimorar a metodologia de previsão do campo de velocidades do vento visando à melhor determinação da distribuição dos equipamentos e aproveitamento do potencial existente. Neste contexto, esta dissertação apresenta um estudo sobre a aplicação da Dinâmica dos Fluidos Computacional - CFD para avaliação do potencial eólico e o comportamento do vento sobre um modelo de área urbana. Resultados numéricos com diferentes alternativas de modelagem do problema são comparados com dados de um experimento em túnel de vento e dados obtidos com CFD publicados na literatura, visando determinar uma metodologia adequada para avaliação do problema proposto. As simulações numéricas do escoamento de ar em área urbana são realizadas com o uso do programa ANSYS Fluent 13.0, que utiliza o método de volumes finitos para a solução das equações de Navier-Stokes com médias de Reynolds (RANS) e Simulação de Grandes Escalas (LES). O estudo é dividido em dois casos. No primeiro caso, a simulação da área em questão é realizada com as médias de Reynolds de forma a validar o modelo criado. No segundo caso, a simulação é realizada com o método de grandes escalas, cujos resultados serão explorados e comparados com o método RANS. Os resultados obtidos mostram essencialmente duas realidades, a primeira diz respeito à intensidade do recurso eólico e a segunda a sua particular complexidade. / The growing demand for today’s energy, combined with the serious environmental problems of its production, has stimulated the demand for alternative forms of energy. Among the various forms of energy, those that are based on renewable sources have been gaining a particularly significant space for discussion. In this context, the wind energy has proved to be a reliable source, energetically, economically and environmentally speaking. Furthermore, the micro production has been appointed as one of the solutions for the aforementioned problems. Thus, the combination of both these areas has particular interest. This work aims to characterize the wind in the urban environment, with particular interest for the possibility to carry out its energy. In this context this work presents a study on the application of computational fluid dynamics to evaluate the wind potential and the behavior of the wind on urban area. Numerical results with different alternatives for modeling the problem are compared with data from an experiment in wind tunnel and CFD simulation to determine one appropriate methodology for evaluation of the problem. The numerical simulations of the air flow over urban area are performed using the ANSYS Fluent 13.0 which uses the finite volumes method for solving the Reynolds Averaged Navier-Stokes (RANS) and Large-Eddy Simulation (LES). The study is divided in two cases. In the first one, the simulation of the area is performed with the Reynolds Averaged Navier-Stokes to validate the model created. In the second case, the simulation is performed with the Large-Eddy Simulation, which will be explored and compared with the results with RANS method. The results obtained show essentially two realities, the first is related to the intensity of wind resource and the second due to its high complexity. It was noted that in the urban areas turbulence intensity is significantly higher than in the rural environment. The non-horizontal flow components showed a high dependence on geometry of the building. The brief comparison of results obtained by simulation software shows an acceptable approach to the wind velocity field. In terms of turbulence intensity, it’s over estimated by the software. It is concluded from the results obtained that the energy use in urban environment provides a delicate task, due to the low wind intensity and complexity of the resource.