Ce travail a été consacré à la simulation numérique des équations de la mécanique des fluides par des méthodes de volumes finis implicites. Tout d’abord, nous avons étudié et mis en place une version implicite du schéma de Roe pour les écoulements monophasiques et diphasiques compressibles. Grâce à la méthode de Newton utilisée pour résoudre les systèmes nonlinéaires, nos schémas sont conservatifs. Malheureusement, la résolution de ces systèmes est très coûteuse. Il est donc impératif d’utiliser des algorithmes de résolution performants. Pour des matrices de grande taille, on utilise souvent des méthodes itératives dont la convergence dépend de leur spectre. Nous avons donc étudié le spectre du système linéaire et proposé une stratégie de Scaling pour améliorer le conditionnement de la matrice. Combinée avec le préconditionneur classique ILU, notre stratégie de Scaling a réduit de façon significative le nombre d’itérations GMRES du système local et le temps de calcul. Nous avons également montré l’intérêt du schéma centré pour la simulation de certains écoulements à faible nombre de Mach. Nous avons ensuite étudié et implémenté la méthode de décomposition de domaine pour les écoulements compressibles. Nous avons proposé une nouvelle variable interface qui rend la méthode du complément de Schur plus facile à construire et nous permet de traiter les termes de diffusion. L’utilisation du solveur itératif GMRES plutôt que Richardson pour le système interface apporte aussi une amélioration des performances par rapport aux autres méthodes. Nous pouvons également découper notre domaine de calcul en un nombre quelconque de sous-domaines. En utilisant la stratégie de Scaling pour le système interface, nous avons amélioré le conditionnement de la matrice et réduit le nombre d’itérations GMRES de ce système. En comparaison avec le calcul distribué classique, nous avons montré que notre méthode est robuste et efficace. / This thesis deals with numerical simulations of compressible fluid flows by implicit finite volume methods. Firstly, we studied and implemented an implicit version of the Roe scheme for compressible single-phase and two-phase flows. Thanks to Newton method for solving nonlinear systems, our schemes are conservative. Unfortunately, the resolution of nonlinear systems is very expensive. It is therefore essential to use an efficient algorithm to solve these systems. For large size matrices, we often use iterative methods whose convergence depends on the spectrum. We have studied the spectrum of the linear system and proposed a strategy, called Scaling, to improve the condition number of the matrix. Combined with the classical ILU preconditioner, our strategy has reduced significantly the GMRES iterations for local systems and the computation time. We also show some satisfactory results for low Mach-number flows using the implicit centered scheme. We then studied and implemented a domain decomposition method for compressible fluid flows. We have proposed a new interface variable which makes the Schur complement method easy to build and allows us to treat diffusion terms. Using GMRES iterative solver rather than Richardson for the interface system also provides a better performance compared to other methods. We can also decompose the computational domain into any number of subdomains. Moreover, the Scaling strategy for the interface system has improved the condition number of the matrix and reduced the number of GMRES iterations. In comparison with the classical distributed computing, we have shown that our method is more robust and efficient.

Esta tese apresenta os resultados da utilização de superfícies com microcanais como método de controle passivo de escoamento visando à redução do arrasto turbulento. Obtém-se essa redução pelo aumento da anisotropia das tensões de Reynolds junto à parede, situação em que o estado da turbulência se torna localmente axissimétrico e a uma componente. Utilizam-se as equações de Navier-Stokes para formular os escoamentos e o modelo de transporte de Tensões de Reynolds para computar as quantidades turbulentas. Aproximam-se essas equações pelo Método dos Volumes Finitos e a solução numérica é computada com o solucionador Star-CCM+ v. 8.06. Propõe-se um modelo de predição de redução de arrasto para uma placa plana com microcanais de seção retangular e dimensões geométricas variáveis, aplicando a formulação e o método numérico para calcular a anisotropia das tensões de Reynolds, o estado da turbulência e a redução de arrasto relativamente a uma placa de superfície lisa e de mesma área molhada. No capítulo de análise e discussão emprega-se o triângulo de Lumley-Pope para determinar o estado da turbulência e a trajetória de retorno à isotropia do escoamento, após verificar a incerteza numérica e validar o resultado com o modelo de predição e o da teoria da placa plana. Conduzem-se as análises quantitativas examinando as tensões de cisalhamento, as tensões de Reynolds e a morfologia do escoamento em pontos do interior e em superfícies lisas e adjacentes aos microcanais. Conclui-se apresentando uma visão geral dos resultados e propondo alternativas de desdobramento e continuidade deste trabalho. / The potential reduction of turbulent drag is investigated for flows over a flat plate with streamwise aligned microgrooves. For this purpose, the connection between the anisotropy of the Reynolds stresses and drag reduction effect is presented, and a model is developed in order to estimate the drag reduction potential according to flow and geometrical settings. The Navier-Stokes transport equations particularized for incompressible flows are used to describe the fluid motion, and the turbulence quantities are evaluated using the linear pressure-strain Reynolds stress transport model. The quantities are estimated using the Finite-Volume Method, which is applied to a set of grids with different refinement levels and groove topologies. After validating the numerical results against the predictions of the proposed model, and the theoretical estimates available in the literature, the author discusses the drag reducing effect by examining the state of turbulence in the microgrooves, also providing an assessment on the anisotropy of the Reynolds Stresses inside, near and outside the grooves. In the final chapter, conclusions are drawn, and outlooks of possible extensions to this work are suggested.

Orientador: Kátia Lucchesi Cavalca Dedini / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-17T15:01:37Z (GMT). No. of bitstreams: 1 Machado_TiagoHenrique_M.pdf: 6488903 bytes, checksum: f0285c2b84eaaa9effae0b3cce163537 (MD5) Previous issue date: 2011 / Resumo:Mancais hidrodinâmicos são amplamente utilizados em máquinas rotativas, sendo o elemento responsável pela interação entre o rotor e a estrutura de fundação. Portanto, a fim de descrever o comportamento dinâmico dos eixos em rotação, é necessário conhecer as características desses mancais. Por esta razão, a proposta deste trabalho é analisar as características dinâmicas e operacionais de diferentes configurações de mancais radiais hidrodinâmicos. O objetivo deste projeto é desenvolver e implementar uma representação consistente para o modelo matemático que descreve o comportamento do mancal hidrodinâmico, através da equação de Reynolds e, assim, obter a distribuição da pressão no mancal, bem como os coeficientes dinâmicos, a fim de avaliar o desempenho de diferentes geometrias em condições severas, altas velocidades de rotação e alta carga aplicada. Um modelo numérico para a solução da equação de Reynolds, com a ajuda da equação de Bernoulli, obtém o campo de pressão gerado dentro do filme de óleo. Os coeficientes dinâmicos são obtidos a partir de um conceito semelhante ao de um sistema mola-amortecedor, a fim de representar a flexibilidade e o amortecimento inerentes ao filme de óleo. A abordagem numérica utilizada neste estudo é feita através do método dos volumes finitos. A escolha deste método foi feita devido à sua ampla utilização na simulação numérica dos fenômenos da dinâmica dos fluidos e transferência de calor. O trabalho apresenta uma análise para três diferentes geometrias de mancais: um cilíndrico, um elíptico e um trilobular. Além das diferentes geometrias, o projeto também propõe um estudo sobre o efeito de algumas descontinuidades geométricas no comportamento dos mancais. As descontinuidades foram introduzidas na forma de ranhuras axiais, e seu efeito é então comparado com o desempenho dos mancais sem a introdução destas descontinuidades / Abstract: Hydrodynamic bearings are widely used in rotating machines, being the element responsible for the interaction between the rotor and the supporting structure. Therefore, in order to describe the dynamic behavior of the rotating shafts, it is necessary to know the journal bearings characteristics. For this reason, the purpose of this work is to analyze the dynamic and operational characteristics of different configurations of radial hydrodynamic bearings. The objective of this project is to develop and implement a consistent mathematical representation for the model of the hydrodynamic bearing through the Reynolds equation and thereby obtain the pressure distribution in the bearing, as well as the dynamic coefficients, in order to evaluate the performance of different geometries in severe conditions, high speeds and high-applied load. A numerical model for the solution of Reynolds equation, with help from the Bernoulli equation, obtains the pressure field generated inside the oil film. For the dynamic coefficients, they are obtained from a similar concept to the spring-damper, in order to represent the inherent flexibility and damping of the oil film. The numerical approach used in this study is done using the finite volume method. The choice of this method was made due to their widespread use in the numerical simulation of phenomena of fluid dynamics and heat transfer. The work presents an analysis for three different geometries of bearings: a cylindrical, an elliptical and a three-lobe. In addition to the different geometries, the project also proposes a study on the effect of some geometric discontinuities in the behavior of the bearings. The discontinuities were introduced in the form of axial grooves, and its effect is then compared with the performance of bearings without these grooves / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador: Kátia Lucchesi Cavalca Dedini / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecânica / Made available in DSpace on 2018-08-17T14:52:53Z (GMT). No. of bitstreams: 1 Vieira_LeonardoCarpinetti_M.pdf: 3383681 bytes, checksum: 7b1d6a8f8eb035c38aaa391c784305bb (MD5) Previous issue date: 2011 / Resumo: A geração de pressão no fluido lubrificante presente na folga entre um mancal axial e o colar do eixo em rotação é de fundamental importância para evitar contato entre partes sólidas com movimento relativo. Qualquer contato existente poderá causar fricção, desgaste e, consequentemente, falha dos elementos de máquinas rotativas. Sendo assim, no projeto de um mancal eficaz, é importante conhecer o mecanismo pelo qual a pressão é gerada no filme de óleo e a magnitude da carga axial a ser transmitida do colar para o mancal através do fluido. Para isso, é necessário resolver a Equação de Reynolds, obtendo-se a distribuição de pressão no óleo. A partir disso, diversos parâmetros operacionais, tais como capacidade de carga, fluxos de fluido, posição do pico de pressão entre outros, podem ser obtidos. De maneira a avaliar o problema de lubrificação hidrodinâmica proposto, é obtida uma solução utilizando o Método dos Volumes Finitos em coordenadas polares. Utiliza-se este método para resolver o balanço de fluxo em cada volume de controle com auxílio da Equação de Bernoulli, permitindo a análise de descontinuidades do filme de óleo. Segmentos individuais e mancais compostos, incluindo as suas ranhuras, são analisados. A influência da gama de parâmetros envolvidos, como a área plana ao fim de cada segmento de mancal, espessura de fluido e viscosidade do fluido foi analisada, buscando-se valores ótimos destes parâmetros. Finalmente, foram calculados e analisados o coeficiente de rigidez e o coeficiente de amortecimento para este tipo de mancal / Abstract: The pressure generation within the lubricant fluid present in the clearance between a thrust bearing and the collar attached to the shaft has a fundamental importance to avoid contact between solid parts with axial relative motion. Any existing contact can lead to friction, wear and, as a consequence, failure of elements on a rotating machine. Therefore, in order to design an effective bearing, it is important to know how the pressure is generated in the oil film and the magnitude of the load capacity transmitted from the collar to the bearing throughout the fluid. Thus, it is necessary to solve the Reynolds' Equation to obtain the pressure distribution on the sections under Hydrodynamic Lubrication. Afterwards, several operational parameters can be obtained, such as the total load capacity, lubricant fluid flow, position of the maximum pressure and so on. In order to evaluate this hydrodynamic lubrication problem, the solution of the Reynolds Equation by using the Finite Volume Method in polar coordinates is obtained. This method is used to solve the balance of fluid flow in each control volume with help from the Bernoulli Equation, allowing the analysis of the film thickness discontinuities. Individual pads and complete axial bearings, including its grooves, are analyzed. The influence of the set of parameters involved, such as area of the flat part of the pad, film thickness and viscosity of the fluid, on the results of the distribution of pressure and the calculated load capacity was evaluated, as well as the evaluation of the optimum dimensions for bearings. Finally, the stiffness and damping coefficients for such type of bearings were obtained and the obtained results were analyzed / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador: Katia Lucchesi Cavalca / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-14T04:10:09Z (GMT). No. of bitstreams: 1 Ferreira_RafaelBittencourt_M.pdf: 2098484 bytes, checksum: e3c06e9d3444778c194664385fd94c5a (MD5) Previous issue date: 2009 / Resumo: Desde os primeiros estudos numéricos sobre lubrificação em mancais hidrodinâmicos até a atualidade, existiram relativamente poucos trabalhos abordando o problema através do método dos volumes finitos. Na verdade, a maioria dos estudos, na área de lubrificação, utiliza a equação de Reynolds e o método das diferenças finitas para estabelecer o campo de pressão gerado no filme de óleo, mesmo sendo este tipo de abordagem inviável em certos casos, pois para geometrias mais complexas, ou com certas condições dinâmicas do sistema, a qualidade da discretização e a consistência dos resultados gerados pela equação de Reynolds se mostram insuficientes. Sendo assim, este trabalho trata da análise do campo de pressão em um filme gerado através do efeito hidrodinâmico em um mancal radial cilíndrico plano. Para tanto, o fluido de lubrificação é colocado em condição isotérmica e em regime de escoamento laminar. Para avaliar o campo de pressão é desenvolvido um algoritmo que determina a solução através do método dos volumes finitos para a equação simplificada de Reynolds em um domínio com duas dimensões, e segue com as seguintes etapas: construção da malha, integração das equações dentro dos volumes, definição das equações discretizadas, estabelecimento de condições iniciais e de contorno e avaliação do campo de pressão no domínio de cálculo. Outro ponto analisado é a influência de ranhuras nos mancais, as quais trazem novos elementos para o comportamento do mancal, tais como a redução da pressão no filme de óleo e a maior retenção do fluido entre o eixo e o mancal. Ao final, é possível concluir que a grande importância desta abordagem está na similaridade do método dos volumes finitos com a tradicional representação diferencial do escoamento em volumes de controle, tornando mais claros os conceitos físicos envolvidos, outro ponto abordado é a comparação com o método das diferenças finitas e a solução analítica de Ocvirk, com esses resultados é possível avaliar as principais características e vantagens de cada método. Com respeito ao campo de pressão, também é possível chegar aos valores das forças de sustentação do eixo e ainda avaliar os coeficientes de rigidez do mancal, que são fundamentais no estudo da dinâmica de máquinas rotativas. / Abstract: Since the first studies on hydrodynamic journal bearings using numerical methods until today, there has been relatively few studies addressing the problem through the Finite Volume Method. Indeed, the majority of studies in the area of lubrication use the equation of Reynolds and method of finite differences to establish the field of pressure generated in the oil film, even if this type of approach is not viable in certain cases, specially for complex geometries or for certain dynamic conditions of the system. Thus, this work deals with the analysis of the pressure distribution in a film generated by the hydrodynamic effect in a radial radial journal bearing. Wherefore, the fluid lubrication is placed in isotherm condition and laminar flow. To evaluate the pressure distribution an algorithm is developed to determine the solution through the Finite Volume Method for the simplified Reynolds equation in a field with two dimensions, and follows with the steps: construction of the mesh, integration of equations within the volumes, the definition of discrete equations, establishment of initial conditions and boundary conditions and finally the evaluation of pressure distribution in the calculation domain. Another point analyzed is the influence of grooves in the bearings, these grooves bring new elements for the performance of journal bearings, such as to reduce the pressure in the oil film and the greater retention of the fluid between rotor and bearing, this study is highlight also the behavior of the stiffness of the oil film in certain conditions of operation, and this dynamic parameter has high importance in studies of rotative machines response. At the end, it is possible to conclude that the great importance of this approach is the similarity of the Finite Volume Method with the traditional differential representation in studies with fluid mechanics, making clearer the physical concepts involved, another point raised is the comparison with the Ocvirk's analytical solution, allowing to evaluate the main characteristics and the results quality for this numerical method. / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientador: Milton Mori / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-09T09:36:46Z (GMT). No. of bitstreams: 1 Vegini_AtilanoAntonio_D.pdf: 28017713 bytes, checksum: b5c4951104b5d515ebe084a6388a6ea2 (MD5) Previous issue date: 2007 / Resumo: Neste trabalho é apresentado o código computacional CYCLO, programa desenvolvido especificamente para predizer o escoamento multifásico em ciclones. Os ciclones são equipamentos utilizados na indústria, principalmente, para a separação de material particulado de correntes gasosas. Desde o início da industrialização até a atualidade, os ciclones são investigados cientificamente na tentativa de elucidar o real funcionamento deste complexo sistema de escoamento que ocorre no interior destes equipamentos. Com o avanço computacional, as modernas técnicas da fluidodinâmica computacional (CFD), têm se mostrado capaz de auxiliar nas investigações científicas sobre os ciclones. O código CYCLO, é um destes programas baseado nas técnicas da fluidodinâmica computacional (CFD). Consiste de um modelo multifásico num referencial euleriano, e composto de equações diferenciais parciais que são resolvidas numericamente através da aplicação do método dos volumes finitos, com acoplamento pressão-velocidade e malha numérica deslocada. A hipótese de simetria no eixo é assumida, tornando possível a utilização de um modelo 3DSimétrico em coordenadas cilíndricas. O código CYCLO torna possível à utilização de até três fases sólidas, cada uma caracterizada por seu tamanho de partícula, massa específica e fração volumétrica. O escoamento da fase sólida é considerado invíscido e a transferência de quantidade de movimento interfase predito pela adoção de um modelo de força de arraste. Para representar o comportamento anisotrópico dos tensores de Reynolds na fase gasosa, é utilizado um modelo de turbulência híbrido composto por uma combinação do modelo k-e padrão para os componentes radiais e axiais do tensor de Reynolds, e, para os componentes tangenciais por intermédio da teoria do comprimento de mistura de Prandtl. Os objetivos deste trabalho vão desde a verificação e validação do código CYCLO com seu modelo multifásico euleriano-euleriano e todas suas hipóteses simplificadoras, até a apresentação de aplicações práticas do uso das técnicas de CFD. Para a verificação e validação do código, os resultados obtidos numericamente são comparados com valores experimentais para o escoamento monofásico e para o escoamento multifásico. Como exemplo de aplicação prática das técnicas de CFD, é apresentado a utilização do código na análise e projeto de ciclones industriais / Abstract: In this work it is presented the computational code CYCLO, which is a program developed specifically to predict the dynamic behavior of the gas-solid flow in cyclones. The cyclones have been used extensively in many industries for a long time for separating the particles from gaseous streams. Since the beginning of the industrialization until the present time, the cyclones are investigated scientifically in order to elucidate the real functioning of this complex flow behavior inside of these equipments. With the computational advance, the modern techniques of the computational fluid dynamics (CFD), has been capable to assist in the scientific inquiries concerning cyclones. The CYCLO code is one of these programs based on computational fluid dynamics (CFD) techniques. The program consists of a multiphase flow model based on the Eulerian-Eulerian approach and it is composed of several partial differential equations, which are solved using the finite volume method with pressure-velocity coupling and staggered grids. The axial symmetry hypothesis is assumed, which makes possible to apply the 3D symmetric model in cylindrical coordinate. The CYCLO code makes possible the use of up to three solid phases, each one represented by a particle diameter, density and specific volumetric fraction. The solid-phase flow is considered inviscid and drag forces between phases are responsible for the gas-solid interaction. To represent an anisotropic behavior of the Reynolds stress in the gas phase, it is used a hybrid model composed of the combination of the (k-e) standard model for the radial and axial components and Prandt¿s longitudinal mixing model for the tangential components. The purpose of the present study is to verify and validate the model used by the CYCLO code and to show practical applications of the use of CFD techniques. In order to verify and validate the code, numerical results were compared with experimental data for gas and gas-solid flows. As practical examples of application of CFD techniques, the code was used to the analysis and design of industrial cyclones / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química

Orientadores: Maria Cristina de Castro Cunha, Persio Leister A. Barros / Dissertação (mestrado profissional) - Universidade Estadual de Campinas, Instituto de Matematica, Estatistica e Computação Cientifica / Made available in DSpace on 2018-08-11T16:59:15Z (GMT). No. of bitstreams: 1 Santos_UilbiranChaves_M.pdf: 537649 bytes, checksum: d5123c50fbd2d334c5620c343887ad48 (MD5) Previous issue date: 2008 / Resumo: Inicialmente é apresentado um capítulo sobre a modelagem matemática e sua importância nos dias atuais. Discutimos o papel das equações diferenciais, ordinárias e parciais, na representação matemática de processos reais que envolvam leis de conservação usadas na física dando enfoque às equações de difusão, que resultam do princípio de conservação de energia. No segundo capítulo é apresentado o método dos volumes finitos, uma técnica atual e bastante útil usada na discretização de equações diferenciais parciais parabólicas e hiperbólicas. A seguir o método dos volumes finitos é aplicado para encontrar aproximações para a solução da equação do calor transiente e espacialmente bidimensional. / Abstract: Initially it presented a chapter on mathematical modeling and its relevance today. We discussed the role of differential equations, ordinary and partial, the mathematical representation of actual cases involving conservation laws of physics focusing used in the equations of delivery, that result from the principle of conservation of energy. The second chapter is presented the method of finite volume, a current and very useful technique used in discretization of partial differential equations and parabolas hyperbolic. Following the method finite volume is applied to find approaches to the solution of the equation of heat and spatially transient two-dimensional. / Mestrado / Mestre em Matemática

[pt] Uma grande parcela da população depende de águas subterrâneas para o consumo. Assim, vazamentos dos tanques de gasolina representam um risco significativo de saúde pública, além dos riscos de incêndio e explosão que o vazamento de vapores destes tanques pode causar nas linhas de esgotos dos edifícios. Estes tanques estão em geral localizados na zona não saturada dos solos, acima do nível freático. Torna-se assim importante avaliar a mobilidade destes compostos nestas condições e em particular de prever a sua chegada no lençol freático. Esta avaliação pode ser feita através de procedimentos numéricos de análise. O principal objetivo deste trabalho é a implementação numérica de um programa que simule o fluxo transiente de NAPLs (non aqueous phase liquids) em meios porosos não saturados. A solução numérica é obtida pela utilização do método dos volumes finitos. Para representar as características hidráulicas dos solos, a função Brooks e Corey e os modelos de Brooks e Corey-Burdine são utilizados. O trabalho apresenta detalhes do programa desenvolvido e o seu uso para simular ensaios de fluxo de gasolina e gasolina misturada com etanol realizados na PUC em estudos anteriores. / [en] A large proportion of the population depends on groundwater for consumption. Thus, leakage of gasoline tanks represents a significant risk to public health besides the risk of fire and explosion in the sewer lines, due to the flammable gases. These tanks are generally located in the unsaturad soil above the groundwater level. It is therefore important to evaluate the mobility of these compounds under these conditions and in particular to predict its arrival in the water table. This assessment can be done by numerical procedures for analysis. The main objective of this work is the numerical implementation of a program that simulates the transient flow of NAPLs (non aqueous phase liquids) in unsaturated porous media. The numerical solution is obtained by using the method of finite volumes. To represent the hydraulic characteristics of soils, the function of Brooks and Corey, as well as the Brooks and Corey-Burdine’s model are used. The essay presents details of the program developed and its use to simulate the flow testing of gasoline and gasoline blended with ethanol made in the PUC in previus studies.

The objective of this dissertation is to develop robust and accurate numerical methods for solving the compressible, non-linear Euler equations of gas dynamics in one and two space dimensions. In theory, solutions of the Euler equations can display various characteristics including shock waves, rarefaction waves and contact discontinuities. To capture these features correctly, highly accurate numerical schemes are designed. In this thesis, two different projects have been studied to show the accuracy and utility of these numerical schemes. Firstly, the compressible, non-linear Euler equations of gas dynamics in one space dimension are considered. Since the non-linear partial differential equations (PDEs) can develop discontinuities (shock waves), the numerical code is designed to obtain stable numerical solutions of the Euler equations in the presence of shocks. Discontinuous solutions are defined in a weak sense, which means that there are many different solutions of the initial value problems of PDEs. To choose the physically relevant solution among the others, the entropy condition was applied to the problem. This condition is then used to derive a bound on the solution in order to satisfy L2-stability. Also, it provides information on how to add an adequate amount of diffusion to smooth the numerical shock waves. Furthermore, numerical solutions are obtained using far-field and no penetration (wall) boundary conditions. Grid interfaces were also included in these numerical computations. Secondly, the two dimensional compressible, non-linear Euler equations are considered. These equations are used to obtain numerical solutions for compressible ow in a shock tube with a 90° circular bend for two channels of different curvatures. The cell centered finite volume numerical scheme is employed to achieve these numerical solutions. The accuracy of this numerical scheme is tested using two different methods. In the first method, manufactured solutions are used to the test the convergence rate of the code. Then, Sod's shock tube test case is implemented into the numerical code to show the correctness of the code in both ow directions. The numerical method is then used to obtain numerical solutions which are compared with experimental data available in the literature. It is found that the numerical solutions are in a good agreement with these experimental results.

A representação de malhas por meio de estrutura de dados e operadores topológicos e um dos focos principais da modelagem geométrica, onde permite uma implementação robusta e eficiente de mecanismos de refinamento adaptativo, alinhamento de células e acesso as relações de incidência e adjacência entre os elementos da malha, o que é de grande importância na maioria das aplicações em mecânica dos fluidos. No caso de malhas não estruturadas, a não uniformidade da decomposição celular e melhor representada por uma estrategia mais sofisticada, que são as estruturas de dados topológicas. As estruturas de dados topológicas indexam elementos de uma malha representando relações de incidência e adjacência entre elementos, garantindo acesso rápido às informações. Um dos aspectos mais comuns aos problemas tratados pela mecânica dos fluidos computacional é a complexidade da geometria do domínio onde ocorre o escoamento. O uso de estruturas de dados para manipular malhas computacionais e de grande importância pois realiza de modo eficiente as consultas às informações da malha e centraliza todas as operações sobre a malha em um único módulo, possibilitando sua extensão e adaptação em diversas situações. Este trabalho visou explorar o acoplamento de uma estrutura de dados topológica, a Mate Face, em um módulo simulador existente, de modo a gerenciar todos os acessos à malha e dispor operações e iteradores para pesquisas complexas nas vizinhanças de cada elemento na malha. O módulo simulador resolve as equações governantes da mecânica dos fluidos através da técnica de volumes finitos. Foi utilizada uma formulação que atribui os valores das propriedades aos centroides dos volumes de controle, utiliza métodos de alta ordem, os esquemas ENO e WENO, que tem a finalidade de capturar com eficiência descontinuidades presentes em problemas governados por equações diferenciais parciais hiperbólicas. As equações de Euler em duas dimensões representam os escoamentos de interesse no presente trabalho. O acoplamento da estrutura de dados Mate Face ao simulador foi realizada através da criação de uma biblioteca desenvolvida que atua como uma interface de comunicação entre os dois módulos, a estrutura de dados e o simulador, que foram implementados em diferentes linguagens de programação. Deste modo, todas as funcionalidades existentes na Mate Face tornaram-se acessíveis ao simulador na forma de procedimentos. Um estudo sobre malhas dinâmicas foi realizado envolvendo o método das molas para movimentação de malhas simulando-se operações de arfagem. A idéia foi verificar a aplicabilidade deste método para auxiliar simulações de escoamentos não estacionarios. Uma outra vertente do trabalho foi estender a estrutura Mate Face de forma a representar elementos não suportados a priori, de modo a flexibilizar o seu uso em simulações de escoamentos baseados no método de volumes finitos espectrais. O método dos volumes espectrais e utilizado para se obter alta resolução espacial do domínio computacional, que também atribui valores das propriedades aos centroides dos volumes de controle, porém, os volumes de controle são particionados em volumes menores de variadas topologias. Assim, uma extensão da Mate Face foi desenvolvida para representar a nova malha para a aplicação do método, representando-se cada particionamento localmente em cada volume espectral. Para todas as etapas deste trabalho, realizaram-se experimentos que validaram a utilizaação da estrutura de dados Mate Face junto a métodos numéricos. Desta forma, a estrutura pode auxiliar as ferramentas de simulações de escoamentos de fluidos no gerenciamento e acesso à malha computacional / The storage and access of grid files by data structures and topologic operators is one of the most important goals of geometric modeling research field, which allows an efficient and stable implementation of adaptive refinement mechanisms, cells alignment and access to incidence and adjacency properties from grid elements, representing great concernment in the majority of applications from fluid mechanics. In the case of non-structured grids, the cellular decomposition if non-uniform and is better suited by a more sophisticated strategy - the topologic data structs. The topologic data structs index grid elements representing incidence and adjacency properties from grid elements, ensuring quick access to information. One of most common aspects from problems solved by computational fluid mechanic is the complexity of the domain geometry where the fluid ows. The usage of data structures to manipulate computational grids is of great importance because it performs efficiently queries on grid information and centers all operations to the grid on a unique module, allowing its extension and flexible usage on many problems. This work aims at exploring the coupling of a topologic data structure, the Mate Face, on a solver module, by controlling all grid access providing operators and iterators that perform complex neighbor queries at each grid element. The solver module solves the governing equations from fluid mechanics though the finite volume technique with a formulation that sets the property values to the control volume centroids, using high order methods - the ENO and WENO schemes, which have the purpose of efficiently capture the discontinuities appearing in problems governed by hyperbolic conservation laws. The two dimensional Euler equations are considered to represent the flows of interest. The coupling of the Mate Face data structure to the solver module was achieved by a creation of a library that acts as an interface layer between both modules, the Mate Face and the solver, which had been implemented using different programming languages. Therefore, all Mate Face class methods are available to the solver module though the interface library in the form of procedures. A study of dynamic grids was made by using spring methods for the moving grid under pitch movement case. The goal was to analyze the applicability of such method to aid non stationary simulations. Another contribution of this work was to show how the Mate Face can be extended in order to deal with non-supported types of elements, allowing it to aid numeric simulations using the spectral finite volume method. The spectral nite volume method is used to obtain high spatial resolution, also by setting the property values to the control volume centroids, but here the control volumes are partitioned into smaller volumes of different types, from triangles to hexagons. Then, an extension of the Mate Face was developed in order to hold the new generated grid by the partitioning specfied by the spectral finite volume method. The extension of Mate Face represents all partitioned elements locally for each original control volume. For all implementations and proposals from this work, experiments were performed to validate the usage of the Mate Face along with numeric methods. Finally, the data structure can aid the fluid flow simulation tools by managing the grid file and providing efficient query operators