skip to main content


This content will become publicly available on October 10, 2018

Title: Discontinuous Galerkin algorithms for fully kinetic plasmas

Here, we present a new algorithm for the discretization of the non-relativistic Vlasov–Maxwell system of equations for the study of plasmas in the kinetic regime. Using the discontinuous Galerkin finite element method for the spatial discretization, we obtain a high order accurate solution for the plasma's distribution function. Time stepping for the distribution function is done explicitly with a third order strong-stability preserving Runge–Kutta method. Since the Vlasov equation in the Vlasov–Maxwell system is a high dimensional transport equation, up to six dimensions plus time, we take special care to note various features we have implemented to reduce the cost while maintaining the integrity of the solution, including the use of a reduced high-order basis set. A series of benchmarks, from simple wave and shock calculations, to a five dimensional turbulence simulation, are presented to verify the efficacy of our set of numerical methods, as well as demonstrate the power of the implemented features.
ORCiD logo [1] ;  [2] ;  [3] ; ORCiD logo [2] ;  [1]
  1. Univ. of Maryland, College Park, MD (United States)
  2. Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
  3. Univ. of Maryland, College Park, MD (United States); Princeton Univ., Princeton, NJ (United States)
Publication Date:
Grant/Contract Number:
AGS-1622306; FG02-93ER54197; AC02-09CH11466; FA9550-15-1-0193; ACI-1548562
Accepted Manuscript
Journal Name:
Journal of Computational Physics
Additional Journal Information:
Journal Volume: 353; Journal Issue: C; Journal ID: ISSN 0021-9991
Research Org:
Princeton Plasma Physics Lab. (PPPL), Princeton, NJ (United States)
Sponsoring Org:
Country of Publication:
United States
70 PLASMA PHYSICS AND FUSION TECHNOLOGY; Discontinuous Galerkin; Vlasov–Maxwell
OSTI Identifier: