Showing 1–50 of 76 results for author: Calvo, I

Piecewise omnigenous stellarators

Authors: J. L. Velasco, I. Calvo, F. J. Escoto, E. Sánchez, H. Thienpondt, F. I. Parra

Abstract: 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 contou… ▽ More 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. △ Less

Submitted 3 September, 2024; v1 submitted 13 May, 2024; originally announced May 2024.

Reduced turbulent transport in the quasi-isodynamic stellarator configuration CIEMAT-QI4

Authors: J. M. García-Regaña, I. Calvo, E. Sánchez, H. Thienpondt, J. L. Velasco, J. A. Capitán

Abstract: CIEMAT-QI4 is a quasi-isodynamic stellarator configuration that simultaneously features very good fast-ion confinement in a broad range of $β$ values, low neoclassical transport and bootstrap current, and ideal magnetohydrodynamic stability up to $β=5\%$. In this paper it is shown that CIEMAT-QI4 also exhibits reduced turbulent transport. This is demonstrated through nonlinear electrostatic simula… ▽ More CIEMAT-QI4 is a quasi-isodynamic stellarator configuration that simultaneously features very good fast-ion confinement in a broad range of $β$ values, low neoclassical transport and bootstrap current, and ideal magnetohydrodynamic stability up to $β=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. △ Less

Submitted 10 September, 2024; v1 submitted 25 April, 2024; originally announced April 2024.

Comments: 12 pages, 9 figures

Influence of the density gradient on turbulent heat transport at ion-scales: an inter-machine study with the gyrokinetic code stella

Authors: H Thienpondt, JM García-Regaña, I Calvo, G Acton, M Barnes

Abstract: Efficient control of turbulent heat transport is crucial for magnetic confinement fusion reactors. This work discusses the complex interplay between density gradients and micro-instabilities, shedding light on their impact on turbulent heat transport in different fusion devices. In particular, the influence of density gradients on turbulent heat transport is investigated through an extensive inter… ▽ More Efficient control of turbulent heat transport is crucial for magnetic confinement fusion reactors. This work discusses the complex interplay between density gradients and micro-instabilities, shedding light on their impact on turbulent heat transport in different fusion devices. In particular, the influence of density gradients on turbulent heat transport is investigated through an extensive inter-machine study, including various stellarators such as W7-X, LHD, TJ-II and NCSX, along with the Asdex-Upgrade tokamak and the tokamak geometry of the Cyclone Base Case (CBC). Linear and nonlinear simulations are performed employing the $δ$f-gyrokinetic code stella across a wide range of parameters to explore the effects of density gradients, temperature gradients, and kinetic electrons. A strong reduction in ion heat flux with increasing density gradients is found in NCSX and W7-X due to the stabilization of temperature-gradient-driven modes without significantly destabilizing density-gradient-driven modes. In contrast, the tokamaks exhibit an increase in ion heat flux with density gradients. Notably, the behavior of ion heat fluxes in stellarators does not align with that of linear growth rates. Additionally, this study provides physical insights into the micro-instabilities, emphasizing the dominance of trapped-electron-modes in CBC, AUG, TJ-II, LHD and NCSX, while both the trapped-electron-mode and the passing-particle-driven universal instability contribute significantly in W7-X. △ Less

Submitted 15 April, 2024; originally announced April 2024.

Comments: 26 pages, 30 figures

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

Authors: F. J. Escoto, J. L. Velasco, I. Calvo, M. Landreman, F. I. Parra

Abstract: 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,… ▽ More 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/ △ Less

Submitted 20 December, 2023; v1 submitted 19 December, 2023; originally announced December 2023.

Electrostatic microturbulence in W7-X: comparison of local gyrokinetic simulations with Doppler reflectometry measurements

Authors: A. González-Jerez, J. M. García-Regaña, I. Calvo, D. Carralero, T. Estrada, E. Sánchez, M. Barnes, the W7-X Team

Abstract: The first experimental campaigns of Wendelstein 7-X (W7-X) have shown that turbulence plays a decisive role in the performance of neoclassically optimized stellarators. This stresses the importance of understanding microturbulence from the theoretical and experimental points of view. To this end, this paper addresses a comprehensive characterization of the turbulent fluctuations by means of nonlin… ▽ More The first experimental campaigns of Wendelstein 7-X (W7-X) have shown that turbulence plays a decisive role in the performance of neoclassically optimized stellarators. This stresses the importance of understanding microturbulence from the theoretical and experimental points of view. To this end, this paper addresses a comprehensive characterization of the turbulent fluctuations by means of nonlinear gyrokinetic simulations performed with the code stella in two W7-X scenarios. In the first part of the paper, the amplitude of the density fluctuations is calculated and compared with measurements obtained by Doppler reflectometry (DR) in the OP1 experimental campaigns. It is found that the trend of the fluctuations along the radius is explained by the access of the DR system to different regions of the turbulence wavenumber spectrum. In the second part of the article, frequency spectra of the density fluctuations and the zonal component of the turbulent flow are numerically characterized for comparisons against future experimental analyses. Both quantities feature broad frequency spectra with dominant frequencies of O(1)-O(10) kHz. △ Less

Submitted 15 December, 2023; originally announced December 2023.

Comments: 27 pages, 16 figures

A polynomial quantum computing algorithm for solving the dualization problem

Authors: Mauro Mezzini, Fernando Cuartero Gomez, Fernando Pelayo, Jose Javier Paulet Gonzales, Hernan Indibil de la Cruz Calvo, Vicente Pascual

Abstract: Given two prime monotone boolean functions $f:\{0,1\}^n \to \{0,1\}$ and $g:\{0,1\}^n \to \{0,1\}$ the dualization problem consists in determining if $g$ is the dual of $f$, that is if $f(x_1, \dots, x_n)= \overline{g}(\overline{x_1}, \dots \overline{x_n})$ for all $(x_1, \dots x_n) \in \{0,1\}^n$. Associated to the dualization problem there is the corresponding decision problem: given two monoton… ▽ More Given two prime monotone boolean functions $f:\{0,1\}^n \to \{0,1\}$ and $g:\{0,1\}^n \to \{0,1\}$ the dualization problem consists in determining if $g$ is the dual of $f$, that is if $f(x_1, \dots, x_n)= \overline{g}(\overline{x_1}, \dots \overline{x_n})$ for all $(x_1, \dots x_n) \in \{0,1\}^n$. Associated to the dualization problem there is the corresponding decision problem: given two monotone prime boolean functions $f$ and $g$ is $g$ the dual of $f$? In this paper we present a quantum computing algorithm that solves the decision version of the dualization problem in polynomial time. △ Less

Submitted 28 August, 2023; originally announced August 2023.

Robust stellarator optimization via flat mirror magnetic fields

Authors: J. L. Velasco, I. Calvo, E. Sánchez, F. I. Parra

Abstract: 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… ▽ More 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 $β$ 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. △ Less

Submitted 30 June, 2023; originally announced June 2023.

Journal ref: Nucl. Fusion 63 126038 (2023)

Reduction or enhancement of stellarator turbulence by impurities

Authors: J. M. García-Regaña, I. Calvo, F. I. Parra, H. Thienpondt

Abstract: A systematic study of the impact of impurities on the turbulent heat fluxes is presented for the stellarator Wendelstein 7-X (W7-X) and, for comparison, the Large Helical Device and ITER. By means of nonlinear multispecies gyrokinetic simulations, it is shown that impurities, depending on the sign of their density gradient, can significantly enhance or reduce turbulent ion heat losses. For the rel… ▽ More A systematic study of the impact of impurities on the turbulent heat fluxes is presented for the stellarator Wendelstein 7-X (W7-X) and, for comparison, the Large Helical Device and ITER. By means of nonlinear multispecies gyrokinetic simulations, it is shown that impurities, depending on the sign of their density gradient, can significantly enhance or reduce turbulent ion heat losses. For the relevant scenario of turbulence reduction, an optimal impurity concentration that minimizes the ion heat diffusivity emerges as a universal feature. This result demonstrates the potential of impurities for controlling turbulence and accessing enhanced confinement regimes in fusion plasmas and, in particular, in W7-X. △ Less

Submitted 6 September, 2024; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 8 pages, 8 figures

Neoclassical transport in strong gradient regions of large aspect ratio tokamaks

Authors: Silvia Trinczek, Felix I. Parra, Peter J. Catto, Iván Calvo, Matt Landreman

Abstract: We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionalit… ▽ More We present a new neoclassical transport model for large aspect ratio tokamaks where the gradient scale lengths are of the size of the poloidal gyroradius. Previous work on neoclassical transport across transport barriers assumed large density and potential gradients but a small temperature gradient, or neglected the gradient of the mean parallel flow. Using large aspect ratio and low collisionality expansions, we relax these restrictive assumptions. We define a new set of variables based on conserved quantities, which simplifies the drift kinetic equation whilst keeping strong gradients, and derive equations describing the transport of particles, parallel momentum and energy by ions in the banana regime. The poloidally varying parts of density and electric potential are included. Studying contributions from both passing and trapped particles, we show that the resulting transport is dominated by trapped particles. We find that a non-zero neoclassical particle flux requires parallel momentum input which could be provided through interaction with turbulence or impurities. We derive upper and lower bounds for the energy flux across a transport barrier in both temperature and density and present example profiles and fluxes. △ Less

Submitted 19 April, 2023; v1 submitted 17 January, 2023; originally announced January 2023.

Comments: 52 pages, 13 figures, submitted to Journal of Plasma Physics

A quasi-isodynamic configuration with good confinement of fast ions at low plasma $β$

Authors: E. Sánchez, J. L. Velasco, I. Calvo, S. Mulas

Abstract: A new quasi-isodynamic stellarator configuration optimized for the confinement of energetic ions at low plasma $β$ 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, lo… ▽ More A new quasi-isodynamic stellarator configuration optimized for the confinement of energetic ions at low plasma $β$ 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 $β= 5%$, reduced effective ripple, with $ε_{eff} < 0.5%$ in the plasma core. Even at low $β$, the configuration approximately satisfies the maximum-$J$ property, and the confinement of fast ions is good at $β\sim 1.5%$ and becomes excellent at reactor values, $β\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. △ Less

Submitted 9 December, 2022; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: 21 pages, 16 figures

Prevention of core particle depletion in stellarators by turbulence

Authors: H. Thienpondt, J. M. García-Regaña, I. Calvo, J. A. Alonso, J. L. Velasco, A. González-Jerez, M. Barnes, K. Brunner, O. Ford, G. Fuchert, J. Knauer, E. Pasch, L. Vanó, the Wendelstein 7-X team

Abstract: 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… ▽ More 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. △ Less

Submitted 8 November, 2023; v1 submitted 9 September, 2022; originally announced September 2022.

Comments: 8 pages, 6 figures, updated the title and incorporated the appendix into the main letter

Journal ref: Physical Review Research 5 (2023) L022053

Finite orbit width effects in large aspect ratio stellarators

Authors: Vincent d'Herbemont, Felix I. Parra, Ivan Calvo, Jose Luis Velasco

Abstract: New orbit averaged equations for low collisionality neoclassical fluxes in large aspect ratio stellarators with mirror ratios close to unity are derived. The equations retain finite orbit width effects by employing the second adiabatic invariant $J$ as a velocity space coordinate and they have been implemented in the orbit-averaged neoclassical code KNOSOS. The equations are used to study the… ▽ More New orbit averaged equations for low collisionality neoclassical fluxes in large aspect ratio stellarators with mirror ratios close to unity are derived. The equations retain finite orbit width effects by employing the second adiabatic invariant $J$ as a velocity space coordinate and they have been implemented in the orbit-averaged neoclassical code KNOSOS. The equations are used to study the $1/ν$ regime and the lower collisionality regimes. For generic large aspect ratio stellarators with mirror ratios close to unity, as the collision frequency decreases, the $1/ν$ regime transitions directly into the $ν$ regime, without passing through a $\sqrtν$ regime. An explicit formula for the neoclassical fluxes in the $ν$ regime is obtained. The formula includes the effect of particles that transition between different types of wells. While these transitions produce stochastic scattering independent of the value of the collision frequency in velocity space, the diffusion in real space remains proportional to the collision frequency. The $\sqrtν$ regime is only recovered in large aspect ratio stellarators close to omnigeneity: large aspect ratio stellarators with large mirror ratios and optimized large aspect ratio stellarators with mirror ratios close to unity. Neoclassical transport in large aspect ratio stellarators with large mirror ratios can be calculated with the orbit-averaged equations derived by \cite{calvo17}. In these stellarators, the $\sqrtν$ regime exists in the collisionality interval $ε|\ln ε| \ll ν_{ii} R a/ρ_i v_{ti} \ll 1/ε$. In optimized large aspect ratio stellarators with mirror ratios close to unity, the $\sqrtν$ regime occurs in an interval of collisionality that depends on the deviation from omnigeneity $δ$: $δ^2 |\ln δ| \ll ν_{ii} R a/ρ_i v_{ti} \ll 1$. △ Less

Submitted 2 September, 2022; v1 submitted 2 February, 2022; originally announced February 2022.

Comments: Accepted in the Journal of Plasma Physics. 79 pages, 15 figures

On the role of density fluctuations in the core turbulent transport of Wendelstein 7-X

Authors: D. Carralero, T. Estrada, E. Maragkoudakis, T. Windisch, J. A. Alonso, J. L. Velasco, O. Ford, M. Jakubowski, S. Lazerson, M. Beurskens, S. Bozhenkov, I. Calvo, H. Damm, G. Fuchert, J. M. García-Regaña, U. Höfel, N. Marushchenko, N. Pablant, E. Sánchez, H. M. Smith, E. Pasch, T. Stange

Abstract: 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 $η_i = L_n/L_{T_i}$ [D. Carralero et al., Nucl. Fusion, 2021]. In order to confirm the role of… ▽ More 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 $η_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 $η_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. △ Less

Submitted 1 October, 2021; originally announced October 2021.

Comments: Submitted to PPCF, 47th EPS special issue

Physics design point of high-field stellarator reactors

Authors: J. A. Alonso, I. Calvo, D. Carralero, J. L. Velasco, J. M. García-Regaña, I. Palermo, D. Rapisarda

Abstract: 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.… ▽ More 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 △ Less

Submitted 25 December, 2021; v1 submitted 30 September, 2021; originally announced September 2021.

Electrostatic gyrokinetic simulations in Wendelstein 7-X geometry: benchmark between the codes stella and GENE

Authors: A. González-Jerez, P. Xanthopoulos, J. M. García-Regaña, I. Calvo, J. Alcusón, A. Bañón-Navarro, M. Barnes, F. I. Parra, J. Geiger

Abstract: The first experimental campaigns have proven that, due to the optimization of the magnetic configuration with respect to neoclassical transport, the contribution of turbulence is essential to understand and predict the total particle and energy transport in Wendelstein 7-X (W7-X). This has spurred much work on gyrokinetic modelling for the interpretation of the available experimental results and f… ▽ More The first experimental campaigns have proven that, due to the optimization of the magnetic configuration with respect to neoclassical transport, the contribution of turbulence is essential to understand and predict the total particle and energy transport in Wendelstein 7-X (W7-X). This has spurred much work on gyrokinetic modelling for the interpretation of the available experimental results and for the preparation of the next campaigns. At the same time, new stellarator gyrokinetic codes have just been or are being developed. It is therefore desirable to have a sufficiently complete, documented and verified set of gyrokinetic simulations in W7-X geometry against which new codes or upgrades of existing codes can be tested and benchmarked. This paper attemps to provide such a set of simulations in the form of a comprehensive benchmark between the recently developed code stella and the well-established code GENE. The benchmark consists of electrostatic gyrokinetic simulations in W7-X magnetic geometry and includes different flux tubes, linear ion-temperature-gradient (ITG) and trapped-electron-mode (TEM)} stability analyses, computation of linear zonal flow responses and calculation of ITG-driven heat fluxes. △ Less

Submitted 13 July, 2021; originally announced July 2021.

A model for the fast evaluation of prompt losses of energetic ions in stellarators

Authors: J. L. Velasco, I. Calvo, S. Mulas, E. Sánchez, F. I. Parra, Á. Cappa, the W7-X team

Abstract: 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… ▽ More 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. △ Less

Submitted 6 July, 2021; v1 submitted 10 June, 2021; originally announced June 2021.

Turbulent transport of impurities in 3D devices

Authors: J. M. García-Regaña, M. Barnes, I. Calvo, A. González-Jerez, H. Thienpondt, E. Sánchez, F. I. Parra, D. St. -Onge

Abstract: A large diffusive turbulent contribution to the radial impurity transport in Wendelstein 7-X (W7-X) plasmas has been experimentally inferred during the first campaigns and numerically confirmed by means of gyrokinetic simulations with the code stella. In general, the absence of strong impurity accumulation during the initial W7-X campaigns is attributed to this diffusive term. In the present work… ▽ More A large diffusive turbulent contribution to the radial impurity transport in Wendelstein 7-X (W7-X) plasmas has been experimentally inferred during the first campaigns and numerically confirmed by means of gyrokinetic simulations with the code stella. In general, the absence of strong impurity accumulation during the initial W7-X campaigns is attributed to this diffusive term. In the present work the diffusive contribution is also calculated in other stellarator plasmas. In particular, the diffusion (D) and convection (V) coefficients of carbon and iron impurities produced by ion-temperature-gradient (ITG) turbulence are obtained for W7-X, LHD, TJ-II and NCSX. The results show that, although the size of D and V can differ across the four devices, inward convection is found for all of them. For W7-X, TJ-II and NCSX the two coefficients are comparable and the turbulent peaking factor is surprisingly similar. In LHD, appreciably weaker diffusive and convective impurity transport and significantly larger turbulent peaking factor are predicted. All this suggests that ITG turbulence, although not strongly, would lead to negative impurity density gradients in stellarators. Then, considering mixed ITG/Trapped Electron Mode (TEM) turbulence for the specific case of W7-X, it has been quantitatively assessed to what degree pellet fueled reduced turbulence scenarios feature reduced turbulent transport of impurities as well. The results for trace iron impurities show that, although their turbulent transport is not entirely suppressed, a significant reduction of V and a stronger decrease of D are found. Although the diffusion is still above neoclassical levels, the neoclassical convection would gain under such conditions a greater specific weight on the dynamics of impurities in comparison with standard ECRH scenarios. △ Less

Submitted 22 October, 2021; v1 submitted 9 June, 2021; originally announced June 2021.

Comments: 16 pages, 8 figures

Journal ref: Nucl. Fusion 61 116019 2021

Gyrokinetic simulations in stellarators using different computational domains

Authors: E. Sánchez, J. M. García-Regaña, A. Bañón Navarro, J. H. E. Proll, C. Mora Moreno, A. González-Jerez, I. Calvo, R. Kleiber, J. Riemann, J. Smoniewski, M. Barnes, F. I. Parra

Abstract: In this work, we compare gyrokinetic simulations in stellarators using different computational domains, namely, flux tube, full-flux-surface, and radially global domains. Two problems are studied: the linear relaxation of zonal flows and the linear stability of ion temperature gradient (ITG) modes. Simulations are carried out with the codes EUTERPE, GENE, GENE-3D, and stella in magnetic configurat… ▽ More In this work, we compare gyrokinetic simulations in stellarators using different computational domains, namely, flux tube, full-flux-surface, and radially global domains. Two problems are studied: the linear relaxation of zonal flows and the linear stability of ion temperature gradient (ITG) modes. Simulations are carried out with the codes EUTERPE, GENE, GENE-3D, and stella in magnetic configurations of LHD and W7-X using adiabatic electrons. The zonal flow relaxation properties obtained in different flux tubes are found to differ with each other and with the radially global result, except for sufficiently long flux tubes, in general. The flux tube length required for convergence is configuration-dependent. Similarly, for ITG instabilities, different flux tubes provide different results, but the discrepancy between them diminishes with increasing flux tube length. Full-flux-surface and flux tube simulations show good agreement in the calculation of the growth rate and frequency of the most unstable modes in LHD, while for W7-X differences in the growth rates are found between the flux tube and the full-flux-surface domains. Radially global simulations provide results close to the full-flux-surface ones. The radial scale of unstable ITG modes is studied in global and flux tube simulations finding that in W7-X, the radial scale of the most unstable modes depends on the binormal wavenumber, while in LHD no clear dependency is found. △ Less

Submitted 5 June, 2021; originally announced June 2021.

Comments: submitted to Nuclear Fusion

Fast simulations for large aspect ratio stellarators with the neoclassical code KNOSOS

Authors: J. L. Velasco, I. Calvo, F. I. Parra, V. d'Herbemont, H. M. Smith, D. Carralero, T. Estrada, the W7-X team

Abstract: 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 ve… ▽ More 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. △ Less

Submitted 3 June, 2021; originally announced June 2021.

Comments: IAEA paper, submitted to Nuclear Fusion

An experimental characterization of core turbulence regimes in Wendelstein 7-X

Authors: D. Carralero, T. Estrada, E. Maragkoudakis, T. Windisch, J. A. Alonso, M. Beurskens, S. Bozhenkov, I. Calvo, H. Damm, O. Ford, G. Fuchert, J. M. García-Regaña, N. Pablant, E. Sánchez, E. Pasch, J. L. Velasco, the Wendelstein 7-X team

Abstract: 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 cla… ▽ More 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ρ_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 $η_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. △ Less

Submitted 11 May, 2021; originally announced May 2021.

Comments: Manuscript submitted to Nuclear Fusion in May 2021

Study on impurity hole plasmas by global neoclassical simulation

Authors: Keiji Fujita, Shinsuke Satake, Masanori Nunami, José Manuel García-Regaña, José Luis Velasco, Iván Calvo

Abstract: An impurity hole observed in the Large Helical Device (LHD) is a hollow density profile of an impurity ion species formed in the core plasma where the negative (inward-pointing) ambipolar radial electric field ($E_r$) exists. Although local neoclassical models have predicted the sign of $E_r$ in impurity hole plasmas is negative for the entire minor radius, an experimental measurement of an impuri… ▽ More An impurity hole observed in the Large Helical Device (LHD) is a hollow density profile of an impurity ion species formed in the core plasma where the negative (inward-pointing) ambipolar radial electric field ($E_r$) exists. Although local neoclassical models have predicted the sign of $E_r$ in impurity hole plasmas is negative for the entire minor radius, an experimental measurement of an impurity hole plasma has shown that the $E_r$ changes the sign from negative to positive along the minor radius. In the present work, we investigate neoclassical impurity transport in an impurity hole plasma using a global neoclassical simulation code FORTEC-3D. The variation of electrostatic potential on flux surface ($Φ_1$) is evaluated from the quasi-neutrality condition in multi-ion-species plasma by the global simulation. The ambipolar $E_r$ and neoclassical fluxes are determined in solving a global drift-kinetic equation including the effect of $Φ_1$. By the global simulation, we show that an $E_r$ which changes the sign along the radius is obtained as a solution of the ambipolar condition and with such an $E_r$, impurity carbon flux can be outwardly directed even where $E_r<0$ and the carbon density profile is hollow around the magnetic axis. Furthermore, it is found that the outward carbon flux is only a factor 2-3 from balancing the modeled inward turbulent flux. Our result indicates that we have moved one step closer to reproducing the impurity transport in impurity hole plasmas by kinetic simulation. △ Less

Submitted 25 July, 2021; v1 submitted 27 April, 2021; originally announced April 2021.

Journal ref: Nucl. Fusion 61 086025 (2021)

Comparison of local and global gyrokinetic calculations of collisionless zonal flow damping in quasi-symmetric stellarators

Authors: J. Smoniewski, E. Sánchez, I. Calvo, M. J. Pueschel, J. N. Talmadge

Abstract: The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (N… ▽ More The linear collisionless damping of zonal flows is calculated for quasi-symmetric stellarator equilibria in flux-tube, flux-surface, and full-volume geometry. Equilibria are studied from the quasi-helical symmetry configuration of the Helically Symmetric eXperiment (HSX), a broken symmetry configuration of HSX, and the quasi-axial symmetry geometry of the National Compact Stellarator eXperiment (NCSX). Zonal flow oscillations and long-time damping affect the zonal flow evolution, and the zonal flow residual goes to zero for small radial wavenumber. The oscillation frequency and damping rate depend on the bounce-averaged radial particle drift in accordance with theory. While each flux tube on a flux surface is unique, several different flux tubes in HSX or NCSX can reproduce the zonal flow damping from a flux-surface calculation given an adequate parallel extent. The flux-surface or flux-tube calculations can accurately reproduce the full-volume long-time residual for moderate $k_x$, but the oscillation and damping time scales are longer in local representations, particularly for small $k_x$ approaching the system size. △ Less

Submitted 13 March, 2021; v1 submitted 22 December, 2020; originally announced December 2020.

Comments: The following article has been accepted by Physics of Plasmas. After it is published, it will be found at https://aip.scitation.org/journal/php. 33 pages, 18 figures

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

Authors: J. M. García-Regaña, M. Barnes, I. Calvo, F. I. Parra, J. Alcusón, R. Davies, A. González-Jerez, A. Mollén, E. Sánchez, J. L. Velasco, A. Zocco

Abstract: 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;… ▽ More 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. △ Less

Submitted 17 August, 2020; originally announced August 2020.

Comments: 19 pages, 10 figures, 2 tables

Journal ref: J. Plasma Phys. 87 (2021) 855870103

KNOSOS: a fast orbit-averaging neoclassical code for stellarator geometry

Authors: J. L. Velasco, I. Calvo, F. I. Parra, J. M. García-Regaña

Abstract: 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 featur… ▽ More 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. △ Less

Submitted 22 June, 2020; v1 submitted 30 August, 2019; originally announced August 2019.

Impact of main ion pressure anisotropy on stellarator impurity transport

Authors: Ivan Calvo, Felix I. Parra, J. L. Velasco, J. M. García-Regaña

Abstract: 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… ▽ More 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üter, plateau and $1/ν$ 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. △ Less

Submitted 12 December, 2019; v1 submitted 19 July, 2019; originally announced July 2019.

Comments: 40 pages, 4 figures. Published version

Journal ref: Nucl. Fusion 60, 016035 (2020)

Electrostatic potential variations on stellarator magnetic surfaces in low collisionality regimes

Authors: Ivan Calvo, J. L. Velasco, Felix I. Parra, J. Arturo Alonso, J. M. García-Regaña

Abstract: 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 dema… ▽ More 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/ν$, $\sqrtν$ 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. △ Less

Submitted 21 August, 2018; v1 submitted 30 April, 2018; originally announced April 2018.

Comments: 30 pages, 12 figures, 1 table. Published version

Journal ref: J. Plasma Phys. 84, 905840407 (2018)

On-surface potential and radial electric field variations in electron root stellarator plasmas

Authors: J. M. García-Regaña, T. Estrada, I. Calvo, J. L. Velasco, J. A. Alonso, D. Carralero, R. Kleiber, M. Landreman, A. Mollén, E. Sánchez, C. Slaby

Abstract: 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… ▽ More 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 $Φ_1$ as a possible correction term to the total radial electric field. Here $Φ_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), $Φ_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 $Φ_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. △ Less

Submitted 27 April, 2018; originally announced April 2018.

Comments: 23 pages, 10 figures

Stellarator impurity flux driven by electric fields tangent to magnetic surfaces

Authors: Ivan Calvo, Felix I. Parra, J. L. Velasco, J. Arturo Alonso, J. M. García-Regaña

Abstract: 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. Le… ▽ More 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. △ Less

Submitted 16 October, 2018; v1 submitted 15 March, 2018; originally announced March 2018.

Comments: 10 pages, 1 figure. Published version

Journal ref: Nucl. Fusion 58, 124005 (2018)

Oscillatory relaxation of zonal flows in a multi-species stellarator plasma

Authors: E. Sánchez, I. Calvo, J. L. Velasco, F. Medina, A. Alonso, P. Monreal, R. Kleiber, the TJ-II team

Abstract: 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. I… ▽ More 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. △ Less

Submitted 29 January, 2018; originally announced January 2018.

Comments: 11 figures, 22 pages

Large tangential electric fields in plasmas close to temperature screening

Authors: J. L. Velasco, I. Calvo, J. M. García-Regaña, F. I. Parra, S. Satake, J. A. Alonso, the LHD team

Abstract: 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… ▽ More 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). △ Less

Submitted 11 December, 2017; originally announced December 2017.

Semianalytical calculation of the zonal-flow oscillation frequency in stellarators

Authors: Pedro Monreal, Edilberto Sánchez, Iván Calvo, Andrés Bustos, Félix I. Parra, Alexey Mishchenko, Axel Könies, Ralf Kleiber

Abstract: Due to their capability to reduce turbulent transport in magnetized plasmas, understanding the dynamics of zonal flows is an important problem in the fusion programme. Since the pioneering work by Rosenbluth and Hinton in axisymmetric tokamaks, it is known that studying the linear and collisionless relaxation of zonal flow perturbations gives valuable information and physical insight. Recently, th… ▽ More Due to their capability to reduce turbulent transport in magnetized plasmas, understanding the dynamics of zonal flows is an important problem in the fusion programme. Since the pioneering work by Rosenbluth and Hinton in axisymmetric tokamaks, it is known that studying the linear and collisionless relaxation of zonal flow perturbations gives valuable information and physical insight. Recently, the problem has been investigated in stellarators and it has been found that in these devices the relaxation process exhibits a characteristic feature: a damped oscillation. The frequency of this oscillation might be a relevant parameter in the regulation of turbulent transport, and therefore its efficient and accurate calculation is important. Although an analytical expression can be derived for the frequency, its numerical evaluation is not simple and has not been exploited systematically so far. Here, a numerical method for its evaluation is considered, and the results are compared with those obtained by calculating the frequency from gyrokinetic simulations. This "semianalytical" approach for the determination of the zonal-flow frequency reveals accurate and faster than the one based on gyrokinetic simulations. △ Less

Submitted 4 May, 2017; v1 submitted 10 January, 2017; originally announced January 2017.

Comments: 30 pages, 14 figures

The effect of tangential drifts on neoclassical transport in stellarators close to omnigeneity

Authors: Ivan Calvo, Felix I. Parra, J. L. Velasco, J. Arturo Alonso

Abstract: 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… ▽ More 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/ν$ 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/ν$ regime. Then, it is shown that neoclassical transport is determined by two layers in phase space. One of the layers corresponds to the $\sqrtν$ 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ν$ 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. △ Less

Submitted 10 March, 2017; v1 submitted 19 October, 2016; originally announced October 2016.

Comments: 36 pages. Version to be published in Plasma Physics and Controlled Fusion

Journal ref: Plasma Physics and Controlled Fusion 59, 055014 (2017)

Observation of oscillatory radial electric field relaxation in a helical plasma

Authors: J. A. Alonso, E. Sanchez, I. Calvo, J. L. Velasco, S. Perfilov, A. Chmyga, L. G. Eliseev, L. I. Krupnik, T. Estrada, R. Kleiber, K. J. McCarthy, A. V. Melnikov, P. Monreal, F. I. Parra, A. I. Zhezhera, the TJ-II Team

Abstract: 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… ▽ More 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. △ Less

Submitted 1 September, 2016; originally announced September 2016.

Journal ref: Phys. Rev. Lett. 118, 185002 (2017)

Moderation of neoclassical impurity accumulation in high temperature plasmas of helical devices

Authors: J. L. Velasco, I. Calvo, S. Satake, A. Alonso, M. Nunami, M. Yokoyama, M. Sato, T. Estrada, J. M. Fontdecaba, M. Liniers, K. J. McCarthy, F. Medina, B. Ph Van Milligen, M. Ochando, F. Parra, H. Sugama, A. Zhezhera, the LHD experimental team, the TJ-II team

Abstract: 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 st… ▽ More 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. △ Less

Submitted 28 March, 2017; v1 submitted 26 July, 2016; originally announced July 2016.

Journal ref: Nuclear Fusion 57, 016016 (2017)

Parallel impurity dynamics in the TJ-II stellarator

Authors: J A Alonso, J L Velasco, I Calvo, T Estrada, J M Fontdecaba, J M García-Regaña, J Geiger, M Landreman, K J McCarthy, F Medina, B Ph Van Milligen, M A Ochando, F I Parra, the TJ-II Team, the W7-X Team

Abstract: 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,… ▽ More 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. △ Less

Submitted 4 July, 2016; v1 submitted 20 April, 2016; originally announced April 2016.

Journal ref: 2016 Plasma Phys. Control. Fusion 58 074009

Particle transport after pellet injection in the TJ-II stellarator

Authors: J. L. Velasco, K. J. McCarthy, N. Panadero, S. Satake, D. López-Bruna, A. Alonso, I. Calvo, T. Estrada, J. M. Fontdecaba, J. Hernández, R. García, F. Medina, M. Ochando, I. Pastor, S. Perfilov, E. Sánchez, A. Soleto, B. Ph. Van Milligen, A. Zhezhera, the TJ-II team

Abstract: 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 ca… ▽ More 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. △ Less

Submitted 28 March, 2017; v1 submitted 29 January, 2016; originally announced January 2016.

Journal ref: Plasma Physics and Controlled Fusion 58, 084004 (2016)

Residual zonal flows in tokamaks and stellarators at arbitrary wavelengths

Authors: P. Monreal, I. Calvo, E. Sánchez, F. I. Parra, A. Bustos, A. Könies, R. Kleiber, T. Görler

Abstract: 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, a… ▽ More 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. △ Less

Submitted 13 October, 2015; v1 submitted 12 May, 2015; originally announced May 2015.

Comments: 31 pages, 15 figures

Journal ref: Plasma Physics and Controlled Fusion 58, 045018 (2016)

Radial transport of toroidal angular momentum in tokamaks

Authors: Ivan Calvo, Felix I. Parra

Abstract: The radial flux of toroidal angular momentum is needed to determine tokamak intrinsic rotation profiles. Its computation requires knowledge of the gyrokinetic distribution functions and turbulent electrostatic potential to second-order in $ε= ρ/L$, where $ρ$ is the ion Larmor radius and $L$ is the variation length of the magnetic field. In this article, a complete set of equations to calculate the… ▽ More The radial flux of toroidal angular momentum is needed to determine tokamak intrinsic rotation profiles. Its computation requires knowledge of the gyrokinetic distribution functions and turbulent electrostatic potential to second-order in $ε= ρ/L$, where $ρ$ is the ion Larmor radius and $L$ is the variation length of the magnetic field. In this article, a complete set of equations to calculate the radial transport of toroidal angular momentum in any tokamak is presented. In particular, the $O(ε^2)$ equations for the turbulent components of the distribution functions and electrostatic potential are given for the first time without assuming that the poloidal magnetic field over the magnetic field strength is small. △ Less

Submitted 29 April, 2015; v1 submitted 2 December, 2014; originally announced December 2014.

Comments: 39 pages. Version to be published in Plasma Physics and Controlled Fusion

Journal ref: Plasma Physics and Controlled Fusion 57, 075006 (2015)

Less constrained omnigeneous stellarators

Authors: Felix I. Parra, Ivan Calvo, Per Helander, Matt Landreman

Abstract: A stellarator is said to be omnigeneous if all particles have vanishing average radial drifts. In omnigeneous stellarators, particles are perfectly confined in the absence of turbulence and collisions, whereas in non-omnigeneous configurations, particle can drift large radial distances. One of the consequences of omnigeneity is that the unfavorable inverse scaling with collisionality of the stella… ▽ More A stellarator is said to be omnigeneous if all particles have vanishing average radial drifts. In omnigeneous stellarators, particles are perfectly confined in the absence of turbulence and collisions, whereas in non-omnigeneous configurations, particle can drift large radial distances. One of the consequences of omnigeneity is that the unfavorable inverse scaling with collisionality of the stellarator neoclassical fluxes disappears. In the pioneering and influential article [Cary~J~R and Shasharina~S~G 1997 {\it Phys. Plasmas} {\bf 4} 3323], the conditions that the magnetic field of a stellarator must satisfy to be omnigeneous are derived. However, reference [Cary~J~R and Shasharina~S~G 1997 {\it Phys. Plasmas} {\bf 4} 3323] only considered omnigeneous stellarators in which all the minima of the magnetic field strength on a flux surface must have the same value. The same is assumed for the maxima. We show that omnigenenous magnetic fields can have local minima and maxima with different values. Thus, the parameter space in which omnigeneous stellarators are possible is larger than previously expected. The analysis presented in this article is only valid for orbits with vanishing radial width, and in principle it is not applicable to energetic particles. However, one would expect that improving neoclassical confinement would improve energetic particle confinement. △ Less

Submitted 7 February, 2015; v1 submitted 21 October, 2014; originally announced October 2014.

Comments: 10 pages, 5 figures

Journal ref: Nucl. Fusion 55, 033005 (2015)

Equivalence of two independent calculations of the higher order guiding center Lagrangian

Authors: F. I. Parra, I. Calvo, J. W. Burby, J. Squire, H. Qin

Abstract: The difference between the guiding center phase-space Lagrangians derived in [J.W. Burby, J. Squire, and H. Qin, Phys. Plasmas {\bf 20}, 072105 (2013)] and [F.I. Parra, and I. Calvo, Plasma Phys. Control. Fusion {\bf 53}, 045001 (2011)] is due to a different definition of the guiding center coordinates. In this brief communication the difference between the guiding center coordinates is calculated… ▽ More The difference between the guiding center phase-space Lagrangians derived in [J.W. Burby, J. Squire, and H. Qin, Phys. Plasmas {\bf 20}, 072105 (2013)] and [F.I. Parra, and I. Calvo, Plasma Phys. Control. Fusion {\bf 53}, 045001 (2011)] is due to a different definition of the guiding center coordinates. In this brief communication the difference between the guiding center coordinates is calculated explicitly. △ Less

Submitted 7 February, 2015; v1 submitted 24 July, 2014; originally announced July 2014.

Comments: 6 pages

Journal ref: Phys. Plasmas 21, 104506 (2014)

Flow damping in stellarators close to quasisymmetry

Authors: Ivan Calvo, Felix I. Parra, J. L. Velasco, J. Arturo Alonso

Abstract: 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… ▽ More 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/ν$ regime. Denoting by $α$ a parameter that gives the size of the deviation from quasisymmetry, it was proven that particle fluxes do not scale with $α^{3/2}$, as typically claimed, but with $α$. 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. △ Less

Submitted 13 October, 2014; v1 submitted 4 July, 2014; originally announced July 2014.

Comments: 14 pages, 2 figures. To appear in Plasma Physics and Controlled Fusion

Journal ref: Plasma Physics and Controlled Fusion 57, 014014 (2015)

Optimizing stellarators for large flows

Authors: Ivan Calvo, Felix I. Parra, J. Arturo Alonso, J. L. Velasco

Abstract: 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… ▽ More 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 $α$ 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 $α^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(|α|)$ deviation from quasisymmetry. Heuristic estimations in the literature incorrectly predicted an $O(|α|^{3/2})$ deviation. △ Less

Submitted 5 May, 2014; originally announced May 2014.

Comments: 24 pages, 2 figures. To appear in Plasma Physics and Controlled Fusion

Journal ref: Plasma Physics and Controlled Fusion 56, 094003 (2014)

When omnigeneity fails

Authors: Felix I. Parra, Ivan Calvo, Jose Luis Velasco, J. Arturo Alonso

Abstract: A generic non-symmetric magnetic field does not confine magnetized charged particles for long times due to secular magnetic drifts. Stellarator magnetic fields should be omnigeneous (that is, designed such that the secular drifts vanish), but perfect omnigeneity is technically impossible. There always are small deviations from omnigeneity that necessarily have large gradients. The amplification of… ▽ More A generic non-symmetric magnetic field does not confine magnetized charged particles for long times due to secular magnetic drifts. Stellarator magnetic fields should be omnigeneous (that is, designed such that the secular drifts vanish), but perfect omnigeneity is technically impossible. There always are small deviations from omnigeneity that necessarily have large gradients. The amplification of the energy flux caused by a deviation of size $ε$ is calculated and it is shown that the scaling with $ε$ of the amplification factor can be as large as linear. In opposition to common wisdom, most of the transport is not due to particles trapped in ripple wells, but to the perturbed motion of particles trapped in the omnigeneous magnetic wells around their bounce points. △ Less

Submitted 4 April, 2014; originally announced April 2014.

Comments: 6 pages, 2 figures

Electrostatic potential variations along flux surfaces in stellarators

Authors: M. A. Pedrosa, J. A. Alonso, J. M. García-Regaña, C. Hidalgo, J. L. Velasco, I. Calvo, C. Silva, P. Helander, R. Kleiber

Abstract: 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 el… ▽ More 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. △ Less

Submitted 12 December, 2014; v1 submitted 3 April, 2014; originally announced April 2014.

Journal ref: 2015 Nucl. Fusion 55 052001

Stellarators close to quasisymmetry

Authors: Ivan Calvo, Felix I. Parra, J. L. Velasco, J. Arturo Alonso

Abstract: 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 + α\mathbf{B}_1$, where $\mathbf{B}_0$ is quasisymmetric, $α\mathbf{B}_1$ is a deviation from quasisymmetry, and $α\ll 1$, the stellarator can rotate at high v… ▽ More 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 + α\mathbf{B}_1$, where $\mathbf{B}_0$ is quasisymmetric, $α\mathbf{B}_1$ is a deviation from quasisymmetry, and $α\ll 1$, the stellarator can rotate at high velocities if $α< ε^{1/2}$, with $ε$ 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(ε^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(ε^2)$ pieces. △ Less

Submitted 6 February, 2014; v1 submitted 12 July, 2013; originally announced July 2013.

Comments: 59 pages, 2 figures. arXiv admin note: text overlap with arXiv:1204.1509

Journal ref: Plasma Physics and Controlled Fusion 55, 125014 (2013)

Damping of radial electric field fluctuations in the TJ-II stellarator

Authors: J L Velasco, J A Alonso, I Calvo, J Arévalo, E Sánchez, L Eliseev, S Perfilov, T Estrada, A López-Fraguas, C Hidalgo, the TJ-II team

Abstract: 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 fluctua… ▽ More 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. △ Less

Submitted 28 March, 2017; v1 submitted 5 July, 2013; originally announced July 2013.

Journal ref: Plasma Physics and Controlled Fusion, 55, 124044 (2013)

A DDS-Based Scalable and Reconfigurable Framework for Cyber-Physical Systems

Authors: Ismael Etxeberria-Agiriano, Isidro Calvo, Liliana Montero, Ivan Alonso

Abstract: Cyber-Physical Systems (CPSs) involve the interconnection of heterogeneous computing devices which are closely integrated with the physical processes under control. Often, these systems are resource-constrained and require specific features such as the ability to adapt in a timeliness and efficient fashion to dynamic environments. Also, they must support fault tolerance and avoid single points of… ▽ More Cyber-Physical Systems (CPSs) involve the interconnection of heterogeneous computing devices which are closely integrated with the physical processes under control. Often, these systems are resource-constrained and require specific features such as the ability to adapt in a timeliness and efficient fashion to dynamic environments. Also, they must support fault tolerance and avoid single points of failure. This paper describes a scalable framework for CPSs based on the OMG DDS standard. The proposed solution allows reconfiguring this kind of systems at run-time and managing efficiently their resources. △ Less

Submitted 10 August, 2012; originally announced August 2012.

Comments: 13 pages, 6 figures

Journal ref: International Journal of Software Engineering & Applications (IJSEA), vol. 3, no. 4 (2012) 25-37

Vanishing neoclassical viscosity and physics of the shear layer in stellarators

Authors: J. L. Velasco, J. A. Alonso, I. Calvo, J. Arevalo

Abstract: 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. T… ▽ More 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. △ Less

Submitted 16 August, 2012; v1 submitted 5 June, 2012; originally announced June 2012.

Comments: 5 pages, 5 figures, accepted for publication in Physical Review Letters

Journal ref: Physical Review Letters 109, 135003 (2012)

Long-wavelength limit of gyrokinetics in a turbulent tokamak and its intrinsic ambipolarity

Authors: Ivan Calvo, Felix I. Parra

Abstract: Recently, the electrostatic gyrokinetic Hamiltonian and change of coordinates have been computed to order $ε^2$ in general magnetic geometry. Here $ε$ is the gyrokinetic expansion parameter, the gyroradius over the macroscopic scale length. Starting from these results, the long-wavelength limit of the gyrokinetic Fokker-Planck and quasineutrality equations is taken for tokamak geometry. Employing… ▽ More Recently, the electrostatic gyrokinetic Hamiltonian and change of coordinates have been computed to order $ε^2$ in general magnetic geometry. Here $ε$ is the gyrokinetic expansion parameter, the gyroradius over the macroscopic scale length. Starting from these results, the long-wavelength limit of the gyrokinetic Fokker-Planck and quasineutrality equations is taken for tokamak geometry. Employing the set of equations derived in the present article, it is possible to calculate the long-wavelength components of the distribution functions and of the poloidal electric field to order $ε^2$. These higher-order pieces contain both neoclassical and turbulent contributions, and constitute one of the necessary ingredients (the other is given by the short-wavelength components up to second order) that will eventually enter a complete model for the radial transport of toroidal angular momentum in a tokamak in the low flow ordering. Finally, we provide an explicit and detailed proof that the system consisting of second-order gyrokinetic Fokker-Planck and quasineutrality equations leaves the long-wavelength radial electric field undetermined; that is, the turbulent tokamak is intrinsically ambipolar. △ Less

Submitted 23 July, 2014; v1 submitted 6 April, 2012; originally announced April 2012.

Comments: 70 pages. Typos in equations (63), (90), (91), (92) and (129) corrected

Journal ref: Plasma Physics and Controlled Fusion 54, 115007 (2012)

Intrinsic rotation with gyrokinetic models

Authors: Felix I. Parra, Michael Barnes, Ivan Calvo, Peter J. Catto

Abstract: The generation of intrinsic rotation by turbulence and neoclassical effects in tokamaks is considered. To obtain the complex dependences observed in experiments, it is necessary to have a model of the radial flux of momentum that redistributes the momentum within the tokamak in the absence of a preexisting velocity. When the lowest order gyrokinetic formulation is used, a symmetry of the model pre… ▽ More The generation of intrinsic rotation by turbulence and neoclassical effects in tokamaks is considered. To obtain the complex dependences observed in experiments, it is necessary to have a model of the radial flux of momentum that redistributes the momentum within the tokamak in the absence of a preexisting velocity. When the lowest order gyrokinetic formulation is used, a symmetry of the model precludes this possibility, making small effects in the gyroradius over scale length expansion necessary. These effects that are usually small become important for momentum transport because the symmetry of the lowest order gyrokinetic formulation leads to the cancellation of the lowest order momentum flux. The accuracy to which the gyrokinetic equation needs to be obtained to retain all the physically relevant effects is discussed. △ Less

Submitted 22 March, 2012; originally announced March 2012.

Journal ref: Phys. Plasmas 19, 056116 (2012)