Nonlinear force-free magnetic field extrapolations: comparison of the Grad-Rubin and Wheatland-Sturrock-Roumeliotis algorithm

Abstract

We compare the performance of two alternative algorithms which aim to construct a force-free magnetic field given suitable boundary conditions. For this comparison, we have implemented both algorithms on the same finite element grid which uses Whitney forms to describe the fields within the grid cells. The additional use of conjugate gradient and multigrid iterations result in quite effective codes. The Grad-Rubin and Wheatland-Sturrock-Roumeliotis algorithms both perform well for the reconstruction of a known analytic force-free field. For more arbitrary boundary conditions the Wheatland-Sturrock-Roumeliotis approach has some difficulties because it requires overdetermined boundary information which may include inconsistencies. The Grad-Rubin code on the other hand loses convergence for strong current densities. For the example we have investigated, however, the maximum possible current density seems to be not far from the limit beyond which a force free field cannot exist anymore for a given normal magnetic field intensity on the boundary.

21 pages, 13 figures

References in corpus (3)

• Optimization code with weighting function for the reconstruction of coronal magnetic fields
• Preprocessing of vector magnetograph data for a non-linear force-free magnetic field reconstruction
• Magnetic modelling and tomography: First steps towards a consistent reconstruction of the solar corona

Cited by in corpus (6)

• Computing nonlinear force-free fields in spherical geometry
• An optimization principle for the computation of MHD equilibria in the solar corona
• Optimization approach for the computation of magnetohydrostatic coronal equilibria in spherical geometry
• Segmentation of Loops from Coronal EUV Images
• Magnetic Stereoscopy
• Testing non-linear force-free coronal magnetic field extrapolations with the Titov-Demoulin equilibrium