MATHEMATICA notebook for computing tetrahedral finite element shape functions and matrices for the Helmholtz equation

A MATHEMATICA notebook to compute the elements of the matrices which arise in the solution of the Helmholtz equation by the finite element method (nodal approximation) for tetrahedral elements of any approximation order is presented. The results of the notebook enable a fast computational implementation of finite element codes for high order simplex 3D elements reducing the overheads due to implementation and test of the complex mathematical expressions obtained from the analytical integrations. These matrices can be used in a large number of applications related to physical phenomena described by the Poisson, Laplace and Schrodinger equations with anisotropic physical properties.

Edge detection and noise removal by use of a partial differential equation with automatic selection of parameters

This work deals with noise removal by the use of an edge preserving method whose parameters are automatically estimated, for any application, by simply providing information about the standard deviation noise level we wish to eliminate. The desired noiseless image u(x), in a Partial Differential Equation based model, can be viewed as the solution of an evolutionary differential equation u t(x) = F(u xx, u x, u, x, t) which means that the true solution will be reached when t ® ¥. In practical applications we should stop the time ''t'' at some moment during this evolutionary process. This work presents a sufficient condition, related to time t and to the standard deviation s of the noise we desire to remove, which gives a constant T such that u(x, T) is a good approximation of u(x). The approach here focused on edge preservation during the noise elimination process as its main characteristic. The balance between edge points and interior points is carried out by a function g which depends on the initial noisy image u(x, t0), the standard deviation of the noise we want to eliminate and a constant k. The k parameter estimation is also presented in this work therefore making, the proposed model automatic. The model's feasibility and the choice of the optimal time scale is evident through out the various experimental results.

A Decomposition and Noise Removal Method Combining Diffusion Equation and Wave Atoms for Textured Images

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

Solutions of the Fokker-Planck equation for a Morse isospectral potential

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

A General Solution of the Fokker-Planck Equation

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

Particular solution for anomalous diffusion equation with source term

In literature the phenomenon of diffusion has been widely studied, however for nonextensive systems which are governed by a nonlinear stochastic dynamic, there are a few soluble models. The purpose of this study is to present the solution of the nonlinear Fokker-Planck equation for a model of potential with barrier considering a term of absorption. Systems of this nature can be observed in various chemical or biological processes and their solution enriches the studies of existing nonextensive systems.

Discussion on the limit cycles of the Lev Ginzburg equation by M. Bellamy and R.E. Mickens, Journal of Sound and Vibration 308 (2007) 337-342

The authors M. Bellamy and R.E. Mickens in the article "Hopf bifurcation analysis of the Lev Ginzburg equation" published in Journal of Sound and Vibration 308 (2007) 337-342, claimed that this differential equation in the plane can exhibit a limit cycle. Here we prove that the Lev Ginzburg differential equation has no limit cycles. (C) 2012 Elsevier Ltd. All rights reserved.

Global dynamics of stationary solutions of the extended Fisher-Kolmogorov equation

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

On the symplectic two-form of gravity in terms of Dirac eigenvalues

The Dirac eigenvalues form a subset of observables of the Euclidean gravity. The symplectic two-form in the covariant phase space could be expressed, in principle, in terms of the Dirac eigenvalues. We discuss the existence of the formal solution of the equations defining the components of the symplectic form in this framework. (C) 2002 Published by Elsevier B.V. B.V.

Vortices in nonlocal Gross-Pitaevskii equation

We consider vortices in the nonlocal two-dimensional Gross-Pitaevskii equation with the interaction potential having Lorentz-shaped dependence on the relative momentum. It is shown that in the Fourier series expansion with respect to the polar angle, the unstable modes of the axial n-fold vortex have orbital numbers l satisfying 0 < \l\ < 2\n\, as in the local model. Numerical simulations show that nonlocality slightly decreases the threshold rotation frequency above which the nonvortex state ceases to be the global energy minimum and decreases the frequency of the anomalous mode of the 1-vortex. In the case of higher axial vortices, nonlocality leads to instability against splitting with the creation of antivortices and gives rise to additional anomalous modes with higher orbital numbers. Despite new instability channels with the creation of antivortices, for a stationary solution comprised of vortices and antivortices there always exists another vortex solution, composed solely of vortices, with the same total vorticity but with a lower energy.

Relaxation algorithm to hyperbolic states in Gross-Pitaevskii equation

A new version of the relaxation algorithm is proposed in order to obtain the stationary ground-state solutions of nonlinear Schrodinger-type equations, including the hyperbolic solutions. In a first example, the method is applied to the three-dimensional Gross-Pitaevskii equation, describing a condensed atomic system with attractive two-body interaction in a non-symmetrical trap, to obtain results for the unstable branch. Next, the approach is also shown to be very reliable and easy to be implemented in a non-symmetrical case that we have bifurcation, with nonlinear cubic and quintic terms. (c) 2006 Elsevier B.V. All rights reserved.

Quadratic gravity in (2+1)D with a topological Chern-Simons term

Three-dimensional quadratic gravity, unlike general relativity in (2+1)D, is dynamically nontrivial and has a well behaved nonrelativistic potential. Here we analyse the changes that occur when a topological Chem-Simons term is added to this theory. It is found that the harmless massive scalar mode of the latter gives rise to a troublesome massive spin-0 ghost, while the massive spin-2 ghost is replaced by two massive physical particles both of spin 2. We also found that light deflection does not have the 'wrong sign' such as in the framework of three-dimensional quadratic gravity.

Propagator, tree-level unitarity and effective nonrelativistic potential for higher-derivative gravity theories in D dimensions

A prescription for computing the propagator for D-dimensional higher-derivative gravity theories, based on the Barnes-Rivers operators, is presented. A systematic study of the tree-level unitarity of these theories is developed and the agreement of their linearized versions with Newton's law is investigated by computing the corresponding effective nonrelativistic potential. Three-dimensional quadratic gravity with a gravitational Chern-Simons term is also analyzed. A discussion on the issue of light bending within the framework of both D-dimensional quadratic gravity and three-dimensional quadratic gravity with a Chern-Simons term is provided as well. (C) 2002 American Institute of Physics.

Can one cure the nommitarity of three-dimensional quadratic gravity using a topological Chern-Simons term?

Quadratic gravity in (2+1)D is nonunitarity at the tree level. When a topological Chern-Simons term is added to this theory, the harmless massive scalar mode of the former gives rise to a troublesome massive spin-0 ghost, while the massive spin-2 ghost is replaced by two massive physical particles both of spin-2. Therefore, unlike what it is claimed in the literature, quadratic Chern-Simons gravity in (2+1)D is nonunitary at the tree level.

Whitham method for the Benjamin-Ono-Burgers equation and dispersive shocks

The Whitham modulation equations for the parameters of a periodic solution are derived using the generalized Lagrangian approach for the case of the damped Benjamin-Ono equation. The structure of the dispersive shock is considered in this method.