Simulador computacional para poços de petróleo com método de elevação artificial por bombeio mecânico

This work proposes a computer simulator for sucker rod pumped vertical wells. The simulator is able to represent the dynamic behavior of the systems and the computation of several important parameters, allowing the easy visualization of several pertinent phenomena. The use of the simulator allows the execution of several tests at lower costs and shorter times, than real wells experiments. The simulation uses a model based on the dynamic behavior of the rod string. This dynamic model is represented by a second order partial differencial equation. Through this model, several common field situations can be verified. Moreover, the simulation includes 3D animations, facilitating the physical understanding of the process, due to a better visual interpretation of the phenomena. Another important characteristic is the emulation of the main sensors used in sucker rod pumping automation. The emulation of the sensors is implemented through a microcontrolled interface between the simulator and the industrial controllers. By means of this interface, the controllers interpret the simulator as a real well. A "fault module" was included in the simulator. This module incorporates the six more important faults found in sucker rod pumping. Therefore, the analysis and verification of these problems through the simulator, allows the user to identify such situations that otherwise could be observed only in the field. The simulation of these faults receives a different treatment due to the different boundary conditions imposed to the numeric solution of the problem. Possible applications of the simulator are: the design and analysis of wells, training of technicians and engineers, execution of tests in controllers and supervisory systems, and validation of control algorithms

Avaliacao de dados de secagem de suspensoes de polpas de frutas em leito de jorro com alimentacao intermitente

The objective of this work was the development and improvement of the mathematical models based on mass and heat balances, representing the drying transient process fruit pulp in spouted bed dryer with intermittent feeding. Mass and energy balance for drying, represented by a system of differential equations, were developed in Fortran language and adapted to the condition of intermittent feeding and mass accumulation. Were used the DASSL routine (Differential Algebraic System Solver) for solving the differential equation system and used a heuristic optimization algorithm in parameter estimation, the Particle Swarm algorithm. From the experimental data food drying, the differential models were used to determine the quantity of water and the drying air temperature at the exit of a spouted bed and accumulated mass of powder in the dryer. The models were validated using the experimental data of drying whose operating conditions, air temperature, flow rate and time intermittency, varied within the limits studied. In reviewing the results predicted, it was found that these models represent the experimental data of the kinetics of production and accumulation of powder and humidity and air temperature at the outlet of the dryer

Evolução top-hat hidrodinâmica em campos de velocidades peculiares primordiais

We investigate the cosmology of the vacuum energy decaying into cold dark matter according to thermodynamics description of Alcaniz & Lima. We apply this model to analyze the evolution of primordial density perturbations in the matter that gave rise to the first generation of structures bounded by gravity in the Universe, called Population III Objects. The analysis of the dynamics of those systems will involve the calculation of a differential equation system governing the evolution of perturbations to the case of two coupled fluids (dark matter and baryonic matter), modeled with a Top-Hat profile based in the perturbation of the hydrodynamics equations, an efficient analytical tool to study the properties of dark energy models such as the behavior of the linear growth factor and the linear growth index, physical quantities closely related to the fields of peculiar velocities at any time, for different models of dark energy. The properties and the dynamics of current Universe are analyzed through the exact analytical form of the linear growth factor of density fluctuations, taking into account the influence of several physical cooling mechanisms acting on the density fluctuations of the baryonic component of matter during the evolution of the clouds of matter, studied from the primordial hydrogen recombination. This study is naturally extended to more general models of dark energy with constant equation of state parameter in a flat Universe

Métodos numéricos para resolução de equações diferenciais ordinárias lineares baseados em interpolação por spline

In this work we have elaborated a spline-based method of solution of inicial value problems involving ordinary differential equations, with emphasis on linear equations. The method can be seen as an alternative for the traditional solvers such as Runge-Kutta, and avoids root calculations in the linear time invariant case. The method is then applied on a central problem of control theory, namely, the step response problem for linear EDOs with possibly varying coefficients, where root calculations do not apply. We have implemented an efficient algorithm which uses exclusively matrix-vector operations. The working interval (till the settling time) was determined through a calculation of the least stable mode using a modified power method. Several variants of the method have been compared by simulation. For general linear problems with fine grid, the proposed method compares favorably with the Euler method. In the time invariant case, where the alternative is root calculation, we have indications that the proposed method is competitive for equations of sifficiently high order.

Scattering and bound states of spinless particles in a mixed vector-scalar smooth step potential

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

Colombeau's theory and shock wave solutions for systems of PDEs

In this article we study the existence of shock wave solutions for systems of partial differential equations of hydrodynamics with viscosity in one space dimension in the context of Colombeau's theory of generalized functions. This study uses the equality in the strict sense and the association of generalized functions (that is the weak equality). The shock wave solutions are given in terms of generalized functions that have the classical Heaviside step function as macroscopic aspect. This means that solutions are sought in the form of sequences of regularizations to the Heaviside function that have to satisfy part of the equations in the strict sense and part of the equations in the sense of association.

Dynamic stationary response of reinforced plates by the boundary element method

A direct version of the boundary element method (BEM) is developed to model the stationary dynamic response of reinforced plate structures, such as reinforced panels in buildings, automobiles, and airplanes. The dynamic stationary fundamental solutions of thin plates and plane stress state are used to transform the governing partial differential equations into boundary integral equations (BIEs). Two sets of uncoupled BIEs are formulated, respectively, for the in-plane state ( membrane) and for the out-of-plane state ( bending). These uncoupled systems are joined to formamacro-element, in which membrane and bending effects are present. The association of these macro-elements is able to simulate thin-walled structures, including reinforced plate structures. In the present formulation, the BIE is discretized by continuous and/or discontinuous linear elements. Four displacement integral equations are written for every boundary node. Modal data, that is, natural frequencies and the corresponding mode shapes of reinforced plates, are obtained from information contained in the frequency response functions (FRFs). A specific example is presented to illustrate the versatility of the proposed methodology. Different configurations of the reinforcements are used to simulate simply supported and clamped boundary conditions for the plate structures. The procedure is validated by comparison with results determined by the finite element method (FEM).

Traveling waves in the Lethargic Crab Disease

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

Zeros of Jacobi functions of second kind

The number of zeros in (- 1, 1) of the Jacobi function of second kind Q(n)((alpha, beta)) (x), alpha, beta > - 1, i.e. The second solution of the differential equation(1 - x(2))y (x) + (beta - alpha - (alpha + beta + 2)x)y' (x) + n(n + alpha + beta + 1)y(x) = 0,is determined for every n is an element of N and for all values of the parameters alpha > - 1 and beta > - 1. It turns out that this number depends essentially on alpha and beta as well as on the specific normalization of the function Q(n)((alpha, beta)) (x). Interlacing properties of the zeros are also obtained. As a consequence of the main result, we determine the number of zeros of Laguerre's and Hermite's functions of second kind. (c) 2005 Elsevier B.V. All rights reserved.

New steps on Sobolev orthogonality in two variables

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

Superfluid Fermi-Fermi mixture: Phase diagram, stability, and soliton formation

We study the phase diagram for a dilute Bardeen-Cooper-Schrieffer superfluid Fermi-Fermi mixture (of distinct mass) at zero temperature using energy densities for the superfluid fermions in one (1D), two (2D), and three (3D) dimensions. We also derive the dynamical time-dependent nonlinear Euler-Lagrange equation satisfied by the mixture in one dimension using this energy density. We obtain the linear stability conditions for the mixture in terms of fermion densities of the components and the interspecies Fermi-Fermi interaction. In equilibrium there are two possibilities. The first is that of a uniform mixture of the two components, the second is that of two pure phases of two components without any overlap between them. In addition, a mixed and a pure phase, impossible in 1D and 2D, can be created in 3D. We also obtain the conditions under which the uniform mixture is stable from an energetic consideration. The same conditions are obtained from a modulational instability analysis of the dynamical equations in 1D. Finally, the 1D dynamical equations for the system are solved numerically and by variational approximation (VA) to study the bright solitons of the system for attractive interspecies Fermi-Fermi interaction in 1D. The VA is found to yield good agreement to the numerical result for the density profile and chemical potential of the bright solitons. The bright solitons are demonstrated to be dynamically stable. The experimental realization of these Fermi-Fermi bright solitons seems possible with present setups.

Homoclinic bifurcations in reversible Hamiltonian systems

We study the existence of homoclic solutions for reversible Hamiltonian systems taking the family of differential equations u(iv) + au - u +f(u, b) = 0 as a model, where fis an analytic function and a, b real parameters. These equations are important in several physical situations such as solitons and in the existence of finite energy stationary states of partial differential equations, but no assumptions of any kind of discrete symmetry is made and the analysis here developed can be extended to others Hamiltonian systems and successfully employed in situations where standard methods fail. We reduce the problem of computing these orbits to that of finding the intersection of the unstable manifold with a suitable set and then apply it to concrete situations. We also plot the homoclinic values configuration in parameters space, giving a picture of the structural distribution and a geometrical view of homoclinic bifurcations. (c) 2005 Published by Elsevier B.V.

Mean-field equations for cigar- and disc-shaped Bose and Fermi superfluids

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

Correlation times in stochastic equations with delayed feedback and multiplicative noise

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

Negative Even Grade mKdV Hierarchy and its Soliton Solutions

In this paper we discuss the algebraic construction of the mKdV hierarchy in terms of an affine Lie algebra (s) over capl(2). An interesting novelty araises from the negative even grade sector of the affine algebra leading to nonlinear integro-differential equations admiting non-trivial vacuum configuration. These solitons solutions are constructed systematically from generalization of the dressing method based on non zero vacua. The sub-hierarchies admiting such class of solutions are classified.