Real-time plasma state monitoring and supervisory control on TCV

In ITER and DEMO, various control objectives related to plasma control must be simultaneously achieved by the plasma control system (PCS), in both normal operation as well as off-normal conditions. The PCS must act on off-normal events and deviations from the target scenario, since certain sequences (chains) of events can precede disruptions. It is important that these decisions are made while maintaining a coherent prioritization between the real-time control tasks to ensure high-performance operation. In this paper, a generic architecture for task-based integrated plasma control is proposed. The architecture is characterized by the separation of state estimation, event detection, decisions and task execution among different algorithms, with standardized signal interfaces. Central to the architecture are a plasma state monitor and supervisory controller. In the plasma state monitor, discrete events in the continuous-valued plasma state are modeled using finite state machines. This provides a high-level representation of the plasma state. The supervisory controller coordinates the execution of multiple plasma control tasks by assigning task priorities, based on the finite states of the plasma and the pulse schedule. These algorithms were implemented on the TCV digital control system and integrated with actuator resource management and existing state estimation algorithms and controllers. The plasma state monitor on TCV can track a multitude of plasma events, related to plasma current, rotating and locked neoclassical tearing modes, and position displacements. In TCV experiments on simultaneous control of plasma pressure, safety factor profile and NTMs using electron cyclotron heating (ECH) and current drive (ECCD), the supervisory controller assigns priorities to the relevant control tasks. The tasks are then executed by feedback controllers and actuator allocation management. This work forms a significant step forward in the ongoing integration of control capabilities in experiments on TCV, in support of tokamak reactor operation

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

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

Testing a prediction model for the H-mode density pedestal against JET-ILW pedestals

The neutral ionisation model proposed by Groebner et al (2002 Phys. Plasmas 9 2134) to determine the plasma density profile in the H-mode pedestal, is extended to include charge exchange processes in the pedestal stimulated by the ideas of Mahdavi et al (2003 Phys. Plasmas 10 3984). The model is then tested against JET H-mode pedestal data, both in a 'standalone' version using experimental temperature profiles and also by incorporating it in the Europed version of EPED. The model is able to predict the density pedestal over a wide range of conditions with good accuracy. It is also able to predict the experimentally observed isotope effect on the density pedestal that eludes simpler neutral ionization models

Overview of JET results for optimising ITER operation

The JET 2019–2020 scientific and technological programme exploited the results of years of concerted scientific and engineering work, including the ITER-like wall (ILW: Be wall and W divertor) installed in 2010, improved diagnostic capabilities now fully available, a major neutral beam injection upgrade providing record power in 2019–2020, and tested the technical and procedural preparation for safe operation with tritium. Research along three complementary axes yielded a wealth of new results. Firstly, the JET plasma programme delivered scenarios suitable for high fusion power and alpha particle (α) physics in the coming D–T campaign (DTE2), with record sustained neutron rates, as well as plasmas for clarifying the impact of isotope mass on plasma core, edge and plasma-wall interactions, and for ITER pre-fusion power operation. The efficacy of the newly installed shattered pellet injector for mitigating disruption forces and runaway electrons was demonstrated. Secondly, research on the consequences of long-term exposure to JET-ILW plasma was completed, with emphasis on wall damage and fuel retention, and with analyses of wall materials and dust particles that will help validate assumptions and codes for design and operation of ITER and DEMO. Thirdly, the nuclear technology programme aiming to deliver maximum technological return from operations in D, T and D–T benefited from the highest D–D neutron yield in years, securing results for validating radiation transport and activation codes, and nuclear data for ITER

Interpretative and predictive modelling of Joint European Torus collisionality scans

Transport modelling of Joint European Torus (JET) dimensionless collisionality scaling experiments in various operational scenarios is presented. Interpretative simulations at a fixed radial position are combined with predictive JETTO simulations of temperatures and densities, using the TGLF transport model. The model includes electromagnetic effects and collisions as well as □(→┬E ) X □(→┬B ) shear in Miller geometry. Focus is on particle transport and the role of the neutral beam injection (NBI) particle source for the density peaking. The experimental 3-point collisionality scans include L-mode, and H-mode (D and H and higher beta D plasma) plasmas in a total of 12 discharges. Experimental results presented in (Tala et al 2017 44th EPS Conf.) indicate that for the H-mode scans, the NBI particle source plays an important role for the density peaking, whereas for the L-mode scan, the influence of the particle source is small. In general, both the interpretative and predictive transport simulations support the experimental conclusions on the role of the NBI particle source for the 12 JET discharges

Impact of fast ions on density peaking in JET: fluid and gyrokinetic modeling

The effect of fast ions on turbulent particle transport, driven by ion temperature gradient (ITG)/ trapped electron mode turbulence, is studied. Two neutral beam injection (NBI) heated JET discharges in different regimes are analyzed at the radial position ρ$_{t}$=0.6, one of them an L-mode and the other one an H-mode discharge. Results obtained from the computationally efficient fluid model EDWM and the gyro-fluid model TGLF are compared to linear and nonlinear gyrokinetic GENE simulations as well as the experimentally obtained density peaking. In these models, the fast ions are treated as a dynamic species with a Maxwellian background distribution. The dependence of the zero particle flux density gradient (peaking factor) on fast ion density, temperature and corresponding gradients, is investigated. The simulations show that the inclusion of a fast ion species has a stabilizing influence on the ITG mode and reduces the peaking of the main ion and electron density profiles in the absence of sources. The models mostly reproduce the experimentally obtained density peaking for the L-mode discharge whereas the H-mode density peaking is significantly underpredicted, indicating the importance of the NBI particle source for the H-mode density profile

Current Research into Applications of Tomography for Fusion Diagnostics

Retrieving spatial distribution of plasma emissivity from line integrated measurements on tokamaks presents a challenging task due to ill-posedness of the tomography problem and limited number of the lines of sight. Modern methods of plasma tomography therefore implement a-priori information as well as constraints, in particular some form of penalisation of complexity. In this contribution, the current tomography methods under development (Tikhonov regularisation, Bayesian methods and neural networks) are briefly explained taking into account their potential for integration into the fusion reactor diagnostics. In particular, current development of the Minimum Fisher Regularisation method is exemplified with respect to real-time reconstruction capability, combination with spectral unfolding and other prospective tasks

The effect of beryllium oxide on retention in JET ITER-like wall tiles

Preliminary results investigating the microstructure, bonding and effect of beryllium oxide formation on retention in the JET ITER-like wall beryllium tiles, are presented. The tiles have been investigated by several techniques: Scanning Electron Microscopy (SEM) equipped with Energy Dispersive X-ray (EDX), Transmission Electron microscopy (TEM) equipped with EDX and Electron Energy Loss Spectroscopy (EELS), Raman Spectroscopy and Thermal Desorption Spectroscopy (TDS). This paper focuses on results from melted materials of the dump plate tiles in JET. From our results and the literature, it is concluded, beryllium can form micron deep oxide islands contrary to the nanometric oxides predicted under vacuum conditions. The deepest oxides analyzed were up to 2-micron thicknesses. The beryllium Deuteroxide (BeOxDy) bond was found with Raman Spectroscopy. Application of EELS confirmed the oxide presence and stoichiometry. Literature suggests these oxides form at temperatures greater than 700 °C where self-diffusion of beryllium ions through the surface oxide layer can occur. Further oxidation is made possible between oxygen plasma impurities and the beryllium ions now present at the wall surface. Under Ultra High Vacuum (UHV) nanometric Beryllium oxide layers are formed and passivate at room temperature. After continual cyclic heating (to the point of melt formation) in the presence of oxygen impurities from the plasma, oxide growth to the levels seen experimentally (approximately two microns) is proposed. This retention mechanism is not considered to contribute dramatically to overall retention in JET, due to low levels of melt formation. However, this mechanism, thought the result of operation environment and melt formation, could be of wider concern to ITER, dependent on wall temperatures

Modelling of tungsten erosion and deposition in the divertor of JET-ILW in comparison to experimental findings

The erosion, transport and deposition of tungsten in the outer divertor of JET-ILW has been studied for an HMode discharge with low frequency ELMs. For this specific case with an inter-ELM electron temperature at the strike point of about 20 eV, tungsten sputtering between ELMs is almost exclusively due to beryllium impurity and self-sputtering. However, during ELMs tungsten sputtering due to deuterium becomes important and even dominates. The amount of simulated local deposition of tungsten relative to the amount of sputtered tungsten in between ELMs is very high and reaches values of 99% for an electron density of 5E13 cm$^{-3}$ at the strike point and electron temperatures between 10 and 30 eV. Smaller deposition values are simulated with reduced electron density. The direction of the B-field significantly influences the local deposition and leads to a reduction if the E×B drift directs towards the scrape-off-layer. Also, the thermal force can reduce the tungsten deposition, however, an ion temperature gradient of about 0.1 eV/mm or larger is needed for a significant effect. The tungsten deposition simulated during ELMs reaches values of about 98% assuming ELM parameters according to free-streaming model. The measured WI emission profiles in between and within ELMs have been reproduced by the simulation. The contribution to the overall net tungsten erosion during ELMs is about 5 times larger than the one in between ELMs for the studied case. However, this is due to the rather low electron temperature in between ELMs, which leads to deuterium impact energies below the sputtering threshold for tungsten

Impact of fast ions on density peaking in JET : fluid and gyrokinetic modeling

The effect of fast ions on turbulent particle transport, driven by ion temperature gradient (ITG)/trapped electron mode turbulence, is studied. Two neutral beam injection (NBI) heated JET discharges in different regimes are analyzed at the radial position rho(t) = 0.6, one of them an L-mode and the other one an H-mode discharge. Results obtained from the computationally efficient fluid model EDWM and the gyro-fluid model TGLF are compared to linear and nonlinear gyrokinetic GENE simulations as well as the experimentally obtained density peaking. In these models, the fast ions are treated as a dynamic species with a Maxwellian background distribution. The dependence of the zero particle flux density gradient (peaking factor) on fast ion density, temperature and corresponding gradients, is investigated. The simulations show that the inclusion of a fast ion species has a stabilizing influence on the ITG mode and reduces the peaking of the main ion and electron density profiles in the absence of sources. The models mostly reproduce the experimentally obtained density peaking for the L-mode discharge whereas the H-mode density peaking is significantly underpredicted, indicating the importance of the NBI particle source for the H-mode density profile