Solver comparison for Poisson-like equations on tokamak geometries
Bourne E, Leleux P, Kormann K, Kruse C, Grandgirard V, Güçlü Y, Kühn MJ, Rüde U, Sonnendrücker E, Zoni E (2023)
Publication Type: Journal article
Publication year: 2023
Journal
Book Volume: 488
Pages Range: 112249
Article Number: 112249
DOI: 10.1016/j.jcp.2023.112249
Abstract
The solution of Poisson-like equations defined on a complex geometry is required for gyrokinetic simulations, which are important for the modelling of plasma turbulence in nuclear fusion devices such as the ITER tokamak. In this paper, we compare three existing solvers finely tuned to solve this problem, in terms of the accuracy of the solution, and their computational efficiency. We also consider practical implementation aspects, including the parallel efficiency of the code, potentially enabling an integration of the solvers in a state-of-the-art first-principle gyrokinetic simulation framework. The first, the Spline FEM solver, uses polar splines to construct a finite elements method which solves the equation on curvilinear coordinates. The resulting linear system is solved using a conjugate gradient method. The second, the GMGPolar solver, uses a symmetric finite difference method to discretise the differential equation. The resulting linear system is solved using a tailored geometric multigrid scheme, with a combination of zebra circle and radial line smoothers, together with an implicit extrapolation scheme. The third, the Embedded Boundary solver, uses a finite volumes method on Cartesian coordinates with an embedded boundary scheme. The resulting linear system is solved using a multigrid scheme. The Spline FEM solver is shown to be the most accurate. The GMGPolar solver is shown to use the least memory. The Embedded Boundary solver is shown to be the fastest in most cases. All three solvers are shown to be capable of solving the equation on a realistic non-analytical geometry. The Embedded Boundary solver is additionally used to attempt to solve an X-point geometry.
Authors with CRIS profile
Additional Organisation(s)
Involved external institutions
Atomic Energy and Alternative Energies Commission / Commissariat à l'Énergie Atomique et aux Énergies Alternatives (CEA)
France (FR)
Max-Planck-Institut für Plasmaphysik / Max Planck Institute for Plasma Physics (MPG IPP)
Germany (DE)
Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) / German Aerospace Center
Germany (DE)
Centre Européen de Recherche et de Formation Avancée en Calcul Scientifique (CERFACS)
France (FR)
Ruhr-Universität Bochum (RUB)
Germany (DE)
Lawrence Berkeley National Laboratory (LBNL)
United States (USA) (US)
France (FR)
Max-Planck-Institut für Plasmaphysik / Max Planck Institute for Plasma Physics (MPG IPP)
Germany (DE)
Deutsches Zentrum für Luft- und Raumfahrt e.V. (DLR) / German Aerospace Center
Germany (DE)
Centre Européen de Recherche et de Formation Avancée en Calcul Scientifique (CERFACS)
France (FR)
Ruhr-Universität Bochum (RUB)
Germany (DE)
Lawrence Berkeley National Laboratory (LBNL)
United States (USA) (US)
How to cite
APA:
Bourne, E., Leleux, P., Kormann, K., Kruse, C., Grandgirard, V., Güçlü, Y.,... Zoni, E. (2023). Solver comparison for Poisson-like equations on tokamak geometries. Journal of Computational Physics, 488, 112249. https://doi.org/10.1016/j.jcp.2023.112249
MLA:
Bourne, Emily, et al. "Solver comparison for Poisson-like equations on tokamak geometries." Journal of Computational Physics 488 (2023): 112249.
BibTeX: Download