Made available in DSpace on 2016-12-12T20:25:12Z (GMT). No. of bitstreams: 1 Maico Jeremias da Silva.pdf: 5502697 bytes, checksum: 292822b79334828b19d2d99087aae96b (MD5) Previous issue date: 2015-07-28 / Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / This work discusses the behaviour of the convection heat transfer coefficient in plane fins with variation of air flow velocity and fin spacing. The differential equation that governs heat transfer in fins is discretized by the finite volume method, which enables computation of the fin thermal characteristics, such as temperature distribution, heat transfer and fin efficiency. The determination of the convection coefficient, important parameter for thermal analysis, is performed by applying a heuristic optimization method, known as Particle Swarm Optimization, which combines data measured from experimental analysis conducted in wind tunnel with the aforementioned heat transfer numerical approximation of fins. / Este trabalho analisa os efeitos da velocidade do escoamento de ar e do espaçamento entre aletas no coeficiente de convecção forçada sobre as superfícies de aletas planas. As equações diferenciais que definem a transferência de calor em aletas são discretizadas pela técnica de volumes finitos possibilitando a determinação das características térmicas da aleta, tais como perfil de temperatura, calor transferido e eficiência térmica. A determinação do coeficiente de convecção, parâmetro fundamental para análises térmicas, é realizada mediante a aplicação de um método de otimização heurístico, conhecido como Método do Enxame de Partículas, que combina os dados obtidos da análise experimental realizada em túnel de vento com a aproximação numérica da troca de calor na aleta descrita acima

O presente trabalho trata de construir um modelo computacional utilizando o método dos volumes finitos para malhas não-estruturadas, a fim de se calcular a carga hidráulica num meio poroso, considerando este meio não homogêneo e anisotrópico. A anisotropia é uma característica de muitos materiais encontrados na natureza e depende da propriedade estudada no meio. Primeiramente apresenta-se a dedução da equação do transporte advectivo dispersivo e a formulação matemática para a equação de Laplace, esta última utilizada para o cálculo da carga hidráulica. Em seguida, apresenta-se o algoritmo de solução de um programa computacional em linguagem C++ que permite calcular a velocidade do fluxo em cada face de um volume de controle. Finalmente são feitos vários testes para validação do código computacional utilizado, o que levou a crer que o método utilizado é eficaz para os tipos de malhas testados, apresentando algumas diferenças quanto ao erro da solução. / The present work build a computational model using the finite volumes method for unstructured meshes, with the purpose of calculating the hydraulic load in a porous medium, considering it material non - homogeneous and anisotropic. The Anisotropy is a characteristic of many materials found in the nature and it depends on the property studied in this material. First, we present the deduction of the equation of advective-dispersive transport and the mathematical formulation for the Laplaces equation, this last one used for the calculation of the hydraulic load. Soon afterwards, we present the solution algorithm of a computational program in the C++ language that allows to calculate the speed of the flow in each face of the control volume. Finally several tests for validation of the code are made, which makes it that the plausible to assume method is effective for the types of meshes tested, presenting some differences for the wrong solution.

Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro / A Amazônia exibe uma variedade de cenários que se complementam. Parte desse ecossistema sofre anualmente severas alterações em seu ciclo hidrológico, fazendo com que vastos trechos de floresta sejam inundados. Esse fenômeno, entretanto, é extremamente importante para a manutenção de ciclos naturais. Neste contexto, compreender a dinâmica das áreas alagáveis amazônicas é importante para antecipar o efeito de ações não sustentáveis. Sob esta motivação, este trabalho estuda um modelo de escoamento em áreas alagáveis amazônicas, baseado nas equações de Navier-Stokes, além de ferramentas que possam ser aplicadas ao modelo, favorecendo uma nova abordagem do problema. Para a discretização das equações é utilizado o Método dos Volumes Finitos, sendo o Método do Gradiente Conjugado a técnica escolhida para resolver os sistemas lineares associados. Como técnica de resolução numérica das equações, empregou-se o Método Marker and Cell, procedimento explícito para solução das equações de Navier-Stokes. Por fim, as técnicas são aplicadas a simulações preliminares utilizando a estrutura de dados Autonomous Leaves Graph, que tem recursos adaptativos para manipulação da malha que representa o domínio do problema / Amazon exhibits a variety of scenarios that complement each other. Yearly, part of this ecosystem suffers severe changes in its hydrological cycle, causing flood through vast stretches of the forest. This phenomenon, however, is extremely important for the maintenance of natural cycles. In this context, understanding the Amazonian floodplains dynamics is important to anticipate the effect of unsustainable actions. Under this motivation, this work presents a study of an Amazonian floodplains flow model, based on Navier-Stokes equations, besides other tools that can be applied to the model, thus leading to a new approach to the problem. To discretize the model equations the Finite Volume Method was employed. To solve the associated linear systems the Conjugate Gradient Method technique was chosen. For numerical solution of the equations, the Marker and Cell Method, a explicit solution procedure for Navier-Stokes equations, was employed. All these techniques are applied to preliminary simulations using the Autonomous Leaves Graph algorithm which has adaptive features to manipulate the domain grid problem

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards. / Geotechnical Engineering is the area of Civil Engineering that studies the interaction between constructions carried out by man or natural phenomena with geological environment, which most of times is partially saturated soil. In this sense, work developing as stabilization, dam containing, retaining walls, foundations and highways are conditioned to a right prediction of water flow into the soil. However, considering the water flow, the studied region areas are commonly on the order of square kilometers, mathematical models solutions require computational meshes of large proportions, causing serious limitations linked to computational memory requirements and processing time. In order to overcome these limitations, efficient numerical methods must be used in the solution of the considered problem. Hence iterative methods for solving nonlinear and large sparse linear systems must be used in this type of application. In short, this study approached a solution to the Richard partial differential equation by the two dimensions finite volume method, bringing Picard and Newton method with greater efficiency. Linear system resolution iterative techniques based on Krylov space with pre-conditioners matrix were used. Portable Extensible Toolkit for Scientific Computation (PETSc) numerical library was a tool used during the task. The results indicate when a Richards equation is solved considering thr PICARD-KRYLOV method, no matter the soil evaluation model, the best combination for solving linear systems is the stabilized double gradient method and the SOR preconditioning. On the other hand, when the van Genuchten equations are used the gradients methods with the SOR preconditioning must be chosen. Adopting the NEWTON-KRYLOV method, the stabilized double gradient method is more efficient in soling Newton linear system, in relation to the preconditioning it must be giving preference to the Jacob block. Finally, there are strong indications that the PICARDKRYLOV method can be more effective than the NEWTON-KRYLOV one, when used for solving Richards partial differential equation.

Made available in DSpace on 2016-06-02T19:56:42Z (GMT). No. of bitstreams: 1 3506.pdf: 2271432 bytes, checksum: fc1d05123fb3be32395a9e84b2da5ad8 (MD5) Previous issue date: 2011-02-25 / Universidade Federal de Sao Carlos / This work investigated numerical methods in the solution of mathematical models of fixed-bed reactors. For the reactors modeling and simulation, two numerical methods were used: sequencing method (SM) and finite volume method (FVM). There were also proposed two mathematical models: the pseudo-homogeneous model and the dimensionless one, which is based on the Peclet (Pe) and Biot (Bi) numbers. A horizontal-flow anaerobic immobilized sludge (HAIS) reactor developed in bench scale and after a scale-up, reducing the COD in the wastewater treatment was simulated by sequencing method, varying the numbers of mesh; a tubular fixed-bed reactor with biomass immobilized for the startup period of lactic acid fermentation, also simulated by sequencing method and compared with experimental data; and was also evaluated the precision of sequencing and finite volume methods over the reactor s profile, varying the Peclet e Biot numbers. The models development was based on studies about hydrodynamics and biochemistry kinetics. Both methods described satisfactorily the behavior of the reactors in the performed simulations, but in high values of Peclet, the finite volume method generated inadequacies such as oscillatory responses and over the limit. This paper is an elucidation to sequencing method, which besides its huge range and simplicity, still is not so studied neither known, because it s a recent method. / Este trabalho investigou metodos numericos na solucao de modelos matematicos para reatores de leito fixo. Para a modelagem e simulacao dos reatores, foram utilizados dois metodos numericos: metodo da sequencia (SM) e metodo dos volumes finitos (FVM). Foram propostos dois modelos matematicos: o pseudo-homogeneo e o adimensional, sendo este ultimo baseado nos numeros de Peclet (Pe) e Biot (Bi). Um reator anaerobio horizontal de leito fixo (RAHLF) desenvolvido inicialmente em escala de bancada e posterior aumento de escala na reducao da DQO de aguas residuarias foi simulado pelo metodo da sequencia, variando-se o numero de malhas; um reator tubular de leito fixo com biomassa anaerobia imobilizada no periodo da partida da fermentacao de acido latico, tambem simulado pelo metodo da sequencia e comparado com dados experimentais; e avaliou-se a precisao dos metodos da sequencia e dos volumes finitos sobre o perfil da concentracao de um reator, variando-se os valores de Peclet e Biot. O desenvolvimento dos modelos foi baseado em estudos sobre caracteristicas hidrodinamicas do sistema e de cinetica bioquimica. Ambos os metodos descreveram satisfatoriamente o comportamento dos reatores nas simulacoes realizadas, porem em valores elevados de Peclet, o metodo dos volumes finitos gerou inadequacoes como respostas oscilatorias e superiores ao limite. Este trabalho foi uma elucidacao ao metodo da sequencia, que apesar da sua grande abrangencia e simplicidade, por ser um metodo recente, ainda e pouco estudado e conhecido.

O presente trabalho trata de construir um modelo computacional utilizando o método dos volumes finitos para malhas não-estruturadas, a fim de se calcular a carga hidráulica num meio poroso, considerando este meio não homogêneo e anisotrópico. A anisotropia é uma característica de muitos materiais encontrados na natureza e depende da propriedade estudada no meio. Primeiramente apresenta-se a dedução da equação do transporte advectivo dispersivo e a formulação matemática para a equação de Laplace, esta última utilizada para o cálculo da carga hidráulica. Em seguida, apresenta-se o algoritmo de solução de um programa computacional em linguagem C++ que permite calcular a velocidade do fluxo em cada face de um volume de controle. Finalmente são feitos vários testes para validação do código computacional utilizado, o que levou a crer que o método utilizado é eficaz para os tipos de malhas testados, apresentando algumas diferenças quanto ao erro da solução. / The present work build a computational model using the finite volumes method for unstructured meshes, with the purpose of calculating the hydraulic load in a porous medium, considering it material non - homogeneous and anisotropic. The Anisotropy is a characteristic of many materials found in the nature and it depends on the property studied in this material. First, we present the deduction of the equation of advective-dispersive transport and the mathematical formulation for the Laplaces equation, this last one used for the calculation of the hydraulic load. Soon afterwards, we present the solution algorithm of a computational program in the C++ language that allows to calculate the speed of the flow in each face of the control volume. Finally several tests for validation of the code are made, which makes it that the plausible to assume method is effective for the types of meshes tested, presenting some differences for the wrong solution.

Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro / A Amazônia exibe uma variedade de cenários que se complementam. Parte desse ecossistema sofre anualmente severas alterações em seu ciclo hidrológico, fazendo com que vastos trechos de floresta sejam inundados. Esse fenômeno, entretanto, é extremamente importante para a manutenção de ciclos naturais. Neste contexto, compreender a dinâmica das áreas alagáveis amazônicas é importante para antecipar o efeito de ações não sustentáveis. Sob esta motivação, este trabalho estuda um modelo de escoamento em áreas alagáveis amazônicas, baseado nas equações de Navier-Stokes, além de ferramentas que possam ser aplicadas ao modelo, favorecendo uma nova abordagem do problema. Para a discretização das equações é utilizado o Método dos Volumes Finitos, sendo o Método do Gradiente Conjugado a técnica escolhida para resolver os sistemas lineares associados. Como técnica de resolução numérica das equações, empregou-se o Método Marker and Cell, procedimento explícito para solução das equações de Navier-Stokes. Por fim, as técnicas são aplicadas a simulações preliminares utilizando a estrutura de dados Autonomous Leaves Graph, que tem recursos adaptativos para manipulação da malha que representa o domínio do problema / Amazon exhibits a variety of scenarios that complement each other. Yearly, part of this ecosystem suffers severe changes in its hydrological cycle, causing flood through vast stretches of the forest. This phenomenon, however, is extremely important for the maintenance of natural cycles. In this context, understanding the Amazonian floodplains dynamics is important to anticipate the effect of unsustainable actions. Under this motivation, this work presents a study of an Amazonian floodplains flow model, based on Navier-Stokes equations, besides other tools that can be applied to the model, thus leading to a new approach to the problem. To discretize the model equations the Finite Volume Method was employed. To solve the associated linear systems the Conjugate Gradient Method technique was chosen. For numerical solution of the equations, the Marker and Cell Method, a explicit solution procedure for Navier-Stokes equations, was employed. All these techniques are applied to preliminary simulations using the Autonomous Leaves Graph algorithm which has adaptive features to manipulate the domain grid problem

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / A engenharia geotécnica é uma das grandes áreas da engenharia civil que estuda a interação entre as construções realizadas pelo homem ou de fenômenos naturais com o ambiente geológico, que na grande maioria das vezes trata-se de solos parcialmente saturados. Neste sentido, o desempenho de obras como estabilização, contenção de barragens, muros de contenção, fundações e estradas estão condicionados a uma correta predição do fluxo de água no interior dos solos. Porém, como a área das regiões a serem estudas com relação à predição do fluxo de água são comumente da ordem de quilômetros quadrados, as soluções dos modelos matemáticos exigem malhas computacionais de grandes proporções, ocasionando sérias limitações associadas aos requisitos de memória computacional e tempo de processamento. A fim de contornar estas limitações, métodos numéricos eficientes devem ser empregados na solução do problema em análise. Portanto, métodos iterativos para solução de sistemas não lineares e lineares esparsos de grande porte devem ser utilizados neste tipo de aplicação. Em suma, visto a relevância do tema, esta pesquisa aproximou uma solução para a equação diferencial parcial de Richards pelo método dos volumes finitos em duas dimensões, empregando o método de Picard e Newton com maior eficiência computacional. Para tanto, foram utilizadas técnicas iterativas de resolução de sistemas lineares baseados no espaço de Krylov com matrizes pré-condicionadoras com a biblioteca numérica Portable, Extensible Toolkit for Scientific Computation (PETSc). Os resultados indicam que quando se resolve a equação de Richards considerando-se o método de PICARD-KRYLOV, não importando o modelo de avaliação do solo, a melhor combinação para resolução dos sistemas lineares é o método dos gradientes biconjugados estabilizado mais o pré-condicionador SOR. Por outro lado, quando se utiliza as equações de van Genuchten deve ser optar pela combinação do método dos gradientes conjugados em conjunto com pré-condicionador SOR. Quando se adota o método de NEWTON-KRYLOV, o método gradientes biconjugados estabilizado é o mais eficiente na resolução do sistema linear do passo de Newton, com relação ao pré-condicionador deve-se dar preferência ao bloco Jacobi. Por fim, há evidências que apontam que o método PICARD-KRYLOV pode ser mais vantajoso que o método de NEWTON-KRYLOV, quando empregados na resolução da equação diferencial parcial de Richards. / Geotechnical Engineering is the area of Civil Engineering that studies the interaction between constructions carried out by man or natural phenomena with geological environment, which most of times is partially saturated soil. In this sense, work developing as stabilization, dam containing, retaining walls, foundations and highways are conditioned to a right prediction of water flow into the soil. However, considering the water flow, the studied region areas are commonly on the order of square kilometers, mathematical models solutions require computational meshes of large proportions, causing serious limitations linked to computational memory requirements and processing time. In order to overcome these limitations, efficient numerical methods must be used in the solution of the considered problem. Hence iterative methods for solving nonlinear and large sparse linear systems must be used in this type of application. In short, this study approached a solution to the Richard partial differential equation by the two dimensions finite volume method, bringing Picard and Newton method with greater efficiency. Linear system resolution iterative techniques based on Krylov space with pre-conditioners matrix were used. Portable Extensible Toolkit for Scientific Computation (PETSc) numerical library was a tool used during the task. The results indicate when a Richards equation is solved considering thr PICARD-KRYLOV method, no matter the soil evaluation model, the best combination for solving linear systems is the stabilized double gradient method and the SOR preconditioning. On the other hand, when the van Genuchten equations are used the gradients methods with the SOR preconditioning must be chosen. Adopting the NEWTON-KRYLOV method, the stabilized double gradient method is more efficient in soling Newton linear system, in relation to the preconditioning it must be giving preference to the Jacob block. Finally, there are strong indications that the PICARDKRYLOV method can be more effective than the NEWTON-KRYLOV one, when used for solving Richards partial differential equation.

Orientadores: Jose Roberto Nunhez, Edson Tomaz / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-06T01:51:29Z (GMT). No. of bitstreams: 1 Machado_MarcioBezerra_D.pdf: 3569468 bytes, checksum: 4c69c5d76d7d79aa717808b84c4701d2 (MD5) Previous issue date: 2006 / Resumo: Estudos têm mostrado que a humanidade enfrentará severa falta de água nas próximas décadas. Muitos esforços têm sido direcionados para o desenvolvimento de novas ferramentas computacionais a fim de se garantir uma melhor utilização dos recursos hídricos. Diversos estudos estão sendo realizados utilizando ferramentas de CFD (Computational Fluid Dynamics) para obtenção de novas formas de gerenciamento destes recursos. Neste contexto, é de suma importância o desenvolvimento de novas técnicas para predizer o impacto ambiental causado por emissões industriais em rios de modo que estratégias possam ser planejadas para diminuir os efeitos desta poluição. Este trabalho apresenta um modelo Fluidodinâmico Computacional tridimensional para simular a dispersão de substâncias solúveis em rios. O método dos volumes finitos foi utilizado para aproximar as equações de conservação de momento, de massa e de espécie química. O sistema de coordenadas cartesianas foi escolhido para representar o sistema. Foi utilizado um modelo algébrico de turbulência de ordem zero. O modelo de StreeterPhelps foi usado para predizer a concentração de substâncias orgânicas e de oxigênio dissolvido ao longo do rio. O modelo pode também predizer o impacto causado pela ocorrência de múltiplos pontos de emissão no trecho estudado. O modelo matemático foi desenvolvido em linguagem Fortran. Os resultados mostram que a metodologia proposta é uma boa ferramenta para a avaliação do impacto ambiental causado pela emissão de efluentes em rios. O software é bastante rápido, especialmente quando comparado com outros pacotes de CFD disponíveis comercialmente. Foram feitas comparações entre os resultados numéricos e dados experimentais coletados no rio Atibaia. Os resultados numéricos apresentaram uma boa concordância com os dados coletados experimentalmente / Abstract: A future lack of water in the next decades has been observed by many studies. Much effort has been devoted to find strategies which will help to manage proper1y water resources. Theoretical studies have been used recent1y since the scope of computational fluid dynamics (CFD) has increased, allowing its use in the issue of water quality. In this scenario, it is important to develop new techniques to predict the environmental impact of emissions in rivers so that strategies can be devised to decrease the effects of pollution. This work presents a three-dimensional Computational Fluid Dynamics (CFD) in house model to simulate the dispersion of soluble substances in a river. The finite volume method is used to approximate the momentum, mass and species conservation equations. A Cartesian coordinate system has been chosen to represent the river. Turbulence is taken into account by a zero-order equation model. The Streeter-Phelps model has been used to predict the concentration of organic substances and dissolved oxygen along the river. The model can also predict the impact of multiple effluents discharges. Results show that the proposed methodology is a good tool for the evaluation of the environmental impact caused for pollutants emissions in rivers. The software has been developed from the model and use the Fortran language. It is very fast, especially when compared to available commercial CFD packages. Experimental comparisons for soluble substances dispersion have been made for the Atibaia River. The results show good agreement with experimental data / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química

Orientador: Jose Roberto Nunhez / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-11T08:15:27Z (GMT). No. of bitstreams: 1 Lopes_GabrielaCantarelli_M.pdf: 1710566 bytes, checksum: efab0c6308734536b15c2ef91d4d8337 (MD5) Previous issue date: 2008 / Resumo: Processos que envolvem escoamentos de fluidos no interior de dutos estão presentes em muitas aplicações industriais. Na indústria petroquímica um desses processos que vem se tornando cada vez mais importante é o craqueamento catalítico de frações pesadas do petróleo, já que as converte em frações leves e mais nobres. Por causa da crescente demanda mundial por gasolina e GLP e à sua alta rentabilidade para uma refinaria de petróleo, unidades de craqueamento catalítico em leito fluidizado (FCC) estão presentes em todo o mundo. Testes experimentais costumam ser usados no estudo dos fenômenos envolvidos nesse processo. Porém esse tipo de análise possui alto custo, que pode ser reduzido com o uso de simulações computacionais em seu estudo inicial. Assim, o objetivo deste trabalho foi desenvolver um modelo tridimensional e elíptico, em linguagem Fortran, capaz de fornecer dados para a análise preliminar de escoamentos no interior de reatores de FCC. Na modelagem desses problemas são usadas equações diferenciais parciais, e essas não possuem solução analítica conhecida, sendo necessário o emprego de métodos numéricos para esse fim. Neste trabalho foi usado o Método dos Volumes Finitos, que tem a função de substituir as equações diferenciais parciais por equações algébricas aplicadas a pequenos volumes de controle finitos pertencentes ao domínio. Uma das maiores dificuldades encontradas no tratamento numérico de escoamentos incompressíveis é a determinação de um campo de pressão que satisfaça a Equação da Continuidade. Esse problema foi resolvido fazendo-se uso da abordagem acoplp.da de solução. Para análise do modelo foram obtidos perfis numéricos de velocidade e pressão para fluidos escoando em regimes laminar e turbulento, que foram validados usando-se os dados obtidos com a solução analítica das equações, por correlações (semi-) empíricas ou por dados experimentais, conforme cada um dos casos. Notou-se que o modelo representa muito bem casos laminares, e gue 11.oS casos turbulentos foi necessário um maior refino da malha próximo a parede do tubo. Também foram feitas simulações para que se pudessem observar as características tridimensionais, elípticas e transientes da modelagem do escoamento. De maneira geral o modelo se mostrou bastante rápido, convergindo em poucas iterações. Palavras-Chave: Fluidodinâmica computacional; dutos cilíndricos; modelo tridimensional e elíptico; método dos Volumes Finitos; solução acoplada; turbulência / Abstract: Processes involving fluid flow in tubes are present in many industrial applications. In petrochemical industry one of these processes that are becoming more and more important is the fluid catalytic cracking of heavy petroleum fractions. This fact is due to the process capacity to convert heavy fractions in light and valuable ones. Because of the increasingly worldwide demand for gasoline and LPG and its high yield for a petroleum refinery, fluid catalytic cracking (FCC) units are present in the whole world. Experimental tests are used in the study of the phenomena involved in this processo However this kind of analysis has high cost which can be reduced by using computational simulations in its initial study. Thus, the aim of this work was the development of a three-dimensional and elliptical mo dei in Fortran language in order to provide data for fluid flow preliminary analysis in FCC reactors. Partial differential equations were used in the modeling of these problems. These equations do not have known analytical solution, being necessary therefore the use of numerical methods. In this work the Finite Volume Method were applied with this purpose. This method has as a role to substitute the partial differential equations of the mo dei for algebric equations applied to small finite control volumes of the domain. One of the biggest difficulties found in the numerical treatment of incompressible fluid flows is the determination of apressure field that satisfies the Continuity Equation. This problem was solved using the coupled solution approach. For model analysis, numerical velocity and pressure proJ:iles for laminar and turbulent flows were obtained, that had been validated using the data obtained through the analytical solution of the equations, by empirical correlations or by experimental data, according to each one of the cases. The model represented well laminar cases, and in the turbulent ones the mesh had to be more refined near the tube wall. Other simulations were performed, in aQalyzing the three-dimensional, elliptical and transient model characteristics. In general, the mo dei was very fast, converging in a few interations. Keywords: Computational fluid dyn~mic; cylindrical ducts; three-dimensional and elliptical model; Finite Volume Method; coupled solution; turbulence / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química