2014/05/05 by Simon Candelaresi, Candelaresi, Simon, D. I. Pontin +3 · 1 citation
Engineering · Physics and Astronomy · #47F05 #65D25 #65M06 #68W40 #76W05 #85-08 #85A30 #Advanced Numerical Methods in Computational Mathematics #Computational Physics (physics.comp-ph) #Electromagnetic Scattering and Analysis #FOS: Physical sciences #Magnetic confinement fusion research #Plasma Physics (physics.plasm-ph) #Solar and Stellar Astrophysics (astro-ph.SR)
paper · pdf · doi:10.48550/arxiv.1405.0942
openalex publication_date 2014/05/05 · openalex created_date 2025/10/10 · openalex updated_date 2026/07/28
We present a new code that performs a relaxation of a magnetic field towards a force-free state (Beltrami field) using a Lagrangian numerical scheme. Beltrami fields are of interest for the dynamics of many technical and astrophysical plasmas as they are the lowest energy states that the magnetic field can reach. The numerical method strictly preserves the magnetic flux and the topology of magnetic field lines. In contrast to other implementations we use mimetic operators for the spatial derivatives in order to improve accuracy for high distortions of the grid. Compared with schemes using direct derivatives we find that the final state of the simulation approximates a force-free state with a significantly higher accuracy. We implement the scheme in a code which runs on graphical processing units (GPU), which leads to an enhanced computing speed compared to previous relaxation codes.