Generalised collisional-radiative model for light elements. A: Data for the Li isonuclear sequence

A first stage collision database is assembled which contains electron-impact effective collision strengths, and ionization and recombination rate coefficients for Li, Li+, and Li2+. The first stage database is constructed using the R-matrix with pseudo-states, time-dependent close-coupling, converged close-coupling, and perturbative distorted-wave methods. A second stage collision database is then assembled which contains generalized collisional-radiative and radiated power loss coefficients. The second stage database is constructed by solution of collisional-radiative equations in the quasi-static equilibrium approximation using the first stage database. Both collision database stages reside in electronic form at the ORNL Controlled Fusion Atomic Data Center and in the ADAS database, and are easily accessed over the worldwide internet. ?? 2006 Elsevier Inc. All rights reserved.

DK-1D: a drift-kinetic simulation tool for modelling the shear Alfvén wave and its interaction with collisionless plasma

We present a highly accurate tool for the simulation of shear Alfven waves (SAW) in collisionless plasma. SAW are important in space plasma environments because for small perpendicular scale lengths they can support an electric field parallel to the ambient magnetic field. Electrons can be accelerated by the parallel electric field and these waves have been implicated as the source of vibrant auroral displays. However, the parallel electric field carried by SAW is small in comparison to the perpendicular electric field of the wave, making it difficult to measure directly in the laboratory, or by satellites in the near-Earth plasma environment. In this paper, we present a simulation code that provides a means to study in detail the SAW-particle interaction in both space and laboratory plasma. Using idealised, small-amplitude propagating waves with a single perpendicular wavenumber, the simulation code accurately reproduces the damping rates and parallel electric field amplitudes predicted by linear theory for varying temperatures and perpendicular scale lengths. We present a rigorous kinetic derivation of the parallel electric field strength for small-amplitude SAW and show that commonly-used inertial and kinetic approximations are valid except for where the ratio of thermal to Alfv\'{e}n speed is between 0.7 and 1.0. We also present nonlinear simulations of large-amplitude waves and show that in cases of strong damping, the damping rates and parallel electric field strength deviate from linear predictions when wave energies are greater than only a few percent of the plasma kinetic energy, a situation which is often observed in the magnetosphere. The drift-kinetic code provides reliable, testable predictions of the parallel electric field strength which can be investigated directly in the laboratory, and will help to bridge the gap between studies of SAW in man-made and naturally occuring plasma.

Identification of geodesic chirping Alfven modes and q-factor estimation in hot core tokamak plasmas in ASDEX Upgrade

Alfven eigenmodes (AE) driven by ion cyclotron resonance heating are usually registered by different diagnostic channels in the hot core plasmas of large tokamaks like JET and ASDEX Upgrade. These AE appear very near to the extremum points of Alfven wave continuum, which is modified by the geodesic effect due to poloidal mode coupling. It is shown that the AE spectrum may be explored as the magnetic spectroscopy (like Alfven cascades by Sharapov et al 2001 Phys. Lett. A 289 127) to determine the q-factor minimum and geodesic frequency at the magnetic axis in standard sawtoothed discharges without reversed shear.

Rotation effect on geodesic and zonal flow modes in tokamak plasmas with isothermal magnetic surfaces

The electrostatic geodesic mode oscillations are investigated in rotating large aspect ratio tokamak plasmas with circular isothermal magnetic surfaces. The analysis is carried out within the magnetohydrodynamic model including heat flux to compensate for the non-adiabatic pressure distribution along the magnetic surfaces in plasmas with poloidal rotation. Instead of two standard geodesic modes, three geodesic continua are found. The two higher branches of the geodesic modes have a small frequency up-shift from ordinary geodesic acoustic and sonic modes due to rotation. The lower geodesic continuum is a newzonal flowmode (geodesic Doppler mode) in plasmas with mainly poloidal rotation. Limits to standard geodesic modes are found. Bifurcation of Alfven continuum by geodesic modes at the rational surfaces is also discussed. Due to that, the frequency of combined geodesic continuum extends from the poloidal rotation frequency to the ion-sound band that can have an important role in suppressing plasma turbulence.

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.

Shearless transport barriers in magnetically confined plasmas

Shearless transport barriers appear in confined plasmas due to non-monotonic radial profiles and cause localized reduction of transport even after they have been broken. In this paper we summarize our recent theoretical and experimental research on shearless transport barriers in plasmas confined in toroidal devices. In particular, we discuss shearless barriers in Lagrangian magnetic field line transport caused by non-monotonic safety factor profiles. We also discuss evidence of particle transport barriers found in the TCABR Tokamak (University of Sao Paulo) and the Texas Helimak (University of Texas at Austin) in biased discharges with non-monotonic plasma flows.

Divertor map with freedom of geometry and safety factor profile

An explicit, area-preserving and integrable magnetic field line map for a single-null divertor tokamak is obtained using a trajectory integration method to represent equilibrium magnetic surfaces. The magnetic surfaces obtained from the map are capable of fitting different geometries with freely specified position of the X-point, by varying free model parameters. The safety factor profile of the map is independent of the geometric parameters and can also be chosen arbitrarily. The divertor integrable map is composed of a nonintegrable map that simulates the effect of external symmetry-breaking resonances, so as to generate a chaotic region near the separatrix passing through the X-point. The composed field line map is used to analyze escape patterns (the connection length distribution and magnetic footprints on the divertor plate) for two equilibrium configurations with different magnetic shear profiles at the plasma edge.

Modelling of spectral properties and population kinetics studies of inertial fusión and laboratory astrophysical plasmas

Fundamental research and modelling in plasma atomic physics continue to be essential for providing basic understanding of many different topics relevant to high-energy-density plasmas. The Atomic Physics Group at the Institute of Nuclear Fusion has accumulated experience over the years in developing a collection of computational models and tools for determining the atomic energy structure, ionization balance and radiative properties of, mainly, inertial fusion and laser-produced plasmas in a variety of conditions. In this work, we discuss some of the latest advances and results of our research, with emphasis on inertial fusion and laboratory-astrophysical applications.

Supersonic regime of the Hall-magnetohydrodynamics resistive tearing instability

An earlier analysis of the Hall-magnetohydrodynamics (MHD) tearing instability [E. Ahedo and J. J. Ramos, Plasma Phys. Controlled Fusion 51, 055018 (2009)] is extended to cover the regime where the growth rate becomes comparable or exceeds the sound frequency. Like in the previous subsonic work, a resistive, two-fluid Hall-MHD model with massless electrons and zero-Larmor-radius ions is adopted and a linear stability analysis about a force-free equilibrium in slab geometry is carried out. A salient feature of this supersonic regime is that the mode eigenfunctions become intrinsically complex, but the growth rate remains purely real. Even more interestingly, the dispersion relation remains of the same form as in the subsonic regime for any value of the instability Mach number, provided only that the ion skin depth is sufficiently small for the mode ion inertial layer width to be smaller than the macroscopic lengths, a generous bound that scales like a positive power of the Lundquist number

Inverse bremsstrahlung absorption in spherical laser targets

Inverse bremsstrahlung has been incorporated into an analytical model of the expanding corona of a laser-irradiated spherical target. Absorption decreases slowly with increasing intensity, in agreement with some numerical simulations, and contrary to estimates from simple models in use up to now, which are optimistic at low values of intensity and very pessimistic at high values. Present results agree well with experimental data from many laboratories; substantial absorption is found up to moderate intensities,say below IOl5 W cm-2 for 1.06 pm light. Anomalous absorption, wher, included in the analysis, leaves practically unaffected the ablation pressure and mass ablation rate, for given absorbed intensity. Universal results are given in dimensionless fom.

Resonant absorption in a plasma step profile

Resonant absorption of p-polarized light shined on a plane-layered plasma with a step profile, is discussed as a function of wavelength (or critical density n,) of the light: for simplicity the incidence angle is assumed small. If n, lies within or above the step, the absorption A is given by Ginzburg’s result modified by strong reflections at the foot and top of the step. The absorption above is total for particular values of nc and U. For n, crossing the top of the density step the absorption is not monotonical: it exhibits a minimum that vanishes for zero radius of curvature U there and zero collision frequency 1’ (A - Iln VI-’). The results are applied to the profile produced by irradiating a solid target with a high-intensity pulse that steepens the plasma by radiation pressure.

Plasma kinetics issues in an ESA study for a plasma laboratory in space

A study supported by the European Space Agency (ESA), in the context of its General Studies Programme, performed an investigation of the possible use of space for studies in pure and applied plasma physics, in areas not traditionally covered by ‘space plasma physics’. A set of experiments have been identified that can potentially provide access to new phenomena and to allow advances in several fields of plasma science. These experiments concern phenomena on a spatial scale (101–104 m) intermediate between what is achievable on the ground and the usual solar system plasma observations. Detailed feasibility studies have been performed for three experiments: active magnetic experiments, largescale discharges and long tether–plasma interactions. The perspectives opened by these experiments are discussed for magnetic reconnection, instabilities, MHD turbulence, atomic excited states kinetics, weakly ionized plasmas,plasma diagnostics, artificial auroras and atmospheric studies. The discussion is also supported by results of numerical simulations and estimates.

Electron cooling in probe collection from magnetized plasmas with anomalous transport

The electron-retarding range of the current-voltage characteristic of a flat Langmuir probe perpendicular to a strong magnetic field in a fully ionized plasma is analysed allowing for anomalous (Bohm) cross-field transport and temperature changes in the collection process. With probe size and ion thermal gyroradius comparable, and smaller than the electron mean free path, there is an outer quasineutral region with ion viscosity determinant in allowing nonambipolar parallel and cross flow. A potential overshoot lying either at the base or inside the quasineutral region both makes ions follow Boltzmann's law at negative bias and extends the electron-retarding range to probe bias e(j)p ~ +2Too. Electron heating and cooling occur roughly at positive and negative bias, with a re-minimum around efa ~ - 2 7 ^ ; far from the probe heat conduction cools and heats electrons at and radially away from the probe axis, respectively. The potential overshoot with no thermal effects would reduce the electron current Ie, making the In Ie versus 4>p graph downwards-concave,but cooling further reduces Ie substantially, and may tilt the slope upwards past the temperature minimum. The domain of strict validity of our analysis is narrow in case of low ion mass (deuterium), breaking down with the ion Boltzmann law.

Self-similar deflagration in laser half-plasmas

The self-similar motion of a half-space plasma, generated by a linear pulse of laser radiation absorbed anomalously at the critical density, has been studied. The resulting plasma structure has been completely determined for [pulse duration (critical density)maximum irradiation] large enough