Orientador: Osvair Vidal Trevisan / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica e Instituto de Geociencias / Made available in DSpace on 2018-08-12T18:20:10Z (GMT). No. of bitstreams: 1 Marques_JoaoBoscoDias_M.pdf: 2275485 bytes, checksum: 7008d88db0703cad4c59a10ca5e7b368 (MD5) Previous issue date: 2008 / Resumo: Este trabalho disserta sobre a análise de teste transitório de produção e de crescimento de pressão em poço vertical sob condição de escoamento bifásico de óleo e gás. Os resultados obtidos são o fator de película e pontos das curvas das permeabilidades efetivas ao óleo e ao gás em função da saturação de óleo. A permeabilidade efetiva em função da saturação do óleo é obtida da combinação da permeabilidade efetiva em função da pressão de teste de poço com a solução de uma equação diferencial ordinária que relaciona a saturação do óleo com a pressão. Um dos métodos, denominado método p2, é desenvolvido a partir da representação da curva da função de pseudopressão, kro/(µoBo) versus a pressão de teste, por meio de duas linhas retas. O problema físico analisado consiste de um reservatório homogêneo operando a partir da pressão de bolha, água conata imóvel e efeitos de capilaridade e de gravidade desprezíveis. Os dados para validação dos métodos são obtidos com testes de poços realizados em modelos de reservatórios preparados num simulador comercial. O método p2 utiliza basicamente pontos e inclinações de retas de curvas típicas, construídas com dados de teste de poço, para determinar o fator de película e as curvas das permeabilidades efetivas aos fluidos. As curvas obtidas pelo método p2 são comparadas com as obtidas por outro método, chamado de método de determinação in situ / Abstract: This thesis provides an analysis of drawdown and buildup tests from vertical well under solution gas-drive conditions. The results obtained include the skin factor and the points of the effective permeability curves for the oil and gas as a function of oil saturation. The effective permeability as a function of oil saturation is the result of a combination of the effective permeability as a function of wellbore pressure and the solution of an ordinary differential equation which relates oil saturation to the pressure. One of the methods, named square-pressure method, is developed using two straight lines to approximate the curve of the pseudo-pressure function, kro/(µoBo) versus the wellbore pressure. The physical problem consists of a homogeneous reservoir operating from bubble-point pressure with immobile connate water, neglecting the capillarity and gravity effects. The data for the validation of the methods was obtained using well tests reservoirs models built in a commercial simulator. The p2 method basically utilizes the points and the straight lines slope of typical curves built from well test data to determinate the effective permeability curves and the skin factor. Then, the results of these curves are compared to those obtained with another method, named in situ determination method / Mestrado / Reservatórios e Gestão / Mestre em Ciências e Engenharia de Petróleo

Orientadores: Milton Mori, Henry França Meier / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-02T16:32:40Z (GMT). No. of bitstreams: 1 Peres_AlexandredePaula_D.pdf: 4006667 bytes, checksum: cea470fc549f4246fc4881f41e89e239 (MD5) Previous issue date: 2002 / Resumo: Este trabalho se propôs a estudar o escoamento do gás em ciclones em duas frentes: por intermédio do mapeamento do campo de velocidades em um ciclone instalado em uma planta semi-piloto de uma unidade da PETROBRAS-SIX e por meio do desenvolvimento de um modelo tluidodinâmico computacional (CFD) tridimensional, amparado pelos resultados experimentais obtidos nesse trabalho e de resultados disponíveis na literatura. Os ciclones têm sido empregados na indústria de processos como separadores de partículas desde o início do século passado, sendo que inúmeros trabalhos de investigação científica visando a obtenção de técnicas que predissessem a queda de pressão e a curva de eficiência de separação são encontrados na literatura. Contudo, somente modernas técnicas, como as da fluido dinâmica computacional (CFD), têm se mostrado capazes de predizer a contento os complexos fenômenos que ocorrem no interior ciclônico. Desta forma, utilizou-se o código de CFD comercial CFX - versões 4.3 e 4.4, da ABA Technology, o qual utiliza o método dos volumes finitos com variáveis co-localizadas resolvendo o sistema em coordenadas generalizadas para a obtenção das soluções das equações de transporte. As malhas numéricas foram construídas por intermédio da técnica dos multi-blocos, utilizando o pré processador "BUILD" presente no código. Os principais resultados mostraram que a metodologia empregada se mostrou adequada e apta a reproduzir o comportamento tluidodinâmico do escoamento do gás (monofásico) e do gás com a presença de particulados (bifásico, em uma abordagem Euleriana-Euleriana), em ciclones tanto em escala de laboratório quanto em escala industrial / Abstract: The aim of this work was to analyze the gas flow in cyclones experimentally and numerically. The experimental study was conducted in an acrylic cyclone settled in a pilot unity belonging to SIX/PETROBRAS in São Mateus do Sul, PR. Numerical solutions were obtained by the development of a 3-D computational fluid dynamics (CFD) model, sustained by the experimental resuIts obtained in this work and results from the literature. Cyclones have been used as solid particle separators in large-scale chemical processes since the beginning of the past century. Several studies of scientific investigation in order to obtain prediction techniques of pressure drop and collection efficiency are found in the literature. However, the design of new cyclones and the analysis of the actual equipment can be best achieved using computational fluid dynamics (CFD) techniques in order to obtain higher collection efficiency and a lower pressure drop. The 3-D simulations of the turbulent gas flow in the cyclone was carried out using the computational fluid-dynamics code CFX, versions 4.3 and 4.4, by AEA Technology. The numerical solutions were obtained with a finite-volume method and body fitted grid generation, using the pre-processor BUILD to generate the grids with multi-block techniques. The methodology showed great capability to represent the swirling flow in small cyclones as in industrial cyclones / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química

Orientador : Marcio Luiz de Souza-Santos / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-02T20:13:42Z (GMT). No. of bitstreams: 1 Costa_AlexandreMarconideSouzada_D.pdf: 1045987 bytes, checksum: 017ef3657be12ae9dedd4df55fd326b1 (MD5) Previous issue date: 2002 / Doutorado

Orientador: Leonardo Goldstein Jr / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-03T18:09:40Z (GMT). No. of bitstreams: 1 Lopes_MarcosDonizete_M.pdf: 5077011 bytes, checksum: c1b9066b877bcfe1af8be5b8a935e84d (MD5) Previous issue date: 2003 / Resumo: O emprego de modelos constitui uma solução prática para estudar a influência dos fatores operacionais no desempenho dos ciclones de porte industrial. A construção dos modelos deve obedecer a fatores de escala apropriados que determinam a existência de semelhança entre o modelo e o protótipo. Modelo e protótipo apresentarão comportamento fluidodinâmico semelhante se forem geometricamente semelhantes e se os valores dos números adimensionais correspondentes forem iguais. Um conjunto de números adimensionais representativo dos fenômenos envolvidos no escoamento gás-sólidos em ciclones deve levar em conta as forças interparticulares, como a de van der Waals e a eletrostática. A dificuldade para obedecer às condições de semelhança levando-as em conta seria, no entanto, muito grande. Assim sendo, desconsiderou-se a eventual participação das forças interparticulares na determinação do grupo de números adimensionais a partir dos quais se gerou as leis de escala. Estas leis simplificadas foram aplicadas para a construção de dois ciclones, modelo e protótipo, em escala, que foram testados para determinação da faixa de aplicação da simplificação introduzida. Adicionalmente abriu-se mão de operar com uma semelhança dinâmica completa entre modelo e protótipo, tendo-se a chamada Semelhança Incompleta, pois, por motivos práticos, operou-se o modelo com ar em condições ambientais, o que equivale a desprezar as forças viscosas em relação às forças de inércia, subtraindo-se o requisito da igualdade do número de Reynolds das relações de escala. O presente trabalho mostrou experimentalmente a existência de concentrações a partir das quais a simplificação proposta perde a validade, implicando que as forças desconsideradas são de magnitude tal que não podem deixar de ser consideradas na geração dos adimensionais. Verificouse que, abaixo destas concentrações limite houve semelhança quanto ao rendimento de coleta, como se constatou pela existência de igualdade dos rendimentos de coleta dos ciclones modelo e protótipo, quando em condições operacionais de velocidades correspondentes e com a mesma concentração de sólidos na entrada dos ciclones. Verificou-se também que a perda de pressão, Dp, adimensionalizada em relação à perda de pressão no ciclone vazio mostrou-se mais adequada que a razão Dp/romo2; para definição da semelhança quanto à perda de pressão, possivelmente porque a velocidade na entrada do ciclone, mo, não seja a velocidade mais adequada a ser utilizada na adimensionalização da perda de pressão / Abstract: The use of models in scale is a practical way of studying how the operational parameters influence the performance of industrial cyclones. The model and the prototype cyclones should be built according to a scaling law - given by a set of dimensionless parameters --, so that they are similar. It there will be fluid-dynamical similarity between the cyclones if they are geometrically similar and the sets of corresponding dimensionless parameters are equal. A group of dimensionless parameters which represent the gas-solid flow phenomenon in cyclones should consider interpartic1e forces, such as the van der Waals and the electrostatic forces. However, as it would be quite difficult to comply with a scaling law determined from taking into consideration these forces, it was decided not to include them and check the range of applicability of the simplified set. Moreover, during the tests, all the dimensionless parameters were made equal, with the exception of the Reynolds number, because, for practical reasons, atmospheric air was used in the tests for both prototype and model cyc1ones, what is equivalent to neglect the viscous forces as compared to the inertial forces. The present work showed that when the model and prototype cyclones were operating at the corresponding velocities and with the same solids concentration at the inlet, similarity was obtained between the collection efficiencies when the solids concentration was less than 1.0 kg/kg, meaning that the proposed methodology was appropriate, that is, that the forces not included were not significant and could be neglected. Concerning the pressure drop in the cyclones, it was observed that the ratio of the pressure drop, Dp, by the pressure drop in the cyclone with no solids, was more appropriate than the ratio Dp/romo2; for expressing similarity, maybe because the entrance velocity, mo, is not the best velocity to be used as a scaling parameter / Mestrado / Termica e Fluidos / Mestre em Engenharia Mecânica

Orientadores : Marco Aurelio Cremasco, Henry França Meier / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-03T22:11:25Z (GMT). No. of bitstreams: 1 Ropelato_Karolline_M.pdf: 3593278 bytes, checksum: 1a6dd755239ec9a10ccf68cc2acb4232 (MD5) Previous issue date: 2004 / Resumo: O tema central desta investigação consiste no desenvolvimento de um modelo matemático multifásico (gás-sólido) tridimensional e transiente, como uma aproximação da dinâmica do escoamento em reatores downer de FCC, ampliando horizontes, com a corroboração ou validação do modelo com dados experimentais da literatura, nos estudos de incremento de performance do processo de FCC. Foram testadas quatro configurações geométricas para a validação do modelo proposto: downer com entrada axial anular; downer com entrada radial cruzada; downer com entrada radial tangencial e o downer com entrada transversal. As principais conclusões obtidas a partir de um estudo de verificação, observando a influência do número de elementos da malha numérica para uma geometria 3-D e um modelo matemático multifásico, transiente e turbulento, bem como do esquema de interpolação, determinando o número adequado de elementos para a resolução das equações diferenciais parciais. Além do estudo de validação, avaliando a influência das condições de alimentação na distribuição de partículas, determinando a melhor configuração geométrica entre as estudadas, para um perfil de distribuição uniforme de partículas no leito do reator, validando o modelo matemático multifásico (gás-sólido) tridimensional e transiente, para a predição da dinâmica do escoamento em reatores downer de FCC, com a corroboração a partir de dados experimentais disponíveis na literatura. Convencionalmente a tecnologia de FCC usa reatores de leito móvel ascendente conhecidos por reator riser. Os chineses, fortemente na década de 90, resolveram analisar o processo de "ponta cabeça", ou seja, em vez de fluxos ascendentes de partículas e gás, estudaram esse mesmo contato considerando-o descendente, o downer. Como resultado, verificaram o aumento da produção de gasolina e diminuição de produtos secundários. Eis o exemplo clássico da modificação conceitual sobre o já existente / Abstract: The present work provides an unsteady, three-dimensional, two-phase, inviscid model of the solid phase, as an approach of the flow dynamics inside downer reactors of FCC processes, and the validation of the proposed model against experimental data available on the literature. Four geometrical configurations were evaluated to validate the proposed model: downer with an annular axial inlet, downer with two cross radial inlets, downer with one tangential inlet and a downer with four transversal inlets. The main conc1usions about the grid size analysis for the unsteady, three-dimensional and multi-phase flow shows the correct grid independence results to solve partial differential equations system. Studies involving the advection scheme were also performed. By comparison between the numerical predictions on this work and experimental data from literature it is possible to conc1ude that the results present good agreement against each other. In fact, the inviscid model of the solid phase, despite of its simplicity modeling, with 3D transient features, could reproduce quite well the main phenomenological characteristics of the turbulent gas-solid flow. In the 90's decades, Chinese's started to study the FCC process gravity concurrent1y, studying down flow (downer reactors) instead of upflow behaviors (riser reactors). The experimental studies show fuel increases production and decreases of secondary products formation. This a c1assical example of a change on conceptual between the downer and riser approaches in industrial facilities / Mestrado / Mestre em Engenharia Química

Orientador: Sergio Nascimento Bordalo / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-24T22:34:31Z (GMT). No. of bitstreams: 1 SevericheMaldonado_JuanCarlos_M.pdf: 5161429 bytes, checksum: 8172a7b4b32287da4014bb30b23addab (MD5) Previous issue date: 1999 / Resumo: o objetivo deste estudo é pesquisar experimentalmente a separação de um fluxo bifásico sob condições de padrão de fluxo estratificado ondulado e anular em uma junta "tê" horizontal de configuração regular e irregular. Os dados da separação foram adquiridos para uma extensa série de condições de fluxo estratificado ondulado e anular; os dados experimentais incluem medições de razões de separação gás-líquido e distribuição da pressão ao longo dos tres ramos. Analisando os dados se observa uma forte dependencia do fenômeno da separação sob o padrão de fluxo existente a montante da junta "tê". Comparando os dados de separação entre uma tê de configuração regular e irregular sob as mesmas condições de fluxo de entrada observa-se que, para condições de fluxo estratificado, ocorre uma substancial redução da fração de líquido que se desvia para o ramo lateral. Para fluxo anular, no entanto, a redução do líquido desviado para o ramo lateral é muito pequena. Grande queda de pressão ocorre no ramo lateral no "tê" de configuração irregular. Isto é mais pronunciado para fluxo anular. O modelo geométrico estudado (Shoham et al. 1987) para predizer o fenômeno da separação em um "tê" regular e em um "tê" irregular comparado com os dados experimentais não apresenta uma boa concordáncia com os resultados experimentais / Abstract: The objective of this study is to investigate experimentally the splitting of two-phase flow under conditions of stratified wavy and annular pattems in horizontal tee-junctions of regular and irregular configurations. Data were acquired for a wide range of stratified wavy and annular flow conditions; experimental data include measurements of the gas and liquid splitting ratios and the pressure distribution along the three arms of the junction. The results show a strong dependence of the splitting phenomenon on the flow pattern existing upstream of the tee. For the stratified wavy flow condition. Comparison of the irregular tee and the regular tee splitting under the same inlet flow condition reveals a substancial reduction in the fraction of the liquid intake that flows to the lateral branch. For the annular flow, however, that reduction is very small. Higher pressure drop occurs in the irregular branch armoThis is . more pronounced for annular flow. The geometrical model estudied (Shoam et al. 1987) for prediction of the splitting phenomenon through regular and irregular tees does not show a reasonable agreement with the data / Mestrado / Mestre em Engenharia de Petróleo

Orientador: Luiz Felipe Mendes de Moura / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-25T04:14:28Z (GMT). No. of bitstreams: 1 Tambellini_Flavio_D.pdf: 10882636 bytes, checksum: 4cd40929cd4ea82a33b7db9c8c968bf0 (MD5) Previous issue date: 1999 / Resumo: o objetivo principal deste trabalho é fazer uma análise da tomografia de capacitância elétrica aplicada ao sistema bifásico, onde os componentes são ar e água, visando obter a reconstrução de imagem deste sistema. A tomografia de capacitância é constituída de um sensor com Seletrodos localizados em tomo de um tubo de acrílico e um transdutor eletrônico baseado em corrente alternada, cuja frequência de excitação pode atingir até 16 MHz. Foi feita a calibração do transdutor eletrônico e esta comparada com as curvas teóricas. Testes experimentaise simulações numéricas foram realizados, sendo os resultados comparados para frações de líquido variando de 0,1 a 0,95. A imagem foi reconstruída usando o algoritmo de retro-projeção linear, juntamente com a utilização da técnica de filtro para melhorar a imagem reconstruída e com uma estimativa da fração de líquido, obtida a partir de uma correlação entre a fração de líquido volumétrica e as capacitâncias experimentais. As informações necessárias para o algoritmo são as medidas experimentaisde capacitância, a fração de líquido estimada e a distribuição de sensibilidade, sendo esta última, obtida a partir de simulações numéricas. Os resultados são apresentados, as limitações descritas e propostas de melhoria para trabalhos futuros apresentadas / Abstract: The main objective ofthis work is the analysis ofthe electrical capacitance tomography for imaging of two-phase system, where the components are water and air. The capacitance tomography is consists by a sensor built with eight electrodes around an insulating pipe. An electronic transducer based in alternated current was developed to support excitation frequency up to 1611Hz. The calibration of the electronic transducer was made and compared with the theoretical curves. Numerical and experimental tests were made and the results compared for liquid fractions ranging from 0,1 to 0,95. The image was reconstructed using the linear back projection algoritm. A filtering technique was applied to the algoritm to improve the image quality. A liquid fraction estimative was obtained from the correlation of volumetric liquid fraction and experimental capacitances. The necessary information to the algoritm are the experimental measurements of capacitance, the estimated liquid fraction and sensitivitydistribuition. The results are shown, the limitationsare described and the improvement ofthis technique to the future works are indicated. / Doutorado / Termica e Fluidos / Doutor em Engenharia Mecânica

Orientador: Niederauer Mastelari / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-08-14T13:06:41Z (GMT). No. of bitstreams: 1 Correa_FernandaCristina_M.pdf: 2774513 bytes, checksum: 26936aae7fde336d3017925056dfbe4e (MD5) Previous issue date: 2009 / Resumo: Escoamentos bifásicos gás-líquido são comumente encontrados na natureza e em diferentes atividades industriais. Neste tipo de escoamento, as fases presentes podem escoar dispostas em diferentes configurações espaciais no interior do duto, denominadas de regimes ou padrões de escoamento, cuja determinação é crítica para o projeto de sistemas de escoamento bifásico. No presente trabalho, a indicação dos padrões de escoamento é realizada a partir da medida de fração de gás utilizando um instrumento microcontrolado composto por um transdutor e um condicionador de sinais. Um sistema de reconhecimento de padrões é composto por diferentes módulos, que consistem em: adquirir o sinal; extrair informações relevantes do sinal; e por fim a classificação propriamente dita. Para a extração das informações relevantes do sinal, aplicou-se funções estatísticas de média e desvio padrão por ser interessante sob o ponto de vista computacional. Na etapa classificatória, a utilização de um conjunto de regras através de técnicas de Sistemas Especialistas é utilizado pelo desempenho deste tipo de classificação e também pouca necessidade computacional. Desta forma, com a junção das funções estatísticas de média e desvio padrão e um conjunto de regras com Sistemas Especialistas, é desenvolvido neste trabalho um sistema de reconhecimento de padrões de escoamento bifásico gás-líquido embarcado em um instrumento de medição de fração de gás microcontrolado. / Abstract: Two-phase flow gas-liquid are commonly found in nature and in various industrial activities. In this type of flow, the phases present can escape arranged in different spatial configurations in the duct, called regimes or flow patterns, the determination is critical for the design of two-phase flow systems. In this study, an indication of the flow pattern is achieved from the measurement of gas fraction microcontrolled using an instrument composed of a transducer and signal conditioner. A system of pattern recognition consists of different modules, which consist of: acquiring the signal, extracting relevant information from the signal, and finally the classification itself. Statistical functions of mean and standard deviation are applied to extract the relevant information from the signal, for computacional purposes. In the classification stage, a set of rules using techniques of expert systems is used by the performance of this type of classification and also little need computing. Thus, with the addition of the statistical functions of mean and standard deviation and a set of rules with expert systems is developed in this work a system of pattern recognition of two-phase flow gas-liquid loaded into an microcontrolled instrument for measuring gas fraction. / Mestrado / Mecanica dos Sólidos e Projeto Mecanico / Mestre em Engenharia Mecânica

Orientadores: Fernando de Almeida França, Edson Yoshihito Nekagawa / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-20T08:46:19Z (GMT). No. of bitstreams: 1 Sousa_JoaoTadeuVidalde_M.pdf: 30145690 bytes, checksum: 24119006d22dc1551983013469607047 (MD5) Previous issue date: 1995 / Resumo: O presente trabalho apresenta os resultados de uma investigação experimental realizada para avaliar a influência de um meio não-newtoniano e da inclinação da tubulação sobre a velocidade de ascensão de bolhas alongadas. o aparato experimental consiste de um tubo de acrílico transparente de 9,44 m, com diâmetro interno de 77,2 mm, que pode se fixado em qualquer inclinação entre as posições horizontal e vertical. Foram utilizadas válvulas de fechamento rápido para a medição da fração de vazio volumétrica, e sondas de condutividade elétrica para a determinação da velocidade e comprimento das bolhas alongadas. Foram realizados 622 testes, tendo ar como a fase gasosa e uma solução de goma xantana como o fluido de perfuração com comportamento não-newtoniano, para diversas inclinações do tubo. Adicionalmente, foram realizados 619 testes com a mistura ar-água. Os resultados experimentais foram comparados com trabalhos já publicados / Abstract; This work shows the results of an experimental investigation carried out to evaluate the influence of a non-Newtonian medium and pipe inclination on the rise velocity of elongated bubbles. The experimental apparatus consists of a 9.44 m long transparent pipe with 77.2 mm inside diameter that can be set at any inclination from vertical to horizontal. It contains ball valves that allows measurements of volumetric void fractions, and electrical conductivity probes for length and velocity determination of elongated bubbles. For several pipe inclinations, 622 tests were run using air, as the gas phase, and Xanthum Gum solution, as the drilling fluid with non-Newtonian behavior. Additional 619 tests were run for air-water mixture. The experimental results were compared with previous works / Mestrado / Mestre em Engenharia de Petróleo

Orientadores: Luiz Felipe M. Moura, Ibere Nascentes Alves / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica / Made available in DSpace on 2018-07-19T04:52:19Z (GMT). No. of bitstreams: 1 FagundesNetto_JoseRoberto_M.pdf: 3664146 bytes, checksum: d63485aec9f8f87a0d50ebc9f4a03e99 (MD5) Previous issue date: 1993 / Resumo: Uma modelagem para o escoamento bifásico vertical a bolhas e pistonado é proposta, através de uma abordagem Lagrangeana para a fase discreta e um tratamento Euleriano para a fase contínua. Cada bolha é tratada individualmente através de três equações: as de balanço de massa e de quantidade de movimento e a que determina a sua posição em relação a um referencÍal fixo na tubulação. A fase contínua é resolvida através da discretização com malhas defasadas das duas equações de balanço: de massa e de quantidade de movimento. O acoplamento entre as equações das duas fases é realizado através da fração volumétrica e das forças interfaciais. Foi desenvolvido um simulador computacional baseado neste modelo que, devido à abordagem individualizada dada às bolhas, permite que se faça um tratamento estatístico dos resultados. Pode-se determinar, em cada ponto da tubulação, o valor médio e o desvio padrão da altura das bolhas, do comprimento dos pistões de líquido, das velocidades das fases e do gradiente de pressão. Os valores calculados foram comparados com dados pesquisados na literatura e apresentaram um resultado satisfatório tanto para a altura das bolhas como para o comprimento dos pistões de líquido. As velocidades e o gradiente de pressão previstos obtiveram boa correlação com os valores medidos. Foi realizada ainda uma análise de sensibilidade para avaliar como os resultados são influenciados pela variação da vazão volumétrica das fases, do diâmetro do tubo e da pressão na extremidade jusante da tubulação / Abstract: A model to predict the behaviour of vertical two-phase flow regimes, such as bubbles and slug flow has been proposed, using a Lagrangian representation of the discrete phase and a Eulerian representation of the continuous phase. Each bubble is modeled individually through three equations: the mass and the momentum balance equations and the bubble position equation. The continuous phase is solved using a staggered mesh for the mass and momentum equations. Volumetric fraction and interfaces forces provide the coupling of the phase models. A computer program has been developed based on this model and it allows a statistical treatment of the results, due to the individual bubbles approach. It's possible to calculate the average and standard deviation of the bubbles heights, liquid slugs lengths, phase velocities and pressure gradient, for each pipe position. The results were checked against existing data and a reasonable agreement for the lengths was obtained. The velocities and pressure gradient show a good correlation with measured data. It was also performed a sensitive analysis to check the influence of the phases volumetric flow rates, pipe diameter and outlet pressure on the simulation results / Mestrado / Mestre em Engenharia de Petróleo