As Vibrações Induzidas por Escoamentos (VIEs) são dos problemas dinâmicos mais importantes no projeto de trocadores de calor de tipo carcaça e tubos. Esta interação fluido - estrutura pode gerar vibrações de grande amplitude dos tubos ou partes internas do equipamento, o que pode provocar o impacto entre componentes e até falha por fadiga. Várias bancadas experimentais têm sido construídas para o estudo desta interação fluido - estrutura, no entanto, alguns mecanismos de vibração, principalmente aqueles relacionados com escoamentos bifásicos, não são completamente compreendidos. Portanto, neste trabalho, o projeto de uma estrutura dinâmica para o estudo das vibrações induzidas por escoamento bifásico é apresentado. Esta estrutura é composta por um sistema de fios de aço tensionados que permitem que a primeira frequência de ressonância do tubo seja sintonizada. O tubo instrumentado foi instalado numa seção de testes de tubos rígidos em configuração triangular normal com uma razão de passo transversal e diâmetro de 1,26. Os resultados experimentais foram contrastados com uma modelagem teórica desenvolvida para esta estrutura dinâmica. Pôde-se verificar que a primeira frequência de ressonância e seu fator de amortecimento coincidem com os valores estimados pelo modelo, enquanto que a segunda frequência de ressonância apresenta valor teórico 6,6% maior do que o experimental. A seguir, um estudo experimental no ar foi desenvolvido para o reconhecimento dos modos de vibração e as frequências de ressonância da estrutura dinâmica. O fator de amortecimento no ar também foi estimado, o método de Kennedy - Pancu combinado com o algoritmo ERA foi usado para este fim. Assim que a primeira frequência de ressonância foi reconhecida e sintonizada, testes em escoamentos de água foram realizados. A influência da velocidade do escoamento monofásico de água na amplitude de vibração e nas características da Densidade Espectral (PSD) de aceleração foi analisada. Finalmente, testes em escoamento bifásico foram feitos, a sensibilidade do tubo ao escoamento bifásico fica evidente quando analisadas as amplitudes nas PSDs. Além disso, a influência da fração de vazio na amplitude de vibração e no fator de amortecimento foi revisada, no entanto, resultados de estudos anteriores não foram conferidos pelo alto amortecimento concentrado introduzido por alguns elementos construtivos utilizados. / Flow - Induced Vibration (FIV) is the most critical dynamic issue in the design of heat exchangers. This fluid - structure phenomenon may generate high amplitude vibration of tubes or structural parts, which may lead to impacts between internal components or even failure due to fatigue. Many test benches have been constructed to study this fluid - structure interactions, however, some vibration mechanisms, mostly those related to multiphase flow, are not yet fully understood. Therefore, in this work, the design of a dynamic structure for the study of multiphase flow induced vibration is presented. This structure is composed by a system of tensioned piano wires that allow the first natural frequency of the instrumented tube to be calibrated. The tube was installed in a rigid bundle configured in a normal triangular pattern with pitch-todiameter ratio of 1.26. Then, an experimental study in air environment was performed aiming at addressing the mode shapes and resonance frequencies of the dynamic structure. The damping ratio in air environment was also estimated, the Kennedy - Pancu method in combination with the Eigensystem Realization Algorithm have been used to that purpose. Experimental results in air were compared with a theoretical model. The first resonance frequency and its damping factor presented a good agreement with the model, while the second resonance frequency deviated some 6,6% from the predicted value. Tests in water and multiphase flow were also performed. The influence of water flow velocity on vibration amplitudes and on the Power Spectrum Densities (PSDs) has been analyzed. Finally, tests for two - phase water/air flow were carried out; the tube sensitivity to multiphase flow is evident when the vibration level in acceleration PSDs is by analyzed. Also, the influence of void fraction on vibration amplitude and damping ratio have been addressed, however, previous studies could not be validated due to the high damping introduced by some constructive elements of the structure.

Estudou-se a aplicação de uma modelagem RANS, com a utilização do modelo de turbulência k-? na modelagem do escoamento no interior de um sistema de lavagem de gases, composto por um lavador Venturi e uma coluna de spray, por técnicas de fluidodinâmica computacional (CFD), com a fase líquida resolvida pela abordagem Lagrangiana e a atomização secundária dessas por diferentes modelos de quebra. Utilizou-se o pacote comercial de CFD da ANSYS CFX 15.0 com as simulações realizadas em estado transiente. Verificaram-se os resultados obtidos pelo simulador para um lavador Venturi em escala piloto com dados da literatura, obtendo-se boa adequação para a perda de pressão ao longo do equipamento. Comparou-se a influência do tamanho das gotas de líquido da injeção na distribuição de tamanho de gotas ao longo do equipamento utilizando o modelo de quebra CAB (Cascade Atomization and drop Breakup), obtendo-se uma melhor dispersão das gotas de líquido na secção transversal do lavador ao se utilizar na injeção, uma distribuição de tamanho de gotas ao invés de gotas com mesmo diâmetro. Simularam-se diferentes condições de vazão de líquido e gás em um lavador Venturi em escala industrial e compararam-se os resultados de perda de pressão obtidos com correlações da literatura, obtendo-se boa adequação para as regiões iniciais do equipamento. Utilizando o modelo de quebra CAB e a abordagem Euler-Lagrange, comparou-se a distribuição de tamanho de gotas ao longo do lavador Venturi industrial para diferentes vazões de gás. Com a simulação do lavador Venturi, utilizou-se os resultados obtidos neste e simulou-se o sistema de lavagem de gases completo (coluna de Spray e lavador Venturi). Comparou-se ao final a retirada de poluentes por meio da impactação inercial e da difusão, obtendo-se uma melhor probabilidade de retirada de poluentes com o aumento da vazão de líquido. / The flow field inside a wet scrubber was simulated with a RANS model using Computational Fluid Dynamics (CFD) in which the multiphase flow was solved by an Euler-Lagrange approach. Turbulence was modeled by the k-? model and droplet breakup was assumed to occur. The transient flow equations were solved using a commercial CFD package (ANSYS CFX 15.0). Simulations of a pilot plant Venturi scrubber were compared with literature data, in which a good agreement level was achieved for pressure loss through the scrubber. Droplet size distribution was evaluated throughout the scrubber by changing droplet diameter of injected liquid and using a Cascade Atomization and drop Breakup model (CAB). A better liquid dispersion inside the scrubber was achieved for a droplet size distribution injection when compared with a single diameter droplet liquid injection. After the results verification by comparing with data from a pilot scale Venturi, pressure loss for different flow conditions was compared between numerical simulations and literature correlations in an industrial size Venturi scrubber, achieving a good agreement for the entrance section of the scrubber. Droplet size distribution throughout the industrial Venturi scrubber for different gas mass flow conditions were also compared using the CAB model and Euler-Lagrange approach. Removal of contaminants was evaluated by two parameters, being inertial impaction of particulate and diffusion to droplet surface, by which a better contaminant removal was achieved for higher liquid flow rates.

Sensores capacitivos são amplamente utilizados em vários sistemas de medida, tanto na indústria como em laboratórios de pesquisa. Estes dispositivos medem variadas grandezas físicas, como deslocamento, força, pressão, densidade e concentração de substâncias em escoamentos bifásicos, entre outros. Sua aplicação vêm sendo utilizada há algum tempo de forma relativamente elementar, necessitando ainda de estudos específicos para entender a fundo o fenômeno de transferência de carga. A sonda proposta neste trabalho visa melhorar a aquisição do sinal capacitivo e obter uma medida precisa de fração volumétrica em condições bifásicas e trifásicas. Portanto, é necessário um circuito transdutor que transmita de forma precisa a capacitância registrada entre os eletrodos dos sensores até uma central de aquisição. Neste trabalho buscamos desenvolver um modelo de sonda para medir frações volumétricas in situ. Para tanto, apresentamos as etapas de calibração dos sensores, através de válvulas de fechamento rápido, e validamos a técnica sob determinadas condições e variados padrões de escoamento. Dois sensores de geometrias distintas (anéis e hélices) foram utilizados na composição da sonda capacitiva não intrusiva. Propõe-se a solução do sistema linear para medição direta de fração volumétrica em escoamento trifásico. / Capacitive sensors are widely employed in many measurement systems in both industry and laboratory activities. These devices measure several physical quantities such as displacement, forces, pressure, density and concentration of substances in two-phase flows. Their application in industry is still insipient, in part because specific studies are still necessary to properly understand the phenomenon of charge transfer. The goal of this research is the development of a non-intrusive double-sensor probe to on-line measurement of in-situ volumetric fractions in two and three-phase flows. Therefore, we present the calibration work, through quick-closing valves, and validation techniques under certain conditions and for different flows patterns. A transducer circuit to transmit the capacitance sign accurately from the sensors to the central acquisition was developed. Two sensors with different geometries (rings and helix) comprise the double-sensor capacitive probe. Finally, we propose the solution of a linear system for the direct measurement of in-situ volumetric fraction in three-phase flow.

The study in formation and evolution of sprays is essential for developing more detailed physical models and new injection strategies for direct injection internal combustion engines. In the present work, the sprays from a multi-hole injector are evaluated in an effort to characterize the effects of the spray development in a constant surrounding air flow. These interactions are studied in terms of the air turbulence characteristics, the inlet air mass flow and the fuel injection pressure. Ethanol sprays are injected in a constant-flow chamber. The apparatus purpose is to isolate the experiment from the fluid flow properties intrinsic to engine operations, such as instabilities and moving walls. The factors that affects the air-spray interactions were assessed with the air velocity fields in the presence of the ethanol spray. The two-phase particle image velocimetry technique was enhanced to allow measuring in the required experimental conditions. In all conditions, the interaction is based in a pressure gradient formed between the inner and outer regions of the spray. The results indicates a different mechanism when compared to quiescent conditions. The recirculation vortex at the spray border is present only in the initial injection stages. However, the end of injection transient, an instability initiated with the injector needle closing, is still present for these conditions. The interaction mechanism accelerates the velocity distributions towards the spray main boundary. The experiments indicate that the increase of the air mass flow modifies the air penetration velocity but without altering the interaction mechanism characteristics. Higher injection pressures suggests a lower degree of air interaction at the initial instants of the spray development. Turbulent intensity distributions are calculated for the air flow during the injection event. The distributions indicate that the sprays attenuate the turbulent intensity in all conditions, consistent with the observations of the velocity fields. To assess the effects of air turbulence, sets of interchangeable perforated plates are used to limit the integral scales of turbulence. The spectrogram analyses indicate turbulence is reduced not only in the integral scales, but also in all the measured frequency scales. The inlet turbulence integral scales of the air flow have little influence in the spray development. In the turbulence field, the power levels at the end of injection were similar regardless of the inlet turbulence integral scales. / O estudo da formação e evolução dos sprays é essencial para o desenvolvimento de modelos físicos mais detalhados e novas estratégias de injeção para motores de combustão interna de injeção direta. No presente trabalho, os sprays de um injetor multi-furos são avaliados em um esforço para caracterizar os efeitos do desenvolvimento do spray em fluxo de ar constante. Estas interações são estudadas em termos de características de turbulência do ar,do fluxo mássico de ar e da pressão de injeção de combustível. Sprays de etanol são injetados em uma câmara de fluxo constante. O objetivo do aparato é isolar o experimento de propriedades do escoamento intrínsecas ao funcionamento de motores de combustão interna, tais como instabilidades e geometrias móveis. Os fatores que afetam as interações de arspray foram avaliados com os campos de velocidade do ar obtidos na presença de spray. A técnica de velocimetria por imagem de partículas de duas fases foi aprimorada para permitir a medição nas condições experimentais. Em todas as condições, a interação é baseada em um diferencial de pressão formado entre as regiões interna e externa do spray. Os resultados indicam um mecanismo diferente quando comparado com condições quiescentes. O vortex formado na fronteira do spray é observado apenas nos estágios iniciais de injeção. No entanto, o transiente de fim de injeção ainda está presente para essas condições. O mecanismo de interação acelera as distribuições de velocidade em direção à fronteira do jato. Os experimentos indicam que o aumento do fluxo de massa de ar modifica a velocidade de penetração do ar, mas sem alterar as características do mecanismo de interação. Distribuições de intensidade turbulenta são calculadas para o fluxo de ar durante o evento de injeção. As distribuições indicam que os sprays atenuam a intensidade turbulenta em todas as condições, consistente com as observações dos campos de velocidade. Para avaliar os efeitos da turbulência do ar, conjuntos de placas perfuradas intercambiáveis são utilizadas para limitar as escalas integrais de turbulência. As análises do espectrograma indicam que a turbulência é reduzida não apenas nas escalas integrais, mas também em todas as escalas de freqüência medidas. Estas escalas integrais de turbulência do fluxo de ar de entrada têm pouca influência no desenvolvimento do spray. No campo de turbulência, os níveis de potência ao final da injeção foram semelhantes, independentemente das escalas integrais de turbulência de entrada.

Orientador: Milton Mori, Marcos A. d'Avila / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-06T05:47:36Z (GMT). No. of bitstreams: 1 Santos_CelsoMurilodos_M.pdf: 3396606 bytes, checksum: 6f04eb0860468401aa6369dba2b5c839 (MD5) Previous issue date: 2005 / Resumo: Colunas de bolhas são equipamentos de simples operação, mas que apresentam uma fluidodinâmica complexa. Seu funcionamento básico é dado pela injeção de um gás (fase dispersa) na base da coluna em um meio líquido (fase contínua). Elas são empregadas nas mais diversas áreas industriais, como processos de química fina, reações de oxidação, reações de alquilação, síntese de Fischer-Tropsch, tratamento de efluentes, reações de fermentação e produção de proteínas, e mais recentemente, em cultura de células. Neste trabalho foram realizadas simulações tridimensionais do escoamento transiente em coluna de bolhas utilizando o pacote computacional de CFD CFX® 5.7 (ANSYS®). Este pacote utiliza o método dos volumes finitos para solucionar as equações governantes de transporte em escoamento multifásico, utilizando a abordagem Euleriana-Euleriana, que mostrou-se adequada para este estudo. Os estudos de casos foram feitos para verificar a influência das propriedades da fase contínua (viscosidade e tensão superficial), da fase dispersa (diâmetro das bolhas) e dos modelos de arraste ... Observação: O resumo, na íntegra, poderá ser visualizado no texto completo da tese digital / Abstract: Bubble columns are equipments of simple operation, but presents a complex fluidodynamic. Its basic operation consists of the injection of a gas (disperses phase) in the bottom of the column in a liquid medium (continuous phase). As reactors, they are used in a variety of chemical processes, such as Fischer-Tropsch synthesis, oxidation reactions, alkylation reactions, effluent treatment, fermentation, reactions, protein production, and more recently, in cell cultures. In this work, three-dimensional simulations of the transient flow in bubble columns was performed using the computational package CFD CFX® 5.7 (ANSYS®). This package uses the finite volumes method to solve the governing equations of transport in multiphase flow using the Eulerian-Eulerian approach, which showed to be adequate for this study. The case studies were carried out in order to verify the influence of the properties of the continuous phase (viscosity and superficial tension), the dispersed phase (diameter of the bubbles) and the drag models ... Note: The complete abstract is available with the full electronic digital thesis or dissertations / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Orientador: Milton Mori / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-12T20:09:08Z (GMT). No. of bitstreams: 1 Hodapp_MaximilianJoachim_M.pdf: 3485880 bytes, checksum: 79b4777bac2b444262795c5ddaaa44ac (MD5) Previous issue date: 2009 / Resumo: O objetivo deste trabalho foi o estudo comparativo de duas modelagens para representação de um escoamento gás-sólido. Na primeira modelagem avaliou-se uma correlação de arraste apresentada recentemente na literatura, baseada em simulações lattice Boltzmann, aplicada a escoamentos gás-sólido. Na segunda observou-se o efeito da variação do coeficiente de especularidade na condição de contorno na parede sobre os perfis de velocidade da fase particulada. Grande atenção tem sido dada a modelagem matemática de escoamentos multifásicos, em especial ao gás-sólido, uma vez que vários processos industriais utilizam-se desta operação. Porém o desenvolvimento de modelos analíticos que incluam todos os fenômenos de transferência de massa, energia e quantidade de movimento não encontra-se disponível, devido principalmente a grande complexidade dos fenômenos envolvidos. Neste aspecto a fluidodinâmica computacional (CFD) tem demonstrado ser uma boa alternativa para o estudo de sistemas complexos, sendo diversos estudos, não somente de engenharia, aplicando esta técnica, publicados todos os anos. Como forma de validar os resultados obtidos por este método numérico, escolheu-se um caso de estudo em escala de laboratório. Os softwares comerciais ANSYS CFX 10 e FLUENT 6.3 foram utilizados para a definição e resolução do problema, além do pós-processamento. Os resultados obtidos foram comparados com os dados numéricos de um trabalho de mestrado do PQGe, bem como com dados experimentais da literatura. Pode-se perceber que os resultados, para a primeira abordagem não apresentaram melhoras em relação a outras modelagens das forças entre e intra-partículas, além do maior tempo computacional requerido. A segunda abordagem demonstrou valores adequados para o coeficiente estudado / Abstract: The aim of this work was a comparative study of two different modeling to represent a gas-solid flow. The first approach consists of a new drag correlation presented in the literature. This relation was obtained through lattice Boltzmann simulations of gas-solid flow, thus not depending on empirical data. The second looked for the effects of the variation of the specularity coefficient at the wall boundary condition. Multiphase flow modeling is gathering great attention, especially to gas-solid flows, due to its importance in industrial processes. However analytical models that take into account the mass, momentum and energy transfers are not available, mainly because of the complexities evolved in such systems. Therefore Computational Fluid Dynamics (CFD) has proved to be a viable alternative, having a large number of scientific works been published in recent years. In order to validate the results a comparison with other simulations using different modeling, done by another member of the PQGe laboratory, and with experimental data was carried out. The commercial softwares ANSYS CFX 10 and FLUENT 6.3 were used to define and numerically solve the problem, also to post process the results. For the first approach, the comparison showed that the studied drag correlation gave no improvement upon the other two models analyzed. Also a longer computational time was required, which can not be ignored as an important parameter in CFD simulations. As for the second approach, it was possible to obtain adequate values for the specularity coefficient / Mestrado / Desenvolvimento de Processos Químicos / Mestre em Engenharia Química

Orientador: Antonio Carlos Bannwart / Dissertação (mestrado) - Universidade Estadual de Campinas, Faculdade de Engenharia Mecanica, Instituto de Geociencias / Made available in DSpace on 2018-08-14T16:09:37Z (GMT). No. of bitstreams: 1 Lima_IvanNovilleRochaCorrea_M.pdf: 2854489 bytes, checksum: fec392ec12651861ca4547ea5ce279c9 (MD5) Previous issue date: 2009 / Resumo: O presente trabalho tem por objetivo central a medição experimental do coeficiente de transferência de calor bifásico no padrão gás-líquido horizontal intermitente, por ser este um padrão muito comum nas operações de produção de petróleo. Para alcançar esse objetivo, foi construído um aparato experimental consistindo de uma seção de testes alimentada por uma mistura pré-aquecida de ar e água, a diferentes pares de vazões e temperaturas de entrada. Para realizar a troca térmica, o tubo foi circundado por uma jaqueta de água gelada submetida a diversas vazões e temperaturas de entrada, formando, com o tubo central, uma seção de troca térmica de correntes paralelas. A metodologia experimental consistiu em medir a taxa de transferência de calor do fluido interno para o fluido externo e, a partir da diferença média de temperatura entre o fluido quente e a parede, determinar o coeficiente de transferência de calor entre o fluido quente e a parede, para cada elenco de condições de operação. O procedimento foi realizado tanto para os ensaios em escoamento monofásico como para o escoamento bifásico. Os resultados para o coeficiente de transferência de calor do fluido quente foram comparados às correlações e aos modelos existentes e apresentaram resultados satisfatórios. / Abstract: The main objective of the present work is the experimental measurement of the two-phase heat transfer coefficient in the intermittent horizontal gas-liquid flow pattern. This flow pattern is very common in oil production pipelines. In order to reach this objective, an experimental apparatus was built, consisting of a test section fed with an air-water mixture pre-heated at different pairs of inlet flow rates and temperatures. In order to realize the heat transfer, the pipe was surrounded by a jacket containing cold water under various inlet flow rates and temperatures, forming with the central pipe a heat transfer section of parallel current. The experimental methodology consisted of measuring the heat transfer rate from the internal to the external fluid, and from the average difference of the temperature between the hot fluid and the pipe wall, determine the heat transfer coefficient between the hot fluid and the pipe wall for each cast of operating condition. The procedure was carried through for both single phase and two-phase tests. The results of the heat transfer coefficients between the hot fluid were compared with some correlations and existing models and presented satisfactory results. / Mestrado / Explotação / Mestre em Ciências e Engenharia de Petróleo

Orientador: Milton Mori / Tese (doutorado) - Universidade Estadual de Campinas, Faculdade de Engenharia Quimica / Made available in DSpace on 2018-08-10T23:15:40Z (GMT). No. of bitstreams: 1 Rosa_LeonardoMachadoda_D.pdf: 1367270 bytes, checksum: 1819e68c66a949c1c8fae024725df552 (MD5) Previous issue date: 2008 / Resumo: A aplicação de Fluidodinâmica Computacional (CFD) em operações industriais vem aumentando sensivelmente, em razão desta auxiliar no projeto de novos equipamentos, bem como a otimização dos já existentes. O presente trabalho apresenta uma série de simulações tridimensionais, bifásicas e adiabáticas, utilizando-se o modelo de dois-fluidos para o cálculo do escoamento gás-sólido em risers de diferentes tamanhos ¿ um com as dimensões de 10m de altura e 76mm de diâmetro interno, e outro com 18m de altura e 50mm de diâmetro interno. Utilizou-se o modelo k-epsilon para tratar a fase gasosa como turbulenta. Foram considerados os modelos cinéticos de 4 lumps e 10 lumps para o cálculo das taxas de craqueamento catalítico presentes no reator. Dados experimentais para comparação do escoamento foram obtidos em literatura, e os resultados obtidos mostram que, para a correta predição dos fenômenos presentes no interior de um riser, deve-se levar em consideração também as reações de craqueamento, visto que a alteração da composição química da fase gasosa interfere diretamente no comportamento dinâmico do escoamento / Abstract: The application of Computational Fluid Dynamics (CFD) in industrial operations has been increased not only in the project of new equipments but also in the optimization of the existing ones. The present work shows a set of tridimensional two-phase and adiabatic conditions, using the two-fluids model for the calculation of gas-solids flow in a high-flux risers of two different geometries, one with 10m of height and 76mm of internal diameter, and the second with 18m of height and 50mm of internal diameter. The k-epsilon model was applied to consider the influence of turbulence over the gas phase. The four-lump and ten-lump schemes kinetic models were used to predict the catalytic cracking reactions in the two different reactors. The riser with 10m of height and 76 mm of internal diameter had two entrances, one for particles inlet and one for the organic materials. The size of mesh for this case was 584.000 control volumes. The second riser studied with geometry of 18m of height and 50mm of internal iameter had one entrance for vapor, one for organic materials, and a third for particles inlet. The mesh generated over this geometry was considered 650.000 control volumes. The numerical results of gas-solids dispersion flow in a high-flux were compared with the xperimental data from the work of P¨arssinen and Zhu (2001), the four lump scheme for he mass concentration of chemical species was compared with the work of (ALI et al., 997), and the ten lump scheme was compares with the work of NAYAK et al. (2005). The catalytic cracking reactions models show an important role in the prediction of all henomena present in the reactor. The simulated results show that the kinetics models are important not only because of the variation of the gas phase and the temperature profiles inside the reactor but also because these parameters interfere directly in the dynamic behavior of the gas-solid flows / Doutorado / Desenvolvimento de Processos Químicos / Doutor em Engenharia Química

Neste trabalho, é apresentada uma análise de estabilidade linear para um sistema flowline-riser. O modelo utiliza as equações de continuidade para as fases líquidas e gasosas e uma equação de momento para a mistura onde são considerados os efeitos de atrito. A mistura água-óleo é modelada de forma homogênea. Para a determinação da fração de vazio, adota-se o modelo de fluxo de deriva, baseado em várias correlações drift. Se o escoamento é estratificado a fração de vazio é modelada através do modelo de equilíbrio local de Taitel e Dukler (1976). Para a caracterização dos fluidos é adotado um modelo de equilíbrio de fase black-oil onde a transferência de massa é considerada entre as correntes óleo e gás em função das condições locais de pressão e temperatura. É considerada a abordagem de parâmetros distribuídos, onde diversas geometrias com discretização variável podem ser consideradas tais como riser em catenária e riser lazy wave. Para realizar a análise de estabilidade linear, as equações do modelo são linearizadas em torno do estado estacionário e discretizadas pelo método das diferenças finitas, onde foi utilizado um programa escrito em Matlab. A partir do sistema linearizado, é avaliada a estabilidade do estado estacionário pelas raízes do polinômio característico, que são solução do problema de autovalores e autovetores. É avaliada a convergência numérica do modelo e mapas de estabilidade são apresentados para vários sistemas de produção de petróleo. Os resultados obtidos numericamente são comparados com pontos operacionais de sistemas reais de produção de petróleo. Observou-se que o modelo convergiu para todos os casos avaliados e apresentou ótima concordância com os dados de campo. / This work presents a linear stability analysis for a flowline-riser system. The model considers continuity equations for the liquid and gas phases, and a momentum equation for the mixture that accounts for friction effects. The water-oil mixture is modeled as being homogeneous. The void fraction is determined by using the drift-flux model based on several drift correlations. For stratified flow, the void fraction is expressed by the local equilibrium model of Taitel and Dukler (1976). Fluid characterization uses a black-oil model that considers mass transfer between the oil and gas flows as a function of the local pressure and temperature conditions. In the proposed approach with distributed parameters, several geometries with variable discretization can be considered, such as catenary and lazy wave risers. To perform the linear stability analysis, the equations of the model are linearized around the stationary state and discretized using the finite difference method, implemented using custom-written code in Matlab. From the linearized system, the stability of the steady state is evaluated by computing the roots of the characteristic polynomial equation of the eigenvalues and eigenvectors problem. Convergence of the numerical model is evaluated and stability maps for several oil production systems were presented. Numerical results were compared with the actual measurements of oil production systems. The model converged in all cases tested and presented an excellent agreement with field data.

O trabalho tem como objetivo realizar um estudo experimental na bancada do Laboratório Multipropósito de Escoamento Multifásico com o intuito de estudar o fenômeno de intermitência severa em um sistema pipeline-riser com os fluidos água e ar. A intermitência severa pode ocorrer em sistemas pipeline-riser onde o pipeline é descendente seguido de um riser vertical, além de ser necessário a presença de baixas vazões de fluidos. Esse fenômeno é caracterizado por ser cíclico em que há acúmulo de líquido na base do riser e por acabar causando perdas na produção de petróleo e gás devido a grandes flutuações de pressão e vazão que podem durar horas a depender do comprimento do sistema. Os picos de pressão e vazão também podem causar o desligamento do sistema de separação na plataforma. O estudo se dividiu em várias etapas onde inicialmente foram definidos os procedimentos experimentais a serem utilizados. A calibração de placas de orifício foi necessária para o controle de vazão mássica de gás. Os resultados experimentais foram divididos em casos estáveis e instáveis onde a instabilidade é caracterizada pela presença de ciclos de pressão que podem ser observados em históricos de pressão na base do riser. Mapas de estabilidade foram criados e a região instável obtida experimentalmente foi comparada com a curva de estabilidade obtida pelo uso da teoria de estabilidade linear. Os históricos de pressão na base do riser para os casos instáveis obtidos foram comparados com dois modelos numéricos. A variação da pressão no separador foi usada para verificar a mitigação da intermitência severa e/ou da condição instável obtida no sistema para alguns casos instáveis. / The objective of this work is to do an experimental study of the severe slugging phenomenon in the pipeline-riser system of the Multipurpose Multiphase Flow Laboratory by using the fluids air and water. Severe slugging may occur for low flow rates in pipeline-riser systems where a downward pipeline is followed by a vertical riser. In this phenomenon there is liquid accumulation at the bottom of the riser resulting in production losses due to the great fluctuations of pressure and flow rate during its cycles which may last for hours depending on the length of the system. The high pressure values can also cause shutdown of the platform separation system. The first stage of this study was to define the experimental procedure to be adopted. It was necessary to perform a calibration of the orifice plates in order to have a precise control of the gas mass flow rate. The experimental results were divided in stable and unstable cases. The instability is defined by the presence of pressure oscillations at the bottom of the riser. Stability maps were created to compare the stabiliy curve obtained by the stability linear theory with the experimental results. The experimental pressure oscillations were compared with two numerical models. The pressure variation at the separator was studied to verify the mitigation effects during unstable and/or severe slugging conditions.