BZPEERPapers, people, evidence and research paths
Menu

Discover

SearchFind papers, people, topics, institutions and opportunities from one place.GraphSee how papers, authors, topics and evidence connect.

Research tools

HomeContinue from the papers, people and topics that matter now.CollectionsTurn saved papers into reading lists, evidence maps and shareable context.ResearchReconstruct the path that led you from search to paper to question.CollaborationMove from a paper, topic or expert request into a focused thread.PublishCreate full-text research articles with metadata, formulas and versions.

Network

ExpertsFind specialists whose work and reviews explain why they can help.CommunitiesJoin research groups organized around papers, topics and shared evidence.

Opportunities

ListingsBrowse roles, calls and collaborations connected to real research context.

Account

Log inReturn to your saved papers, graph views and active threads.Create accountStart saving papers, building a research trail and joining teams.
Log inCreate account

Source author record

Felix I. Parra

Felix I. Parra appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record
physics.plasm-phphysics.acc-ph

Catalog footprint

What is connected

23works
2topics
4close collaborators

Actions

Connect this record

Log in to claim
Open graphBrowse works

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

23 published item(s)

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.

0 citations0 reviewsNo code linkedOpen access
preprint2018arXiv

Solution to a collisionless shallow-angle magnetic presheath with kinetic ions

Using a kinetic model for the ions and adiabatic electrons, we solve a steady state, electron-repelling magnetic presheath in which a uniform magnetic field makes a small angle $α\ll 1$ (in radians) with the wall. The presheath characteristic thickness is the typical ion gyroradius $ρ_{\text{i}}$. The Debye length $λ_{\text{D}}$ and the collisional mean free path of an ion $λ_{\text{mfp}}$ satisfy the ordering $λ_{\text{D}} \ll ρ_{\text{i}} \ll αλ_{\text{mfp}}$, so a quasineutral and collisionless model is used. We assume that the electrostatic potential is a function only of distance from the wall, and it varies over the scale $ρ_{\text{i}}$. Using the expansion in $α\ll 1$, we derive an analytical expression for the ion density that only depends on the ion distribution function at the entrance of the magnetic presheath and the electrostatic potential profile. Importantly, we have added the crucial contribution of the orbits in the region near the wall. By imposing the quasineutrality equation, we derive a condition that the ion distribution function must satisfy at the magnetic presheath entrance --- the kinetic equivalent of the Chodura condition. Using an ion distribution function at the entrance of the magnetic presheath that satisfies the kinetic Chodura condition, we find numerical solutions for the self-consistent electrostatic potential, ion density and flow across the magnetic presheath for several values of $α$. Our numerical results also include the distribution of ion velocities at the Debye sheath entrance. We find that at small values of $α$ there are substantially fewer ions travelling with a large normal component of the velocity into the wall.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Effect of the Shafranov shift and the gradient of $β$ on intrinsic momentum transport in up-down asymmetric tokamaks

Tokamaks with up-down asymmetric poloidal cross-sections spontaneously rotate due to turbulent transport of momentum. In this work, we investigate the effect of the Shafranov shift on this intrinsic rotation, primarily by analyzing tokamaks with tilted elliptical flux surfaces. By expanding the Grad-Shafranov equation in the large aspect ratio limit we calculate the magnitude and direction of the Shafranov shift in tilted elliptical tokamaks. The results show that, while the Shafranov shift becomes up-down asymmetric and depends strongly on the tilt angle of the flux surfaces, it is insensitive to the shape of the current and pressure profiles (when the geometry, total plasma current, and average pressure gradient are kept fixed). Next, local nonlinear gyrokinetic simulations of these MHD equilibria are performed with GS2, which reveal that the Shafranov shift can significantly enhance the momentum transport. However, to be consistent, the effect of $β'$ (i.e. the radial gradient of $β$) on the magnetic equilibrium was also included, which was found to significantly reduce momentum transport. Including these two competing effects broadens the rotation profile, but leaves the on-axis value of the rotation roughly unchanged. Consequently, the shape of the $β$ profile has a significant effect on the rotation profile of an up-down asymmetric tokamak.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Poloidal tilting symmetry of high order tokamak flux surface shaping in gyrokinetics

A poloidal tilting symmetry of the local nonlinear $δf$ gyrokinetic model is demonstrated analytically and verified numerically. This symmetry shows that poloidally rotating all the flux surface shaping effects with large poloidal mode number by a single tilt angle has an exponentially small effect on the transport properties of a tokamak. This is shown using a generalization of the Miller local equilibrium model to specify an arbitrary flux surface geometry. With this geometry specification we find that, when performing an expansion in large flux surface shaping mode number, the governing equations of gyrokinetics are symmetric in the poloidal tilt of the high order shaping effects. This allows us to take the fluxes from a single configuration and calculate the fluxes in any configuration that can be produced by tilting the large mode number shaping effects. This creates a distinction between tokamaks with mirror symmetric flux surfaces and tokamaks without mirror symmetry, which is expected to have important consequences for generating toroidal rotation using up-down asymmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2016arXiv

Scaling of up-down asymmetric turbulent momentum flux with poloidal shaping mode number in tokamaks

Breaking the up-down symmetry of tokamaks removes a constraint limiting intrinsic momentum transport, and hence toroidal rotation, to be small. Using gyrokinetic theory, we study the effect of different up-down asymmetric flux surface shapes on the turbulent transport of momentum. This is done by perturbatively expanding the gyrokinetic equation in large flux surface shaping mode number. It is found that the momentum flux generated by shaping that lacks mirror symmetry (which is necessarily up-down asymmetric) has a power law scaling with the shaping mode number. However, the momentum flux generated by mirror symmetric flux surface shaping (even if it is up-down asymmetric) decays exponentially with large shaping mode number. These scalings are consistent with nonlinear local gyrokinetic simulations and indicate that low mode number shaping effects (e.g. elongation, triangularity) are optimal for creating rotation. Additionally it suggests that breaking the mirror symmetry of flux surfaces may generate significantly more toroidal rotation.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Equivalence of two different approaches to global $δf$ gyrokinetic simulations

A set of flux tube gyrokinetic equations that includes the effect of the spatial variation of the density, temperature and rotation gradients on the turbulence is derived. This new set of equations uses periodic boundary conditions. In the limit $l_\bot/L \ll 1$, where $l_\bot$ is the characteristic perpendicular length of turbulent structures and $L$ is the characteristic size of the device, this new set of flux tube gyrokinetic equations is shown to be equivalent to the traditional global $δf$ gyrokinetic equations to an order higher in $l_\bot/L$ than the usual flux tube formulations.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Intrinsic rotation in tokamaks: theory

Self-consistent equations for intrinsic rotation in tokamaks with small poloidal magnetic field $B_p$ compared to the total magnetic field $B$ are derived. The model gives the momentum redistribution due to turbulence, collisional transport and energy injection. Intrinsic rotation is determined by the balance between the momentum redistribution and the turbulent diffusion and convection. Two different turbulence regimes are considered: turbulence with characteristic perpendicular lengths of the order of the ion gyroradius, $ρ_i$, and turbulence with characteristic lengths of the order of the poloidal gyroradius, $(B/B_p) ρ_i$. Intrinsic rotation driven by gyroradius scale turbulence is mainly due to the effect of neoclassical corrections and of finite orbit widths on turbulent momentum transport, whereas for the intrinsic rotation driven by poloidal gyroradius scale turbulence, the slow variation of turbulence characteristics in the radial and poloidal directions and the turbulent particle acceleration can be become as important as the neoclassical and finite orbit width effects. The magnetic drift is shown to be indispensable for the intrinsic rotation driven by the slow variation of turbulence characteristics and the turbulent particle acceleration. The equations are written in a form conducive to implementation in a flux tube code, and the effect of the radial variation of the turbulence is included in a novel way that does not require a global gyrokinetic formalism.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Less constrained omnigeneous stellarators

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.

0 citations0 reviewsNo code linkedOpen access
preprint2015arXiv

Radial penetration of flux surface shaping in tokamaks

Using analytic calculations, the effects of the edge flux surface shape and the toroidal current profile on the penetration of flux surface shaping are investigated in a tokamak. It is shown that the penetration of shaping is determined by the poloidal variation of the poloidal magnetic field on the surface. This fact is used to investigate how different flux surface shapes penetrate from the edge. Then, a technique to separate the effects of magnetic pressure and tension in the Grad-Shafranov equation is presented and used to calculate radial profiles of strong elongation for nearly constant current profiles. Lastly, it is shown that more hollow toroidal current profiles are significantly better at conveying shaping from the edge to the core.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

A current driven electromagnetic mode in sheared and toroidal configurations

The induced electric field in a tokamak drives a parallel electron current flow. In an inhomogeneous, finite beta plasma, when this electron flow is comparable to the ion thermal speed, the Alfven mode wave solutions of the electromagnetic gyrokinetic equation can become nearly purely growing kink modes. Using the new "low-flow" version of the gyrokinetic code GS2 developed for momentum transport studies [Barnes et al 2013 Phys. Rev. Lett. 111, 055005], we are able to model the effect of the induced parallel electric field on the electron distribution to study the destabilizing influence of current on stability. We identify high mode number kink modes in GS2 simulations and make comparisons to analytical theory in sheared magnetic geometry. We demonstrate reassuring agreement with analytical results both in terms of parametric dependences of mode frequencies and growth rates, and regarding the radial mode structure.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Intrinsic momentum transport in up-down asymmetric tokamaks

Recent work demonstrated that breaking the up-down symmetry of tokamak flux surfaces removes a constraint that limits intrinsic momentum transport, and hence toroidal rotation, to be small. We show, through MHD analysis, that ellipticity is most effective at introducing up-down asymmetry throughout the plasma. We detail an extension to GS2, a local $δf$ gyrokinetic code that self-consistently calculates momentum transport, to permit up-down asymmetric configurations. Tokamaks with tilted elliptical poloidal cross-sections were simulated to determine nonlinear momentum transport. The results, which are consistent with experiment in magnitude, suggest that a toroidal velocity gradient, $(\partial u_{ζi} / \partial ρ) / v_{th i}$, of 5% of the temperature gradient, $(\partial T_{i} / \partial ρ) / T_{i}$, is sustainable. Here $v_{th i}$ is the ion thermal speed, $u_{ζi}$ is the ion toroidal mean flow, $ρ$ is the minor radial coordinate normalized to the tokamak minor radius, and $T_{i}$ is the ion temperature. Since other intrinsic momentum transport mechanisms scale poorly to larger machines, these results indicate that up-down asymmetry is the most feasible method to generate the current experimentally-measured rotation levels in reactor-sized devices.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

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

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.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

Phase-space Lagrangian derivation of electrostatic gyrokinetics in general geometry

Gyrokinetic theory is based on an asymptotic expansion in the small parameter $ε$, defined as the ratio of the gyroradius and the characteristic length of variation of the magnetic field. In this article, this ordering is strictly implemented to compute the electrostatic gyrokinetic phase-space Lagrangian in general magnetic geometry to order $ε^2$. In particular, a new expression for the complete second-order gyrokinetic Hamiltonian is provided, showing that in a rigorous treatment of gyrokinetic theory magnetic geometry and turbulence cannot be dealt with independently. The new phase-space gyrokinetic Lagrangian gives a Vlasov equation accurate to order $ε^2$ and a Poisson equation accurate to order $ε$. The final expressions are explicit and can be implemented into any simulation without further computations.

0 citations0 reviewsNo code linkedOpen access
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.

0 citations0 reviewsNo code linkedOpen access
preprint2014arXiv

When omnigeneity fails

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.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Conditions for up-down asymmetry in the core of tokamak equilibria

A local magnetic equilibrium solution is sought around the magnetic axis in order to identify the key parameters defining the magnetic-surface's up-down asymmetry in the core of tokamak plasmas. The asymmetry is found to be determined essentially by the ratio of the toroidal current density flowing on axis to the fraction of the external field's odd perturbation that manages to propagate from the plasma boundary into the core. The predictions are tested and illustrated first with an analytical Solovev equilibrium and then using experimentally relevant numerical equilibria. Hollow current-density distributions, and hence reverse magnetic shear, are seen to be crucial to bring into the core asymmetry values that are usually found only near the plasma edge.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Extension of gyrokinetics to transport time scales

Gyrokinetic simulations have greatly improved our theoretical understanding of turbulent transport in fusion devices. Most gyrokinetic models in use are delta-f simulations in which the slowly varying radial profiles of density and temperature are assumed to be constant for turbulence saturation times, and only the turbulent electromagnetic fluctuations are calculated. New massive simulations are being built to self-consistently determine the radial profiles of density and temperature. However, these new codes have failed to realize that modern gyrokinetic formulations, composed of a gyrokinetic Fokker-Planck equation and a gyrokinetic quasineutrality equation, are only valid for delta-f simulations that do not reach the longer transport time scales necessary to evolve radial profiles. In tokamaks, due to axisymmetry, the evolution of the axisymmetric radial electric field is a challenging problem requiring substantial modifications to gyrokinetic treatments. In this thesis, I study the effect of turbulence on the global electric field and plasma flows. By studying the current conservation equation, or vorticity equation, I prove that the long wavelength, axisymmetric flow must remain neoclassical and I show that the tokamak is intrinsically ambipolar, i.e., the radial current is zero to a very high order for any long wavelength radial electric field. Intrinsic ambipolarity is the origin of the problems with the modern gyrokinetic approach since the lower order gyrokinetic quasineutrality (if properly evaluated) is effectively independent of the radial electric field. I propose a new gyrokinetic formalism to solve for the global radial electric field.

0 citations0 reviewsNo code linkedOpen access
preprint2013arXiv

Turbulent momentum pinch of diamagnetic flows in a tokamak

The ion toroidal rotation in a tokamak consists of an $E\times B$ flow due to the radial electric field and a diamagnetic flow due to the radial pressure gradient. The turbulent pinch of toroidal angular momentum due to the Coriolis force studied in previous work is only applicable to the $E\times B$ flow. In this Letter, the momentum pinch for the rotation generated by the radial pressure gradient is calculated and is compared with the Coriolis pinch. This distinction is important for subsonic flows or the flow in the pedestal where the two types of flows are similar in size and opposite in direction. In the edge, the different pinches due to the opposite rotations can result in intrinsic momentum transport that gives significant rotation peaking.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Analytic fluid theory of beam spiraling in high-intensity cyclotrons

Using a two-dimensional fluid description, we investigate the nonlinear radial-longitudinal dynamics of intense beams in storage rings and cyclotrons. With a multiscale analysis separating the time scale associated with the betatron motion and the slower time scale associated with space-charge effects, we show that the longitudinal-radial vortex motion can be understood in the frame moving with the charged beam as the nonlinear advection of the beam by the $\mathbf{E}\times\mathbf{B}$ velocity field, where $\mathbf{E}$ is the electric field due to the space charge and $\mathbf{B}$ is the external magnetic field. This interpretation provides simple explanations for the stability of round beams and for the development of spiral halos in elongated beams. By numerically solving the nonlinear advection equation for the beam density, we find that it is also in quantitative agreement with results obtained in PIC simulations.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Perpendicular momentum injection by lower hybrid wave in a tokamak

The injection of lower hybrid waves for current drive into a tokamak affects the profile of intrinsic rotation. In this article, the momentum deposition by the lower hybrid wave on the electrons is studied. Due to the increase in the poloidal momentum of the wave as it propagates into the tokamak, the parallel momentum of the wave increases considerably. The change of the perpendicular momentum of the wave is such that the toroidal angular momentum of the wave is conserved. If the perpendicular momentum transfer via electron Landau damping is ignored, the transfer of the toroidal angular momentum to the plasma will be larger than the injected toroidal angular momentum. A proper quasilinear treatment proves that both perpendicular and parallel momentum are transferred to the electrons. The toroidal angular momentum of the electrons is then transferred to the ions via different mechanisms for the parallel and perpendicular momentum. The perpendicular momentum is transferred to ions through an outward radial electron pinch, while the parallel momentum is transferred through collisions.

0 citations0 reviewsNo code linkedOpen access
preprint2012arXiv

Sources of intrinsic rotation in the low flow ordering

A low flow, $δf$ gyrokinetic formulation to obtain the intrinsic rotation profiles is presented. The momentum conservation equation in the low flow ordering contains new terms, neglected in previous first principles formulations, that may explain the intrinsic rotation observed in tokamaks in the absence of external sources of momentum. The intrinsic rotation profile depends on the density and temperature profiles and on the up-down asymmetry.

0 citations0 reviewsNo code linkedOpen access
preprint2011arXiv

Up-down symmetry of the turbulent transport of toroidal angular momentum in tokamaks

Two symmetries of the local nonlinear delta-f gyrokinetic system of equations in tokamaks in the high flow regime are presented. The turbulent transport of toroidal angular momentum changes sign under an up-down reflection of the tokamak and a sign change of both the rotation and the rotation shear. Thus, the turbulent transport of toroidal angular momentum must vanish for up-down symmetric tokamaks in the absence of both rotation and rotation shear. This has important implications for the modeling of spontaneous rotation.

0 citations0 reviewsNo code linkedOpen access