LEAST ACTION PRINCIPLE and THE INCOMPRESSIBLE EULER EQUATIONS WITH VARIABLE DENSITY

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

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.

Classes of exact Klein-Gordon equations with spatially dependent masses: Regularizing the one-dimensional inversely linear potential

In this work we consider the effect of a spatially dependent mass over the solution of the Klein-Gordon equation in 1 + 1 dimensions, particularly the case of inversely linear scalar potential, which usually presents problems of divergence of the ground-state wave function at the origin, and possible nonexistence of the even-parity wave functions. Here we study this problem, showing that for a certain dependence of the mass with respect to the coordinate, this problem disappears. (c) 2006 Elsevier B.V. All rights reserved.

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.

Difference equations with delays depending on time

For linear difference equations with coefficients and delays varying in time, sufficient conditions are given, in the scalar case, the zero solution to be stable. © 1990 Sociedade Brasileira de Matemática.

Analytical/numerical hybrid solution for one-dimensional ablation problem

Ablation is a thermal protection process with several applications in engineering, mainly in the field of airspace industry. The use of conventional materials must be quite restricted, because they would suffer catastrophic flaws due to thermal degradation of their structures. However, the same materials can be quite suitable once being protected by well-known ablative materials. The process that involves the ablative phenomena is complex, could involve the whole or partial loss of material that is sacrificed for absorption of energy. The analysis of the ablative process in a blunt body with revolution geometry will be made on the stagnation point area that can be simplified as a one-dimensional plane plate problem, hi this work the Generalized Integral Transform Technique (GITT) is employed for the solution of the non-linear system of coupled partial differential equations that model the phenomena. The solution of the problem is obtained by transforming the non-linear partial differential equation system to a system of coupled first order ordinary differential equations and then solving it by using well-established numerical routines. The results of interest such as the temperature field, the depth and the rate of removal of the ablative material are presented and compared with those ones available in the open literature.

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.

Detecção de bordas em imagens digitais através de um processo de difusão anisotrópica não linear

The edges detection model by a non-linear anisotropic diffusion, consists in a mathematical model of smoothing based in Partial Differential Equation (PDE), alternative to the conventional low-pass filters. The smoothing model consists in a selective process, where homogeneous areas of the image are smoothed intensely in agreement with the temporal evolution applied to the model. The level of smoothing is related with the amount of undesired information contained in the image, i.e., the model is directly related with the optimal level of smoothing, eliminating the undesired information and keeping selectively the interest features for Cartography area. The model is primordial for cartographic applications, its function is to realize the image preprocessing without losing edges and other important details on the image, mainly airports tracks and paved roads. Experiments carried out with digital images showed that the methodology allows to obtain the features, e.g. airports tracks, with efficiency.

Differential mixer with NMOS/PMOS stack at switching stage

This paper proposes a novel differential mixer topology. The traditional stage of switching is replaced by a stack of NMOS and PMOS transistors combined. A design is given of a 900 MHz down-conversion mixer using a 0.35 μm CMOS process. Comparison with conventional mixer shows that the topology leads to a better performance in terms of conversion gain and linearity. ©2012 IEEE.