VELOCITY FIELD OF COMPRESSIBLE MAGNETOHYDRODYNAMIC TURBULENCE: WAVELET DECOMPOSITION AND MODE SCALINGS

We study compressible magnetohydrodynamic turbulence, which holds the key to many astrophysical processes, including star formation and cosmic-ray propagation. To account for the variations of the magnetic field in the strongly turbulent fluid, we use wavelet decomposition of the turbulent velocity field into Alfven, slow, and fast modes, which presents an extension of the Cho & Lazarian decomposition approach based on Fourier transforms. The wavelets allow us to follow the variations of the local direction of the magnetic field and therefore improve the quality of the decomposition compared to the Fourier transforms, which are done in the mean field reference frame. For each resulting component, we calculate the spectra and two-point statistics such as longitudinal and transverse structure functions as well as higher order intermittency statistics. In addition, we perform a Helmholtz-Hodge decomposition of the velocity field into incompressible and compressible parts and analyze these components. We find that the turbulence intermittency is different for different components, and we show that the intermittency statistics depend on whether the phenomenon was studied in the global reference frame related to the mean magnetic field or in the frame defined by the local magnetic field. The dependencies of the measures we obtained are different for different components of the velocity; for instance, we show that while the Alfven mode intermittency changes marginally with the Mach number, the intermittency of the fast mode is substantially affected by the change.

A geometric mass-preserving redistancing scheme for the level set function

In this paper we describe and evaluate a geometric mass-preserving redistancing procedure for the level set function on general structured grids. The proposed algorithm is adapted from a recent finite element-based method and preserves the mass by means of a localized mass correction. A salient feature of the scheme is the absence of adjustable parameters. The algorithm is tested in two and three spatial dimensions and compared with the widely used partial differential equation (PDE)-based redistancing method using structured Cartesian grids. Through the use of quantitative error measures of interest in level set methods, we show that the overall performance of the proposed geometric procedure is better than PDE-based reinitialization schemes, since it is more robust with comparable accuracy. We also show that the algorithm is well-suited for the highly stretched curvilinear grids used in CFD simulations. Copyright (C) 2010 John Wiley & Sons, Ltd.

Numerical study of the atmospheric flow over a coastal cliff

This work presents a numerical method suitable for the study of the development of internal boundary layers (IBL) and their characteristics for flows over various types of coastal cliffs. The IBL is an important meteorological occurrence for flows with surface roughness and topographical step changes. A two-dimensional flow program was used for this study. The governing equations were written using the vorticity-velocity formulation. The spatial derivatives were discretized by high-order compact finite differences schemes. The time integration was performed with a low storage fourth-order Runge-Kutta scheme. The coastal cliff (step) was specified through an immersed boundary method. The validation of the code was done by comparison of the results with experimental and observational data. The numerical simulations were carried out for different coastal cliff heights and inclinations. The results show that the predominant factors for the height of the IBL and its characteristics are the upstream velocity, and the height and form (inclination) of the coastal cliff. Copyright (C) 2010 John Wiley & Sons, Ltd.

Fractal structures in nonlinear plasma physics

Fractal structures appear in many situations related to the dynamics of conservative as well as dissipative dynamical systems, being a manifestation of chaotic behaviour. In open area-preserving discrete dynamical systems we can find fractal structures in the form of fractal boundaries, associated to escape basins, and even possessing the more general property of Wada. Such systems appear in certain applications in plasma physics, like the magnetic field line behaviour in tokamaks with ergodic limiters. The main purpose of this paper is to show how such fractal structures have observable consequences in terms of the transport properties in the plasma edge of tokamaks, some of which have been experimentally verified. We emphasize the role of the fractal structures in the understanding of mesoscale phenomena in plasmas, such as electromagnetic turbulence.

Robust tori in a double-waved Hamiltonian model

A Hamiltonian system perturbed by two waves with particular wave numbers can present robust tori, which are barriers created by the vanishing of the perturbed Hamiltonian at some defined positions. When robust tori exist, any trajectory in phase space passing close to them is blocked by emergent invariant curves that prevent the chaotic transport. Our results indicate that the considered particular solution for the two waves Hamiltonian model shows plenty of robust tori blocking radial transport. (C) 2010 Elsevier B.V. All rights reserved.

Upgrading a TCABR data analysis and acquisition system for remote participation using Java, XML, RCP and modern client/server communication/authentication

The TCABR data analysis and acquisition system has been upgraded to support a joint research programme using remote participation technologies. The architecture of the new system uses Java language as programming environment. Since application parameters and hardware in a joint experiment are complex with a large variability of components, requirements and specification solutions need to be flexible and modular, independent from operating system and computer architecture. To describe and organize the information on all the components and the connections among them, systems are developed using the extensible Markup Language (XML) technology. The communication between clients and servers uses remote procedure call (RPC) based on the XML (RPC-XML technology). The integration among Java language, XML and RPC-XML technologies allows to develop easily a standard data and communication access layer between users and laboratories using common software libraries and Web application. The libraries allow data retrieval using the same methods for all user laboratories in the joint collaboration, and the Web application allows a simple graphical user interface (GUI) access. The TCABR tokamak team in collaboration with the IPFN (Instituto de Plasmas e Fusao Nuclear, Instituto Superior Tecnico, Universidade Tecnica de Lisboa) is implementing this remote participation technologies. The first version was tested at the Joint Experiment on TCABR (TCABRJE), a Host Laboratory Experiment, organized in cooperation with the IAEA (International Atomic Energy Agency) in the framework of the IAEA Coordinated Research Project (CRP) on ""Joint Research Using Small Tokamaks"". (C) 2010 Elsevier B.V. All rights reserved.

The non-twist standard map with robust tori

The non-twist standard map occurs frequently in many fields of science specially in modelling the dynamics of the magnetic field lines in tokamaks. Robust tori, dynamical barriers that impede the radial transport among different regions of the phase space, are introduced in the non-twist standard map in a conservative fashion. The resulting non-twist standard map with robust tori is an improved model to study transport barriers in plasmas confined in tokamaks.

Integrable maps with non-trivial topology: application to divertor configurations

We explore a method for constructing two-dimensional area-preserving, integrable maps associated with Hamiltonian systems, with a given set of fixed points and given invariant curves. The method is used to find an integrable Poincare map for the field lines in a large aspect ratio tokamak with a poloidal single-null divertor. The divertor field is a superposition of a magnetohydrodynamic equilibrium with an arbitrarily chosen safety factor profile, with a wire carrying an electric current to create an X-point. This integrable map is perturbed by an impulsive perturbation that describes non-axisymmetric magnetic resonances at the plasma edge. The non-integrable perturbed map is applied to study the structure of the open field lines in the scrape-off layer, reproducing the main transport features obtained by integrating numerically the magnetic field line equations, such as the connection lengths and magnetic footprints on the divertor plate.

Temporal behaviour of toroidal rotation velocity in the TCABR tokamak

A new method for determining the temporal evolution of plasma rotation is reported in this work. The method is based upon the detection of two different portions of the spectral profile of a plasma impurity line, using a monochromator with two photomultipliers installed at the exit slits. The plasma rotation velocity is determined by the ratio of the two detected signals. The measured toroidal rotation velocities of C III (4647.4 angstrom) and C VI (5290.6 angstrom), at different radial positions in TCABR discharges, show good agreement, within experimental uncertainty, with previous results (Severo et al 2003 Nucl. Fusion 43 1047). In particular, they confirm that the plasma core rotates in the direction opposite to the plasma current, while near the plasma edge (r/a > 0.9) the rotation is in the same direction. This technique was also used to investigate the dependence of toroidal rotation on the poloidal position of gas puffing. The results show that there is no dependence for the plasma core, while for plasma edge (r/a > 0.9) some dependence is observed.

Results of Joint Experiments and other IAEA activities on research using small tokamaks

This paper presents an overview of the results obtained during the Joint Experiments organized in the framework of the IAEA Coordinated Research Project on `Joint Research Using Small Tokamaks` that have been carried out on the tokamaks CASTOR at IPP Prague, Czech Republic (2005), T-10 at RRC `Kurchatov Institute`, Moscow, Russia (2006), and the most recent one at ISTTOK at IST, Lisbon, Portugal, in 2007. Experimental programmes were aimed at diagnosing and characterizing the core and the edge plasma turbulence in a tokamak in order to investigate correlations between the occurrence of transport barriers, improved confinement, electric fields and electrostatic turbulence using advanced diagnostics with high spatial and temporal resolution. On CASTOR and ISTTOK, electric fields were generated by biasing an electrode inserted into the edge plasma and an improvement of the global particle confinement induced by the electrode positive biasing has been observed. Geodesic acoustic modes were studied using heavy ion beam diagnostics on T-10 and ISTTOK and correlation reflectometry on T-10. ISTTOK is equipped with a gallium jet injector and the technical feasibility of gallium jets interacting with plasmas has been investigated in pulsed and ac operation. The first Joint Experiments have clearly demonstrated that small tokamaks are suitable for broad international cooperation to conduct dedicated joint research programmes. Other activities within the IAEA Coordinated Research Project on Joint Research Using Small Tokamaks are also overviewed.

Ideal internal kink modes in a differentially rotating cylindrical plasma

The Velikhov effect leading to magnetorotational instability (MRI) is incorporated into the theory of ideal internal kink modes in a differentially rotating cylindrical plasma column. It is shown that this effect can play a stabilizing role for suitably organized plasma rotation profiles, leading to suppression of MHD (magnetohydrodynamic) instabilities in magnetic confinement systems. The role of this effect in the problem of the Suydam and the m = 1 internal kink modes is elucidated, where m is the poloidal mode number.

Chaotic transport in reversed shear tokamaks

For tokamak models using simplified geometries and reversed shear plasma profiles, we have numerically investigated how the onset of Lagrangian chaos at the plasma edge may affect the plasma confinement in two distinct but closely related problems. Firstly, we have considered the motion of particles in drift waves in the presence of an equilibrium radial electric field with shear. We have shown that the radial particle transport caused by this motion is selective in phase space, being determined by the resonant drift waves and depending on the parameters of both the resonant waves and the electric field profile. Moreover, we have shown that an additional transport barrier may be created at the plasma edge by increasing the electric field. In the second place, we have studied escape patterns and magnetic footprints of chaotic magnetic field lines in the region near a tokamak wall, when there are resonant modes due to the action of an ergodic magnetic limiter. A non-monotonic safety factor profile has been used in the analysis of field line topology in a region of negative magnetic shear. We have observed that, if internal modes are perturbed, the distributions of field line connection lengths and magnetic footprints exhibit spatially localized escape channels. For typical physical parameters of a fusion plasma, the two Lagrangian chaotic processes considered in this work can be effective in usual conditions so as to influence plasma confinement. The reversed shear effects discussed in this work may also contribute to evaluate the transport barrier relevance in advanced confinement scenarios in future tokamak experiments.

Decomposition of spectral density in individual eigenvalue contributions

The eigenvalue densities of two random matrix ensembles, the Wigner Gaussian matrices and the Wishart covariant matrices, are decomposed in the contributions of each individual eigenvalue distribution. It is shown that the fluctuations of all eigenvalues, for medium matrix sizes, are described with a good precision by nearly normal distributions.

Trojan Horse Method: a useful tool for electron screening effect investigation

Direct measurements in the last decades have highlighted a new problem related to the lowering of the Coulomb barrier between the interacting nuclei due to the presence of the ""electron screening"" in the laboratory measurements. It was systematically observed that the presence of the electronic cloud around the interacting ions in measurements of nuclear reactions cross sections at astrophysical energies gives rise to an enhancement of the astrophysical S(E)-factor as lower and lower energies are explored [1]. Moreover, at present Such an effect is not well understood as the value of the potential for screening extracted from these measurements is higher than the tipper limit of theoretical predictions (adiabatic limit). On the other hand, the electron screening potential in laboratory measurement is different from that occurring in stellar plasmas thus the quantity of interest in astrophysics is the so-called ""bare nucleus cross section"". This quantity can only be extrapolated in direct measurements. These are the reasons that led to a considerable growth on interest in indirect measurement techniques and in particular the Trojan Horse Method (THM) [2,3]. Results concerning the bare nucleus cross sections measurements will be shown in several cases of astrophysical interest. In those cases the screening potential evaluated by means of the THM will be compared with the adiabatic limit and results arising from extrapolation in direct measurements.

OSCILLATOR STRENGTHS AND LIFETIMES FOR THE P XII

Semi-empirical weighted oscillator strengths (gf) and lifetimes presented in this work for all experimentally known electric dipole P XII spectral lines and energy levels were computed within a multiconfiguration Hartree-Fock relativistic approach. In this calculation, the electrostatic parameters were optimized by a least-squares procedure in order to improve the adjustment to experimental energy levels. The method produces lifetime and gf values that are in agreement with intensity observations used for the interpretation of spectrograms of solar and laboratory plasmas.