Abstract
In this work we present a new Eulerian kinetic 1D3V code devoted to the study of plasma-wall interactions and sheath structure. The code solves the Vlasov-Poisson system for two or more kinetic species on a one-dimensional spatial grid between two limiting plates, in the presence of a uniform magnetic field tilted with respect to the normal to the plates. A source-sink term in the Vlasov equation leads to a stationary solution comprising a magnetic and a collisional presheath, as well as the usual Debye sheath in front of the wall. Thanks to a nonuniform spatial grid the sheath structure can be resolved accurately. Several advection schemes are implemented and the code is parallelized. Here, the code is used to illustrate the problems of a stationary sheath-presheath structure in front of a negatively biased wall, with and without a tilted magnetic field.
Highlights
Despite the constant improvement of power and accessibility of computing resources, the accurate numerical simulation of plasma-wall interactions remains a challenging task
Numerical results We tested our code on two well-documented steady-state problems: the plasma-wall transition with and without a tilted magnetic field
In order to check that the chosen grid and interpolation scheme do not introduce any artificial widening of the distribution function, either by cell averaging or numerical diffusion, the simulation results are compared against an analytical ballistic model in the sheath (Fig. 5)
Summary
Despite the constant improvement of power and accessibility of computing resources, the accurate numerical simulation of plasma-wall interactions remains a challenging task. 3. Numerical results We tested our code on two well-documented steady-state problems: the plasma-wall transition with and without a tilted magnetic field. Numerical results We tested our code on two well-documented steady-state problems: the plasma-wall transition with and without a tilted magnetic field In both cases, we considered a half-infinite plasma between an absorbing plate (located at x = 0) and a bulk plasma at equilibrium (x = L). In order to check that the chosen grid and interpolation scheme do not introduce any artificial widening of the distribution function, either by cell averaging or numerical diffusion, the simulation results are compared against an analytical ballistic model in the sheath (Fig. 5). A slight erosion of the IVDF maximum can be seen (mostly due to the action of the limiters), the IVDF geometry, and most importantly its asymmetry, are well preserved by the PFC3 scheme, and spurious oscillations are controlled
Sign-in/Register to view highlights & summary Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Schedule a callDisclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.
