Grid-free tree-code simulations of the plasma-material interaction region

Joint Varenna-Lausanne International Workshop 2014

Salmagne, C.^1,2 ; Reiter, D.^2 ; Gibbon, P.^3

German Research School for Simulation Sciences GmbH, Jülich, Germany^1

Institute of Energy and Climate Research, IEK-4, Forschungszentrum Jülich, Jülich, Germany^2

Institute for Advanced Simulation, Jülich Supertcomputing Centre, Forschungszentrum Jülich, Jülich, Germany^3

关键词: Astrophysical applications;    Electrostatic sheath;    Inclined magnetic fields;    Magnetic field angle;    Massively parallels;    Particle simulations;    Plasma-material interactions;    Plasma-wall transitions;   

Others  :  https://iopscience.iop.org/article/10.1088/1742-6596/561/1/012014/pdf
DOI  :  10.1088/1742-6596/561/1/012014

来源: IOP

【 摘 要 】

A fully kinetic grid-free model based on a Barnes-Hut tree code is used to selfconsistently simulate a collisionless plasma bounded by two floating walls. The workhorse for simulating such plasma wall transition layers is currently the PIC method. However, the present grid-free formulation provides a powerful independent tool to test it [1] and to possibly extend particle simulations towards collisional regimes in a more internally consistent way. Here, we use the grid-free massively parallel Barnes-Hut tree-code PEPC - a well established tool for simulations of Laser-plasmas and astrophysical applications - to develop a 3D ab initio plasma target interaction model. With our approach an electrostatic sheath naturally builds up within the first couple of Debye lengths close to the wall rather than being imposed as a prescribed boundary condition. We verified the code using analytic results [2] as well as 1D PIC simulations [3]. The model was then used to investigate the influence of inclined magnetic fields on the plasma material interface. We used the code to study the correlation between the magnetic field angle and the angular distribution of incident particles.