Demonstrating improved confinement of energetic ions is one of the key goals of the Wendelstein 7-X (W7-X) stellarator. In the past campaigns, measuring confined fast ions has proven to be challenging. Future deuterium campaigns would open up the option of using fusion-produced neutrons to indirectly observe confined fast ions. There are two neutron populations: 2.45 MeV neutrons from thermonuclear and beam-target fusion, and 14.1 MeV neutrons from DT reactions between tritium fusion products and bulk deuterium. The 14.1 MeV neutron signal can be measured using a scintillating fiber neutron detector, whereas the overall neutron rate is monitored by common radiation safety detectors, for instance fission chambers. The fusion rates are dependent on the slowing-down distribution of the deuterium and tritium ions, which in turn depend on the magnetic configuration via fast ion orbits. In this work, we investigate the effect of magnetic configuration on neutron production rates in W7-X. The neutral beam injection, beam and triton slowing-down distributions, and the fusion reactivity are simulated with the ASCOT suite of codes. The results indicate that the magnetic configuration has only a small effect on the production of 2.45 MeV neutrons from DD fusion and, particularly, on the 14.1 MeV neutron production rates. Despite triton losses of up to 50 %, the amount of 14.1 MeV neutrons produced might be sufficient for a time-resolved detection using a scintillating fiber detector, although only in high-performance discharges.

This paper presents the current state of the global gyrokinetic code Orb5 as an update of the previous reference [Jolliet et al., Comp. Phys. Commun. 177 409 (2007)]. The Orb5 code solves the electromagnetic Vlasov-Maxwell system of equations using a PIC scheme and also includes collisions and strong flows. The code assumes multiple gyrokinetic ion species at all wavelengths for the polarization density and drift-kinetic electrons. Variants of the physical model can be selected for electrons such as assuming an adiabatic response or a "hybrid" model in which passing electrons are assumed adiabatic and trapped electrons are drift-kinetic. A Fourier filter as well as various control variates and noise reduction techniques enable simulations with good signal-to-noise ratios at a limited numerical cost. They are completed with different momentum and zonal flow-conserving heat sources allowing for temperature-gradient and flux-driven simulations. The code, which runs on both CPUs and GPUs, is well benchmarked against other similar codes and analytical predictions, and shows good scalability up to thousands of nodes.

The collisionless time evolution of zonal flows in stellarator systems is investigated. An analytical solution of the kinetic and quasineutrality equations describing the residual zonal flow is derived for arbitrary three-dimensional systems without approximations in the magnetic geometry. The theory allows for an arbitrary number of particle species. It has been found that in stellarators the residual zonal flows are not in general steady but oscillate with a certain frequency. This frequency is determined by the speed of the bounce-averaged radial drifts of the particles trapped in the magnetic field and vanishes in tokamaks, where such net drifts are absent. A reduction of the bounce-averaged radial drifts in configurations optimized with respect to neoclassical transport results in a smaller zonal-flow frequency. * alexey.mishchenko@ipp.mpg.de

The linear response of a collisionless stellarator plasma to an applied radial electric field is calculated, both analytically and numerically. Unlike in a tokamak, the electric field and associated zonal flow develop oscillations before settling down to a stationary state, the so-called Rosenbluth-Hinton flow residual. These oscillations are caused by locally trapped particles with radially drifting bounce orbits. These particles also cause a kind of Landau damping of the oscillations that depends on the magnetic configuration. The relative importance of geodesic acoustic modes and zonal-flow oscillations therefore varies among different stellarators.

No abstract