Professor Felix Parra Diaz

Abstract:

In this comment, we show that by using the discrete particle distribution function the changes of the phase-space volume of gyrocenter coordinates due to the fluctuating E×B velocity do not explicitly appear in the Poisson equation and the [Sosenko, Phys. Scr. 64, 264 (2001)] result is recovered. It is demonstrated that there is no contradiction between the work presented by Sosenko and the work presented by [Wang, Phys. Plasmas 17, 082304 (2010)]. © 2010 American Institute of Physics.

Abstract:

The effect of flow shear on turbulent transport in tokamaks is studied numerically in the experimentally relevant limit of zero magnetic shear. It is found that the plasma is linearly stable for all nonzero flow shear values, but that subcritical turbulence can be sustained nonlinearly at a wide range of temperature gradients. Flow shear increases the nonlinear temperature gradient threshold for turbulence but also increases the sensitivity of the heat flux to changes in the temperature gradient, except over a small range near the threshold where the sensitivity is decreased. A bifurcation in the equilibrium gradients is found: for a given input of heat, it is possible, by varying the applied torque, to trigger a transition to significantly higher temperature and flow gradients. © 2010 The American Physical Society.

Abstract:

Gyrokinetic theory is based on an asymptotic expansion in the small parameter $\epsilon$, 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 $\epsilon^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 $\epsilon^2$ and a Poisson equation accurate to order $\epsilon$. The final expressions are explicit and can be implemented into any simulation without further computations.

Abstract:

Quasineutral particle-in-cell models of ions must fulfill the kinetic Bohm criterion, in its inequality form, at the domain boundary in order to match correctly with solutions of the Debye sheaths tied to the walls. The simple, fluid form of the Bohm criterion is shown to be a bad approximation of the exact, kinetic form when the ion velocity distribution function has a significant dispersion and involves different charge numbers. The fulfillment of the Bohm criterion is measured by a weighting algorithm at the boundary, but linear weighting algorithms have difficulties to reproduce the nonlinear behavior around the sheath edge. A surface weighting algorithm with an extended temporal weighting is proposed and shown to behave better than the standard volumetric weighting. Still, this must be supplemented by a forcing algorithm of the kinetic Bohm criterion. This postulates a small potential fall in a supplementary, thin, transition layer. The electron-wall interaction is shown to be of little relevance in the fulfillment of the Bohm criterion. © 2010 American Institute of Physics.

Abstract:

Full f electrostatic gyrokinetic formulations employ two gyrokinetic equations, one for ions and the other for electrons, and quasineutrality to obtain the ion and electron distribution functions and the electrostatic potential. We demonstrate with several examples that the long wavelength radial electric field obtained with full f approaches is extremely sensitive to errors in the ion and electron density since small deviations in density give rise to large, nonphysical deviations in the conservation of toroidal angular momentum. For typical tokamak values, a relative error of 10-7 in the ion or electron densities is enough to obtain the incorrect toroidal rotation. Based on the insights gained with the examples considered, three simple tests to check transport of toroidal angular momentum in full f simulations are proposed. © 2010 American Institute of Physics.