Related papers: ORB5: a global electromagnetic gyrokinetic code us…
Modeling multi-scale collisionless magnetized processes constitutes an important numerical challenge. By treating electrons as a fluid and ions kinetically, the so-called hybrid Particle-In-Cell (PIC) codes represent a promising…
We develop a gyrokinetic treatment for ions in the magnetic presheath, close to the plasma-wall boundary. We focus on magnetic presheaths with a small magnetic field to wall angle, $\alpha \ll 1$ (in radians). Characteristic lengths…
We describe a novel Godunov-type numerical method for solving the equations of resistive relativistic magnetohydrodynamics. In the proposed approach, the spatial components of both magnetic and electric fields are located at zone interfaces…
High-temperature superconducting (HTS) magnets and other advances have led to renewed interest in magnetic mirrors for fusion energy. The non-Maxwellian nature of mirror plasmas necessitates kinetic modeling to predict, optimize and design…
Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in…
We report on the algorithms and numerical methods used in Viriato, a novel fluid-kinetic code that solves two distinct sets of equations: (i) the Kinetic Reduced Electron Heating Model (KREHM) equations [Zocco & Schekochihin, Phys. Plasmas…
We have performed fully-kinetic simulations of X-B and O-X-B mode conversion in one and two dimensional setups using the PIC code EPOCH. We have recovered the linear dispersion relation for electron Bernstein waves by employing relatively…
An analytic equilibrium, the Toroidal Bessel Function Model, is used in conjunction with the gyrokinetic code GYRO to investigate the nature of microinstabilities in a reversed field pinch (RFP) plasma. The effect of the normalized electron…
In this work, we present the first global gyrokinetic simulations of the ITER baseline scenario operating at 15 MA using GENE-Tango electrostatic and electromagnetic simulations. The modeled radial region spans close to the magnetic axis up…
A new mechanism for toroidal momentum transport in a tokamak is investigated using the gyro-kinetic model. First, an analytic model is developed through the use of the ballooning transform. The terms that generate the momentum transport are…
This paper presents a hybrid numerical method for linear collisional kinetic equations with diffusive scaling. The aim of the method is to reduce the computational cost of kinetic equations by taking advantage of the lower dimensionality of…
For the simulation of disruptions in tokamak fusion plasmas, a fluid model describing the evolution of relativistic runaway electrons and their interaction with the background plasma is presented. The overall aim of the model is to…
Ion orbit loss has been used to model the formation of a strong negative radial electric field $E_r$ in the tokamak edge, as well as edge momentum transport and toroidal rotation. To quantitatively measure ion orbit loss, an orbit-flux…
Ion-gyroradius-scale microinstabilities typically have a frequency comparable to the ion transit frequency. Due to the small electron-to-ion mass ratio and the large electron transit frequency, it is conventionally assumed that passing…
Turbulent dynamics in the scrape-off layer (SOL) of magnetic fusion devices is intermittent with large fluctuations in density and pressure. Therefore, a model is required that allows perturbations of similar or even larger magnitude to the…
Long global gyrokinetic turbulence simulations are particularly challenging in situations where the system deviates strongly from its initial state and when fluctuation level are high e.g. in strong gradient regions. For Particle-in-Cell…
Linear kinetic Monte Carlo particle transport models are frequently employed in fusion plasma simulations to quantify atomic and surface effects on the main plasma flow dynamics. Separate codes are used for transport of neutral particles…
A Continuous Galerkin method-based approach is presented to compute the seismic normal modes of rotating planets. Special care is taken to separate out the essential spectrum in the presence of a fluid outer core using a polynomial…
In this paper we present a newly developed GRPIC code framework called FPIC, providing a detailed description of the Maxwell-equations solver, of the particle ``pushers'', and of the other algorithms that are needed in this approach. We…
The Eulerian variational formulation is presented to obtain governing equations of the electromagnetic turbulent gyrokinetic system. A local momentum balance in the system is derived from the invariance of the Lagrangian of the system under…