Simulations of microturbulence in magnetised plasmas using a delta-f gyrokinetic approach with an evolving background Maxwellian
Open Access
- 1 February 2021
- journal article
- research article
- Published by IOP Publishing in Journal of Physics: Conference Series
- Vol. 1785 (1), 012003
- https://doi.org/10.1088/1742-6596/1785/1/012003
Abstract
The gyrokinetic delta-f particle-in-cell (PIC) approach is known to be successful for simulating turbulence in the core of magnetic fusion plasmas, where fluctuations are relatively small and therefore the unperturbed particle distribution function, usually represented by a stationary Maxwellian f0, remains a good choice of a control variate for reducing statistical sampling noise. However, towards the plasma edge, characterized by low density and temperature and strong gradients, relative deviation amplitudes typically become large, so that the essential assumption of |δf/f0| << 1 underlying the delta-f PIC approach will not be valid, where δf is the fluctuating part of distribution. This motivates the study of the limits of the delta-f approach in a simplified system mimicking the plasma edge. To this end, simulations are run using GK-engine, which is a delta-f PIC code that solves the nonlinear gyrokinetic equation in a sheared slab geometry, using B-spline finite elements to represent the self-consistent electrostatic field. Initial radial density and ion temperature profiles exhibiting high logarithmic gradients representing plasma edge conditions are used. In order to avoid practical problems of particles exiting the simulation domain as the ion temperature profile relaxes, all profiles are mirrored at domain-centre and periodic boundary conditions are imposed. The validity of the delta-f approach is measured by statistical noise estimates, while monitoring relative deviation levels of temperature via the kinetic energy. In particular, the effect of background profile gradients on these measures is investigated. As a first step towards reducing the amplitude of the deviation δf, an adaptive Maxwellian f0 is implemented, whose time dependent temperature profiles are obtained by locally relaxing kinetic energy accumulating in δf into f0.This publication has 10 references indexed in Scilit:
- Global turbulence features across marginality and non-local pedestal-core interactionsPlasma Physics and Controlled Fusion, 2018
- A new hybrid-Lagrangian numerical scheme for gyrokinetic simulation of tokamak edge plasmaJournal of Computational Physics, 2016
- Monte Carlo particle-in-cell methods for the simulation of the Vlasov–Maxwell gyrokinetic equationsJournal of Plasma Physics, 2015
- Gyrokinetic simulations of turbulent transportNuclear Fusion, 2010
- Long global gyrokinetic simulations: Source terms and particle noise controlPhysics of Plasmas, 2008
- A guiding-center Fokker–Planck collision operator for nonuniform magnetic fieldsPhysics of Plasmas, 2004
- A revised δf algorithm for nonlinear PIC simulationComputer Physics Communications, 2003
- Collisional delta-f scheme with evolving background for transport time scale simulationsPhysics of Plasmas, 1999
- A unified Monte Carlo interpretation of particle simulations and applications to non-neutral plasmasPhysics of Plasmas, 1994
- Nonlinear gyrokinetic equations for tokamak microturbulencePhysics of Fluids, 1988