Limit cycles of cubic polynomial differential systems with rational first integrals of degree 2

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

Secular series and renormalization group for amplitude equations

We have developed a technique that circumvents the process of elimination of secular terms and reproduces the uniformly valid approximations, amplitude equations, and first integrals. The technique is based on a rearrangement of secular terms and their grouping into the secular series that multiplies the constants of the asymptotic expansion. We illustrate the technique by deriving amplitude equations for standard nonlinear oscillator and boundary-layer problems. © 2008 The American Physical Society.

Integrability and singularity structure of predator‐prey system

The ‘‘extended’’ ARS (Ablowitz, Ramani, and Segur) algorithm is introduced to characterize a dynamical system as Painlevé or otherwise; to that end, it is required that the formal series—the Laurent series, logarithmic, algebraic psi series about a movable singularity—are shown to converge in the deleted neighborhood of the singularity. The determinations thus obtained are compared with those following from the α method of Painlevé. An attempt is made to relate the structure of solutions about a movable singularity with that of first integrals (when they exist). All these ideas are illustrated by a comprehensive analysis of the general two‐dimensional predator‐prey system.

Determinação de integrais primeiras liouvillianas em equações diferenciais ordinárias de segunda ordem

Nesta Tese desenvolvemos várias abordagens "Darbouxianas"para buscar integrais primeiras (elementares e Liouvillianas) de equações diferenciais ordinárias de segunda ordem (2EDOs) racionais. Os algoritmos (semi-algoritmos) que desenvolvemos seguem a linha do trabalho de Prelle e Singer. Basicamente, os métodos que buscam integrais primeiras elementares são uma extensão da técnica desenvolvida por Prelle e Singer para encontrar soluções elementares de equações diferenciais ordinárias de primeira ordem (1EDOs) racionais. O procedimento que lida com 2EDOs racionais que apresentam integrais primeiras Liouvillianas é baseado em uma extensão ao nosso método para encontrar soluções Liouvillianas de 1EDOs racionais. A ideia fundamental por tras do nosso trabalho consiste em que os fatores integrantes para 1-formas polinomiais geradas pela diferenciação de funções elementares e Liouvillianas são formados por certos polinômios denominados polinômios de Darboux. Vamos mostrar como combinar esses polinômios de Darboux para construir fatores integrantes e, de posse deles, determinar integrais primeiras. Vamos ainda discutir algumas implementações computacionais dos semi-algoritmos.

Classification analytique de systèmes différentiels linéaires déployant une singularité irrégulière de rang de Poincaré 1

Cette thèse traite de la classification analytique du déploiement de systèmes différentiels linéaires ayant une singularité irrégulière. Elle est composée de deux articles sur le sujet: le premier présente des résultats obtenus lors de l'étude de la confluence de l'équation hypergéométrique et peut être considéré comme un cas particulier du second; le deuxième contient les théorèmes et résultats principaux. Dans les deux articles, nous considérons la confluence de deux points singuliers réguliers en un point singulier irrégulier et nous étudions les conséquences de la divergence des solutions au point singulier irrégulier sur le comportement des solutions du système déployé. Pour ce faire, nous recouvrons un voisinage de l'origine (de manière ramifiée) dans l'espace du paramètre de déploiement $\epsilon$. La monodromie d'une base de solutions bien choisie est directement reliée aux matrices de Stokes déployées. Ces dernières donnent une interprétation géométrique aux matrices de Stokes, incluant le lien (existant au moins pour les cas génériques) entre la divergence des solutions à $\epsilon=0$ et la présence de solutions logarithmiques autour des points singuliers réguliers lors de la résonance. La monodromie d'intégrales premières de systèmes de Riccati correspondants est aussi interprétée en fonction des éléments des matrices de Stokes déployées. De plus, dans le second article, nous donnons le système complet d'invariants analytiques pour le déploiement de systèmes différentiels linéaires $x^2y'=A(x)y$ ayant une singularité irrégulière de rang de Poincaré $1$ à l'origine au-dessus d'un voisinage fixé $\mathbb{D}_r$ dans la variable $x$. Ce système est constitué d'une partie formelle, donnée par des polynômes, et d'une partie analytique, donnée par une classe d'équivalence de matrices de Stokes déployées. Pour chaque valeur du paramètre $\epsilon$ dans un secteur pointé à l'origine d'ouverture plus grande que $2\pi$, nous recouvrons l'espace de la variable, $\mathbb{D}_r$, avec deux secteurs et, au-dessus de chacun, nous choisissons une base de solutions du système déployé. Cette base sert à définir les matrices de Stokes déployées. Finalement, nous prouvons un théorème de réalisation des invariants qui satisfont une condition nécessaire et suffisante, identifiant ainsi l'ensemble des modules.

Global dynamics of the Rikitake system

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

On the global dynamics of the Rabinovich system

In this paper by using the Poincare compactification in R(3) make a global analysis of the Rabinovich system(x) over dot = hy - v(1)x + yz, (y) over dot = hx - v(2)y - xz, (z) over dot = -v(3)z + xy,with (x, y, z) is an element of R(3) and ( h, v(1), v(2), v(3)) is an element of R(4). We give the complete description of its dynamics on the sphere at infinity. For ten sets of the parameter values the system has either first integrals or invariants. For these ten sets we provide the global phase portrait of the Rabinovich system in the Poincare ball (i.e. in the compactification of R(3) with the sphere S(2) of the infinity). We prove that for convenient values of the parameters the system has two families of singularly degenerate heteroclinic cycles. Then changing slightly the parameters we numerically found a four wings butterfly shaped strange attractor.

Conjuntos minimais de sistemas lineares por partes

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

On symmetries of singular implicit ODEs

We study implicit ODEs, cubic in derivative, with infinitesimal symmetry at singular points. Cartan showed that even at regular points the existence of nontrivial symmetry imposes restrictions on the ODE. Namely, this algebra has the maximal possible dimension 3 iff the web of solutions is flat. For cubic ODEs with flat 3-web of solutions we establish sufficient conditions for the existence of nontrivial symmetries at singular points and show that under natural assumptions such a symmetry is semi-simple, i.e. is a scaling is some coordinates. We use this symmetry to find first integrals of the ODE.

Equações diferenciais implícitas com descontinuidade

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

On the chaotic rotation of a liquid-filled gyrostat via the Melnikov-Holmes-Marsden integral

The non-linear motions of a gyrostat with an axisymmetrical, fluid-filled cavity are investigated. The cavity is considered to be completely filled with an ideal incompressible liquid performing uniform rotational motion. Helmholtz theorem, Euler's angular momentum theorem and Poisson equations are used to develop the disturbed Hamiltonian equations of the motions of the liquid-filled gyrostat subjected to small perturbing moments. The equations are established in terms of a set of canonical variables comprised of Euler angles and the conjugate angular momenta in order to facilitate the application of the Melnikov-Holmes-Marsden (MHM) method to investigate homoclinic/heteroclinic transversal intersections. In such a way, a criterion for the onset of chaotic oscillations is formulated for liquid-filled gyrostats with ellipsoidal and torus-shaped cavities and the results are confirmed via numerical simulations. (c) 2006 Elsevier Ltd. All rights reserved.

First results for the Coulomb gauge integrals using NDIM

The Coulomb gauge has at least two advantages over other gauge choices in that bound states between quarks and studies of confinement are easier to understand in this gauge. However, perturbative calculations, namely Feynman loop integrations, are not well defined (there are the so-called energy integrals) even within the context of dimensional regularization. Leibbrandt and Williams proposed a possible cure to such a problem by splitting the space-time dimension into D = ω + ρ, i.e., introducing a specific parameter ρ to regulate the energy integrals. The aim of our work is to apply the negative dimensional integration method (NDIM) to the Coulomb gauge integrals using the recipe of split-dimension parameters and present complete results - finite and divergent parts - to the one- and two-loop level for arbitrary exponents of the propagators and dimension.

An example of cognitive obstacles in advanced integration: the case of scalar line integrals

Cognitive obstacles that arise in the teaching and learning of scalar line integrals, derived from cognitive aids provided to students when first learning about integration of single variable functions are described. A discussion of how and why the obstacles cause students problems is presented and possible strategies to overcome the obstacles are outlined.

Boundedness and convergence of singular integrals on fractal type sets

The topic of this dissertation lies in the intersection of harmonic analysis and fractal geometry. We particulary consider singular integrals in Euclidean spaces with respect to general measures, and we study how the geometric structure of the measures affects certain analytic properties of the operators. The thesis consists of three research articles and an overview. In the first article we construct singular integral operators on lower dimensional Sierpinski gaskets associated with homogeneous Calderón-Zygmund kernels. While these operators are bounded their principal values fail to exist almost everywhere. Conformal iterated function systems generate a broad range of fractal sets. In the second article we prove that many of these limit sets are porous in a very strong sense, by showing that they contain holes spread in every direction. In the following we connect these results with singular integrals. We exploit the fractal structure of these limit sets, in order to establish that singular integrals associated with very general kernels converge weakly. Boundedness questions consist a central topic of investigation in the theory of singular integrals. In the third article we study singular integrals of different measures. We prove a very general boundedness result in the case where the two underlying measures are separated by a Lipshitz graph. As a consequence we show that a certain weak convergence holds for a large class of singular integrals.

Beyond fractional derivatives: local approximation of other convolution integrals

Dynamic systems involving convolution integrals with decaying kernels, of which fractionally damped systems form a special case, are non-local in time and hence infinite dimensional. Straightforward numerical solution of such systems up to time t needs O(t(2)) computations owing to the repeated evaluation of integrals over intervals that grow like t. Finite-dimensional and local approximations are thus desirable. We present here an approximation method which first rewrites the evolution equation as a coupled in finite-dimensional system with no convolution, and then uses Galerkin approximation with finite elements to obtain linear, finite-dimensional, constant coefficient approximations for the convolution. This paper is a broad generalization, based on a new insight, of our prior work with fractional order derivatives (Singh & Chatterjee 2006 Nonlinear Dyn. 45, 183-206). In particular, the decaying kernels we can address are now generalized to the Laplace transforms of known functions; of these, the power law kernel of fractional order differentiation is a special case. The approximation can be refined easily. The local nature of the approximation allows numerical solution up to time t with O(t) computations. Examples with several different kernels show excellent performance. A key feature of our approach is that the dynamic system in which the convolution integral appears is itself approximated using another system, as distinct from numerically approximating just the solution for the given initial values; this allows non-standard uses of the approximation, e. g. in stability analyses.