Comprehensive study of Kinetic Trajectory Simulation method for multi-component magnetized plasma-wall interaction process

Abstract

For a wide range of plasma applications across diverse fields, a comprehensive understanding of the plasma-wall interaction mechanism is indispensable due to its inherent connection with confined plasma. This work compilation delves into the Kinetic Trajectory Simulation (KTS) method for the interaction of multi-component magnetized plasma with wall, specifically focusing on its implications for the tungsten wall sputtering model. In the evolution of the 1d3v (one dimensional spatial coordinates and three dimensional velocity coordinates) KTS method, the coupled set of kinetic equations has been solved under specified boundary conditions which yields results of higher accuracy. At the particle injection boundary, we have assumed the velocity distribution function of particle species to be cut-off Maxwellians, meeting essential requirements for plasma-wall transition processes: quasineutrality, sheath edge singularity, continuity of macroscopic fluid variables, and the kinetic Bohm sheath condition. The kinetic Bohm sheath condition, a fundamental criterion for plasma sheath formation, is extended for multi-component plasmas, accounting for the cut-off Maxwellian distribution of negatively charged particles. A comparative study of the kinetic Bohm sheath condition for cut-off and Boltzmann distributions reveals a deviation of less than 2.0% in magnitude. The concentration ratio of positive or negative ion species and the presheath side electron temperature influence various plasma-wall transition characteristics, including wall potential, Debye sheath thickness, particle densities, potential distribution, particle fluxes towards the surface, particle drift velocity, phase-space trajectory evolution, and physical sputtering of the tungsten surface. Although lighter ions possess higher energy when striking the surface, the physical sputtering yield of the tungsten surface is greater for heavier ions due to their lower threshold energy and larger collision cross-section. Furthermore, a comparative study of plasma-wall transition properties using kinetic and fluid approaches demonstrates qualitative similarities, with a notable deviation of approximately 4.0% in the magnitude in the vicinity of the material surface.