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.

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.

We present an ultrafast neural network (NN) model, QLKNN, which predicts core tokamak transport heat and particle fluxes. QLKNN is a surrogate model based on a database of 300 million flux calculations of the quasilinear gyrokinetic transport model QuaLiKiz. The database covers a wide range of realistic tokamak core parameters. Physical features such as the existence of a critical gradient for the onset of turbulent transport were integrated into the neural network training methodology. We have coupled QLKNN to the tokamak modelling framework JINTRAC and rapid control-oriented tokamak transport solver RAPTOR. The coupled frameworks are demonstrated and validated through application to three JET shots covering a representative spread of H-mode operating space, predicting turbulent transport of energy and particles in the plasma core. JINTRAC-QLKNN and RAPTOR-QLKNN are able to accurately reproduce JINTRAC-QuaLiKiz T i,e and n e profiles, but 3 to 5 orders of magnitude faster. Simulations which take hours are reduced down to only a few tens of seconds. The discrepancy in the final source-driven predicted profiles between QLKNN and QuaLiKiz is on the order 1%-15%. Also the dynamic behaviour was well captured by QLKNN, with differences of only 4%-10% compared to JINTRAC-QuaLiKiz observed at mid-radius, for a study of density buildup following the L-H transition. Deployment of neural network surrogate models in multi-physics integrated tokamak modelling is a promising route towards enabling accurate and fast tokamak scenario optimization, Uncertainty Quantification, and control applications.

Detailed experimental studies of ion heat transport have been carried out in JET exploiting the upgrade of Active Charge Exchange Spectroscopy and the availability of multi-frequency ICRH with 3 He minority. The determination of ion threshold and stiffness offers unique opportunities for validation of the well-established theory of Ion Temperature Gradient driven modes. Ion stiffness is observed to decrease strongly in presence of toroidal rotation when the magnetic shear is sufficiently low. This effect is dominant with respect to the well-known w ExB threshold up-shift and plays a major role in enhancing core confinement in Hybrid regimes and Ion Internal Transport Barriers. The effects of T e /T i and s/q on ion threshold are found rather weak in the domain explored. Quasi-linear fluid/gyro-fluid and linear/non-linear gyro-kinetic simu lations have been carried out. Whilst threshold predictions show good match with experimental observations, some significant discrepancies are found on the stiffness behaviour.