Discontinuous Galerkin Methods for Relativistic Vlasov–Maxwell System

He Yang, Fengyan Li

Research output: Contribution to journalArticle

2 Citations (Scopus)

Abstract

The relativistic Vlasov–Maxwell (RVM) system is a kinetic model that describes the dynamics of plasma when the charged particles move in the relativistic regime and their collisions are not important. In this paper, we formulate and investigate discontinuous Galerkin (DG) methods to solve the RVM system. When standard piecewise polynomial functions are used to define trial and test spaces, the methods conserve mass as expected. However the energy conservation does not hold due to the specific form of the total energy of the system. In order to obtain provable mass and energy conservation, we take advantage of the flexibility of DG discretizations and enrich the discrete spaces with some non-polynomial function. For the semi-discrete DG methods with standard and enriched spaces, stability and error estimates are established together with their properties in conservation. In actual implementation with the enriched space, special care is needed to reduce the loss of significance for better numerical stability. Numerical experiments, including streaming Weibel instability and wakefield acceleration, are presented to demonstrate the performance of the methods. Positivity-preserving limiter is also used in simulating wakefield acceleration to obtain physically more relevant solutions.

Original languageEnglish (US)
Pages (from-to)1216-1248
Number of pages33
JournalJournal of Scientific Computing
Volume73
Issue number2-3
DOIs
StatePublished - Dec 1 2017

Fingerprint

Discontinuous Galerkin Method
Convergence of numerical methods
Galerkin methods
Energy conservation
Limiters
Energy Conservation
Charged particles
Conservation
Polynomials
Plasmas
Kinetics
Limiter
Mass Conservation
Discontinuous Galerkin
Stability Estimates
Conserve
Piecewise Polynomials
Numerical Stability
Kinetic Model
Polynomial function

Keywords

  • Conservation
  • Discontinuous Galerkin methods
  • Non-polynomial space
  • Relativistic Vlasov–Maxwell system
  • Streaming Weibel instability
  • Wakefield acceleration

ASJC Scopus subject areas

  • Software
  • Theoretical Computer Science
  • Numerical Analysis
  • Engineering(all)
  • Computational Theory and Mathematics
  • Computational Mathematics
  • Applied Mathematics

Cite this

Discontinuous Galerkin Methods for Relativistic Vlasov–Maxwell System. / Yang, He; Li, Fengyan.

In: Journal of Scientific Computing, Vol. 73, No. 2-3, 01.12.2017, p. 1216-1248.

Research output: Contribution to journalArticle

@article{0a53b7e400584fcaa1db2d63c318ebdf,
title = "Discontinuous Galerkin Methods for Relativistic Vlasov–Maxwell System",
abstract = "The relativistic Vlasov–Maxwell (RVM) system is a kinetic model that describes the dynamics of plasma when the charged particles move in the relativistic regime and their collisions are not important. In this paper, we formulate and investigate discontinuous Galerkin (DG) methods to solve the RVM system. When standard piecewise polynomial functions are used to define trial and test spaces, the methods conserve mass as expected. However the energy conservation does not hold due to the specific form of the total energy of the system. In order to obtain provable mass and energy conservation, we take advantage of the flexibility of DG discretizations and enrich the discrete spaces with some non-polynomial function. For the semi-discrete DG methods with standard and enriched spaces, stability and error estimates are established together with their properties in conservation. In actual implementation with the enriched space, special care is needed to reduce the loss of significance for better numerical stability. Numerical experiments, including streaming Weibel instability and wakefield acceleration, are presented to demonstrate the performance of the methods. Positivity-preserving limiter is also used in simulating wakefield acceleration to obtain physically more relevant solutions.",
keywords = "Conservation, Discontinuous Galerkin methods, Non-polynomial space, Relativistic Vlasov–Maxwell system, Streaming Weibel instability, Wakefield acceleration",
author = "He Yang and Fengyan Li",
year = "2017",
month = "12",
day = "1",
doi = "10.1007/s10915-016-0332-4",
language = "English (US)",
volume = "73",
pages = "1216--1248",
journal = "Journal of Scientific Computing",
issn = "0885-7474",
publisher = "Springer New York",
number = "2-3",

}

TY - JOUR

T1 - Discontinuous Galerkin Methods for Relativistic Vlasov–Maxwell System

AU - Yang, He

AU - Li, Fengyan

PY - 2017/12/1

Y1 - 2017/12/1

N2 - The relativistic Vlasov–Maxwell (RVM) system is a kinetic model that describes the dynamics of plasma when the charged particles move in the relativistic regime and their collisions are not important. In this paper, we formulate and investigate discontinuous Galerkin (DG) methods to solve the RVM system. When standard piecewise polynomial functions are used to define trial and test spaces, the methods conserve mass as expected. However the energy conservation does not hold due to the specific form of the total energy of the system. In order to obtain provable mass and energy conservation, we take advantage of the flexibility of DG discretizations and enrich the discrete spaces with some non-polynomial function. For the semi-discrete DG methods with standard and enriched spaces, stability and error estimates are established together with their properties in conservation. In actual implementation with the enriched space, special care is needed to reduce the loss of significance for better numerical stability. Numerical experiments, including streaming Weibel instability and wakefield acceleration, are presented to demonstrate the performance of the methods. Positivity-preserving limiter is also used in simulating wakefield acceleration to obtain physically more relevant solutions.

AB - The relativistic Vlasov–Maxwell (RVM) system is a kinetic model that describes the dynamics of plasma when the charged particles move in the relativistic regime and their collisions are not important. In this paper, we formulate and investigate discontinuous Galerkin (DG) methods to solve the RVM system. When standard piecewise polynomial functions are used to define trial and test spaces, the methods conserve mass as expected. However the energy conservation does not hold due to the specific form of the total energy of the system. In order to obtain provable mass and energy conservation, we take advantage of the flexibility of DG discretizations and enrich the discrete spaces with some non-polynomial function. For the semi-discrete DG methods with standard and enriched spaces, stability and error estimates are established together with their properties in conservation. In actual implementation with the enriched space, special care is needed to reduce the loss of significance for better numerical stability. Numerical experiments, including streaming Weibel instability and wakefield acceleration, are presented to demonstrate the performance of the methods. Positivity-preserving limiter is also used in simulating wakefield acceleration to obtain physically more relevant solutions.

KW - Conservation

KW - Discontinuous Galerkin methods

KW - Non-polynomial space

KW - Relativistic Vlasov–Maxwell system

KW - Streaming Weibel instability

KW - Wakefield acceleration

UR - http://www.scopus.com/inward/record.url?scp=85006381096&partnerID=8YFLogxK

UR - http://www.scopus.com/inward/citedby.url?scp=85006381096&partnerID=8YFLogxK

U2 - 10.1007/s10915-016-0332-4

DO - 10.1007/s10915-016-0332-4

M3 - Article

AN - SCOPUS:85006381096

VL - 73

SP - 1216

EP - 1248

JO - Journal of Scientific Computing

JF - Journal of Scientific Computing

SN - 0885-7474

IS - 2-3

ER -