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.

Modelling the dynamics of dengue real epidemics

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

Dynamics in dumbbell domains II. The limiting problem

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

Riccati-type equations, generalised WZNW equations, and multidimensional Toda systems

We associate to an arbitrary Z-gradation of the Lie algebra of a Lie group a system of Riccati-type first order differential equations. The particular cases under consideration are the ordinary Riccati and the matrix Riccati equations. The multidimensional extension of these equations is given. The generalisation of the associated Redheffer-Reid differential systems appears in a natural way. The connection between the Toda systems and the Riccati-type equations in lower and higher dimensions is established. Within this context the integrability problem for those equations is studied. As an illustration, some examples of the integrable multidimensional Riccati-type equations related to the maximally nonabelian Toda systems are given.

Bose-Einstein condensation dynamics in three dimensions by the pseudospectral and finite-difference methods

We suggest a pseudospectral method for solving the three-dimensional time-dependent Gross-Pitaevskii (GP) equation, and use it to study the resonance dynamics of a trapped Bose-Einstein condensate induced by a periodic variation in the atomic scattering length. When the frequency of oscillation of the scattering length is an even multiple of one of the trapping frequencies along the x, y or z direction, the corresponding size of the condensate executes resonant oscillation. Using the concept of the differentiation matrix, the partial-differential GP equation is reduced to a set of coupled ordinary differential equations, which is solved by a fourth-order adaptive step-size control Runge-Kutta method. The pseudospectral method is contrasted with the finite-difference method for the same problem, where the time evolution is performed by the Crank-Nicholson algorithm. The latter method is illustrated to be more suitable for a three-dimensional standing-wave optical-lattice trapping potential.

C programs for solving the time-dependent Gross-Pitaevskii equation in a fully anisotropic trap

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

On Lugre Friction Model to Mitigate Nonideal Vibrations

In this paper, a nonideal mechanical system with the LuGre friction damping model is considered. The mechanical model of the system is an oscillator not necessarily linear connected with an unbalanced motor of excitation with limited power supply. The control of motion and the attenuation of the Sommerfeld effect of the considered nonideal system are analyzed in this paper The mathematical model of the system is represented by coupled non-linear differential equations. The identification of some interesting nonlinear phenomenon in the transient and steady state motion of the system during the passage through resonance (using applied voltages at dc motor as control parameter) is investigated in detail using numerical simulation. [DOI: 10.1115/1.3124783]

Uniform exponential dichotomy and continuity of attractors for singularly perturbed damped wave equations

For eta >= 0, we consider a family of damped wave equations u(u) + eta Lambda 1/2u(t) + au(t) + Lambda u = f(u), t > 0, x is an element of Omega subset of R-N, where -Lambda denotes the Laplacian with zero Dirichlet boundary condition in L-2(Omega). For a dissipative nonlinearity f satisfying a suitable growth restrictions these equations define on the phase space H-0(1)(Omega) x L-2(Omega) semigroups {T-eta(t) : t >= 0} which have global attractors A(eta) eta >= 0. We show that the family {A(eta)}(eta >= 0), behaves upper and lower semi-continuously as the parameter eta tends to 0(+).

Dynamics in dumbbell domains I. Continuity of the set of equilibria

We analyze the dynamics of a reaction-diffusion equation with homogeneous Neumann boundary conditions in a dumbbell domain. We provide an appropriate functional setting to treat this problem and, as a first step, we show in this paper the continuity of the set of equilibria and of its linear unstable manifolds. (c) 2006 Elsevier B.V. All rights reserved.

Dissipative three-wave structures in stimulated backscattering. I. A subluminous solitary attractor

We present a solitary solution of the three-wave nonlinear partial differential equation (PDE) model - governing resonant space-time stimulated Brillouin or Raman backscattering - in the presence of a cw pump and dissipative material and Stokes waves. The study is motivated by pulse formation in optical fiber experiments. As a result of the instability any initial bounded Stokes signal is amplified and evolves to a subluminous backscattered Stokes pulse whose shape and velocity are uniquely determined by the damping coefficients and the cw-pump level. This asymptotically stable solitary three-wave structure is an attractor for any initial conditions in a compact support, in contrast to the known superluminous dissipative soliton solution which calls for an unbounded support. The linear asymptotic theory based on the Kolmogorov-Petrovskii-Piskunov assertion allows us to determine analytically the wave-front slope and the subluminous velocity, which are in remarkable agreement with the numerical computation of the nonlinear PDE model when the dynamics attains the asymptotic steady regime. © 1997 The American Physical Society.

Stable two-dimensional dispersion-managed soliton

The existence of a dispersion-managed soliton in two-dimensional nonlinear Schrodinger equation with periodically varying dispersion has been explored. The averaged equations for the soliton width and chirp are obtained which successfully describe the long time evolution of the soliton. The slow dynamics of the soliton around the fixed points for the width and chirp are investigated and the corresponding frequencies are calculated. Analytical predictions are confirmed by direct partial differential equation (PDE) and ordinary differential equation (ODE) simulations. Application to a Bose-Einstein condensate in optical lattice is discussed. The existence of a dispersion-managed matter-wave soliton in such system is shown.

Substitution-newton-Raphson method for the solution of electric network equations

The power flow problem, in transmission networks, has been well solved, for most cases, using Newton-Raphson method (NR) and its decoupled versions. Generally speaking, the solution of a non-linear system of equations refers to two methods: NR and Successive Substitution. The proposal of this paper is to evaluate the potential of the Substitution-Newton-Raphson Method (SNR), which combines both methods, on the solution of the power flow problem. Simulations were performed using a two-bus test network in order to observe the characteristics of these methods. It was verified that the NR is faster than SNR, in terms of convergence, considering non-stressed scenarios. For those cases where the power flow in the network is closed to the limits (stressed system), the SNR converges faster. This paper presents the power flow formulation of the SNR and describes its potential for its application in special cases such as stressed scenarios. © 2006 IEEE.

Transmission problem for waves with frictional damping April 22, 2007

In this paper we consider the transmission problem, in one space dimension, for linear dissipative waves with frictional damping. We study the wave propagation in a medium with a component with attrition and another simply elastic. We show that for this type of material, the dissipation produced by the frictional part is strong enough to produce exponential decay of the solution, no matter how small is its size. ©2007 Texas State University.

Consequences of quadratic frictional force on the one dimensional bouncing ball model

Some dynamical properties of the one dimensional Fermi accelerator model, under the presence of frictional force are studied. The frictional force is assumed as being proportional to the square particle's velocity. The problem is described by use of a two dimensional non linear mapping, therefore obtained via the solution of differential equations. We confirm that the model experiences contraction of the phase space area and in special, we characterized the behavior of the particle approaching an attracting fixed point. © 2007 American Institute of Physics.

Necessary conditions of optimality for state constrained infinite horizon differential inclusions

This article presents and discusses necessary conditions of optimality for infinite horizon dynamic optimization problems with inequality state constraints and set inclusion constraints at both endpoints of the trajectory. The cost functional depends on the state variable at the final time, and the dynamics are given by a differential inclusion. Moreover, the optimization is carried out over asymptotically convergent state trajectories. The novelty of the proposed optimality conditions for this class of problems is that the boundary condition of the adjoint variable is given as a weak directional inclusion at infinity. This improves on the currently available necessary conditions of optimality for infinite horizon problems. © 2011 IEEE.