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

Simulação computacional da interação fluido-estrutura em bombas de cavidades progressivas

The pumping through progressing cavities system has been more and more employed in the petroleum industry. This occurs because of its capacity of elevation of highly viscous oils or fluids with great concentration of sand or other solid particles. A Progressing Cavity Pump (PCP) consists, basically, of a rotor - a metallic device similar to an eccentric screw, and a stator - a steel tube internally covered by a double helix, which may be rigid or deformable/elastomeric. In general, it is submitted to a combination of well pressure with the pressure generated by the pumping process itself. In elastomeric PCPs, this combined effort compresses the stator and generates, or enlarges, the clearance existing between the rotor and the stator, thus reducing the closing effect between their cavities. Such opening of the sealing region produces what is known as fluid slip or slippage, reducing the efficiency of the PCP pumping system. Therefore, this research aims to develop a transient three-dimensional computational model that, based on single-lobe PCP kinematics, is able to simulate the fluid-structure interaction that occurs in the interior of metallic and elastomeric PCPs. The main goal is to evaluate the dynamic characteristics of PCP s efficiency based on detailed and instantaneous information of velocity, pressure and deformation fields in their interior. To reach these goals (development and use of the model), it was also necessary the development of a methodology for generation of dynamic, mobile and deformable, computational meshes representing fluid and structural regions of a PCP. This additional intermediary step has been characterized as the biggest challenge for the elaboration and running of the computational model due to the complex kinematic and critical geometry of this type of pump (different helix angles between rotor and stator as well as large length scale aspect ratios). The processes of dynamic generation of meshes and of simultaneous evaluation of the deformations suffered by the elastomer are fulfilled through subroutines written in Fortan 90 language that dynamically interact with the CFX/ANSYS fluid dynamic software. Since a structural elastic linear model is employed to evaluate elastomer deformations, it is not necessary to use any CAE package for structural analysis. However, an initial proposal for dynamic simulation using hyperelastic models through ANSYS software is also presented in this research. Validation of the results produced with the present methodology (mesh generation, flow simulation in metallic PCPs and simulation of fluid-structure interaction in elastomeric PCPs) is obtained through comparison with experimental results reported by the literature. It is expected that the development and application of such a computational model may provide better details of the dynamics of the flow within metallic and elastomeric PCPs, so that better control systems may be implemented in the artificial elevation area by PCP

Metodologia, Sistema de Codificação e Ferramenta de software para Modelagem de Móveis Modulados

On this research we investigated how new technologies can help the process of design and manufacturing of furniture in such small manufacturers in Rio Grande do Norte state. Google SketchUp, a 3D software tool, was developed in such a way that its internal structures are opened and can be accessed using SketchUp s API for Ruby and programs written in Ruby language (plugins). Using the concepts of the so-called Group Technology and the flexibility that enables adding new functionalities to this software, it was created a Methodology for Modeling of Furniture, a Coding System and a plugin for Google s tool in order to implement the Methodology developed. As resulted, the following facilities are available: the user may create and reuse the library s models over-and-over; reports of the materials manufacturing process costs are provided and, finally, detailed drawings, getting a better integration between the furniture design and manufacturing process

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

Medição e modelagem da resposta de um sensor de pig perfilométrico sob diferentes solicitações dinâmicas

Ensure the integrity of the pipeline network is an extremely important factor in the oil and gas industry. The engineering of pipelines uses sophisticated robotic inspection tools in-line known as instrumented pigs. Several relevant factors difficult the inspection of pipelines, especially in offshore field which uses pipelines with multi-diameters, radii of curvature accentuated, wall thickness of the pipe above the conventional, multi-phase flow and so on. Within this context, appeared a new instrumented Pig, called Feeler PIG, for detection and sizing of thickness loss in pipelines with internal damage. This tool was developed to overcome several limitations that other conventional instrumented pigs have during the inspection. Several factors influence the measurement errors of the pig affecting the reliability of the results. This work shows different operating conditions and provides a test rig for feeler sensors of an inspection pig under different dynamic loads. The results of measurements of the damage type of shoulder and holes in a cyclic flat surface are evaluated, as well as a mathematical model for the sensor response and their errors from the actual behavior

Obtenção de antioxidante a partir de derivados do LCC

Cashew-nut-shell-liquid (CNSL) is a phenolic oil that hás been due its their antioxirsion properties for use in fuels. The present work develops a method to the conversion of hidrogenated cardanol, that is the main component of the CNSL, in a compound with similar chacteristics to antioxidants used in products from petroleum. The antioxidants wasd obtained by exhaustive alkylation of the compound with tert-butyl chloride. After completing the optimization of several reaction steps, the product 2,4,6 tri-tert-butyl (pentadecylphenol) was obtained for the first tima. Characteeization and determination of physico-chemical properties were realized too, as well as wasd developed a study for check your application as an oxidative inhibitor by the molecular modeling. Estimation of process evalution was executed as well, where a rapid and practical computational methodology was utilizated in projects of the fine chemistry. The research showed satisfactory results and it could be concluded that the commercialization of this chemical products is feasible

Avaliação e modelagem reológica de fluido de perfuração base água

The determination of the rheology of drilling fluids is of fundamental importance to select the best composition and the best treatment to be applied in these fluids. This work presents a study of the rheological behavior of some addictives used as viscosifiers in water-based drilling fluids. The evaluated addictives were: Carboxymethylcellulose (CMC), Xanthan gum (GX), and Bentonite. The main objective was to rheologically characterize suspensions composed by these addictives, by applying mathematical models for fluid flow behavior, in order to determine the best flow equation to represent the system, as well as the model parameters. The mathematical models applied in this research were: the Bingham Model, the Ostwald de Wale Model, and the Herschel-Bulkley Model. A previous study of hydration time for each used addictive was accomplished seeking to evaluate the effect of polymer and clay hydration on rheological behavior of the fluid. The rheological characterization was made through typical rheology experiments, using a coaxial cylinder viscosimeter, where the flow curves and the thixotropic magnitude of each fluid was obtained. For each used addictive the rheological behavior as a function of temperature was also evaluated as well as fluid stability as a function of the concentration and kind of addictive used. After analyses of results, mixtures of polymer and clay were made seeking to evaluate the rheological modifications provided by the polymer incorporation in the water + bentonite system. The obtained results showed that the Ostwald de Waale model provided the best fit for fluids prepared using CMC and for fluids with Xanthan gum and Bentonite the best fit was given by the Herschel-Bulkley one

Determinação experimental e modelagem termodinâmica do ponto de turbidez de sistemas aquosos com tensoativos nonilfenolpolietoxilados

The nonionic surfactants are composed of substances whose molecules in solution, does not ionize. The solubility of these surfactants in water due to the presence of functional groups that have strong affinity for water. When these surfactants are heated is the formation of two liquid phases, evidenced by the phenomenon of turbidity. This study was aimed to determine the experimental temperature and turbidity nonilfenolpoliethoxyled subsequently perform a thermodynamic modeling, considering the models of Flory-Huggins and the empirical solid-liquid equilibrium (SLE). The method used for determining the turbidity point was the visual method (Inoue et al., 2008). The experimental methodology consisted of preparing synthetic solutions of 0,25%, 0,5%, 1%, 2%, 3%, 4%, 5%, 6%, 7%, 8%, 9%, 10%, 12,5%, 15%, 17% and 20% by weight of surfactant. The nonionic surfactants used according to their degree of ethoxylation (9.5, 10, 11, 12 and 13). During the experiments the solutions were homogenized and the bath temperature was gradually increased while the turbidity of the solution temperature was checked visually Inoue et al. (2003). These temperature data of turbidity were used to feed the models evaluated and obtain thermodynamic parameters for systems of surfactants nonilfenolpoliethoxyled. Then the models can be used in phase separation processes, facilitating the extraction of organic solvents, therefore serve as quantitative and qualitative parameters. It was observed that the solidliquid equilibrium model (ESL) was best represented the experimental data.

Determinação de dados de equilíbrio líquido-vapor a altas pressões para sistemas de hidrocarbonetos assimétricos

Crude oil is a complex liquid mixture of organic and inorganic compounds that are dominated by hydrocarbons. It is a mixture of alkanes from the simplest to more complex aromatic compounds that are present derivatives such as gasoline, diesel, alcohol, kerosene, naphtha, etc.. These derivatives are extracted from any oil, however, only with a very high quality, in other words, when the content of hydrocarbons of low molecular weight is high means that production of these compounds is feasible. The American Petroleum Institute (API) developed a classification system for the various types of oil. In Brazil, the quality of most of the oil taken from wells is very low, so it is necessary to generate new technology to develop best practices for refining in order to produce petroleum products of higher commercial value. Therefore, it is necessary to study the thermodynamic equilibrium properties of its derivative compounds of interest. This dissertation aims to determine vapor-liquid equilibrium (VLE) data for the systems Phenilcyclohexane - CO2, and Cyclohexane - Phenilcyclohexane - CO2 at high pressure and temperatures between 30 to 70oC. Furthermore, comparisons between measured VLE experimental data from this work and from the literature in relation to the Peng- Robinson molecular thermodynamic model, using a simulation program SPECS IVCSEP v5.60 and two adjustable interaction parameters, have been performed for modeling and simulation purposes. Finally, the developed apparatus for determination of phase equilibrium data at high pressures is presented

Modelagem e simulação da deposição de parafinas em escoamento turbulento

The flow assurance has become one of the topics of greatest interest in the oil industry, mainly due to production and transportation of oil in regions with extreme temperature and pressure. In these operations the wax deposition is a commonly problem in flow of paraffinic oils, causing the rising costs of the process, due to increased energy cost of pumping, decreased production, increased pressure on the line and risk of blockage of the pipeline. In order to describe the behavior of the wax deposition phenomena in turbulent flow of paraffinic oils, under different operations conditions, in this work we developed a simulator with easy interface. For that we divided de work in four steps: (i) properties estimation (physical, thermals, of transport and thermodynamics) of n-alkanes and paraffinic mixtures by using correlations; (ii) obtainment of the solubility curve and determination the wax appearance temperature, by calculating the solid-liquid equilibrium of parafinnic systems; (iii) modelling wax deposition process, comprising momentum, mass and heat transfer; (iv) development of graphic interface in MATLAB® environment for to allow the understanding of simulation in different flow conditions as well as understand the matter of the variables (inlet temperature, external temperature, wax appearance temperature, oil composition, and time) on the behavior of the deposition process. The results showed that the simulator developed, called DepoSim, is able to calculate the profile of temperature, thickness of the deposit, and the amount of wax deposited in a simple and fast way, and also with consistent results and applicable to the operation

Modelagem da transferência de massa na remoção do óleo da água produzida por flotação em coluna com uso de tensoativos de origem vegetal

The treatment of oil produced water and its implications are continually under investigation and several questions are related to this subject. In the Northeast Region Brazil, the onshore reservoirs are, in its majority, mature oil fields with high production of water. As this oil produced water has high levels of oil, it cannot be directly discarded into the environment because it represents a risk for contamination of soil, water, and groundwater, or even may cause harm to living bodies. Currently, polyelectrolytes that promote the coalescence of the oil droplets are used to remove the dispersed oil phase, enhancing the effectiveness of the flotation process. The non-biodegradability and high cost of polyelectrolytes are limiting factors for its application. On this context, it is necessary to develop studies for the search of more environmentally friendly products to apply in the flotation process. In this work it is proposed the modeling of the flotation process, in a glass column, using surfactants derived from vegetal oils to replace the polyelectrolytes, as well as to obtain a model that represents the experimental data. In addition, it was made a comparative study between the models described in the literature and the one developed in this research. The obtained results showed that the developed model presented high correlation coefficients when fitting the experimental data (R2 > 0.98), thus proving its efficiency in modeling the experimental data.

Modelagem e simulação de um separador bifásico com inversão de fases induzida

The present work has the main goal to study the modeling and simulation of a biphasic separator with induced phase inversion, the MDIF, with the utilization of the finite differences method for the resolution of the partial differencial equations which describe the transport of contaminant s mass fraction inside the equipment s settling chamber. With this aim, was developed the deterministic differential model AMADDA, wich was admensionalizated and then semidiscretizated with the method of lines. The integration of the resultant system of ordinary differential equations was realized by means of a modified algorithm of the Adam-Bashfort- Moulton method, and the sthocastic optimization routine of Basin-Hopping was used in the model s parameter estimation procedure . With the aim to establish a comparative referential for the results obtained with the model AMADDA, were used experimental data presented in previous works of the MDIF s research group. The experimental data and those obtained with the model was assessed regarding its normality by means of the Shapiro-Wilk s test, and validated against the experimental results with the Student s t test and the Kruskal-Wallis s test, depending on the result. The results showed satisfactory performance of the model AMADDA in the evaluation of the MDIF s separation efficiency, being possible to determinate that at 1% significance level the calculated results are equivalent to those determinated experimentally in the reference works

Modelagem e simulação da solubilização da parafina em tubos circulares e em regime laminar

O óleo produzido nos novos campos de petróleo está cada vez mais parafínico e viscoso, com isso, à medida que o óleo é escoado, parafinas são depositadas sobre as paredes internas do tubo, e ao longo do tempo, tendem a reduzir drasticamente a área transversal ao escoamento. Visando estudar o processo de solubilização da parafina em dutos, esse trabalho objetiva desenvolver modelos matemáticos que represente o processo, com base nos fenômenos envolvidos no mesmo tais como transferência de massa, transferência de energia e equilíbrio sólido-líquido, implementando-os em um ambiente de desenvolvimento VBA (Visual Basic) for Excel ®. O presente trabalho foi realizado em quatro etapas: i) modelagem dos fenômenos de transferência de calor e massa, ii) modelagem da rotina dos coeficientes de atividade através do modelo UNIFAC e modelagem do sistema de equilíbrio sólido-líquido; iii) modelagem matemática do processo de solubilização e cálculo da espessura da parafina ao longo do tempo; iv) implementação dos modelos em um ambiente de desenvolvimento VBA for Excel® e criação de um simulador com uma interface gráfica, para simular o processo de solubilização da parafina depositada em dutos e sua otimização. O simulador conseguiu produzir soluções bastante adequadas, mantendo continuidade das equações diferenciáveis do balanço de energia e de massa, com uma interpretação física viável, sem a presença de dissipação de oscilações nos perfis de temperatura e massa. Além disso, esse simulador visa permitir a simulação nas diversas condições de escoamento, bem como compreender a importância das variáveis (vazão, temperatura de entrada, temperatura externa, cadeia carbônica do solvente). Através dos resultados foram possíveis verificar os perfis de temperatura, fração molar e o de solubilização

Recuperação e purificação de quitosanases usando adsorção em leito expandido com streamline DEAE com modelagem e simulação usando redes neurais

Expanded Bed Adsorption (EBA) is an integrative process that combines concepts of chromatography and fluidization of solids. The many parameters involved and their synergistic effects complicate the optimization of the process. Fortunately, some mathematical tools have been developed in order to guide the investigation of the EBA system. In this work the application of experimental design, phenomenological modeling and artificial neural networks (ANN) in understanding chitosanases adsorption on ion exchange resin Streamline® DEAE have been investigated. The strain Paenibacillus ehimensis NRRL B-23118 was used for chitosanase production. EBA experiments were carried out using a column of 2.6 cm inner diameter with 30.0 cm in height that was coupled to a peristaltic pump. At the bottom of the column there was a distributor of glass beads having a height of 3.0 cm. Assays for residence time distribution (RTD) revelead a high degree of mixing, however, the Richardson-Zaki coefficients showed that the column was on the threshold of stability. Isotherm models fitted the adsorption equilibrium data in the presence of lyotropic salts. The results of experiment design indicated that the ionic strength and superficial velocity are important to the recovery and purity of chitosanases. The molecular mass of the two chitosanases were approximately 23 kDa and 52 kDa as estimated by SDS-PAGE. The phenomenological modeling was aimed to describe the operations in batch and column chromatography. The simulations were performed in Microsoft Visual Studio. The kinetic rate constant model set to kinetic curves efficiently under conditions of initial enzyme activity 0.232, 0.142 e 0.079 UA/mL. The simulated breakthrough curves showed some differences with experimental data, especially regarding the slope. Sensitivity tests of the model on the surface velocity, axial dispersion and initial concentration showed agreement with the literature. The neural network was constructed in MATLAB and Neural Network Toolbox. The cross-validation was used to improve the ability of generalization. The parameters of ANN were improved to obtain the settings 6-6 (enzyme activity) and 9-6 (total protein), as well as tansig transfer function and Levenberg-Marquardt training algorithm. The neural Carlos Eduardo de Araújo Padilha dezembro/2013 9 networks simulations, including all the steps of cycle, showed good agreement with experimental data, with a correlation coefficient of approximately 0.974. The effects of input variables on profiles of the stages of loading, washing and elution were consistent with the literature