Um método de linearização local com passo adaptativo para solução numérica de equações diferenciais estocásticas com ruído aditivo

Neste trabalho apresentamos um novo método numérico com passo adaptativo baseado na abordagem de linearização local, para a integração de equações diferenciais estocásticas com ruído aditivo. Propomos, também, um esquema computacional que permite a implementação eficiente deste método, adaptando adequadamente o algorítimo de Padé com a estratégia “scaling-squaring” para o cálculo das exponenciais de matrizes envolvidas. Antes de introduzirmos a construção deste método, apresentaremos de forma breve o que são equações diferenciais estocásticas, a matemática que as fundamenta, a sua relevância para a modelagem dos mais diversos fenômenos, e a importância da utilização de métodos numéricos para avaliar tais equações. Também é feito um breve estudo sobre estabilidade numérica. Com isto, pretendemos introduzir as bases necessárias para a construção do novo método/esquema. Ao final, vários experimentos numéricos são realizados para mostrar, de forma prática, a eficácia do método proposto, e compará-lo com outros métodos usualmente utilizados.

Symmetry analysis of an integrable reaction-diffusion equation

In this paper, we investigate the invariance and integrability properties of an integrable two-component reaction-diffusion equation. We perform Painleve analysis for both the reaction-diffusion equation modelled by a coupled nonlinear partial differential equations and its general similarity reduced ordinary differential equation and confirm its integrability. Further, we perform Lie symmetry analysis for this model. Interestingly our investigations reveals a rich variety of particular solutions, which have not been reported in the literature, for this model. (C) 2000 Elsevier B.V. Ltd. All rights reserved.

General soliton matrices in the Riemann-Hilbert problem for integrable nonlinear equations

We derive the soliton matrices corresponding to an arbitrary number of higher-order normal zeros for the matrix Riemann-Hilbert problem of arbitrary matrix dimension, thus giving the complete solution to the problem of higher-order solitons. Our soliton matrices explicitly give all higher-order multisoliton solutions to the nonlinear partial differential equations integrable through the matrix Riemann-Hilbert problem. We have applied these general results to the three-wave interaction system, and derived new classes of higher-order soliton and two-soliton solutions, in complement to those from our previous publication [Stud. Appl. Math. 110, 297 (2003)], where only the elementary higher-order zeros were considered. The higher-order solitons corresponding to nonelementary zeros generically describe the simultaneous breakup of a pumping wave (u(3)) into the other two components (u(1) and u(2)) and merger of u(1) and u(2) waves into the pumping u(3) wave. The two-soliton solutions corresponding to two simple zeros generically describe the breakup of the pumping u(3) wave into the u(1) and u(2) components, and the reverse process. In the nongeneric cases, these two-soliton solutions could describe the elastic interaction of the u(1) and u(2) waves, thus reproducing previous results obtained by Zakharov and Manakov [Zh. Eksp. Teor. Fiz. 69, 1654 (1975)] and Kaup [Stud. Appl. Math. 55, 9 (1976)]. (C) 2003 American Institute of Physics.

Noise and ultraviolet divergences in the dynamics of the chiral condensate in QCD

The time evolution of the matter produced in high energy heavy-ion collisions seems to be well described by relativistic viscous hydrodynamics. In addition to the hydrodynamic degrees of freedom related to energy-momentum conservation, degrees of freedom associated with order parameters of broken continuous symmetries must be considered because they are all coupled to each other. of particular interest is the coupling of degrees of freedom associated with the chiral symmetry of QCD. Quantum and thermal fluctuations of the chiral fields act as noise sources in the classical equations of motion, turning them into stochastic differential equations in the form of Ginzburg-Landau-Langevin (GLL) equations. Analytic solutions of GLL equations are attainable only in very special circumstances and extensive numerical simulations are necessary, usually by discretizing the equations on a spatial lattice. However, a not much appreciated issue in the numerical simulations of GLL equations is that ultraviolet divergences in the form of lattice-spacing dependence plague the solutions. The divergences are related to the well-known Rayleigh-Jeans catastrophe in classical field theory. In the present communication we present a systematic lattice renormalization method to control the catastrophe. We discuss the implementation of the method for a GLL equation derived in the context of a model for the QCD chiral phase transition and consider the nonequilibrium evolution of the chiral condensate during the hydrodynamic flow of the quark-gluon plasma.

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

ALPHA-CONTRACTIONS AND ATTRACTORS FOR DISSIPATIVE SEMILINEAR HYPERBOLIC-EQUATIONS AND SYSTEMS

In this paper we discuss the existence of compact attractor for the abstract semilinear evolution equation u = Au + f (t, u); the results are applied to damped partial differential equations of hyperbolic type. Our approach is a combination of Liapunov method with the theory of alpha-contractions.

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.

Modelos matemáticos e aplicações ao ensino médio

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

Teoria de estabilidade de equações diferenciais ordinárias e aplicações: modelos presa-predador e competição entre espécies

Pós-graduação em Matemática Universitária - IGCE

Description of the attainable sets of one-dimensional differential inclusions

The role played by the attainable set of a differential inclusion, in the study of dynamic control systems and fuzzy differential equations, is widely acknowledged. A procedure for estimating the attainable set is rather complicated compared to the numerical methods for differential equations. This article addresses an alternative approach, based on an optimal control tool, to obtain a description of the attainable sets of differential inclusions. In particular, we obtain an exact delineation of the attainable set for a large class of nonlinear differential inclusions.

Soluções numéricas de Equações Diferenciais Fracionária e Ordinária para o cálculo da distribuição da concentração em um sistema de compartimento

The purpose of this work was the study of numerical methods for differential equations of fractional order and ordinary. These methods were applied to the problem of calculating the distribution of the concentration of a given substance over time in a given physical system. The two compartment model was used for representation of this system. Comparison between numerical solutions obtained were performed and, in particular, also compared with the analytical solution of this problem. Finally, estimates for the error between the solutions were calculated

Differential orthogonality: Laguerre and Hermite cases with applications

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

Differential cross sections for elastic and inelastic positronium-hydrogen-atom scattering

A time reversal symmetric regularized electron exchange model was used to elastic scattering, target elastic Ps excitations and target inelastic excitation of hydrogen in a five state coupled model. A singlet Ps-H-S-wave resonance at 4.01 eV of width 0.15 eV and a P-wave resonance at 5.08 eV of width 0.004 eV were obtained using this model. The effect on the convergence of the coupled-channel scheme due to the inclusion of the excited Ps and H states was also analyzed.

Estabilidade de pontos de equilíbrio e existência de soluções periódicas em alguns modelos bidimensionais

Pós-graduação em Matemática Universitária - IGCE

Continuity of Dynamical Structures for Nonautonomous Evolution Equations Under Singular Perturbations

In this paper we study the continuity of invariant sets for nonautonomous infinite-dimensional dynamical systems under singular perturbations. We extend the existing results on lower-semicontinuity of attractors of autonomous and nonautonomous dynamical systems. This is accomplished through a detailed analysis of the structure of the invariant sets and its behavior under perturbation. We prove that a bounded hyperbolic global solutions persists under singular perturbations and that their nonlinear unstable manifold behave continuously. To accomplish this, we need to establish results on roughness of exponential dichotomies under these singular perturbations. Our results imply that, if the limiting pullback attractor of a nonautonomous dynamical system is the closure of a countable union of unstable manifolds of global bounded hyperbolic solutions, then it behaves continuously (upper and lower) under singular perturbations.