Search

Piecewise omnigenous fields are stellarator magnetic fields that are optimized with respect to radial neoclassical transport thanks to a second adiabatic invariant that is piecewisely constant on the flux-surface. They are qualitatively different from omnigenous fields (including quasi-isodynamic or quasisymmetric fields), for which the second adiabatic invariant is a flux-surface constant. Piecewise omnigenous fields thus open an alternative path towards stellarator reactors. In this work, piecewise omnigenous fields are characterized and parametrized in a systematic manner. This is a step towards including piecewise omnigenity as an explicit design criterion in stellarator optimization, and towards a systematic study of the properties of nearly piecewise omnigenous stellarator configurations.

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

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

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

People with diabetes must carefully monitor their blood glucose levels, especially after eating. Blood glucose regulation requires a proper combination of food intake and insulin boluses. Glucose prediction is vital to avoid dangerous post-meal complications in treating individuals with diabetes. Although traditional methods, such as artificial neural networks, have shown high accuracy rates, sometimes they are not suitable for developing personalised treatments by physicians due to their lack of interpretability. In this study, we propose a novel glucose prediction method emphasising interpretability: Interpretable Sparse Identification by Grammatical Evolution. Combined with a previous clustering stage, our approach provides finite difference equations to predict postprandial glucose levels up to two hours after meals. We divide the dataset into four-hour segments and perform clustering based on blood glucose values for the twohour window before the meal. Prediction models are trained for each cluster for the two-hour windows after meals, allowing predictions in 15-minute steps, yielding up to eight predictions at different time horizons. Prediction safety was evaluated based on Parkes Error Grid regions. Our technique produces safe predictions through explainable expressions, avoiding zones D (0.2% average) and E (0%) and reducing predictions on zone C (6.2%). In addition, our proposal has slightly better accuracy than other techniques, including sparse identification of non-linear dynamics and artificial neural networks. The results demonstrate that our proposal provides interpretable solutions without sacrificing prediction accuracy, offering a promising approach to glucose prediction in diabetes management that balances accuracy, interpretability, and computational efficiency.

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

People with diabetes have to manage their blood glucose level to keep it within an appropriate range. Predicting whether future glucose values will be outside the healthy threshold is of vital importance in order to take corrective actions to avoid potential health damage. In this paper we describe our research with the aim of predicting the future behavior of blood glucose levels, so that hypoglycemic events may be anticipated. The approach of this work is the application of transformation functions on glucose time series, and their use in convolutional neural networks. We have tested our proposed method using real data from 4 different diabetes patients with promising results., Comment: 7 pages, 7 figures, 3 tables

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

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

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

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

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

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

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

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

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

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

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

This paper presents a study conducted to allow urban farmers to remotely monitor their farm through the design and development of an Internet of Things-based (IoT) microfarm prototype which utilized wick system as planting method. The system involves the detection of three environmental parameters namely, light intensity, soil moisture and temperature through the use of respective sensors which were connected to the Arduino microcontroller, the sensor node of the system. Irregularities in the aforementioned parameters were neutralized through the use of parameter regulators such as LED growlight strips, water pump and air cooler. The data collected by these sensors were gathered by the Arduino microcontroller and were sent to the Web database through the IoT gateway which was the Raspberry Pi computer chip. These data were also sent to an Android unit installed with the Microfarm Companion application which was capable of monitoring and controlling the environmental parameters observed in the microfarm. The application allows the user to view the current value of the parameter involved and to choose whether to control the parameter regulators automatically or manually. The microfarm system runs autonomously which reduces the labor required to produce healthy plants and crops. Mustard greens samples were used in testing the system. After a month of monitoring the height of the samples, it was observed that the average height of the samples is about 0.23 cm taller than the standard height. The proponents has also tested the system functionality by evaluating the sensor data log that provides the values gathered by the sensors and the turn-on times of the parameter regulators. From these data, it can be observed that whenever the values obtained by the sensors fall outside the threshold range, the parameter regulators turns on, indicating that the system is working properly.

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

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

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

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

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

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

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