Bereiche | Tage | Auswahl | Suche | Aktualisierungen | Downloads | Hilfe
P: Fachverband Plasmaphysik
P 15: Helmholtz Graduate School II
P 15.10: Poster
Dienstag, 14. März 2017, 16:30–18:30, HS Foyer
Semi-Lagrangian drift-kinetic simulations: field-aligned interpolation and splitting in complex geometries — •Edoardo Zoni1,2, Yaman Güçlü1, Michel Mehrenberger3, and Eric Sonnendrücker1,2 — 1Max-Planck-Institut für Plasmaphysik, Germany — 2Zentrum Mathematik, TU München, Germany — 3Institut de Recherche Mathématique Avancée, Université de Strasbourg, France
Global turbulence simulations of magnetic fusion devices based on the solution of the gyrokinetic Vlasov-Maxwell equations are computationally very expensive because the thermal ion Larmor radius must be resolved.
The computational burden may be reduced by aligning one grid coordinate with the local magnetic field line, along which the gradients are known to be small. Unfortunately, this methodology poses restrictions on the poloidal mesh, and cannot easily handle complex magnetic field configurations found in diverted Tokamaks and Stellarators.
An alternative and more flexible approach was developed by Ottaviani and Hariri, where local field-aligned differentiation (or interpolation) was performed between adjacent poloidal planes.
Such a method was adapted to the semi-Lagrangian context and combined with dimensional splitting of the transport equation by Latu et al.. We now extend it to general curvilinear coordinates.
We describe here the mathematical formulation and the details of our field-aligned interpolation algorithm. Our code is verified with a linear dispersion analysis for the ITG instability in screw-pinch configuration.