FPPAC: A two-dimensional multispecies nonlinear Fokker-Planck package

Abstract

Energetic electrons are of interest in many types of plasmas, however previous modeling of their properties has been restricted to the use of linear Fokker–Planck collision operators or non-relativistic formulations. Here, we describe a fully non-linear kinetic-equation solver, capable of handling large electric-field strengths (compared to the Dreicer field) and relativistic temperatures. This tool allows modeling of the momentum-space dynamics of the electrons in cases where strong departures from Maxwellian distributions may arise. As an example, we consider electron runaway in magnetic-confinement fusion plasmas and describe a transition to electron slide-away at field strengths significantly lower than previously predicted.Program title: NORSEProgram Files doi: http://dx.doi.org/10.17632/86wmgj758w.1 Licensing provisions: GPLv3Programming language: MatlabNature of problem: Solves the Fokker–Planck equation for electrons in 2D momentum space in a homogeneous plasma (allowing for magnetization), using a relativistic non-linear electron–electron collision operator. Electric-field acceleration, synchrotron-radiation-reaction losses, as well as heat and particle sources are included. Scenarios with time-dependent plasma parameters can be studied.Solution method: The kinetic equation is represented on a non-uniform 2D finite-difference grid and is evolved using a linearly implicit time-advancement scheme. A mixed finite-difference-Legendre-mode representation is used to obtain the relativistic potentials (analogous to the non-relativistic Rosenbluth potentials) from the distribution.