https://doi.org/10.1140/epjd/e2014-50180-9
Regular Article
A new fully two-dimensional conservative semi-Lagrangian method: applications on polar grids, from diocotron instability to ITG turbulence*
1 Inria-Rennes (IPSO team) &
IRMAR, Rennes,
France
2 Inria-Nancy (TONUS team) &
IRMA, Strasbourg,
France
3 Institut Jean Lamour,
Nancy,
France
4 Max-Planck-Institut für Plasmaphysik,
Garching,
Germany
5 Observatoire Astronomique de
Strasbourg, Strasbourg, France
a
e-mail: mehrenbe@math.unistra.fr
Received:
5
March
2014
Received in final form:
26
May
2014
Published online:
18
September
2014
While developing a new semi-Lagrangian solver, the gap between a linear Landau run in 1D × 1D and a 5D gyrokinetic simulation in toroidal geometry is quite huge. Intermediate test cases are welcome for testing the code. A new fully two-dimensional conservative semi-Lagrangian (CSL) method is presented here and is validated on 2D polar geometries. We consider here as building block, a 2D guiding-center type equation on an annulus and apply it on two test cases. First, we revisit a 2D test case previously done with a PIC approach [J. Pétri, A&A 503, 1 (2009)] and detail the boundary conditions. Second, we consider a 4D drift-kinetic slab simulation (see [V. Grandgirard, M. Brunetti, P. Bertrand, N. Besse, X. Garbet, P. Ghendrih, G. Manfredi, Y. Sarazin, O. Sauter, E. Sonnendrücker, J. Vaclavik, L. Villard, J. Comput. Phys. 217, 395 (2006)]). In both cases, the new method appears to be a good alternative to deal with this type of models since it improves the lack of mass conservation of the standard semi-Lagrangian (BSL) method.
© EDP Sciences, Società Italiana di Fisica, Springer-Verlag 2014