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.

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior / Em uma grande gama de problemas físicos, governados por equações diferenciais, muitas vezes é de interesse obter-se soluções para o regime transiente e, portanto, deve-se empregar técnicas de integração temporal. Uma primeira possibilidade seria a de aplicar-se métodos explícitos, devido à sua simplicidade e eficiência computacional. Entretanto, esses métodos frequentemente são somente condicionalmente estáveis e estão sujeitos a severas restrições na escolha do passo no tempo. Para problemas advectivos, governados por equações hiperbólicas, esta restrição é conhecida como a condição de Courant-Friedrichs-Lewy (CFL). Quando temse a necessidade de obter soluções numéricas para grandes períodos de tempo, ou quando o custo computacional a cada passo é elevado, esta condição torna-se um empecilho. A fim de contornar esta restrição, métodos implícitos, que são geralmente incondicionalmente estáveis, são utilizados. Neste trabalho, foram aplicadas algumas formulações implícitas para a integração temporal no método Smoothed Particle Hydrodynamics (SPH) de modo a possibilitar o uso de maiores incrementos de tempo e uma forte estabilidade no processo de marcha temporal. Devido ao alto custo computacional exigido pela busca das partículas a cada passo no tempo, esta implementação só será viável se forem aplicados algoritmos eficientes para o tipo de estrutura matricial considerada, tais como os métodos do subespaço de Krylov. Portanto, fez-se um estudo para a escolha apropriada dos métodos que mais se adequavam a este problema, sendo os escolhidos os métodos Bi-Conjugate Gradient (BiCG), o Bi-Conjugate Gradient Stabilized (BiCGSTAB) e o Quasi-Minimal Residual (QMR). Alguns problemas testes foram utilizados a fim de validar as soluções numéricas obtidas com a versão implícita do método SPH. / In a wide range of physical problems governed by differential equations, it is often of interest to obtain solutions for the unsteady state and therefore it must be employed temporal integration techniques. One possibility could be the use of an explicit methods due to its simplicity and computational efficiency. However, these methods are often only conditionally stable and are subject to severe restrictions for the time step choice. For advective problems governed by hyperbolic equations, this restriction is known as the Courant-Friedrichs-Lewy (CFL) condition. When there is the need to obtain numerical solutions for long periods of time, or when the computational cost for each time step is high, this condition becomes a handicap. In order to overcome this restriction implicit methods can be used, which are generally unconditionally stable. In this study, some implicit formulations for time integration are used in the Smoothed Particle Hydrodynamics (SPH) method to enable the use of larger time increments and obtain a strong stability in the time evolution process. Due to the high computational cost required by the particles tracking at each time step, the implementation will be feasible only if efficient algorithms were applied for this type of matrix structure such as Krylov subspace methods. Therefore, we carried out a study for the appropriate choice of methods best suited to this problem, and the methods chosen were the Bi-Conjugate Gradient (BiCG), the Bi-Conjugate Gradient Stabilized (BiCGSTAB) and the Quasi-Minimal Residual(QMR). Some test problems were used to validate the numerical solutions obtained with the implicit version of the SPH method.

CNPq / O bloqueio de linhas de produção devido à formação de hidratos é uma das principais preocupações na garantia do escoamento em operações de produção de óleo e gás devido aos altos custos associados às paradas de produção. Os hidratos são formados pelo aprisionamento de moléculas de gás em uma gaiola de moléculas de água formada por ligações de hidrogênio. As condições de alta pressão e baixa temperatura necessárias para a formação de hidratos são frequentemente encontradas em cenários de produção no mar, sendo que o escoamento em golfadas é considerado o padrão de escoamento predominante. O presente trabalho utiliza uma abordagem em regime permanente para modelar a formação de dispersões homogêneas de hidrato-em-água durante o escoamento em golfadas em tubulações horizontais. O consumo das fases para a formação de hidratos é estimado em função da área interfacial gás-água e da temperatura de subresfriamento do sistema. O modelo de escoamento em golfadas é acoplado com: (i) termos de transferência de massa para recalcular as velocidades das estruturas das golfadas e (ii) termos de geração de calor, visto que a formação de hidratos é um processo exotérmico. O modelo foi implementado em linguagem Fortran90, utilizando lógica nodal de marcha. Comparações do modelo com dados experimentais da literatura apresentam desvios de aproximadamente ±20%. O modelo foi utilizado para analisar os efeitos da formação de hidratos sobre a hidrodinâmica e transferência de calor do escoamento em golfadas, analisando: velocidades (superficiais e das estruturas do escoamento em golfadas), pressão, temperatura, coeficiente de transferência de calor e geometria da célula unitária do escoamento em golfadas (comprimento das regiões características e frações de fase). O modelo também foi utilizado para analisar a influência das condições de entrada do escoamento (velocidade superficial da mistura, razão de líquido/mistura, pressão, uso de inibidores químicos) e de geometria da tubulação (diâmetro interno, espessura da parede, condutividade da parede) que retardam a formação de hidratos. Ao final, a influência da deposição de camadas de hidratos sobre o escoamento em golfadas (pressão, temperatura, velocidades, geometria da geometria da célula unitária) é analisado. / Hydrate formation is one of the main flow assurance concerns in offshore oil and gas production due to the high cost of production interruptions or impairments. Hydrates are formed by the imprisonment of gas molecules into hydrogen bonded cages of water molecules. The high pressure and low temperature conditions needed for hydrate formation are frequently found in offshore production scenarios, where slug flow is often the prevailing flow regime. The present work uses a steady-state approach for modeling homogeneous hydrate-in-water dispersion formation on horizontal slug flow in pipelines. The consumption of the phases during hydrate formation is estimated by a kinetic model in terms of the gaswater interfacial surface and the subcooling of the system. The slug flow model is coupled with: (i) mass transfer terms so as to recalculate the velocities of the slug flow structures and (ii) heat generation terms, since hydrate formation is an exothermic process. Comparisons of the model with experimental data present an average deviation of ±20%. The model was used to analyze the effects brought by hydrate formation over the slug flow hydrodynamics and heat transfer, analyzing: velocities (superficial and structure velocities), pressure, temperature, heat transfer coefficient and unit cell geometry (region lengths and phase fractions). The model was also used to analyze the influence of the pipe inlet conditions (mixture superficial velocity, liquid loading, pressure, use of chemical inhibitors) and of the pipeline geometry (internal diameter, wall width, wall thermal conductivity) in delaying hydrate formation. Finally, the influence of hydrate deposits on the slug flow behavior (pressure, temperature, velocities, unit cell geometry) is analyzed.

O fenômeno de perda de circulação é consequência da perfuração de poços de petróleo e gás em substratos permeáveis e é caracterizado pelo influxo de fluido de perfuração da região anular para a formação. Além dos custos adicionais relacionados com a necessidade constante de reposição do fluido perdido, danos à formação nas adjacências do poço são de ocorrência comum. Uma potencial solução é a adição de materiais de perda de circulação (LCM’s) ao fluido de perfuração para que uma camada de reboco se forme na parede do poço em filtração dinâmica. Neste trabalho, o escoamento bifásico líquido-sólido é simulado numericamente através de uma abordagem de Euler-Lagrange para representar o crescimento do reboco na interface poço-formação. A região anular do poço é idealizada como um canal vertical de seção transversal retangular em contato com um meio poroso anisotrópico, o qual é concebido na escala do poro (heterogêneo) como um arranjo de cilindros alternados. O escoamento particulado é ascendente e eventualmente as partículas, as quais representam o LCM, abandonam o canal adentrando o meio poroso onde pode ocorrer a geração do reboco. A simulação é realizada através do Método de Fase Discreta Densa (DDPM) no qual as equações da fase sólida e fluida são solucionadas separadamente. O Método dos Elementos Discretos (DEM) é empregado para calcular as interações de contato (colisão e atrito). A caracterização do combate à perda de circulação é realizada através da redução da vazão de fuga ao longo do tempo. Resultados mostram o efeito da variação do número de Reynolds no canal (125, 250, 500), da vazão inicial de fuga (5, 10, 20%), do diâmetro das partículas (0,50; 0,75; 1,00 mm), da razão de massa específica partícula-fluido (1,5; 2,5) e da configuração do meio poroso (porosidade, número de cilindros e garganta de poro). A eficiência do reboco na redução da vazão de fuga recai na formação de um plugue de partículas que seja capaz de obturar toda a altura do meio poroso, não permitindo que haja gargantas de poros desobstruídas. Tal aspecto é favorecido pela redução do diâmetro e da razão de massa específica, bem como pelo aumento do Reynolds e da vazão inicial de fuga. / The lost circulation is a consequence of the well drilling in a permeable substrate, being characterized by the drilling fluid influx to the porous formation. Despite the additional costs associated with the need of continuous fluid replacement, another misfortune is the irreversible formation damage due to the fluid invasion, which may reflect negatively throughout the productive life of the well. Therefore, it is of utter importance to make use of preventive and corrective techniques, like the addition of lost circulation materials (LCM) to the drilling fluid. The LCM particles deposit over the porous formation under dynamic filtration and create a mud cake (filter cake) that helps diminishing the fluid invasion flow rate. In this work, the liquid-solid two-phase flow is numerically simulated via an Euler-Lagrange approach to represent the mud cake growth. The well annular region is considered as a vertical channel bounded by an anisotropic porous formation. The porous medium is conceived in the pore-scale as a periodic array of staggered cylinders. The fluid flows upward through the channel carrying the solid particles that mimic the LCM’s. The particles might eventually get into the porous formation and deposit, creating the mud cake. The equations for the fluid flow and for the particles movement are solved separately via the Dense Discrete Phase Model (DDPM). The particle-particle and particle-wall interactions like collision and friction are calculated via the Discrete Element Methods that is coupled to the DDPM. The effectiveness of the lost circulation mitigation is evaluated mainly by the decrease in the fluid invasion flow rate over time. The results demonstrate the influence of the variation of the Reynolds number on the channel (125, 250, 500), the initial fluid invasion flow rate (5, 10, 20%), the particles diameter (0.50, 0.75, 1.00 mm), the particle-fluid specific mass ratio (1.5, 2.5) and the configuration of the porous medium (porosity, number of cylinders and pore throat).

Escoamentos bifásicos líquido-gás estão presentes na natureza e em muitas atividades industriais. Neste tipo de escoamento, as fases líquida e gasosa podem assumir diferentes configurações espaciais dentro da tubulação, chamadas padrões de escoamento. O escoamento bifásico líquido-gás em golfadas é o padrão de escoamento mais frequente nas aplicações industriais, ocorrendo em uma ampla faixa de velocidades das fases segundo os estudos de diversos autores. A modelagem matemática para o escoamento em golfadas compreende desde modelos simples em regime estacionário até modelos mais complexos, em regime transiente. E, para solução destes modelos são necessárias correlações empíricas e distribuições estatísticas dos parâmetros característicos do escoamento. Assim, no presente trabalho, vários modelos baseados em redes neurais artificiais são apresentados como suporte à caracterização dos parâmetros do escoamento bifásico em golfadas em função das séries temporais de fração de vazio obtidas experimentalmente. As séries temporais de fração de vazio são medidas com um par de sensores de malha de eletrodos instalado na seção de testes de uma planta experimental do NUEMUTFPR e descritas em Castillo (2015). A partir das séries temporais de fração de vazio medidas são calculados os parâmetros de interesse para o escoamento em estudo: comprimento da bolha alongada de gás, comprimento do pistão de líquido, velocidade de translação da bolha alongada e desvios padrões para essas variáveis. Essas variáveis medidas e calculadas são utilizadas para a obtenção de um conjunto de modelos baseados em rede neural artificial. Após obtenção dos modelos é realizado um estudo de simulação no qual esses modelos são usados para estimar os parâmetros que caracterizam o escoamento bifásico em golfadas. Análises detalhadas dos resultados mostraram que as variáveis relacionadas à fase gasosa são estimadas com maior acurácia que as variáveis relacionadas à fase líquida. Como aplicação imediata do modelo obtido, apresenta-se sua utilização como uma ferramenta de cálculo das condições iniciais para um modelo matemático fenomenológico de escoamento bifásico em golfadas com leve mudança de inclinação baseado no método de seguimento de pistões. O diferencial do presente trabalho está na predição da característica intermitente do escoamento bifásico líquido-gás em golfadas a partir do modelo neural, além da estimação de parâmetros médios para as variáveis de interesse com taxas de incerteza variando entre 10% e 16%. / Gas-liquid two-phase flows are present in nature and in different industrial activities. In this type of flow, the liquid and gas phases assume different spatial configurations inside the pipe, called flow patterns. Slug flow is one of the most frequent flow patterns in industrial applications, occurring over a wide range of phase velocities according to studies presented by several authors. The mathematical modelling of slug flow comprises from simple steady state models to more complex models for transient regimes. Those models require closure relationships, e.g. empirical correlations and statistical distributions of characteristic flow parameters. In this work, several models based on artificial neural networks are presented as a support to the characterization of the two-phase slug flow parameters that depend on experimentally obtained void fraction time series. The void fraction time series are measured with a pair of wiremesh sensors installed in a test section of an experimental rig in the premises of the NUEM/UTFPR labs and described in Castillo (2015). From the time series of void fraction measurements relevant parameters to the flow under consideration are computed: the length of the elongated gas bubble, the liquid slug length, the translational velocity of the elongated bubble and the standard deviations for those variables. Those measured and calculated variables are used to obtain a set of artificial neural network-based models. After obtaining such models, a simulation study in which those models are used to estimate the parameters that characterize the two-phase slug flows is carried out. Detailed analysis of the results showed that the variables related to the gas phase are estimated with greater accuracy than the ones related to the liquid phase. As an immediate application of the obtained model, its use as a tool to calculate the initial conditions for a phenomenological mathematical model of twophase slug flow with a slight change of inclination based on a slug tracking method is presented. The differential of this study is to predict the intermittent features of the twophase slug flow by means of a neural model, as well as the estimation of average parameters for the variables of interest with uncertainly rates ranging between 10% and 16%.

Petrobras / No processo de perfuração de poços de petróleo e gás o controle da pressão é uma importante tarefa. Se a pressão no interior do poço estiver abaixo da pressão de poros, um influxo da formação pode ocorrer, fenômeno denominado de kick. Se o influxo não for controlado pode ocorrer um fluxo descontrolado da formação para a superfície (blowout). As pressões de fechamento são utilizadas para calcular a pressão do reservatório, entretanto, o caráter viscoplástico do fluido de perfuração atenua as pressões medidas. No presente trabalho, é apresentada uma modelagem matemática e numérica para prever a propagação de pressão ao longo do poço durante um influxo de gás (kick). Considerou-se a compressibilidade e o comportamento viscoplástico do fluido de perfuração, modelando-o como Plástico de Bingham. O escoamento é considerado unidimensional, laminar, transiente, isotérmico e homogêneo. A solubilidade do gás no fluido de perfuração é desconsiderada e o gás é modelado como um gás real. O fluxo da formação para o interior do poço é tratado através da lei de Darcy. Os balanços de massa e de quantidade de movimento para a mistura formam um sistema de equações diferenciais parciais hiperbólicas, tendo como incógnitas a pressão e a velocidade, as quais são resolvidas pelo método das características. Consideram-se dois casos de estudo: o kick dinâmico e o kick estático. A obtenção dos campos de velocidade e pressão ao longo do poço é realizada através de um programa computacional desenvolvido em linguagem FORTRAN. Posteriormente, aferiu-se os resultados obtidos pelo modelo desenvolvido comparando-os com solução analítica e dados experimentais. Os resultados apresentam que quanto maior a tensão limite de escoamento, menor é a transmissão de pressão e que se a pressão do reservatório não for grande o suficiente, não há ganho de volume de fluido nos tanques de lama. / In the process of drilling oil and gas wells the pressure control is an important task. If the pressure inside the well is smaller than the formation pressure, an influx of the formation may occur, phenomenon denominated as kick. If the influx is not controlled there may be an uncontrolled flux from the formation to the surface (blowout). The closing pressures are used to calculate the formation pressure, however, the viscoplastic character of the drilling fluid reduces the pressures measured on the surface. This current work presents a mathematical and numerical modeling to predict the pressure propagation along the well during a gas influx (kick). The compressibility and the viscoplastic behavior of the drilling fluid were considered, modeling it as a Bingham Plastic. The flow is considered as one-dimensional, laminar, transient, isothermal and homogeneous. The solubility of the gas in the drilling fluid is disregarded and the gas is modeled as a real gas. The flux of the formation into the wellbore is treated by Darcy’s law. The balance equations of mass and momentum for the mixture result in a system of hyperbolic partial differential equations, having as unknowns the pressure and the velocity, which are solved by the method of characteristics. Two study cases were considered: the dynamic kick and the static kick. The pressure and velocity fields along the well are obtained by a computer program developed in FORTRAN language. Afterwards, the results obtained from the model developed were assessed by comparing to analytic solution and experimental data. The results present that the bigger the yield stress is, the smaller is the pressure transmission and that if the formation pressure is not high enough, there is no pit-gain in the mud tanks.

Petrobras / Petróleos parafínicos são complexas misturas de hidrocarbonetos saturados, aromáticos, naftênicos e, em menor quantidade, de componentes que apresentam heteroátomos como resinas e asfaltenos. Em baixas temperaturas, a solubilidade das parafinas em solução diminui e ocorre a precipitação de hidrocarbonetos no petróleo. Devido à cristalização nota-se não apenas o aumento considerável da viscosidade dos petróleos em baixas temperaturas, mas também a deposição de parafina nas paredes internas da tubulação e a gelificação do petróleo nos oleodutos em eventuais paradas da produção. Os trabalhos disponíveis na literatura mostram que os históricos térmico e de cisalhamento influenciam o comportamento mecânico de petróleos parafínicos em baixas temperaturas. No presente trabalho são apresentados resultados da influência da temperatura inicial de resfriamento na temperatura de cristalização, na temperatura de gelificação e na tensão crítica necessária para reiniciar o escoamento de petróleos parafínicos. A temperatura em que o material é submetido antes do resfriamento influencia nas propriedades reológicas do petróleo. A partir da revisão de conceitos de cristalização define-se que para a precipitação do primeiro cristal é necessário a ocorrência de super-resfriamento na solução. Em outras palavras, não é possível que ocorra a cristalização da parafina exatamente na temperatura de equilíbrio sólido-líquido, Teq,SL. É proposta uma metodologia experimental para determinar a temperatura de saturação de óleos parafínicos. A partir de resultados reométricos pode-se notar que a magnitude do superresfriamento influencia a taxa de nucleação e o tamanho final dos cristais precipitados. O tamanho da região metaestável é influenciado pela taxa de resfriamento e pela temperatura inicial do teste. A partir de microscopia de campo claro observa-se que quanto maior o superresfriamento maior o número de cristais formados. Pode-se adiantar que a morfologia e o tamanho dos cristais formados influenciam o comportamento reológico do óleo parafínico em baixas temperaturas. / Waxy crude oil is a complex mixture of hydrocarbons consisting of paraffins, aromatics, naphthenes, asphaltenes and resins. At low temperatures the solubility of high molecular weight components in the oil is decreased and especially the n-paraffins tend to precipitate as crystal structures. These crystals provide a non-Newtonian behavior to the fluid, tend to deposit in the inner surface of pipelines and are responsible for the gelation of the oil when the flow is interrupted. In such cases, the pressure needed to restart the oil flow in subsea pipelines can be much larger than the usual steady-state pressure as the temperature in such environment can be as low as 4 °C. The literature has shown that not only the temperature itself but also the fluid shear and thermal histories have significant influence on the yield stress of waxy crude oils. The current work shows the effect of the initial cooling temperature on the waxy crude oil viscosity, gelation temperature and critical stress. It is showed that the initial cooling temperature affects the mechanical behavior of crude oils at low temperature. It is presented that crystallization is the process where an ordered solid structure is formed from a disordered phase. This process involves two main stages, namely nucleation and crystals growth. As paraffins nucleation is a stochastic process, a supercooling is needed to the onset of this phenomenon. In other words precipitation of crystals could not take place at the solidliquid equilibrium thermodynamic temperature, Teq,SL. Therefore, during the cooling, it is required that the fluid reaches a certain temperature, below the Teq,SL, to initiate the nucleation process. This work proposes an experimental procedure to determine consistently the highest solid-liquid thermodynamic equilibrium temperature. By means of rheometric tests, it is showed that the supercooling affects the nucleation rate and the crystals morphology, and also it is influenced by the cooling rate and the initial cooling temperature. One can conclude that the higher the supercooling the higher the number of crystals and, as consequence, the lower the crystals length. Finally, the crystals morphology and the length affect the mechanical properties of the waxy oil.

PETROBRAS / O presente trabalho contempla a análise do desenvolvimento de filme de líquido em escoamento ascendente sob a ação dos campos centrífugo e gravitacional, localizado em uma câmara ciclônica de um sistema de distribuição. No estudo foi analisada a dinâmica de um escoamento bifásico líquido-gás a partir de simulações numéricas, tridimensionais e transientes. Na modelagem numérica do escoamento foi utilizado o modelo de dois fluidos euleriano-euleriano junto com o esquema compressivo para a captura da interface líquido-gás e o modelo de turbulência SST (Shear Stress Transport). As equações de balanço foram discretizadas utilizado o método dos volumes finitos baseados em elementos finitos. As simulações numéricas foram realizadas no software comercial ANSYS-CFX 15.0. O modelo numérico desenvolvido foi validado com base em dados experimentais fornecidos pelo NUEM (Núcleo de Escoamento Multifásico). Os resultados obtidos neste trabalho possibilitaram o entendimento da dinâmica do escoamento (formação, espalhamento, estabilidade e desenvolvimento) e como a variação das velocidades superficiais do líquido e do gás (0,5; 1,0; 1,5 e 2,0 m/s para ambas as fases) influenciam na estabilidade e no comportamento do escoamento de filme de líquido ascendente. Foi elaborado ainda um modelo matemático capaz de fornecer condições de comportamento médio de um escoamento no padrão golfadas para os casos analisados. / Flows with free surface compose an important class of two-phase flows in the fluids mechanics. In various applications a thin film of liquid builds up on horizontal or vertical surfaces. This type of flow is denominated as liquid film flow. The present work contemplates the analysis of the development of liquid film in ascendant flow under the action of the centrifugal and gravitational fields, located in a cyclonic chamber of a distribution system. In the study, the dynamics of a liquid film flow were analyzed with the aid of three-dimensional and transient numerical simulations. The numerical model developed was validated based on experimental data provided by NUEM (Centre for Multiphase Flows). In the numerical modeling of the flow, were used the EulerianEulerian two fluid model coupled with the compressive discretization scheme for the capture of the liquid-gas interface and the Shear Stress Transport turbulence model. The balance equations were discretized using the finite volume based on finite elements method. Numerical simulations were supported by the commercial software ANSYS-CFX 15.0. The results obtained in this work enabled the understanding of flow dynamics (formation, spreading, stability and development) and how the superficial velocities of liquid and gas (0.5, 1.0, 1.5 and 2.0 m/s for both phases) influence on the stability and behavior of the ascendant liquid film flow. A mathematical model capable of providing conditions that represent the average behavior of a slug flow pattern was also elaborate.

CAPES / Bombas centrífugas são vastamente utilizadas em muitas aplicações industriais. O conhecimento de como estes componentes se comportam nas diversas aplicações é determinante para o desenvolvimento de sistemas de bombeamento eficientes e econômicos. A combinação de múltiplos rotores, difusores aletados e uma voluta podem apresentar características complexas no escoamento que em grande parte desviam-se da teoria de escoamento invíscido. Uma alternativa para extrair informações detalhadas do escoamento em bombas centrífugas, que forneçam subsídios para entender os fenômenos físicos envolvidos, é o uso da técnica numérica de Dinâmica dos Fluidos Computacional (CFD). Neste cenário, no presente trabalho é desenvolvido um estudo numérico do escoamento monofásico em uma bomba centrífuga de dois estágios (Imbil ITAP 65-330/2) com difusor aletado e uma voluta. As equações de conservação da massa e balanço da quantidade de movimento, que modelam o escoamento em bombas centrífugas, são resolvidas numericamente utilizando o programa comercial de CFD Ansys CFX, usando uma metodologia multi-bloco e transiente. A malha utilizada é estruturada em todo o domínio computacional. As simulações numéricas foram realizadas para o escoamento de água de um fluidos viscosificado para viscosidades variando entre 87 e 720 cP e para diferentes condições de vazão e rotação. Os resultados numéricos referentes ao desempenho global foram validados com dados experimentais obtidos por Amaral (2007) mostrando uma boa concordância, onde o erro médio foi de 6,8% para água. Foram analisados os comportamentos dos campos de escoamento, de velocidades e de pressão em função das condições de operação da bomba, no ponto de melhor eficiência e fora dela. De maneira geral os resultados obtidos atendem aos objetivos propostos e representam uma significativa contribuição ao estado da arte, devendo servir como referência para o desenvolvimento de estudos futuros. / Centrifugal pumps are vastly used in many industrial applications. Knowledge of how these components behave in several circumstances is crucial for the development of more efficient and, therefore, less expensive pumping installations. The combination of multiple impellers, vaned diffusers and a volute might introduce several complex flow characteristics that largely deviate from regular inviscid pump flow theory. Computational Fluid Dynamics can be very helpful to extract information about which physical phenomena are involved in such flows. In this sense, this work performs a numerical study of the flow in a two-stage centrifugal pump (Imbil ITAP 65-330/2) with a vaned diffuser and a volute. The flow in the pump is modeled using the software Ansys CFX, by means of a multi-block, transient rotor-stator technique, with structured grids for all pump parts. The simulations were performed using water and a mixture of water and glycerin as work fluids. Several viscosities were considered, in a range between 87 and 720 cP. Comparisons between experimental data obtained by Amaral (2007) and numerical head curves showed a good agreement, with an average deviation of 6.8% for water. The behavior of velocity, pressure and turbulence kinetic energy fields was evaluated for several operational conditions. In general, the results obtained by this work achieved the proposed goals and are a significant contribution to the understanding of the flow studied.

CAPES / Escoamentos bifásicos gás-líquido são muito comuns em aplicações industriais, especialmente nas indústrias de petróleo e gás, química e nuclear. Conforme variam as condições operacionais tais como vazões das fases, dimensão das tubulações e propriedades físicas dos fluidos, ocorrem diferentes configurações geométricas conhecidas como padrões de escoamento. No caso da produção de petróleo, o padrão encontrado com maior frequência é o de golfadas, onde se alternam regiões contínuas de líquido (pistões) e regiões de preponderância gasosa (bolhas alongadas). Em operações de produção de petróleo em águas profundas é comum encontrar cenários onde a tubulação ajusta-se ao relevo do leito marinho, apresentando leves mudanças de direção. Nesse contexto, no presente trabalho é desenvolvido um estudo experimental do escoamento bifásico líquido-gás no padrão golfadas em um duto com leve mudança de direção, composto por um trecho horizontal seguido por um trecho inclinado descendente. São produzidos escoamentos em condições controladas, dos quais são extraídos alguns parâmetros característicos. Os fluidos usados são água e ar e os experimentos são realizados no circuito de testes presente no NUEM (Núcleo de Escoamentos Multifásicos da UTFPR). A extração dos parâmetros é feita através de sensores resistivos, instalados em quatro seções da tubulação. Adicionalmente, duas câmeras de alta taxa de aquisição de imagens são utilizadas. Com os resultados obtidos, avalia-se a influência de uma leve mudança de direção nos parâmetros característicos e na transição do escoamento em golfadas para estratificado na região descendente. / Gas-liquid two-phase flow is very common in industrial applications, especially in the oil and gas, chemical, and nuclear industries. As operating conditions change such as the flow rates of the phases, the pipe diameter and physical properties of the fluids, different configurations called flow patterns take place. In the case of oil production, the most frequent pattern found is slug flow, in which continuous liquid plugs (liquid slugs) and gas-dominated regions (elongated bubbles) alternate. Offshore scenarios where the pipe lies onto the seabed with slight changes of direction are extremely common. With those scenarios and issues in mind, this work presents an experimental study of two-phase gas-liquid slug flows in a duct with a slight change of direction, represented by a horizontal section followed by a downward sloping pipe stretch. The experiments were carried out at NUEM (Núcleo de Escoamentos Multifásicos UTFPR). The flow initiated and developed under controlled conditions and their characteristic parameters were measured with resistive sensors installed at four pipe sections. Two high-speed cameras were also used. With the measured results, it was evaluated the influence of a slight direction change on the slug flow structures and on the transition between slug flow and stratified flow in the downward section.