Search

In fusion plasma physics, the large-scale trajectories of energetic particles in magnetic confinement devices are known as orbits. To effectively and efficiently be able to work with orbits, the Orbit Weight Computational Framework (OWCF) was developed. The OWCF constitutes a set of scripts, functions and applications capable of computing, visualizing and working with quantities related to fast-ion (FI) orbits in toroidally symmetric fusion devices. The current version is highly integrated with the DRESS code, which enables the OWCF to compute and analyze the orbit sensitivity for arbitrary neutron- and gammadiagnostics. However, the framework is modular in the sense that any future codes (e.g. FIDASIM) can be easily integrated. The OWCF can also compute projected velocity spectra for FI orbits, which play a key role in many FI diagnostics. Via interactive applications, the OWCF can function both as a tool for investigative research but also for teaching. The OWCF will be used to analyze and simulate the diagnostic results of current and future fusion experiments such as ITER. The orbit weight functions computed with the OWCF can be used to reconstruct the FI distribution in terms of FI orbits from experimental measurements using tomographic inversion.

In fusion plasma physics, the large-scale trajectories of energetic particles in magnetic confinement devices are known as orbits. To effectively and efficiently be able to work with orbits, the Orbit Weight Computational Framework (OWCF) was developed. The OWCF constitutes a set of scripts, functions and applications capable of computing, visualizing and working with quantities related to fast-ion (FI) orbits in toroidally symmetric fusion devices. The current version is highly integrated with the DRESS code, which enables the OWCF to compute and analyze the orbit sensitivity for arbitrary neutron- and gammadiagnostics. However, the framework is modular in the sense that any future codes (e.g. FIDASIM) can be easily integrated. The OWCF can also compute projected velocity spectra for FI orbits, which play a key role in many FI diagnostics. Via interactive applications, the OWCF can function both as a tool for investigative research but also for teaching. The OWCF will be used to analyze and simulate the diagnostic results of current and future fusion experiments such as ITER. The orbit weight functions computed with the OWCF can be used to reconstruct the FI distribution in terms of FI orbits from experimental measurements using tomographic inversion.

In fusion plasma physics, the large-scale trajectories of energetic particles in magnetic confinement devices are known as orbits. To effectively and efficiently be able to work with orbits, the Orbit Weight Computational Framework (OWCF) was developed. The OWCF constitutes a set of scripts, functions and applications capable of computing, visualizing and working with quantities related to fast-ion (FI) orbits in toroidally symmetric fusion devices. The current version is highly integrated with the DRESS code, which enables the OWCF to compute and analyze the orbit sensitivity for arbitrary neutron- and gammadiagnostics. However, the framework is modular in the sense that any future codes (e.g. FIDASIM) can be easily integrated. The OWCF can also compute projected velocity spectra for FI orbits, which play a key role in many FI diagnostics. Via interactive applications, the OWCF can function both as a tool for investigative research but also for teaching. The OWCF will be used to analyze and simulate the diagnostic results of current and future fusion experiments such as ITER. The orbit weight functions computed with the OWCF can be used to reconstruct the FI distribution in terms of FI orbits from experimental measurements using tomographic inversion.

The fast-ion phase-space distribution function in axisymmetric tokamak plasmas is completely described by the three constants of motion: energy, magnetic moment and toroidal canonical angular momentum. In this work, the observable regions of constants-of-motion phase-space, given a diagnostic setup, are identified and explained using projected velocities of the fast ions along the diagnostic lines-of-sight as a proxy for several fast-ion diagnostics, such as fast-ion D alpha spectroscopy, collective Thomson scattering, neutron emission spectroscopy and gamma-ray spectroscopy. The observable region in constants-of-motion space is given by a position condition and a velocity condition, and the diagnostic sensitivity is given by a gyro-orbit and a drift-orbit weighting. As a practical example, 3D orbit weight functions quantifying the diagnostic sensitivity to each point in phase-space are computed and investigated for the future COMPASS-Upgrade and MAST-Upgrade tokamaks.

A new scintillator based fast-ion loss detector (FILD) system has been designed for the Wendelstein 7-X (W7- X) stellarator. The mechanical design of the system is presented here along with engineering analyses of the system. This includes an assessment of the structural loads which also considers electromagnetic forces calculated for possible disruption scenarios as well as a mode frequency analysis of the system. Furthermore, an analysis of the thermal characteristics of the actively cooled probe head under prescribed steady state and transient heat loads is presented. Finally, the optical relay system is described and its properties used to forward model synthetic camera signals using the FILDSIM code and based on theoretical fast-ion losses incident on the probe head calculated using the ASCOT5 code.

Tomographic reconstructions of a 3D fast-ion constants-of-motion phase-space distribution function are computed by inverting synthetic signals based on projected velocities of the fast ions along the diagnostic lines of sight. A spectrum of projected velocities is a key element of the spectrum formation in fast-ion D-alpha spectroscopy, collective Thomson scattering, and gamma-ray and neutron emission spectroscopy, and it can hence serve as a proxy for any of these. The fast-ion distribution functions are parameterised by three constants of motion, the kinetic energy, the magnetic moment and the toroidal canonical angular momentum. The reconstructions are computed using both zeroth-order and first-order Tikhonov regularisation expressed in terms of Bayesian inference to allow uncertainty quantification. In addition to this, a discontinuity appears to be present in the solution across the trapped-passing boundary surface in the three-dimensional phase space due to a singularity in the Jacobian of the transformation from position and velocity space to phase space. A method to allow for this apparent discontinuity while simultaneously penalising large gradients in the solution is demonstrated. Finally, we use our new methods to optimise the diagnostic performance of a set of six fans of sightlines by finding where the detectors contribute most complementary diagnostic information for the future COMPASS-Upgrade tokamak.

Sawtooth oscillations in tokamaks frequently lead to the redistribution of energetic ions, mainly on passing orbits, causing their expulsion from the core. This paper discusses the first measurements of the interaction of fast ions and sawteeth in the Large Helical Device. The crashes were caused by the plasma current induced by Electron Cyclotron Current Drive and Neutral Beam Current Drive. Despite these crashes, there was no detectable redistribution effect on fast ions in either the core or at the edge of the plasma.

Spherical tokamaks are key to the successful design of operating scenarios of future fusion reactors in the areas of divertor physics, neutral beam current drive and fast ion physics. MAST Upgrade, which has successfully concluded its first experimental campaign, was specifically designed to address the role of the radial gradient of the fast ion distribution in driving the excitation of magneto-hydrodynamic (MHD) instabilities, such as toroidal Alfven eigenmodes, fish-bones and long-lived mode, thanks to its two tangential neutral beam injection systems, one on the equatorial plane and one that is vertically shifted 65 cm above the equatorial plane. To study the fast ion dynamics in the presence of such instabilities, as well as of sawteeth and neo-classical tearing modes, several fast ion diagnostics were upgraded and new ones added. Among them, the MAST prototype neutron camera (NC) has been upgraded to six, equatorial sight-lines. The first observations of the confined fast ion behavior with the upgraded NC in a wide range of plasma scenarios characterized by on-axis and/or off-axis heating and different MHD instabilities are presented here. The observations presented in this study confirm previous results on MAST but with a higher level of detail and highlight new physics observations unique to the MAST Upgrade. The results presented here confirm the improved performance of the NC Upgrade, which thus becomes one of the key elements, in combination with the rich set of fast ion diagnostics available on the MAST Upgrade, for a more constrained modeling of the fast ion dynamics in fusion reactor relevant scenarios.

We compute reconstructions of 4D and 5D fast-ion phase-space distribution functions in fusion plasmas from synthetic projections of these functions. The fast-ion phase-space distribution functions originating from neutral beam injection (NBI) at TCV and Wendelstein 7-X (W7-X) at full, half, and one-third injection energies can be distinguished and particle densities of each component inferred based on 20 synthetic spectra of projected velocities at TCV and 680 at W7-X. Further, we demonstrate that an expansion into a basis of slowing-down distribution functions is equivalent to regularization using slowing-down physics as prior information. Using this technique in a Tikhonov formulation, we infer the particle density fractions for each NBI energy for each NBI beam from synthetic measurements, resulting in six unknowns at TCV and 24 unknowns at W7-X. Additionally, we show that installing 40 LOS in each of 17 ports at W7-X, providing full beam coverage and almost full angle coverage, produces the highest quality reconstructions.

Spherical tokamaks are key to the successful design of operating scenarios of future fusion reactors in the areas of divertor physics, neutral beam current drive and fast ion physics. MAST Upgrade, which has successfully concluded its first experimental campaign, was specifically designed to address the role of the radial gradient of the fast ion distribution in driving the excitation of magneto-hydrodynamic (MHD) instabilities, such as toroidal Alfven eigenmodes, fish-bones and long-lived mode, thanks to its two tangential neutral beam injection systems, one on the equatorial plane and one that is vertically shifted 65 cm above the equatorial plane. To study the fast ion dynamics in the presence of such instabilities, as well as of sawteeth and neo-classical tearing modes, several fast ion diagnostics were upgraded and new ones added. Among them, the MAST prototype neutron camera (NC) has been upgraded to six, equatorial sight-lines. The first observations of the confined fast ion behavior with the upgraded NC in a wide range of plasma scenarios characterized by on-axis and/or off-axis heating and different MHD instabilities are presented here. The observations presented in this study confirm previous results on MAST but with a higher level of detail and highlight new physics observations unique to the MAST Upgrade. The results presented here confirm the improved performance of the NC Upgrade, which thus becomes one of the key elements, in combination with the rich set of fast ion diagnostics available on the MAST Upgrade, for a more constrained modeling of the fast ion dynamics in fusion reactor relevant scenarios.

We compute reconstructions of 4D and 5D fast-ion phase-space distribution functions in fusion plasmas from synthetic projections of these functions. The fast-ion phase-space distribution functions originating from neutral beam injection (NBI) at TCV and Wendelstein 7-X (W7-X) at full, half, and one-third injection energies can be distinguished and particle densities of each component inferred based on 20 synthetic spectra of projected velocities at TCV and 680 at W7-X. Further, we demonstrate that an expansion into a basis of slowing-down distribution functions is equivalent to regularization using slowing-down physics as prior information. Using this technique in a Tikhonov formulation, we infer the particle density fractions for each NBI energy for each NBI beam from synthetic measurements, resulting in six unknowns at TCV and 24 unknowns at W7-X. Additionally, we show that installing 40 LOS in each of 17 ports at W7-X, providing full beam coverage and almost full angle coverage, produces the highest quality reconstructions.

We compute reconstructions of 4D and 5D fast-ion phase-space distribution functions in fusion plasmas from synthetic projections of these functions. The fast-ion phase-space distribution functions originating from neutral beam injection (NBI) at TCV and Wendelstein 7-X (W7-X) at full, half, and one-third injection energies can be distinguished and particle densities of each component inferred based on 20 synthetic spectra of projected velocities at TCV and 680 at W7-X. Further, we demonstrate that an expansion into a basis of slowing-down distribution functions is equivalent to regularization using slowing-down physics as prior information. Using this technique in a Tikhonov formulation, we infer the particle density fractions for each NBI energy for each NBI beam from synthetic measurements, resulting in six unknowns at TCV and 24 unknowns at W7-X. Additionally, we show that installing 40 LOS in each of 17 ports at W7-X, providing full beam coverage and almost full angle coverage, produces the highest quality reconstructions.

Spherical tokamaks are key to the successful design of operating scenarios of future fusion reactors in the areas of divertor physics, neutral beam current drive and fast ion physics. MAST Upgrade, which has successfully concluded its first experimental campaign, was specifically designed to address the role of the radial gradient of the fast ion distribution in driving the excitation of magneto-hydrodynamic (MHD) instabilities, such as toroidal Alfven eigenmodes, fish-bones and long-lived mode, thanks to its two tangential neutral beam injection systems, one on the equatorial plane and one that is vertically shifted 65 cm above the equatorial plane. To study the fast ion dynamics in the presence of such instabilities, as well as of sawteeth and neo-classical tearing modes, several fast ion diagnostics were upgraded and new ones added. Among them, the MAST prototype neutron camera (NC) has been upgraded to six, equatorial sight-lines. The first observations of the confined fast ion behavior with the upgraded NC in a wide range of plasma scenarios characterized by on-axis and/or off-axis heating and different MHD instabilities are presented here. The observations presented in this study confirm previous results on MAST but with a higher level of detail and highlight new physics observations unique to the MAST Upgrade. The results presented here confirm the improved performance of the NC Upgrade, which thus becomes one of the key elements, in combination with the rich set of fast ion diagnostics available on the MAST Upgrade, for a more constrained modeling of the fast ion dynamics in fusion reactor relevant scenarios.

We compute reconstructions of 4D and 5D fast-ion phase-space distribution functions in fusion plasmas from synthetic projections of these functions. The fast-ion phase-space distribution functions originating from neutral beam injection (NBI) at TCV and Wendelstein 7-X (W7-X) at full, half, and one-third injection energies can be distinguished and particle densities of each component inferred based on 20 synthetic spectra of projected velocities at TCV and 680 at W7-X. Further, we demonstrate that an expansion into a basis of slowing-down distribution functions is equivalent to regularization using slowing-down physics as prior information. Using this technique in a Tikhonov formulation, we infer the particle density fractions for each NBI energy for each NBI beam from synthetic measurements, resulting in six unknowns at TCV and 24 unknowns at W7-X. Additionally, we show that installing 40 LOS in each of 17 ports at W7-X, providing full beam coverage and almost full angle coverage, produces the highest quality reconstructions.

Spherical tokamaks are key to the successful design of operating scenarios of future fusion reactors in the areas of divertor physics, neutral beam current drive and fast ion physics. MAST Upgrade, which has successfully concluded its first experimental campaign, was specifically designed to address the role of the radial gradient of the fast ion distribution in driving the excitation of magneto-hydrodynamic (MHD) instabilities, such as toroidal Alfven eigenmodes, fish-bones and long-lived mode, thanks to its two tangential neutral beam injection systems, one on the equatorial plane and one that is vertically shifted 65 cm above the equatorial plane. To study the fast ion dynamics in the presence of such instabilities, as well as of sawteeth and neo-classical tearing modes, several fast ion diagnostics were upgraded and new ones added. Among them, the MAST prototype neutron camera (NC) has been upgraded to six, equatorial sight-lines. The first observations of the confined fast ion behavior with the upgraded NC in a wide range of plasma scenarios characterized by on-axis and/or off-axis heating and different MHD instabilities are presented here. The observations presented in this study confirm previous results on MAST but with a higher level of detail and highlight new physics observations unique to the MAST Upgrade. The results presented here confirm the improved performance of the NC Upgrade, which thus becomes one of the key elements, in combination with the rich set of fast ion diagnostics available on the MAST Upgrade, for a more constrained modeling of the fast ion dynamics in fusion reactor relevant scenarios.

The reference ion cyclotron resonance frequency (ICRF) heating schemes for ITER deuterium-tritium (D-T) plasmas at the full magnetic field of 5.3 T are second harmonic heating of T and 3He minority heating. The wave-particle resonance location for these schemes coincide and are central at a wave frequency of 53 MHz at 5.3 T. Experiments have been carried out in the second major D-T campaign (DTE2) at JET, and in its prior D campaigns, to integrate these ICRF scenarios in JET high-performance plasmas and to compare their performance with the commonly used hydrogen (H) minority heating. In 50:50 D:T plasmas, up to 35% and 5% larger fusion power and diamagnetic energy content, respectively, were obtained with second harmonic heating of T as compared to H minority heating at comparable total input powers and gas injection rates. The core ion temperature was up to 30% and 20% higher with second harmonic T and 3He minority heating, respectively, with respect to H minority heating. These are favourable results for the use of these scenarios in ITER and future fusion reactors. According to modelling, adding ICRF heating to neutral beam injection using D and T beams resulted in a 10%-20% increase of on-axis bulk ion heating in the D-T plasmas due to its localisation in the plasma core. Central power deposition was confirmed with the break-in-slope and fast Fourier transform analysis of ion and electron temperature in response to ICRF modulation. The tail temperature of fast ICRF-accelerated tritons, their enhancement of the fusion yield and time behaviour as measured by an upgraded magnetic proton recoil spectrometer and neutral particle analyser were found in agreement with theoretical predictions. No losses of ICRF-accelerated ions were observed by fast ion detectors, which was as expected given the high plasma density of n e approximate to 7-8 x 1019 m-3 in the main heating phase that limited the formation of ICRF-accelerated fast ion tails. 3He was introduced in the mach

The reference ion cyclotron resonance frequency (ICRF) heating schemes for ITER deuterium-tritium (D-T) plasmas at the full magnetic field of 5.3 T are second harmonic heating of T and 3He minority heating. The wave-particle resonance location for these schemes coincide and are central at a wave frequency of 53 MHz at 5.3 T. Experiments have been carried out in the second major D-T campaign (DTE2) at JET, and in its prior D campaigns, to integrate these ICRF scenarios in JET high-performance plasmas and to compare their performance with the commonly used hydrogen (H) minority heating. In 50:50 D:T plasmas, up to 35% and 5% larger fusion power and diamagnetic energy content, respectively, were obtained with second harmonic heating of T as compared to H minority heating at comparable total input powers and gas injection rates. The core ion temperature was up to 30% and 20% higher with second harmonic T and 3He minority heating, respectively, with respect to H minority heating. These are favourable results for the use of these scenarios in ITER and future fusion reactors. According to modelling, adding ICRF heating to neutral beam injection using D and T beams resulted in a 10%-20% increase of on-axis bulk ion heating in the D-T plasmas due to its localisation in the plasma core. Central power deposition was confirmed with the break-in-slope and fast Fourier transform analysis of ion and electron temperature in response to ICRF modulation. The tail temperature of fast ICRF-accelerated tritons, their enhancement of the fusion yield and time behaviour as measured by an upgraded magnetic proton recoil spectrometer and neutral particle analyser were found in agreement with theoretical predictions. No losses of ICRF-accelerated ions were observed by fast ion detectors, which was as expected given the high plasma density of n e approximate to 7-8 x 1019 m-3 in the main heating phase that limited the formation of ICRF-accelerated fast ion tails. 3He was introduced in the mach

We compute reconstructions of 4D and 5D fast-ion phase-space distribution functions in fusion plasmas from synthetic projections of these functions. The fast-ion phase-space distribution functions originating from neutral beam injection (NBI) at TCV and Wendelstein 7-X (W7-X) at full, half, and one-third injection energies can be distinguished and particle densities of each component inferred based on 20 synthetic spectra of projected velocities at TCV and 680 at W7-X. Further, we demonstrate that an expansion into a basis of slowing-down distribution functions is equivalent to regularization using slowing-down physics as prior information. Using this technique in a Tikhonov formulation, we infer the particle density fractions for each NBI energy for each NBI beam from synthetic measurements, resulting in six unknowns at TCV and 24 unknowns at W7-X. Additionally, we show that installing 40 LOS in each of 17 ports at W7-X, providing full beam coverage and almost full angle coverage, produces the highest quality reconstructions.

We compute reconstructions of 4D and 5D fast-ion phase-space distribution functions in fusion plasmas from synthetic projections of these functions. The fast-ion phase-space distribution functions originating from neutral beam injection (NBI) at TCV and Wendelstein 7-X (W7-X) at full, half, and one-third injection energies can be distinguished and particle densities of each component inferred based on 20 synthetic spectra of projected velocities at TCV and 680 at W7-X. Further, we demonstrate that an expansion into a basis of slowing-down distribution functions is equivalent to regularization using slowing-down physics as prior information. Using this technique in a Tikhonov formulation, we infer the particle density fractions for each NBI energy for each NBI beam from synthetic measurements, resulting in six unknowns at TCV and 24 unknowns at W7-X. Additionally, we show that installing 40 LOS in each of 17 ports at W7-X, providing full beam coverage and almost full angle coverage, produces the highest quality reconstructions.

In recent years tokamak experiments and modelling have increasingly indicated that the interaction between suprathermal (fast) ions and thermal plasma can lead to a reduction of turbulence and an improvement of confinement. The regimes in which this stabilization occurs are relevant to burning plasmas, and their understanding will inform reactor scenario optimization. This review summarizes observations, simulations, theoretical understanding, and open questions on this emerging topic.

Spherical tokamaks are key to the successful design of operating scenarios of future fusion reactors in the areas of divertor physics, neutral beam current drive and fast ion physics. MAST Upgrade, which has successfully concluded its first experimental campaign, was specifically designed to address the role of the radial gradient of the fast ion distribution in driving the excitation of magneto-hydrodynamic (MHD) instabilities, such as toroidal Alfven eigenmodes, fish-bones and long-lived mode, thanks to its two tangential neutral beam injection systems, one on the equatorial plane and one that is vertically shifted 65 cm above the equatorial plane. To study the fast ion dynamics in the presence of such instabilities, as well as of sawteeth and neo-classical tearing modes, several fast ion diagnostics were upgraded and new ones added. Among them, the MAST prototype neutron camera (NC) has been upgraded to six, equatorial sight-lines. The first observations of the confined fast ion behavior with the upgraded NC in a wide range of plasma scenarios characterized by on-axis and/or off-axis heating and different MHD instabilities are presented here. The observations presented in this study confirm previous results on MAST but with a higher level of detail and highlight new physics observations unique to the MAST Upgrade. The results presented here confirm the improved performance of the NC Upgrade, which thus becomes one of the key elements, in combination with the rich set of fast ion diagnostics available on the MAST Upgrade, for a more constrained modeling of the fast ion dynamics in fusion reactor relevant scenarios.

The sawtooth instability is known for inducing transport and loss of energetic particles (EPs), and for generating seed magnetic islands that can trigger tearing modes. Both effects degrade the overall plasma performance. Several theories and numerical models have been previously developed to quantify the expected EP transport caused by sawteeth, with various degrees of sophistication to differentiate the response of EPs at different energies and on different orbits (e.g. passing vs. trapped), although the analysis is frequently limited to a single time slice during a tokamak discharge. This work describes the development and initial benchmark of a framework that enables a reduced model for EP transport by sawteeth retaining the full EP phase-space information. The model, implemented in the ORBIT hamiltonian particle-following code, can be used either as a standalone post-processor taking input data from codes such as TRANSP, or as a pre-processor to compute transport coefficients that can be fed back to TRANSP for time-dependent simulations including the effects of sawteeth on EPs. The advantage of the latter approach is that the evolution of the EP distribution can be simulated quantitatively for sawtoothing discharges, thus enabling a more accurate modeling of sources, sinks and overall transport properties of EP and thermal plasma species for comprehensive physics studies that require detailed information of the fast-ion distribution function and its evolution over time.

Thailand Tokamak-1 (TT-1) will be a small research tokamak being operated by the Thailand Institute of Nuclear Technology (TINT) in Nakhonnayok province, Thailand. The device partly uses the infrastructure of the former HT-6M alongside new hardware. Auxiliary heating sources, such as ion cyclotron range of frequency (ICRH) and/or neutral beam injection (NBI), maybe also included in subsequent operations. This study aims to investigate the characteristics and behavior of fast ions produced by NBI heating in various operation scenarios. This work employs the collisionless Lorentz-orbit (LORBIT) code for simulating the motion of ions in the TT-1. When a hydrogen ion (H+) with an energy of 20 keV travels in the plasma with the plasma current of 100 kA, the toroidal magnetic field of 1.0 T, the average Larmor radius is found to approximately be 0.48 cm for a passing transit particle and 2.17 cm for a trapped particle. The Larmor radius reduces as the plasma current is increased or a stronger magnetic field is used. Furthermore, when the ion has pitch angles in the range of 85°-100°, the ion orbit clearly demonstrates a banana shape with a width and height of 7.0 and 18.8 cm, respectively. The size of the banana orbit shows an inverse correlation with the plasma current. In the last part of this work, we investigate the motion of ions that had initial positions along the paths of a co-current and counter-current NBI. It is shown that the number of lost ions can be reduced if the initial radial positions of the ions are in the range of 0.64–0.68 m., source:https://doi.org/10.1016/j.fusengdes.2022.113254, identifier:0000-0002-0801-5106

0000-0002-0801-5106, Thailand Tokamak-1 (TT-1) will be a small research tokamak being operated by the Thailand Institute of Nuclear Technology (TINT) in Nakhonnayok province, Thailand. The device partly uses the infrastructure of the former HT-6M alongside new hardware. Auxiliary heating sources, such as ion cyclotron range of frequency (ICRH) and/or neutral beam injection (NBI), maybe also included in subsequent operations. This study aims to investigate the characteristics and behavior of fast ions produced by NBI heating in various operation scenarios. This work employs the collisionless Lorentz-orbit (LORBIT) code for simulating the motion of ions in the TT-1. When a hydrogen ion (H+) with an energy of 20 keV travels in the plasma with the plasma current of 100 kA, the toroidal magnetic field of 1.0 T, the average Larmor radius is found to approximately be 0.48 cm for a passing transit particle and 2.17 cm for a trapped particle. The Larmor radius reduces as the plasma current is increased or a stronger magnetic field is used. Furthermore, when the ion has pitch angles in the range of 85°-100°, the ion orbit clearly demonstrates a banana shape with a width and height of 7.0 and 18.8 cm, respectively. The size of the banana orbit shows an inverse correlation with the plasma current. In the last part of this work, we investigate the motion of ions that had initial positions along the paths of a co-current and counter-current NBI. It is shown that the number of lost ions can be reduced if the initial radial positions of the ions are in the range of 0.64–0.68 m.

Both fast ions and runaway electrons are described by distribution functions, the understanding of which are of critical importance for the success of future fusion devices such as ITER. Typically, energetic particle diagnostics are only sensitive to a limited subsection of the energetic particle phase-space which is often insufficient for model validation. However, previous publications show that multiple measurements of a single spatially localized volume can be used to reconstruct a distribution function of the energetic particle velocity-space by using the diagnostics' velocity-space weight functions, i.e. velocity-space Tomography. In this work we use the recently formulated orbit weight functions to remove the restriction of spatially localized measurements and present orbit tomography, which is used to reconstruct the 3D phase-space distribution of all energetic particle orbits in the plasma. Through a transformation of the orbit distribution, the full energetic particle distribution function can be determined in the standard {energy, pitch, r, z}-space. We benchmark the technique by reconstructing the fast-ion distribution function of an MHD-quiescent DIII-D discharge using synthetic and experimental FIDA measurements. We also use the method to study the redistribution of fast ions during a sawtooth crash at ASDEX upgrade using FIDA measurements. Finally, a comparison between the orbit tomography and velocity-space tomography is shown.

Collisionless physics primarily determines the transport of fusion-born alpha particles in 3D equilibria. Several transport mechanisms have been implicated in stellarator configurations, including stochastic diffusion due to class transitions, ripple trapping, and banana drift-convective orbits. Given the guiding center dynamics in a set of six quasihelical and quasiaxisymmetric equilibria, we perform a classification of trapping states and transport mechanisms. In addition to banana drift convection and ripple transport, we observe substantial non-conservation of the parallel adiabatic invariant which can cause losses through diffusive banana tip motion. Furthermore, many lost trajectories undergo transitions between trapping classes on longer time scales, either with periodic or irregular behavior. We discuss possible optimization strategies for each of the relevant transport mechanisms. We perform a comparison between fast ion losses and metrics for the prevalence of mechanisms such as banana-drift convection (Velasco et al 2021 Nucl. Fusion 61 116059), transitioning orbits, and wide orbit widths. Quasihelical configurations are found to have natural protection against ripple-trapping and diffusive banana tip motion leading to a reduction in prompt losses.

A conceptual design of a scintillator-based fast-ion loss detector (FILD) has been developed for the Wendelstein 7-X stellarator (W7-X). Simulations using the Monte Carlo codes ASCOT5 and BEAMS3D have been performed to calculate the expected flux of neutral beam injection (NBI)-generated fast hydrogen ions onto the conceptual detector probe head. These fast-ion loss fluxes have been calculated for several magnetic field configurations as well as probe insertion positions. At the maximum insertion position, both co- and counter-going losses with high incident pitch angles are observed; however, at retracted positions, only co-going fast ions reach the probe head. The FILDSIM code has been used to optimize the geometry of the detector entrance and collimating elements to achieve a wide velocity space coverage as well as a high velocity-space resolution. A synthetic FILD signal is calculated for the expected loss distribution via forward modeling using the instrument response function. The synthetic signal is found to largely retain the velocity space features of the loss distribution.

Using the PION ICRH modelling code and comparisons against JET tokamak experiments, the effect of including pitch angle dependence within the RF diffusion operator on the fast ion particle distribution functions is quantified. It is found to be of greatest importance in cases of higher harmonic heating and lower heating ion mass, resulting in faster drop-off of the distribution's high energy tail. We see differences of several orders of magnitude in the high-energy range and significant non-linear alterations by several tens of percent to ion species power partition. ITER scenario operational parameters are also considered, and this improved treatment is shown to benefit anticipated ITER scenarios with second harmonic hydrogen heating, according to our predictions. PION's combination of benchmarked simplified wave physics and Fokker-Planck treatment offers modelling advantages. Since including the pitch angle dependence in the RF diffusion operator has not led to a significant increase in the required computing time when modelling different ICRF schemes in JET discharges, it has been made available within the production code., The CCFE part of 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 and 2019–2020 under Grant Agreement No. 633053. The views and opinions expressed herein do not necessarily reflect those of the European Commission. The BSC part of this project is co-financed by the European Union Regional Development Fund within the framework of the ERDF Operational Program of Catalonia 2014–2020 with a grant of 50% of total cost eligible. The authors are grateful to Jacob Eriksson for assistance with experimental data, to Lars-Göran Eriksson for discussions on the implementation of the new features, and to Colin Roach and Michael Fitzgerald for valuable comments on the manuscript., Peer Reviewed, Postprint (published version)

Fast ion wall loads can result in excessively high heat fluxes to the plasma-facing components (PFCs). To allow for the development of mitigation strategies, and thereby protect the PFCs, the fast ion losses have to be predicted by faithful models. To ensure that fast ion models are an accurate representation of the real world, they need to be verified. The neutral-beam experiments performed in Wendelstein 7-X (W7-X) allow to investigate and verify models of the fast ion losses in the stellarator configuration. Infrared thermographic measurements were used to obtain the heat flux to both the baffle plates and the divertor. We found evidence of fast ion wall loads on the baffle plates, with loads between 100 kW m−2 and 1 MW m−2. The loads are attributed to fast ions which escape the main plasma via magnetic ripples. The fast ion wall loads on the baffle plates show up-down and toroidal asymmetry. The experimental results were compared to numerical simulations performed by the BEAMS3D and ASCOT codes. Qualitative agreement in up-down asymmetry is found, but the magnitude and toroidal asymmetry are not yet well predicted by the simulations. The asymmetries of the strike lines on the divertor suggest that fast ions also play a role here. Specifically, a second strike line emerged consistently in the high-iota configuration on the horizontal divertor. The shape and magnitude of the strike lines changed considerably during the neutral-beam injection (NBI) operation phase. Although no damage to steel components of W7-X was found, fast ion loads to the baffle plates could possibly limit the NBI operation in the upcoming campaigns of W7-X.

A discrepancy between predicted and measured neutron rates on MAST using TRANSP/NUBEAM has previously been observed and a correction factor of about 0.6 was needed to match the two: this correction factor could not be accounted for by the experimental uncertainties in the plasma kinetic profiles nor in the NBI energy and power (Cecconello et al 2019 Nucl. Fusion 59 016006). Further causes of this discrepancy are here studied by means of TRANSP/NUBEAM and ASCOT/BBNBI simulations. Different equilibria, toroidal field ripples, uncertainties on the NBI divergence value and gyro-orbit effects were studied and simulations were performed with both transport codes. It was found that the first three effects accounted for only a 5% variation in the fast ion density. On the other hand, full gyro-orbit simulations of the fast ions dynamics carried out in ASCOT/BBNBI resulted in an approximately 20% reduction of the fast ion population compared to TRANSP/NUBEAM. A detailed analysis of the fast ion distributions showed how the drop occurred regardless of the energy at pitch values <=-0.4. The DRESS code was then used to calculate the neutron rate at the neutron camera detector's location showing that the discrepancy is considerably reduced when the full gyro-orbit fast ion distribution is used, with now the correction factor, used to match experimental and predicted neutron rates, being around 0.9.

An integrated modelling methodology for the calculation of realistic plasma neutron sources for the JET tokamak has been developed. The computational chain comprises TRANSP plasma transport and DRESS neutron spectrum calculations, and their coupling to the MCNP neutron transport code, bridging plasma physics and neutronics. In the paper we apply the developed methodology to the analysis of neutron emission properties of deuterium and helium plasmas at JET, and validate individual modelling steps against neutron diagnostic measurements. Two types of JET discharges are modelled-baseline-like and three-ion radio-frequency scenarios-due to their diversity in plasma heating, characteristics of the induced fast ion population, and the imprint of these on neutron emission properties. The neutron emission modelling results are quantitatively compared to the total neutron yield from fission chambers, neutron emissivity profiles from the neutron camera, neutron spectra from the time-of-flight spectrometer, and neutron activation measurements. The agreement between measured and calculated quantities is found to be satisfactory for all four diagnostic systems within the estimated experimental and computational uncertainties. Additionally, the effect of neutrons not originating from the dominating D(D, n)He-3 reactions is studied through modelling of triton burnup DT neutrons, and, in mixed D-He-3 plasmas, neutrons produced in the Be-9(D, n gamma)B-10 reaction on impurities. It is found that these reactions can contribute up to several percent to the total neutron yield and dominate the neutron activation of samples. The effect of MeV-range fast ions on the neutron activation of In-115 and Al-27 samples is measured and computationally validated.

An understanding of fast (supra-thermal) ion behavior is of great importance in tokamak physics and is a subject studied from both theoretical and experimental perspectives. This paper investigates the spatial energy and density distributions of RF accelerated deuterons using the neutron camera at the tokamak JET. Using the 19 liquid scintillator detectors available in the neutron camera system, we obtain spatial information that cannot be accessed with a single sightline. We present a spectroscopic analysis method in which a spatially resolved model of the fast ion distribution is fitted to the pulse height spectra from all 19 detectors simultaneously. The fast ion distribution is parameterized in such a way that the density, energy, and pitch-angle parts are uncoupled. The energy part is composed of a Maxwellian distribution, characterized by an 'RF tail temperature,' and the spatial dependence is modeled as a two-dimensional Gaussian distribution on the poloidal plane of the tokamak. From this parameterized model, we can calculate the spectrum of fusion born neutrons originating from reactions involving RF accelerated deuterons, and by fitting this model to the measured neutron camera pulse height spectra, we obtain an estimate of the spatially resolved distribution of the fast deuterons. The method has been applied to three JET pulses using different RF heating schemes and is shown to identify several features of the fast ion distribution produced in the various scenarios. Hence, this method is able to provide quantitative information about the fast ion distribution resulting from different heating schemes, and can also be useful, e.g., to validate simulation results from RF modeling codes.

Recurring bursts of chirping Alfvén modes that were observed in JT-60U tokamak plasmas driven by negative-ion-based neutral beams (N-NB) are reproduced in first-principle simulations performed with an extended version of the hybrid code MEGA. This code simulates the interactions between gyrokinetic fast ions and magnetohydrodynamic (MHD) modes in the presence of a realistic fast ion source and collisions, so that it self-consistently captures dynamics across a wide range of time scales (0.01–100 ms). The simulation confirms that the experimentally observed phenomena known as 'fast frequency sweeping (fast FS) modes' are caused by bursts of energetic particle modes (EPM) with dominant toroidal mode number n = 1. On the long time scale (1–10 ms), the simulation reproduces the chirping range (40–60 kHz), the burst duration (few ms) and intervals (5–10 ms). On the short time scale (0.01–0.1 ms), it reproduces pulsations and phase jumps, which we interpret as the result of beating between multiple resonant wave packets. Having reproduced at multiple levels of detail the dynamics of low-amplitude long-wavelength Alfvén modes driven by N-NB ions, the next goal is to reproduce and explain abrupt large-amplitude events (ALE) that were seen in the same experiments at longer time intervals (10–100 ms)., source:Citation A. Bierwage et al 2017 Nucl. Fusion 57 016036, source:http://dx.doi.org/10.1088/1741-4326/57/1/016036, identifier:0000-0003-1243-0502

Modelling neutral beam injection (NBI) in fusion reactors requires computing the trajectories of large ensembles of particles. Slowing down times of up to one second combined with nanosecond time steps make these simulations computationally very costly. This paper explores the performance of BGSDC, a new numerical time stepping method, for tracking ions generated by NBI in the DIII-D and JET reactors. BGSDC is a high-order generalisation of the Boris method, combining it with spectral deferred corrections and the Generalized Minimal Residual method GMRES. Without collision modelling, where numerical drift can be quantified accurately, we find that BGSDC can deliver higher quality particle distributions than the standard Boris integrator at comparable cost or comparable distributions at lower cost. With collision models, quantifying accuracy is difficult but we show that BGSDC produces stable distributions at larger time steps than Boris.

JT-60SA, the largest tokamak that will operate before ITER, has been designed and built jointly by Japan and Europe, and is due to start operation in 2020. Its main missions are to support ITER exploitation and to contribute to the demonstration fusion reactor machine and scenario design. Peculiar properties of JT-60SA are its capability to produce long-pulse, high-β, and highly shaped plasmas. The preparation of the JT-60SA Research Plan, plasma scenarios, and exploitation are producing physics results that are not only relevant to future JT-60SA experiments, but often constitute original contributions to plasma physics and fusion research. Results of this kind are presented in this paper, in particular in the areas of fast ion physics, highbeta plasma properties and control, and non-linear edge localised mode stability studies.

A comparison between fast ion measurements and sawtooth models in the Mega Amp Spherical Tokamak (MAST) is extended to include fast ion deuterium-alpha (FIDA) data. It is concluded that FIDA data cannot be used to distinguish between three alternative models used in the plasma transport/fast particle code TRANSP/NUBEAM to simulate fast ion redistribution during sawteeth. For FIDA lines-of-sight that probe the sawtoothing region, at each sawtooth crash there is an overall drop in the emission of up to 60%. Data from passive FIDA lines-of-sight (i.e.with emission resulting from neutralisation by thermal neutrals in the plasma periphery rather than beam neutrals) show a sudden increase in the emission following sawtooth crashes. The subsequent decay in the emission in these passive channels indicates that redistributed passing fast ions are rapidly lost from the edge region, probably as a result of charge-exchange reactions with edge neutrals.

JT-60SA, the largest tokamak that will operate before ITER, has been designed and built jointly by Japan and Europe, and is due to start operation in 2020. Its main missions are to support ITER exploitation and to contribute to the demonstration fusion reactor machine and scenario design. Peculiar properties of JT-60SA are its capability to produce long-pulse, high-ß, and highly shaped plasmas. The preparation of the JT-60SA Research Plan, plasma scenarios, and exploitation are producing physics results that are not only relevant to future JT-60SA experiments, but often constitute original contributions to plasma physics and fusion research. Results of this kind are presented in this paper, in particular in the areas of fast ion physics, high-beta plasma properties and control, and non-linear edge localised mode stability studies., Postprint (published version)

We demonstrate 2D reconstructions of the fast-ion velocity distribution from two-view fast-ion D-alpha (FIDA) measurements at the EAST tokamak. By expressing the distribution in a basis relying on the fast-ion slowing-down process in fusion plasmas, the smoothness and velocity-space resolution of reconstructions are improved. We reconstruct distributions of fast ions born from simultaneous co- and counter-current neutral beam injection and detect the expected distinct change in fast-ion birth pitch when comparing discharges utilizing different neutral beam injectors. For purely co-current injection, we find a good agreement between TRANSP-predicted and reconstructed fast-ion densities, pressures and current densities for energies above 20 keV. We furthermore illustrate the improvement of the reconstructed high-energy range (> 40 keV) of the distribution by combining FIDA with neutron emission spectroscopy measurements with the compact single-plate EJ301 scintillator.

Understanding the effect of Alfvén eigenmodes (AEs) and neoclassical tearing modes (NTMs) on fast ions is highly important for fusion reactors due to potentially strong resonant interactions between the fast ions and the modes. Here, we use the four-view fast-ion D-alpha (FIDA) diagnostic installed in the DIII-D tokamak to reconstruct the fast-ion velocity distribution at two radial positions during two sequential discharges with strong and weak mode activity, respectively. The velocity-space coverage of the diagnostics, however, only allows reliable reconstructions of fast ions with positive pitches. Therefore, we suggest new tomographic inversion methods relying on prior information outside the well-diagnosed region. We find that within the population of fast ions with positive pitches, ions, at all energies, are transported away from the measurement volumes. Comparisons between the reconstructions and kick model simulations, where the mode activity is considered, reveal that low-frequency modes such as the NTMs and low-frequency AEs contribute significantly to the positive-pitch fast-ion transport in the central measurement volume, whereas TAEs and EAEs become important farther out and are responsible for decreased fast-ion confinement.

In plasmas heated with deuterium beams a deficit of the expected fusion neutron rate is an indicator of the deterioration of the fast-ion confinement, caused, for instance, by magnetohydrodynamic instabilities. The capability of predicting this deficit during the discharge relies on the availability of real-time estimates of the neutron rate from NBI codes which must be fast and accurate at the same time. Therefore, the recently developed real-time RABBIT code for neutral beam injection (NBI) simulations has been extended to output the distribution function and calculate the neutron emission. After the description of this newly installed diagnostics in RABBIT, benchmarks with NUBEAM, a massively used and validated Monte Carlo NBI solver, are discussed on ASDEX-Upgrade and JET cases. A first application for control-room intershot analysis on DIII-D is presented, and the results are compared on a large database with a slower NUBEAM analysis. Further application possibilities, e.g. for real-time control of Alfven eigenmodes, are outlined., For complete list of authors see http://dx.doi.org/10.1088/1741-4326/ab1edd

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., For complete list of authors see http://dx.doi.org/10.1088/1361-6587/ab1e65

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., For complete list of authors see http://dx.doi.org/10.1088/1741-4326/ab2d4f

Fast ions in fusion plasmas often leave characteristic signatures in the neutron emission from the plasma. In this paper, we show how neutron measurements can be used to study fast ions and give examples of physics results obtained on present day tokamaks. The focus is on measurements with dedicated neutron spectrometers and with compact neutron detectors used in each channel of neutron profile monitors. A measured neutron spectrum can be analyzed in several different ways, depending on the physics scenario under consideration. Gross features of a fast ion energy distribution can be studied by applying suitably chosen thresholds to the measured spectrum, thus probing ions with different energies. With this technique it is possible to study the interaction between fast ions and MHD activity, such as toroidal Alfven eigenmodes (TAEs) and sawtooth instabilities. Quantitative comparisons with modeling can be performed by a direct computation of the neutron emission expected from a given fast ion distribution. Within this framework it is also possible to determine physics parameters, such as the supra-thermal fraction of the neutron emission, by fitting model parameters to the data. A detailed, model-independent estimate of the fast ion distribution can be obtained by analyzing the data in terms of velocity space weight functions. Using this method, fast ion distributions can be resolved in both energy and pitch by combining neutron and gamma-ray measurements obtained along several different sightlines. Fast ion measurements of the type described in this paper will also be possible at ITER, provided that the spectrometers have the dynamic range required to resolve the fast ion spectral features in the presence of the dominating thermonuclear neutron emission. A dedicated high-resolution neutron spectrometer has been designed for this purpose.

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., QC 20200313

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., QC 20200311

In plasmas heated with deuterium beams a deficit of the expected fusion neutron rate is an indicator of the deterioration of the fast-ion confinement, caused, for instance, by magnetohydrodynamic instabilities. The capability of predicting this deficit during the discharge relies on the availability of real-time estimates of the neutron rate from NBI codes which must be fast and accurate at the same time. Therefore, the recently developed real-time RABBIT code for neutral beam injection (NBI) simulations has been extended to output the distribution function and calculate the neutron emission. After the description of this newly installed diagnostics in RABBIT, benchmarks with NUBEAM, a massively used and validated Monte Carlo NBI solver, are discussed on ASDEX-Upgrade and JET cases. A first application for control-room intershot analysis on DIII-D is presented, and the results are compared on a large database with a slower NUBEAM analysis. Further application possibilities, e.g. for real-time control of Alfven eigenmodes, are outlined., QC 20200304

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., For complete list of authors see http://dx.doi.org/10.1088/1361-6587/ab1e65

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., For complete list of authors see http://dx.doi.org/10.1088/1361-6587/ab1e65

The interaction between fast-ions and edge localized modes (ELMs) is investigated by means of fast-ion loss detector measurements. Fast-ion losses are increased during ELMs exhibiting a 3D filamentary-like behaviour. An accelerated beam ion population has been observed during ELMs in a tokamak for the first time. Tomographic inversion of the measured fast-ion losses reveal multiple velocity-space structures. Attending to the experimental observations, an acceleration mechanism is proposed based on a resonant interaction between the beam ions and parallel electric fields emerging during the ELM crash. The key experimental observations can be qualitatively reproduced by full-orbit following simulations of fast-ions in the presence of the ELM magnetic and electric perturbation fields. Our findings may shed light on the possible contribution of fast-ions to the ELM stability and the transient heat loads on plasma facing components.

The interaction between fast-ions and edge localized modes (ELMs) is investigated by means of fast-ion loss detector measurements. Fast-ion losses are increased during ELMs exhibiting a 3D filamentary-like behaviour. An accelerated beam ion population has been observed during ELMs in a tokamak for the first time. Tomographic inversion of the measured fast-ion losses reveal multiple velocity-space structures. Attending to the experimental observations, an acceleration mechanism is proposed based on a resonant interaction between the beam ions and parallel electric fields emerging during the ELM crash. The key experimental observations can be qualitatively reproduced by full-orbit following simulations of fast-ions in the presence of the ELM magnetic and electric perturbation fields. Our findings may shed light on the possible contribution of fast-ions to the ELM stability and the transient heat loads on plasma facing components.