Modelling third harmonic ion cyclotron acceleration of deuterium beams for JET fusion product studies experiments

Recent JET experiments have been dedicated to the studies of fusion reactions between deuterium (D) and Helium-3 (3He) ions using neutral beam injection (NBI) in synergy with third harmonic ion cyclotron radio-frequency heating (ICRH) of the beam. This scenario generates a fast ion deuterium tail enhancing DD and D3He fusion reactions. Modelling and measuring the fast deuterium tail accurately is essential for quantifying the fusion products. This paper presents the modelling of the D distribution function resulting from the NBI+ICRF heating scheme, reinforced by a comparison with dedicated JET fast ion diagnostics, showing an overall good agreement. Finally, a sawtooth activity for these experiments has been observed and interpreted using SPOT/RFOF simulations in the framework of Porcelli's theoretical model, where NBI+ICRH accelerated ions are found to have a strong stabilizing effect, leading to monster sawteeth

Fusion product studies via fast ion D-D and D-3He fusion on JET

Dedicated fast ion D-D and D-3He fusion experiments were performed on JET with carbon wall (2008) and ITER-like wall (2014) for testing the upgraded neutron and energetic ion diagnostics of fusion products. Energy spectrum of D-D neutrons was the focus of the studies in pure deuterium plasmas. A significant broadening of the energy spectrum of neutrons born in D-D fast fusion was observed, and dependence of the maximum D and D-D neutron energies on plasma density was established. Diagnostics of charged products of aneutronic D-3He fusion reactions, 3.7 MeV alpha-particles similar to those in D-T fusion, and 14.6 MeV protons, were the focus of the studies in D-3He plasmas. Measurements of 16.4 MeV gamma-rays born in the weak secondary branch of D(3He, γ)5Li reaction were used for assessing D-3He fusion power. For achieving high yield of D-D and D-3He reactions at relatively low levels of input heating power, an acceleration of D beam up to the MeV energy range was used employing 3rd harmonic () ICRH technique. These results were compared to the techniques of D beam injection into D-3He mixture, and 3He-minority ICRH in D plasmas

Overview of the JET results

Since the installation of an ITER-like wall, the JET programme has focused on the consolidation of ITER design choices and the preparation for ITER operation, with a specific emphasis given to the bulk tungsten melt experiment, which has been crucial for the final decision on the material choice for the day-one tungsten divertor in ITER. Integrated scenarios have been progressed with the re-establishment of long-pulse, high-confinement H-modes by optimizing the magnetic configuration and the use of ICRH to avoid tungsten impurity accumulation. Stationary discharges with detached divertor conditions and small edge localized modes have been demonstrated by nitrogen seeding. The differences in confinement and pedestal behaviour before and after the ITER-like wall installation have been better characterized towards the development of high fusion yield scenarios in DT. Post-mortem analyses of the plasma-facing components have confirmed the previously reported low fuel retention obtained by gas balance and shown that the pattern of deposition within the divertor has changed significantly with respect to the JET carbon wall campaigns due to the absence of thermally activated chemical erosion of beryllium in contrast to carbon. Transport to remote areas is almost absent and two orders of magnitude less material is found in the divertor

Micro ion beam analysis for the erosion of beryllium marker tiles in a tokamak limiter

Beryllium limiter marker tiles were exposed to plasma in the Joint European Torus to diagnose the erosion of main chamber wall materials. A limiter marker tile consists of a beryllium coating layer (7-9 mu m) on the top of bulk beryllium, with a nickel interlayer (2-3 mu m) between them. The thickness variation of the beryllium coating layer, after exposure to plasma, could indicate the erosion measured by ion beam analysis with backscattering spectrometry. However, interpretations from broad beam backscattering spectra were limited by the non-uniform surface structures. Therefore, micro-ion beam analysis (mu-IBA) with 3 MeV proton beam for Elastic back scattering spectrometry (EBS) and PIXE was used to scan samples. The spot size was in the range of 3-10 mu m. Scanned areas were analysed with scanning electron microscopy (SEM) as well. Combining results from mu-IBA and SEM, we obtained local spectra from carefully chosen areas on which the surface structures were relatively uniform. Local spectra suggested that the scanned area (approximate to 600 mu m x 1200 mu m) contained regions with serious erosion with only 2-3 mu m coating beryllium left, regions with intact marker tile, and droplets with 90% beryllium. The nonuniform erosion, droplets mainly formed by beryllium, and the possible mixture of beryllium and nickel were the major reasons that confused interpretation from broad beam EBS

Long-lived coupled peeling ballooning modes preceding ELMs on JET

In some JET discharges, type-I edge localised modes (ELMs) are preceded by a class of low-frequency oscillations (Perez et al 2004 Nucl. Fusion 44 609). While in many cases the ELM is triggered during the growth phase of this oscillation, it is also observed that this type of oscillation can saturate and last for several tens of ms until an ELM occurs. In order to identify the nature of these modes, a wide pre-ELM oscillation database, including detailed pedestal profile information, has been assembled and analysed in terms of MHD stability parameters. The existence domain of these pre-ELM oscillations and the statistical distribution of toroidal mode numbers (n) up to n = 16 have been mapped in ballooning alpha (alpha(ball)) and either edge current density (J(edge)) or pedestal collisionality (nu(ee,ped)*) coordinates and compared to linear MHD stability predictions. The pre-ELM oscillations are reliably observed when the J/alpha ratio is high enough for the pedestal to access the coupled peeling-ballooning (PB) domain (aka stability nose). Conversely, when the pedestal is found to be in or near the high-n ballooning domain (which is at low J/alpha ratio), ELMs are usually triggered promptly, i.e. with no detectable pre-ELM oscillations, or with pre-ELM oscillations only observable on ECE whose n appears to be too high to be resolved by the magnetics. Individual discharges can sometimes exhibit a fairly wide range of pre-ELM mode numbers, but for a wider database, the statistical n-number domains are found to be well ordered along the J - alpha stability boundary and behave as expected from PB theory: the higher the J/alpha ratio, the lower the mode's measured n tends to be. Within the measurement uncertainties, the measured n is usually found to be compatible with the most unstable n predicted by the linear stability code MISHKA1. These results confirm the earlier hypothesis that these modes are coupled peeling-ballooning modes, and extend and generalise to higher-mode numbers the work by Huysmans et al (1998 Nucl. Fusion 38 179), who identified the lowest n modes as pure external kink modes. Since the destabilisation of PB modes is widely accepted to give rise to ELMs, the mode saturation and delayed ELM triggering that is sometimes observed is rather unexpected. Possibilities to reconcile the extended lifetime of these modes with current ELM models are briefly discussed, but will require further investigation

Diagnostic of fast-ion energy spectra and densities in magnetized plasmas

The measurement of the energy spectra and densities of alpha-particles and other fast ions are part of the ITER measurement requirements, highlighting the importance of energy-resolved energetic-particle measurements for the mission of ITER. However, it has been found in recent years that the velocity-space interrogation regions of the foreseen energetic-particle diagnostics do not allow these measurements directly. We will demonstrate this for gamma-ray spectroscopy (GRS), collective Thomson scattering (CTS), neutron emission spectroscopy and fast-ion D-alpha spectroscopy by invoking energy and momentum conservation in each case, highlighting analogies and differences between the different diagnostic velocity-space sensitivities. Nevertheless, energy spectra and densities can be inferred by velocity-space tomography which we demonstrate using measurements at JET and ASDEX Upgrade. The measured energy spectra agree well with corresponding simulations. At ITER, alpha-particle energy spectra and densities can be inferred for energies larger than 1.7 MeV by velocity-space tomography based on GRS and CTS. Further, assuming isotropy of the alpha-particles in velocity space, their energy spectra and densities can be inferred by 1D inversion of spectral single-detector measurements down to about 300 keV by CTS. The alpha-particle density can also be found by fitting a model to the CTS measurements assuming the alpha-particle distribution to be an isotropic slowing-down distribution

EDGE2D-EIRENE simulations of the influence of isotope effects and anomalous transport coefficients on near scrape-off layer radial electric field

EDGE2D-EIRENE (the 'code') simulations show that radial electric field, Er, in the near scrape-off layer (SOL) of tokamaks can have large variations leading to a strong local E x B shear greatly exceeding that in the core region. This was pointed out in simulations of JET plasmas with varying divertor geometry, where the magnetic configuration with larger predicted near SOL E-r was found to have lower H-mode power threshold, suggesting that turbulence suppression in the SOL by local E. x. B shear can be a player in the L-H transition physics (Delabie et al 2015 42nd EPS Conf. on Plasma Physics (Lisbon, Portugal, 22-26 June 2015) paper O3.113 (http://ocs.ciemat.es/EPS2015PAP/pdf/O3.113.pdf), Chankin et al 2017 Nucl. Mater. Energy 12 273). Further code modeling of JET plasmas by changing hydrogen isotopes (H-D-T) showed that the magnitude of the near SOL E-r is lower in H cases in which the H-mode threshold power is higher (Chankin et al 2017 Plasma Phys. Control. Fusion 59 045012). From the experiment it is also known that hydrogen plasmas have poorer particle and energy confinement than deuterium plasmas, consistent with the code simulation results showing larger particle diffusion coefficients at the plasma edge, including SOL, in hydrogen plasmas (Maggi et al 2018 Plasma Phys. Control. Fusion 60 014045). All these experimental observations and code results support the hypothesis that the near SOL E x B shear can have an impact on the plasma confinement. The present work analyzes neutral ionization patterns of JET plasmas with different hydrogen isotopes in L-mode cases with fixed input power and gas puffing rate, and its impact on target electron temperature, T-e, and SOL E-r. The possibility of a self-feeding mechanism for the increase in the SOL E-r via the interplay between poloidal E x B drift and target T-e is discussed. It is also shown that reducing anomalous turbulent transport coefficients, particle diffusion and electron and ion heat conductivities, leads to higher peak target T-e and larger E-r, suggesting the possibility of a positive feedback loop, under an implicitly made assumption that the E x B shear in the SOL is capable of suppressing turbulence

A power-balance model of the density limit in fusion plasmas: application to the L-mode tokamak

A power-balance model, with radiation losses from impurities and neutrals, gives a unified description of the density limit (DL) of the stellarator, the L-mode tokamak, and the reversed field pinch (RFP). The model predicts a Sudo-like scaling for the stellarator, a Greenwald- like scaling, alpha I-p(8/9), for the RFP and the ohmic tokamak, a mixed scaling, alpha (PIp4/9)-I-4/9, for the additionally heated L-mode tokamak. In a previous paper (Zanca et al 2017 Nucl. Fusion 57 056010) the model was compared with ohmic tokamak, RFP and stellarator experiments. Here, we address the issue of the DL dependence on heating power in the L-mode tokamak. Experimental data from high-density disrupted L-mode discharges performed at JET, as well as in other machines, arc taken as a term of comparison. The model fits the observed maximum densities better than the pure Greenwald limit