MONKES: a fast neoclassical code for the evaluation of monoenergetic transport coefficients

MONKES is a new neoclassical code for the evaluation of monoenergetic transport coefficients in stellarators. By means of a convergence study and benchmarks with other codes, it is shown that MONKES is accurate and efficient. The combination of spectral discretization in spatial and velocity coordinates with block sparsity allows MONKES to compute monoenergetic coefficients at low collisionality, in a single core, in approximately one minute. MONKES is sufficiently fast to be integrated into stellarator optimization codes for direct optimization of the bootstrap current and to be included in predictive transport suites. The code and data from this paper are available at https://github.com/JavierEscoto/MONKES

Comparison of particle trajectories and collision operators for collisional transport in nonaxisymmetric plasmas

In this work, we examine the validity of several common simplifying assumptions used in numerical neoclassical calculations for nonaxisymmetric plasmas, both by using a new continuum drift-kinetic code and by considering analytic properties of the kinetic equation. First, neoclassical phenomena are computed for the LHD and W7-X stellarators using several versions of the drift-kinetic equation, including the commonly used incompressible-ExB-drift approximation and two other variants, corresponding to different effective particle trajectories. It is found that for electric fields below roughly one third of the resonant value, the different formulations give nearly identical results, demonstrating the incompressible ExB-drift approximation is quite accurate in this regime. However, near the electric field resonance, the models yield substantially different results. We also compare results for various collision operators, including the full linearized Fokker-Planck operator. At low collisionality, the radial transport driven by radial gradients is nearly identical for the different operators, while in other cases it is found to be important that collisions conserve momentum

Energetic particle loss mechanisms in reactor-scale equilibria close to quasisymmetry

Collisionless physics primarily determines the transport of fusion-born alpha particles in 3D equilibria. Several transport mechanisms have been implicated in stellarator configurations, including stochastic diffusion due to class transitions, ripple trapping, and banana drift-convective orbits. Given the guiding center dynamics in a set of six quasihelical and quasiaxisymmetric equilibria, we perform a classification of trapping states and transport mechanisms. In addition to banana drift convection and ripple transport, we observe substantial non-conservation of the parallel adiabatic invariant which can cause losses through diffusive banana tip motion. Furthermore, many lost trajectories undergo transitions between trapping classes on longer time scales, either with periodic or irregular behavior. We discuss possible optimization strategies for each of the relevant transport mechanisms. We perform a comparison between fast ion losses and metrics for the prevalence of mechanisms such as banana-drift convection [1], transitioning orbits, and wide orbit widths. Quasihelical configurations are found to have natural protection against ripple-trapping and diffusive banana tip motion leading to a reduction in prompt losses

Local and global Fokker-Planck neoclassical calculations showing flow and bootstrap current modification in a pedestal

In transport barriers, particularly H-mode edge pedestals, radial scale lengths can become comparable to the ion orbit width, causing neoclassical physics to become radially nonlocal. In this work, the resulting changes to neoclassical flow and current are examined both analytically and numerically. Steep density gradients are considered, with scale lengths comparable to the poloidal ion gyroradius, together with strong radial electric fields sufficient to electrostatically confine the ions. Attention is restricted to relatively weak ion temperature gradients (but permitting arbitrary electron temperature gradients), since in this limit a delta-f (small departures from a Maxwellian distribution) rather than full-f approach is justified. This assumption is in fact consistent with measured inter-ELM H-Mode edge pedestal density and ion temperature profiles in many present experiments, and is expected to be increasingly valid in future lower collisionality experiments. In the numerical analysis, the distribution function and Rosenbluth potentials are solved for simultaneously, allowing use of the exact field term in the linearized Fokker-Planck collision operator. In the pedestal, the parallel and poloidal flows are found to deviate strongly from the best available conventional neoclassical prediction, with large poloidal variation of a different form than in the local theory. These predicted effects may be observable experimentally. In the local limit, the Sauter bootstrap current formulae appear accurate at low collisionality, but they can overestimate the bootstrap current near the plateau regime. In the pedestal ordering, ion contributions to the bootstrap and Pfirsch-Schluter currents are also modified

A New Optimized Quasihelically SymmetricStellarator

A new optimized quasihelically symmetric configuration is described that has the desir-able properties of improved energetic particle confinement, reduced turbulent transportby 3D shaping, and non-resonant divertor capabilities. The configuration presented in thispaper is explicitly optimized for quasihelical symmetry, energetic particle confinement,neoclassical confinement, and stability near the axis. Post optimization, the configurationwas evaluated for its performance with regard to energetic particle transport, idealmagnetohydrodynamic (MHD) stability at various values of plasma pressure, and iontemperature gradient instability induced turbulent transport. The effect of discrete coilson various confinement figures of merit, including energetic particle confinement, aredetermined by generating single-filament coils for the configuration. Preliminary divertoranalysis shows that coils can be created that do not interfere with expansion of thevessel volume near the regions of outgoing heat flux, thus demonstrating the possibilityof operating a non-resonant divertor.Comment: 13 figures, 3 table

Main results of the first experimental campaign in the stellarator W7-X

A summary of the first operational phase (OP1.1) at the stellarator W7-X is given. The operational setup of heating and diagnostics as well the results of experiments are briefly described. Plasma parameters and confinement are better than expected: Te > 8 keV and Ti > 2 keV at ne ≈ 3×1019 m-3 yielding β0 ≈ 2.5 %. The results for ECR heating with X2-mode as well the ECCD are in good agreement with the theory predictions. The heating scenario with the O2-mode alone was successfully first time performed. Stellarator specific regime of core “electron root” confinement was obtained

Towards a new image processing system at Wendelstein 7-X: From spatial calibration to characterization of thermal events

Wendelstein 7-X (W7-X) is the most advanced fusion experiment in the stellarator line and is aimed at proving that the stellarator concept is suitable for a fusion reactor. One of the most important issues for fusion reactors is the monitoring of plasma facing components when exposed to very high heat loads, through the use of visible and infrared (IR) cameras. In this paper, a new image processing system for the analysis of the strike lines on the inboard limiters from the first W7-X experimental campaign is presented. This system builds a model of the IR cameras through the use of spatial calibration techniques, helping to characterize the strike lines by using the information given by real spatial coordinates of each pixel. The characterization of the strike lines is made in terms of position, size, and shape, after projecting the camera image in a 2D grid which tries to preserve the curvilinear surface distances between points. The description of the strike-line shape is made by means of the Fourier Descriptors