Search

The first divertor was installed in the JET machine between 1992 and 1994 and was operated with carbon tiles and then beryllium tiles in 1994-5. Post-mortem studies after these first experiments demonstrated that most of the impurities deposited in the divertor originate in the main chamber, and that asymmetric deposition patterns generally favouring the inner divertor region result from drift in the scrape-off layer. A new monolithic divertor structure was installed in 1996 which produced heavy deposition at shadowed areas in the inner divertor corner, which is where the majority of the tritium was trapped by co-deposition during the deuterium-tritium experiment in 1997. Different divertor geometries have been tested since such as the Gas-Box and High-Delta divertors; a principle objective has been to predict plasma behaviour, transport and tritium retention in ITER. Transport modelling experiments were carried out at the end of four campaigns by puffing C-13-labelled methane, and a range of diagnostics such as quartz-microbalance and rotating collectors have been installed to add time resolution to the post-mortem analyses. The study of material migration after D-D and D-T campaigns clearly revealed important consequences of fuel retention in the presence of carbon walls. They gave a strong impulse to make a fundamental change of wall materials. In 2010 the carbon divertor and wall tiles were removed and replaced with tiles with Be or W surfaces for the ITER-Like Wall Project. For complete list of authors see http://dx.doi.org/10.1016/j.fusengdes.2018.10.002

In this paper we validate the finite element code RAPLICASOL, which models radiofrequency wave propagation in edge plasmas near ICRF antennas, against calculations with the TOPICA code. We compare the output of both codes for the ASDEX Upgrade 2-strap antenna, and for a 4-strap WEST-like antenna. Although RAPLICASOL requires considerably fewer computational resources than TOPICA, we find that the predicted quantities of experimental interest (including reflection coefficients, coupling resistances, S- and Z-matrix entries, optimal matching settings, and even radiofrequency electric fields) are in good agreement provided we are careful to use the same geometry in both codes.

The MHD response and the penetration of a deuterium shattered pellet into a JET plasma is investigated via the non-linear reduced MHD code JOREK with the neutral gas shielding (NGS) ablation model. The dominant MHD destabilizing mechanism by the injection is identified as the local helical cooling at each rational surface, as opposed to the global current profile contraction. Thus the injected fragments destabilize each rational surface as they pass through them. The injection penetration is found to be much better compared to MGI, with the convective transport caused by core MHD instabilities (e.g. 1/1 kink) contributing significantly to the core penetration. Moreover, the injection with realistic JET SPI system configurations is simulated in order to provide some insights into future operations, and the impact on the total assimilation and penetration depth of varying injection parameters such as the injection velocity or fineness of shattering is assessed. Further, the effect of changing the target equilibrium temperature or q profile on the assimilation and penetration is also investigated. Such analysis will form the basis of further investigation into a desirable configuration for the future SPI system in ITER. EURATOM 633053

Tungsten (W) melting due to transient power loads, for example those delivered by edge localised modes (ELMs), is a major concern for next step fusion devices. A series of experiments has been performed on JET to investigate the dynamics of Type-I ELM-induced transient melting. Following initial exposures in 2013 of a W-lamella with sharp leading edge in the bulk W outer divertor, new experiments have been performed in 2016–2017 on a protruding W-lamella with a 15° slope, allowing direct and spatially resolved (0.85 mm/pixel) observation of the top surface using the IR thermography system viewing from the top of the poloidal cross-section. Thermal and IR analysis have already been conducted assuming the geometrical projection of the parallel heat flux on the W-lamellas, thus ignoring the gyro-radius orbit of plasma particles. Although it is well justified during L-mode or inter-ELM period, the hypothesis becomes questionable during ELM when the ion Larmor radius is larger. The goal of this paper is to extend the previous analysis based on the forward approach to the H-mode discharges and investigate in particular the gyro-radius effect during the Type-I ELMs, those used to achieve transient melting on the slope of the protruding W-lamella. Surface temperatures measured by the IR camera are compared with reconstructed synthetic data from 3D thermal modelling using heat loads derived from optical projection of the parallel heat flux (ignoring the gyro-radius orbit), 2D gyro-radius orbit and particle-in-cell (PIC) simulations describing the influence of finite Larmor-radius effects and electrical potential on the deposited power flux. Results show that the ELM power deposition behaves differently than the optical projection of the parallel heat flux, contrary to the L-mode observations, and may thus be due to the much larger gyro-orbits of the energetic ELM ions in comparison to L-mode or inter-ELM conditions. Keywords: ELM heat flux analysis, Optical projection, PIC modelling, Ion orbit modelling, IR data, Spatial resolution effect

Micro/nanoscopic observations on the surface of the divertor tiles used in the first campaign (2011-2012) of the JET tokamak with ITER-like Wall (JET ILW) have been carried out by means of several material analysis techniques. Previous results from the inner divertor were reported for a single poloidal section of the tile numbers 1, 3 and 4, i.e., upper, vertical and horizontal targets, respectively. The formation of the thick stratified mixed-material deposition layer on tiles 1 and 4, and erosion on tile 3 were identified. This study is mostly focused on the outer divertor: tiles 6, 7 and 8. In contrast to the inner tile, remarkable surface modifications have not been observed on the vertical target (tiles 7 and 8) where sputtering erosion and impurity deposition would have been almost balanced. Only a specific part of tile 6 (horizontal target) located near the exhaust channel was covered with a stratified ("geological-like") mixed-material deposition layer which mainly included Be and Ni with the thickness of similar to 2 mu m. Special feature of this mixed layer was that a certain amount of nitrogen (N) was clearly detected in the layer. Since the concentration of N varied with the depth position, it could be depended on the amount of that gas puffed for plasma edge cooling during the JET experimental campaign. In addition to the outer divertor tiles, a very interesting feature of the local erosion and deposition effects is reported in this paper.

International audience; In Ion-Temperature-Gradient (ITG) driven turbulence, the resonance condition leads to ion particle turbulent transport coefficients significantly larger than electron particle turbulent transport coefficients. This is shown in non-linear gyrokinetic simulations and explained by an analytical quasilinear model. It is then illustrated by JETTO-QuaLiKiz integrated modelling. Large ion particle transport coefficients implies that the ion density profiles are uncorrelated to the corresponding ion source, allowing peaked isotope density profiles even in the absence of core source. This also implies no strong core accumulation of He ash. Furthermore, the relaxation time of the individual ion profiles in a multi-species plasma can be significantly faster than the total density profile relaxation time which is constrained by the electrons. This leads to fast isotope mixing and fast impurity transport in ITG regimes. In Trapped-Electron-Mode (TEM) turbulence , in presence of electron heating about twice the ion heating, the situation is the inverse: ion particle turbulent transport coefficients are smaller than their electron counterpart.

The EUROfusion Work Package JET3 programme, established to enable the technological exploitation of the JET experiments over the next years, includes, within the NEXP subproject, a novel Shutdown Dose Rate (SDR) experiment. Considering its ITER-relevance, SDR experiments at JET represent a unique opportunity to validate the numerical tools for ITER nuclear analysis, through the comparison between numerical predictions and measured quantities (C/E). Within this framework, two active gamma dosimeters based on spherical air-vented ionization chambers (ICs) have been installed in ex-vessel positions close to the horizontal ports of the tokamak in Octants 1 and 2. The first JET campaign exploited in the novel SDR experiment is the latest 5-week Deuterium-Deuterium campaign (c36b), which achieved the best results in recent years in terms of high power operation. The present work is dedicated to the analysis of dose rate measurements carried out during this campaign and after shutdown. Proper correction factors are evaluated and applied to the instrument reading, while influence quantities and error sources are analyzed in order to calculate the overall experimental uncertainty.

Plasma tomography consists of reconstructing a two-dimensional radiation profile of a poloidal cross section of a fusion device based on line-integrated measurements along several lines of sight. The reconstruction process is computationally intensive, and in practice, only a few reconstructions are usually computed per pulse. In this work, we trained a deep neural network based on a large collection of sample tomograms that have been produced at JET over several years. Once trained, the network is able to reproduce those results with high accuracy. More importantly, it can compute all the tomographic reconstructions for a given pulse in just a few seconds. This makes it possible to visualize several phenomena—such as plasma heating, disruptions, and impurity transport—over the course of the entire pulse.

Beryllium oxide (BeO) and deuteroxide (BeOxDy) have been found on the melted zone of a beryllium tile extracted from the upper dump plate of JET-ILW (2011-2012 campaign). Results have been obtained using Raman microscopy, which is sensitive to both the chemical bond and crystal structure, with a micrometric lateral resolution. BeO is found with a wurtzite crystal structure. BeOxDy is found as three different types which are not the beta-phase but behaves as molecular species like Be(OD)(2), O(Be-D)(2) and DBeOD. The presence of a small amount of trapped D2O is also suspected. Our results therefore strongly suggest that D trapping occurs after melting through the formation of deuteroxides. The temperature increase favors the formation of crystal BeO which favors deuterium trapping through OD bonding. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2018.11.008

Measurements of neutron streaming through penetrations in biological shields are being carried out at JET fusion device by means of thermoluminescence (TL) detectors with the objective to validate the neutronics codes and nuclear data widely applied in ITER nuclear analyses in a real fusion environment. TLDs response due to the neutron component of the radiation field is related to the neutron fluence in a well-defined neutron energy spectrum. Therefore, a TLDs calibration in real fusion radiation fields is necessary to allow neutron fluence from TL measurements at JET to be more precisely calculated. Hence, a MCP-N and MCP-7 TLDs produced at the IFJ PAN in Kraków were calibrated at the ENEA facilities of Frascati and Casaccia laboratories. The obtained results have been analysed and new calibration factors are proposed. The detection system based on TLDs developed and calibrated for JET experiments can then be generally applied not only to fusion neutron fields, but also to ITER to monitor the neutron fluence outside the biological shield. © 2018

Understanding the details of the correlation between time series is an essential step on the route to assessing the causal relation between systems. Traditional statistical indicators, such as the Pearson correlation coefficient and the mutual information, have some significant limitations. More recently, transfer entropy has been proposed as a powerful tool to understand the flow of information between signals. In this paper, the comparative advantages of transfer entropy, for determining the time horizon of causal influence, are illustrated with the help of synthetic data. The technique has been specifically revised for the analysis of synchronization experiments. The investigation of experimental data from thermonuclear plasma diagnostics proves the potential and limitations of the developed approach. For complete list of authors see http://dx.doi.org/10.3390/e20090627

This paper analyses the final phase of hybrid scenario discharges at JET, the reduction of auxiliary heating towards finally the Ohmic phase. The here considered Ohmic phase is mostly still in the current flattop but may also be in the current ramp down. For this purpose a database is created of 54 parameters in 7 phases distributed in time of the discharge. It is found that the occurrence of a locked mode is in most cases preceded by a radiation peaking after the main heating phase either in a low power phase and/or in the Ohmic phase. To gain insight on the importance of different parameters in this process a correlation analysis to the radiation peaking in the Ohmic phase is done. The first finding is that the further away in time the analysed phases are the less the correlation is. This means in the end that a good termination scenario might also be able to terminate unhealthy plasmas safely. The second finding is that remaining impurities in the plasma after reducing the heating power in the termination phase are the most important reason for generating a locked mode which can lead to a disruption. © EURATOM 2018.

A novel Shutdown dose rate benchmark experiment has been performed at Joint European Torus (JET) machine during the last high performance Deuterium-Deuterium (DD) campaign in preparation of future ...

In order to contribute to the understanding of runaway electron generation mechanisms during tokamak disruptions, a test particle tracker is introduced in the JOREK 3D non-linear MHD code, able to compute both full and guiding center relativistic orbits. Tests of the module show good conservation of the invariants of motion and consistency between full orbit and guiding center solutions. A first application is presented where test electron confinement properties are investigated in a massive gas injection-triggered disruption simulation in JET-like geometry. It is found that electron populations initialised before the thermal quench (TQ) are typically not fully deconfined in spite of the global stochasticity of the magnetic field during the TQ. The fraction of ‘survivors’ decreases from a few tens down to a few tenths of percent as the electron energy varies from 1 keV to 10 MeV. The underlying mechanism for electron ‘survival’ is the prompt reformation of closed magnetic surfaces at the plasma core and, to a smaller extent, the subsequent reappearance of a magnetic surface at the edge. It is also found that electrons are less deconfined at 10 MeV than at 1 MeV, which appears consistent with a phase averaging effect due to orbit shifts at high energy.

Heavy impurities, such as tungsten (W), can exhibit strongly poloidally asymmetric density profiles in rotating or radio frequency heated plasmas. In the metallic environment of JET, the poloidal asymmetry of tungsten enhances its neoclassical transport up to an order of magnitude, so that neoclassical convection dominates over turbulent transport in the core. Accounting for asymmetries in neoclassical transport is hence necessary in the integrated modeling framework. The neoclassical drift kinetic code, NEO [E. Belli and J. Candy, Plasma Phys. Controlled Fusion P50, 095010 (2008)], includes the impact of poloidal asymmetries on W transport. However, the computational cost required to run NEO slows down significantly integrated modeling. A previous analytical formulation to describe heavy impurity neoclassical transport in the presence of poloidal asymmetries in specific collisional regimes [C. Angioni and P. Helander, Plasma Phys. Controlled Fusion 56, 124001 (2014)] is compared in this work to numerical results from NEO. Within the domain of validity of the formula, the factor for reducing the temperature screening due to poloidal asymmetries had to be empirically adjusted. After adjustment, the modified formula can reproduce NEO results outside of its definition domain, with some limitations: When main ions are in the banana regime, the formula reproduces NEO results whatever the collisionality regime of impurities, provided that the poloidal asymmetry is not too large. However, for very strong poloidal asymmetries, agreement requires impurities in the Pfirsch-Schlüter regime. Within the JETTO integrated transport code, the analytical formula combined with the poloidally symmetric neoclassical code NCLASS [W. A. Houlberg et al., Phys. Plasmas 4, 3230 (1997)] predicts the same tungsten profile as NEO in certain cases, while saving a factor of one thousand in computer time, which can be useful in scoping studies. The parametric dependencies of the temperature screening reduction due to poloidal asymmetries would need to be better characterised for this faster model to be extended to a more general applicability.

The impact of the finite grid size in SOLPS-ITER edge plasma simulations is assessed for JET H-mode discharges with a metal wall. For a semi-horizontal divertor configuration it is shown that the separatrix density is at least 30% higher when a narrow scrape-off layer (SOL) grid width is chosen in SOLPS-ITER compared to the case for which the SOL grid width is maximised. The density increase is caused by kinetic neutrals being not confined inside the divertor region because of the reduced extent of the plasma grid. In this case, an enhanced level of reflections of energetic neutrals at the low-field side (LFS) metal divertor wall is observed. This leads to a shift of the ionisation source further upstream which must be accounted for as a numerical artefact. An overestimate in the cooling at the divertor entrance is observed in this case, identified by a reduced heat flux decay parameters lambda(div)(q). Otherwise and further upstream the mid-plane heat decay length lambda(q) parameter is not affected by any change in divertor dissipation. This confirms the assumptions made for the ITER divertor design studies, i.e. that lambda(q) upstream is essentially set by the assumptions for the ratio radial to parallel heat conductivity. It is also shown that even for attached conditions the decay length relations lambda(ne)>lambda(Te)>lambda(q) hold in the near-SOL upstream. Thus for interpretative edge plasma simulations one must take the (experimental) value of lambda(ne) into account, rather than lambda(q), as the former actually defines the required minimum upstream SOL grid extent. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2018.10.013

ASCOT Fusion Source Integrator AFSI, an efficient tool for calculating fusion reaction rates and characterizing the fusion products, based on arbitrary reactant distributions, has been developed and is reported in this paper. Calculation of reactor-relevant D-D, D-T and D-3He fusion reactions has been implemented based on the Bosch-Hale fusion cross sections. The reactions can be calculated between arbitrary particle populations, including Maxwellian thermal particles and minority energetic particles. Reaction rate profiles, energy spectra and full 4D phase space distributions can be calculated for the non-isotropic reaction products. The code is especially suitable for integrated modelling in self-consistent plasma physics simulations as well as in the Serpent neutronics calculation chain. Validation of the model has been performed for neutron measurements at the JET tokamak and the code has been applied to predictive simulations in ITER.

The kinetic code for plasma periphery (KIPP) was used to assess the importance of the kinetic effects of parallel electron transport in the scrape-off layer (SOL) and divertor of JET high radiative H-mode inter-ELM plasma conditions with the ITER-like wall and strong nitrogen (N2) injection. Plasma parameter profiles along a magnetic field from one of the EDGE2D-EIRENE simulation cases were used as an input for KIPP runs. Profiles were maintained by particle and power sources. KIPP generated electron distribution functions, f e, parallel power fluxes, electron-ion thermoforces, Debye sheath potential drops and electron sheath transmission factors at divertor targets. For heat fluxes in the main SOL, KIPP results showed deviations from classical (e.g. Braginskii) fluxes by factors typically of ∼1.5, sometimes up to 2, with the flux limiting for more upstream positions and flux enhancement near entrances to the divertor. In the divertor, at the same time, for radial positions closer to the separatrix, very large heat flux enhancement factors of up to ten or even higher, indicative of a strong nonlocal heat transport, were found at the outer target, with heat power flux density exhibiting bump-on-tail features at high energies. Under such extreme conditions, however, contributions of conductive power fluxes to total power fluxes were strongly reduced, with convective power fluxes becoming comparable, or sometimes exceeding, conductive power fluxes. Electron-ion thermoforce, on the other hand, which is known to be determined mostly by thermal and subthermal electrons, was found to be in good agreement with Braginskii formulas, including the Z eff dependence. Overall, KIPP results indicate, at least for the plasma conditions used in this modelling, a sizable, but not dominant, effect of kinetics on parallel electron transport.

We report the first experimental observation of the depolarization of the Thomson scattering (TS) radiation, a relativistic effect expected to occur in very high T e plasmas and never observed so far in a fusion machine. A set of unused optical fibers in the collection optics of the high resolution Thomson scattering system of JET has been used to detect the depolarized TS radiation during a JET campaign with T e 8keV . A linear polarizer with the axis perpendicular to the direction of the incident E -field was placed in front of a fiber optic pair observing a region close to the plasma core, while another fiber pair with no polariser simultaneously observed an adjacent plasma region. The measured intensity ratio was found to be consistent with the theory, taking into account sensitivity coefficients of the two measurement channels determined with post-experiment calibrations and Raman scattering. This depolarization effect is at the basis of polarimetric TS, a different and complementary method for the analysis of TS spectra that can provide significant advantages for T e measurements in very hot plasmas such as in ITER ( T e 40keV ).

A new compact gamma-ray spectrometer was developed in order to optimise the measurement of bremsstrahlung radiation emitted from runaway electrons in the MeV range. The detector is based on a cerium doped lutetium-yttrium oxyorthosilicate (LYSO:Ce) scintillator coupled to a silicon photomultiplier and is insensitive to magnetic fields. A dedicated electronic board was developed to optimise the signal readout as well as for online control of the device. The detector combines a dynamic range up to 10 MeV with moderate energy non-linearity, counting rate capabilities in excess of 1 MHz, and an energy resolution that extrapolates to a few % in the MeV range, thus meeting the requirements for its application to runaway electron studies by bremsstrahlung measurements in the gamma-ray energy range.

Neutron emission spectroscopy is a diagnostic technique that allows for energy measurements of neutrons born in nuclear reactions. The JET tokamak fusion experiment (Culham, UK) has a special role in this respect as advanced spectrometers for 2.5 MeV and 14 MeV neutrons have been developed here for the first time for measurements of the neutron emission spectrum from D and DT plasmas with unprecedented accuracy. Twin liquid scintillating neutron spectrometers were built and calibrated at the Physikalisch-Technische Bundesanstalt (PTB) (Braunschweig, Germany) and installed on JET in the recent years with tangential-equatorial (KM12) and vertical-radial (KM13) view lines, with the latter only recently operational. This article reports on the performance of KM12 and on the development of the data analysis methods in order to extract physics information upon D ions kinematics in JET auxiliary-heated D plasmas from 2.5 MeV neutron measurements. The comparison of these results with the correspondents from other JET neutron spectrometers is also presented: their agreement allows for JET unique capability of multi-lines of sight neutron spectroscopy and for benchmarking other 14 MeV neutron spectrometers installed on the same lines of sight in preparation for the DT experimental campaign at JET. © 2018 EURATOM.

Compact DT neutron generators based on accelerators are often built on the principle of a mixed beam operation, meaning that deuterium (D) and tritium (T) are both present in the ion beam and in the target. Moreover, the beam consists of a mixture of ions and ionized molecules (D, T ions, and ionized D-D, T-T and D-T molecules) so the relevant source components come from T(d, n), D(t, n), D(d, n) and T(t, 2n) reactions at different ion energies. The method for assessing the relative amplitudes of different source components (DD, DT, TT) is presented. The assessment relies on the measurement of the neutron spectrum of different DT components (T(d, n) and D(t, n) at different energies) using a high resolution neutron spectrometer, e.g. a diamond detector, fusion reaction cross-sections, and simulations of neutron generation in the target. Through this process a complete description of the neutron source properties of the mixed beam neutron generator can be made and a neutron source description card, in a format suitable for Monte Carlo code MCNP, produced.

Neutral particle analysis (NPA) is one of the few methods for diagnosing fast ions inside a plasma by measuring neutral atom fluxes emitted due to charge exchange reactions. The JET tokamak features an NPA diagnostic which measures neutral atom fluxes and energy spectra simultaneously for hydrogen, deuterium and tritium species. A synthetic NPA diagnostic has been developed and used to interpret these measurements to diagnose energetic particles in JET plasmas with neutral beam injection (NBI) heating. The synthetic NPA diagnostic performs a Monte Carlo calculation of the neutral atom fluxes in a realistic geometry. The 4D fast ion distributions, representing NBI ions, were simulated using the Monte Carlo orbit-following code ASCOT. Neutral atom density profiles were calculated using the FRANTIC neutral code in the JINTRAC modelling suite. Additionally, for rapid analysis, a scan of neutral profiles was precalculated with FRANTIC for a range of typical plasma parameters. These were taken from the JETPEAK database, which includes a comprehensive set of data from the flat-top phases of nearly all discharges in recent JET campaigns. The synthetic diagnostic was applied to various JET plasmas in the recent hydrogen campaign where different hydrogen/deuterium mixtures and NBI configurations were used. The simulated neutral fluxes from the fast ion distributions were found to agree with the measured fluxes, reproducing the slowing-down profiles for different beam isotopes and energies and quantitatively estimating the fraction of hydrogen and deuterium fast ions.

The ASCOT Fusion Source Integrator (AFSI) has been used to calculate neutron production rates and spectra corresponding to the JET 19-channel neutron camera (KN3) and the time-of-flight spectrometer (TOFOR) as ideal diagnostics, without detector-related effects. AFSI calculates fusion product distributions in 4D, based on Monte Carlo integration from arbitrary reactant distribution functions. The distribution functions were calculated by the ASCOT Monte Carlo particle orbit following code for thermal, NBI and ICRH particle reactions. Fusion cross-sections were defined based on the Bosch-Hale model and both DD and DT reactions have been included. Neutrons generated by AFSI-ASCOT simulations have already been applied as a neutron source of the Serpent neutron transport code in ITER studies. Additionally, AFSI has been selected to be a main tool as the fusion product generator in the complete analysis calculation chain: ASCOT - AFSI - SERPENT (neutron and gamma transport Monte Carlo code) - APROS (system and power plant modelling code), which encompasses the plasma as an energy source, heat deposition in plant structures as well as cooling and balance-of-plant in DEMO applications and other reactor relevant analyses. This conference paper presents the first results and validation of the AFSI DD fusion model for different auxiliary heating scenarios (NBI, ICRH) with very different fast particle distribution functions. Both calculated quantities (production rates and spectra) have been compared with experimental data from KN3 and synthetic spectrometer data from ControlRoom code. No unexplained differences have been observed. In future work, AFSI will be extended for synthetic gamma diagnostics and additionally, AFSI will be used as part of the neutron transport calculation chain to model real diagnostics instead of ideal synthetic diagnostics for quantitative benchmarking.

The fuel outgassing efficiency from plasma-facing components exposed in JET-ILW has been studied at ITER-relevant baking temperatures. Samples retrieved from the W divertor and Be main chamber were annealed at 350 and 240 °C, respectively. Annealing was performed with thermal desoprtion spectrometry (TDS) for 0, 5 and 15 h to study the deuterium removal effectiveness at the nominal baking temperatures. The remained fraction was determined by emptying the samples fully of deuterium by heating W and Be samples up to 1000 and 775 °C, respectively. Results showed the deposits in the divertor having an increasing effect to the remaining retention at temperatures above baking. Highest remaining fractions 54 and 87% were observed with deposit thicknesses of 10 and 40 μm, respectively. Substantially high fractions were obtained in the main chamber samples from the deposit-free erosion zone of the limiter midplane, in which the dominant fuel retention mechanism is via implantation: 15 h annealing resulted in retained deuterium higher than 90%. TDS results from the divertor were simulated with TMAP7 calculations. The spectra were modelled with three deuterium activation energies resulting in good agreement with the experiments EURATOM 633053

Future devices like JT-60SA, ITER and DEMO require quantitative predictions of pedestal density and temperature levels, as well as inter-ELM and ELM divertor heat fluxes, in order to improve global confinement capabilities while preventing divertor erosion/melting in the planning of future experiments. Such predictions can be obtained from dedicated pedestal models like EPED, and from non-linear MHD codes like JOREK, for which systematic validation against current experiments is necessary. In this paper, we show progress in the quantitative validation of the JOREK code using JET simulations. Results analyse the impact of diamagnetic terms on the dynamics and size of the ELMs, and evidence is provided that the onset of type-I ELMs is not governed by linear MHD stability alone, but that a nonlinear threshold could be responsible for large MHD events at the plasma edge. This work has been carried out within the framework of the EUROfusion Consortium and has received funding from the Euratom research and training programme 2014–2018 under grant agreement No 633053, and from the RCUK Energy Programme (grant number EP/I501045). To obtain further information on the data and models underlying this paper please contact PublicationsManagerccfe.ac.uk. This work used the HELIOS supercomputer (IFERC-CSC), Japan, under the Broader Approach collaboration, implemented by Fusion for Energy and JAEA. The views and opinions expressed herein do not necessarily reflect those of the European Commission or the ITER Organization. The HEC ARCHER computer (UK), as part of the Plasma HEC Consortium EPSRC grant EP/L000237/1, and the MARCONI computer at CINECA in Italy, were also used.

In metallic machines ICRH heating is playing an increasingly important role. One of its most recent applications on the Joint Europena Torus (JET) is sawtooth control by ICRH modulation, for avoiding triggering dangerous neo-classical tearing modes (NTMs) and counteracting impurity accumulation. Some of the main difficulties of these experiments are the assessment of the synchronization efficiency and the understanding of the main physical mechanisms at play. In this paper, three independent classes of statistical indicators are introduced to address these issues: Recurrence Plots, Convergent Cross Mapping and Transfer Entropy. The application to JET experiments with the ILW shows that the proposed indicators agree quite well among themselves and provide sound estimates of the efficiency of the synchronisation scheme investigated. They also support, with a shot to shot basis analysis and an estimate of the uncertainties, the interpretation that the fast ions play a fundamental role in the stabilization of the sawteeth, in both L and H mode. Proposals for experiments to be carried out in the future to consolidate the interpretation of the results are discussed.

In this work, an investigation of the neutral gas flow in the JET sub-divertor area is presented, with respect to the interaction between the plasma side and the pumping side. The edge plasma side is simulated with the SOLPS code, while the sub-divertor area is modeled by means of the direct simulation Monte Carlo (DSMC) method, which in the last few years has proved well able to describe rarefied, collisional flows in tokamak sub-divertor structures. Four different plasma scenarios have been selected, and for each of them a user-defined, iterative procedure between SOLPS and DSMC has been established, using the neutral flux as the key communication term between the two codes. The goal is to understand and quantify the mutual influence between the two regions in a self-consistent manner, that is to say, how the particle exhaust pumping system controls the upstream plasma conditions. Parametric studies of the flow conditions in the sub-divertor, including additional flow outlets and variations of the cryopump capture coefficient, have been performed as well, in order to understand their overall impact on the flow field. The DSMC analyses resulted in the calculation of both the macroscopic quantities—i.e. temperature, number density and pressure—and the recirculation fluxes towards the plasma chamber. The consistent values for the recirculation rates were found to be smaller than those according to the initial standard assumption made by SOLPS.

The JET ITER-Like Wall experiment, with its all-metal plasma-facing components, provides a unique environment for plasma and plasma-wall interaction studies. These studies are of great importance in understanding the underlying phenomena taking place during the operation of a future fusion reactor. Present work summarizes and reports the plasma fuel retention in the divertor resulting from the two first experimental campaigns with the ITER-Like Wall. The deposition pattern in the divertor after the second campaign shows same trend as was observed after the first campaign: highest deposition of 10-15 μm was found on the top part of the inner divertor. Due to the change in plasma magnetic configurations from the first to the second campaign, and the resulted strike point locations, an increase of deposition was observed on the base of the divertor. The deuterium retention was found to be affected by the hydrogen plasma experiments done at the end of second experimental campaign.

The isotope content of the plasma reflects on the dynamics of isotope changeover experiments, efficiency of wall conditioning and the performance of a fusion device in the active phase of operation ...

ERO is a Monte-Carlo code for modeling plasma-wall interaction and 3D plasma impurity transport for applications in fusion research. The code has undergone a significant upgrade (ERO2.0) which allows increasing the simulation volume in order to cover the entire plasma edge of a fusion device, allowing a more self-consistent treatment of impurity transport and comparison with a larger number and variety of experimental diagnostics. In this contribution, the physics-relevant technical innovations of the new code version are described and discussed. The new capabilities of the code are demonstrated by modeling of beryllium (Be) erosion of the main wall during JET limiter discharges. Results for erosion patterns along the limiter surfaces and global Be transport including incident particle distributions are presented. A novel synthetic diagnostic, which mimics experimental wide-angle 2D camera images, is presented and used for validating various aspects of the code, including erosion, magnetic shadowing, non-local impurity transport, and light emission simulation. EURATOM 633053

In this work, discriminant analysis is used as the main approach for building a physics based automated classifier for the discrimination of the edge-localized mode (ELM) plasma instability. The classifier is then applied for distinguishing type I and type III ELMs from a set of carbon-wall plasmas at JET. This provides a fast, standardized classification of ELM types which is expected to significantly reduce the effort of ELM experts in identifying ELM types. Further, the classifier yields a separation hyperplane in terms of global plasma parameters, which provides an insight into the range of conditions under which specific ELM behaviors occur. EURATOM 633053

Since 2011 the JET tokamak has been operated with a metal ITER-like wall (JET-ILW) including castellated beryllium limiters and lamellae-type bulk tungsten tiles in the divertor. This has allowed for a large scale test of castellated plasma-facing components (PFC). Procedures for sectioning the limiters into single blocks of castellation have been developed. This facilitated morphology studies of morphology of surfaces inside the grooves for limiters after experimental campaigns 2011–2012 and 2013–2014. The deposition in the 0.4–0.5mm wide grooves of the castellation is ‘shallow’. It reaches 1–2mm into the 12mm deep gap. Deuterium concentrations are small (mostly below 1 × 1018 cm−2 ). The estimated total amount of deuterium in all the castellated limiters does not exceed the inventory of the plasma-facing surfaces (PFS) of the limiters. There are only traces of Ni, Cr and Fe deposited in the castellation gaps. The same applies to the carbon content. Also low deposition of D, Be and C has been measured on the sides of the bulk tungsten lamellae pieces. Modelling clearly reflects: (a) a sharp decrease in the measured deposition profiles and (b) an increase in deposition with the gap width. Both experimental and modelling data give a strong indication and information to ITER that narrow gaps in the castellated PFC are essential. X-ray diffraction on PFS has clearly shown two distinct composition patterns: Be with an admixture of Be–W intermetallic compounds (e.g. Be22W) in the deposition zone, whilst only pure Be has been detected in the erosion zone. The lack of compound formation in the erosion zone indicates that no distinct changes in the thermo-mechanical properties of the Be PFC might be expected. EURATOM 633053 Swedish Research Council (VR) 2015–04844

Recent JET (ITER-Like Wall) experiments have shown that the fueling gas puffed from different locations of the vessel can result in different scrape-off layer (SOL) density profiles and therefore different radio frequency (RF) coupling. To reproduce the experimental observations, to understand the associated physics and to optimize the gas puff methods, we have carried out three-dimensional (3D) simulations with the EMC3-EIRENE code in JET-ILW including a realistic description of the vessel geometry and the gas injection modules (GIMs) configuration. Various gas puffing methods have been investigated, in which the location of gas fueling is the only variable parameter. The simulation results are in quantitative agreement with the experimental measurements. They confirm that compared to divertor gas fueling, mid-plane gas puffing increases the SOL density most significantly but locally, while top gas puffing increases it uniformly in toroidal direction but to a lower degree. Moreover, the present analysis corroborates the experimental findings that combined gas puff scenarios—based on distributed main chamber gas puffing—can be effective in increasing the RF coupling for multiple antennas simultaneously. The results indicate that the spreading of the gas, the local ionization and the transport of the ionized gas along the magnetic field lines connecting the local gas cloud in front of the GIMs to the antennas are responsible for the enhanced SOL density and thus the larger RF coupling. EURATOM 633053

A new compact gamma-ray spectrometer based on a Silicon Photo-Multiplier (SiPM) coupled to a LaBr3(Ce) crystal has been developed for the upgrade of the Gamma Camera (GC) of JET, where it must operate in a high intensity neutron/gamma-ray admixed field. The work presents the results of an experiment aimed at characterizing the effect of 14 MeV neutron irradiation on both LaBr3(Ce) and SiPM that compose the full detector. The pulse height spectrum from neutron interactions with the crystal has been measured and is successfully reproduced by MCNP simulations. It is calculated that about 8% of the impinging neutrons leave a detectable signal of which less than < 4% of the events occur in the energy region above 3 MeV, which is of interest for gamma-ray spectroscopy applications. Neutron irradiation also partly degrades the performance of the SiPM and this is mostly manifested as an increase of the dark current versus the neutron fluence. However, it was found that the SiPM can be still operated up to a fluence of 4×1010 n/cm2, which is the highest value we experimentally tested. Implications of these results for GC measurements at JET are discussed. © 2017 Istituto di Fisica del Plasma - CNR. Published by IOP Publishing Ltd on behalf of Sissa Medialab.

In the past at JET, with the MkI divertor, a systematic study of the influence of X-point height and poloidal flux expansion has been conducted [1,2] showing minor differences in the radiation distribution, whereas in [3] experiment and simulations have shown enhancement of detachment as the flux expansion was increased. More recently at JET, equipped with the ITER-like wall (ILW), radiative seeded scenarios have been studied and a maximum radiation fraction 75% has been achieved. EDGE2D-EIRENE [4-6] simulations [7,8] have already shown that the divertor heat fluxes can be reduced with N2 injection, qualitatively consistent with experimental observations [9], by adjusting the impurity injection rate to reproduce the measured divertor radiation. In this paper we will present edge predictive simulations on the assessment of effects of poloidal flux expansion and recycling on radiation distribution and X-point peaking on JET-ILW nitrogen seeded plasmas. (C) 2017 The Authors. Published by Elsevier Ltd. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2017.07.004

To understand the fuel retention mechanism correlation of surface chemical states and hydrogen isotope retention behavior determined by XPS (X-ray photoelectron spectroscopy) and TDS (Thermal desorption spectroscopy), respectively, for JET ITER-Like Wall samples from operational period 2011–2012 were investigated. It was found that the deposition layer was formed on the upper part of the inner vertical divertor area. At the inner plasma strike point region, the original surface materials, W or Mo, were found, indicating to an erosion-dominated region, but deposition of impurities was also found. Higher heat load would induce the formation of metal carbide. At the outer horizontal divertor tile, mixed material layer was formed with iron as an impurity. TDS showed the H and D desorption behavior and the major D desorption temperature for the upper part of the inner vertical tile was located at 370 °C and 530 °C. At the strike point region, the D desorption temperature was clearly shifted toward higher release temperatures, indicating the stabilization of D trapping by higher heat load

A calibration of the JET neutron detectors was carried out prior to the upcoming deuterium-tritium experimental campaign. Two Compact DT neutron generators (NGs) were purchased for this purpose from VNIIA, Russia. These generators are capable of producing approximately 2 × 108 neutrons/s with a DT fusion energy spectrum. Preceding the in-vessel calibration measurements, these compact generators were tested and fully characterised at the UK's National Physical Laboratory (NPL). In order to support the characterisation measurements, detailed neutronics models were developed of the NGs, monitoring detectors and remote handling (RH) apparatus. Neutron spectra calculated from these models have been used to help determine NPL long counter efficiencies and effective centres, as well as NPL reference iron and aluminium activation foil reaction rates. The neutron emission rate has been measured for both generators as a function of angle using absolutely calibrated long counters and the relative emission rate by monitoring single crystal diamond detectors. The measured anisotropy profile is shown to be reproducible with a detailed NG MCNP model. Consequently, the neutron source routine and the MCNP model of the NGs can be reliably used for the analysis of the in-vessel calibration at JET.

At the Joint European Torus (JET) the ex-vessel fission chambers and in-vessel activation detectors are used as the neutron production rate and neutron yield monitors respectively. In order to ensure that these detectors produce accurate measurements they need to be experimentally calibrated. A new calibration of neutron detectors to 14 MeV neutrons, resulting from deuterium–tritium (DT) plasmas, is planned at JET using a compact accelerator based neutron generator (NG) in which a D/T beam impinges on a solid target containing T/D, producing neutrons by DT fusion reactions. This paper presents the analysis that was performed to model the neutron source characteristics in terms of energy spectrum, angle–energy distribution and the effect of the neutron generator geometry. Different codes capable of simulating the accelerator based DT neutron sources are compared and sensitivities to uncertainties in the generator's internal structure analysed. The analysis was performed to support preparation to the experimental measurements performed to characterize the NG as a calibration source. Further extensive neutronics analyses, performed with this model of the NG, will be needed to support the neutron calibration experiments and take into account various differences between the calibration experiment and experiments using the plasma as a source of neutrons.

International audience; The Materials Detritiation Facility has been designed to thermally treat solid non-combustible radioactive waste produced during operations of the Joint European Torus (JET) that is classified as Intermediate Level Waste in the UK due to its tritium inventory (>12kBqg). The waste will be thermally treated in a retort furnace at temperatures up to 1000DC under a flowing air atmosphere for approximately 6 hours to reduce its tritium inventory sufficiently to allow its disposal at a lower waste category via existing disposal routes. The gaseous flow from the furnace will be processed via a bubbler system, where released tritium will be trapped in water.Commissioning of the facility will be divided into two main parts inactive and active. The main purpose of the inactive commissioning is to verify that all components and safety systems of the facility are installed, tested and operated properly and within their operational limits. Several trials of the furnace with non-radioactive materials will be performed to verify its temperature profile, and to verify operation of the gaseous process line.During the active commissioning, small amounts of tritium-contaminated material will be introduced into the facility and used for active trials. The tritium inventory in this material has been selected based on the As low as reasonably practicable (ALARP) principle, to ensure that the activity levels are sufficient to fully test the control instrumentation and pose minimal risk to operators during commissioning. Overall, four active trials will be performed with carbon-based and Inconel materials with total tritium inventories of 1MBq, 3GBq, 20GBq and 26GBq. Tritium levels in the bubblers as well as in aerial discharge from the facility will be monitored. Furthermore, all materials used in the active trials will be sampled and analyzed to verify the performance of the process and confirm that a major part of tritium inventory can be removed from materials by the process.

In order to extend the investigation of the plasma cleaning of ITER first mirrors, a set of molybdenum mirrors was coated in a laboratory with ITER-relevant contaminants, namely beryllium and tungsten. Different coating techniques as well as several contaminant compositions were used to ensure a large variety of films to clean, completing a previous study conducted on mirrors exposed in the JET ITER-like wall (tokamak deposits) [1]. Due to the toxicity of beryllium, the samples were treated in a vacuum chamber specially built for this purpose. The cleaning was performed using capacitively coupled RF plasma and evaluated by performing reflectivity measurements, scanning electron microscopy, x-ray photoelectron spectroscopy and ion beam analysis. The removal of all types of contaminants was achieved by using different plasma compositions (argon, helium and mixtures of the two) with various ion energies (from 200-600 eV) and in some cases the mirror's reflectivity was restored towards initial values. Pure helium discharges were capable of removing mixed beryllium/tungsten layers and oxidized molybdenum. In addition, no significant increase in the diffuse reflectivity of the mirrors was observed for the helium cleaning, though this was the case for some samples cleaned with argon. Helium is therefore appropriate for cleaning all mirrors in ITER leading to a possible cleaning regime where the entire vessel is filled with He and all mirrors are cleaned simultaneously without damaging their surfaces.

One crucial requirement for the development of fusion power is to know where, and how much, impurities collect in the machine, and how much of the fuelling isotope tritium will be trapped therein. ...

Halo currents generated during unmitigated tokamak disruptions are known to develop rotating asymmetric features that are of great concern to ITER because they can dynamically amplify the mechanical stresses on the machine. This paper presents a multi-machine analysis of these phenomena. More specifically, data from C-Mod, NSTX, ASDEX Upgrade, DIII-D, and JET are used to develop empirical scalings of three key quantities: (1) the machinespecific minimum current quench time, τCQ; (2) the halo current rotation duration, trot; and (3) the average halo current rotation frequency, fh . These data reveal that the normalized rotation duration, trot/τCQ, and the average rotation velocity, vh , are surprisingly consistent from machine to machine. Furthermore, comparisons between carbon and metal wall machines show that metal walls have minimal impact on the behavior of rotating halo currents. Finally, upon projecting to ITER, the empirical scalings indicate that substantial halo current rotation above fh = 20 Hz is to be expected. More importantly, depending on the projected value of τCQ in ITER, substantial rotation could also occur in the resonant frequency range of 6–20 Hz. As such, the possibility of damaging halo current rotation during unmitigated disruptions in ITER cannot be ruled out. EURATOM 633053 RCUK Energy Programme EP/ I501045 Princeton University DE-AC02-09CH11466

Pellet injection is a likely fuelling method of reactor grade plasmas. When the pellet ablates, it will transiently perturb the density and temperature profiles of the plasma. This will in turn change dimensionless parameters such as a/Ln,a/LT and plasma β. The microstability properties of the plasma then changes which influences the transport of heat and particles. In this paper, gyrokinetic simulations of a JET L-mode pellet fuelled discharge are performed. The ion temperature gradient/trapped electron mode turbulence is compared at the time point when the effect from the pellet is the most pronounced with a hollow density profile and when the profiles have relaxed again. Linear and nonlinear simulations are performed using the gyrokinetic code GENE including electromagnetic effects and collisions in a realistic geometry in local mode. Furthermore, global nonlinear simulations are performed in order to assess any nonlocal effects. It is found that the positive density gradient has a stabilizing effect that is partly counteracted by the increased temperature gradient in the this region. The effective diffusion coefficients are reduced in the positive density region region compared to the intra pellet time point. No major effect on the turbulent transport due to nonlocal effects are observed.

A novel approach to uncertainty quantification in codes simulating fusion plasma, deterministic sampling (DS), is evaluated. This method uses a few carefully selected samples and can be used to pro ...

The erosion of tungsten (W), induced by the bombardment of plasma and impurity particles, determines the lifetime of plasma-facing components as well as impacting on plasma performance by the influx of W into the confined region. The screening of W by the divertor and the transport of W in the plasma determines largely the W content in the plasma core, but the W source strength itself has a vital impact on this process. The JET tokamak experiment provides access to a large set of W erosion-determining parameters and permits a detailed description of the W source in the divertor closest to the ITER one: (i) effective sputtering yields and fluxes as function of impact energy of intrinsic (Be, C) and extrinsic (Ne, N) impurities as well as hydrogenic isotopes (H, D) are determined and predictions for the tritium (T) isotope are made. This includes the quantification of intra- and inter-edge localised mode (ELM) contributions to the total W source in H-mode plasmas which vary owing to the complex flux compositions and energy distributions in the corresponding phases. The sputtering threshold behaviour and the spectroscopic composition analysis provides an insight in the dominating species and plasma phases causing W erosion. (ii) The interplay between the net and gross W erosion source is discussed considering (prompt) re-deposition, thus, the immediate return of W ions back to the surface due to their large Larmor radius, and surface roughness, thus, the difference between smooth bulk-W and rough W-coating components used in the JET divertor. Both effects impact on the balance equation of local W erosion and deposition. (iii) Post-mortem analysis reveals the net erosion/deposition pattern and the W migration paths over long periods of plasma operation identifying the net W transport to remote areas. This W transport is related to the divertor plasma regime, e.g. attached operation with high impact energies of impinging particles or detached operation, as well as to the applied magnetic configuration in the divertor, e.g. close divertor with good geometrical screening of W or open divertor configuration with poor screening. JET equipped with the ITER-like wall (ILW) provided unique access to the net W erosion rate within a series of 151 subsequent H-mode discharges (magnetic field: T, plasma current: MA, auxiliary power MW) in one magnetic configuration accumulating 900 s of plasma with particle fluences in the range of 5– in the semi-detached inner and attached outer divertor leg. The comparison of W spectroscopy in the intra-ELM and inter-ELM phases with post-mortem analysis of W marker tiles provides a set of gross and net W erosion values at the outer target plate. ERO code simulations are applied to both divertor leg conditions and reproduce the erosion/deposition pattern as well as confirm the high experimentally observed prompt W re-deposition factors of more than 95% in the intra- and inter-ELM phase of the unseeded deuterium H-mode plasma. Conclusions to the expected divertor conditions in ITER as well as to the JET operation in the DT plasma mixture are drawn on basis of this unique benchmark experiment.

Experiments at JET showed locally enhanced, asymmetric beryllium (Be) erosion at outer wall limiters when magnetically connected ICRH antennas were in operation. A first modeling effort using the 3D erosion and scrape-off layer impurity transport modeling code ERO reproduced qualitatively the experimental outcome. However, local plasma parameters - in particular when 3D distributions are of interest - can be difficult to determine from available diagnostics and so erosion / impurity transport modeling input relies on output from other codes and simplified models, increasing uncertainties in the outcome. In the present contribution, we introduce and evaluate the impact of improved models and parameters with largest uncertainties of processes that impact impurity production and transport across the scrape-off layer, when simulated in ERO: (i) the magnetic geometry has been revised, for affecting the separatrix position (located 50-60 mm away from limiter surface) and thus the background plasma profiles; (ii) connection lengths between components, which lead to shadowing of ion fluxes, are also affected by the magnetic configuration; (iii) anomalous transport of ionized impurities, defined by the perpendicular diffusion coefficient, has been revisited; (iv) erosion yields that account for energy and angular distributions of background plasma ions under the present enhanced sheath potential and oblique magnetic field, have been introduced; (v) the effect of additional erosion sources, such as charge-exchange neutral fluxes, which are dominant in recessed areas like antennas, has been evaluated; (vi) chemically assisted release of Be in molecular form has been included. Sensitivity analysis highlights a qualitative effect (i.e. change in emission patterns) of magnetic shadowing, anomalous diffusion, and inclusion of neutral fluxes and molecular release of Be. The separatrix location, and energy and angular distribution of background plasma fluxes impact erosion quantitatively. ERO simulations that include all features described above match experimentally measured Be I (457.3 nm) and Be II (467.4 nm) signals, and erosion increases with varying ICRH antenna's RF power. However, this increase in erosion is only partially captured by ERO's emission measurements, as most contributions from plasma wetted surfaces fall outside the volume observed by sightlines.

Équipe 107 : Physique des plasmas chauds; International audience; The ITER baseline scenario, with 500 MW of DT fusion power and Q = 10, will rely on a Type I ELMy H-mode, with Delta W = 0.7 MJ mitigated edge localized modes (ELMs). Tungsten (W) is the material now decided for the divertor plasma-facing components from the start of plasma operations. W atoms sputtered from divertor targets during ELMs are expected to be the dominant source under the partially detached divertor conditions required for safe ITER operation. W impurity concentration in the plasma core can dramatically degrade its performance and lead to potentially damaging disruptions. Understanding the physics of plasma-wall interaction during ELMs is important and a primary input for this is the energy of incoming ions during an ELM event. In this paper, coupled Infrared thermography and Langmuir Probe (LP) measurements in JET-ITER-Like-Wall unseeded H-mode experiments with ITER relevant ELM energy drop have been used to estimate the impact energy of deuterium ions (D+) on the divertor target. This analysis gives an ion energy of several keV during ELMs, which makes D+ responsible for most of the W sputtering in unseeded H-mode discharges. These LP measurements were possible because of the low electron temperature (T-e) during ELMs which allowed saturation of the ion current. Although at first sight surprising, the observation of low T-e at the divertor target during ELMs is consistent with the `Free-Streaming' kinetic model which predicts a near-complete transfer of parallel energy from electrons to ions in order to maintain quasi-neutrality of the ELM filaments while they are transported to the divertor targets.

This paper presents a diagnostic system, gamma-ray alpha-particle monitor (GRAM), for continuous monitoring of deuterium-tritium fusion alpha-particles in the MeV energy range escaped from the plas ...