[pt] Neste trabalho é apresentado o estudo da interação entre um jato supersônico e uma superfície plana, com o objetivo de analisar o comportamento do campo de velocidade, pressão e temperatura do escoamento. Este estudo encontra sua motivação no processo de descamação térmica de rochas duras, a qual pode resultar da iteração entre um jato a alta pressão e temperatura e a rocha. Este processo, que pode ser útil na perfuração de rochas duras e profundas, ocorre devido ao acúmulo de tensões térmicas na rocha, o qual pode acarretar sua fratura. Este tipo de processo também envolve diversos mecanismos aerodinâmicos e termodinâmicos, que são isoladamente fenômenos abertos. No desenvolvimento deste trabalho o escoamento foi modelado pelas equações de Navier - Stokes bidimensionais para uma mistura de gases perfeitos em um sistema de coordenadas cilíndrico. O modelo considerado para descrever o transporte turbulento é o modelo de uma equação de Spalart - Allmaras, o qual envolve a solução de uma equação diferencial para a viscosidade turbulenta. Estas equações são resolvidas utilizando-se uma metodologia de volumes finitos adaptada a escoamentos compressíveis. A descrição dos escoamentos transientes obtidos necessitou de diversas modificações ao código computacional existente. Estas modificações trataram, em particular, das condições de contorno, que utilizam a noção de características, e do modelo de turbulência. A estrutura do escoamento resultante da interação entre o jato supersônico e a parede é estudada, avaliando-se a influência (i) da distância entre a saída do jato e a parede, (ii) da razão de pressões entre o jato e o ambiente. Além disso, é examinada a evolução transiente do escoamento. Os resultados obtidos são analisados com vista a obter as melhores condições aerodinâmicas para o processo de descamação térmica. / [en] I in this work a study of the interaction between a supersonic jet and a planar surface is presented, with the aim to analyze the behavior of the velocity, pressure and temperature flowfield. This study finds its motivation in the process of thermal spallation of hard rocks, which may result from the interaction between a high pressure and high temperature jet and the rock. This process, that can be used in the drilling of hard and deep rocks, occurs due to the accumulation of thermal stresses in the rock, which can cause its fracture. This type of process also involves several aerodynamic and thermodynamic mechanisms, which are still open phenomena. In the development of this work the flow was modeled by the two-dimensional Navier-Stokes equation for a mixture of perfect gases in a cylindrical coordinates system. The model considered to describe the turbulent transport is the one equation of Spalart - Allmaras model, which involves the solution of a differential equation for the turbulent viscosity. These equations are solved using a finite volumes methodology which is adapted to compressible flows. The description of the obtained transient flow required several modifications in the existing computational code. These modifications involved, in particular, the choice of boundary conditions, that use the notion of characteristics, and the turbulence model. The structure of the flow resulting from the interaction between the supersonic jet and the wall is studied. In particular, are examined the influence (i) the distance between the jet and wall, (II) of the pressures ratio between the jet and the environment. Moreover, the transient evolution of the flow is examined. The obtained results are examined to determine the best aerodynamic conditions for the process of thermal spallation to occur.

Este trabalho tem como propósito investigar como a variação geométrica de determinados parâmetros envolvidos na construção de uma geometria bifurcada de seção circular, em forma de Y, afeta a resistência ao escoamento, tanto de fluidos newtonianos como não newtonianos. As geometrias estudadas foram construídas utilizando-se o princípio do Design Construtal. Os parâmetros variados foram a relação entre os comprimentos dos dutos pais e filhos, a relação entre os diâmetros dos mesmos dutos, e o ângulo central da estrutura em forma de Y. Para as relações geométricas lineares foram utilizados os valores de 0,5; 0,6; 0,7; 0,8; 0,9 e 1, enquanto para os ângulos, foram utilizados os valores de 155°, 135°, 115°, 95°, 75°, 45°, 25° e 10°. Os fluidos utilizados foram do tipo newtoniano e não newtoniano, dentre estes últimos, foram estudados fluidos pseudoplásticos e dilatantes. O trabalho foi realizado através de simulações numéricas, implementadas com a utilização do software comercial Ansys Fluent, o qual resolve as equações governantes através do método dos volumes finitos. As malhas utilizadas foram do tipo poliédrica. Os resultados indicam que há uma diferença em relação ao que se espera da literatura para as relações entre os diâmetros e os comprimentos. A Lei Hess-Murray indica que estas relações ótimas seriam de 2-1/3 para as relações entre os diâmetros e comprimentos. No presente trabalho, foram determinadas relações entre os diâmetros próximas de 0,6, e entre os comprimentos, iguais a 1. Os ângulos ótimos ficaram localizados no intervalo entre 100° e 135°. / This work aims to investigate how the geometric variation of certain parameters involved in the construction of a bifurcated Y-shaped circular cross-section geometry affects the flow resistance of both Newtonian and non-Newtonian fluids. The geometries studied were constructed using the Constructal Design principle. The parameters were the relationship between the lengths of the daughter and parent ducts, the relationship between the diameters of the same ducts, and the central angle of the Y-shaped structure. For the linear geometric relations, values of 0.5; 0.6; 0.7; 0.8; 0.9 and 1 where used, for the angles, the values of 155 °, 135 °, 115 °, 95°, 75 °, 45 °, 25 ° and 10 ° were used. The fluids used were of the Newtonian and non-Newtonian type, among the latter, pseudo plastic and dilatant fluids were studied. The work was carried out through numerical simulations, implemented with the commercial software Ansys Fluent, which solves the governing equations through the finite volume method. The meshes used were of the polyhedral type. The results indicate that there is a difference in relation to what is expected from the literature for the relationships between diameters and lengths. The Hess-Murray Law indicates that these optimal relations would be 2-1/3 for the relationships between diameters and lengths. In the present work, relationships between the diameters close to 0,6 were found and s equal to 1 between the lengths. The optimum angles were located in the range between 100 ° and 135 °.

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.

Neste trabalho desenvolveu-se a modelagem matemática e computacional na análise de Funções de Interpolação nos Métodos de Volumes Finitos. As equações foram integradas no volume de controle onde nas fronteiras devem ser avaliados os fluxos advectivos e difusivos. Para calcular estes fluxos, utilizamos a função de interpolação cujo, principal objetivo é conectar estes pontos nodais. Utilizamos três esquemas de função de interpolação: CDS, UDS, QUICK estes esquemas foram escolhidos, pois apresentaram os menores erros de truncamentos possíveis. Em síntese, os testes computacionais mostraram as limitações dos esquemas CDS e UDS, bem como também revelaram os ótimos resultados do esquema QUICK. / In this work, the mathematical and computer modeling in the analysis of the interpolation functions in the Finite Volume Method was developed. The equations were integrated in the control of volume in which frontiers the advective and diffuse flows must be evaluated. The interpolation function, which main objective is to connect the nodal points, was used to calculate those flows. Three interpolation function schemes were used: CDS, UDS and QUICK. These schemes were chosen due to the fact that they present less errors and blockings. To sum up, the computer tests showed the limitations of the CDS and UDS schemes and the excellent results of the QUICK scheme.

Made available in DSpace on 2015-05-08T14:59:37Z (GMT). No. of bitstreams: 1 arquivototal.pdf: 3644722 bytes, checksum: e20bf7ab6319ef9aed02f134cd90a808 (MD5) Previous issue date: 2012-01-23 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior - CAPES / Food drying is a process used to ensure the quality and stability of foods, considering that the decrease of water in the material reduces the biological activity and the chemical and physical changes that occur during storage. The objective of this research was to study the heat and mass transfer in drying process of slices of tomatoes and to evaluate the characteristics of drying kinetics in two resistive equipments and a heat pump. The operating conditions for dryers resistive were set at temperatures of 40 and 50 and 60°C with an air velocity of 2.5m/s to the heat pump and a dryer resistive temperature of 55°C and air velocity of 5.7m/s. The experimental results of drying kinetics were represented using mathematical models and are shown in terms of effective diffusivity and convective of the mass. The values found for the effective diffusivity of mass are in agreement with the literature. The convective coefficients mass are the order of 10-1. Those coefficients were not found references in literature. A comparative study between the heat pump and the dryer resistive showed that the first one was more efficient to perform the drying of tomato. We studied the sorption isotherm for tomatoes at temperatures of 50, 60, 70 and 80°C. The results were represented using mathematical models and the GAB model showed the best fit to the experimental data. A numerical study for heat transfer and mass during drying of the product was studied, through the development of a computational code, in configuring an infinite slab. In the modeling process were studied the Fourier Law for heat transfer and Law of Fick for mass transfer. For code validation were used results from the literature and from this present research. The results were satisfactory, but it is necessary to study the variables involved in the drying process for a better adjustment. / Secagem de alimentos é um processo utilizado para assegurar a qualidade e a estabilidade dos alimentos, considerando que a diminuição da quantidade de água do material reduz a atividade biológica e as mudanças químicas e físicas que ocorrem durante o armazenamento. O objetivo deste trabalho foi estudar a transferência de calor e massa na secagem de tomate cortado em fatias e avaliar as características da cinética de secagem em dois equipamentos resistivos e uma bomba de calor. As condições operacionais para os secadores resistivos foram fixadas nas temperaturas de 40, 50 e 60ºC, na velocidade do ar de 2,5m/s e para a bomba de calor e um secador resistivo na temperatura de 55ºC e velocidade do ar de 5,7m/s. Os resultados experimentais da cinética de secagem foram representados através de modelos matemáticos e são mostrados em termos de difusividade efetiva e convectiva de massa. Os valores encontrados para a difusividade efetiva de massa estão de acordo com os encontrados na literatura. Os coeficientes convectivos de massa estão na ordem de 10-1, não sendo encontrado referencial da literatura. Um estudo comparativo entre a bomba de calor e o secador resistivo e o primeiro se mostrou mais eficiente para realizar a secagem do tomate. Foi estudada a isoterma de sorção para o tomate nas temperaturas de 50, 60, 70 e 80ºC. Os resultados foram representados através de modelos matemáticos e o modelo de GAB mostrou o melhor ajuste aos dados experimentais. Um estudo numérico para a transferência de calor e massa durante a secagem do produto foi estudado, através do desenvolvimento de um código computacional, na configuração de uma placa plana infinita. Na modelagem foram estudadas a Lei de Fourier para a transferência de calor e a Lei de Fick para a transferência de massa. Para validação do código foram utilizados resultados obtidos na literatura e os do presente trabalho. Os resultados obtidos mostrou-se satisfatório, porém é necessário estudar as variáveis envolvidas na secagem em busca de uma melhor ajuste.

Neste trabalho desenvolveu-se a modelagem matemática e computacional na análise de Funções de Interpolação nos Métodos de Volumes Finitos. As equações foram integradas no volume de controle onde nas fronteiras devem ser avaliados os fluxos advectivos e difusivos. Para calcular estes fluxos, utilizamos a função de interpolação cujo, principal objetivo é conectar estes pontos nodais. Utilizamos três esquemas de função de interpolação: CDS, UDS, QUICK estes esquemas foram escolhidos, pois apresentaram os menores erros de truncamentos possíveis. Em síntese, os testes computacionais mostraram as limitações dos esquemas CDS e UDS, bem como também revelaram os ótimos resultados do esquema QUICK. / In this work, the mathematical and computer modeling in the analysis of the interpolation functions in the Finite Volume Method was developed. The equations were integrated in the control of volume in which frontiers the advective and diffuse flows must be evaluated. The interpolation function, which main objective is to connect the nodal points, was used to calculate those flows. Three interpolation function schemes were used: CDS, UDS and QUICK. These schemes were chosen due to the fact that they present less errors and blockings. To sum up, the computer tests showed the limitations of the CDS and UDS schemes and the excellent results of the QUICK scheme.

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