Implementação da rotina de unfolding para determinação de distribuição de tamanho de grãos esféricos via distribuição de interceptos lineares e de área de seção

Na unfolding method of linear intercept distributions and secction área distribution was implemented for structures with spherical grains. Although the unfolding routine depends on the grain shape, structures with spheroidal grains can also be treated by this routine. Grains of non-spheroidal shape can be treated only as approximation. A software was developed with two parts. The first part calculates the probability matrix. The second part uses this matrix and minimizes the chi-square. The results are presented with any number of size classes as required. The probability matrix was determined by means of the linear intercept and section area distributions created by computer simulation. Using curve fittings the probability matrix for spheres of any sizes could be determined. Two kinds of tests were carried out to prove the efficiency of the Technique. The theoretical tests represent ideal cases. The software was able to exactly find the proposed grain size distribution. In the second test, a structure was simulated in computer and images of its slices were used to produce the corresponding linear intercept the section area distributions. These distributions were then unfolded. This test simulates better reality. The results show deviations from the real size distribution. This deviations are caused by statistic fluctuation. The unfolding of the linear intercept distribution works perfectly, but the unfolding of section area distribution does not work due to a failure in the chi-square minimization. The minimization method uses a matrix inversion routine. The matrix generated by this procedure cannot be inverted. Other minimization method must be used

Otimização e análise mecânica de pastas geopoliméricas para uso em poços sujeitos à injeção cíclica de vapor

Oil wells subjected to cyclic steam injection present important challenges for the development of well cementing systems, mainly due to tensile stresses caused by thermal gradients during its useful life. Cement sheath failures in wells using conventional high compressive strength systems lead to the use of cement systems that are more flexible and/or ductile, with emphasis on Portland cement systems with latex addition. Recent research efforts have presented geopolymeric systems as alternatives. These cementing systems are based on alkaline activation of amorphous aluminosilicates such as metakaolin or fly ash and display advantageous properties such as high compressive strength, fast setting and thermal stability. Basic geopolymeric formulations can be found in the literature, which meet basic oil industry specifications such as rheology, compressive strength and thickening time. In this work, new geopolymeric formulations were developed, based on metakaolin, potassium silicate, potassium hydroxide, silica fume and mineral fiber, using the state of the art in chemical composition, mixture modeling and additivation to optimize the most relevant properties for oil well cementing. Starting from molar ratios considered ideal in the literature (SiO2/Al2O3 = 3.8 e K2O/Al2O3 = 1.0), a study of dry mixtures was performed,based on the compressive packing model, resulting in an optimal volume of 6% for the added solid material. This material (silica fume and mineral fiber) works both as an additional silica source (in the case of silica fume) and as mechanical reinforcement, especially in the case of mineral fiber, which incremented the tensile strength. The first triaxial mechanical study of this class of materials was performed. For comparison, a mechanical study of conventional latex-based cementing systems was also carried out. Regardless of differences in the failure mode (brittle for geopolymers, ductile for latex-based systems), the superior uniaxial compressive strength (37 MPa for the geopolymeric slurry P5 versus 18 MPa for the conventional slurry P2), similar triaxial behavior (friction angle 21° for P5 and P2) and lower stifness (in the elastic region 5.1 GPa for P5 versus 6.8 GPa for P2) of the geopolymeric systems allowed them to withstand a similar amount of mechanical energy (155 kJ/m3 for P5 versus 208 kJ/m3 for P2), noting that geopolymers work in the elastic regime, without the microcracking present in the case of latex-based systems. Therefore, the geopolymers studied on this work must be designed for application in the elastic region to avoid brittle failure. Finally, the tensile strength of geopolymers is originally poor (1.3 MPa for the geopolymeric slurry P3) due to its brittle structure. However, after additivation with mineral fiber, the tensile strength became equivalent to that of latex-based systems (2.3 MPa for P5 and 2.1 MPa for P2). The technical viability of conventional and proposed formulations was evaluated for the whole well life, including stresses due to cyclic steam injection. This analysis was performed using finite element-based simulation software. It was verified that conventional slurries are viable up to 204ºF (400ºC) and geopolymeric slurries are viable above 500ºF (260ºC)

Análise gráfica de estruturas porosas sobre a ótica da estereologia

In this work we developed a computer simulation program for physics porous structures based on programming language C + + using a Geforce 9600 GT with the PhysX chip, originally developed for video games. With this tool, the ability of physical interaction between simulated objects is enlarged, allowing to simulate a porous structure, for example, reservoir rocks and structures with high density. The initial procedure for developing the simulation is the construction of porous cubic structure consisting of spheres with a single size and with varying sizes. In addition, structures can also be simulated with various volume fractions. The results presented are divided into two parts: first, the ball shall be deemed as solid grains, ie the matrix phase represents the porosity, the second, the spheres are considered as pores. In this case the matrix phase represents the solid phase. The simulations in both cases are the same, but the simulated structures are intrinsically different. To validate the results presented by the program, simulations were performed by varying the amount of grain, the grain size distribution and void fraction in the structure. All results showed statistically reliable and consistent with those presented in the literature. The mean values and distributions of stereological parameters measured, such as intercept linear section of perimeter area, sectional area and mean free path are in agreement with the results obtained in the literature for the structures simulated. The results may help the understanding of real structures.

Simulação computacional de medidas estereológicas em estruturas de metal duro (WC-Co)

This work focuses on the creation and applications of a dynamic simulation software in order to study the hard metal structure (WC-Co). The technological ground used to increase the GPU hardware capacity was Geforce 9600 GT along with the PhysX chip created to make games more realistic. The software simulates the three-dimensional carbide structure to the shape of a cubic box where tungsten carbide (WC) are modeled as triangular prisms and truncated triangular prisms. The program was proven effective regarding checking testes, ranging from calculations of parameter measures such as the capacity to increase the number of particles simulated dynamically. It was possible to make an investigation of both the mean parameters and distributions stereological parameters used to characterize the carbide structure through cutting plans. Grounded on the cutting plans concerning the analyzed structures, we have investigated the linear intercepts, the intercepts to the area, and the perimeter section of the intercepted grains as well as the binder phase to the structure by calculating the mean value and distribution of the free path. As literature shows almost consensually that the distribution of the linear intercepts is lognormal, this suggests that the grain distribution is also lognormal. Thus, a routine was developed regarding the program which made possible a more detailed research on this issue. We have observed that it is possible, under certain values for the parameters which define the shape and size of the Prismatic grain to find out the distribution to the linear intercepts that approach the lognormal shape. Regarding a number of developed simulations, we have observed that the distribution curves of the linear and area intercepts as well as the perimeter section are consistent with studies on static computer simulation to these parameters.

Modelagem e simulação de um sistema de bombeio mecânico em poços direcionais utilizando parâmetros concentrados

This work aims presenting the development of a model and computer simulation of a sucker rod pumping system. This system take into account the well geometry, the flow through the tubing, the dynamic behavior of the rod string and the use of a induction motor model. The rod string were modeled using concentrated parameters, allowing the use of ordinary differential equations systems to simulate it s behavior

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

Perdas de solo e nutrientes por erosão num Argissolo com resíduos vegetais de cana-de-açúcar

O presente trabalho teve como objetivo determinar as perdas de solo, matéria orgânica (MO) e nutrientes (P, K, Ca e Mg) por erosão, em área cultivada com cana-de-açúcar, cuja palhada é mantida sobre a superfície do solo, localizada em Catanduva - SP, sob um Argissolo Vermelho--Amarelo. em parcelas experimentais, com 0 (CS0), 50 (CS50) e 100% (CS100) de cobertura sobre a superfície do solo, aplicou-se chuva simulada com intensidade de 60 mm h-1, durante 65 minutos. Análises do sedimento erodido indicaram taxas de enriquecimento da seguinte ordem: 2,7 a 1,9 (MO), 3,8 a 2,7 (P), 1,3 a 1,7 (K), 3,9 a 3,6 (Ca) e 2,9 a 2,6 (Mg) vezes em relação ao solo original para CS0 e CS50, respectivamente. A CS50 propiciou controle significativo da erosão de 69%, mas não reduziu a concentração de MO e nutrientes no sedimento erodido. A CS100 foi significativamente eficiente no controle da erosão (89%) e na redução das concentrações de MO (69%), P (88%), K (23%), Ca (74%) e Mg (75%) no sedimento.

Desenvolvimento de uma ferramenta computacional para aplicação em métodos de elevação por bombeio centrífugo submerso

The Electrical Submersible Pump (ESP) has been one of the most appropriate solutions for lifting method in onshore and offshore applications. The typical features for this application are adverse temperature, viscosity fluids and gas environments. The difficulties in equipments maintenance and setup contributing to increasing costs of oil production in deep water, therefore, the optimization through automation can be a excellent approach for decrease costs and failures in subsurface equipment. This work describe a computer simulation related with the artificial lifting method ESP. This tool support the dynamic behavior of ESP approach, considering the source and electric energy transmission model for the motor, the electric motor model (including the thermal calculation), flow tubbing simulation, centrifugal pump behavior simulation with liquid nature effects and reservoir requirements. In addition, there are tri-dimensional animation for each ESP subsytem (transformer, motor, pump, seal, gas separator, command unit). This computer simulation propose a improvement for monitoring oil wells for maximization of well production. Currenty, the proprietaries simulators are based on specific equipments manufactures. Therefore, it is not possible simulation equipments of another manufactures. In the propose approach there are support for diverse kinds of manufactures equipments

Controle vetorial para velocidade de um motor de indução trifásico utilizando estimador filtro de Kalman

This paper describes the study, computer simulation and feasibility of implementation of vector control speed of an induction motor using for this purpose the Extended Kalman Filter as an estimator of rotor flux. The motivation for such work is the use of a control system that requires no sensors on the machine shaft, thus providing a considerable cost reduction of drives and their maintenance, increased reliability, robustness and noise immunity as compared to control systems with conventional sensors

The succession dynamics of a macroalgal community after a flood disturbance in a tropical stream from São Paulo State, southeastern Brazil

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

Propagação de danos em sistemas cooperativos: propriedades termodinâmicas

High-precision calculations of the correlation functions and order parameters were performed in order to investigate the critical properties of several two-dimensional ferro- magnetic systems: (i) the q-state Potts model; (ii) the Ashkin-Teller isotropic model; (iii) the spin-1 Ising model. We deduced exact relations connecting specific damages (the difference between two microscopic configurations of a model) and the above mentioned thermodynamic quanti- ties which permit its numerical calculation, by computer simulation and using any ergodic dynamics. The results obtained (critical temperature and exponents) reproduced all the known values, with an agreement up to several significant figures; of particular relevance were the estimates along the Baxter critical line (Ashkin-Teller model) where the exponents have a continuous variation. We also showed that this approach is less sensitive to the finite-size effects than the standard Monte-Carlo method. This analysis shows that the present approach produces equal or more accurate results, as compared to the usual Monte Carlo simulation, and can be useful to investigate these models in circumstances for which their behavior is not yet fully understood

Propagação de danos no modelo de ising em redes de bravais e em fractais

In this work we have studied, by Monte Carlo computer simulation, several properties that characterize the damage spreading in the Ising model, defined in Bravais lattices (the square and the triangular lattices) and in the Sierpinski Gasket. First, we investigated the antiferromagnetic model in the triangular lattice with uniform magnetic field, by Glauber dynamics; The chaotic-frozen critical frontier that we obtained coincides , within error bars, with the paramegnetic-ferromagnetic frontier of the static transition. Using heat-bath dynamics, we have studied the ferromagnetic model in the Sierpinski Gasket: We have shown that there are two times that characterize the relaxation of the damage: One of them satisfy the generalized scaling theory proposed by Henley (critical exponent z~A/T for low temperatures). On the other hand, the other time does not obey any of the known scaling theories. Finally, we have used methods of time series analysis to study in Glauber dynamics, the damage in the ferromagnetic Ising model on a square lattice. We have obtained a Hurst exponent with value 0.5 in high temperatures and that grows to 1, close to the temperature TD, that separates the chaotic and the frozen phases

Effect of target bias on magnetic field enhanced plasma immersion ion implantation

Recent studies have demonstrated that sheath dynamics in plasma immersion ion implantation (PIII) is significantly affected by an external magnetic field, especially in the case when the magnetic field is parallel to the workpiece surface or intersects it at small angles. In this work we report the results from two-dimensional, particle-in-cell (PIC) computer simulations of magnetic field enhanced plasma immersion implantation system at different bias voltages. The simulations begin with initial low-density nitrogen plasma, which extends with uniform density through a grounded cylindrical chamber. Negative bias voltage is applied to a cylindrical target located on the axis of the vacuum chamber. An axial magnetic field is created by a solenoid installed inside the target holder. A set of simulations at a fixed magnetic field of 0.0025 T at the target surface is performed. Secondary electron emission from the target subjected to ion bombardment is also included. It is found that the plasma density around the cylindrical target increases because of intense background gas ionization by the electrons drifting in the crossed E x B fields. Suppression of the sheath expansion and increase of the implantation current density in front of the high-density plasma region are observed. The effect of target bias on the sheath dynamics and implantation current of the magnetic field enhanced PIII is discussed. (C) 2007 Elsevier B.V. All rights reserved.

Effect of electron magnetic trapping in a plasma immersion ion implantation system

In this work we describe a two-dimensional computer simulation of magnetic field enhanced plasma immersion implantation system. Negative bias voltage of 10.0 kV is applied to a cylindrical target located on the axis of a grounded vacuum chamber filled with uniform nitrogen plasma. A pair of external coils creates a static magnetic field with main vector component along the axial direction. Thus, a system of crossed ExB field is generated inside the vessel forcing plasma electrons to rotate in azimuthal direction. In addition, the axial variation of the magnetic field intensity produces magnetic mirror effect that enables axial particle confinement. It is found that high-density plasma regions are formed around the target due to intense background gas ionization by the trapped electrons. Effect of the magnetic field on the sheath dynamics and the implantation current density of the PIII system is investigated. By changing the magnetic field axial profile (varying coils separation) an enhancement of about 30% of the retained dose can be achieved. The results of the simulation show that the magnetic mirror configuration brings additional benefits to the PIII process, permitting more precise control of the implanted dose.

A planar flexible robotic manipulator

Presents the dynamic modelling 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, which is used to discretize the displacements such that the small motion is represented in terms of nodal displacements. Computer simulation results are presented to illustrate this study. The dynamic model becomes necessary for use in future design and control applications.