期刊
MULTISCALE MODELING & SIMULATION
卷 7, 期 2, 页码 865-887出版社
SIAM PUBLICATIONS
DOI: 10.1137/070704939
关键词
Coulomb collisions; plasma; simulation; hybrid method; bump on tail; thermalization
资金
- U.S. DOE by the University of California, Los Angeles [DE-FG02-05ER25710]
- Lawrence Livermore National Laboratory [DE-AC52-07NA27344]
- Office of Advanced Scientific Computing Research, DOE Office of Science
If the collisional time scale for Coulomb collisions is comparable to the characteristic time scales for a plasma, then simulation of Coulomb collisions may be important for computation of kinetic plasma dynamics. This can be a computational bottleneck because of the large number of simulated particles and collisions (or phase-space resolution requirements in continuum algorithms), as well as the wide range of collision rates over the velocity distribution function. This paper considers Monte Carlo simulation of Coulomb collisions using the binary collision models of Takizuka and Abe and of Nanbu. It presents a hybrid method for accelerating the computation of Coulomb collisions. The hybrid method represents the velocity distribution function as a combination of a thermal component (a Maxwellian distribution) and a kinetic component (a set of discrete particles). Collisions between particles from the thermal component preserve the Maxwellian; collisions between particles from the kinetic component are performed using the method of Takizuka and Abe or of Nanbu. Collisions between the kinetic and thermal components are performed by sampling a particle from the thermal component and selecting a particle from the kinetic component. Particles are also transferred between the two components according to thermalization and dethermalization probabilities, which are functions of phase space.
作者
我是这篇论文的作者
点击您的名字以认领此论文并将其添加到您的个人资料中。
推荐
暂无数据