Analytical calculation of the transition to complete phase synchronization in coupled oscillators

Here we present a system of coupled phase oscillators with nearest neighbors coupling, which we study for different boundary conditions. We concentrate at the transition to the total synchronization. We are able to develop exact solutions for the value of the coupling parameter when the system becomes completely synchronized, for the case of periodic boundary conditions as well as for a chain with fixed ends. We compare the results with those calculated numerically.

Determination of the critical coupling for oscillators in a ring

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

The canonical structure of Podolsky's generalized electrodynamics on the null-plane

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

Numerical simulation of a radial diffuser turbulent airflow

In the present work are presented results from numerical simulations performed with the ANSYS-CFX (R) code. We have studied a radial diffuser flow case, which is the main academic problem used to study the flow behavior on flat plate valves. The radial flow inside the diffuser has important behavior such as the turbulence decay downstream and recirculation regions inside the valve flow channel due to boundary layer detachment. These flow structures are present in compressor reed valve configurations, influencing to a greater extent the compressor efficiency. The main target of the present paper was finding the simulation set-up (computational domain, boundary conditions and turbulence model) that better fits with experimental data published by Tabatabai and Pollard. The local flow turbulence and velocity profiles were investigated using four different turbulence models, two different boundary conditions set-up, two different computational domains and three different flow conditions (Re-in - Reynolds number at the diffuser inlet). We used the Reynolds stress (BSL); the k-epsilon; the RNG k-epsilon; and the shear stress transport (SST) k-omega turbulence models. The performed analysis and comparison of the computational results with experimental data show that the choice of the turbulence model, as well as the choice of the other computational conditions, plays an important role in the results physical quality and accuracy. (c) 2007 Elsevier B.V. All rights reserved.

Energy spectrum and eigenfunctions through the Quantum Section Method

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

PULLBACK ATTRACTORS FOR A SINGULARLY NONAUTONOMOUS PLATE EQUATION

We consider the family of singularly nonautonomous plate equation with structural dampingu(tt) + a(t, x)u(t) - Delta u(t) + (-Delta)(2)(u) + lambda u = f(u),in a bounded domain Omega subset of R(n), with Navier boundary conditions. When the nonlinearity f is dissipative we show that this problem is globally well posed in H(0)(2)(Omega) x L(2)(Omega) and has a family of pullback attractors which is upper-semicontinuous under small perturbations of the damping a.

Lower bounds on blow up solutions of the three-dimensional Navier-Stokes equations in homogeneous Sobolev spaces

Suppose that u(t) is a solution of the three-dimensional Navier-Stokes equations, either on the whole space or with periodic boundary conditions, that has a singularity at time T. In this paper we show that the norm of u(T - t) in the homogeneous Sobolev space (H)over dot(s) must be bounded below by c(s)t(-(2s-1)/4) for 1/2 < s < 5/2 (s not equal 3/2), where c(s) is an absolute constant depending only on s; and by c(s)parallel to u(0)parallel to((5-2s)/5)(L2)t(-2s/5) for s > 5/2. (The result for 1/2 < s < 3/2 follows from well-known lower bounds on blowup in Lp spaces.) We show in particular that the local existence time in (H)over dot(s)(R-3) depends only on the (H)over dot(s)-norm for 1/2 < s < 5/2, s not equal 3/2. (C) 2012 American Institute of Physics. [http://dx.doi.org/10.1063/1.4762841]

Effects of varying curvature and width on the electronic states of GaAs quantum rings

Stationary states of an electron in thin GaAs elliptical quantum rings are calculated within the effective-mass approximation. The width of the ring varies smoothly along the centerline, which is an ellipse. The solutions of the Schrödinger equation with Dirichlet boundary conditions are approximated by a product of longitudinal and transversal wave functions. The ground-state probability density shows peaks: (i) where the curvature is larger in a constant-with ring, and (ii) in thicker parts of a circular ring. For rings of typical dimensions, it is shown that the effects of a varying width may be stronger than those of the varying curvature. Also, a width profile which compensates the main localization effects of the varying curvature is obtained.

Interpolating EFGM for computing continuous and discontinuous electromagnetic fields

Purpose - This paper proposes an interpolating approach of the element-free Galerkin method (EFGM) coupled with a modified truncation scheme for solving Poisson's boundary value problems in domains involving material non-homogeneities. The suitability and efficiency of the proposed implementation are evaluated for a given set of test cases of electrostatic field in domains involving different material interfaces.Design/methodology/approach - the authors combined an interpolating approximation with a modified domain truncation scheme, which avoids additional techniques for enforcing the Dirichlet boundary conditions and for dealing with material interfaces usually employed in meshfree formulations.Findings - the local electric potential and field distributions were correctly described as well as the global quantities like the total potency and resistance. Since, the treatment of the material interfaces becomes practically the same for both the finite element method (FEM) and the proposed EFGM, FEM-oriented programs can, thus, be easily extended to provide EFGM approximations.Research limitations/implications - the robustness of the proposed formulation became evident from the error analyses of the local and global variables, including in the case of high-material discontinuity.Practical implications - the proposed approach has shown to be as robust as linear FEM. Thus, it becomes an attractive alternative, also because it avoids the use of additional techniques to deal with boundary/interface conditions commonly employed in meshfree formulations.Originality/value - This paper reintroduces the domain truncation in the EFGM context, but by using a set of interpolating shape functions the authors avoided the use of Lagrange multipliers as well Mathematics in Engineering high-material discontinuity.

Application of a meshless method in electromagnetics

An improved meshless method is presented with an emphasis on the detailed description of this new computational technique and its numerical implementations by investigating the usefulness of a commonly neglected parameter in this paper. Two approaches to enforce essential boundary conditions are also thoroughly investigated. Numerical tests on a mathematical function is carried out as a means of validating the proposed method. It will be seen that the proposed method is more robust than the conventional ones. Applications in solving electromagnetic problems are also presented.

NONLINEAR DIFFUSION PROCESS IN A BENARD SYSTEM AT THE CRITICAL-POINT FOR THE ONSET OF CONVECTION

The evolution equation governing surface perturbations of a shallow fluid heated from below at the critical Rayleigh number for the onset of convective motion, and with boundary conditions leading to zero critical wave number, is obtained. A solution for negative or cooling perturbations is explicitly exhibited, which shows that the system presents sharp propagating fronts.

The electrochemical behavior of pyrite-pyrrhotite mixtures

To assess the response of common sulfide minerals to oxidizing conditions, a methodology to immobilize mechanically solid particles on carbon surfaces (voltammetry of microparticles, VMP) was employed, to define the influence of the pyrrhotite content in pyrite-pyrrhotite mixtures. The influence of the galvanic interactions and local pH on the oxidation reaction of pyrite was also investigated. With this purpose, artificial two-mineral electrodes were constructed, ranging in weight from 20 to 80% pyrrhotite. The resulting cyclic voltammograms were analyzed and relative quantities of oxidation products were evaluated. The goal of this work was to define the boundary conditions, in terms of pyrrhotite content in the mixture, that determine the SO42-/S ratio obtained and to describe some parameters which influence this ratio: local pH and galvanic interactions. (C) 2003 Elsevier B.V. All rights reserved.

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.

A MODEL OF THE DEUTERON BASED ON A BREIT TYPE EQUATION

A relativistic treatment of the deuteron and its observables based on a two-body Dirac (Breit) equation, with phenomenological interactions, associated to one-boson exchanges with cutoff masses, is presented. The 16-component wave function for the deuteron (J(pi) = 1+) solution contains four independent radial functions which obey a system of four coupled differential equations of first order. This radial system is numerically integrated, from infinity to the origin, by fixing the value of the deuteron binding energy and using appropriate boundary conditions at infinity. Specific examples of mixtures containing scalar, pseudoscalar and vector like terms are discussed in some detail and several observables of the deuteron are calculated. Our treatment differs from more conventional ones in that nonrelativistic reductions of the order c-2 are not used.

Self-consistent equilibrium calculation through a direct variational technique in tokamak plasmas

A self-consistent equilibrium calculation, valid for arbitrary aspect ratio tokamaks, is obtained through a direct variational technique that reduces the equilibrium solution, in general obtained from the 2D Grad-Shafranov equation, to a 1D problem in the radial flux coordinate rho. The plasma current profile is supposed to have contributions of the diamagnetic, Pfirsch-Schluter and the neoclassical ohmic and bootstrap currents. An iterative procedure is introduced into our code until the flux surface averaged toroidal current density (J(T)), converges to within a specified tolerance for a given pressure profile and prescribed boundary conditions. The convergence criterion is applied between the (J(T)) profile used to calculate the equilibrium through the variational procedure and the one that results from the equilibrium and given by the sum of all current components. The ohmic contribution is calculated from the neoclassical conductivity and from the self-consistently determined loop voltage in order to give the prescribed value of the total plasma current. The bootstrap current is estimated through the full matrix Hirshman-Sigmar model with the viscosity coefficients as proposed by Shaing, which are valid in all plasma collisionality regimes and arbitrary aspect ratios. The results of the self-consistent calculation are presented for the low aspect ratio tokamak Experimento Tokamak Esferico. A comparison among different models for the bootstrap current estimate is also performed and their possible Limitations to the self-consistent calculation is analysed.