English

On hp-Streamline Diffusion and Nitsche schemes for the Relativistic Vlasov-Maxwell System

Numerical Analysis 2017-11-02 v1

Abstract

We study stability and convergence of hphp-streamline diffusion (SD) finite element, and Nitsche's schemes for the three dimensional, relativistic (3 spatial dimension and 3 velocities), time dependent Vlasov-Maxwell system and Maxwell's equations, respectively. For the hphp scheme for the Vlasov-Maxwell system, assuming that the exact solution is in the Sobolev space Hs+1(Ω)H^{s+1}(\Omega), we derive global {\sl a priori} error bound of order O(h/p)s+1/2{\mathcal O}(h/p)^{s+1/2}, where h(=maxKhK)h (= \max_K h_K) is the mesh parameter and p(=maxKpK)p (= \max_K p_K) is the spectral order. This estimate is based on the local version with hK=\mboxdiamKh_K=\mbox{ diam } K being the diameter of the {\sl phase-space-time} element KK and pKp_K is the spectral order (the degree of approximating finite element polynomial) for KK. As for the Nitsche's scheme, by a simple calculus of the field equations, first we convert the Maxwell's system to an {\sl elliptic type} equation. Then, combining the Nitsche's method for the spatial discretization with a second order time scheme, we obtain optimal convergence of O(h2+k2){\mathcal O}(h^2+k^2), where hh is the spatial mesh size and kk is the time step. Here, as in the classical literature, the second order time scheme requires higher order regularity assumptions. Numerical justification of the results, in lower dimensions, is presented and is also the subject of a forthcoming computational work [20].

Keywords

Cite

@article{arxiv.1711.00271,
  title  = {On hp-Streamline Diffusion and Nitsche schemes for the Relativistic Vlasov-Maxwell System},
  author = {Mohammad Asadzadeh and Piotr Kowalczyk and Christoffer Standar},
  journal= {arXiv preprint arXiv:1711.00271},
  year   = {2017}
}

Comments

24 pages, 4 figures