Source author record

J. L. Velasco

J. L. Velasco appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

14works
1topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

14 published item(s)

preprint2021arXiv

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

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.

preprint2021arXiv

Physics design point of high-field stellarator reactors

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

preprint2020arXiv

Characterization of the radial electric field and edge velocity shear in Wendelstein 7-X

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 $ι$, 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.

preprint2020arXiv

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

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, φ_1. Furthermore, the equation solved is linear in φ_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.

preprint2020arXiv

Turbulent impurity transport simulations in Wendelstein 7-X plasmas

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.

preprint2019arXiv

Impact of main ion pressure anisotropy on stellarator impurity transport

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.

preprint2015arXiv

Electrostatic potential variation on the flux surface and its impact on impurity transport

The particle transport of impurities in magnetically confined plasmas under some conditions does not find, neither quantitatively nor qualitatively, a satisfactory theory-based explanation. This compromise the successful realization of thermo-nuclear fusion for energy production since its accumulation is known to be one of the causes that leads to the plasma breakdown. In standard reactor-relevant conditions this accumulation is in most stellarators intrinsic to the lack of toroidal symmetry, that leads to the neoclassical electric field to point radially inwards. This statement, that the standard theory allows to formulate, has been contradicted by some experiments that showed weaker or no accumulation under such conditions \cite{Ida_pop_16_056111_2009, Yoshinuma_nf_49_062002_2009}. The charge state of the impurities makes its transport more sensitive to the electric fields. Thus, the short length scale turbulent electrostatic potential or its long wave-length variation on the flux surface $Φ_{1}$ -- that the standard neoclassical approach usually neglects -- might possibly shed some light on the experimental findings. In the present work the focus is put on the second of the two, and investigate its influence of the radial transport of C$^{6+}$. We show that in LHD it is strongly modified by $Φ_{1}$, both resulting in mitigated/enhanced accumulation at internal/external radial positions; for Wendelstein 7-X, on the contrary, $Φ_{1}$ is expected to be considerably smaller and the transport of C$^{6+}$ not affected up to an appreciable extent; and in TJ-II the potential shows a moderate impact despite of the large amplitude of $Φ_1$ for the parameters considered.

preprint2014arXiv

Electrostatic potential variations along flux surfaces in stellarators

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.

preprint2014arXiv

Optimizing stellarators for large flows

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.

preprint2014arXiv

Stellarators close to quasisymmetry

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.

preprint2012arXiv

Dynamics of zonal flow-like structures in the edge of the TJ-II stellarator

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

preprint2012arXiv

Vanishing neoclassical viscosity and physics of the shear layer in stellarators

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.

preprint2011arXiv

Calculation of the bootstrap current profile for the TJ-II stellarator

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