Collisions are an important dissipation mechanism in plasmas. In one-dimensional modelling, a commonly used collision operator is the Lenard-Bernstein operator, or its modified energy- and momentum-conserving counterpart. When approximating such operators numerically, it is important to respect their structure in order to satisfy the laws of thermodynamics. It is, however, challenging to discretise such operators in a structure-preserving way when considering particle methods. In this work, we present a macro-particle discretisation of the Lenard-Bernstein collision operator that is energy and momentum preserving.
Comment: 17 pages, 9 figures