Dedico este trabalho a todas as pessoas que passaram pela minha vida e que fazem parte dela, e com as quais aprendi lições valiosas que seguirão norteando meu caminho para sempre.Epígrafe "Aprender é a única coisa de que a mente nunca se cansa, nunca tem medo e nunca se arrepende." -Leonardo da Vinci v 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 acoplada 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 que nos 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.