Search

In omnigeneous magnetic fields, charged particles are perfectly confined in the absence of collisions and turbulence. For this reason, the magnetic configuration is optimized to be close to omnigenity in any candidate for a stellarator fusion reactor. However, approaching omnigenity imposes severe constraints on the spatial variation of the magnetic field. In particular, the topology of the contours of constant magnetic-field-strength on each magnetic surface must be such that there are no particles transitioning between different types of wells. This, in turn, usually leads to complicated plasma shapes and coils. This Letter presents a new family of optimized fields that display tokamak-like collisional energy transport while having transitioning particles. This result radically broadens the space of accessible reactor-relevant configurations.

CIEMAT-QI4 is a quasi-isodynamic stellarator configuration that simultaneously features very good fast-ion confinement in a broad range of $\beta$ values, low neoclassical transport and bootstrap current, and ideal magnetohydrodynamic stability up to $\beta=5\%$. In this paper it is shown that CIEMAT-QI4 also exhibits reduced turbulent transport. This is demonstrated through nonlinear electrostatic simulations with the gyrokinetic code stella, including kinetic ions and electrons. The relation between reduced turbulent transport and the fact that CIEMAT-QI4 very approximately satisfies the so-called maximum-J property is discussed., Comment: 12 pages, 9 figures

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/

Stellarator magnetic configurations need to be optimized in order to meet all the required properties of a fusion reactor. In this work, it is shown that a flat-mirror quasi-isodynamic configuration (i.e. a quasi-isodynamic configuration with sufficiently small radial variation of the mirror term) can achieve small radial transport of energy and good confinement of bulk and fast ions even if it is not very close to perfect omnigeneity, and for a wide range of plasma scenarios, including low $\beta$ and small radial electric field. This opens the door to constructing better stellarator reactors. On the one hand, they would be easier to design, as they would be robust against error fields. On the other hand, they would be easier to operate since, both during startup and steady-state operation, they would require less auxiliary power, and the damage to plasma-facing components caused by fast ion losses would be reduced to acceptable levels.

A new quasi-isodynamic stellarator configuration optimized for the confinement of energetic ions at low plasma $\beta$ is obtained. The numerical optimization is carried out using the STELLOPT suite of codes. New proxies to measure closeness to quasi-isodynamicity and quality of fast ion confinement have been included. The new configuration has poloidally closed contours of magnetic field strength, low magnetic shear and a rotational transform profile allowing an island divertor. It shows ideal and ballooning magnetohydrodynamic stability up to $\beta = 5%$, reduced effective ripple, with $\epsilon_{eff} < 0.5%$ in the plasma core. Even at low $\beta$, the configuration approximately satisfies the maximum-$J$ property, and the confinement of fast ions is good at $\beta \sim 1.5%$ and becomes excellent at reactor values, $\beta \sim 4%$. An evaluation of the $D_{31}$ neoclassical mono-energetic coefficient supports the expectation of a reduced bootstrap current for plasmas confined in quasi-isodynamic configurations. A set of filamentary coils that preserve the good confinement of fast ions in the core is presented., Comment: 21 pages, 16 figures

In the present Letter, first-of-its-kind computer simulations predicting plasma profiles for modern optimized stellarators -- while self-consistently retaining neoclassical transport, turbulent transport with 3D effects, and external physical sources -- are presented. These simulations exploit a newly developed coupling framework involving the global gyrokinetic turbulence code GENE-3D, the neoclassical transport code KNOSOS, and the 1D transport solver TANGO. This framework is used to analyze the recently observed degradation of energy confinement in electron-heated plasmas in the Wendelstein 7-X stellarator, where the central ion temperature was "clamped" to $T_i \approx 1.5$ keV regardless of the external heating power. By performing first-principles based simulations, we provide key evidence to understand this effect, namely the inefficient thermal coupling between electrons and ions in a turbulence-dominated regime, which is exacerbated by the large $T_e/T_i$ ratios, and show that a more efficient ion heat source, such as direct ion heating, will increase the on-axis ion temperature. This work paves the way towards the use of high-fidelity models for the development of the next generation of stellarators, in which neoclassical and turbulent transport are optimized simultaneously.

In reactor-relevant plasmas, neoclassical transport drives an outward particle flux in the core of large stellarators and predicts strongly hollow density profiles. However, this theoretical prediction is contradicted by experiments. In particular, in Wendelstein 7-X, the first large optimized stellarator, flat or weakly peaked density profiles are generally measured, indicating that neoclassical theory is not sufficient and that an inward contribution to the particle flux is missing in the core. In this Research Letter, it is shown that the turbulent contribution to the particle flux can explain the difference between experimental measurements and neoclassical predictions. The results of this Research Letter also prove that theoretical and numerical tools are approaching the level of maturity needed for the prediction of equilibrium density profiles in stellarator plasmas, which is a fundamental requirement for the design of operation scenarios of present devices and future reactors., Comment: 8 pages, 6 figures, updated the title and incorporated the appendix into the main letter

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.

A recent characterization of core turbulence carried out with a Doppler reflectometer in the optimized stellarator Wendelstein 7-X (W7-X) found that discharges achieving high ion temperatures at the core featured an ITG-like suppression of density fluctuations driven by a reduction of the gradient ratio $\eta_i = L_n/L_{T_i}$ [D. Carralero et al., Nucl. Fusion, 2021]. In order to confirm the role of ITG turbulence in this process, we set out to establish experimentally the relation between core density fluctuations, turbulent heat flux and global confinement. With this aim, we consider the scenarios found in the previous work and carry out power balance analysis for a number of representative ones, including some featuring high ion temperature. As well, we evaluate the global energy confinement time and discuss it in the context of the ISS04 inter-stellarator scaling. We find that, when turbulence is suppressed as a result of a reduction of $\eta_i$, there is a reduction of ion turbulent transport, and global performance is improved as a result. This is consistent with ITG turbulence limiting the ion temperature at the core of W7-X. In contrast, when turbulence is reduced following a decrease in collisionality, no changes are observed in transport or confinement. This could be explained by ITG modes being combined with TEM turbulence when the later is destabilized at low collisionalities., Comment: Submitted to PPCF, 47th EPS special issue

The ongoing development of electromagnets based on High Temperature Superconductors has led to the conceptual exploration of high-magnetic-field fusion reactors of the tokamak type, operating at on-axis fields above 10 T. In this work we explore the consequences of the potential future availability of high-field three-dimensional electromagnets on the physics design point of a stellarator reactor. We find that, when an increase in the magnetic field strength $B$ is used to maximally reduce the device linear size $R\sim B^{-4/3}$ (with otherwise fixed magnetic geometry), the physics design point is largely independent of the chosen field strength/device size. A similar degree of optimization is to be imposed on the magnetohydrodynamic, transport and fast ion confinement properties of the magnetic configuration of that family of reactor design points. Additionally, we show that the family shares an invariant operation map of fusion power output as a function of the auxiliary power and relative density variation. The effects of magnetic field over-engineering and the $R(B)$ scaling of design points with constant neutron wall loading are also inspected. In this study we use geometric parameters characteristic of the helias reactor, but most results apply to other stellarator configurations

A good understanding of the confinement of energetic ions in non-axisymmetric magnetic fields is key for the design of reactors based on the stellarator concept. In this work, we develop a model that, based on the radially-local bounce-averaged drift-kinetic equation, classifies orbits and succeeds in predicting configuration-dependent aspects of the prompt losses of energetic ions in stellarators. Such a model could in turn be employed in the optimization stage of the design of new devices.

In this work, a new version of KNOSOS is presented. KNOSOS is a low-collisionality radially-local, bounce-averaged neoclassical code that is extremely fast, and at the same time, includes physical effects often neglected by more standard codes: the component of the magnetic drift that is tangent to the flux-surface and the variation of the electrostatic potential on the flux-surface. An earlier version of the code could only describe configurations that were sufficiently optimized with respect to neoclassical transport. KNOSOS can now be applied to any large aspect ratio stellarator, and its performance is demonstrated by means of detailed simulations in the configuration space of Wendelstein 7-X., Comment: IAEA paper, submitted to Nuclear Fusion

First results from the optimized helias Wendelstein 7-X stellarator (W7-X) have shown that core transport is no longer mostly neoclassical, as is the case in previous kinds of stellarators. Instead, turbulent transport poses a serious limitation to the global performance of the machine. Several studies have found this particularly relevant for ion transport, with core ion temperatures becoming clamped at relatively low values of $T_{i} \simeq 1.7$ keV, except in the few scenarios in which turbulence can be suppressed. In order to understand turbulent mechanisms at play, it is important to have a clear understanding of the parametric dependencies of turbulent fluctuations, and the relation between them and turbulent transport. In this work we use Doppler reflectometry measurements carried out during a number of relevant operational scenarios to provide a systematic characterization of ion-scale ($k_\perp\rho_i\simeq 1$) density fluctuations in the core of W7-X. Then, we study the relation between fluctuation amplitude and plasma profiles and show how distinct regimes can be defined for the former, depending on normalized gradients $a/L_{ne}$ and $a/L_{Ti}$. Furthermore, we discuss the importance of other potentially relevant parameters such as $T_e/T_i$, $E_r$ or collisionality. Comparing the different regimes, we find that turbulence amplitude depends generally on the gradient ratio $\eta_i=L_{ne}/L_{Ti}$, as would be expected for ITG modes, with the exception of a range of discharges, for which turbulence suppression may be better explained by an ITG to TEM transition triggered by a drop in collisionality. Finally, we show a number of scenarios under which $T_{i,core} > 1.7$ keV is achieved and how core fluctuations are suppressed in all of them, thus providing experimental evidence of microturbulence being the main responsible for the limited ion confinement in W7-X., Comment: Manuscript submitted to Nuclear Fusion in May 2021

A study of turbulent impurity transport by means of quasilinear and nonlinear gyrokinetic simulations is presented for Wendelstein 7-X (W7-X). The calculations have been carried out with the recently developed gyrokinetic code stella. Different impurity species are considered in the presence of various types of background instabilities: ITG, TEM and ETG modes for the quasilinear part of the work; ITG and TEM for the nonlinear results. While the quasilinear approach allows one to draw qualitative conclusions about the sign or relative importance of the various contributions to the flux, the nonlinear simulations quantitatively determine the size of the turbulent flux and check the extent to which the quasilinear conclusions hold. Although the bulk of the nonlinear simulations are performed at trace impurity concentration, nonlinear simulations are also carried out at realistic effective charge values, in order to know to what degree the conclusions based on the simulations performed for trace impurities can be extrapolated to realistic impurity concentrations. The presented results conclude that the turbulent radial impurity transport in W7-X is mainly dominated by ordinary diffusion, which is close to that measured during the recent W7-X experimental campaigns. It is also confirmed that thermo-diffusion adds a weak inward flux contribution and that, in the absence of impurity temperature and density gradients, ITG- and TEM-driven turbulence push the impurities inwards and outwards, respectively., Comment: 19 pages, 10 figures, 2 tables

In this work we present the first measurements obtained by the V-band Doppler reflectometer during the second operation phase of Wendelstein 7-X to discuss the influence in the velocity shear layer and the radial electric field, E$_r$, of several plasma parameters such as magnetic configuration, rotational transform or degree of detachment. In the first place, we carry out a systematic characterization of the turbulence rotation velocity profile in order to describe the influence of density and heating power on E$_r$ under the four most frequent magnetic configurations. The $|$E$_r|$ value in the edge is found to increase with configurations featuring higher $\iota$, although this does not apply for the high mirror configuration, KJM. As well, the E$_r$ value in the SOL and the velocity shear near the separatrix are found to display a clear dependence on heating power and density for all configurations. For a number of relevant cases, these results are assessed by comparing them to neoclassical predictions obtained from the codes DKES and KNOSOS, finding generally good agreement with experimental results. Finally, the evolution of E$_r$ at the edge is evaluated throughout the island-divertor detachment regime achieved for the first time in the 2018 campaign. After detachment, $|$E$_r|$ is reduced both at the SOL and edge, and the plasma column shrinks, with the shear layer seemingly moving radially inwards from the separatrix., Comment: Sent for publication to Nuclear Fusion

KNOSOS (KiNetic Orbit-averaging SOlver for Stellarators) is a freely available, open-source code (\href{https://github.com/joseluisvelasco/KNOSOS}{https://github.com/joseluisvelasco/KNOSOS}) that calculates neoclassical transport in low-collisionality plasmas of three-dimensional magnetic confinement devices by solving the radially local drift-kinetic and quasineutrality equations. The main feature of KNOSOS is that it relies on orbit-averaging to solve the drift-kinetic equation very fast. KNOSOS treats rigorously the effect of the component of the magnetic drift that is tangent to magnetic surfaces, and of the component of the electrostatic potential that varies on the flux surface, {\varphi}_1. Furthermore, the equation solved is linear in {\varphi}_1, which permits an efficient solution of the quasineutrality equation. As long as the radially local approach is valid, KNOSOS can be applied to the calculation of neoclassical transport in stellarators (helias, heliotrons, heliacs, etc.) and tokamaks with broken axisymmetry. In this paper, we show several calculations for the stellarators W7-X, LHD, NCSX and TJ-II that provide benchmark with standard local codes and demonstrate the advantages of this approach.

Main ions influence impurity dynamics through a variety of mechanisms; in particular, via impurity-ion collisions. To lowest order in an expansion in the main ion mass over the impurity mass, the impurity-ion collision operator only depends on the component of the main ion distribution that is odd in the parallel velocity. These lowest order terms give the parallel friction of the impurities with the main ions, which is typically assumed to be the main cause of collisional impurity transport. Next-order terms in the mass ratio expansion of the impurity-ion collision operator, proportional to the component of the main ion distribution that is even in the parallel velocity, are usually neglected. However, in stellarators, the even component of the main ion distribution can be very large. In this article, such next-order terms in the mass ratio expansion of the impurity-ion collision operator are retained, and analytical expressions for the neoclassical radial flux of trace impurities are calculated in the Pfirsch-Schl\"uter, plateau and $1/\nu$ regimes. The new terms provide a drive for impurity transport that is physically very different from parallel friction: they are associated to anisotropy in the pressure of the main ions, which translates into impurity pressure anisotropy. It is argued that main ion pressure anisotropy must be taken into account for a correct description of impurity transport in certain realistic stellarator plasmas. Examples are given by numerically evaluating the analytical expressions for the impurity flux., Comment: 40 pages, 4 figures. Published version

The component of the neoclassical electrostatic potential that is non-constant on the magnetic surface, that we denote by $\tilde\varphi$, can affect radial transport of highly charged impurities, and this has motivated its inclusion in some modern neoclassical codes. The number of neoclassical simulations in which $\tilde\varphi$ is calculated is still scarce, partly because they are usually demanding in terms of computational resources, especially at low collisionality. In this paper the size, the scaling with collisionality and with aspect ratio, and the structure of $\tilde\varphi$ on the magnetic surface are analytically derived in the $1/\nu$, $\sqrt{\nu}$ and superbanana-plateau regimes of stellarators close to omnigeneity; i. e. stellarators that have been optimized for neoclassical transport. It is found that the largest $\tilde\varphi$ that the neoclassical equations admit scales linearly with the inverse aspect ratio and with the size of the deviation from omnigeneity. Using a model for a perturbed omnigeneous configuration, the analytical results are verified and illustrated with calculations by the code KNOSOS. The techniques, results and numerical tools employed in this paper can be applied to neoclassical transport problems in tokamaks with broken axisymmetry., Comment: 30 pages, 12 figures, 1 table. Published version

In the present work we report recent radial electric field measurements carried out with the Doppler reflectometry system in the TJ-II stellarator. The study focuses on the fact that, under some conditions, the radial electric field measured at different points over the same flux surface shows significantly different values. A numerical analysis is carried out considering the contribution arising from the radial dependence of $\Phi_1$ as a possible correction term to the total radial electric field. Here $\Phi_1$ is the neoclassical electrostatic potential variation over the surface. The comparison shows good agreement in some aspects, like the conditions under which this correction is large (electron-root conditions) or negligible (ion-root conditions). But it disagrees in others like the sign of the correction. The results are discussed together with the underlying reasons of this partial disagreement. In addition, motivated by the recent installation of the dual Doppler reflectometry system in Wendelstein 7-X (W7-X), $\Phi_1$ estimations for W7-X are revisited considering Core-Electron-Root-Plasma (CERC) plasmas from its first experimental campaign. The simulations show larger values of $\Phi_1$ under electron-root conditions than under ion root ones. The contribution from the kinetic electron response is shown to become important at some radii. All this results in a potential variation size noticeably larger than estimated in our previous work in W7-X \cite{Regana_nf_57_056004_2017} for other plasma parameters and another configuration., Comment: 23 pages, 10 figures

The control of impurity accumulation is one of the main challenges for future stellarator fusion reactors. The standard argument to explain this accumulation relies on the, in principle, large inward pinch in the neoclassical impurity flux caused by the typically negative radial electric field in stellarators. This simplified interpretation was proven to be flawed by Helander et al. [Phys. Rev. Lett. 118, 155002 (2017)], who showed that in a relevant regime (low-collisionality main ions and collisional impurities) the radial electric field does not drive impurity transport. In that reference, the effect of the component of the electric field that is tangent to the magnetic surface was not included. In this Letter, an analytical calculation of the neoclassical radial impurity flux incorporating such effect is given, showing that it can be very strong for highly charged impurities and that, once it is taken into account, the dependence of the impurity flux on the radial electric field reappears. Realistic examples are provided in which the inclusion of the tangential electric field leads to impurity expulsion., Comment: 10 pages, 1 figure. Published version

The low frequency oscillatory relaxation of zonal potential perturbations is studied numerically in the TJ-II stellarator (where it was experimentally detected for the first time). It is studied in full global gyrokinetic simulations of multi-species plasmas. The oscillation frequency obtained is compared with predictions based on single-species simulations using simplified analytical relations. It is shown that the frequency of this oscillation for a multi-species plasma can be accurately obtained from single-species calculations using extrapolation formulas. The damping of the oscillation and the influence of the different inter-species collisions is studied in detail. It is concluded that taking into account multiple kinetic ions and electrons with impurity concentrations realistic for TJ-II plasmas allows to account for the values of frequency and damping rate in zonal flows relaxations observed experimentally., Comment: 11 figures, 22 pages

Low-collisionality stellarator plasmas usually display a large negative radial electric field that has been expected to cause accumulation of impurities due to their high charge number. In this paper, two combined effects that can potentially modify this scenario are discussed. First, it is shown that, in low collisionality plasmas, the kinetic contribution of the electrons to the radial electric field can make it negative but small, bringing the plasma close to impurity temperature screening (i.e., to a situation in which the ion temperature gradient is the main drive of impurity transport and causes outward flux); in plasmas of very low collisionality, such as those of the Large Helical Device displaying impurity hole, screening may actually occur. Second, the component of the electric field that is tangent to the flux surface (in other words, the variation of the electrostatic potential on the flux surface), although smaller than the radial component, has recently been suggested to be an additional relevant drive for radial impurity transport. Here, it is explained that, especially when the radial electric field is small, the tangential magnetic drift has to be kept in order to correctly compute the tangential electric field, that can be larger than previously expected. This can have a strong impact on impurity transport, as we illustrate by means of simulations using the newly-developed code KNOSOS (KiNetic Orbit-averaging-SOlver for Stellarators).

In general, the orbit-averaged radial magnetic drift of trapped particles in stellarators is non-zero due to the three-dimensional nature of the magnetic field. Stellarators in which the orbit-averaged radial magnetic drift vanishes are called omnigeneous, and they exhibit neoclassical transport levels comparable to those of axisymmetric tokamaks. However, the effect of deviations from omnigeneity cannot be neglected in practice. For sufficiently low collision frequencies (below the values that define the $1/\nu$ regime), the components of the drifts tangential to the flux surface become relevant. This article focuses on the study of such collisionality regimes in stellarators close to omnigeneity when the gradient of the non-omnigeneous perturbation is small. First, it is proven that closeness to omnigeneity is required to preserve radial locality in the drift-kinetic equation for collisionalities below the $1/\nu$ regime. Then, it is shown that neoclassical transport is determined by two layers in phase space. One of the layers corresponds to the $\sqrt{\nu}$ regime and the other to the superbanana-plateau regime. The importance of the superbanana-plateau layer for the calculation of the tangential electric field is emphasized, as well as the relevance of the latter for neoclassical transport in the collisionality regimes considered in this paper. In particular, the tangential electric field is essential for the emergence of a new subregime of superbanana-plateau transport when the radial electric field is small. A formula for the ion energy flux that includes the $\sqrt{\nu}$ regime and the superbanana-plateau regime is given. The energy flux scales with the square of the size of the deviation from omnigeneity. Finally, it is explained why below a certain collisionality value the formulation presented in this article ceases to be valid., Comment: 36 pages. Version to be published in Plasma Physics and Controlled Fusion

Measurements of the relaxation of a zonal electrostatic potential perturbation in a non-axisymmetric magnetically confined plasma are presented. A sudden perturbation of the plasma equilibrium is induced by the injection of a cryogenic hydrogen pellet in the TJ-II stellarator, which is observed to be followed by a damped oscillation in the electrostatic potential. The waveform of the relaxation is consistent with theoretical calculations of zonal potential relaxation in a non-axisymmetric magnetic geometry. The turbulent transport properties of a magnetic confinement configuration are expected to depend on the features of the collisionless damping of zonal flows, of which the present letter is the first direct observation.

Achieving impurity and helium ash control is a crucial issue in the path towards fusion-grade magnetic confinement devices, and this is particularly the case of helical reactors, whose low-collisionality ion-root operation scenarios usually display a negative radial electric field which is expected to cause inwards impurity pinch. In these work we discuss, based on experimental measurements and standard predictions of neoclassical theory, how plasmas of very low ion collisionality, similar to those observed in the impurity hole of the Large Helical Device, can be an exception to this general rule, and how a negative radial electric field can coexist with an outward impurity flux. This interpretation is supported by comparison with documented discharges available in the International Stellarator-Heliotron Profile Database, and it can be extrapolated to show that achievement of high ion temperature in the core of helical devices is not fundamentally incompatible with low core impurity content.

We review in a tutorial fashion some of the causes of impurity density variations along field lines and radial impurity transport in the moment approach framework. An explicit and compact form of the parallel inertia force valid for arbitrary toroidal geometry and magnetic coordinates is derived and shown to be non-negligible for typical TJ-II plasma conditions. In the second part of the article, we apply the fluid model including main ion-impurity friction and inertia to observations of asymmetric emissivity patterns in neutral beam heated plasmas of the TJ-II stellarator. The model is able to explain qualitatively several features of the radiation asymmetry, both in stationary and transient conditions, based on the calculated in-surface variations of the impurity density.

We study radial particle transport in stellarator plasmas using cryogenic pellet injection. By means of perturbative experiments, we estimate the experimental particle flux and compare it with neoclassical simulations. Experimental evidence is obtained of the fact that core depletion in helical devices can be slowed-down even by pellets that do not reach the core region. This phenomenon is well captured by neoclassical predictions with DKES and FORTEC-3D.

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 \cite{Ida_pop_16_056111_2009, Yoshinuma_nf_49_062002_2009}. 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 wave-length variation on the flux surface $\Phi_{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 $\Phi_{1}$, both resulting in mitigated/enhanced accumulation at internal/external radial positions; for Wendelstein 7-X, on the contrary, $\Phi_{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 $\Phi_1$ for the parameters considered.

Quasisymmetric stellarators are a type of optimized stellarators for which flows are undamped to lowest order in an expansion in the normalized Larmor radius. However, perfect quasisymmetry is impossible. Since large flows may be desirable as a means to reduce turbulent transport, it is important to know when a stellarator can be considered to be sufficiently close to quasisymmetry. The answer to this question depends strongly on the size of the spatial gradients of the deviation from quasisymmetry and on the collisionality regime. Recently, formal criteria for closeness to quasisymmetry have been derived in a variety of situations. In particular, the case of deviations with large gradients was solved in the $1/\nu$ regime. Denoting by $\alpha$ a parameter that gives the size of the deviation from quasisymmetry, it was proven that particle fluxes do not scale with $\alpha^{3/2}$, as typically claimed, but with $\alpha$. It was also shown that ripple wells are not necessarily the main cause of transport. This paper reviews those works and presents a new result in another collisionality regime, in which particles trapped in ripple wells are collisional and the rest are collisionless., Comment: 14 pages, 2 figures. To appear in Plasma Physics and Controlled Fusion

Plasma flow is damped in stellarators because they are not intrinsically ambipolar, unlike tokamaks, in which the flux-surface averaged radial electric current vanishes for any value of the radial electric field. Only quasisymmetric stellarators are intrinsically ambipolar, but exact quasisymmetry is impossible to achieve in non-axisymmetric toroidal configurations. By calculating the violation of intrinsic ambipolarity due to deviations from quasisymmetry, one can derive criteria to assess when a stellarator can be considered quasisymmetric in practice, i.e. when the flow damping is weak enough. Let us denote by $\alpha$ a small parameter that controls the size of a perturbation to an exactly quasisymmetric magnetic field. Recently, it has been shown that if the gradient of the perturbation is sufficiently small, the flux-surface averaged radial electric current scales as $\alpha^2$ for any value of the collisionality. It was also argued that when the gradient of the perturbation is large, the quadratic scaling is replaced by a more unfavorable one. In this paper, perturbations with large gradients are rigorously treated. In particular, it is proven that for low collisionality a perturbation with large gradient yields, at best, an $O(|\alpha|)$ deviation from quasisymmetry. Heuristic estimations in the literature incorrectly predicted an $O(|\alpha|^{3/2})$ deviation., Comment: 24 pages, 2 figures. To appear in Plasma Physics and Controlled Fusion

First observations of electrostatic potential variations within the flux surfaces of a toroidal magnetic confinement device are presented. Measurements are taken in the TJ-II stellarator with two distant Langmuir probe arrays. The edge floating potentials display differences of several tens of Volts in electron-root wave-heated plasmas. The differences are reduced for higher densities and lower electron temperatures after the ion-root electric field forms at the plasma edge. Neoclassical Monte Carlo simulations estimate the correct order of magnitude for the overall variation in potential and predict the trend observed with the radial electric field. However, for the specific location of the probes, the simulations give differences smaller than those observed experimentally.

Fully-ionised carbon impurity flow is studied in ion-root, neutral beam heated plasmas by means of Charge Exchange Recombination Spectroscopy (CXRS) in the TJ-II stellarator. Perpendicular flows are found to be in reasonable agreement with neoclassical calculations of the radial electric field. The parallel flow of the impurity is obtained at two locations of the same flux surface after subtraction of the calculated Pfirsch-Schl\"uter parallel velocity. For the medium density plasmas studied, $\bar{n}_{\rm e}\in(1.2-2.4)\times 10^{19}$ m$^{-3}$, the measured impurity flow is found to be inconsistent with a total incompressible flow, i.e. $\nabla\cdot{\bf u}_z\ne0$, thus contradicting the usual assumption of a constant density on each flux surface. The experimentally observed velocity deviations are compared with the parallel return flow calculated from a modelled impurity density redistribution driven by ion-impurity friction. Although the calculated return flow substantially modifies the incompressible velocity pattern, the modifications do not explain the in-surface variations of impurity parallel mass flow at the precise locations of the CXRS measurements. Small inhomogeneities of the electrostatic potential in a surface are also shown to affect the impurity redistribution but do not provide a better understanding of the measurements., Comment: 22 pages, 11 figures

Rotation is favorable for confinement, but a stellarator can rotate at high speeds if and only if it is sufficiently close to quasisymmetry. This article investigates how close it needs to be. For a magnetic field $\mathbf{B} = \mathbf{B}_0 + \alpha \mathbf{B}_1$, where $\mathbf{B}_0$ is quasisymmetric, $\alpha\mathbf{B}_1$ is a deviation from quasisymmetry, and $\alpha\ll 1$, the stellarator can rotate at high velocities if $\alpha < \epsilon^{1/2}$, with $\epsilon$ the ion Larmor radius over the characteristic variation length of $\mathbf{B}_0$. The cases in which this result may break down are discussed. If the stellarator is sufficiently quasisymmetric in the above sense, the rotation profile, and equivalently, the long-wavelength radial electric field, are not set neoclassically; instead, they can be affected by turbulent transport. Their computation requires the $O(\epsilon^2)$ pieces of both the turbulent and the long-wavelength components of the distribution function. This article contains the first step towards a formulation to calculate the rotation profile by providing the equations determining the long-wavelength components of the $O(\epsilon^2)$ pieces., Comment: 59 pages, 2 figures. arXiv admin note: text overlap with arXiv:1204.1509

The drift kinetic equation is solved for low density TJ-II plasmas employing slowly varying, time-dependent profiles. This allows to simulate density ramp-up experiments and to describe from first principles the formation and physics of the radial electric field shear, which is associated to the transition from electron to ion root. We show that the range of frequencies of plasma potential fluctuations in which zonal flows are experimentally observed is neoclassically undamped in a neighborhood of the transition. This makes the electron root regime of stellarators, close to the transition to ion root, a propitious regime for the study of zonal-flow evolution. We present simulations of collisionless relaxation of zonal flows, in the sense of the Rosenbluth and Hinton test, that show an oscillatory behaviour in qualitative agreement with the experiment close to the transition.

Poloidal and toroidal velocities of fully-ionised carbon are measured by means of charge exchange recombination spectroscopy (CXRS) in the TJ-II stellarator. We present a detailed treatment of the 3D geometry and show that flow measurements, performed at different locations of the same flux surface, are compatible with flow incompressibility for the low density plasmas under study (line averaged electron densities $\bar{n}_e\,\le\,10^{19}$ m$^{-3}$). Furthermore, comparison with neoclassical calculations shows quantitative agreement with the measured radial electric field and ion bootstrap parallel flow in the absence of an external momentum input., Comment: 17 pages, 10 figures

The dynamics of fluctuating electric field structures in the edge of the TJ-II stellarator, that display zonal flow-like traits, is studied. These structures have been shown to be global and affect particle transport dynamically [J.A. Alonso et al., Nucl. Fus. 52 063010 (2012)]. In this article we discuss possible drive (Reynolds stress) and damping (Neoclassical viscosity, geodesic transfer) mechanisms for the associated ExB velocity. We show that: (a) while the observed turbulence-driven forces can provide the necessary perpendicular acceleration, a causal relation could not be firmly established, possibly because of the locality of the Reynolds stress measurements, (b) the calculated neoclassical viscosity and damping times are comparable to the observed zonal flow relaxation times, and (c) although an accompanying density modulation is observed to be associated to the zonal flow, it is not consistent with the excitation of pressure side-bands, like those present in geodesic acoustic oscillations, caused by the compression of the ExB flow field.

The drift kinetic equation is solved for low density TJ-II plasmas employing slowly varying, time-dependent profiles. This allows to simulate density ramp-up experiments and describe from first principles the formation and physics of the radial electric field shear layer. The main features of the transition are perfectly captured by the calculation, and good quantitative agreement is also found. The results presented here, that should be valid for other non-quasisymmetric stellarators, provide a fundamental explanation for a wealth of experimental observations connected to the shear layer emergence in TJ-II. The key quantity is the neoclassical viscosity, which is shown to go smoothly to zero when the critical density is approached from below. This makes it possible for turbulence-related phenomena, and particularly zonal flows, to arise in the neighborhood of the transition., Comment: 5 pages, 5 figures, accepted for publication in Physical Review Letters

Calculations of the monoenergetic radial diffusion coefficients are presented for several configurations of the TJ-II stellarator usually explored in operation. The neoclassical radial fluxes and the ambipolar electric field for the standard configuration are then studied for three different collisionality regimes, obtaining precise results in all cases.

Calculations of the bootstrap current for the TJ-II stellarator are presented. DKES and NEO-MC codes are employed; the latter has allowed, for the first time, the precise computation of the bootstrap transport coefficient in the long mean free path regime of this device. The low error bars allow a precise convolution of the monoenergetic coefficients, which is confirmed by error analysis. The radial profile of the bootstrap current is presented for the first time for the 100_44_64 configuration of TJ-II for three different collisionality regimes. The bootstrap coefficient is then compared to that of other configurations of TJ-II regularly operated. The results show qualitative agreement with toroidal current measurements; precise comparison with real discharges is ongoing.

We perform numerical simulations, including parallel tempering, on the Potts glass model with binary random quenched couplings using the JANUS application-oriented computer. We find and characterize a glassy transition, estimating the location of the transition and the value of the critical exponents. We show that there is no ferromagnetic transition in a large temperature range around the glassy critical temperature. We also compare our results with those obtained recently on the "random permutation" Potts glass., Comment: 7 pages and 3 figures. Corrected minor typos

We study numerically the nonequilibrium dynamics of the Ising Spin Glass, for a time that spans eleven orders of magnitude, thus approaching the experimentally relevant scale (i.e. {\em seconds}). We introduce novel analysis techniques that allow to compute the coherence length in a model-independent way. Besides, we present strong evidence for a replicon correlator and for overlap equivalence. The emerging picture is compatible with non-coarsening behavior., Comment: 4 pages, 4 eps color figures. Version accepted for publication in Phys. Rev. Lett

We investigate the critical properties of the four-state commutative random permutation glassy Potts model in three and four dimensions by means of Monte Carlo simulation and of a finite size scaling analysis. Thanks to the use of a field programmable gate array we have been able to thermalize a large number of samples of systems with large volume. This has allowed us to observe a spin-glass ordered phase in d=4 and to study the critical properties of the transition. In d=3, our results are consistent with the presence of a Kosterlitz-Thouless transition, but we cannot exclude transient effects due to a value of the lower critical dimension slightly below 3., Comment: 9 pages, 8 Postscript figures

This paper describes JANUS, a modular massively parallel and reconfigurable FPGA-based computing system. Each JANUS module has a computational core and a host. The computational core is a 4x4 array of FPGA-based processing elements with nearest-neighbor data links. Processors are also directly connected to an I/O node attached to the JANUS host, a conventional PC. JANUS is tailored for, but not limited to, the requirements of a class of hard scientific applications characterized by regular code structure, unconventional data manipulation instructions and not too large data-base size. We discuss the architecture of this configurable machine, and focus on its use on Monte Carlo simulations of statistical mechanics. On this class of application JANUS achieves impressive performances: in some cases one JANUS processing element outperfoms high-end PCs by a factor ~ 1000. We also discuss the role of JANUS on other classes of scientific applications., Comment: 11 pages, 6 figures. Improved version, largely rewritten, submitted to Computing in Science & Engineering

We describe the hardwired implementation of algorithms for Monte Carlo simulations of a large class of spin models. We have implemented these algorithms as VHDL codes and we have mapped them onto a dedicated processor based on a large FPGA device. The measured performance on one such processor is comparable to O(100) carefully programmed high-end PCs: it turns out to be even better for some selected spin models. We describe here codes that we are currently executing on the IANUS massively parallel FPGA-based system., Comment: 19 pages, 8 figures; submitted to Computer Physics Communications

We investigate by means of Monte Carlo simulation and Finite-Size Scaling analysis the critical properties of the three dimensional O(5) non linear sigma model and of the antiferromagnetic RP2 model, both of them regularized on a lattice. High accuracy estimates are obtained for the critical exponents, universal dimensionless quantities and critical couplings. It is concluded that both models belong to the same Universality Class, provide that rather non standard identifications are made for the momentum-space propagator of the RP2 model. We have also investigated the phase diagram of the RP2 model extended by a second-neighbor interaction. A rich phase diagram is found, where most phase transitions are first order., Comment: 3 postscript figures

Dedicated machines designed for specific computational algorithms can outperform conventional computers by several orders of magnitude. In this note we describe {\it Ianus}, a new generation FPGA based machine and its basic features: hardware integration and wide reprogrammability. Our goal is to build a machine that can fully exploit the performance potential of new generation FPGA devices. We also plan a software platform which simplifies its programming, in order to extend its intended range of application to a wide class of interesting and computationally demanding problems. The decision to develop a dedicated processor is a complex one, involving careful assessment of its performance lead, during its expected lifetime, over traditional computers, taking into account their performance increase, as predicted by Moore's law. We discuss this point in detail.