Controlling the weights of simulation particles: adaptive particle management using kd trees
Abstract
In particle simulations, the weights of particles determine how many physical particles they represent. Adaptively adjusting these weights can greatly improve the efficiency of the simulation, without creating severe nonphysical artifacts. We present a new method for the pairwise merging of particles, in which two particles are combined into one. To find particles that are ‘close' to each other, we use a kd tree data structure. With a kd tree, close neighbors can be searched for efficiently, and independently of the mesh used in the simulation. The merging can be done in different ways, conserving for example momentum or energy. We introduce probabilistic schemes, which set properties for the merged particle using random numbers. The effect of various merge schemes on the energy distribution, the momentum distribution and the grid moments is compared. We also compare their performance in the simulation of the twostream instability.
 Publication:

Journal of Computational Physics
 Pub Date:
 February 2014
 DOI:
 10.1016/j.jcp.2013.12.005
 arXiv:
 arXiv:1301.1552
 Bibcode:
 2014JCoPh.259..318T
 Keywords:

 Physics  Computational Physics
 EPrint:
 Journal of Computational Physics 259 (2014) 318330