NORSE: A solver for the relativistic non-linear Fokker-Planck equation for electrons in a homogeneous plasma

详细信息    查看全文

• 作者：A. Stahl^a ; ^{stahla@chalmers.se} ; M. Landreman^b ; O. Embré ; us^a ; T. Fü ; lö ; p^a
• 关键词：Non-linear relativistic Fokker&ndash ; Planck equation ; Kinetic plasma theory ; Energetic electrons ; Runaway electrons
• 刊名：Computer Physics Communications
• 出版年：2017
• 出版时间：March 2017
• 年：2017
• 卷：212
• 期：Complete
• 页码：269-279
• 全文大小：819 K
• 卷排序：212

文摘

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 summaryProgram title: NORSEProgram Files doi:http://dx.doi.org/10.17632/86wmgj758w.1Licensing 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.