Search

Accurate measurement of the nuclear environment within a test tritium breeding-blanket module of a fusion reactor is crucial to determine tritium production rates which are relevant to self-suffici ...

Numerical scan at constant shows that core and scrape-off layer (SOL) radiations do not depend on the plasma current (I-p). Whereas the SOL radiation increases with seeding, the core radiation, how ...

Activities under the EURO fusion work package JET3 programme have been established to enable the technological exploitation of the planned JET experiments over the next few years, which culminates in a D-T experimental campaign, DTE-2. In the areas of nuclear technology and nuclear safety the programme offers a unique opportunity to provide experimental data that is relevant to ITER. The key purpose of the collected data will be to support bench marking and validation activities relating to neutronics and activation codes, and associated nuclear data, that are used to predict the nuclear behavior of ITER component and materials, during and after operations. This paper details the status and key issues of the ongoing ACT sub-project under work package JET3, which aims to take advantage of the large 14 MeV neutron fluence expected during JET DTE-2 to irradiate samples of real ITER materials used in the manufacturing of the main in-vessel tokamak components. The materials considered, with specified minor elemental impurity levels, include: Nb3Sn, SS316L steels from a range of manufacturers, SS304B, Alloy 660, W, CuCrZr, OF-Cu, XM-19, Al bronze, NbTi and EUROFER. The activities include provision for measurement of nuclide activities for each material and comparison against the predicted quantities through calculation with the FISPACT-II inventory code. Included here are key pre-analysis results for the selected ITER irradiation samples, and corresponding optimization of dosimetry foils (Ti, Mn, Co, Ni, Y, Fe, Co, Sc, Ta) that will be irradiated at selected positions inside JET irradiation stations in order to determine the neutron spectrum. Preliminary experimental activation results through recent JET D–D operations are discussed.

For simulation of plasma-facing component erosion in fusion experiments, an analytical expression for the ion velocity just before the surface impact including the local electric field and an optional surface biasing effect is suggested. Energy and angular impact distributions and the resulting effective sputtering yields were produced for several experimental scenarios at JET ILW mostly involving PFCs exposed to an oblique magnetic field. The analytic solution has been applied as an improvement to earlier ERO modelling of localized, Be outer limiter, RF-enhanced erosion, modulated by toggling of a remote, however magnetically connected ICRH antenna. The effective W sputtering yields due to D and Be ion impact in Type-I and Type-III ELMs and inter-ELM conditions were also estimated using the analytical approach and benchmarked by spectroscopy. The intra-ELM W sputtering flux increases almost 10 times in comparison to the inter-ELM flux. (C) 2017 Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2017.03.031

Beryllium film deposition was studied with cavity samples in remote areas of the inner and outer JET divertor and below divertor tile 5 during the 2011-2012 campaign with the ITER-like wall. Predominantly beryllium films were formed inside the cavities with some additional carbon, the ratio Be/C was > 2. These deposited layers had high D/(Be+C) ratios of about 0.3. The formation of these films is mainly due to sticking of beryllium-containing particles with low sticking coefficients < 0.5. The observed surface loss probabilities depend on the position in the divertor. The particles responsible for film deposition originated from the location of in the divertor strike points. (C) 2016 Elsevier Ltd. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2016.12.005

In the frame of the WPJET3-DT Technology project within the EUROfusion Consortium program, neutronics experiments are in preparation for the future deuterium–tritium campaign on JET (DTE2). The experiments will be conducted with the purpose to validate the neutronics codes and tools used in ITER, thus reducing the related uncertainties and the associated risks in the machine operation. This paper summarizes the status of previous shutdown dose rate benchmarks experiments and analyses performed at JET and focuses on the computational and experimental efforts conducted in preparation of the future DTE2 experiments. In particular, preliminary calculations and studies to select detectors and positions aimed to reduce uncertainties in the shutdown dose rate experiment are presented and discussed.

A new Safety Case is required to permit tritium operations on JET during the forthcoming DTE2 campaign. The outputs, benefits and lessons learned associated with the production of this Safety Case are presented. The changes that have occurred to the Safety Case methodology since the last JET tritium Safety Case are reviewed. Consideration is given to the effects of modifications, particularly ITER related changes, made to the JET and the impact these have on the hazard assessments as well as normal operations. Several specialized assessments, including recent MELCOR modelling, have been undertaken to support the production of this Safety Case and the impact of these assessments is outlined. Discussion of the preliminary actions being taken to progress implementation of this Safety Case is provided, highlighting new methods to improve the dissemination of the key Safety Case results to the plant operators. Finally, the work required to complete this Safety Case, before the next tritium campaign, is summarized.

The experimental method developed at ASDEX Upgrade for the determination of the intrinsic tungsten (W) density profile coupling data from the soft X-ray (SXR) diagnostic and vacuum-ultra-violet (VUV) spectroscopy has been upgraded for application to JET plasmas. The strong poloidal asymmetries in the SXR emission are modeled assuming a ln(ϵ(ρ,R)/ϵ(ρ,R0))=λ(ρ)(R2−R02) distribution, where ρ is the flux coordinate, R is the major radius, and λ is the fit parameter. The W density is calculated from the resulting 2D SXR emissivity maps accounting for contributions from a low-Z impurity (typically beryllium) and main ion with the assumption that their contributions are poloidally symmetric. Comparing the result with the independent W concentration measurement of VUV spectroscopy, a recalibration factor for the SXR emissivity is calculated making the method robust against the decrease in the sensitivity of the SXR diodes which has been observed across multiple campaigns. The final 2D W density map is checked for consistency versus the time-evolution of the W concentration measurement from VUV spectroscopy, toroidal rotation measurements from charge exchange recombination spectroscopy, and tomographic reconstructions of bolometry data. The method has been found to be robust for W concentrations above a few 10−5 and in cases where the contributions from other medium-Z impurities such as Ni are negligible.The experimental method developed at ASDEX Upgrade for the determination of the intrinsic tungsten (W) density profile coupling data from the soft X-ray (SXR) diagnostic and vacuum-ultra-violet (VUV) spectroscopy has been upgraded for application to JET plasmas. The strong poloidal asymmetries in the SXR emission are modeled assuming a ln(ϵ(ρ,R)/ϵ(ρ,R0))=λ(ρ)(R2−R02) distribution, where ρ is the flux coordinate, R is the major radius, and λ is the fit parameter. The W density is calculated from the resulting 2D SXR emissivity maps accounting for contributions from a low-Z impurity (typically beryllium) and main ion with the assumption that their contributions are poloidally symmetric. Comparing the result with the independent W concentration measurement of VUV spectroscopy, a recalibration factor for the SXR emissivity is calculated making the method robust against the decrease in the sensitivity of the SXR diodes which has been observed across multiple campaigns. The final 2D W density map is checked for...

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.

There is experimental evidence that the pedestal dynamics in type-I ELMy H-mode discharges is significantly affected by a change in the recycling conditions at the tungsten plasma-facing components ...

In this paper we present the comparison of simulations with the numerical codes COREDIV and SOLPS5.0 for JET L-mode discharges with ITER like wall (ILW). The simulations have been performed for L-m ...

In this paper numerical simulations with the self-consistent COREDIV code of the planned JET DT experiments have been performed. First, record shot from the 1997 experiments was simulated and good ...

A new analytical approximation for the electric potential profile in the presence of an oblique magnetic field and the analytical solution for the particle motion just before the impact with a plas ...

The JET tokamak is operated with beryllium limiter tiles in the main chamber and tungsten coated carbon fiber composite tiles and solid W tiles in the divertor. One important issue is how wall materials are migrating during plasma operation. To study beryllium redistribution in the main chamber and in the divertor, a 10Be enriched limiter tile was installed prior to plasma operations in 2011-2012. Methods to take surface samples have been developed, an abrasive method for bulk Be tiles in the main chamber, which permits reuse of the tiles, and leaching with hot HCl to remove all Be deposited at W coated surfaces in the divertor. Quantitative analysis of the total amount of Be in cm2 sized samples was made with inductively coupled plasma atomic emission spectroscopy (ICP-AES). The 10Be/9Be ratio in the samples was measured with accelerator mass spectrometry (AMS). The experimental setup and methods are described in detail, including sample preparation, measures to eliminate contributions in AMS from the 10B isobar, possible activation due to plasma generated neutrons and effects of diffusive isotope mixing. For the first time marker concentrations are measured in the divertor deposits. They are in the range 0.4-1.2% of the source concentration, with moderate poloidal variation.

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.

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

This paper deals with a divertor coil currents optimized procedure to design High Flux Expansion (HFE) configurations in the JET tokamak aimed to study the effects of flux expansion variation on the radiation fraction and radiated power re-distribution. A number of benefits of HFE configuration have been experimentally demonstrated on TCV, EAST, NSTX and DIII-D tokamaks and are under investigation for next generation devices, as DEMO and DTT. The procedure proposed here exploits the linearized relation between the plasma-wall gaps and the Poloidal Field (PF) coil currents. Once the linearized model is provided by means of CREATE-NL code, the divertor coils currents are calculated using a constrained quadratic programming optimization procedure, in order to achieve HFE configuration. Flux expanded configurations have been experimentally realized both in ohmic and heated plasma with and without nitrogen seeding. Preliminary results on the effects of the flux expansion variation on total power radiation increase will be also briefly discussed. © 2018 Elsevier B.V.

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.

Fast ions have recently been found to significantly impact and partially suppress plasma turbulence both in experimental and numerical studies in a number of scenarios. Understanding the underlying physics and identifying the range of their beneficial effect is an essential task for future fusion reactors, where highly energetic ions are generated through fusion reactions and external heating schemes. However, in many of the gyrokinetic codes fast ions are, for simplicity, treated as equivalent-Maxwellian-distributed particle species, although it is well known that to rigorously model highly non-thermalised particles, a non-Maxwellian background distribution function is needed. To study the impact of this assumption, the gyrokinetic code GENE has recently been extended to support arbitrary background distribution functions which might be either analytic, e.g. slowing down and bi-Maxwellian, or obtained from numerical fast ion models. A particular JET plasma with strong fast-ion related turbulence suppression is revised with these new code capabilities both with linear and nonlinear gyrokinetic simulations. It appears that the fast ion stabilization tends to be less strong but still substantial with more realistic distributions, and this improves the quantitative power balance agreement with experiments., Submitted to Physics of Plasmas

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 ...

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.

This paper further analyses the EDGE2D-EIRENE simulations presented by [Chankin et al 2017 Nucl. Mater. Energy 12 273], of L-mode JET plasmas in vertical-vertical (VV) and vertical-horizontal (VH) divertor configurations. As expected, the simulated outer divertor ionisation source peaks near the separatrix in VV and radially further out in VH. We identify the reflections of recycled neutrals from lower divertor tiles as the primary mechanism by which ionisation is concentrated on the outer divertor separatrix in the VV configuration. These lower tile reflection pathways (of neutrals from the outer divertor, and to an even greater extent from the inner divertor) dominate the outer divertor separatrix ionisation. In contrast, the lower-tile-reflection pathways are much weaker in the VH simulation and its outer divertor ionisation is dominated by neutrals which do not reflect from any surfaces. Interestingly, these differences in neutral pathways give rise to strong differences in the heat flux density width λq at the outer divertor entrance: λq = 3.2 mm in VH compared to λq = 11.8 mm in VV. In VH, a narrow channel exists in the near scrape-off-layer (SOL) where the convected heat flux, driven by strong Er × B flow and thermoelectric current, dominates over the conducted heat flux. The width of this channel sets λq and is determined by the radial distance between the separatrix and the ionisation peak in the outer divertor EURATOM 633053 RCUK Energy Programme P012450/1

In the preparation for the Deuterium-Tritium campaign, JET will operate with a tritium plasma. The T + T reaction consists of two notable channels: (1) T + T -> He-4 + 2n, (2) T + T -> He-5 + ...

In this paper, the dynamic characteristics of type-I ELM time-series from the JET tokamak, the world's largest magnetic confinement plasma physics experiment, have been investigated. The dynamic analysis has been focused on the detection of nonlinear structure in D a radiation time series. Firstly, the method of surrogate data has been applied to evaluate the statistical significance of the null hypothesis of static nonlinear distortion of an underlying Gaussian linear process. Several nonlinear statistics have been evaluated, such us the time delayed mutual information, the correlation dimension and the maximal Lyapunov exponent. The obtained results allow us to reject the null hypothesis, giving evidence of underlying nonlinear dynamics. Moreover, no evidence of low-dimensional chaos has been found; indeed, the analysed time series are better characterized by the power law sensitivity to initial conditions which can suggest a motion at the 'edge of chaos', at the border between chaotic and regular non-chaotic dynamics. This uncertainty makes it necessary to further investigate about the nature of the nonlinear dynamics. For this purpose, a second surrogate test to distinguish chaotic orbits from pseudoperiodic orbits has been applied. In this case, we cannot reject the null hypothesis which means that the ELM time series is possibly pseudo-periodic. In order to reproduce pseudo-periodic dynamical properties, a periodic state-of-the-art model, proposed to reproduce the ELM cycle, has been corrupted by a dynamical noise, obtaining time series qualitatively in agreement with experimental time series.

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.

International audience; The paper is dedicated to the characterisation of multi-strap ICRH antenna coupling to plasma. Relevance of traditional concept of coupling resistance to antennas with mutually coupled straps is revised and the importance of antenna port excitation consistency for application of the concept is highlighted. A method of antenna S-matrix measurement in presence of plasma is discussed allowing deeper insight into the problem of antenna-plasma coupling. The method is based entirely on the RF plant hardware and control facilities available at JET and it involves application of variable phasing between the antenna straps during the RF plant operations at >100 kW. Unlike traditional techniques relying on low-power (~10 mW) network analysers, the applied antenna voltage amplitudes are relevant to practical conditions of ICRH operations; crucially, they are high enough to minimise possible effects of antenna loading non-linearity due to the RF sheath effects and other phenomena which could affect low-power measurements. The method has been successfully applied at JET to conventional 4-port ICRH antennas energised at frequencies of 33 MHz, 42 MHz and 51 MHz during L-mode plasma discharges while different gas injection modules (GIMs) were used to maintain comparable plasma densities during the pulses. The S-matrix assessment and its subsequent processing yielding `global' antenna coupling resistances in conditions of equalised port maximum voltages allowed consistent description of antenna coupling to plasma at different strap phasing, operational frequencies and applied GIMs. Comprehensive experimental characterisation of mutually coupled antenna straps in presence of plasma also provided a unique opportunity for in-depth verification of TOPICA computer simulations.

For the first time, over five confinement times, the self-consistent flux driven time evolution of heat, momentum transport and particle fluxes of electrons and multiple ions including Tungsten (W) is modeled within the integrated modeling platform JETTO (Romanelli et al 2014 Plasma Fusion Res. 9 1–4), using first principle-based codes: namely, QuaLiKiz (Bourdelle et al 2016 Plasma Phys. Control. Fusion 58 014036) for turbulent transport and NEO (Belli and Candy 2008 Plasma Phys. Control. Fusion 50 95010) for neoclassical transport. For a JET-ILW pulse, the evolution of measured temperatures, rotation and density profiles are successfully predicted and the observed W central core accumulation is obtained. The poloidal asymmetries of the W density modifying its neoclassical and turbulent transport are accounted for. Actuators of the W core accumulation are studied: removing the central particle source annihilates the central W accumulation whereas the suppression of the torque reduces significantly the W central accumulation. Finally, the presence of W slightly reduces main ion heat turbulent transport through complex nonlinear interplays involving radiation, effective charge impact on ITG and collisionality. EURATOM 633053

Reliable algorithms for disruption avoidance and prediction are foreseen to play a fundamental role in the JET control system for the successful operation of the machine in the upcoming deuterium-Tritium campaigns. The integration of such algorithms is expected to be a key part also in the implementation of the ITER plasma control system. So far, most of the effort has been devoted to the prediction of disruptions, which is required to mitigate the effects of these transient events, protecting the integrity of in-vessel components. Nevertheless, in order to put in place recover strategies or to have the possibility of a soft landing for the plasma current, the paradigm must be shifted to avoiding disruptions. In this paper, plasma profile-based indicators will be statistically analyzed showing their potential in such a perspective, where warning times and reliability of detection are crucial.

Experimental observations in JET tokamak plasmas and gyrokinetic simulations point to an important role, for electron heat transport, of electron-scale instabilities and of their interaction with ion-scale instabilities. Since these effects are maximized for strong electron heating and ion-scale modes close to marginal stability, these findings are of high relevance for ITER plasmas, featuring both conditions. Gyrokinetic and quasi-linear transport models accounting for multi-scale effects are assessed against JET experimental results.

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.

The connection between plasma physics and neutronics is crucial for the understanding of the operation and performance of modern and future tokamak devices. Neutrons are one of the primary carriers ...

The design and operation of future fusion devices relying on H-mode plasmas requires reliable modelling of edge-localized modes (ELMs) for precise prediction of divertor target conditions. An exten ...

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 present paper describes the findings of the calibration of the neutron yield monitors on the Joint European Torus (JET) performed in 2013 using a 252Cf source deployed inside the torus by the remote handling system, with particular regard to the calibration of fission chambers which provide the time resolved neutron yield from JET plasmas. The experimental data obtained in toroidal, radial, and vertical scans are presented. These data are first analysed following an analytical approach adopted in the previous neutron calibrations at JET. In this way, a calibration function for the volumetric plasma source is derived which allows us to understand the importance of the different plasma regions and of different spatial profiles of neutron emissivity on fission chamber response. Neutronics analyses have also been performed to calculate the correction factors needed to derive the plasma calibration factors taking into account the different energy spectrum and angular emission distribution of the calibrating (point) 252Cf source, the discrete positions compared to the plasma volumetric source, and the calibration circumstances. All correction factors are presented and discussed. We discuss also the lessons learnt which are the basis for the on-going 14 MeV neutron calibration at JET and for ITER. © 2017 EURATOM.

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.

The divertor solutions obtained with the plasma edge modelling tool SOLPS5.0 are discussed. The code results are benchmarked against carefully analysed L-mode discharges at various density levels with and without impurity seeding in the full-metal tokamaks ASDEX Upgrade and JET. The role of the cross-field drifts and currents in the solutions is analysed in detail, and the improvements achieved by fully activating the drift and current terms in view of matching the experimental signals are addressed. The persisting discrepancies are also discussed.

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.

Collection and ex situ studies of dust generated in controlled fusion devices during plasma operation are regularly carried out after experimental campaigns. Herewith results of the dust survey per ...

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 ...

Results are presented for the dust survey performed at JET after the second experimental campaign with the ITER-Like Wall: 2013-2014. Samples were collected on adhesive stickers from several different positions in the divertor both on the tiles and on the divertor carrier. Brittle dust-forming deposits on test mirrors from the inner divertor wall were also studied. Comprehensive characterization accomplished by a wide range of high-resolution microscopy techniques, including focused ion beam, has led to the identification of several classes of particles: (i) beryllium flakes originating either from the Be coatings from the inner wall cladding or Be-rich mixed co-deposits resulting from material migration; (ii) beryllium droplets and splashes; (iii) tungsten and nickel-rich (from Inconel) droplets; (iv) mixed material layers with a various content of small (8-200 nm) W-Mo and Ni-based debris. A significant content of nitrogen from plasma edge cooling has been identified in all types of co-deposits. A comparison between particles collected after the first and second experimental campaign is also presented and discussed. (C) 2016 Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2016.11.027

A strong effect of divertor configuration on the threshold power for the L-H transition (P-LH) was observed in recent JET experiments in the new ITER-like Wall (ILW) [1-3]. Following a series of EDGE2D-EIRENE code simulations with Be impurity and drifts a possible mechanism for the P-LH variation with the divertor geometry is proposed. Both experiment and code simulations show that in the configuration with lower neutral recycling near the outer strike point (OSP), electron temperature (T-e) peaks near the OSP prior to the L-H transition, while in the configuration with higher OSP recycling T-e peaks further out in the scrape-offlayer (SOL) and the plasma stays in the L-mode at the same input power. Code results show large positive radial electric field (E-r) in the near SOL under lower recycling conditions leading to a large E x B shear across the separatrix which may trigger earlier (at lower input power) edge turbulence suppression and lower P-LH. Suppressed T-e's at OSP in configurations with strike points on vertical targets (VT) were observed earlier and explained by a geometrical effect of neutral recycling near this particular position, whereas in configurations with strike points on horizontal targets (HT) the OSP appears to be more open for neutrals (see e.g. review paper [4]). (C) 2016 Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2016.10.004

A series of neon seeded JET ELMy H-mode pulses is considered from the modeling as well as from the experimental point of view. For two different Ne seeding rates and two different D puffing gas levels the heating power, P-heat, is in the range 22-29.5 MW. The main focus is on the numerical reconstruction of the total radiated power (which mostly depends on the W concentration) and its distribution between core and divertor and of Z(eff)(which mostly depends on the Ne concentration). To model self-consistently the core and the SOL two input parameters had to be adjusted case by case: the SOL diffusivity, D SOL, and the core impurity inward pinch, v(pinch). D-SOL had to be increased with increasing Gamma(Ne) and the level of v(pinch) had to be changed, for any given Gamma(Ne), according to the level of P-heat : it decreases with increasing P-heat. Since the ELM frequency, f(ELM), is experimentally correlated with P-heat, (it increases with P-heat) the impurity inward pinch can be seen as to depend on f(ELM). Therefore, to maintain a low v(pinch) level (i.e. high f(ELM)) Gamma(Ne)/P-heat should not exceed a certain threshold, which slightly increases with the Gamma(D) puffing rate. This might lead to a limitation in the viability of reducing the target heat load by Ne seeding at moderate Gamma(D), while keeping Z(eff) at acceptably low level. (C) 2016 Published by Elsevier Ltd. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2016.10.012

A set of quartz crystal microbalances (QMB) was used at JET with full carbon wall to monitor mass erosion/ deposition rates in the remote areas of the divertor. After introduction of the ITER-like wall (ILW) in JET with beryllium main wall and tungsten divertor, strong reduction of the material deposition and accompanied fuel retention was observed. Therefore the existing QMB electronics have been modified to improve the accuracy of frequency measurements by a factor of ten down to 0.1 Hz which corresponds to 1.4 ng cm(-2). The averaged deposition rates of 1.2-3 ng cm(-2) s(-1) and erosion rates of 5.6-8.1 ng cm(-2) s(-1) were observed in the inner divertor of JET-ILW with the inner strike point positions close to the bottom edge of vertical tile 3 and at the horizontal tile 4 respectively. The erosion with averaged rates of approximate to 2.1 ng cm(-2) s(-1) and approximate to 120 ng cm(-2) s(-1) were observed in the outer divertor for the outer strike point positions at tile 5 and tile 6 respectively. (C) 2017 The Authors. Published by Elsevier Ltd. This is an open access article under the CC BY-NC-ND license. For complete list of authors see http://dx.doi.org/10.1016/j.nme.2017.03.007