Error estimates of Runge–Kutta discontinuous galerkin methods for the Vlasov–Maxwell system

Yang, He¹ ; Li, Fengyan¹

¹ Department of Mathematical Sciences, Rensselaer Polytechnic Institute, Troy, 12180-3590 New York, USA.

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 49 (2015) no. 1, pp. 69-99

Résumé

In this paper, error analysis is established for Runge–Kutta discontinuous Galerkin (RKDG) methods to solve the Vlasov–Maxwell system. This nonlinear hyperbolic system describes the time evolution of collisionless plasma particles of a single species under the self-consistent electromagnetic field, and it models many phenomena in both laboratory and astrophysical plasmas. The methods involve a third order TVD Runge–Kutta discretization in time and upwind discontinuous Galerkin discretizations of arbitrary order in phase domain. With the assumption that the exact solutions have sufficient regularity, the $L^2$ errors of the particle number density function as well as electric and magnetic fields at any given time $T$ are bounded by $C{h}^{k+1/2}+C\tau {}^3$ under a CFL condition $\tau /h\le \gamma$. Here $k$ is the polynomial degree used in phase space discretization, satisfying $k>\frac{d_x+1}{2}$ (with $d_x$ being the dimension of spatial domain), $\tau$ is the time step, and $h$ is the maximum mesh size in phase space. Both $C$ and $\gamma$ are positive constants independent of $h$ and $\tau$, and they may depend on the polynomial degree $k$, time $T$, the size of the phase domain, certain mesh parameters, and some Sobolev norms of the exact solution. The analysis can be extended to RKDG methods with other numerical fluxes and to RKDG methods solving relativistic Vlasov–Maxwell equations.

Détail

Reçu le : 2013-12-03
DOI : 10.1051/m2an/2014025

Classification : 65M15, 65M60, 65M06, 35Q83, 35L50
Keywords: Vlasov–Maxwell system, Runge–Kutta discontinuous Galerkin methods, error estimates