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 (15)

• Relativistic Reconnection: an Efficient Source of Non-Thermal Particles
• SMILEI: a collaborative, open-source, multi-purpose particle-in-cell code for plasma simulation
• Formation of Hard Power-laws in the Energetic Particle Spectra Resulting from Relativistic Magnetic Reconnection
• First-Principles Plasma Simulations of Black-Hole Jet Launching
• The serendipity family of finite elements
• Vlasov methods in space physics and astrophysics
• Discontinuous Galerkin algorithms for fully kinetic plasmas
• Multi-dimensional simulations of ergospheric pair discharges around black holes
• Electron Weibel Instability in Relativistic Counter-Streaming Plasmas with Flow-Aligned External Magnetic Fields
• Collisionless accretion onto black holes: dynamics and flares
• On the Manifestly Covariant Juttner Distribution and Equipartition Theorem
• Conservative Discontinuous Galerkin Schemes for Nonlinear Fokker-Planck Collision Operators
• Loading Relativistic Maxwell Distributions in Particle Simulations
• Noise-Induced Magnetic Field Saturation in Kinetic Simulations
• A positivity-preserving and conservative high-order flux reconstruction method for the polyatomic Boltzmann--BGK equation