Hartman-Grobman decomposition when in presence of jumps in the relative degree of the dynamics in an energy harvesting device

This work deals with the nonlinear piezoelectric coupling in vibration-based energy harvesting, done by A. Triplett and D.D. Quinn in J. of Intelligent Material Syst. and Structures (2009). In that paper the first order nonlinear fundamental equation has a three dimensional state variable. Introducing both observable and control variables in such a way the controlled system became a SISO system, we can obtain as a corollary that for a particular choice of the observable variable it is possible to present an explicit functional relation between this variable one, and the variable representing the charge harvested. After-by observing that the structure in the Input-Output decomposition essentially changes depending on the relative degree changes, presenting bifurcation branches in its zero dynamics-we are able in to identify this type of bifurcation indicating its close relation with the Hartman - Grobman theorem telling about decomposition into stable and the unstable manifolds for hyperbolic points.

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).

On Complementary Root Locus of Biproper Transfer Functions

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

Um estudo algorítmico de problemas logísticos na indústria de petróleo e gás natural

This work consists on the study of two important problems arising from the operations of petroleum and natural gas industries. The first problem the pipe dimensioning problem on constrained gas distribution networks consists in finding the least cost combination of diameters from a discrete set of commercially available ones for the pipes of a given gas network, such that it respects minimum pressure requirements at each demand node and upstream pipe conditions. On its turn, the second problem the piston pump unit routing problem comes from the need of defining the piston pump unit routes for visiting a number of non-emergent wells in on-shore fields, i.e., wells which don t have enough pressure to make the oil emerge to surface. The periodic version of this problem takes into account the wells re-filling equation to provide a more accurate planning in the long term. Besides the mathematical formulation of both problems, an exact algorithm and a taboo search were developed for the solution of the first problem and a theoretical limit and a ProtoGene transgenetic algorithm were developed for the solution of the second problem. The main concepts of the metaheuristics are presented along with the details of their application to the cited problems. The obtained results for both applications are promising when compared to theoretical limits and alternate solutions, either relative to the quality of the solutions or to associated running time

Lobatto deferred correction for stiff two-point boundary value problems

An iterated deferred correction algorithm based on Lobatto Runge-Kutta formulae is developed for the efficient numerical solution of nonlinear stiff two-point boundary value problems. An analysis of the stability properties of general deferred correction schemes which are based on implicit Runge-Kutta methods is given and results which are analogous to those obtained for initial value problems are derived. A revised definition of symmetry is presented and this ensures that each deferred correction produces an optimal increase in order. Finally, some numerical results are given to demonstrate the superior performance of Lobatto formulae compared with mono-implicit formulae on stiff two-point boundary value problems. (C) 1998 Elsevier B.V. Ltd. All rights reserved.

Formation of bright solitons and soliton trains in a fermion-fermion mixture by modulational instability

We employ a time- dependent mean- field- hydrodynamic model to study the generation of bright solitons in a degenerate fermion - fermion mixture in a cigar- shaped geometry using variational and numerical methods. Due to a strong Pauli- blocking repulsion among identical spin- polarized fermions at short distances there cannot be bright solitons for repulsive interspecies interactions. Employing a linear stability analysis we demonstrate the formation of stable solitons due to modulational instability of a constant-amplitude solution of the model equations for a sufficiently attractive interspecies interaction. We perform a numerical stability analysis of these solitons and also demonstrate the formation of soliton trains by jumping the effective interspecies interaction from repulsive to attractive. These fermionic solitons can be formed and studied in laboratory with present technology.

An exact equation for the free surface of a fluid in a porous medium

We study the problem of the evolution of the free surface of a fluid in a saturated porous medium, bounded from below by a. at impermeable bottom, and described by the Laplace equation with moving-boundary conditions. By making use of a convenient conformal transformation, we show that the solution to this problem is equivalent to the solution of the Laplace equation on a fixed domain, with new variable coefficients, the boundary conditions. We use a kernel of the Laplace equation which allows us to write the Dirichlet-to-Neumann operator, and in this way we are able to find an exact differential-integral equation for the evolution of the free surface in one space dimension. Although not amenable to direct analytical solutions, this equation turns out to allow an easy numerical implementation. We give an explicit illustrative case at the end of the article.

Effect of a frictional force on the Fermi-Ulam model

The dynamical properties of a classical particle bouncing between two rigid walls, in the presence of a drag force, are studied for the case where one wall is fixed and the other one moves periodically in time. The system is described in terms of a two-dimensional nonlinear map obtained by solution of the relevant differential equations. It is shown that the structure of the KAM curves and the chaotic sea is destroyed as the drag force is introduced. At high energy, the velocity of the particle decreases linearly with increasing iteration number, but with a small superimposed sinusoidal modulation. If the motion passes near enough to a fixed point, the particle approaches it exponentially as the iteration number evolves, with a speed of approach that depends on the strength of the drag force. For a simplified version of the model it is shown that, at low energies corresponding to the region of the chaotic sea in the non-dissipative model, the particle wanders in a chaotic transient that depends on the strength of the drag coefficient. However, the KAM islands survive in the presence of dissipation. It is confirmed that the fixed points and periodic orbits go over smoothly into the orbits of the well-known (non-dissipative) Fermi-Ulam model as the drag force goes to zero.

Can Drag Force Suppress Fermi Acceleration in a Bouncer Model?

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

Scaling Properties of a Hybrid Fermi-Ulam-Bouncer Model

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

Phase Transition in Dynamical Systems: Defining Classes of Universality for Two-Dimensional Hamiltonian Mappings via Critical Exponents

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

Variation of the Equator due to a Highly Inclined and Eccentric Disturber

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

Rapidly varying boundaries in equations with nonlinear boundary conditions. The case of a Lipschitz deformation

We analyze the behavior of solutions of nonlinear elliptic equations with nonlinear boundary conditions of type partial derivative u/partial derivative n + g( x, u) = 0 when the boundary of the domain varies very rapidly. We show that the limit boundary condition is given by partial derivative u/partial derivative n+gamma(x) g(x, u) = 0, where gamma(x) is a factor related to the oscillations of the boundary at point x. For the case where we have a Lipschitz deformation of the boundary,. is a bounded function and we show the convergence of the solutions in H-1 and C-alpha norms and the convergence of the eigenvalues and eigenfunctions of the linearization around the solutions. If, moreover, a solution of the limit problem is hyperbolic, then we show that the perturbed equation has one and only one solution nearby.

VERY RAPIDLY VARYING BOUNDARIES IN EQUATIONS WITH NONLINEAR BOUNDARY CONDITIONS. THE CASE of A NON UNIFORMLY LIPSCHITZ DEFORMATION

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