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

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.

References in corpus (4)

• Particle acceleration in relativistic collisionless shocks: Fermi process at last?
• Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection
• One-to-one direct modeling of experiments and astrophysical scenarios: pushing the envelope on kinetic plasma simulations
• Statistical kinetic treatment of relativistic binary collisions

Cited by in corpus (3)

• Data-driven discovery of reduced plasma physics models from fully-kinetic simulations
• Exploring ultra-high-intensity wakefields in carbon nanotube arrays: an effective plasma-density approach
• Implementation of higher-order velocity mapping between marker particles and grid in the particle-in-cell code XGC