Radial variation of heat transport in L-mode JET discharges

In this paper, we analyze heat transport in the JET tokamak using data from its high resolution ECE diagnostic and analyses based on the transfer entropy (TE). The analysis reveals that heat transport is not smooth and continuous, but is characterized by 'trapping regions' separated by `minor transport barriers'. Meat may 'jump over' these barriers and when the heating power is raised, this 'jumping' behavior becomes more prominent. To check that our results are relevant for global heat transport, we deduced an effective diffusion coefficient from the TE results. Both its value and overall radial variation are consistent with heat diffusivities reported in literature. The detailed radial structure of the effective diffusion coefficient was shown to be linked to the mentioned minor transport barriers

14 MeV calibration of JET neutron detectors-phase 1: Calibration and characterization of the neutron source

In view of the planned DT operations at JET, a calibration of the JET neutron monitors at 14 MeV neutron energy is needed using a 14 MeV neutron generator deployed inside the vacuum vessel by the JET remote handling system. The target accuracy of this calibration is 10% as also required by ITER, where a precise neutron yield measurement is important, e.g. for tritium accountancy. To achieve this accuracy, the 14 MeV neutron generator selected as the calibration source has been fully characterised and calibrated prior to the in-vessel calibration of the JET monitors. This paper describes the measurements performed using different types of neutron detectors, spectrometers, calibrated long counters and activation foils which allowed us to obtain the neutron emission rate and the anisotropy of the neutron generator, i.e.The neutron flux and energy spectrum dependence on emission angle, and to derive the absolute emission rate in 4π sr. The use of high resolution diamond spectrometers made it possible to resolve the complex features of the neutron energy spectra resulting from the mixed D/T beam ions reacting with the D/T nuclei present in the neutron generator target. As the neutron generator is not a stable neutron source, several monitoring detectors were attached to it by means of an ad hoc mechanical structure to continuously monitor the neutron emission rate during the in-vessel calibration. These monitoring detectors, two diamond diodes and activation foils, have been calibrated in terms of neutrons/counts within ± 5% total uncertainty. A neutron source routine has been developed, able to produce the neutron spectra resulting from all possible reactions occurring with the D/T ions in the beam impinging on the Ti D/T target. The neutron energy spectra calculated by combining the source routine with a MCNP model of the neutron generator have been validated by the measurements. These numerical tools will be key in analysing the results from the in-vessel calibration and to derive the response of the JET neutron detectors to DT plasma neutrons starting from the response to the generator neutrons, and taking into account all the calibration circumstances

Investigation of deuterium trapping and release in the JET divertor during the third ILW campaign using TDS

Selected set of samples from JET ITER-Like Wall (JET-ILW) divertor tiles exposed in 2015-2016 has been analysed using Thermal Desorption Spectrometry (TDS). The deuterium (D) amounts obtained with TDS were compared with Nuclear Reaction Analysis (NRA). The highest amount of D was found on the top part of inner divertor which has regions with the thickest deposited layers as for divertor tiles removed in 2014. This area resides deep in the scrape-off layer and plasma configurations for the second (ILW-2, 2013-2014) and the third (ILW-3, 2015-2016) JET-ILW campaigns were similar. Agreement between TDS and NRA is good on the apron of Tile 1 and on the upper vertical region whereas on the lower vertical region of Tile 1 the NRA results are clearly smaller than the TDS results. Inner divertor Tile 3 has somewhat less D than Tiles 0 and 1, and the D amount decreases towards the lower part of the tile. The D retention at the divertor inner and outer corner regions is not symmetric as there is more D retention poloidally at the inner than at the outer divertor corner. In most cases the TDS spectra for the ILW-3 samples are different from the corresponding ILW-2 spectra because HD and D-2 release occurs at higher temperatures than from the ILW-2 samples indicating that the low energy traps have been emptied during the plasma operations and that D is either in the energetically deep traps or located deeper in the sample

Direct gyrokinetic comparison of pedestal transport in JET with carbon and ITER-like walls

This paper compares the gyrokinetic instabilities and transport in two representative JET pedestals, one (pulse 78697) from the JET configuration with a carbon wall (C) and another (pulse 92432) from after the installation of JET's ITER-like Wall (ILW). The discharges were selected for a comparison of JET-ILW and JET-C discharges with good confinement at high current (3 MA, corresponding also to low rho(*)) and retain the distinguishing features of JET-C and JET-ILW, notably, decreased pedestal top temperature for JET-ILW. A comparison of the profiles and heating power reveals a stark qualitative difference between the discharges: the JET-ILW pulse (92432) requires twice the heating power, at a gas rate of 1.9 x 10(22) e s(-1), to sustain roughly half the temperature gradient of the JET-C pulse (78697), operated at zero gas rate. This points to heat transport as a central component of the dynamics limiting the JET-ILW pedestal and reinforces the following emerging JET-ILW pedestal transport paradigm, which is proposed for further examination by both theory and experiment. ILW conditions modify the density pedestal in ways that decrease the normalized pedestal density gradient a/L-n, often via an outward shift in relation to the temperature pedestal. This is attributable to some combination of direct metal wall effects and the need for increased fueling to mitigate tungsten contamination. The modification to the density profile increases eta = L-n/L-T, thereby producing more robust ion temperature gradient (ITG) and electron temperature gradient driven instability. The decreased pedestal gradients for JET-ILW (92432) also result in a strongly reduced E x B shear rate, further enhancing the ion scale turbulence. Collectively, these effects limit the pedestal temperature and demand more heating power to achieve good pedestal performance. Our simulations, consistent with basic theoretical arguments, find higher ITG turbulence, stronger stiffness, and higher pedestal transport in the ILW plasma at lower rho(*)

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

Approximate analytic expressions using Stokes model for tokamak polarimetry and their range of validity

The analysis of the polarimetry measurements has the aim of validating models (De Marco and Segre 1972 Plasma Phys. 14 245), with a careful attention to the clarification of their limits of application. In this paper a new approximation method is introduced, the so-called special constant Omega direction (SCOD), which gives an analytical solution to the polarimetry exact Stokes model equations. The available approximate solutions (including SCOD) of the polarimetry propagation equations are presented, compared and their application limits determined, using a reference tokamak configuration, which is a simplified equilibrium for a circular tokamak. The SCOD approximation is compared successfully to the Stokes model in the context also of equilibria evaluated for two JET discharges. The approximation methods are analytical or very simple mathematical expressions which can also be used in equilibrium codes for their optimization

Modelling of the effect of ELMs on fuel retention at the bulk W divertor of JET

Effect of ELMs on fuel retention at the bulk W target of JET ITER-Like Wall was studied with multi-scale calculations. Plasma input parameters were taken from ELMy H-mode plasma experiment. The energetic intra-ELM fuel particles get implanted and create near-surface defects up to depths of few tens of nm, which act as the main fuel trapping sites during ELMs. Clustering of implantation-induced vacancies were found to take place. The incoming flux of inter-ELM plasma particles increases the different filling levels of trapped fuel in defects. The temperature increase of the W target during the pulse increases the fuel detrapping rate. The inter-ELM fuel particle flux refills the partially emptied trapping sites and fills new sites. This leads to a competing effect on the retention and release rates of the implanted particles. At high temperatures the main retention appeared in larger vacancy clusters due to increased clustering rate

Energetic ion losses 'channeling' mechanism and strategy for mitigation

Results from two different sets of JET experiments are presented. In experiments in which toroidicity-induced Alfven eigenmodes (TAEs) localized at different radial locations had the same frequencies and toroidal mode numbers, the occurrence of enhanced losses after the excitation of TAEs in the core of the plasma was observed. On the contrary, enhanced losses were not observed if the TAEs localized at different radial locations had different frequencies and toroidal mode numbers. Numerical modeling indicates that, in the first set of experiments, the enhanced losses were caused by a combined effect of the TAEs localized at different radial locations. The TAEs localized in the plasma core transported energetic ions from the core to outer regions of the plasma. Then, the TAEs localized in outer regions of the plasma interacted with these ions just transported by the core-localized TAEs causing a further radial displacement of the ions to the plasma edge. This process eventually ends up causing the loss of the resonant ions. In the second set of experiments, it was found that TAEs localized in the plasma core and in outer regions did not interact with the same ions and so no enhanced losses were measured. Sheared profiles of the safety factor combined with flat mass density profiles lead to larger differences on the frequencies of the TAEs localized at different radial locations, eventually avoiding loss of energetic ions through the described mechanism

Radial variation of heat transport in L-mode JET discharges

In this paper, we analyze heat transport in the JET tokamak using data from its high resolution ECE diagnostic and analyses based on the transfer entropy (TE). The analysis reveals that heat transport is not smooth and continuous, but is characterized by 'trapping regions' separated by `minor transport barriers'. Meat may 'jump over' these barriers and when the heating power is raised, this 'jumping' behavior becomes more prominent. To check that our results are relevant for global heat transport, we deduced an effective diffusion coefficient from the TE results. Both its value and overall radial variation are consistent with heat diffusivities reported in literature. The detailed radial structure of the effective diffusion coefficient was shown to be linked to the mentioned minor transport barriers