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.

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.

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.

An alternative formulation for guided electromagnetic fields in grounded chiral slabs - Summary

An alternative formulation for guided electromagnetic fields in grounded chiral slabs is presented. This formulation is formally equivalent to the double Fourier transform method used by the authors to calculate the spectral fields in open chirostrip structures. In this paper, we have addressed the behavior of the electromagnetic fields in the vicinity of the ground plane and at the interface between the chiral substrate and the free space region. It was found that the boundary conditions for the magnetic field, valid for achiral media, are not completely satisfied when we deal with chiral material. Effects of chirality on electromagnetic field distributions and on surface wave dispersion curves were also analyzed.

Vortex lattice and matching fields for a long superconducting wire

We investigate the flux penetration patterns and matching fields of a long cylindrical wire of circular cross section in the presence of an external magnetic field. For this study we write the London theory for a long cylinder both for the mixed and Meissner states, with boundary conditions appropriate for this geometry. Using the Monte Carlo simulated annealing method, the free energy of the mixed state is minimized with respect to the vortex position and we obtain the ground state of the vortex lattice for N=3 up to 18 vortices. The free energy of the Meissner and mixed states provides expressions for the matching fields. We find that, as in the case of samples of different geometry, the finite-size effect provokes a delay on the vortex penetration and a vortex accumulation in the center of the sample. The vortex patterns obtained are in good agreement with experimental results.

Two-body Dirac equation approach to the deuteron

The two-body Dirac(Breit) equation with potentials associated to one-boson-exchanges with cutoff masses is solved for the deuteron and its observables calculated. The 16-component wave-function for the Jπ = 1+ state contains four independent radial functions which satisfy a system of four coupled differential equations of first order. This system is numerically integrated, from infinity towards the origin, by fixing the value of the deuteron binding energy and imposing appropriate boundary conditions at infinity. For the exchange potential of the pion, a mixture of direct plus derivative couplings to the nucleon is considered. We varied the pion-nucleon coupling constant, and the best results of our calculations agree with the lower values recently determined for this constant.

Size effects in the magnetization of a superconducting wire

The size effects in the magnetization of a long cylindrical wire of circular cross section in the presence of an external magnetic field are investigated. For this study the London theory is used with boundary conditions appropriate for this geometry. Using the Monte Carlo simulated annealing method, the free energy of the mixed state is minimized with respect to the vortex positions. The ground state of the vortex lattice for n = 1 up to 18 vortices for a given radius of the cylinder is obtained. It is found that the finite size of the sample provokes a matching effect in the magnetization, as found in experiments with superconducting samples of finite size but different geometry. © 1999 American Institute of Physics.

Dynamical calculation of direct muon-transfer rates from thermalized muonic hydrogen to C6+ and O8+

Moun-transfer reactions from muonic hydrogen to carbon and oxygen nuclei employing a full quantum-mechanical few-body description of rearrangement scattering were studied by solving the Faddeev-Hahn-type equations using close-coupling approximation. The application of a close-coupling-type ansatz led to satisfactory results for direct muon-transfer reactions from muonic hydrogen to C6+ and O8+.