Numerical Thermalization in One- and Two-Dimensional Particle-In-Cell Simulations with Monte-Carlo Collisions

P. Y. Lai,Y. R. Lin-Liu,L. Chen,S. H. Chen
DOI: https://doi.org/10.1109/plasma.2016.7534193
2016-01-01
Abstract:Summary form only given. Numerical thermalization induced by discrete-particle effects was observed in particle-in-cell (PIC) simulations 1 , which show that the thermal relaxation time scales with N D 2 and N D for one-dimensional (1D) and two-dimensional (2D) models 2 , respectively. The parameter N D denotes the number of particles in a Debye length (1D) or a Debye square (2D). However, it was found recently that the thermal relaxation time is anomalously shortened to scale with N D while adding the Monte-Carlo collisions in a 1D PIC simulation. 3,4 In this work, we examine the numerical thermalization in two-dimensional ES PIC simulations with the consideration of electron-ion collision and electron-neutral collision. 4 Our results show that the thermal relaxation time is less sensitive to the collision frequency as compared to the 1D cases, and the thermal relaxation time remains to scale with N D as predicted by the theoretical analysis using the Balescu-Lenard-Landau kinetic theory 5 .
What problem does this paper attempt to address?