Análise da estabilidade de sistemas dinâmicos periódicos usando Teoria de Sinha

Pós-graduação em Matemática - IBILCE

Continuous iteration of dynamical maps

A precise meaning is given to the notion of continuous iteration of a mapping. Usual discrete iterations are extended into a dynamical flow which is a homotopy of them all. The continuous iterate reveals that a dynamic map is formed by independent component modes evolving without interference with each other. An application to turbulent flow suggests that the velocity field assumes nonseparable values. © 1998 American Institute of Physics.

An effortless hybrid method to solve economic load dispatch problem in power systems

This paper proposes a new approach and coding scheme for solving economic dispatch problems (ED) in power systems through an effortless hybrid method (EHM). This novel coding scheme can effectively prevent futile searching and also prevents obtaining infeasible solutions through the application of stochastic search methods, consequently dramatically improves search efficiency and solution quality. The dominant constraint of an economic dispatch problem is power balance. The operational constraints, such as generation limitations, ramp rate limits, prohibited operating zones (POZ), network loss are considered for practical operation. Firstly, in the EHM procedure, the output of generator is obtained with a lambda iteration method and without considering POZ and later in a genetic based algorithm this constraint is satisfied. To demonstrate its efficiency, feasibility and fastness, the EHM algorithm was applied to solve constrained ED problems of power systems with 6 and 15 units. The simulation results obtained from the EHM were compared to those achieved from previous literature in terms of solution quality and computational efficiency. Results reveal that the superiority of this method in both aspects of financial and CPU time. (C) 2011 Elsevier Ltd. All rights reserved.

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.

An improved tabu-based vector optimal algorithm for design optimizations of electromagnetic devices

This paper introduces an improved tabu-based vector optimal algorithm for multiobjective optimal designs of electromagnetic devices. The improvements include a division of the entire search process, a new method for fitness assignment, a novel scheme for the generation and selection of neighborhood solutions, and so forth. Numerical results on a mathematical function and an engineering multiobjective design problem demonstrate that the proposed method can produce virtually the exact Pareto front, in both parameter and objective spaces, even though the iteration number used by it is only about 70% of that required by its ancestor.

Short term hydroelectric scheduling combining network flow and interior point approaches

In this paper, short term hydroelectric scheduling is formulated as a network flow optimization model and solved by interior point methods. The primal-dual and predictor-corrector versions of such interior point methods are developed and the resulting matrix structure is explored. This structure leads to very fast iterations since it avoids computation and factorization of impedance matrices. For each time interval, the linear algebra reduces to the solution of two linear systems, either to the number of buses or to the number of independent loops. Either matrix is invariant and can be factored off-line. As a consequence of such matrix manipulations, a linear system which changes at each iteration has to be solved, although its size is reduced to the number of generating units and is not a function of time intervals. These methods were applied to IEEE and Brazilian power systems, and numerical results were obtained using a MATLAB implementation. Both interior point methods proved to be robust and achieved fast convergence for all instances tested. (C) 2004 Elsevier Ltd. All rights reserved.

Nonlinear Dynamics and Control of an Ideal/Nonideal Load Transportation System With Periodic Coefficients

In this paper, a loads transportation system in platforms or suspended by cables is considered. It is a monorail device and is modeled as an inverted pendulum built on a car driven by a dc motor the governing equations of motion were derived via Lagrange's equations. In the mathematical model we consider the interaction between the dc motor and the dynamical system, that is, we have a so called nonideal periodic problem. The problem is analyzed, qualitatively, through the comparison of the stability diagrams, numerically obtained, for several motor torque constants. Furthermore, we also analyze the problem quantitatively using the Floquet multipliers technique. Finally, we devise a control for the studied nonideal problem. The method that was used for analysis and control of this nonideal periodic system is based on the Chebyshev polynomial exponsion, the Picard iterative method, and the Lyapunov-Floquet transformation (L-F transformation). We call it Sinha's theory.

INSULATOR-TO-METAL TRANSITION IN POLYTHIOPHENE

In the present work, the electronic structure of polythiophene at several doping levels is investigated by the use of the Huckel Hamiltonian with sigma-bond compressibility. Excess charges are assumed to be stored in conformational defects of the bipolaron type. The Hamiltonian matrix elements representative of a bipolaron are obtained from a previous thiophene oligomer calculation, and then transferred to very long chains. Negative factor counting and inverse iteration techniques have been used to evaluate densities of states and wave functions, respectively. Several types of defect distributions were analyzed. Our results are consistent with the following: (i) the bipolaron lattice does not present a finite density of states at the Fermi energy at any doping level; (ii) bipolaron clusters show an insulator-to-metal transition at 8 mol% doping level; (iii) segregation disorder shows an insulator-to-metal transition for doping levels in the range 20-30 mor %.

On a control of a non-ideal mono-rail system with periodics coefficients

In this work, the problem in the loads transport (in platforms or suspended by cables) it is considered. The system in subject is composed for mono-rail system and was modeled through the system: inverted pendulum, car and motor and the movement equations were obtained through the Lagrange equations. In the model, was considered the interaction among of the motor and system dynamics for several potencies motor, that is, the case studied is denominated a non-ideal periodic problem. The non-ideal periodic problem dynamics was analyzed, qualitatively, through the comparison of the stability diagrams, numerically obtained, for several motor torque constants. Furthermore, one was made it analyzes quantitative of the problem through the analysis of the Floquet multipliers. Finally, the non-ideal problem was controlled. The method that was used for analysis and control of non-ideal periodic systems is based on the Chebyshev polynomial expansion, in the Picard iterative method and in the Lyapunov-Floquet transformation (L-F trans formation). This method was presented recently in [3-9].

On non-linear dynamics and a periodic control design applied to the potential of membrane action

The Fitzhugh-Nagumo (fn) mathematical model characterizes the action potential of the membrane. The dynamics of the Fitzhugh-Nagumo model have been extensively studied both with a view to their biological implications and as a test bed for numerical methods, which can be applied to more complex models. This paper deals with the dynamics in the (FH) model. Here, the dynamics are analyzed, qualitatively, through the stability diagrams to the action potential of the membrane. Furthermore, we also analyze quantitatively the problem through the evaluation of Floquet multipliers. Finally, the nonlinear periodic problem is controlled, based on the Chebyshev polynomial expansion, the Picard iterative method and on Lyapunov-Floquet transformation (L-F transformation).

Energy and protein requirements for maintenance and growth of Boer crossbred kids

Meat production by goats has become an important livestock enterprise in several parts of the world. Nonetheless, energy and protein requirements of meat goats have not been defined thoroughly. The objective of this study was to determine the energy and protein requirements for maintenance and growth of 34 3/4 Boer x 1/4 Saanen crossbred, intact male kids (20.5 +/- 0.24 kg of initial BW). The baseline group was 7 randomly selected kids, averaging 21.2 +/- 0.36 kg of BW. An intermediate group consisted of 6 randomly selected kids, fed for ad libitum intake, that were slaughtered when they reached an average BW of 28.2 +/- 0.39 kg. The remaining kids (n = 21) were allocated randomly on d 0 to 3 levels of DMI (treatments were ad libitum or restricted to 70 or 40% of the ad libitum intake) within 7 slaughter groups. A slaughter group contained 1 kid from each treatment, and kids were slaughtered when the ad libitum treatment kid reached 35 kg of BW. Individual body components (head plus feet, hide, internal organs plus blood, and carcass) were weighed, ground, mixed, and subsampled for chemical analyses. Initial body composition was determined using equations developed from the composition of the baseline kids. The calculated daily maintenance requirement for NE was 77.3 +/- 1.05 kcal/kg(0.75) of empty BW (EBW) or 67.4 +/- 1.04 kcal/kg(0.75) of shrunk BW. The daily ME requirement for maintenance (118.1 kcal/g(0.75) of EBW or 103.0 kcal/kg(0.75) of shrunk BW) was calculated by iteration, assuming that the heat produced was equal to the ME intake at maintenance. The partial efficiency of use of ME for NE below maintenance was 0.65. A value of 2.44 +/- 0.4 g of net protein/kg(0.75) of EBW for daily maintenance was determined. Net energy requirements for growth ranged from 2.55 to 3.0 Mcal/kg of EBW gain at 20 and 35 kg of BW, and net protein requirements for growth ranged from 178.8 to 185.2 g/kg of EBW gain. These results suggest that NE and net protein requirements for growing meat goats exceed the requirements previously published for dairy goats. Moreover, results from this study suggest that the N requirement for maintenance for growing goats is greater than the established recommendations.

Optimisation of a fermentation process for butanol production by particle swarm optimisation (PSO)

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

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.

Numerical study of the spherically symmetric Gross-Pitaevskii equation in two space dimensions

The Gross-Pitaevskii equation for Bose-Einstein condensation (BEC) in two space dimensions under the action of a harmonic oscillator trap potential for bosonic atoms with attractive and repulsive interparticle interactions was numerically studied by using time-dependent and time-independent approaches. In both cases, numerical difficulty appeared for large nonlinearity. Nonetheless, the solution of the time-dependent approach exhibited intrinsic oscillation with time iteration which is independent of space and time steps used in discretization.

On the determination of epsilon during discriminative GMM training

Discriminative training of Gaussian Mixture Models (GMMs) for speech or speaker recognition purposes is usually based on the gradient descent method, in which the iteration step-size, ε, uses to be defined experimentally. In this letter, we derive an equation to adaptively determine ε, by showing that the second-order Newton-Raphson iterative method to find roots of equations is equivalent to the gradient descent algorithm. © 2010 IEEE.