Professor Felix Parra Diaz

Abstract:

Refer to readme.pdf in the repository.

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 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.

Abstract:

The magnetic presheath is a boundary layer occurring when magnetized plasma is in contact with a wall and the angle α between the wall and the magnetic field B is oblique. Here, we consider the fusion-relevant case of a shallow-angle, α 1, electron-repelling sheath, with the electron density given by a Boltzmann distribution, valid for α/√τ + 1 √me/mi, where me is the electron mass, mi is the ion mass, τ = Ti/ZTe,Te is the electron temperature, Ti is the ion temperature and Z is the ionic charge state. The thickness of the magnetic presheath is of the order of a few ion sound Larmor radii ρs = √mi(ZTe + Ti)/ZeB, where e is the proton charge and B = |B| is the magnitude of the magnetic field. We study the dependence on τ of the electrostatic potential and ion distribution function in the magnetic presheath by using a set of prescribed ion distribution functions at the magnetic presheath entrance, parameterized by τ . The kinetic model is shown to be asymptotically equivalent to Chodura’s fluid model at small ion temperature, τ 1, for |ln α| > 3|ln τ | 1. In this limit, despite the fact that fluid equations give a reasonable approximation to the potential, ion gyroorbits acquire a spatial extent that occupies a large portion of the magnetic presheath. At large ion temperature, τ 1, relevant because Ti is measured to be a few times larger than Te near divertor targets of fusion devices, ions reach the Debye sheath entrance (and subsequently the wall) at a shallow angle whose size is given by √α or 1/√τ, depending on which is largest.

Abstract:

Here we present details of an operator-split, implicit–explicit numerical scheme for the solution of the gyrokinetic-Poisson system of equations in the local limit. This scheme has been implemented in a new code called stella, which is capable of evolving electrostatic fluctuations with full kinetic electron effects and an arbitrary number of ion species in general magnetic geometry. We demonstrate the advantages of this mixed approach over a fully explicit treatment and provide linear and nonlinear benchmark comparisons for both axisymmetric and non-axisymmetric magnetic equilibria.