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

Numerical simulation of the liquid phase in SnO2 thin film deposition by sol-gel-dip-coating

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

Numerical simulation of drop impact and jet buckling problems using the eXtended Pom-Pom model

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

A continuously differentiable upwinding scheme for the simulation of fluid flow problems

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

NUMERICAL SIMULATION of THE BEHAVIOR of THE BEHAVIOR of REINFORCED CONCRETE BEAMS and COLUMNS BY FEM CAASTEM 2000 PROGRAM

The objective of this paper is the numerical study of the behavior of reinforced concrete beams and columns by non-linear numerical simulations. The numerical analysis is based on the finite element method implemented in CASTEM 2000. This program uses the constitutive elastoplastic perfect model for the steel, the Drucker-Prager model for the concrete and the Newton-Raphson for the solution of non-linear systems. This work concentrates on the determination of equilibrium curves to the beams and force-strain curves to the columns. The numeric responses are confronted with experimental results found in the literature in order to check there liability of the numerical analyses.

Stress Analysis in Simulation Models With or Without Implant Threads Representation

Purpose: This study aimed to evaluate the influence of implants with or without threads representation on the outcome of a two-dimensional finite element (FE) analysis. Materials and Methods: Two-dimensional FE models that reproduced a frontal section of edentulous mandibular posterior bone were constructed using a standard crown/implant/screw system representation. To evaluate the effect of implant threads, two models were created: a model in which the implant threads were accurately simulated (precise model) and a model in which implants with a smooth surface (press-fit implant) were used (simplified model). An evaluation was performed on ANSYS software, in which a load of 133 N was applied at a 30-degree angulation and 2 mm off-axis from the long axis of the implant on the models, The Von Mises stresses were measured. Results: The precise model (1.45 MPa) showed higher maximum stress values than the simplified model (1.2 MPa). Whereas in the cortical bone, the stress values differed by about 36% (292.95 MPa for the precise model and 401.14 MPa for the simplified model), in trabecular bone (19.35 MPa and 20.35 MPa, respectively), the stress distribution and stress values were similar. Stress concentrations occurred around the implant neck and the implant apex. Conclusions: Considering implant and cortical bone analysis, remarkable differences in stress values were found between the models. Although the models showed different absolute stress values, the stress distribution was similar. INT J ORAL MAXILLOFAC IMPLANTS 2009;24:1040-1044

Dynamics of the deconfinement transition of quarks and gluons in a finite volume

We employ finite elements methods for the approximation of solutions of the Ginzburg-Landau equations describing the deconfinement transition in quantum chromodynamics. These methods seem appropriate for situations where the deconfining transition occurs over a finite volume as in relativistic heavy ion collisions. where in addition expansion of the system and flow of matter are important. Simulation results employing finite elements are presented for a Ginzburg-Landau equation based on a model free energy describing the deconfining transition in pure gauge SU(2) theory. Results for finite and infinite system are compared. (C) 2009 Elsevier B.V. All rights reserved.

Statistics and kinetics of single-molecule electron transfer dynamics in complex environments: A simulation model study

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

QCD tests with an infrared finite gluon propagator and coupling constant

Recent progress in the solution of Schwinger-Dyson equations, as well as lattice simulation of pure glue QCD, indicate that the gluon propagator and coupling constant are infrared finite. Such non-perturbative information can be introduced in the QCD perturbative expansion in the scheme named Dynamical Perturbation Theory. We exemplify this procedure with the calculation of some two-body non-leptonic annihilation B meson decays, which show agreement with the experimental data in the case of a gluon propagator characterized by a dynamical gluon mass of 500MeV, compatible with the value found in several processes computed with this method. We give a. preliminary account of the application of this procedure at the loop level in the case of the Bjorken sum rule.

Numerical simulation of Ginzburg-Landau-Langevin equations

This work is concerned with non-equilibrium phenomena, with focus on the numerical simulation of the relaxation of non-conserved order parameters described by stochastic kinetic equations known as Ginzburg-Landau-Langevin (GLL) equations. We propose methods for solving numerically these type of equations, with additive and multiplicative noises. Illustrative applications of the methods are presented for different GLL equations, with emphasis on equations incorporating memory effects.

Diffraction and an infrared finite gluon propagator

We discuss some phenomenological applications of an infrared finite gluon propagator characterized by a dynamically generated gluon mass. In particular we compute the effect of the dynamical gluon mass on pp and ${\bar{p}}p$ diffractive scattering. We also show how the data on gammap photoproduction and hadronic gg reactions can be derived from the pp and ${\bar{p}}p$ forward scattering amplitudes by assuming vector meson dominance and the additive quark model.

A deterministic simulation study of embryo marker-assisted selection for age at first calving in Nellore (Bos indicus) beef cattle

Fundação de Amparo à Pesquisa do Estado de São Paulo (FAPESP)

Study of stress distribution with finite element method in a root with prefabricated threaded split shaft post - Part I

A finite element analysis was carried out to study the role of prefabricated threaded split shaft post (Flexi-Post) on dentinal stress in pulpless tooth. Three dimensional plane strain model of mesio-distal section of a human maxillary central incisor without restoration was analysed with the MSC/NASTRAN (MacNeal/ Schwendler) general purpose finite analysis program was executed on a microcomputer. The model as discretized into 48.954 axisymmetric finite elements defined by 10.355 nodes. Each element was assigned unique elastic properties to represent the materials modeled. Homogeneity, isotropy and linear elasticity were assume for all material. A simulation of static load of 100N was applied to the incisal edge of the post; vertical. Maximal principal stresses and von Mises equivalent stress were calculated. Using the element analysis model employed in this study, the following can be concluded concerning threaded split shaft post (Flexi-Post): Maximum principal stresses in dentin were located at cervical place and at the post apex. The apical threads of the post not redirecting stresses away from the root.

Simulation of free surface flows in a distributed memory environment

A parallel technique, for a distributed memory machine, based on domain decomposition for solving the Navier-Stokes equations in cartesian and cylindrical coordinates in two dimensions with free surfaces is described. It is based on the code by Tome and McKee (J. Comp. Phys. 110 (1994) 171-186) and Tome (Ph.D. Thesis, University of Strathclyde, Glasgow, 1993) which in turn is based on the SMAC method by Amsden and Harlow (Report LA-4370, Los Alamos Scientific Laboratory, 1971), which solves the Navier-Stokes equations in three steps: the momentum and Poisson equations and particle movement, These equations are discretized by explicit and 5-point finite differences. The parallelization is performed by splitting the computation domain into vertical panels and assigning each of these panels to a processor. All the computation can then be performed using nearest neighbour communication. Test runs comparing the performance of the parallel with the serial code, and a discussion of the load balancing question are presented. PVM is used for communication between processes. (C) 1999 Elsevier B.V. B.V. All rights reserved.

Dynamic modeling and simulation of a flexible robotic manipulator

This work focuses on the dynamic modeling of a flexible robotic manipulator with two flexible links and two revolute joints, which rotates in the horizontal plane. The dynamic equations are derived using the Newton-Euler formulation and the finite element method, based on elementary beam theory. Computer simulation results are presented to illustrate this study. The dynamic model becomes necessary for use in future design and control applications.