An optimization principle for the computation of MHD equilibria in the solar corona

Abstract

AIMS: We develop an optimization principle for computing stationary MHD equilibria. METHODS: Our code for the self-consistent computation of the coronal magnetic fields and the coronal plasma uses non-force-free MHD equilibria. Previous versions of the code have been used to compute non-linear force-free coronal magnetic fields from photospheric measurements. The program uses photospheric vector magnetograms and coronal EUV images as input. We tested our reconstruction code with the help of a semi-analytic MHD-equilibrium. The quality of the reconstruction was judged by comparing the exact and reconstructed solution qualitatively by magnetic field-line plots and EUV-images and quantitatively by several different numerical criteria. RESULTS: Our code is able to reconstruct the semi-analytic test equilibrium with high accuracy. The stationary MHD optimization code developed here has about the same accuracy as its predecessor, a non-linear force-free optimization code. The computing time for MHD-equilibria is, however, longer than for force-free magnetic fields. We also extended a well-known class of nonlinear force-free equilibria to the non-force-free regime for purposes of testing the code. CONCLUSIONS: We demonstrate that the code works in principle using tests with analytical equilibria, but it still needs to be applied to real data.

6 pages, 3 figures

References in corpus (7)

• 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
• Comparing magnetic field extrapolations with measurements of magnetic loops
• Nonlinear force-free magnetic field extrapolations: comparison of the Grad-Rubin and Wheatland-Sturrock-Roumeliotis algorithm
• Magnetic modelling and tomography: First steps towards a consistent reconstruction of the solar corona
• Links between magnetic fields and plasma flows in a coronal hole
• How to use magnetic field information for coronal loop identification?

Cited by in corpus (3)

• Can We Improve the Preprocessing of Photospheric Vector Magnetograms by the Inclusion of Chromospheric Observations?
• Optimization approach for the computation of magnetohydrostatic coronal equilibria in spherical geometry
• Coronal Alfven speeds in an isothermal atmosphere I. Global properties