Search

Nitrogen seeded detachment has been achieved in the tokamak a configuration variable (TCV) in alternative divertor configurations (ADCs), namely X-divertor and X-point target, with and without baffles in H-mode plasmas with high core confinement. Both ADCs show a remarkable reduction in the inter-ELM particle and heat fluxes to the target compared to the standard divertor configuration. 95%–98% of the inter-ELM peak heat flux to the target is mitigated as a synergetic effect of ADCs, baffling, and nitrogen seeded detachment. The effect of divertor geometry and baffles on core-divertor compatibility is investigated in detail. The power balance in these experiments is also investigated to explore the physics behind the observed reduction in heat fluxes in the ADCs.

The Tokamak a Configuration Variable (TCV) has recently been equipped with gas baffles to increase its divertor closure for a broad range of divertor magnetic geometries. First experimental results reported in Reimerdes et al. (2021) demonstrated compatibility with a broad range of divertor magnetic geometries and confirmed the main design constraints of the baffles, in particular an increased divertor neutral pressure. The present article presents a more in-depths analysis and extended experiments of this first baffle assessment on the TCV boundary plasma. It is shown that the divertor neutral pressure increased following the installation of the baffles, as predicted by SOLPS-ITER and SolEdge2D-EIRENE simulations. Varying the divertor closure by changing the position of the plasma showed that the plasma equilibrium designed to assess the baffle effect was not far from the optimal trade-off between plasma plugging and recycling on the baffles. The baffles facilitate access to a partially detached regime in both L- and H-modes. In L-Mode, with the ion ∇ B -drift directed from the X-Point to the plasma core, a reduction of the line-averaged density detachment threshold by approximately 20% is observed at the outer target, while inner strike point detachment is only achieved in the presence of baffles. Multispectral imaging shows that the CIII front moves from the outer target towards the X-Point at a lower line-averaged plasma density, indicating a colder outer leg. In H-mode, the CIII front is generally located near the X-Point between the ELMs, while without baffles, N 2 -seeding is required to move the front up to that location.

Neoclassical and turbulent heavy impurity transport in tokamak core plasmas are determined by main ion temperature, density and toroidal rotation profiles. Thus, in order to understand and prevent experimental behaviour of W accumulation, flux-driven integrated modelling of main ion heat and particle transport over multiple confinement times is a vital prerequisite. For the first time, the quasilinear gyrokinetic code QuaLiKiz has been applied for successful predictions of core kinetic profiles in an ASDEX Upgrade H-mode discharge in the turbulence dominated region within the integrated modelling suite JETTO. Neoclassical contributions are calculated by NCLASS; auxiliary heat and particle deposition profiles due to NBI and ECRH prescribed from previous analysis with TRANSP. Turbulent and neoclassical contributions are insufficient in explaining main ion heat and particle transport inside the $q=1$ surface, necessitating the prescription of further transport coefficients to mimic the impact of MHD activity on central transport. The ion to electron temperature ratio at the simulation boundary at $\rho_\mathrm{tor} = 0.85$ stabilizes ion scale modes while destabilizing ETG modes when significantly exceeding unity. Careful analysis of experimental measurements using Gaussian process regression techniques is carried out to explore reasonable uncertainties. In following trace W impurity transport simulations performed with additionally NEO, neoclassical transport under consideration of poloidal asymmetries alone is found to be insufficient to establish hollow central W density profiles. Reproduction of these conditions measured experimentally is found possible only when assuming the direct impact of a saturated $(m,n)=(1,1)$ MHD mode on heavy impurity transport., Comment: 26 pages, 17 figures

A detailed cross-device investigation on the role of filamentary dynamics in high-density regimes has been performed within the EUROfusion framework, comparing the ASDEX Upgrade (AUG) and TCV tokamaks. Both devices run density ramp experiments at different levels of plasma current, keeping the toroidal field or q 95 constant in order to disentangle the role of the parallel connection length and the current. During the scan at a constant toroidal field, in both devices the scrape-off layer (SOL) profiles tend to develop a clear SOL density shoulder at a lower edge density whenever the current is reduced. Different current behaviour is substantially reconciled in terms of the edge density normalized to the Greenwald fraction. During the scan at constant q 95 AUG exhibits similar behaviour, whereas in TCV no upstream profile modification signature has been observed at lower current levels. The latter behaviour has been ascribed to the lack of target density rollover. The relation between the upstream density profile modification and detachment condition has been investigated. For both devices the relation between blob size and the SOL density e-folding length is found independent of the plasma current, with the observation of a clear increase in blob size and the edge density normalized to a Greenwald fraction. ASDEX Upgrade has also explored filamentary behaviour in the H-mode. The experiments in AUG have focused on the role of neutrals, performing discharges with and without cryogenic pumps, highlighting how high neutral pressure, not only in the divertor but also at the midplane, is needed in order to develop an H-mode SOL profile shoulder in AUG.

Synchrotron radiation observed in a quiescent TCV runaway discharge is studied using filtered camera images targeting three distinct wavelength intervals. Through the tomographic SART procedure the high momentum, high pitch angle part of the spatial and momentum distribution of these relativistic particles is reconstructed. Experimental estimates of the distribution are important for verification and refinement of formation-, decay- and transport-models underlying runaway avoidance and mitigation strategy design. Using a test distribution it is demonstrated that the inversion procedure provides estimates accurate to within a few tens of percent in the region of phase-space contributing most to the synchrotron image. We find that combining images filtered around different parts of the emission spectrum widens the probed part of momentum-space and reduces reconstruction errors. Next, the SART algorithm is used to obtain information on the spatiotemporal runaway momentum distribution in a selected TCV discharge. The momentum distribution is found to relax towards an avalanche-like exponentially decaying profile. Anomalously high pitch angles and a radial profile increasing towards the edge are found for the most strongly emitting particles in the distribution. Pitch angle scattering by toroidal magnetic field ripple is consistent with this picture. An alternative explanation is the presence of high frequency instabilities in combination with the formation of a runaway shell at the edge of the plasma., 18 pages, 13 figures, submitted to Nuclear Fusion

Recent experiments on Tokamak à Configuration Variable have made significant progress toward partial detachment of the outer divertor in neutral beam heated H-mode plasmas in conventional and alternative divertor configurations. The heating power required to enter H-mode was measured in a range of divertor configurations. It is found that at the core densities most favourable for H-mode access, the L–H threshold power is largely independent of the poloidal flux expansion and major radius of the outer divertor, and in the snowflake minus configuration. A factor 2 reduction in the outer divertor power load was achieved in ELM-free (using a fuelling and nitrogen seeding) and ELMy (using nitrogen seeding) H-mode plasmas. No significant reduction in the outer divertor particle flux was achieved in the ELM-free scenarios, compared with ~30% reduction in the most strongly detached ELMy cases. The poloidal flux expansion at the outer divertor was not found to significantly alter the cooling of the divertor in the ELM-free scenarios.

Edge localized modes (ELMs) are repetitive instabilities driven by the large pressure gradients and current densities in the edge of H-mode plasmas. Type-I ELMs lead to a fast collapse of the H-mode pedestal within several hundred microseconds to a few milliseconds. Localized transient heat fluxes to divertor targets are expected to exceed tolerable limits for ITER, requiring advanced insights into ELM physics and applicable mitigation methods. This paper describes how non-linear magneto-hydrodynamic (MHD) simulations can contribute to this effort. The JOREK code is introduced, which allows the study of large-scale plasma instabilities in tokamak X-point plasmas covering the main plasma, the scrape-off layer, and the divertor region with its finite element grid. We review key physics relevant for type-I ELMs and show to what extent JOREK simulations agree with experiments and help reveal the underlying mechanisms. Simulations and experimental findings are compared in many respects for type-I ELMs in ASDEX Upgrade. The role of plasma flows and non-linear mode coupling for the spatial and temporal structure of ELMs is emphasized, and the loss mechanisms are discussed. An overview of recent ELM-related research using JOREK is given, including ELM crashes, ELM-free regimes, ELM pacing by pellets and magnetic kicks, and mitigation or suppression by resonant magnetic perturbation coils (RMPs). Simulations of ELMs and ELM control methods agree in many respects with experimental observations from various tokamak experiments. On this basis, predictive simulations become more and more feasible. A brief outlook is given, showing the main priorities for further research in the field of ELM physics and further developments necessary.

Control of the time between sawtooth crashes, necessary for ITER and DEMO, requires real-time detection of the moment of the sawtooth crash. In this paper, estimation of sawtooth crash times is demonstrated using the model-based interacting multiple model (IMM) estimator, based on simplified models for the sawtooth crash. In contrast to previous detectors, this detector uses the spatial extent of the sawtooth crash as detection characteristic. The IMM estimator is tuned and applied to multiple ECE channels at once. A model for the sawtooth crash is introduced, which is used in the IMM algorithm. The IMM algorithm is applied to seven datasets from the ASDEX Upgrade tokamak. Five crash models with different mixing radii are used. All sawtooth crashes that have been identified beforehand by visual inspection of the data, are detected by the algorithm. A few additional detections are made, which upon closer inspection are seen to be sawtooth crashes, which show a partial reconnection. A closer inspection of the detected normal crashes shows that about 42% are not well fitted by any of the full reconnection models and show some characteristics of a partial reconnection. In some case, the measurement time is during the sawtooth crashes, which also results in an incorrect estimate of the mixing radius. For data provided at a sampling rate of 1 kHz, the run time of the IMM estimator is below 1 ms, thereby fulfilling real-time requirements.

In the non-linear phase of a sawtooth, the complete reconnection of field lines around the q = 1 flux surface often occurs resulting in a radial displacement of the plasma core. A complete time-dependent electromagnetic model of this type of reconnection has been developed and implemented in the EBdyna_go code. This contribution aims at studying the behaviour of ions, both impurity and fast particles, in the pattern of reconnecting field lines during sawtoothing plasma experiments in the ASDEX Upgrade tokamak by using the newly developed numerical framework. Simulations of full reconnection with tungsten impurity that include the centrifugal force are achieved and recover the soft x-ray measurements. Based on this full-reconnection description of the sawtooth, a simple tool dedicated to estimate the duration of the reconnection is introduced. This work then studies the redistribution of fast ions during several experimentally observed sawteeth. In some cases of sawteeth at ASDEX Upgrade, full reconnection is not always observed or expected so the code gives an upper estimate of the actual experimental redistribution. The results of detailed simulations of the crashes are compared with measurements from various diagnostics such as collective Thomson scattering and fast-ion D-alpha (FIDA) spectroscopy, including FIDA tomography. A convincing qualitative agreement is found in different parts of velocity space.