Mimetic methods for Lagrangian relaxation of magnetic fields

Candelaresi, S. , Pontin, D. and Hornig, G. (2014) Mimetic methods for Lagrangian relaxation of magnetic fields. SIAM Journal on Scientific Computing, 36(6), B952-B968. (doi: 10.1137/140967404)

Full text not currently available from Enlighten.


We present a new code that performs a relaxation of a magnetic field toward 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, which leads to an enhanced computing speed compared to previous relaxation codes.

Item Type:Articles
Glasgow Author(s) Enlighten ID:Candelaresi, Dr Simon
Authors: Candelaresi, S., Pontin, D., and Hornig, G.
College/School:College of Science and Engineering > School of Mathematics and Statistics > Mathematics
Journal Name:SIAM Journal on Scientific Computing
Publisher:Society for Industrial and Applied Mathematics
ISSN (Online):1095-7197
Published Online:10 December 2014

University Staff: Request a correction | Enlighten Editors: Update this record