Homoclinic bifurcations in reversible Hamiltonian systems

We study the existence of homoclic solutions for reversible Hamiltonian systems taking the family of differential equations u(iv) + au - u +f(u, b) = 0 as a model, where fis an analytic function and a, b real parameters. These equations are important in several physical situations such as solitons and in the existence of finite energy stationary states of partial differential equations, but no assumptions of any kind of discrete symmetry is made and the analysis here developed can be extended to others Hamiltonian systems and successfully employed in situations where standard methods fail. We reduce the problem of computing these orbits to that of finding the intersection of the unstable manifold with a suitable set and then apply it to concrete situations. We also plot the homoclinic values configuration in parameters space, giving a picture of the structural distribution and a geometrical view of homoclinic bifurcations. (c) 2005 Published by Elsevier B.V.

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.

Correlation times in stochastic equations with delayed feedback and multiplicative noise

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

Non-local interactions and the dynamics of dispersal in immature insects

A simple mathematical model is developed to explain the appearance of oscillations in the dispersal of larvae from the food source in experimental populations of certain species of blowflies. The life history of the immature stage in these flies, and in a number of other insects, is a system with two populations, one of larvae dispersing on the soil and the other of larvae that burrow in the soil to pupate. The observed oscillations in the horizontal distribution of buried pupae at the end of the dispersal process are hypothesized to be a consequence of larval crowding at a given point in the pupation substrate. It is assumed that dispersing larvae are capable of perceiving variations in density of larvae buried at a given point in the substrate of pupation, and that pupal density may influence pupation of dispersing larvae. The assumed interaction between dispersing larvae and the larvae that are burrowing to pupate is modeled using the concept of non-local effects. Numerical solutions of integro-partial differential equations developed to model density-dependent immature dispersal demonstrate that variation in the parameter that governs the non-local interaction between dispersing and buried larvae induces oscillations in the final horizontal distribution of pupae. (C) 1997 Academic Press Limited.

The application of interpolating MLS approximations to the analysis of MHD flows

The element-free Galerkin method (EFGM) is a very attractive technique for solutions of partial differential equations, since it makes use of nodal point configurations which do not require a mesh. Therefore, it differs from FEM-like approaches by avoiding the need of meshing, a very demanding task for complicated geometry problems. However, the imposition of boundary conditions is not straightforward, since the EFGM is based on moving-least-squares (MLS) approximations which are not necessarily interpolants. This feature requires, for instance, the introduction of modified functionals with additional unknown parameters such as Lagrange multipliers, a serious drawback which leads to poor conditionings of the matrix equations. In this paper, an interpolatory formulation for MLS approximants is presented: it allows the direct introduction of boundary conditions, reducing the processing time and improving the condition numbers. The formulation is applied to the study of two-dimensional magnetohydrodynamic flow problems, and the computed results confirm the accuracy and correctness of the proposed formulation. (C) 2002 Elsevier B.V. All rights reserved.

On some classes of exactly-solvable Klein-Gordon equations

In this work we discuss some exactly solvable Klein-Gordon equations. We basically discuss the existence of classes of potentials with different nonrelativistic limits, but which shares the intermediate effective Schroedinger differential equation. We comment about the possible use of relativistic exact solutions as approximations for nonrelativistic inexact potentials. (c) 2005 Elsevier B.V. All rights reserved.

Sliding Vector Fields via Slow-Fast Systems

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

Analytical optimization of spread and change of exponential decay of generalized Wannier functions in one dimension

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

Object oriented modeling of piston engines

This work reports a conception phase of a piston engine global model. The model objective is forecast the motor performance (power, torque and specific consumption as a function of rotation and environmental conditions). Global model or Zero-dimensional is based on flux balance through each engine component. The resulting differential equations represents a compressive unsteady flow, in which, all dimensional variables are areas or volumes. A review is presented first. The ordinary differential equation system is presented and a Runge-Kutta method is proposed to solve it numerically. The model includes the momentum conservation equation to link the gas dynamics with the engine moving parts rigid body mechanics. As an oriented to objects model the documentation follows the UML standard. A discussion about the class diagrams is presented, relating the classes with physical model related. The OOP approach allows evolution from simple models to most complex ones without total code rewrite. Copyright © 2001 Society of Automotive Engineers, Inc.

Simulations of incompressible fluid flows by a least squares finite element method

In this work simulations of incompressible fluid flows have been done by a Least Squares Finite Element Method (LSFEM) using velocity-pressure-vorticity and velocity-pressure-stress formulations, named u-p-ω) and u-p-τ formulations respectively. These formulations are preferred because the resulting equations are partial differential equations of first order, which is convenient for implementation by LSFEM. The main purposes of this work are the numerical computation of laminar, transitional and turbulent fluid flows through the application of large eddy simulation (LES) methodology using the LSFEM. The Navier-Stokes equations in u-p-ω and u-p-τ formulations are filtered and the eddy viscosity model of Smagorinsky is used for modeling the sub-grid-scale stresses. Some benchmark problems are solved for validate the numerical code and the preliminary results are presented and compared with available results from the literature. Copyright © 2005 by ABCM.

Integrating open source codes for solid modeling and automatic mesh generation in FEM applications

The applications of the Finite Element Method (FEM) for three-dimensional domains are already well documented in the framework of Computational Electromagnetics. However, despite the power and reliability of this technique for solving partial differential equations, there are only a few examples of open source codes available and dedicated to the solid modeling and automatic constrained tetrahedralization, which are the most time consuming steps in a typical three-dimensional FEM simulation. Besides, these open source codes are usually developed separately by distinct software teams, and even under conflicting specifications. In this paper, we describe an experiment of open source code integration for solid modeling and automatic mesh generation. The integration strategy and techniques are discussed, and examples and performance results are given, specially for complicated and irregular volumes which are not simply connected. © 2011 IEEE.

On nonlinear dynamics in a free fluid surface of a tank excited by a non-ideal power source

We investigate the nonlinear oscillations in a free surface of a fluid in a cylinder tank excited by non-ideal power source, an electric motor with limited power supply. We study the possibility of parametric resonance in this system, showing that the excitation mechanism can generate chaotic response. Additionally, the dynamics of parametrically excited surface waves in the tank can reveal new characteristics of the system. The fluid-dynamic system is modeled in such way as to obtain a nonlinear differential equation system. Numerical experiments are carried out to find the regions of chaotic solutions. Simulation results are presented as phase-portrait diagrams characterizing the resonant vibrations of free fluid surface and the existence of several types of regular and chaotic attractors. We also describe the energy transfer in the interaction process between the hydrodynamic system and the electric motor. Copyright © 2011 by ASME.

Forced oscillations with continuum models of atomic force microscopy

The dynamics of the AFM-atomic force microscope follows a model based in a Timoshenko cantilever beam with a tip attached at the free end and acting with the surface of a sample. General boundary conditions arise when the tip is either in contact or non-contact with the surface. The governing equations are given in matrix conservative form subject to localized loads. The eigenanalysis is done with a fundamental matrix response of a damped second-order matrix differential equation. Forced responses are found by using a Galerkin approximation of the matrix impulse response. Simulations results with harmonic and pulse forcing show the filtering character and the effects of the tip-sample interaction at the end of the beam. © 2012 American Institute of Physics.

Alternate treatments of jacobian singularities in polar coordinates within finite-difference schemes

Jacobian singularities of differential operators in curvilinear coordinates occur when the Jacobian determinant of the curvilinear-to-Cartesian mapping vanishes, thus leading to unbounded coefficients in partial differential equations. Within a finite-difference scheme, we treat the singularity at the pole of polar coordinates by setting up complementary equations. Such equations are obtained by either integral or smoothness conditions. They are assessed by application to analytically solvable steady-state heat-conduction problems.

On certain new exact solutions of a diffusive predator-prey system

We construct exact solutions for a system of two coupled nonlinear partial differential equations describing the spatio-temporal dynamics of a predator-prey system where the prey per capita growth rate is subject to the Allee effect. Using the G'/G expansion method, we derive exact solutions to this model for two different wave speeds. For each wave velocity we report three different forms of solutions. We also discuss the biological relevance of the solutions obtained. © 2012 Elsevier B.V.