Durante a perfuração de poços de petróleo, a determinação do perfil de temperaturas no poço é importante para tomada de decisões relativas ao processo de cimentação, para a seleção de revestimento do poço e equipamentos e, sobretudo, na identificação de zonas de influxo e perda de circulação. Neste trabalho é proposto um modelo matemático da transferência de calor em regime transitório do escoamento de fluido de perfuração em poços fraturados com perda de circulação. O poço é representado de maneira simplificada através de um cilindro anular concêntrico, cuja parede externa (interface poço-formação) apresenta uma ou mais fraturas discretas. Para a obtenção do modelo térmico é realizado um balanço de energia com foco nas trocas de calor entre a coluna de perfuração, região anular e formação rochosa. A característica principal do modelo proposto é a possibilidade de detecção da posição e número de fraturas a partir do perfil do gradiente térmico da região anular ao longo poço. Para tanto, com o código numérico, obtido via método dos volumes finitos, investiga-se a influência de parâmetros: da fratura (profundidade relativa, perda de circulação, número e distância entre fraturas), físicos (tempo de circulação) e do regime de escoamento (número de Reynolds e viscosidade dos fluidos de perfuração). As variáveis-resposta principais analisadas são a temperatura da região anular e o gradiente térmico. Como variáveis-resposta secundárias são utilizadas as evoluções térmicas da temperatura no fundo do poço e na saída da região anular. É constatado que o aumento da profundidade relativa ou número de fraturas diminui a temperatura do fundo do poço, sem causar variação significativa na temperatura de saída do anular. Para a variação da perda de circulação, o efeito na temperatura do fundo do poço é similar ao da variação do aumento da profundidade relativa da fratura, no entanto são observadas diferenças na temperatura de saída. Além disso, é verificado que, conforme se aumenta o número de fraturas distribuídas ao longo da profundidade do poço, a temperatura do poço tende ao caso de poço não fraturado. De maneira similar é evidenciada a tendência de que a diminuição na distância entre fraturas se aproxima dos resultados para um poço com uma única fratura. Finalmente, o aumento da perda de circulação facilita a detecção de fraturas devido a respectiva mudança na descontinuidade do perfil do gradiente térmico da região anular. / During oil drilling operations, the wellbore temperature profile is used when selecting well casing materials, making cementation related decisions, and, most importantly, to identify loss zones. In this work, a transient heat transfer mathematical model for a fractured wellbore is proposed. The well has its geometry simplified to a concentric annular cylinder which has one or more discrete fracture in its external wall (wellformation interface). In order to obtain the thermal model an energy balance is used, focusing the heat transfer between the pipe, the annular region and the formation. The key characteristic of the model is the fracture detection through thermal gradient graphical analysis. The thermal gradient is an output of the solution of the discretized energy equation in the domains, obtained through the finite volume method. The following parameters are investigated: fracture depth, fracture number, fracture interference, loss circulation, circulation time, Reynolds number and drilling fluid viscosity. The analysis is done through the analysis of the annular region temperature profile and its gradient, along with the thermal evolution of both the bottomhole and outlet temperatures. It is verified that increasing the fracture relative depth or number decreases the bottomhole temperature, while having no significant impact in the outlet temperature. The same bottomhole temperature effect is noted when increasing loss rate, however outlet temperature changes are observed. In a similar way, when decreasing the distance between fractures, the temperature profile in the annular region trends to a wellbore with a single fracture. Finally, increasing loss rate favors fracture detection, since the discontinuity in the annular region thermal gradient profile is intensified.

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.

Durante a perfuração de poços de petróleo, a determinação do perfil de temperaturas no poço é importante para tomada de decisões relativas ao processo de cimentação, para a seleção de revestimento do poço e equipamentos e, sobretudo, na identificação de zonas de influxo e perda de circulação. Neste trabalho é proposto um modelo matemático da transferência de calor em regime transitório do escoamento de fluido de perfuração em poços fraturados com perda de circulação. O poço é representado de maneira simplificada através de um cilindro anular concêntrico, cuja parede externa (interface poço-formação) apresenta uma ou mais fraturas discretas. Para a obtenção do modelo térmico é realizado um balanço de energia com foco nas trocas de calor entre a coluna de perfuração, região anular e formação rochosa. A característica principal do modelo proposto é a possibilidade de detecção da posição e número de fraturas a partir do perfil do gradiente térmico da região anular ao longo poço. Para tanto, com o código numérico, obtido via método dos volumes finitos, investiga-se a influência de parâmetros: da fratura (profundidade relativa, perda de circulação, número e distância entre fraturas), físicos (tempo de circulação) e do regime de escoamento (número de Reynolds e viscosidade dos fluidos de perfuração). As variáveis-resposta principais analisadas são a temperatura da região anular e o gradiente térmico. Como variáveis-resposta secundárias são utilizadas as evoluções térmicas da temperatura no fundo do poço e na saída da região anular. É constatado que o aumento da profundidade relativa ou número de fraturas diminui a temperatura do fundo do poço, sem causar variação significativa na temperatura de saída do anular. Para a variação da perda de circulação, o efeito na temperatura do fundo do poço é similar ao da variação do aumento da profundidade relativa da fratura, no entanto são observadas diferenças na temperatura de saída. Além disso, é verificado que, conforme se aumenta o número de fraturas distribuídas ao longo da profundidade do poço, a temperatura do poço tende ao caso de poço não fraturado. De maneira similar é evidenciada a tendência de que a diminuição na distância entre fraturas se aproxima dos resultados para um poço com uma única fratura. Finalmente, o aumento da perda de circulação facilita a detecção de fraturas devido a respectiva mudança na descontinuidade do perfil do gradiente térmico da região anular. / During oil drilling operations, the wellbore temperature profile is used when selecting well casing materials, making cementation related decisions, and, most importantly, to identify loss zones. In this work, a transient heat transfer mathematical model for a fractured wellbore is proposed. The well has its geometry simplified to a concentric annular cylinder which has one or more discrete fracture in its external wall (wellformation interface). In order to obtain the thermal model an energy balance is used, focusing the heat transfer between the pipe, the annular region and the formation. The key characteristic of the model is the fracture detection through thermal gradient graphical analysis. The thermal gradient is an output of the solution of the discretized energy equation in the domains, obtained through the finite volume method. The following parameters are investigated: fracture depth, fracture number, fracture interference, loss circulation, circulation time, Reynolds number and drilling fluid viscosity. The analysis is done through the analysis of the annular region temperature profile and its gradient, along with the thermal evolution of both the bottomhole and outlet temperatures. It is verified that increasing the fracture relative depth or number decreases the bottomhole temperature, while having no significant impact in the outlet temperature. The same bottomhole temperature effect is noted when increasing loss rate, however outlet temperature changes are observed. In a similar way, when decreasing the distance between fractures, the temperature profile in the annular region trends to a wellbore with a single fracture. Finally, increasing loss rate favors fracture detection, since the discontinuity in the annular region thermal gradient profile is intensified.

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.

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.

Os métodos numéricos convencionais, baseados em malhas, têm sido amplamente aplicados na resolução de problemas da Dinâmica dos Fluidos Computacional. Entretanto, em problemas de escoamento de fluidos que envolvem superfícies livres, grandes explosões, grandes deformações, descontinuidades, ondas de choque etc., estes métodos podem apresentar algumas dificuldades práticas quando da resolução destes problemas. Como uma alternativa viável, existem os métodos de partículas livre de malhas. Neste trabalho é feita uma introdução ao método Lagrangeano de partículas, livre de malhas, Smoothed Particle Hydrodynamics (SPH) voltado para a simulação numérica de escoamentos de fluidos newtonianos compressíveis e quase-incompressíveis. Dois códigos numéricos foram desenvolvidos, uma versão serial e outra em paralelo, empregando a linguagem de programação C/C++ e a Compute Unified Device Architecture (CUDA), que possibilita o processamento em paralelo empregando os núcleos das Graphics Processing Units (GPUs) das placas de vídeo da NVIDIA Corporation. Os resultados numéricos foram validados e a eficiência computacional avaliada considerandose a resolução dos problemas unidimensionais Shock Tube e Blast Wave e bidimensional da Cavidade (Shear Driven Cavity Problem). / The conventional mesh-based numerical methods have been widely applied to solving problems in Computational Fluid Dynamics. However, in problems involving fluid flow free surfaces, large explosions, large deformations, discontinuities, shock waves etc. these methods suffer from some inherent difficulties which limit their applications to solving these problems. Meshfree particle methods have emerged as an alternative to the conventional grid-based methods. This work introduces the Smoothed Particle Hydrodynamics (SPH), a meshfree Lagrangian particle method to solve compressible flows. Two numerical codes have been developed, serial and parallel versions, using the Programming Language C/C++ and Compute Unified Device Architecture (CUDA). CUDA is NVIDIAs parallel computing architecture that enables dramatic increasing in computing performance by harnessing the power of the Graphics Processing Units (GPUs). The numerical results were validated and the speedup evaluated for the Shock Tube and Blast Wave one-dimensional problems and Shear Driven Cavity Problem.