Numerical results for the three mono-energetic transport coefficients required for a complete neoclassical description of stellarator plasmas have been benchmarked within an international collaboration. These transport coefficients are flux-surface-averaged moments of solutions to the linearised drift kinetic equation which have been determined using field-line-integration techniques, Monte Carlo simulations, a variational method employing Fourier-Legendre test functions and a finite difference scheme. The benchmarking has been successfully carried out for past, present and future devices which represent different optimisation strategies within the extensive configuration space available to stellarators.A qualitative comparison of the results with theoretical expectations for simple model fields is provided. The behaviour of the results for the mono-energetic radial and parallel transport coefficients can be largely understood from such theoretical considerations but the mono-energetic bootstrap current coefficient exhibits characteristics which have not been predicted.

International collaboration on development of a stellarator confinement database has progressed. More than 3000 data points from nine major stellarator experiments have been compiled. Robust dependences of the energy confinement time on the density and the heating power have been confirmed. Dependences on other operational parameters, i.e. the major and minor radii, magnetic field and the rotational transform , have been evaluated using inter-machine analyses. In order to express the energy confinement in a unified scaling law, systematic differences in each subgroup are quantified. An a posteriori approach using a confinement enhancement factor on ISS95 as a renormalizing configuration-dependent parameter yields a new scaling expression ISS04; . Gyro–Bohm characteristic similar to ISS95 has been confirmed for the extended database with a wider range of plasma parameters and magnetic configurations than in the study of ISS95. It has also been discovered that there is a systematic offset of energy confinement between magnetic configurations, and its measure correlates with the effective helical ripple of the external stellarator field. Full documentation of the International Stellarator Confinement Database is available at http://iscdb.nifs.ac.jp/ and http://www.ipp.mpg.de/ISS.

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.

The neoclassical prediction of the “electron root,” i.e., a strongly positive radial electric field, Er (being the solution of the ambipolarity condition of the particle fluxes), is analyzed for low-density discharges in Wendelstein-7-AS [G. Grieger, W. Lotz, P. Merkel, et al., Phys. Fluids B 4, 2081 (1992)]. In these electron cyclotron resonance heated (ECRH) discharges with highly localized central power deposition, peaked Te profiles [with Te(0) up to 6 keV and with Ti≪Te] and strongly positive Er in the central region are measured. It is shown that this “electron root” feature at W7-AS is driven by ripple-trapped suprathermal electrons generated by the ECRH. The fraction of ripple-trapped particles in the ECRH launching plane, which can be varied at W7-AS, is found to be the most important. After switching off the heating the “electron root” feature disappears nearly immediately, i.e., two different time scales for the electron temperature decay in the central region are observed. Monte Carlo simulations in five-dimensional phase space are presented, clearly indicating that the additional “convective” electron fluxes driven by the ECRH are of the same order as the ambipolar neoclassical prediction for the “ion root” at much lower Er. For the predicted “electron root,” the ion fluxes calculated based on the traditional neoclassical ordering are much too small; shortcomings of the usual approach are indentified and a new ordering scheme is proposed.