In the quest for new energy sources, the research on controlled thermonuclear fusion 1 has been boosted by the start of the construction phase of the International Thermonuclear Experimental Reactor (ITER). ITER is based on the tokamak magnetic configuration 3, which is the best performing one in terms of energy confinement. Alternative concepts are however actively researched, which in the long term could be considered for a second generation of reactors. Here, we show results concerning one of these configurations, the reversed-field pinch 4,5 (RFP). By increasing the plasma current, a spontaneous transition to a helical equilibrium occurs, with a change of magnetic topology. Partially conserved magnetic flux surfaces emerge within residual magnetic chaos, resulting in the onset of a transport barrier. This is a structural change and sheds new light on the potential of the RFP as the basis for a low-magnetic-field ohmic fusion reactor.The main magnetic field configurations studied for the confinement of toroidal fusion-relevant plasmas are the tokamak 3 , the stellarator 6 and the reversed-field pinch 4,5 (RFP). In the tokamak, a strong magnetic field is produced in the toroidal direction by a set of coils approximating a toroidal solenoid, and the poloidal field generated by a toroidal current flowing into the plasma gives the field lines a weak helical twist. This is the configuration that has been most studied and has achieved the best levels of energy confinement time. Thus, it is the natural choice for the International Thermonuclear Experimental Reactor, which has the mission of demonstrating the scientific and technical feasibility of controlled fusion with magnetic confinement.The RFP, like the tokamak, is axisymmetric and exploits the pinch effect due to a current flowing in a plasma embedded in a toroidal magnetic field. The main difference is that, for a given plasma current, the toroidal magnetic field in a RFP is one order of magnitude smaller than in a tokamak, and is mainly generated by currents flowing in the plasma itself. This feature is underlying the main potential advantage of the RFP as a reactor concept, namely the capability of achieving fusion conditions with ohmic heating only in a much simpler and compact device. In the past, this positive feature was overcome by the poorer stability properties, which led to the growth and saturation of several magnetohydrodynamic (MHD) instabilities, eventually downgrading the confinement performance. These instabilities, represented by Fourier modes in the poloidal and toroidal angles θ and φ as exp [i(mθ − nφ) were considered as an unavoidable ingredient of the dynamo self-organization process 4,8,9 , necessary for the sustainment of the configuration in time. The occurrence of several MHD modes resonating on different plasma layers gives rise to overlapping magnetic islands, which result in a chaotic region, extending over most of the plasma volume 10 , where the magnetic surfaces are destroyed and the confinement level is modest. This conditi...

The behaviour of tungsten in the core of hybrid scenario plasmas in JET with the ITER-like wall is analysed and modelled with a combination of neoclassical and gyrokinetic codes. In these discharges, good confinement conditions can be maintained only for the first 2–3 s of the high power phase. Later W accumulation is regularly observed, often accompanied by the onset of magneto-hydrodynamical activity, in particular neoclassical tearing modes (NTMs), both of which have detrimental effects on the global energy confinement. The dynamics of the accumulation process is examined, taking into consideration the concurrent evolution of the background plasma profiles, and the possible onset of NTMs. Two time slices of a representative discharge, before and during the accumulation process, are analysed with two independent methods, in order to reconstruct the W density distribution over the poloidal cross-section. The same time slices are modelled, computing both neoclassical and turbulent transport components and consistently including the impact of centrifugal effects, which can be significant in these plasmas, and strongly enhance W neoclassical transport. The modelling closely reproduces the observations and identifies inward neoclassical convection due to the density peaking of the bulk plasma in the central region as the main cause of the accumulation. The change in W neoclassical convection is directly produced by the transient behaviour of the main plasma density profile, which is hollow in the central region in the initial part of the high power phase of the discharge, but which develops a significant density peaking very close to the magnetic axis in the later phase. The analysis of a large set of discharges provides clear indications that this effect is generic in this scenario. The unfavourable impact of the onset of NTMs on the W behaviour, observed in several discharges, is suggested to be a consequence of a detrimental combination of the effects of neoclassical transport and of the appearance of an island.

Recent developments in theory-based modelling of core heavy impurity transport are presented, and shown to be necessary for quantitative description of present experiments in JET and ASDEX Upgrade. The treatment of heavy impurities is complicated by their large mass and charge, which result in a strong response to plasma rotation or any small background electrostatic field in the plasma, such as that generated by anisotropic external heating. These forces lead to strong poloidal asymmetries of impurity density, which have recently been added to numerical tools describing both neoclassical and turbulent transport. Modelling predictions of the steady-state two-dimensional tungsten impurity distribution are compared with experimental densities interpreted from soft Xray diagnostics. The modelling identifies neoclassical transport enhanced by poloidal asymmetries as the dominant mechanism responsible for tungsten accumulation in the central core of the plasma. Depending on the bulk plasma profiles, neoclassical temperature screening can prevent accumulation, and can be enhanced by externally heated species, demonstrated here in ICRH plasmas.

The dependence of plasma transport and confinement on the main hydrogenic ion isotope mass is of fundamental importance for understanding turbulent transport and, therefore, for accurate extrapolations of confinement from present tokamak experiments, which typically use a single hydrogen isotope, to burning plasmas such as ITER, which will operate in deuterium-tritium mixtures. Knowledge of the dependence of plasma properties and edge transport barrier formation on main ion species is critical in view of the initial, low-activation phase of ITER operations in hydrogen or helium and of its implications on the subsequent operation in deuterium-tritium. The favourable scaling of global energy confinement time with isotope mass, which has been observed in many tokamak experiments, remains largely unexplained theoretically. Moreover, the mass scaling observed in experiments varies depending on the plasma edge conditions. In preparation for upcoming deuterium-tritium experiments in the JET tokamak with the ITER-like Be/W Wall (JET-ILW), a thorough experimental investigation of isotope effects in hydrogen, deuterium and tritium plasmas is being carried out, in order to provide stringent tests of plasma energy, particle and momentum transport models. Recent hydrogen and deuterium isotope experiments in JET-ILW on L-H power threshold, L-mode and H-mode confinement are reviewed and discussed in the context of past and more recent isotope experiments in tokamak plasmas, highlighting common elements as well as contrasting observations that have been reported. The experimental findings are discussed in the context of fundamental aspects of plasma transport models.

Recent progress in the understanding and prediction of the tungsten behaviour in the core of JET H-mode plasmas with ITER-like wall is presented. Particular emphasis is given to the impact of poloidal asymmetries of the impurity density. In particular, it is shown that the predicted reduction of temperature screening induced by the presence of low field side localization of the tungsten density produced by the centrifugal force is consistent with the observed tungsten behaviour in a JET discharge in H-mode baseline scenario. This provides first evidence of the role of poloidal asymmetries in reducing the strength of temperature screening. The main differences between plasma parameters in JET baseline and hybrid scenario discharges which affect the impact of poloidally asymmetric density on the tungsten radial transport are identified. This allows the conditions by which tungsten accumulation can be avoided to be more precisely defined.

This paper describes the behavior of nickel in low confinement (L-mode) and high confinement (H-mode) Joint European Torus (JET) discharges [P. J. Lomas, Plasma Phys. Control. Fusion 31, 1481 (1989)] characterized by the application of radio-frequency (rf) power heating and featuring ITER (International Thermonuclear Experimental Reactor) relevant collisionality. The impurity transport is analyzed on the basis of perturbative experiments (laser blow off injection) and is compared with electron heat and deuterium transport. In the JET plasmas analyzed here, ion cyclotron resonance heating (ICRH) is applied either in mode conversion (MC) to heat the electrons or in minority heating (MH) to heat the ions. The two heating schemes have systematically different effects on nickel transport, yielding flat or slightly hollow nickel density profiles in the case of ICRH in MC and peaked nickel density profiles in the case of rf applied in MH. Accordingly, both diffusion coefficients and pinch velocities of nickel are found to be systematically different. Linear gyrokinetic calculations by means of the code GS2 [M. Kotschenreuther, G. Rewoldt, and W.M. Tang, Comput. Phys. Commun. 88, 128 (1995)] provide a possible explanation of such different behavior by exploring the effects produced by the different microinstabilities present in these plasmas. In particular, trapped electron modes driven by the stronger electron temperature gradients measured in the MC cases, although subdominant, produce a contribution to the impurity pinch directed outwards that is qualitatively in agreement with the pinch reversal found in the experiment. Particle and heat diffusivities appear to be decoupled in MH shots, with χe and DD≫DNi, and are instead quite similar in the MC ones. In the latter case, nickel transport appears to be driven by the same turbulence that drives the electron heat transport and is sensitive to the value of the electron temperature gradient length. These findings give ground to the idea that in ITER it should be possible to find conditions in which the risk of accumulation of metals such as nickel can be contained.

The scan of Ion Cyclotron Resonant Heating power has been used to systematically study the pump out effect of central electron heating on impurities such as Ni and Mo in H mode low collisionality discharges in JET. The transport parameters of Ni and Mo have been measured by introducing a transient perturbation on their densities via the Laser Blow Off technique. Without ICRH, Ni and Mo density profiles are typically peaked. The application of ICRH, induces on Ni and Mo in the plasma center (at normalized poloidal flux r = 0.2) an outward drift approximately proportional to the amount of injected power. Above a threshold, of about 3MW of ICRH power in the specific case, the radial flow of Ni and Mo changes from inward to outward and the impurity profiles, extrapolated to stationary conditions, become hollow. At mid radius the impurity profiles become flat or only slightly hollow. In the plasma centre the variation of the pinch parameter v/D of Ni is particularly well correlated with the change of the ion temperature gradient, in qualitative agreement with the neoclassical theory. However, the experimental radial velocity is larger than the neoclassical one by up to one order of magnitude. Gyrokinetic simulations of the radial impurity fluxes induced by electrostatic turbulence do not foresee a flow reversal in the analyzed discharges.

The evolution of the JET high performance hybrid scenario, including central accumulation of the tungsten (W) impurity, is reproduced with predictive multi-channel integrated modelling over multiple confinement times using first-principle based core transport models. 8 transport channels (𝑇 𝑖 , 𝑇 𝑒 , 𝑗, 𝑛 𝐷 , 𝑛 𝐵𝑒 , 𝑛 𝑁𝑖 , 𝑛 𝑊 , 𝜔) are modelled predictively, with self-consistent sources, radiation and magnetic equilibrium, yielding a system with multiple non-linearities: This system can reproduce the observed radiative temperature collapse after several confinement times. W is transported inward by neoclassical convection driven by the main ion density gradients and enhanced by poloidal asymmetries due to centrifugal acceleration. The slow evolution of the bulk density profile sets the timescale for W accumulation. Modelling this phenomenon requires a turbulent transport model capable of accurately predicting particle and momentum transport (QuaLiKiz) and a neoclassical transport model including the effects of poloidal asymmetries (NEO) coupled to an integrated plasma simulator (JINTRAC). The modelling capability is applied to optimise the available actuators to prevent W accumulation, and to extrapolate in power and pulse length. Central NBI heating is preferred for high performance, but gives central deposition of particles and torque which increase the risk of W accumulation by increasing density peaking and poloidal asymmetry. The primary mechanism for ICRH to control W in JET is via its impact through turbulence in reducing main ion density peaking (which drives inward neoclassical convection), increased temperature screening and turbulent W diffusion. The anisotropy from ICRH also reduces poloidal asymmetry, but this effect is negligible in high rotation JET discharges. High power ICRH near the axis can sensitively mitigate against W accumulation, and dominant ion heating (e.g. He-3 minority) is predicted to provide more resilience to W accumulation than dominant electron heating (e.g. H minority) in the JET hybrid scenario. Extrapolation to DT plasmas finds 17.5MW of fusion power and improved confinement compared to DD, due to reduced ion-electron energy exchange, and increased Ti/Te stabilisation of ITG instabilities. The turbulence reduction in DT increases density peaking and accelerates the arrival of W on axis; this may be mitigated by reducing the penetration of the beam particle source with an increased pedestal density.