10-14 July 2016
Africa/Johannesburg timezone
Paper Submission is Open!! Registration for CCP2016 Conference will be open at the registration desk @ St. George's Hotel reception area on Sunday 10 July / Welcome Dinner at 18:30 SAST.

An asymptotic preserving scheme for the Vlasov-Maxwell system in the classical limit

13 Jul 2016, 11:35


Dr Lukas Einkemmer (University of Innsbruck)

Abstract content <br> &nbsp; (Max 300 words)<br><a href="http://events.saip.org.za/getFile.py/access?resId=0&materialId=0&confId=34" target="_blank">Formatting &<br>Special chars</a>

In order to model non-thermalized plasmas in the fully magnetized and relativistic regime, the relativistic Vlasov-Maxwell (RVM) system is most commonly used. In that context efficient numerical methods have been developed (see, e.g., [1] and [2]). The dimensionless parameter c that determines the relative strength of relativistic effects is the speed of light scaled by the thermal velocity. In some applications, however, the weakly relativistic regime (i.e. where c is large) is of interest. In this case numerical simulations show that the methods that have been developed for the RVM system are extremely inefficient. In addition, many numerical methods are not even able to recover the correct limit behavior (i.e. convergence to the Vlasov-Poisson system).

In this presentation we describe the asymptotic preserving scheme introduced in [3]. This method is robust with respect to the classic limit and imposes no step size restriction. Our approach relies on a time splitting approach for the RVM system and employs an implicit integrator for the linear part of Maxwell's equations. It turns out that the choice of this implicit method is crucial as even L-stable methods can introduce numerical instabilities for large values of c.

To illustrate the efficiency of the numerical scheme in both the relativistic as well as in the classical limit regime we present the results for a number of numerical simulations. We mostly focus on the Weibel instability and on Landau damping. In the relativistic case we compare the asymptotic preserving scheme to the Hamiltonian splitting proposed in [1] and the VALIS method proposed in [2].

[1] N. Crouseilles, L. Einkemmer, E. Faou, J. Comput. Phys. 283, pp.224-240 (2015)
[2] N. J. Sircombe, T.D. Arber, J. Comput. Phys. 228, pp. 4773-4788, (2009)
[3] N. Crouseilles, L. Einkemmer, E. Faou, arXiv 1602.09062, (2016)

Primary author

Dr Lukas Einkemmer (University of Innsbruck)


Prof. Erwan Faou (INRIA-Rennes, Universit\'e de Rennes 1 and IRMAR, Rennes, France) Dr Nicolas Crouseilles (INRIA-Rennes, Universit\'e de Rennes 1 and IRMAR, Rennes, France)

Presentation Materials

There are no materials yet.