A good understanding of the confinement of energetic ions in non-axisymmetric magnetic fields is key for the design of reactors based on the stellarator concept. In this work, we develop a model that, based on the radially-local bounce-averaged drift-kinetic equation, classifies orbits and succeeds in predicting configuration-dependent aspects of the prompt losses of energetic ions in stellarators. Such a model could in turn be employed in the optimization stage of the design of new devices.
In this paper, we analyze the impact of electron cyclotron resonance heating and electron cyclotron current drive on the Alfvénic instabilities driven by neutral beam injection observed in the TJ-II stellarator. An MHD stability analysis of driven Alfvén eigenmodes compatible with the experimental plasma parameters is carried out in order to compare with the data provided by magnetic coils, radiation monitors, and heavy ion beam probes. To this end, the vacuum magnetic configuration modified by the different levels of plasma current, the thermal plasma parameters and the fast ion pressure profiles generated by the co-injected neutral beam, are entered in the FAR3d gyro-fluid code in order to follow the linear evolution of the destabilized plasma equilibrium. Linear growth rates and radial location of the dominant predicted modes coincident in frequency with the observed fluctuations are presented. Despite the uncertainties related to the estimation of the rotational transform profile, the code predictions agree within reasonable accuracy with the experimental results.
The TJ-II stellarator neutral-beam injection (NBI) system, vacuum vessel and magnetic configuration have been included in the orbit-following Monte Carlo code ASCOT5 to simulate neutral-beam heating and current drive for high-density NBI plasmas. Co- and counter-injection beams are simulated separately. A scan in both electron density and temperature is carried out within the range of values corresponding to realistic high-density NBI plasmas, for which a low level of fast-ion losses due to charge-exchange reactions is expected, since the version of ASCOT5 used in the paper does not include such processes. The rest of the kinetic profiles (ion temperature, radial electric field and effective charge) are kept fixed. The initial distribution of markers shows that the amount of available power in the plasma carried by the beam ions depends slightly on the electron temperature and on the injection direction (co/counter). The steady-state fast-ion distribution function is obtained and used to calculate the three-dimensional fast-ion density, the neutral-beam driven current and the amount of power deposited to the plasma in the two injection scenarios. These three quantities are higher in the counter-injected case due to a lower amount of promptly lost particles. The neutral-beam current drive has been calculated using the fast-ion beam current given by ASCOT5 and the electron return current, which is computed with the analytic solution of the drift kinetic equation for electrons in the presence of fast ions in the low-collisionality regime. Neither the calculated fast-ion density nor the neutral-beam current drive are flux functions, in consistency with the fact that fast-ion drift surfaces and flux surfaces are generally not aligned.
A new quasi-isodynamic stellarator configuration optimized for the confinement of energetic ions at low plasma β is obtained. The numerical optimization is carried out using the STELLOPT suite of codes. New proxies to measure closeness to quasi-isodynamicity and quality of fast ion confinement have been included. The new configuration has poloidally closed contours of magnetic field strength, low magnetic shear and a rotational transform profile allowing an island divertor. It shows ideal and ballooning magnetohydrodynamic stability up to β = 5%, reduced effective ripple, with ϵef f < 0.5% in the plasma core. Even at low β, the configuration approximately satisfies the maximum-J property, and the confinement of fast ions is good at β ∼ 1.5% and becomes excellent at reactor values, β ∼ 4%. An evaluation of the D31 neoclassical mono-energetic coefficient supports the expectation of a reduced bootstrap current for plasmas confined in quasi-isodynamic configurations. A set of filamentary coils that preserve the good confinement of fast ions in the core is presented.
scite is a Brooklyn-based organization that helps researchers better discover and understand research articles through Smart Citations–citations that display the context of the citation and describe whether the article provides supporting or contrasting evidence. scite is used by students and researchers from around the world and is funded in part by the National Science Foundation and the National Institute on Drug Abuse of the National Institutes of Health.
hi@scite.ai
10624 S. Eastern Ave., Ste. A-614
Henderson, NV 89052, USA
Copyright © 2024 scite LLC. All rights reserved.
Made with 💙 for researchers
Part of the Research Solutions Family.