Discontinuous Galerkin Representation of the Maxwell-Jüttner Distribution

Abstract

Kinetic simulations of relativistic gases and plasmas are critical for understanding diverse astrophysical and terrestrial systems, but the accurate construction of the relativistic Maxwellian, the Maxwell-Jüttner (MJ) distribution, on a discrete simulation grid is challenging. Difficulties arise from the finite velocity bounds of the domain, which may not capture the entire distribution function, as well as errors introduced by projecting the function onto a discrete grid. Here we present a novel scheme for iteratively correcting the moments of the projected distribution applicable to all grid-based discretizations of the relativistic kinetic equation. In addition, we describe how to compute the needed nonlinear quantities, such as Lorentz boost factors, in a discontinuous Galerkin (DG) scheme through a combination of numerical quadrature and weak operations. The resulting method accurately captures the distribution function and ensures that the moments match the desired values to machine precision.

References in corpus (5)

• Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection
• Electron Weibel Instability in Relativistic Counter-Streaming Plasmas with Flow-Aligned External Magnetic Fields
• Collisionless accretion onto black holes: dynamics and flares
• Loading Relativistic Maxwell Distributions in Particle Simulations
• A positivity-preserving and conservative high-order flux reconstruction method for the polyatomic Boltzmann--BGK equation