Neste trabalho, são inicialmente discutidas as dificuldades referentes à obtenção de resultados numéricos para a operação de uma turbina Savonius independentes do grau de discretização, tamanho de domínio de cálculo e de máximo tempo físico simulado. Também são relatadas as divergências entre as metodologias numéricas e experimentais adotadas por diversos autores, que dificultam análises e comparações dos resultados obtidos por meio dessas metodologias com os resultados próprios obtidos. Devido a esses fatos, no presente trabalho, uma série de procedimentos experimentais e numéricos são realizados para efetuar análises do escoamento sobre uma turbina eólica Savonius. Nos experimentos em canal aerodinâmico, perfis de velocidade e parâmetros da turbulência são obtidos pela técnica de anemometria de fio quente. Medições com o uso de tubos de Pitot e manômetros eletrônicos são efetuadas para avaliar a variação da pressão e os perfis de velocidade média em posições selecionadas. Além de dados para análise, informações úteis para uso como condições de contorno nas simulações numéricas também são obtidas. Os fenômenos são reproduzidos através de simulações numéricas pelo Método de Volumes Finitos, que solucionam as equações da continuidade, de Navier-Stokes com médias de Reynolds e do modelo de turbulência k-ω SST. Análises experimentais e numéricas considerando o escoamento sobre um cilindro, que mantém semelhanças com o escoamento sobre o rotor, também são realizadas. Simulações numéricas do escoamento sobre o cilindro são efetuadas, fornecendo resultados representativos do escoamento real, quando geometrias tridimensionais são aplicadas na modelagem numérica. Nas simulações do escoamento sobre o rotor Savonius em condição estática, resultados representativos do escoamento real são obtidos com o uso de uma modelagem que leva em consideração a rugosidade das pás do rotor, estacionado na posição angular de 90°. Para posições angulares menores, não se obteve uma boa concordância entre os resultados experimentais e numéricos. A realização deste trabalho fornece informações úteis para a análise do fenômeno e tem potencial para contribuir com futuros trabalhos desse tema. / This research work initially presents a discussion about the difficulties related to obtaining numerical results for the operation of a turbine Savonius independent of the degree of discretization, calculation domain size and maximum physical time of the simulation. The differences between the numerical and experimental methodologies adopted by various authors difficult the analysis and comparisons of the results obtained through these methods with the results obtained by the methodology. Due to these facts, in this research work, a series of experimental and numerical procedures are performed to conduct analyzes of flow over a Savonius wind turbine. In the experiments on aerodynamic channel, velocity and turbulence profiles parameters are obtained by the technique of hot wire anemometry. Measurements using Pitot tubes and electronic manometers are made to evaluate the variation of pressure and mean velocity profiles at selected positions. In addition to data analysis, useful information for use as boundary conditions in the numerical simulations are also obtained. The phenomena are reproduced through numerical simulations by Finite Volume Method, that solve the equations of continuity, Reynolds-averaged Navier–Stokes equations and the equation of the turbulence model k-ω SST. Experimental and numerical analyzes considering the flow over a cylinder, which holds similarities with the flow over the rotor, are also performed. Numerical simulations of the flow over the cylinder are made, providing results representative of the actual flow when three-dimensional geometries are applied in numerical modeling. In flow simulations over the Savonius rotor in static condition, at 90°, representative results of the actual flow are obtained using a model that takes into account the roughness of the rotor blades. For smaller angular positions a good agreement between experimental and numerical results was not obtained. This work provides useful information for the analysis of the phenomenon and has the potential to contribute to future work on this theme.

Polímeros sintéticos ocupam uma posição de grande importância no estilo de vida moderno, servindo como matérias-primas para a construção de uma variedade de utensílios. Apesar do grande número de operações de processamento e produtos disponíveis, o planejamento de produtos e a otimização dos processos de produção raramente constituem-se de tarefas triviais. Isso deve-se ao fato da maioria das operações aplicadas na indústria de processamento de polímeros envolverem geometrias e padrões de escoamento complexos, além da dificuldade intrínseca relacionada ao comportamento reológico complexo de polímeros fundidos ou soluções poliméricas. Devido a estes fatores, o desenvolvimento de técnicas de dinâmica de fluido computacional (computational fluid dynamics – CFD) para a simulação de escoamentos de fluidos poliméricos e etapas de operações de processamento tem sido assunto de numerosos estudos durante as últimas décadas. Sob esta perspectiva, o solver viscoelasticFluidFoam, merece destaque. Ele é capaz de resolver simulações de escoamentos de fluidos viscoelásticos utilizando diferentes equações constitutivas. Contudo, apesar de resultados existentes na literatura apresentarem um bom potencial de aplicação, uma análise extensiva de seu desempenho numérico ainda não foi realizada. Neste contexto, a proposta do presente trabalho é a análise da influência de parâmetros de malha, numéricos e constitutivos no comportamento do solver. As bases para os testes compreendem uma geometria simples – escoamento laminar entre duas placas paralelas – o modelo constitutivo de Oldroyd-B e respectivas soluções analíticas para os campos de velocidade e tensão. Mesmo os testes demonstrando a inegável versatilidade do solver, eles revelam limitações em lidar com algumas configurações de malha e parâmetros constitutivos, principalmente com relação ao refinamento na direção perpendicular ao escoamento, diminuição do número de Reynolds e aumento do número de Weisenberg. Estas limitações podem ser parcialmente contornadas com escolha adequada de parâmetros de relaxação das variáveis e da razão de aspecto dos volumes de controle. Tais dificuldades não estão presentes em simulações de escoamentos de fluidos newtonianos em condições semelhantes, sugerindo que trabalhos futuros devem focar em implementações mais robustas do solver viscoelasticFluidFoam. / Synthetic polymers hold a position of great importance in modern lifestyle, serving as raw materials for the construction of a wide variety of appliances. Despite the large number of processing operations and products available, product planning and optimization of production processes rarely constitute a trivial task. This is due to the fact of operations applied in polymer processing industry involve complex geometries and flow patterns, plus the intrinsic difficulty related to the molten polymers or polymer solutions complex rheological behavior. Because of these factors, the development of techniques of computational fluid dynamics (CFD) for the simulation of flows of polymeric fluids and stages of processing operations has been the subject of numerous studies during the last decades. From this perspective, the viscoelasticFluidFoam solver deserves mention. The solver is capable of resolving simulations of viscoelastic fluid flows using different constitutive equations. However, despite the existing results in the literature present a great potential for application, an extensive analysis of their numerical performance has not been performed yet. The purpose of this paper is to examine the influence of mesh, numerical and constitutive parameters in the behavior of the solver. Bases for the tests comprise a simple geometry – laminar flow between two parallel plates – the constitutive model of Oldroyd-B and its analytical solutions for the velocity and stress fields. Although the tests show the undeniable versatility of the solver, they also reveal limitations in dealing with some mesh settings and constitutive parameters, particularly with respect to refinement in the direction perpendicular to the flow, decreasing in the Reynolds number and increasing in the Weisenberg number. This limitation can be partially circumvented with proper choice of variables relaxation parameters and aspect ratio of the control volumes. Such difficulties are not present in simulations of Newtonian fluids flows under similar conditions, suggesting that future works should focus on more robust implementations of the viscoelasticFluidFoam solver.

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.

O presente trabalho apresenta a modelagem do escoamento de ar sobre um terreno complexo empregando, e comparando os resultados, três programas comerciais de previsão de potencial eólico, WaSP, Meteodyn WT e WindSim. Por se tratar de um modelo simplificado, o programa WaSP nem sempre é adequado para previsão de energia em terrenos muito complexos pois não é capaz de prever turbulência, separação de escoamento e efeitos de estabilidade existentes neste tipo de terreno. Para que se consiga modelar a turbulência existente de maneira mais coerente, utilizam-se programas de mecânica dos fluidos computacional como o Meteodyn WT e o WindSim que modelam o escoamento através das equações médias de Reynolds. A finalidade do estudo é identificar a capacidade de previsão das velocidades de vento por meio de cada um dos programas em um parque eólico localizado na Chapada da Diamantina, tendo em vista que a correta previsão da produção de energia depende de um entendimento correto e detalhado dos recursos disponíveis. Avaliando apenas condição atmosférica neutra e comparando os resultados obtidos, conclui-se que o programa WindSim obteve o menor erro na estimativa da velocidade do vento, -3,7% quando comparado aos outros programas WaSP e Meteodyn, -5,1% e -6,5% respectivamente. O maior erro de aproximação foi obtido pelo programa WaSP, +11% quando comparado aos outros programas Meteodyn e WindSim, +10,6% e +7,9% respectivamente. Para o programa Meteodyn foram avaliados casos com diferentes condições de estabilidade, o erro máximo e mínimo foi de +10,6% e -6,5%, respectivamente em atmosfera estável, mantendo-se elevado. O programa WindSim permite a alteração de variáveis que são fixas no Meteodyn e também a simulação dos efeitos de estabilidade através da inclusão da equação de energia inicializada através do comprimento de Monin-Obukhov o que, implicou em uma melhora significativa nos resultados obtidos e em um maior tempo computacional. O erro, menor e maior, apresentado pelo programa WindSim foi reduzido a, -1,8% e +4,4% respectivamente. Destaca-se que as conclusões obtidas se aplicam ao caso estudado. / This work studies the air flow modeling over a complex terrain using and comparing the results of three commercial programs that estimates the wind field in a site: WaSP, Meteodyn WT and WindSim. Being WAsP a simplified model, it is not always suitable to estimate the flow in complex terrain because it is not able to predict the turbulence, flow separations and stability effects present in these kind of terrains. To model the existing turbulence two computational fluid mechanics programs i.e. Meteodyn WT and WindSim, which models the flow through the Reynolds Averaged Navier-Stokes equations have been used. This study aims to identify the wind speed predictive capacity of each program in a windfarm located on Chapada da Diamantina, as a correct production prediction depends on a correct estimation of the wind field. Evaluating only neutral atmospheric setup and comparing it with the measures, WindSim had the lowest wind speed estimation error, -3.7% compared with WaSP and Meteodyn programs, that provided -5.1% and -6.5% respectively. The greatest error was obtained by WaSP program, +11%, followed by Meteodyn, +10.6% and WindSim, +7.9%. WindSim enables changes in some variables that are fixed in Meteodyn and enables also the simulation of stability effects through the introduction of the energy equation with a Monin-Obukhov length initialization. For Meteodyn program different stabilities were evaluated, but the results obtained were not satisfactory. The maximum and minimum error were reduced to +10.6% and -6.5%, respectively, in stable atmosphere. Due to the specific site studied conditions the simulation of stability effects resulted in a significant improvement in the results even though a larger computational time. Enabling the energy equations setting up stable atmospheric stability, both WindSim errors are reduced to -1.8% and +4.4% respectively. It is noteworthy that the conclusions apply to the case studied.

Orientador: Prof. Dr. Rovilson Mafalda / Dissertação (mestrado) - Universidade Federal do ABC, Programa de Pós-Graduação em Engenharia e Gestão da Inovação, 2017. / Presente em vários países, a competição Baja SAE se consolidou como uma importante oportunidade dos estudantes de engenharia de colocar em prática os conhecimentos adquiridos em sala de aula e desenvolver novas habilidades. Apesar disso, após vinte e três edições no Brasil, o cenário atual da competição é de utilização quase exclusiva de soluções de projeto consagradas, grande semelhança entre os veículos, acirramento da disputa e carência de soluções inovadoras. Neste trabalho propõe-se uma abordagem focada na aerodinâmica como uma forma de se buscar uma vantagem competitiva no atual cenário da competição, através da elaboração e disponibilização de um pacote composto de carenagens aerodinâmicas aplicáveis em um veículo tipo Baja SAE típico. Adicionalmente, o impacto da adoção destas peças na velocidade máxima do veículo é avaliado através de análise de dinâmica dos fluidos computacional (do inglês CFD ¿ Computational Fluid Dynamics). De acordo com as estimativas realizadas neste trabalho, a adoção de um pacote aerodinâmico como o proposto seria capaz de reduzir o arrasto aerodinâmico do veículo em 20,45% e aumentar a sua velocidade máxima em 7,92%. / Present in several countries, the Baja SAE competition has consolidated itself as an important opportunity for engineering students to put into practice the knowledge acquired in the classroom and to develop new skills. Despite this, after twenty-three editions in Brazil, the current scenario of the competition is the consolidation of consecrated solutions, great similarity between vehicles, escalation of the dispute and lack of innovative solutions. This work proposes an approach focused on aerodynamics as a way to seek a competitive advantage in the current scenario, through the elaboration and provision of a package composed of aerodynamic fairings possible to be applied in a typical Baja SAE prototype. Additionally, the impact of the adoption of these parts on the top speed of the vehicle is evaluated through Computational Fluid Dynamics (CFD) analysis. According to the estimates made in this work, the adoption of an aerodynamic package like the one proposed would be able to reduce the aerodynamic drag of the vehicle by 20.45% and increase its maximum speed by 7.92%.

Neste trabalho, são inicialmente discutidas as dificuldades referentes à obtenção de resultados numéricos para a operação de uma turbina Savonius independentes do grau de discretização, tamanho de domínio de cálculo e de máximo tempo físico simulado. Também são relatadas as divergências entre as metodologias numéricas e experimentais adotadas por diversos autores, que dificultam análises e comparações dos resultados obtidos por meio dessas metodologias com os resultados próprios obtidos. Devido a esses fatos, no presente trabalho, uma série de procedimentos experimentais e numéricos são realizados para efetuar análises do escoamento sobre uma turbina eólica Savonius. Nos experimentos em canal aerodinâmico, perfis de velocidade e parâmetros da turbulência são obtidos pela técnica de anemometria de fio quente. Medições com o uso de tubos de Pitot e manômetros eletrônicos são efetuadas para avaliar a variação da pressão e os perfis de velocidade média em posições selecionadas. Além de dados para análise, informações úteis para uso como condições de contorno nas simulações numéricas também são obtidas. Os fenômenos são reproduzidos através de simulações numéricas pelo Método de Volumes Finitos, que solucionam as equações da continuidade, de Navier-Stokes com médias de Reynolds e do modelo de turbulência k-ω SST. Análises experimentais e numéricas considerando o escoamento sobre um cilindro, que mantém semelhanças com o escoamento sobre o rotor, também são realizadas. Simulações numéricas do escoamento sobre o cilindro são efetuadas, fornecendo resultados representativos do escoamento real, quando geometrias tridimensionais são aplicadas na modelagem numérica. Nas simulações do escoamento sobre o rotor Savonius em condição estática, resultados representativos do escoamento real são obtidos com o uso de uma modelagem que leva em consideração a rugosidade das pás do rotor, estacionado na posição angular de 90°. Para posições angulares menores, não se obteve uma boa concordância entre os resultados experimentais e numéricos. A realização deste trabalho fornece informações úteis para a análise do fenômeno e tem potencial para contribuir com futuros trabalhos desse tema. / This research work initially presents a discussion about the difficulties related to obtaining numerical results for the operation of a turbine Savonius independent of the degree of discretization, calculation domain size and maximum physical time of the simulation. The differences between the numerical and experimental methodologies adopted by various authors difficult the analysis and comparisons of the results obtained through these methods with the results obtained by the methodology. Due to these facts, in this research work, a series of experimental and numerical procedures are performed to conduct analyzes of flow over a Savonius wind turbine. In the experiments on aerodynamic channel, velocity and turbulence profiles parameters are obtained by the technique of hot wire anemometry. Measurements using Pitot tubes and electronic manometers are made to evaluate the variation of pressure and mean velocity profiles at selected positions. In addition to data analysis, useful information for use as boundary conditions in the numerical simulations are also obtained. The phenomena are reproduced through numerical simulations by Finite Volume Method, that solve the equations of continuity, Reynolds-averaged Navier–Stokes equations and the equation of the turbulence model k-ω SST. Experimental and numerical analyzes considering the flow over a cylinder, which holds similarities with the flow over the rotor, are also performed. Numerical simulations of the flow over the cylinder are made, providing results representative of the actual flow when three-dimensional geometries are applied in numerical modeling. In flow simulations over the Savonius rotor in static condition, at 90°, representative results of the actual flow are obtained using a model that takes into account the roughness of the rotor blades. For smaller angular positions a good agreement between experimental and numerical results was not obtained. This work provides useful information for the analysis of the phenomenon and has the potential to contribute to future work on this theme.

Grades da tomada dágua de usinas hidroelétricas são equipamentos de grande importância porque são responsáveis pela proteção das turbinas hidráulicas contra impacto de corpos flutuantes. Tais estruturas estão submetidas à ação de cargas dinâmicas oriundas do escoamento da água através das barras. O objetivo do presente trabalho é analisar grades submetidas à ação de fluxo de água. Ou seja, analisar as respostas da estrutura e o comportamento do escoamento da água utilizando os cálculos de dinâmica de estruturas acopladas às técnicas de dinâmica de fluidos computacional (CFD). Para a elaboração destas análises, foi utilizado o software comercial CFX versão 14. Tais análises foram elaboradas mediante o processo de interação fluido-estrutura. Primeiramente, um modelo estrutural simplificado das barras verticais das grades é elaborado a partir de dados de projetos conhecidos. A partir deste modelo define-se um volume de controle que representa o escoamento do fluido. Devido ao número de Reynolds calculado, utiliza-se o modelo de turbulência para obtenção dos resultados tais como tensões e deslocamentos nas barras verticais, e perfil de velocidades do escoamento. Considerando o modelo proposto, foram elaboradas análises com sete diferentes valores de velocidade de modo a comparar os dados calculados numericamente através do método de elementos finitos com valores obtidos experimentalmente. A partir do modelo verificado, é apresentada uma análise de uma grade inclinada submetida a um fluxo paralelo. / Trashracks are very important equipment because they are responsible for protecting turbines of hydroelectric plants against floating bodies. These structures are subjected to the action of dynamic loads due to the water flow through the vertical and horizontal bars. The objective of this study was to analyze trashracks submitted by action of water flow. In other words, to analyze the responses of the structure and the behavior of water flow using dynamic of structures calculations coupled with computational fluid dynamics techniques (CFD) for a turbulent regime, through the use of commercial software CFX version 14. This analysis is elaborated by the process of fluid structure interaction. First of all, a simplified structural model of vertical bars is defined from other similar projects. For this model is defined a volume of control that represents fluid flow. Due to the Reynolds number calculated, it is utilized a turbulence model in order to obtain the results. These results are: stresses and displacements of vertical bars; and profile of velocities of flow. The results are analyzed and discussed. After that, considering the simplified model, analyzes with seven different values of speed are executed in order to compare the results between the data calculated numerically by finite elements method, and the values obtained experimentally. Considering the model verified, an analysis, of an inclined trashrack subjected to a parallel flow, is presented.

O presente trabalho apresenta um método de interação fluido estrutura chamado de fronteira virtual para modelagem de escoamentos sobre geometrias complexas. A principal característica do método é que a condição de contorno interna da superfície do corpo imerso é obtida usando uma função que calcula a força de campo das equações de momentum. Este método de cálculo da força de campo é conhecido como método direto. A principal vantagem do método da fronteira virtual é que se trabalha com duas malhas, uma cartesiana (também conhecida como euleriana) para o fluido e outra curvilínea (também conhecida como malha lagrangiana) para o objeto imerso. Como estas malhas são independentes uma da outra, os algoritmos de geração das malhas são simples em comparação com outros métodos de geração de malha, nos quais a esta precisa se adaptar ao corpo que se quer analisar. Nos métodos de body fitted, algoritmos complexos com custo computacional muito elevado são necessários para se trabalhar com geometrias complexas. Neste trabalho, o método de fronteira virtual é desenvolvido para trabalhar escoamentos compressíveis não viscosos em duas dimensões sobre geometrias complexas, e testado sobre um cilindro de seção circular e sobre uma série de aerofólios NACA da série 4. O regime de principal interesse é o regime transônico, em particular para velocidades abaixo da velocidade do som (Mach entre 0,8 e 1). No entanto, resultados para outros regimes de escoamento (subsônico incompressível, subsônico compressível, transônico e supersônico) também são mostrados. / This work presents a fluid-structure interaction method, known as virtual boundary method, developed for flow modeling, over complex geometries. The main characteristic of the method is that, the internal boundary condition at the surface are created using a function, that compute the body force from the momentum equations. This method is known as direct method. The virtual boundary method main advantage is the easy and quick mesh generation. Due the use two different grids, a Cartesian grid for the fluid and a curvilinear grid for the body that interact with the fluid. As these two grids are independent one of other, the grid generation algorithms are easier compared with other methods on CFD as the body fitted. Where the grid has to be adapted to the body of interest. The body fitted methods are implemented with complex algorithms, and high computational cost, especially when complex geometries are analyzed. In this work, the virtual boundary method are developed in order to work with in viscid compressible flow in two dimensions over complex geometries, tested over a cylinder and several NACA series 4 and critical airfoils. The regime of interest is the transonic, below the speed of sound, nevertheless results for several flow regimens (incompressible subsonic, compressible subsonic, transonic and supersonic) are also presented.

Com o atual desenvolvimento industrial e tecnológico da sociedade, a presença de substâncias inflamáveis e/ou tóxicas aumentou significativamente em um grande número de atividades. A possível dispersão de gases perigosos em instalações de armazenamento ou em operações de transporte representam uma grande ameaça à saúde e ao meio ambiente. Portanto, a caracterização de uma nuvem inflamável e/ou tóxica é um ponto crítico na análise quantitativa de riscos. O objetivo principal desta tese foi fornecer novas perspectivas que pudessem auxiliar analistas de risco envolvidos na análise de dispersões em cenários complexos, por exemplo, cenários com barreiras ou semi-confinados. A revisão bibliográfica mostrou que, tradicionalmente, modelos empíricos e integrais são usados na análise de dispersão de substâncias tóxicas / inflamáveis, fornecendo estimativas rápidas e geralmente confiáveis ao descrever cenários simples (por exemplo, dispersão em ambientes sem obstruções sobre terreno plano). No entanto, recentemente, o uso de ferramentas de CFD para simular dispersões aumentou de forma significativa. Estas ferramentas permitem modelar cenários mais complexos, como os que ocorrem em espaços semi-confinados ou com a presença de barreiras físicas. Entre todas as ferramentas CFD disponíveis, consta na bibliografia que o software FLACS® tem bom desempenho na simulação destes cenários. Porém, como outras ferramentas similares, ainda precisa ser totalmente validado. Após a revisão bibliográfica sobre testes de campo já executados ao longo dos anos, alguns testes foram selecionados para realização de um exame preliminar de desempenho da ferramenta CFD utilizado neste estudo. Foram investigadas as possíveis fontes de incertezas em termos de capacidade de reprodutibilidade, de dependência de malha e análise de sensibilidade das variáveis de entrada e parâmetros de simulação. Os principais resultados desta fase foram moldados como princípios práticos a serem utilizados por analistas de risco ao realizar análise de dispersão com a presença de barreiras utilizando ferramentas CFD. Embora a revisão bibliográfica tenha mostrado alguns dados experimentais disponíveis na literatura, nenhuma das fontes encontradas incluem estudos detalhados sobre como realizar simulações de CFD precisas nem fornecem indicadores precisos de desempenho. Portanto, novos testes de campo foram realizados a fim de oferecer novos dados para estudos de validação mais abrangentes. Testes de campo de dispersão de nuvem de propano (com e sem a presença de barreiras obstruindo o fluxo) foram realizados no campo de treinamento da empresa Can Padró Segurança e Proteção (em Barcelona). Quatro testes foram realizados, consistindo em liberações de propano com vazões de até 0,5 kg/s, com duração de 40 segundos em uma área de descarga de 700 m2. Os testes de campo contribuíram para a reavaliação dos pontos críticos mapeados durante as primeiras fases deste estudo e forneceram dados experimentais para serem utilizados pela comunidade internacional no estudo de dispersão e validação de modelos. Simulações feitas utilizando-se a ferramenta CFD foram comparadas com os dados experimentais obtidos nos testes de campo. Em termos gerais, o simulador mostrou bom desempenho em relação às taxas de concentração da nuvem. O simulador reproduziu com sucesso a geometria complexa e seus efeitos sobre a dispersão da nuvem, mostrando claramente o efeito da barreira na distribuição das concentrações. No entanto, as simulações não foram capazes de representar toda a dinâmica da dispersão no que concerne aos efeitos da variação do vento, uma vez que as nuvens simuladas diluíram mais rapidamente do que nuvens experimentais. / With the industrial and technological development of the present-day society, the presence of flammable and toxic substances has increased in a growing number of activities. Dispersion of hazardous gas releases occurring in transportation or storage installations represent a major threat to health and environment. Therefore, forecasting the behaviour of a flammable or toxic cloud is a critical challenge in quantitative risk analysis. The main aim of this dissertation has been to provide new insights that can help technological risks analysts when dealing with complex dispersion modelling problems, particularly those problems involving dispersion scenarios with barriers or semi-confined. A literature survey has shown that, traditionally, empirical and integral models have been used to analyse dispersion of toxic/flammable substances, providing fast estimations and usually reliable results when describing simple scenarios (e.g. unobstructed gas flows over flat terrain). In recent years, however, the use of CFD tools for simulating dispersion accidents has significantly increased, as they allow modelling more complicated gas dispersion scenarios, like those occurring in complex topographies, semi-confined spaces or with the presence of physical barriers. Among all the available CFD tools, FLACS® software is envisaged to have high performance when simulating dispersion scenarios, but, as other codes alike, still needs to be fully validated. This work contributes to the validation of FLACS software for dispersion analysis. After a literature review on historical field tests, some of them have been selected to undertake a preliminary FLACS performance examination, inspecting all possible sources of uncertainties in terms of reproducibility capacity, grid dependence and sensitivity analysis of input variables and simulation parameters. The main outcomes of preliminary FLACS investigations have been shaped as practical guiding principles to be used by risk analysts when performing dispersion analysis with the presence of barriers using CFD tools. Although the literature survey has shown some experimental data available, none of the works include detailed exercises giving new insights of how to perform accurate CFD simulations nor giving precise rates of FLACS performance. Therefore, new experiments have been performed in order to offer new sets of cloud dispersion data for comprehensive validation studies. Propane cloud dispersion field tests (unobstructed and with the presence of a fence obstructing the flow) have been designed and undertaken at Can Padró Security and Safety training site (Barcelona) by which intensive data on concentration has been acquired. Four tests were performed, consisting on releases up to 0.5 kg/s of propane during 40 seconds in a discharge area of 700 m2. The field tests have contributed to the reassessment of the critical points raised in the guiding principles and have provided experimental data to be used by the international community for dispersion studies and models validation exercises. FLACS software has been challenged against the experimental data collected during the field tests. In general terms, the CFD-based simulator has shown good performance when simulating cloud concentration. FLACS reproduces successfully the presence of complex geometry and its effects on cloud dispersion, showing realistic concentration decreases due to cloud dispersion obstruction by the existence of a fence. However, simulated clouds have not represented the whole complex accumulation dynamics due to wind variation, since they have diluted faster than experimental clouds.

A contribuição efetiva do presente trabalho de tese consiste na implementação de esquemas numéricos, MacCormak (1969) e Jameson e Mavriplis (1986), para resolução das equações de Euler no espaço tridimensional , segundo o contexto de volumes finitos e utilizando técnicas estruturadas e não estruturadas de discretização espacial. A aplicação direta seria para simulações de problemas físicos envolvendo, principalmente, escoamentos supersônicos e hipersônicos de "gás frio" em torno de configurações complexas. Em paralelo a este esforço, são apresentados técnicas de aceleração de convergência em uso comum pela comunidade de CFD e alguns operadores de dissipação artificial para esquemas numéricos com discretização espacial simétrica, com o propósito de averigüar suas características no tocante ao desempenho computacional e à qualidade global da solução.