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.

BERRY PHASE, AHARONOV-BOHM EFFECT AND TOPOLOGY

The Berry phase for an electron in a one-dimensional box rotated around a magnetic flux line has contributions from the geometry and the magnetic flux, which gives an Aharonov-Bohm effect. For a circular box enclosing the magnetic flux, the Berry phase depends on the boundary conditions.

GROUND-STATE-ENERGY OF SINGULAR POTENTIALS

It is shown that for singular potentials of the form lambda/r(alpha),the asymptotic form of the wave function both at r --> infinity and r --> 0 plays an important role. Using a wave function having the correct asymptotic behavior for the potential lambda/r(4), it is, shown that it gives the exact ground-state energy for this potential when lambda --> 0, as given earlier by Harrell [Ann. Phys. (NY) 105, 379 (1977)]. For other values of the coupling parameter X, a trial basis;set of wave functions which also satisfy the correct boundary conditions at r --> infinity and r --> 0 are used to find the ground-state energy of the singular potential lambda/r(4) It is shown that the obtained eigenvalues are in excellent agreement with their exact ones for a very large range of lambda values.

Nonadiabatic calculations of the oscillator strengths for the helium atom in the hyperspherical adiabatic approach

Energies and wavefunctions are calculated for the bound states of the helium atom in the hyperspherical adiabatic approach by the full inclusion of nonadiabatic couplings. We show that the use of appropriate asymptotic radial boundary conditions not only allows the efficient calculation of energies accurate up to a few ppm for the ground state but also gives increasingly precise results for high-lying excited states with a unique set of equations. The accuracy of the wavefunctions is demonstrated by the calculation of oscillator strengths in the length form for transitions between stares ii S-1(e) and (n + 1) P-1(0) up to n = 29, in agreement with variational calculations.

Bethe ansatz solutions for Temperley-Lieb quantum spin chains

We solve the spectrum of quantum spin chains based on representations of the Temperley-Lieb algebra associated with the quantum groups U-q(X-n) for X-n = A(1), B-n, C-n and D-n. The tool is a modified version of the coordinate Bethe ansatz through a suitable choice of the Bethe states which give to all models the same status relative to their diagonalization. All these models have equivalent spectra up to degeneracies and the spectra of the lower-dimensional representations are contained in the higher-dimensional ones. Periodic boundary conditions, free boundary conditions and closed nonlocal boundary conditions are considered. Periodic boundary conditions, unlike free boundary conditions, bleak quantum group invariance. For closed nonlocal cases the models are quantum group invariant as well as periodic in a certain sense.

When the Casimir energy is not a sum of zero-point energies

We compute the leading radiative correction to the Casimir force between two parallel plates in the lambdaPhi(4) theory. Dirichlet and periodic boundary conditions are considered. A heuristic approach, in which the Casimir energy is computed as the sum of one-loop corrected zero-point energies, is shown to yield incorrect results, but we show how to amend it. The technique is then used in the case of periodic boundary conditions to construct a perturbative expansion which is free of infrared singularities in the massless limit. In this case we also compute the next-to-leading order radiative correction, which turns out to be proportional to lambda(3/2).

Dynamics in dumbbell domains I. Continuity of the set of equilibria

We analyze the dynamics of a reaction-diffusion equation with homogeneous Neumann boundary conditions in a dumbbell domain. We provide an appropriate functional setting to treat this problem and, as a first step, we show in this paper the continuity of the set of equilibria and of its linear unstable manifolds. (c) 2006 Elsevier B.V. All rights reserved.

An application of the Adomian's decomposition method to the analysis of MHD duct flows

A new approach is proposed in this work for the treatment of boundary value problems through the Adomian's decomposition method. Although frequently claimed as accurate and having fast convergence rates, the original formulation of Adomian's method does not allow the treatment of homogeneous boundary conditions along closed boundaries. The technique here presented overcomes this difficulty, and is applied to the analysis of magnetohydrodynamic duct flows. Results are in good agreement with finite element method calculations and analytical solutions for square ducts. Therefore, new possibilities appear for the application of Adomian's method in electromagnetics.

The Yang-Lee zeros of the 1D Blume-Capel model on connected and non-connected rings

We carry out a numerical and analytic analysis of the Yang-Lee zeros of the ID Blume-Capel model with periodic boundary conditions and its generalization on Feynman diagrams for which we include sums over all connected and nonconnected rings for a given number of spins. In both cases, for a specific range of the parameters, the zeros originally on the unit circle are shown to depart from it as we increase the temperature beyond some limit. The curve of zeros can bifurcate- and become two disjoint arcs as in the 2D case. We also show that in the thermodynamic limit the zeros of both Blume-Capel models on the static (connected ring) and on the dynamical (Feynman diagrams) lattice tend to overlap. In the special case of the 1D Ising model on Feynman diagrams we can prove for arbitrary number of spins that the Yang-Lee zeros must be on the unit circle. The proof is based on a property of the zeros of Legendre polynomials.

Dirac equation in magnetic-solenoid field

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

South Atlantic Ocean cyclogenesis climatology simulated by regional climate model (RegCM3)

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

Quantum rings of arbitrary shape and non-uniform width in a threading magnetic field

The electronic states of quantum rings with centerlines of arbitrary shape and non-uniform width in a threading magnetic field are calculated. The solutions of the Schrodinger equation with Dirichlet boundary conditions are obtained by a variational separation of variables in curvilinear coordinates. We obtain a width profile that compensates for the main effects of the curvature variations in the centerline. Numerical results are shown for circular, elliptical, and limacon-shaped quantum rings. We also show that smooth and tiny variations in the width may strongly affect the Aharonov-Bohm oscillations.

Fourth-order method for solving the Navier-Stokes equations in a constricting channel

A fourth-order numerical method for solving the Navier-Stokes equations in streamfunction/vorticity formulation on a two-dimensional non-uniform orthogonal grid has been tested on the fluid flow in a constricted symmetric channel. The family of grids is generated algebraically using a conformal transformation followed by a non-uniform stretching of the mesh cells in which the shape of the channel boundary can vary from a smooth constriction to one which one possesses a very sharp but smooth corner. The generality of the grids allows the use of long channels upstream and downstream as well as having a refined grid near the sharp corner. Derivatives in the governing equations are replaced by fourth-order central differences and the vorticity is eliminated, either before or after the discretization, to form a wide difference molecule for the streamfunction. Extra boundary conditions, necessary for wide-molecule methods, are supplied by a procedure proposed by Henshaw et al. The ensuing set of non-linear equations is solved using Newton iteration. Results have been obtained for Reynolds numbers up to 250 for three constrictions, the first being smooth, the second having a moderately sharp corner and the third with a very sharp corner. Estimates of the error incurred show that the results are very accurate and substantially better than those of the corresponding second-order method. The observed order of the method has been shown to be close to four, demonstrating that the method is genuinely fourth-order. © 1977 John Wiley & Sons, Ltd.