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

Esta tese apresenta um modelo matemático para simulação numérica do escoamento com turbilhonamento (swirl) em tochas de plasma térmico de arco não transferido que operam em corrente contínua, assim como os resultados obtidos com as simulações para estudo de casos. O modelo magneto-hidrodinâmico (modelo MHD) bidimensional permitiu simular a interação entre o escoamento e o arco elétrico usando uma configuração axissimétrica, que abrange as seguintes regiões: entrada do gás; interior da tocha; jato de plasma livre no ambiente. O modelo foi implementado num código numérico baseado no Método dos Volumes Finitos para a solução numérica das equações governantes. Para os estudos foram simulados casos com diferentes condições operacionais (vazão de gás; intensidade de corrente elétrica; gases plasmogênicos: ar e argônio; intensidade de turbilhonamento). A fim de verificar a qualidade do modelo, alguns resultados foram comparados com a literatura e apresentaram boa concordância: a maior diferença obtida entre valores de temperatura experimentais e valores calculados foi -10%, e a média das diferenças obtidas nas comparações foi de aproximadamente ±3,2%. Os perfis de temperatura e de velocidade obtidos para a região do arco e para o jato de plasma resultante permitiram o estudo do comportamento do escoamento na tocha de plasma em diferentes condições. Conclui-se que o modelo desenvolvido é apto à realização de investigações numéricas do escoamento em tochas de plasma e dos efeitos do turbilhonamento na interação arco/escoamento. / This thesis presents a mathematical model for numerical simulation of swirling flow in DC non-transferred arc thermal plasma torches, as well as the results obtained from simulations to case studies. The two-dimensional magnetohydrodynamic model (MHD model) allowed simulate the interaction between the flow and the electric arc using an axisymmetric configuration, covering the following areas: gas inlet; inside the torch; free jet of plasma in the environment. The model was implemented in a computer code based on the Finite Volume Method (FVM) to enable the numerical solution of the governing equations. For the study, cases were simulated with different operating conditions (gas flow rate; electric current intensity; plasmogenic gases: air and argon; swirl intensity). In order to verify the quality of the model, some results were compared with the literature and showed good agreement: the biggest difference between experimental temperature values and calculated values was 10%, and the average of the differences obtained in the comparisons was approximately ±3.2%. The resulting profiles of temperature and velocity obtained for the region of the arc and the plasma jet allowed the study of the flow behavior in the plasma torch in different conditions. It is concluded that the model developed is able to carry out numerical investigations of the flow in plasma torches and the effects of swirl in the interaction arc/flow.

Orientador: Sônia Maria Gomes, Jorge Stolfi / Tese (doutorado) - Universidade Estadual de Campinas, Instituto de Matemática, Estatística e Computação Científica / Made available in DSpace on 2018-08-19T12:39:30Z (GMT). No. of bitstreams: 1 Castro_DouglasAzevedo_D.pdf: 8872633 bytes, checksum: a17b2761789c6a831631ac143fdf5ca7 (MD5) Previous issue date: 2011 / Resumo: O objetivo desta tese é desenvolver algoritmos baseados em malhas e bases funcionais inovadoras usando técnicas de multiescala para aproximação de funções e resolução de problemas de equações diferenciais. Para certas classes de problemas, é possível incrementar a eficiência dos algoritmos de multiescala usando bases adaptativas, associadas a malhas construídas de forma a se ajustarem com o fenômeno a ser modelado. Nesta abordagem, em cada nível da hierarquia, os detalhes entre a aproximação desse nível e a aproximação definida no próximo nível menos refinado pode ser usada como indicador de regiões que necessitam de mais ou menos refinamento. Desta forma, em regiões onde a solução é suave, basta utilizar os elementos dos níveis menos refinados da hierarquia, enquanto que o maior refinamento é feito apenas onde a solução tiver variações bruscas. Consideramos dois tipos de formulações para representações multiescala, dependendo das bases adotadas: splines diádicos e wavelets. A primeira abordagem considera espaços aproximantes por funções splines sobre uma hierarquia de malhas cuja resolução depende do nível. A outra abordagem considera ferramentas da analise wavelet para representações em multirresolução de médias celulares. O enfoque está no desenvolvimento de algoritmos baseados em dados amostrais d-dimensionais em malhas diádicas que são armazenados em uma estrutura de árvore binária. A adaptatividade ocorre quando o refinamento é interrompido em algumas regiões do domínio, onde os detalhes entre dois níveis consecutivos são suficientemente pequenos. Um importante aspecto deste tipo de representação é que a mesma estrutura de dados é usada em qualquer dimensão, além de facilitar o acesso aos dados nela armezenados. Utilizamos as técnicas desenvolvidas na construção de um método adaptativo de volumes finitos em malhas diádicas para a solução de problemas diferenciais. Analisamos o desempenho do método adaptativo em termos da compressão de memória e tempo de CPU em comparação com os resultados do esquema de referência em malha uniforme no nível mais refinado. Neste sentido, comprovamos a eficiência do método adaptativo, que foi avaliada levando-se em consideração os efeitos da escolha de diferentes tipos de fluxo numérico e dos parâmetros de truncamento / Abstract: The goal of this thesis is to develop algorithms based on innovative meshes and functional bases using multiscale techniques for function approximation and solution of differential equation problems. For certain classes of problems, one can increase the efficiency of multiscale algorithms using hierarchical adaptive bases, associated to meshes whose resolution varies according to the local features of the phenomenon to be modeled. In this approach, at each level of the hierarchy the details-differences between the approximation for that level and that of the next coarser level-can be used as indicators of regions that need more or less refinement. In this way, in regions where the solution is smooth, it suffices to use elements of the less refined levels of the hierarchy, while the maximum refinement is used only where the solution has sharp variations. We consider two classes of formulations for multiscale representations, depending on the bases used: dyadic splines and wavelets. The first approach uses approximation spaces consisting of spline functions defined over a mesh hierarchy whose resolution depends on the level. The other approach uses tools from wavelet analysis for multiresolu-tion representations of cell averages. The focus is on the development of algorithms based on sampled d-dimensional data on dyadic meshes which are stored in a binary tree structures. The adaptivity happens when the refinement is interrupted in certain regions of the domain, where the details between two consecutive levels are sufficiently small. This representation greatly simplifies the access to the data and it can be used in any dimension. We use these techniques to build an adaptive finite volume method on dyadic grids for the solution of differential problems. We analyze the performance of the method in terms of memory compression and CPU time, comparing it with the reference scheme (which uses a uniform mesh at the maximum refinement level). In these tests, we confirmed the efficiency of the adaptive method for various numeric flow formulas and various choices of the thresholding parameters / Doutorado / Matematica Aplicada / Doutor em Matemática Aplicada

Este trabalho tem como principal objetivo o desenvolvimento de um método numérico para simulação das grandes escalas de escoamentos turbulentos tridimensionais utilizando uma modelagem de turbulência baseada em filtragem temporal (denominada TLES - Temporal Large Eddy Simulation). O método desenvolvido combina discretizações temporais com ordem de mínima precisão 2 (Adams-Bashforth, QUICK, Runge-Kutta), um método de projeção de ordem 2, com discretizações espaciais também de ordem 2 obtidas pelo método de volumes finitos. Esta metodologia foi empregada na simulação de problemas teste turbulentos como o canal e a cavidade impulsionada, sendo este último resultado simulado pela primeira vez com modelagem TLES. Os resultados mostram uma excelente concordância quando comparado com resultados de simulações diretas (DNS) e dados experimentais, superando resultados clássicos obtidos com formulação LES com filtragem espacial. / The main objective of this work is to develop a numerical method for large eddy simulation of tridimensional turbulent flows using a model based on temporal filtering (TLES - Temporal Large Eddy Simulation). The developed method combines at least 2nd order temporal discretizations (Adams-Bashforth, QUICK, Runge-Kutta), a 2nd order projection method, and 2nd order spatial discretizations obtained by the finite volume method. This methodology was employed to the simulation of turbulent benchmark problems such as channel and lid-driven cavity flows. The latter is simulated for the first time using a TLES turbulence modelling. Results show excellent agreement when compared to Direct Numerical Simulations (DNS) and experimental data, with better results than classical results produced by standard LES formulation with spatial filtering.

Submitted by Johnny Rodrigues (johnnyrodrigues@ufcg.edu.br) on 2018-03-23T01:44:32Z No. of bitstreams: 1 GUSTAVO GOMES SAMPAIO CURSINO - TESE PPGEQ 2016.pdf: 5939681 bytes, checksum: 518baf4150ea2fb7085706252276b9fa (MD5) / Made available in DSpace on 2018-03-23T01:44:32Z (GMT). No. of bitstreams: 1 GUSTAVO GOMES SAMPAIO CURSINO - TESE PPGEQ 2016.pdf: 5939681 bytes, checksum: 518baf4150ea2fb7085706252276b9fa (MD5) Previous issue date: 2016-04-18 / Este trabalho teve o propósito de determinar o comportamento dos gases na seção de radiação de um combustor de ar que pertence a uma planta industrial. O corpo metálico do equipamento rompeu em seu primeiro ano de operação, devido a um problema conceitual em sua geometria. A fluidodinâmica computacional (CFD), por meio do método dos volumes finitos, foi utilizada para desenvolver um modelo tridimensional que pudesse reproduzir o perfil de temperatura e o comportamento do fluxo do ar de combustão no equipamento. Na simulação, através do uso do software ANSYS CFX, foram utilizados: (i) o modelo de turbulência Reynolds Stress Model (RSM); (ii) as malhas hexaédrica, tetraédrica e prismática; (iii) o modelo de radiação P-1; e (iv) o modelo de combustão Eddy Dissipation Concept (EDC). Como resultado, foram apresentadas quatro possíveis mudanças na geometria do combustor de ar que, caso adotadas, eliminariam os riscos de novas falhas e garantiriam a continuidade operacional da unidade de processo. / This paper has the objective to describe the behavior of the flow and temperature of the flue gas in the radiation section of the vessel used to preheat air in a combustor. The equipment failed in its first operational year, due to a conceptual problem in its geometry. The CFD code based on finite volume method was applied to simulate the physical model of combustor using the ANSYS CFX software, reproducing the main features of the preheater. The simulation had considered: (i) Reynolds Stress Model (RSM) as turbulence model, (ii) The meshes applied were the hexahedral, tetrahedral and prismatic, (iii) P-1 was used as the radiation model and (iv) Eddy Dissipation Concept (EDC) as combustion model. Through the simulation was possible to propose four different kind of combustor geometry modification, that the application of anyone of them would eliminate the risk of new failures, ensuring the unit production availability.

Neste trabalho apresentamos um esquema numérico para a simulação computacional de escoamentos bifásicos, água-óleo, em reservatórios de petróleo. O modelo matemático consiste em um sistema de equações diferenciais parciais não-linear nas incógnitas velocidade, pressão e saturação. Uma quebra de operadores a dois níveis possibilita uma maior eficiência ao método permitindo que a velocidade, fornecida pelo problema de velocidade-pressão, seja atualizada somente para determinados intervalos de tempo associados ao problema de transporte advectivo-difusivo em termos da saturação. O método dos volumes finitos é empregado na resolução numérica do problema de velocidade-pressão e do transporte de massa por advecção e difusão. Na solução do problema de transporte de massa utilizamos limitadores de fluxo na aproximação dos termos advectivos e diferenças centradas para os termos difusivos. O nosso simulador foi validado a partir de confrontações dos seus resultados com as soluções teóricas conhecidas para os problemas unidimensionais, equações de Burgers e de Buckley-Leverett, e com outros resultados numéricos em se tratando do escoamento bifásico água-óleo bidimensional em meios porosos heterogêneos. / A new numerical method is proposed for the solution of two-phase flow problem in petroleum reservoirs. The two-phase (water and oil) flow problem is governed by a pressure-velocity equation coupled to a saturation equation. For computational eficiency an operator spliting technique is used; distinct time steps can be used for the computation of transport and pressure-velocity problems. The finite volume method is used in the numerical solution of the velocity-pressure and mass transport problems. A flux limiter is used for the numerical discretization of the advective terms while centered schemes are employed for the diffusion terms in the mass transport problem. In the validation of our numerical method we compared numerical and theoretical solutions for one dimensional problems, Burgers and Buckley-Leverett equations, and compared our numerical results to others, in the case of oil-water flows in two dimensions for an heterogeneous porous media.