Kinetic simulations of collision-less plasmas in open magnetic geometries

Abstract

Laboratory plasmas in open magnetic geometries can be found in many different applications such as (1) Scrape-Of-Layer (SOL) and divertor regions in toroidal confinement fusion devices (≈ 1 − 10 2 eV), (2) linear divertor simulators (≈ 1 − 10 eV), (3) plasma-based thrusters (≈ 10 eV) and (4) magnetic mirrors (≈ 10 2 − 10 3 eV). A common feature of these plasma systems is the need to resolve, in addition to velocity space, at least one physical dimension (e.g. along flux lines) to capture the relevant physics. In general, this requires a kinetic treatment. Fully kinetic Particle-In-Cell (PIC) simulations can be applied but at the expense of large computational effort. A common way to resolve this is to use a hybrid approach: kinetic ions and fluid electrons. In the present work, the development of a hybrid PIC computational tool suitable for open magnetic geometries is described which includes (1) the effect of non-uniform magnetic fields, (2) finite fully-absorbing boundaries for the particles and (3) volumetric particle sources. Analytical expressions for the momentum transport in the paraxial limit are presented with their underlying assumptions and are used to validate the results from the PIC simulations. A general method is described to construct discrete particle distribution functions in state of mirror-equilibrium. This method is used to obtain the initial state for the PIC simulation. Collisionless simulations in a mirror geometry are performed. Results show that the effect of magnetic compression is correctly described and momentum is conserved. The self-consistent electric field is calculated and is shown to modify the ion velocity distribution function in manner consistent with analytic theory. Based on this analysis, the ion distribution function is understood in terms of a loss-cone distribution and an isotropic Maxwell-Boltzmann distribution driven by a volumetric plasma source. Finally, inclusion of a Monte-Carlo based Fokker-Planck collision operator is discussed in the context of future work. † This manuscript has been authored by UT-Battelle, LLC, under contract DE-AC05-00OR22725 with the US Department of Energy (DOE). The US government retains and the publisher, by accepting the article for publication, acknowledges that the US government retains a nonexclusive, paid-up, irrevocable, worldwide license to publish or reproduce the published form of this manuscript, or allow others to do so, for US government purposes. DOE will provide public access to these results of federally sponsored research in accordance with the DOE Public Access Plan (http://energy.gov/downloads/doepublic-access-plan).