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.

After completing the main construction phase of Wendelstein 7-X (W7-X) and successfully commissioning the device, first plasma operation started at the end of 2015. Integral commissioning of plasma start-up and operation using electron cyclotron resonance heating (ECRH) and an extensive set of plasma diagnostics have been completed, allowing initial physics studies during the first operational campaign. Both in helium and hydrogen, plasma breakdown was easily achieved. Gaining experience with plasma vessel conditioning, discharge lengths could be extended gradually. Eventually, discharges lasted up to 6 s, reaching an injected energy of 4 MJ, which is twice the limit originally agreed for the limiter configuration employed during the first operational campaign. At power levels of 4 MW central electron densities reached 3  ×  1019 m−3, central electron temperatures reached values of 7 keV and ion temperatures reached just above 2 keV. Important physics studies during this first operational phase include a first assessment of power balance and energy confinement, ECRH power deposition experiments, 2nd harmonic O-mode ECRH using multi-pass absorption, and current drive experiments using electron cyclotron current drive. As in many plasma discharges the electron temperature exceeds the ion temperature significantly, these plasmas are governed by core electron root confinement showing a strong positive electric field in the plasma centre.

In a plasma with a population of super-thermal particles generated by heating or fusion processes, kinetic effects can lead to the additional destabilisation of MHD modes or even to additional energetic particle modes. In order to describe these modes, a new linear gyrokinetic MHD code has been developed and tested, LIGKA (linear gyrokinetic shear Alfvén physics) [Ph. Lauber, Linear gyrokinetic description of fast particle effects on the MHD stability in tokamaks, Ph.D. Thesis, TU Mü nchen, 2003; Ph. Lauber, S. Gü nter, S.D. Pinches, Phys. Plasmas 12 (2005) 122501], based on a gyrokinetic model [H. Qin, Gyrokinetic theory and computational methods for electromagnetic perturbations in tokamaks, Ph.D.Thesis, Princeton University, 1998]. A finite Larmor radius expansion together with the construction of some fluid moments and specification to the shear Alfvén regime results in a self-consistent, electromagnetic, non-perturbative model, that allows not only for growing or damped eigenvalues but also for a change in mode-structure of the magnetic perturbation due to the energetic particles and background kinetic effects.Compared to previous implementations [H. Qin, mentioned above], this model is coded in a more general and comprehensive way. LIGKA uses a Fourier decomposition in the poloidal coordinate and a finite element discretisation in the radial direction. Both analytical and numerical equilibria can be treated. Integration over the unperturbed particle orbits is performed with the drift-kinetic HAGIS code [S.D. Pinches, Ph.D. Thesis, The University of Nottingham, 1996; S.D. Pinches et al., CPC 111 (1998) 131] which accurately describes the particles' trajectories. This allows finite-banana-width effects to be implemented in a rigorous way since the linear formulation of the model allows the exchange of the unperturbed orbit integration and the discretisation of the perturbed potentials in the radial direction.Successful benchmarks for toroidal Alfvén eigenmodes (TAEs) and kinetic Alfvén waves (KAWs) with analytical results, ideal MHD codes, drift-kinetic codes and other codes based on kinetic models are reported.

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

No abstract

A global nonlinear simulation code for the time evolution of ion-temperature-gradient-driven modes in θ-pinch geometry as a first approximation to the stellarator Wendelstein 7-X (W7-X) [Grieger et al., Proceedings of the 13th International Conference on Plasma Physics and Controlled Nuclear Fusion Research, Washington, DC, 1990 (International Atomic Energy Agency, Vienna, 1991), Vol. 3, p. 525] has been developed. A δf particle-in-cell (PIC) method is used to solve the coupled system of gyrokinetic equations for the ions, in the electrostatic approximation, and the quasineutrality equation, assuming adiabatically responding electrons. The focus has been on adherence to conservation laws, i.e., particle number and energy conservation. Besides other improvements it has been shown that a well-chosen initial distribution of the markers in reduced phase space makes optimal use of the δf PIC method to reduce the statistical noise for a given number of markers. In a model including all (1351) physically relevant modes, it has been possible to achieve energy conservation beyond the saturation of the instability.

This paper reports self-consistent global linear gyrokinetic particle-in-cell simulations of shear Alfvén waves destabilized by fast particles in tokamak geometry. Resonant excitation of Toroidal Alfvén Eigenmodes by fast particles and their transition to Energetic Particle Modes (when the fast-particle drive is large enough) has been observed in the simulations.

In the linear collisionless limit, a zonal potential perturbation in a toroidal plasma relaxes, in general, to a non-zero residual value. Expressions for the residual value in tokamak and stellarator geometries, and for arbitrary wavelengths, are derived. These expressions involve averages over the lowest order particle trajectories, that typically cannot be evaluated analytically. In this work, an efficient numerical method for the evaluation of such expressions is reported. It is shown that this method is faster than direct gyrokinetic simulations performed with the Gene and EUTERPE codes. Calculations of the residual value in stellarators are provided for much shorter wavelengths than previously available in the literature. Electrons must be treated kinetically in stellarators because, unlike in tokamaks, kinetic electrons modify the residual value even at long wavelengths. This effect, that had already been predicted theoretically, is confirmed by gyrokinetic simulations.