Numerical heating in particle-in-cell simulations with Monte Carlo binary collisions

TL;DR

This paper identifies and analyzes the numerical heating caused by the coupling of Monte Carlo collision algorithms with particle-in-cell simulations, which can distort plasma evolution over long times.

Contribution

It reveals the origin of MC-induced numerical heating in PIC simulations and proposes strategies to mitigate this nonphysical effect.

Findings

MC coupling causes artificial electromagnetic energy production

Numerical heating significantly impacts long-term plasma simulation accuracy

Analytical description of the heating source aids in developing mitigation strategies

Abstract

The binary Monte Carlo (MC) collision algorithm is a standard and robust method to include binary Coulomb collision effects in particle-in-cell (PIC) simulations of plasmas. Here, we show that the coupling between PIC and MC algorithms can give rise to (nonphysical) numerical heating of the system, that significantly exceeds that observed when these algorithms operate independently. We argue that this deleterious effect results from an inconsistency between the particle motion associated with MC-collisions and the work performed by the collective electromagnetic field on the PIC grid. This inconsistency manifests as the (artificial) stochastic production of electromagnetic energy, which ultimately heats the plasma particles. The MC-induced numerical heating can significantly impact the evolution of the simulated system for long simulation times ($\gtrsim 10^3$ collision periods, for typical numerical parameters). We describe the source of the MC-induced numerical heating analytically and discuss strategies to minimize it.