The here presented work studies the dynamics of filaments using 3D fluid simulations in the presence of detached background profiles. It was found that evolving the neutrals on the time-scale of the filament did not have a significant impact on the dynamics of the filament. In general a decreasing filament velocity with increasing plasma background density has been observed, with the exception of detachment onset, where a temporarily increase in radial velocity occurs. The decreasing trend with temporary increase was found for filaments around the critical size and larger, while smaller filaments where less affected by detachment. With detachment the critical filament size increased, as larger filaments were faster in detached conditions. This breaks the trend of attached conditions, where the critical size decreases with increasing density.

This thesis investigates the influence of background profiles and neutrals on scrape-off layer (SOL) filaments in fusion devices using 3D simulations. This was achieved by seeding filaments on steady-state background profiles including neutrals. In terms of the filament motion, both direct and indirect interaction between the neutrals and filaments have been studied. The simulations have been carried out using the STORM module, which is implemented using the BOUT++ library. STORM has been extended to include neutrals, using a new neutrals C++-API which also uses BOUT++. Filaments have been seeded on top of plasma backgrounds with changing plasma density and temperatures. The profiles have been generated by running a reduced set of equations to steady state. The backgrounds cover different regimes, namely low-recycling regime, high recycling as well as detached plasma conditions. In the attached conditions, the filament’s radial velocity increased with increasing target temperature, suggesting that sheath currents and plasma viscosity significantly influence the studied filaments. Interaction with the neutrals reduced the filament’s radial velocity, especially in the high tem- perature and high density conditions. As detachment was reached, the filament’s velocity increased temporar- ily, as sheath currents were strongly suppressed by the high resistivity, caused by a highly dense, cold plasma in front of the target, before the general decreasing trend is recovered for higher densities. In detached con- ditions the neutrals had an increased impact on the filament dynamics, and partially replaced sheath currents by friction in the detachment front.

This paper presents the effect of self-consistent plasma backgrounds including plasma-neutral interactions, on the dynamics of filament propagation. The principle focus is on the influence of the neutrals on the filament through both direct interactions and through their influence on the plasma background. Both direct and indirect interactions influence the motion of filaments. A monotonic increase of filament peak velocity with upstream electron temperature is observed, while a decrease with increasing electron density is observed. If ordered by the target temperature, the density dependence disappears and the filament velocity is only a function of the target temperature. Smaller filaments keep a density dependence, as a result of the density dependence of the plasma viscosity. The critical size δ^*, where filaments are fastest, is shifted to larger sizes for higher densities, due to the plasma viscosity. If the density dependence of the plasma viscosity is removed, δ^* has no temperature dependence, but rather a density dependence.

The present paper offers an overview of the potential of ion cyclotron resonance heating (ICRH) or radio frequency heating for the DEMO machine. It is found that various suitable heating schemes are available. Similar to ITER and in view of the limited bandwidth of about 10 MHz that can be achieved to ensure optimal functioning of the launcher, it is proposed to make core second harmonic tritium heating the key ion heating scheme, assisted by fundamental cyclotron heating He-3 in the early phase of the discharge; for the present design of DEMO-with a static magnetic field strength of B-o = 5.855 T-that places the T and 3He layers in the core for f = 60 MHz and suggests centering the bandwidth around that main operating frequency. In line with earlier studies for hot, dense plasmas in large-size magnetic confinement machines, it is shown that good single pass absorption is achieved but that the size as well as the operating density and temperature of the machine cause the electrons to absorb a non-negligible fraction of the power away from the core when core ion heating is aimed at. Current drive and alternative heating options are briefly discussed and a dedicated computation is done for the traveling wave antenna, proposed for DEMO in view of its compatibility with substantial antenna-plasma distances. The various tasks that ICRH can fulfill are briefly listed. Finally, the impact of transport and the sensitivity of the obtained results to changes in the machine parameters is commented on.

This paper presents results of JET ITER-like wall L-mode experiments in hydrogen and deuterium (D) plasmas, dedicated to the study of the isotope dependence of ion heat transport by determination of the ion critical gradient and stiffness by varying the ion cyclotron resonance heating power deposition. When no strong role of fast ions in the plasma core is expected, the main difference between the two isotope plasmas is determined by the plasma edge and the core behavior is consistent with a gyro-Bohm scaling. When the heating power (and the fast ion pressure) is increased, in addition to the difference in the edge region, also the plasma core shows substantial changes. The stabilization of ion heat transport by fast ions, clearly visible in D plasmas, appears to be weaker in H plasmas, resulting in a higher ion heat flux in H with apparent anti-gyro-Bohm mass scaling. The difference is found to be caused by the different fast ion pressure between H and D plasmas, related to the heating power settings and to the different fast ion slowing down time, and is completely accounted for in non-linear gyrokinetic simulations. The application of the TGLF quasi-linear model to this set of data is also discussed.

The need for predictive capabilities greater than 95% with very limited false alarms are demanding requirements for reliable disruption prediction systems in tokamaks such as JET or, in the near future, ITER.The prediction of an upcoming disruption must be provided sufficiently in advance in order to apply effective disruption avoidance or mitigation actions to prevent the machine from being damaged.In this paper, following the typical machine learning workflow, a generative topographic mapping (GTM) of the operational space of JET has been built using a set of disrupted and regularly terminated discharges.In order to build the predictive model, a suitable set of dimensionless, machine-independent, physics-based features have been synthesized, which make use of 1D plasma profile information, rather than simple zero-D time series.The use of such predicting features, together with the power of the GTM in fitting the model to the data, obtains, in an unsupervised way, a 2D map of the multi-dimensional parameter space of JET, where it is possible to identify a boundary separating the region free from disruption from the disruption region.In addition to helping in operational boundaries studies, the GTM map can also be used for disruption prediction exploiting the potential of the developed GTM toolbox to monitor the discharge dynamics.Following the trajectory of a discharge on the map throughout the different regions, an alarm is triggered depending on the disruption risk of these regions.The proposed approach to predict disruptions has been evaluated on a training and an independent test set and achieves very good performance with only one tardive detection and a limited number of false detections.The warning times are suitable for avoidance purposes and, more important, the detections are consistent with physical causes and mechanisms that destabilize the plasma leading to disruptions.

The First Mirror Test for ITER has been carried out in JET with mirrors exposed during: (i) the third ILW campaign (ILW-3, 2015-2016, 23.33 h plasma) and (ii) all three campaigns, i.e. ILW-1 to ILW-3: 2011-2016, 63,52 h in total. All mirrors from main chamber wall show no significant changes of the total reflectivity from the initial value and the diffuse reflectivity does not exceed 3% in the spectral range above 500 nm. The modified layer on surface has very small amount of impurities such as D, Be, C, N, O and Ni. All mirrors from the divertor (inner, outer, base under the bulk W tile) lost reflectivity by 20-80% due to the beryllium-rich deposition also containing D, C, N, O, Ni and W. In the inner divertor N reaches 5 x 10(17) cm(-2), W is up to 4.3 x 10(17) cm(-2), while the content of Ni is the greatest in the outer divertor: 3.8 x 10(17) cm(-2). Oxygen-18 used as the tracer in experiments at the end of ILW-3 has been detected at the level of 1.1 x 10(16) cm(-2). The thickness of deposited layer is in the range of 90 nm to 900 nm. The layer growth rate in the base (2.7 pm s(-1)) and inner divertor is proportional to the exposure time when a single campaign and all three are compared. In a few cases, on mirrors located at the cassette mouth, flaking of deposits and erosion occurred.

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

This paper outlines an approach towards improved rigour in tokamak turbulence transport model validation within integrated modelling. Gaussian process regression (GPR) techniques were applied for profile fitting during the preparation of integrated modelling simulations allowing for rigourous sensitivity tests of prescribed initial and boundary conditions as both fit and derivative uncertainties are provided. This was demonstrated by a JETTO integrated modelling simulation of the JET ITER-like-wall H-mode baseline discharge #92436 with the QuaLiKiz quasilinear turbulent transport model, which is the subject of extrapolation towards a deuterium-tritium plasma. The simulation simultaneously evaluates the time evolution of heat, particle, and momentum fluxes over similar to 10 confinement times, with a simulation boundary condition at rho(tor) = 0.85. Routine inclusion of momentum transport prediction in multi-channel flux-driven transport modelling is not standard and is facilitated here by recent developments within the QuaLiKiz model. Excellent agreement was achieved between the fitted and simulated profiles for n(e), T-e, T-i, and Omega(tor) within 2 sigma, but the simulation underpredicts the mid-radius Ti and overpredicts the core n(e) and T-e profiles for this discharge. Despite this, it was shown that this approach is capable of deriving reasonable inputs, including derivative quantities, to tokamak models from experimental data. Furthermore, multiple figures-of-merit were defined to quantitatively assess the agreement of integrated modelling predictions to experimental data within the GPR profile fitting framework.

ERO is a 3D Monte-Carlo impurity transport and plasma-surface interaction code. In 2011 it was applied for the ITER first wall (FW) life time predictions [1] (critical blanket module BM11). After that the same code was significantly improved during its application to existing fusion-relevant plasma devices: the tokamak JET equipped with an ITER-like wall and linear plasma device PISCES-B. This has allowed testing the sputtering data for beryllium (Be) and showing that the "ERO-min" fit based on the large (50%) deuterium (D) surface content is well suitable for plasma-wetted areas (D plasma). The improved procedure for calculating of the effective sputtering yields for each location along the plasma-facing surface using the recently developed semi-analytical sheath approach was validated. The re-evaluation of the effective yields for BM11 following the similar revisit of the JET data has indicated significant increase of erosion and motivated the current re-visit of ERO simulations.

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.

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.

This paper presents the first application to real JET data of the new equilibrium code NICE which enables the consistent resolution of the inverse equilibrium reconstruction problem in the framework of non-linear free-boundary equilibrium coupled to the Stokes model equation for polarimetry. The conducted numerical experiments enable first of all to validate NICE by comparing it to the well-established EFIT code on 4 selected high performance shots. Secondly the results indicate that the fit to polarimetry measurements clearly benefits from the use of Stokes vector measurements compared to the classical case of Faraday measurements, and that the reconstructed p' and ff' profiles are better constrained with smaller error bars and are closer to the profiles reconstructed by EFTM, the EFIT JET code using internal MSE constraints.

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.

Progress in understanding lower hybrid current drive (LHCD) at high density has been made through experiments and modeling, which is encouraging given the need for an efficient off-axis current profile control technique in burning plasma. By reducing the wall recycling of neutrals, the edge temperature is increased and the effect of parametric instability (PI) and collisional absorption (CA) is reduced, which is beneficial for increasing the current drive efficiency. Strong single pass absorption is preferred to prevent CA and high LH operating frequency is essential for wave propagation to the core region at high density, presumably to mitigate the effect of PI. The dimensionless parameter that characterizes LH wave accessibility and wave refraction for the experiments in this joint study is shown to bracket the region in parameter space where ITER LHCD experiments will operate in the steady state scenario phase. Further joint experiments and cross modeling are necessary to understand the LHCD physics in weak damping regimes which would increase confidence in predictions for ITER where the absorption is expected to be strong.

Several important ELM control techniques are in large part motivated by the empirically observed inverse relationship between average ELM energy loss and ELM frequency in a plasma. However, to ensure a reliable effect on the energy released by the ELMs, it is important that this relation is verified for individual ELM events. Therefore, in this work the relation between ELM energy loss (WELM) and waiting time (ΔtELM) is investigated for individual ELMs in a set of ITER-like wall plasmas in JET. A comparison is made with the results from a set of carbon-wall and nitrogen-seeded ITER-like wall JET plasmas. It is found that the correlation between W ELM and ΔtELM for individual ELMs varies from strongly positive to zero. Furthermore, the effect of the extended collapse phase often accompanying ELMs from unseeded JET ILW plasmas and referred to as the slow transport event (STE) is studied on the distribution of ELM durations, and on the correlation between W ELM and ΔtELM . A high correlation between W ELM and ΔtELM , comparable to CW plasmas is only found in nitrogen-seeded ILW plasmas. Finally, a regression analysis is performed using plasma engineering parameters as predictors for determining the region of the plasma operational space with a high correlation between W ELM and ΔtELM .

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.

This paper presents the effect of self-consistent plasma backgrounds, on the dynamics of filament propagation. The critical size δ ∗, is an important scaling parameter for filaments. It is here defined as the perpendicular size, where filaments are fastest and is shifted to larger sizes for higher densities, due to the plasma viscosity. If the density dependency of the plasma viscosity is included, δ ∗ does not show a temperature dependency, but instead a density dependency is observed.