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.

Optimal and sub-optimal control in Dengue epidemics

This work concerns the application of the optimal control theory to Dengue epidemics. The dynamics of this insect-borne disease is modelled as a set of non-linear ordinary differential equations including the effect of educational campaigns organized to motivate the population to break the reproduction cycle of the mosquitoes by avoiding the accumulation of still water in open-air recipients. The cost functional is such that it reflects a compromise between actual financial spending (in insecticides and educational campaigns) and the population health (which can be objectively measured in terms of, for instance, treatment costs and loss of productivity). The optimal control problem is solved numerically using a multiple shooting method. However, the optimal control policy is difficult to implement by the health authorities because it is not practical to adjust the investment rate continuously in time. Therefore, a suboptimal control policy is computed assuming, as the admissible set, only those controls which are piecewise constant. The performance achieved by the optimal control and the sub-optimal control policies are compared with the cases of control using only insecticides when Breteau Index is greater or equal to 5 and the case of no-control. The results show that the sub-optimal policy yields a substantial reduction in the cost, in terms of the proposed functional, and is only slightly inferior to the optimal control policy. Copyright (C) 2001 John Wiley & Sons, Ltd.

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.

Euclidean 4d exact solitons in a Skyrme type model

We introduce a Skyrme type, four-dimensional Euclidean field theory made of a triplet of scalar fields n→, taking values on the sphere S2, and an additional real scalar field φ, which is dynamical only on a three-dimensional surface embedded in R4. Using a special ansatz we reduce the 4d non-linear equations of motion into linear ordinary differential equations, which lead to the construction of an infinite number of exact soliton solutions with vanishing Euclidean action. The theory possesses a mass scale which fixes the size of the solitons in way which differs from Derrick's scaling arguments. The model may be relevant to the study of the low energy limit of pure SU(2) Yang-Mills theory. © 2004 Elsevier B.V. All rights reserved.

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.

Limit cycles for singular perturbation problems via inverse integrating factor

In this paper singularly perturbed vector fields Xε defined in ℝ2 are discussed. The main results use the solutions of the linear partial differential equation XεV = div(Xε)V to give conditions for the existence of limit cycles converging to a singular orbit with respect to the Hausdorff distance. © SPM.

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.

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.

A semigroups theory approach to a model of suspension bridges

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

Utilização de equações diferenciais parciais no tratamento de imagens orbitais

Pós-graduação em Ciências Cartográficas - FCT

Modelagem matemática da doença do caranguejo letárgico via ondas viajantes

Pós-graduação em Biometria - IBB

Simulação numérica de escoamentos de fluidos pelo método de elementos finitos de mínimos quadrados

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

Controle exato para a equação de onda: o método HUM

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

Leis de Conservação Hiperbólicas 2D com Termo Fonte Stiff

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