Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems

Simulação da transferência de calor em amostras aquecidas por plasma

The processing of materials through plasma has been growing enough in the last times in several technological applications, more specifically in surfaces treatment. That growth is due, mainly, to the great applicability of plasmas as energy source, where it assumes behavior thermal, chemical and/or physical. On the other hand, the multiplicity of simultaneous physical effects (thermal, chemical and physical interactions) present in plasmas increases the complexity for understanding their interaction with solids. In that sense, as an initial step for the development of that subject, the present work treats of the computational simulation of the heating and cooling processes of steel and copper samples immersed in a plasma atmosphere, by considering two experimental geometric configurations: hollow and plane cathode. In order to reach such goal, three computational models were developed in Fortran 90 language: an one-dimensional transient model (1D, t), a two-dimensional transient model (2D, t) and a two-dimensional transient model (2D, t) which take into account the presence of a sample holder in the experimental assembly. The models were developed based on the finite volume method and, for the two-dimensional configurations, the effect of hollow cathode on the sample was considered as a lateral external heat source. The main results obtained with the three computational models, as temperature distribution and thermal gradients in the samples and in the holder, were compared with those developed by the Laboratory of Plasma, LabPlasma/UFRN, and with experiments available in the literature. The behavior showed indicates the validity of the developed codes and illustrate the need of the use of such computational tool in that process type, due to the great easiness of obtaining thermal information of interest

Simulação computacional do escoamento em bombas de cavidades progressivas

The use of Progressing Cavity Pumps (PCPs) in artificial lift applications in low deep wells is becoming more common in the oil industry, mainly, due to its ability to pump heavy oils, produce oil with large concentrations of sand, besides present high efficiency when compared to other artificial lift methods. Although this system has been widely used as an oil lift method, few investigations about its hydrodynamic behavior are presented, either experimental or numeric. Therefore, in order to increase the knowledge about the BCP operational behavior, this work presents a novel computational model for the 3-D transient flow in progressing cavity pumps, which includes the relative motion between rotor and stator, using an element based finite volume method. The model developed is able to accurately predict the volumetric efficiency and viscous looses as well as to provide detailed information of pressure and velocity fields inside the pump. In order to predict PCP performance for low viscosity fluids, advanced turbulence models were used to treat, accurately, the turbulent effects on the flow, which allowed for obtaining results consistent with experimental values encountered in literature. In addition to the 3D computational model, a simplified model was developed, based on mass balance within cavities and on simplification on the momentum equations for fully developed flow along the seal region between cavities. This simplified model, based on previous approaches encountered in literature, has the ability to predict flow rate for a given differential pressure, presenting exactness and low CPU requirements, becoming an engineering tool for quick calculations and providing adequate results, almost real-time time. The results presented in this work consider a rigid stator PCP and the models developed were validated against experimental results from open literature. The results for the 3-D model showed to be sensitive to the mesh size, such that a numerical mesh refinement study is also presented. Regarding to the simplified model, some improvements were introduced in the calculation of the friction factor, allowing the application fo the model for low viscosity fluids, which was unsuccessful in models using similar approaches, presented in previous works

Estudo numérico de escoamento turbulento em padrão anular gás-líquido em dutos verticais

Multiphase flows in ducts can adopt several morphologies depending on the mass fluxes and the fluids properties. Annular flow is one of the most frequently encountered flow patterns in industrial applications. For gas liquid systems, it consists of a liquid film flowing adjacent to the wall and a gas core flowing in the center of the duct. This work presents a numerical study of this flow pattern in gas liquid systems in vertical ducts. For this, a solution algorithm was developed and implemented in FORTRAN 90 to numerically solve the governing transport equations. The mass and momentum conservation equations are solved simultaneously from the wall to the center of the duct, using the Finite Volumes Technique. Momentum conservation in the gas liquid interface is enforced using an equivalent effective viscosity, which also allows for the solution of both velocity fields in a single system of equations. In this way, the velocity distributions across the gas core and the liquid film are obtained iteratively, together with the global pressure gradient and the liquid film thickness. Convergence criteria are based upon satisfaction of mass balance within the liquid film and the gas core. For system closure, two different approaches are presented for the calculation of the radial turbulent viscosity distribution within the liquid film and the gas core. The first one combines a k- Ɛ one-equation model and a low Reynolds k-Ɛ model. The second one uses a low Reynolds k- Ɛ model to compute the eddy viscosity profile from the center of the duct right to the wall. Appropriate interfacial values for k e Ɛ are proposed, based on concepts and ideas previously used, with success, in stratified gas liquid flow. The proposed approaches are compared with an algebraic model found in the literature, specifically devised for annular gas liquid flow, using available experimental results. This also serves as a validation of the solution algorithm

Modelagem numérica de escoamento em padrão anular com transferência de calor e massa

Annular flow is the prevailing pattern in transport and energy conversion systems and therefore, one of the most important patterns in multiphase flow in ducts. The correct prediction of the pressure gradient and heat transfer coefficient is essential for optimizing the system s capacity. The objective of this work is to develop and implement a numerical algorithm capable of predicting hydrodynamic and thermal characteristics for upflow, vertical, annular flow. The numerical algorithm is then complemented with the physical modeling of phenomena that occurs in this flow pattern. These are, turbulence, entrainment and deposition and phase change. For the development of the numerical model, axial diffusion of heat and momentum is neglected. In this way the time-averaged equations are solved in their parabolic form obtaining the velocity and temperature profiles for each axial step at a time, together with the global parameters, namely, pressure gradient, mean film thickness and heat transfer coefficient, as well as their variation in the axial direction. The model is validated for the following conditions: fully-developed laminar flow with no entrainment; fully developed laminar flow with heat transfer, fully-developed turbulent flow with entrained drops, developing turbulent annular flow with entrained drops, and turbulent flow with heat transfer and phase change

Modelagem dinâmica do escoamento de um sistema de elevação por plunger lift

A critical problem in mature gas wells is the liquid loading. As the reservoir pressure decreases, gas superficial velocities decreases and the drag exerted on the liquid phase may become insufficient to bring all the liquid to the surface. Liquid starts to drain downward, flooding the well and increasing the backpressure which decreases the gas superficial velocity and so on. A popular method to remedy this problem is the Plunger Lift. This method consists of dropping the "plunger"to the bottom of the tubing well with the main production valve closed. When the plunger reaches the well bottom the production valve is opened and the plunger carry the liquid to the surface. However, models presented in literature for predicting the behavior in plunger lift are simplistic, in many cases static (not considering the transient effects). Therefore work presents the development and validation of a numerical algorithm to solve one-dimensional compressible in gas wells using the Finite Volume Method and PRIME techniques for treating coupling of pressure and velocity fields. The code will be then used to develop a dynamic model for the plunger lift which includes the transient compressible flow within the well

Simulação do escoamento monofásico em um estágio de uma bomba centrífuga utilizando técnicas de fluidodinâmica computacional

Oil production and exploration techniques have evolved in the last decades in order to increase fluid flows and optimize how the required equipment are used. The base functioning of Electric Submersible Pumping (ESP) lift method is the use of an electric downhole motor to move a centrifugal pump and transport the fluids to the surface. The Electric Submersible Pumping is an option that has been gaining ground among the methods of Artificial Lift due to the ability to handle a large flow of liquid in onshore and offshore environments. The performance of a well equipped with ESP systems is intrinsically related to the centrifugal pump operation. It is the pump that has the function to turn the motor power into Head. In this present work, a computer model to analyze the three-dimensional flow in a centrifugal pump used in Electric Submersible Pumping has been developed. Through the commercial program, ANSYS® CFX®, initially using water as fluid flow, the geometry and simulation parameters have been defined in order to obtain an approximation of what occurs inside the channels of the impeller and diffuser pump in terms of flow. Three different geometry conditions were initially tested to determine which is most suitable to solving the problem. After choosing the most appropriate geometry, three mesh conditions were analyzed and the obtained values were compared to the experimental characteristic curve of Head provided by the manufacturer. The results have approached the experimental curve, the simulation time and the model convergence were satisfactory if it is considered that the studied problem involves numerical analysis. After the tests with water, oil was used in the simulations. The results were compared to a methodology used in the petroleum industry to correct viscosity. In general, for models with water and oil, the results with single-phase fluids were coherent with the experimental curves and, through three-dimensional computer models, they are a preliminary evaluation for the analysis of the two-phase flow inside the channels of centrifugal pump used in ESP systems

Análise de estruturas de superfície delgadas axissimétricas via método dos elementos finitos com utilização de elemento retilíneo

The present work deals with the linear analysis of bi-dimensional axisymmetric structures, through development and implementation of a Finite Element Method code. The structures are initially studied alone and afterwards compatibilized into coupled structures, that is, assemblages, including tanks and pressure vessels. Examples are analysed and, in order to prove accuracy, the results were compared with those furnished by the analytical solutions

Análise não-linear de treliças pelo método dos elementos finitos posicional

This work presents the positional nonlinear geometric formulation for trusses using different strain measures. The positional formulation presents an alternative approach for nonlinear problems. This formulation considers nodal positions as variables of the nonlinear system instead of displacements (widely found in literature). The work also describes the arc-length method used for tracing equilibrium paths with snap-through and snap-back. Numerical applications for trusses already established in the literature and comparisons with other studies are provided to prove the accuracy of the proposed formulation

Desempenho e efeitos agudos da espirometria de incentivo volume-orientada sobre os volumes pulmonares em pacientes após AVC e sujeitos saudáveis

O Acidente Vascular Cerebral (AVC) é uma síndrome clínica caracterizada por uma perturbação focal da função cerebral. Geralmente ocasiona quadro de disfunção motora acompanhada de prejuízo da função respiratória. Tendo em vista tal quadro e suas possíveis repercussões, a fisioterapia respiratória tem sido amplamente requisitada com o intuito reverter ou minimizar as complicações. Dentre os recursos utilizados para este fim, os espirômetros de incentivo são utilizados com o objetivo de restaurar os volumes pulmonares, modificando o padrão respiratório e de ventilação pulmonar, prevenindo a incidência das complicações pulmonares. O objetivo do presente estudo foi comparar o desempenho na realização da espirometria de incentivo volume-orientada (EI) e seus efeitos agudos sobre os volumes pulmonares em pacientes após AVC e sujeitos saudáveis. Foram selecionados 40 voluntários, de ambos os gêneros, divididos em grupo experimental (GE), composto por 20 pacientes após AVC e grupo controle (GC) composto por 20 sujeitos saudáveis, pareados quanto à idade, gênero e Índice de Massa Corpórea (IMC). A coleta dos dados foi realizada em duas etapas: (1) Avaliação cognitiva e neurofuncional (Mini Exame do Estado Mental, National Institute Stroke Scale, Medida de Independência Funcional, teste de desempenho da aprendizagem) (2) Avaliação Respiratória (espirometria, manovacuometria e cinemática tóraco-abdominal, através da Pletismografia Opto-eletrônica). A análise estatística foi realizada através do software Graphpad Prism 4.0, em que foram utilizados os testes t Student e ANOVA two-way para comparação intergrupos e adotado nível de significância de 5%. Os resultados mostraram que os pacientes apresentam desempenho inferior na aprendizagem da espirometria de incentivo, com uma média de erros maior 2,95 ± 1,39, quando comparados aos sujeitos saudáveis, 1,15 ± 0,98. Em relação à prática observacional utilizada não foi encontrada diferença entre a aprendizagem através do vídeo ou do terapeuta. Em relação aos efeitos agudos da espirometria de incentivo volume-orientada, os pacientes após AVC apresentaram valores de volume corrente 24,7%, 18% e 14,7% inferiores quando comparados aos sujeitos saudáveis nos momentos pré-EI, EI e pós-EI, porém a espirometria de incentivo induziu incrementos de volume similares em ambos os grupos estudados, com 75, 3% para os pacientes e 73,3% para os sujeitos saudáveis. A espirometria de incentivo promove ganhos significativos no volume corrente da parede torácica, tanto em pacientes acometidos por AVC como em sujeitos saudáveis, no entanto o desempenho da aprendizagem é inferior para os pacientes, para ambas as práticas com vídeo ou terapeuta

Efeitos da respiração em freno-labial sobre os volumes pulmonares e o padrão de hiperinsuflação dinâmica em pacientes com asma

Objectives: To evaluate how to develop dynamic hyperinflation (DH) during exercise, the influence of pursed-lip breathing in (PLB) on breathing pattern and operating volume in patients with asthma. Methods: We studied 12 asthmatic patients in three moments: (1) anthropometry and spirometry, (2) submaximal incremental cycle ergometer test in spontaneous breathing and (3), submaximal incremental test on a cycle ergometer with PLB using the Opto-electronic plethysmography. Results: Evaluating the end-expiratory lung volume (EEV) during submaximal incremental test in spontaneous breathing, patients were divided into euvolume and hyperinflated. The RFL has increased significantly, the variation of the EEV group euvolume (1.4L) and decreased in group hyperinflated (0.272L). In group volume observed a significant increase of 140% in Vt at baseline, before exercise, comparing the RFL and spontaneous breathing. Hyperinflated group was observed that the RFL induced significant increases of Vt at all times of the test incremental baseline, 50%, 100% load and 66% recovery, 250%, 61.5% and 66% respectively. Respiratory rate decreased significantly with PLB at all times of the submaximal incremental test in the group euvolume. The speed of shortening of inspiratory muscles (VtRcp/Ti) in the hyperinflated increased from 1.6 ± 0.8L/s vs. 2.55 ± 0.9L/s, whereas in the RFL euvolume group ranged from 0.72 ± 0.31L/s vs. 0.65 ± 0.2L/s. The velocity of shortening of the expiratory muscles (VtAb/Te) showed similarity in response to RFL. In group hyperinflated varied vs. 0.89 ± 0.47 vs. 0.80 ± 0.36 and ± 1.17 ± 1L vs. 0.78 ± 0.6 for group euvolume. Conclusion: Different behavior in relation to EEV in patients with moderate asthma were observed, the HD and decreased EEV in response to exercise. The breathing pattern was modulated by both RFL performance as at home, making it more efficient

Efeitos agudos de diferentes intensidades de pressão expiratória positiva sobre os volumes pulmonares em pacientes com acidente vascular encefálico

observar os efeitos agudos de diferentes intensidades de Pressão Expiratória Positiva (PEP) sobre a cinemática do complexo toracoabdominal de pacientes acometidos por Acidente Vascular Encefálico (AVE). Métodos: Foram selecionados 21 indivíduos com AVE e 16 indivíduos saudáveis pareados por idade sexo e IMC para grupo controle. Avaliamos função pulmonar, pressões dos músculos respiratórios e os volumes pulmonares por meio da Pletismografia Optoeletrônica durante três diferentes intensidades de PEP 10, 15 e 20 cmH2O. Resultados: o efeito da PEP no volume corrente (VC) do grupo AVE em relação ao grupo controle foi diferente. Enquanto o grupo controle aumentou o VC em relação a respiração tranquila em 343%, 395,2% e 431,8% nas PEP10, PEP15 e PEP20 cmH2O o grupo AVE aumento 186%, 218.8% e 209.5% (p < 0.0001). A PEP também influenciou de forma diferente em relação ao Tempo inspiratório com intensidades diferentes no grupo controle e AVE (p < 0.0001). No ciclo de trabalho foi observado um aumento no grupo controle nas PEP10 (p < 0.001) e PEP15 (p < 0.05) e no grupo AVE foi observada uma redução PEP20 (p < 0.01) quando comparada com a respiração tranquila. Os volumes operacionais do grupo AVE foi observado aumento do volume inspiratório final da parede torácica (Vifpt) e do Volume expiratório final da parede torácica (Vefpt) diferente do grupo controle que gerou aumento do Vifpt acompanhado de diminuição do Vefpt durante as três intensidades de PEP. Conclusão: A hiperinsuflação observada no grupo AVE demonstra que essa terapêutica deve ser utilizada com cautela especialmente nas intensidades maiores que 10 cmH2O para essa população

Avaliação dos efeitos agudos de diferentes niveis de PEP sobre o padrão respiratório e volumes operacionais de pacientes com Parkinson

Objetivo: Determinar os efeitos agudos de diferentes intensidades de pressão expiratória positiva sobre o padrão respiratório e volumes operacionais de pacientes com doença de Parkinson. Métodos: Foram eleitos para o estudo 23 pacientes em estágios II ou III da doença, estando na condição ON , sendo 8 excluídos, e 15 controles saudáveis. Os indivíduos foram submetidos a uma avaliação inicial, constando de dados gerais, avaliação antropométrica, função pulmonar e força muscular respiratória. Após esta etapa inicial, realizou-se a avaliação dos volumes pulmonares com pletismografia optoeletrônica associado a utilização de pressão expiratória positiva (PEP) em três intensidades de pressão positiva, 10 cmH2O, 15 cmH2O e 20 cmH2O em ordem aleatória. Resultados: O grupo Parkinson apresentou valores espirométricos e de força muscular respiratória significativamente menores que o grupo controle (p<0.01). Houve diferença nos valores de volume corrente (Vt) do grupo Parkinson em relação ao grupo controle na respiração tranqüila (p<0.001) e aumento no Vt do grupo Parkinson com uso das três intensidades de PEP (p<0.001), sem diferenças estatisticamente significativas entre as intensidades e com diferença na distribuição compartimental do Vt entre os grupos (p<0.001). O fluxo inspiratório médio e o volume minuto do grupo Parkinson foram menores em relação ao grupo controle na respiração tranqüila (p<0.001) e aumentaram com o uso das três intensidades de PEP (p<0.001), sem diferenças estatisticamente significativas entre as intensidades. Houve diferença na variação dos volumes operacionais entre os grupos em todas as intensidades de PEP (p < 0.001). O volume pulmonar expiratório final não aumentou no grupo Parkinson com uso de PEP. O volume pulmonar inspiratório final aumentou no grupo Parkinson em todas intensidades de PEP (p < 0.001), não havendo diferenças entre as intensidades. Conclusão: Os pacientes com doença de Parkinson apresentam alterações em diferentes componentes do padrão respiratório e a terapia com PEP determina modificações clinicamente importantes nos volumes pulmonares, considerando a intensidade de 10 cmH2O suficientes para este objetivo terapêutico

Efeitos das técnicas breath-stacking e espirometria de incentivo nos volumes da caixa torácica em pacientes com Doença de Parkinson

Coordenação de Aperfeiçoamento de Pessoal de Nível Superior