Lagrangian relaxation schemes for calculating force-free magnetic fields, and their limitations

Abstract

Force-free magnetic fields are important in many astrophysical settings. Determining the properties of such force-free fields -- especially smoothness and stability properties -- is crucial to understanding many key phenomena in astrophysical plasmas, for example energy release processes that heat the plasma and lead to dynamic or explosive events. Here we report on a serious limitation on the computation of force-free fields that has the potential to invalidate the results produced by numerical force-free field solvers even for cases in which they appear to converge (at fixed grid resolution) to an equilibrium magnetic field. In the present work we discuss this problem within the context of a Lagrangian relaxation scheme that conserves magnetic flux and div(B) identically. Error estimates are introduced to assess the quality of the calculated equilibrium. We go on to present an algorithm, based on re-writing the curl operation via Stokes' theorem, for calculating the current which holds great promise for improving dramatically the accuracy of the Lagrangian relaxation procedure.

Submitted to ApJ

References in corpus (4)

• General Relativistic Magnetohydrodynamic Simulations of Jet Formation and Large-Scale Propagation from Black Hole Accretion Systems
• Force-Free Models of Magnetically Linked Star-Disk Systems
• Magnetic Braiding and Parallel Electric Fields
• Magnetic Flux Braiding: Force-Free Equilibria and Current Sheets

Cited by in corpus (14)

• Dynamics of braided coronal loops II: Cascade to multiple small-scale reconnection events
• Dynamics of braided coronal loops - I. Onset of magnetic reconnection
• Magnetic Braiding and Quasi-Separatrix Layers
• Heating of braided coronal loops
• An overview of flux braiding experiments
• The structure of current layers and degree of field line braiding in coronal loops
• Braided magnetic fields: equilibria, relaxation and heating
• Formation of current singularity in a topologically constrained plasma
• On the formation of current sheets in response to the compression or expansion of a potential magnetic field
• Magnetic field relaxation and current sheets in an ideal plasma
• Ideal Relaxation of the Hopf Fibration
• Magnetohydrodynamics dynamical relaxation of coronal magnetic fields. I. Parallel untwisted magnetic fields in 2D
• Constructing current singularity in a 3D line-tied plasma
• Simulations of 3D Magnetic Merging: Resistive Scalings for Null Point and QSL Reconnection