TY - JOUR

T1 - Energy and charge conserving semi-implicit particle-in-cell model for simulations of high-pressure plasmas in magnetic traps

AU - Berendeev, E. A.

AU - Timofeev, I. V.

AU - Kurshakov, V. A.

N1 - The work is supported by Russian Science Foundation (grant N0 21-72-10071 ).

PY - 2024/2

Y1 - 2024/2

N2 - The paper presents a new particle-in-cell model for fully kinetic simulations of plasma confinement regimes with high relative pressure. These regimes are considered as the most promising ones for fusion reactor proposals based on field reversed configurations, mirror and multi-cusp magnetic traps. Formation of an equilibrium high-pressure plasma with a fully excluded or reversed magnetic field cannot be investigated using the simplified magnetohydrodynamic and gyrokinetic approaches. A correct description of the electron dynamics under close proximity of zero and strong magnetic field regions can only be achieved within the framework of the kinetic theory that can be implemented most efficiently by the particle-in-cell method. The full-scale particle-in-cell simulations of modern fusion experiments require the use of finite-difference schemes in which temporal steps exceed the period of fastest electron oscillations at the plasma or even cyclotron frequencies and spatial steps do not have to resolve the Debye radius. This requirement is satisfied by implicit schemes capable of conserving the total energy of the system. The particle-in-cell model presented in this paper is based on the energy conserving semi-implicit approach as a predictive step and a new method for suppressing the electrostatic noise inherent in it as a corrective step. This two-step algorithm provides not only accurate energy conservation, but also the exact local fulfillment of the Gauss law.

AB - The paper presents a new particle-in-cell model for fully kinetic simulations of plasma confinement regimes with high relative pressure. These regimes are considered as the most promising ones for fusion reactor proposals based on field reversed configurations, mirror and multi-cusp magnetic traps. Formation of an equilibrium high-pressure plasma with a fully excluded or reversed magnetic field cannot be investigated using the simplified magnetohydrodynamic and gyrokinetic approaches. A correct description of the electron dynamics under close proximity of zero and strong magnetic field regions can only be achieved within the framework of the kinetic theory that can be implemented most efficiently by the particle-in-cell method. The full-scale particle-in-cell simulations of modern fusion experiments require the use of finite-difference schemes in which temporal steps exceed the period of fastest electron oscillations at the plasma or even cyclotron frequencies and spatial steps do not have to resolve the Debye radius. This requirement is satisfied by implicit schemes capable of conserving the total energy of the system. The particle-in-cell model presented in this paper is based on the energy conserving semi-implicit approach as a predictive step and a new method for suppressing the electrostatic noise inherent in it as a corrective step. This two-step algorithm provides not only accurate energy conservation, but also the exact local fulfillment of the Gauss law.

KW - High-β plasma

KW - Mirror traps

KW - Particle-in-cell simulations

UR - https://www.scopus.com/record/display.uri?eid=2-s2.0-85177809827&origin=inward&txGid=e983685e296932ed34d2c6a0d35487f7

UR - https://www.mendeley.com/catalogue/cbd67319-f4b7-3ab6-8420-785bd3fa2f50/

U2 - 10.1016/j.cpc.2023.109020

DO - 10.1016/j.cpc.2023.109020

M3 - Article

VL - 295

JO - Computer Physics Communications

JF - Computer Physics Communications

SN - 0010-4655

M1 - 109020

ER -