Search

The performance of fusion devices relies strongly on the good confinement of energetic particles (EPs). Therefore, the investigation of EP transport by magnetohydrodynamic instabilities is one of the key aspects in the development of plasma scenarios. Alfvénic instabilities in particular can lead to significant losses of alpha particles that are essential for plasma self-heating. A so-called afterglow scheme has been developed to study the destabilization of Alfvén eigenmodes (AEs) by alpha particles and associated EP transport in the JET tokamak. In this work, the linear stability of AEs is discussed for the partial afterglow phase in a JET deuterium plasma discharge and for the full afterglow phase in a projected deuterium–tritium (DT) plasma. Thanks to recent upgrades in the tokamak transport code TRANSP, one can account for the contributions of different EP species to mode stability. Analysis of deuterium plasmas shows that AE growth rates are extremely sensitive to the energy and distribution of fast ions. An increase in fast ion energy can lead to more unstable AEs. In the afterglow phase of projected DT plasmas, it is EPs that mostly drive the AEs. However, the drive by alpha particles is comparable to that by beam ions and their contribution to the net growth rate might be hard to separate. According to the discussed projections, the destabilization of AEs might be ineffective because the background plasma damping significantly exceeds the EP drive. In this case, the development of an alternative plasma scenario that allows us to overcome such damping would be required in future experiments.

Alfvénic instabilities (AEs) are well known to cause enhanced transport of energetic particles (EPs) in fusion devices. Most studies until now have focused on characterizing and understanding AE stability in single-species plasmas heated by neutral beams (NB), where deuterium is typically used as both main plasma species and NB fuel. As the fusion community moves toward fusion reactors that target burning plasma conditions, such as ITER, the single-species picture breaks down. Burning plasmas, which will use a mix of deuterium and tritium (DT) as main fuel, also feature the presence of several supra-thermal fusion products such as alpha particles, protons, helium isotopes and high-energy tritium ions. This work presents the extension of the EP transport kick model implemented in the TRANSP time-dependent tokamak transport code to study the combined effect of multiple EP species on AE stability and, in turn, the response of different EP species to plasma instabilities in terms of their redistribution and losses. Further validation of the enhanced model is planned based on experimental results expected from the JET DT campaign scheduled for 2021, in preparation for ITER plasmas and beyond.

openaire: EC/H2020/633053/EU//EUROfusion The recently developed directional relativistic spectrum simulator (DRESS) code has been validated for the first time against numerical calculations and experimental measurements performed on the Mega Ampere Spherical Tokamak . In this validation, the neutron emissivities and rates computed by DRESS are benchmarked against TRANSP/NUBEAM predictions while the neutron energy spectra provided by DRESS taking as input TRANSP/NUBEAM and ASCOT/BBNBI in Gyro-Orbit mode fast ion distributions are validated against proton pulse height spectra measured by the neutron flux monitor. Excellent agreement was found between DRESS and TRANSP/NUBEAM predictions of local and total neutron emission.

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.

New results from MAST are presented that focus on validating models in order to extrapolate to future devices. Measurements during start-up experiments have shown how the bulk ion temperature rise scales with the square of the reconnecting field. During the current ramp up models are not able to correctly predict the current diffusion. Experiments have been performed looking at edge and core turbulence. At the edge detailed studies have revealed how filament characteristic are responsible for determining the near and far SOL density profiles. In the core the intrinsic rotation and electron scale turbulence have been measured. The role that the fast ion gradient has on redistributing fast ions through fishbone modes has led to a redesign of the neutral beam injector on MAST Upgrade. In H-mode the turbulence at the pedestal top has been shown to be consistent with being due to electron temperature gradient modes. A reconnection process appears to occur during ELMs and the number of filaments released determines the power profile at the divertor. Resonant magnetic perturbations can mitigate ELMs provided the edge peeling response is maximised and the core kink response minimised. The mitigation of intrinsic error fields with toroidal mode number n>1 has been shown to be important for plasma performance., 34 pages, 10 figures. This is an author-created, un-copyedited version of an article submitted for publication in Nuclear Fusion. IoP Publishing Ltd is not responsible for any errors or omissions in this version of the manuscript or any version derived from it

Neutral beam injection (NBI) will be one of the main sources of heating and non-inductive current drive in ITER. Due to high level of injected power the beam induced heat loads present a potential threat to the integrity of the first wall of the device, particularly in the presence of non-axisymmetric perturbations of the magnetic field. Neutral beam injection can also destabilize Alfven eigenmodes and energetic particle modes, and act as a source of plasma rotation. Therefore, reliable and accurate simulation of NBI is important for making predictions for ITER, as well as for any other current or future fusion device. This paper introduces a new beamlet-based neutral beam ionization model called BBNBI. It takes into account the fine structure of the injector, follows the injected neutrals until ionization, and generates a source ensemble of ionized NBI test particles for slowing down calculations. BBNBI can be used as a stand-alone model but together with the particle following code ASCOT it forms a complete and sophisticated tool for simulating neutral beam injection. The test particle ensembles from BBNBI are found to agree well with those produced by PENCIL for JET, and those produced by NUBEAM both for JET and ASDEX Upgrade plasmas. The first comprehensive comparisons of beam slowing down profiles of interest from BBNBI + ASCOT with results from PENCIL and NUBEAM/TRANSP, for both JET and AUG, are presented. It is shown that, for an axisymmetric plasma, BBNBI + ASCOT and NUBEAM agree remarkably well. Together with earlier 3D studies, these results further validate using BBNBI + ASCOT also for studying phenomena that require particle following in a truly three-dimensional geometry.

The National Spherical Torus Experiment (NSTX) is currently being upgraded to operate at twice the toroidal field and plasma current (up to 1 T and 2 MA), with a second, more tangentially aimed neutral beam (NB) for current and rotation control, allowing for pulse lengths up to 5 s. Recent NSTX physics analyses have addressed topics that will allow NSTX-Upgrade to achieve the research goals critical to a Fusion Nuclear Science Facility. These include producing stable, 100% non-inductive operation in high-performance plasmas, assessing plasma–material interface (PMI) solutions to handle the high heat loads expected in the next-step devices and exploring the unique spherical torus (ST) parameter regimes to advance predictive capability. Non-inductive operation and current profile control in NSTX-U will be facilitated by co-axial helicity injection (CHI) as well as radio frequency (RF) and NB heating. CHI studies using NIMROD indicate that the reconnection process is consistent with the 2D Sweet–Parker theory. Full-wave AORSA simulations show that RF power losses in the scrape-off layer (SOL) increase significantly for both NSTX and NSTX-U when the launched waves propagate in the SOL. Toroidal Alfven eigenmode avalanches and higher frequency Alfven eigenmodes can affect NB-driven current through energy loss and redistribution of fast ions. The inclusion of rotation and kinetic resonances, which depend on collisionality, is necessary for predicting experimental stability thresholds of fast growing ideal wall and resistive wall modes. Neutral beams and neoclassical toroidal viscosity generated from applied 3D fields can be used as actuators to produce rotation profiles optimized for global stability. DEGAS-2 has been used to study the dependence of gas penetration on SOL temperatures and densities for the MGI system being implemented on the Upgrade for disruption mitigation. PMI studies have focused on the effect of ELMs and 3D fields on plasma detachment and heat flux handling. Simulations indicate that snowflake and impurity seeded radiative divertors are candidates for heat flux mitigation in NSTX-U. Studies of lithium evaporation on graphite surfaces indicate that lithium increases oxygen surface concentrations on graphite, and deuterium–oxygen affinity, which increases deuterium pumping and reduces recycling. In situ and test-stand experiments of lithiated graphite and molybdenum indicate temperature-enhanced sputtering, although that test-stand studies also show the potential for heat flux reduction through lithium vapour shielding. Non-linear gyro kinetic simulations have indicated that ion transport can be enhanced by a shear-flow instability, and that non-local effects are necessary to explain the observed rapid changes in plasma turbulence. Predictive simulations have shown agreement between a microtearing-based reduced transport model and the measured electron temperatures in a microtearing unstable regime. Two Alfven eigenmode-driven fast ion transport models have been developed and successfully benchmarked against NSTX data. Upgrade construction is moving on schedule with initial physics research operation of NSTX-U planned for mid-2015.

A conceptual, theoretical assessment of using the toroidal fusion reactor, tokamak, for deep space interplanetary and interstellar missions is presented. Toroidal thermonuclear fusion reactors, such as tokamaks and stellarators, are unique for space propulsion, allowing a design with the magnetic configuration localized inside the toroidal magnetic field coils. Plasma energetic ions, including charged fusion products, can escape such closed configuration at certain conditions as a result of vertical drift in the toroidal rippled magnetic field. Escaping particles can be used for direct propulsion (because toroidal drift is directed one way vertically) or to create and heat externally confined plasma, so that the latter can be used for propulsion. In contrast to other fusion concepts proposed for space propulsion, this concept utilizes the natural drift motion of charged particles out of the closed magnetic field configuration. Also, using deuterium-tritium (D-T) plasma is a novel way to use fusion neutrons with the energy of 14 MeV for direct propulsion (neutron rocket) for out of solar system missions. A special design of the blanket of the reactor allows neutrons to escape the device in a preferable direction. This provides a direct (partial) conversion of fusion energy into thrust.

Alfvenic instabilities (AEs) are well known as a potential cause of enhanced fast ion transport in fusion devices. Given a specific plasma scenario, quantitative predictions of (i) expected unstable AE spectrum and (ii) resulting fast ion transport are required to prevent or mitigate the AE-induced degradation in fusion performance. Reduced models are becoming an attractive tool to analyze existing scenarios as well as for scenario prediction in time-dependent simulations. In this work, a neutral beam heated NSTX discharge is used as reference to illustrate the potential of a reduced fast ion transport model, known as kick model, that has been recently implemented for interpretive and predictive analysis within the framework of the time-dependent tokamak transport code TRANSP. Predictive capabilities for AE stability and saturation amplitude are first assessed, based on given thermal plasma profiles only. Predictions are then compared to experimental results, and the interpretive capabilities of the model further discussed. Overall, the reduced model captures the main properties of the instabilities and associated effects on the fast ion population. Additional information from the actual experiment enables further tuning of the model's parameters to achieve a close match with measurements.

Toroidicity-induced Alfven eigenmode (TAE) stability in the National Spherical Torus Experiment (NSTX) [S. M. Kaye, M. Ono, Y.-K. M. Peng et al., Fusion Technol. 36, 16 (1999)] is analyzed using the improved NOVA-K code [N. N. Gorelenkov, C. Z. Cheng, and G. Y. Fu, Phys. Plasmas 6, 2802 (1999)], which includes finite orbit width and Larmor radius effects and is able to predict the saturation amplitude for the mode using the quasilinear theory. A broad spectrum of unstable global TAEs with different toroidal mode numbers is predicted. Due to the strong poloidal field and the presence of the magnetic well in NSTX, better particle confinement in the presence of TAEs in comparison with tokamaks is illustrated making use of the ORBIT code [R. B. White and M. S. Chance, Phys. Fluids 27, 2455 (1984)].

Stability criterion for neoclassical tearing modes is obtained from the drift kinetic equation. A finite amplitude of a magnetic island is required for mode excitation. The threshold is determined by the ratio of the transversal and the parallel transport near the island when the flattening of the pressure profile eliminates the bootstrap current. A number of supershots from the database of the Tokamak Fusion Test Reactor (TFTR) [D. J. Grove and D. M. Meade, Nucl. Fusion 25, 1167 (1985)] are compared with the theory. In cases where the modes were observed in experiment the stability criterion was violated.

Tokamak plasmas can feature a large population of energetic particles (EP) from neutral beam injection or fusion reactions. In turn, energetic particles can drive instabilities, which affect the driving EP population leading to a distortion of the original EP distribution function and of quantities that depend on it. The latter include, for example, neutral beam (NB) current drive and plasma heating through EP thermalization. Those effects must be taken into account to enable reliable and quantitative simulations of discharges for present devices as well as predictions for future burning plasmas. Reduced models for EP transport are emerging as an effective tool for long time-scale integrated simulations of tokamak plasmas, possibly including the effects of instabilities on EP dynamics. Available models differ in how EP distribution properties are modified by instabilities, e.g. in terms of gradients in real or phase space. It is therefore crucial to assess to what extent different assumptions in the transport models affect predicted quantities such as EP profile, energy distribution, NB driven current and energy/momentum transfer to the thermal populations. A newly developed kick model, which includes modifications of the EP distribution by instabilities in both real and velocity space, is used in this work to investigate these issues. Coupled to TRANSP simulations, the kick model is used to analyze NB-heated NSTX and DIII-D discharges featuring unstable Alfvén eigenmodes (AEs). Results show that instabilities can strongly affect the EP distribution function, and modifications propagate to macroscopic quantities such as NB-driven current profile and NB power transferred to the thermal plasma species. Those important aspects are only qualitatively captured by simpler fast ion transport models that are based on radial diffusion of energetic ions only.

Integrated simulations of tokamak discharges typically rely on classical physics to model energetic particle (EP) dynamics. However, there are numerous cases in which energetic particles can suffer additional transport that is not classical in nature. Examples include transport by applied 3D magnetic perturbations and, more notably, by plasma instabilities. Focusing on the effects of instabilities,ad-hocmodels can empirically reproduce increased transport, but the choice of transport coefficients is usually somehow arbitrary. New approaches based on physics-based reduced models are being developed to address those issues in a simplified way, while retaining a more correct treatment of resonant wave-particle interactions. The kick model implemented in the tokamaktransport code TRANSP is an example of such reduced models. It includes modifications of the EP distribution by instabilities in real and velocity space, retaining correlations between transport in energy and space typical of resonant EP transport. The relevance of EP phase space modifications by instabilities is first discussed in terms of predicted fast ion distribution. Results are compared with those from a simple, ad-hoc diffusive model. It is then shown that the phase-space resolved model can also provide additional insight into important issues such as internal consistency of the simulations and mode stability through the analysis of the power exchanged between energetic particles and the instabilities.

A 3D halo neutral code developed at the Princeton Plasma Physics Laboratory and implemented for analysis using the TRANSP code is applied to projected National Spherical Torus eXperiment-Upgrade (NSTX-U plasmas). The legacy TRANSP code did not handle halo neutrals properly since they were distributed over the plasma volume rather than remaining in the vicinity of the neutral beam footprint as is actually the case. The 3D halo neutral code uses a 'beam-in-a-box' model that encompasses both injected beam neutrals and resulting halo neutrals. Upon deposition by charge exchange, a subset of the full, one-half and one-third beam energy components produce first generation halo neutrals that are tracked through successive generations until an ionization event occurs or the descendant halos exit the box. The 3D halo neutral model and neutral particle analyzer (NPA) simulator in the TRANSP code have been benchmarked with the Fast-Ion D-Alpha simulation (FIDAsim) code, which provides Monte Carlo simulations of beam neutral injection, attenuation, halo generation, halo spatial diffusion, and photoemission processes. When using the same atomic physics database, TRANSP and FIDAsim simulations achieve excellent agreement on the spatial profile and magnitude of beam and halo neutral densities and the NPA energy spectrum. The simulations show that the halo neutral density can be comparable to the beam neutral density. These halo neutrals can double the NPA flux, but they have minor effects on the NPA energy spectrum shape. The TRANSP and FIDAsim simulations also suggest that the magnitudes of beam and halo neutral densities are relatively sensitive to the choice of the atomic physics databases.

Neutral beam injection into reversed-magnetic shear DIII-D plasmas produces a variety of Alfvénic activity including toroidicity-induced Alfvén eigenmodes (TAEs) and reversed shear Alfvén eigenmodes (RSAEs). With measured equilibrium profiles as inputs, the ideal MHD code NOVA is used to calculate eigenmodes of these plasmas. The postprocessor code NOVA-K is then used to perturbatively calculate the actual stability of the modes, including finite orbit width and finite Larmor radius effects, and reasonable agreement with the spectrum of observed modes is found. Using experimentally measured mode amplitudes, fast ion orbit following simulations have been carried out in the presence of the NOVA calculated eigenmodes and are found to reproduce the dominant energy, pitch and temporal evolution of the losses measured using a large bandwidth scintillator diagnostic. The same analysis techniques applied to a DT 8 MA ITER steady-state plasma scenario with reversed-magnetic shear and both beam ion and alpha populations show Alfvén eigenmode instability. Both RSAEs and TAEs are found to be unstable with maximum growth rates occurring for toroidal mode number n = 6 and the majority of the drive coming from fast ions injected by the 1 MeV negative ion beams. AE instability due to beam ion drive is confirmed by the non-perturbative code TAEFL. Initial fast ion orbit following simulations using the unstable modes with a range of amplitudes (δ B/B = 10−5–10−3) have been carried out and show negligible fast ion loss. The lack of fast ion loss is a result of loss boundaries being limited to large radii and significantly removed from the actual modes themselves.

Neutral beam injection (NBI) is one of the primary tools foreseen for heating, current drive (CD) and q-profile control in future fusion reactors such as ITER and a Fusion Nuclear Science Facility. However, fast ions from NBI may also provide the drive for energetic particle-driven instabilities (e.g. Alfvenic modes (AEs)), which in turn redistribute fast ions in both space and energy, thus hampering the control capabilities and overall efficiency of NB-driven current. Based on experiments on the NSTX tokamak (M. Ono et al 2000 Nucl. Fusion 40 557), the effects of AEs and other low-frequency magneto-hydrodynamic instabilities on NB-CD efficiency are investigated. A new fast ion transport model, which accounts for particle transport in phase space as required for resonant AE perturbations, is utilized to obtain consistent simulations of NB-CD through the tokamak transport code TRANSP. It is found that instabilities do indeed reduce the NB-driven current density over most of the plasma radius by up to ~50%. Moreover, the details of the current profile evolution are sensitive to the specific model used to mimic the interaction between NB ions and instabilities. Implications for fast ion transport modeling in integrated tokamak simulations are briefly discussed.

Fast ion transport models currently implemented in the tokamak transport code TRANSP (Hawryluk 1980 Physics of Plasmas Close to Thermonuclear Conditions (Brussels: CEC)) are not capturing important aspects of the physics associated with resonant transport caused by instabilities such as toroidal Alfven eigenmodes (TAEs). This work describes the implementation of a fast ion transport model consistent with the basic mechanisms of resonant mode–particle interaction. The model is formulated in terms of a probability distribution function for the particle's steps in phase space, which is consistent with the Monte Carlo approach used in TRANSP. The proposed model is based on the analysis of the fast ion response to TAE modes through the ORBIT code (White and Chance 1984 Phys. Fluids 27 2455), but it can be generalized to higher frequency modes (e.g. compressional and global Alfven eigenmodes) and to other numerical codes or theories.

Energetic ion transport due to microturbulence is investigated in magnetohydrodynamic-quiescent plasmas by way of neutral beam injection in the DIII-D tokamak [J. L. Luxon, Nucl. Fusion 42, 614 (2002)]. A range of on-axis and off-axis beam injection scenarios are employed to vary relevant parameters such as the character of the background microturbulence and the value of E b / T e, where Eb is the energetic ion energy and Te the electron temperature. In all cases, it is found that any transport enhancement due to microturbulence is too small to observe experimentally. These transport effects are modeled using numerical and analytic expectations that calculate the energetic ion diffusivity due to microturbulence. It is determined that energetic ion transport due to coherent fluctuations (e.g., Alfvén eigenmodes) is a considerably larger effect and should therefore be considered more important for ITER. © 2013 AIP Publishing LLC.

Short toroidal Alfven eigenmode (TAE) avalanche bursts in the National Spherical Torus Experiment (NSTX) cause a drop in the neutron rate and sometimes a loss of neutral beam ions at or near the full injection energy over an extended range of pitch angles. The simultaneous loss of wide ranges of pitch angle suggests stochastic transport of the beam ions occurs. When beam ion orbits are followed with a guiding centre code that incorporates the plasma's magnetic equilibrium plus the measured modes, the predicted ranges of lost pitch angle are similar to those seen in the experiment, with distinct populations of trapped and passing orbits lost. These correspond to domains where the stochasticity extends in the orbit phase space from the region of beam ion deposition to the loss boundary and the trajectories along which modes may transport particles extend from the deposition volume to the loss boundary.

Fast beam-ion losses were studied in DIII-D in the presence of a scaled mock-up of two test blanket modules (TBM) for ITER. Heating of the protective tiles on the front of the TBM surface was found when neutral beams were injected and the TBM fields were engaged. The fast-ion core confinement was not significantly affected. Different orbit-following codes predict the formation of a hot spot on the TBM surface arising from beam ions deposited near the edge of the plasma. The codes are in good agreement with each other on the total power deposited at the hot spot, predicting an increase in power with decreasing separation between the plasma edge and the TBM surface. A thermal analysis of the heat flow through the tiles shows that the simulated power can account for the measured tile temperature rise. The thermal analysis, however, is very sensitive to the details of the localization of the hot spot, which is predicted to be different among the various codes. © 2011 IAEA, Vienna.

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.