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.

An overview is given of physics differences between stellarators and tokamaks, including magnetohydrodynamic equilibrium, stability, fast-ion physics, plasma rotation, neoclassical and turbulent transport, and edge physics. Regarding microinstabilities, it is shown that the ordinary, collisionless trapped-electron mode is stable in large parts of parameter space in stellarators that have been designed so that the parallel adiabatic invariant decreases with radius. Also, the first global, electromagnetic, gyrokinetic stability calculations done for Wendelstein 7-X suggest that kinetic ballooning modes are more stable than in a typical tokamak.

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.

Using a global approach for solving an ion gyrokinetic model in three-dimensional geometry the linear stability and structure of ion-temperature-gradient (ITG) modes in the configuration of the stellarator Wendelstein 7-X (W7-X) [G. Grieger et al., in Plasma Physics and Controlled Nuclear Fusion Research 1990 (International Atomic Energy Agency, Vienna, 1991), Vol. 3, p. 525.] is studied. The time evolution of electrostatic perturbations is solved as an initial value problem with a particle-in-cell δf method. The vacuum magnetohydrodynamic equilibrium is calculated by the code VMEC [S. P. Hirshman and D. K. Lee, Comput. Phys. Commun. 39, 161 (1986)]. In this work the most unstable ITG mode in W7-X is presented. This mode has a pronounced ballooning-type structure; however, it is not tokamak-like. A driving mechanism analysis using the energy transfer shows that the contribution of curvature effects is non-negligible. The growth rate and the mixing-length estimate for transport are compared with those for ITG modes found in axisymmetric geometries.

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.

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.

The particle transport of impurities in magnetically confined plasmas under some conditions does not find, neither quantitatively nor qualitatively, a satisfactory theory-based explanation. This compromise the successful realization of thermo-nuclear fusion for energy production since its accumulation is known to be one of the causes that leads to the plasma breakdown. In standard reactor-relevant conditions this accumulation is in most stellarators intrinsic to the lack of toroidal symmetry, that leads to the neoclassical electric field to point radially inwards. This statement, that the standard theory allows to formulate, has been contradicted by some experiments that showed weaker or no accumulation under such conditions [1,2]. The charge state of the impurities makes its transport more sensitive to the electric fields. Thus, the short length scale turbulent electrostatic potential or its long wavelength variation on the flux surface Φ 1 -that the standard neoclassical approach usually neglects -might possibly shed some light on the experimental findings. In the present work the focus is put on the second of the two, and investigate its influence of the radial transport of C 6+ . We show that in LHD it is strongly modified by Φ 1 , both resulting in mitigated/enhanced accumulation at internal/external radial positions; for Wendelstein 7-X, on the contrary, Φ 1 is expected to be considerably smaller and the transport of C 6+ not affected up to an appreciable extent; and in TJ-II the potential shows a moderate impact despite of the large amplitude of Φ 1 for the parameters considered.

Aiming for filling a corresponding lack for sophisticated test cases for global electromagnetic gyrokinetic codes, a new hierarchical benchmark is proposed. Starting from established test sets with adiabatic electrons, fully gyrokinetic electrons and electrostatic fluctuations are taken into account before finally studying global electromagnetic micro-instabilities. Results from up to five codes involving representatives from as different numerical approaches as particle-in-cell methods, Eulerian and Semi-Lagrange are shown. By means of spectrally resolved growth rates and frequencies and mode structure comparisons, agreement can be confirmed on ion-gyro-radius scales thus providing confidence in the correct implementation of the underlying equations.