Processing math: 100%

2506.16676

Total: 1

#1 Fast solvers for Tokamak fluid models with PETSC -- Part I [PDF] [Copy] [Kimi] [REL]

Authors: Mark F. Adams, Jin Chen, Benjamin Sturdevant

This report develops the first step in adding multigrid solvers to scientific and engineering-relevant magnetohydrodynamics (MHD) models of Tokamaks. These models are characterized by a distinguished direction in the toroidal coordinate that is partially aligned with the magnetic guide field, which dominates the plasma dynamics. All Tokamak models exploit this structure, for example, NIMROD (https://nimrodteam.org/) uses 2D, unstructured, high-order finite elements in the poloidal plane with Fourier modes in the toroidal coordinate, and the 3D, extended MHD code M3D-C1 (https://w3.pppl.gov/~nferraro/m3dc1.html) uses 2D, unstructured C1 elements in the poloidal plane with cubic Hermite functions in the toroidal direction. This structure suggests adding toroidal semi-coarsening multigrid to the existing solver and thereby reducing reliance on direct solvers, which do not scale optimally and are not well suited to modern hardware that demands extreme levels of parallelism. This report focuses on the velocity solve in M3D-C1, using the PETSC -- the Portable, Extensible Toolkit for Scientific Computation -- numerical library (https://petsc.org), and shows that with little new application code, one-dimensional multigrid is about 5x faster than the existing one-level method on an MHD disruption, with runaway electrons, test problem.

Subjects: Plasma Physics , Performance

Publish: 2025-06-20 01:36:08 UTC