Well posed reconstruction of the solar coronal magnetic field

We present and compare two methods for the reconstruction of the solar coronal magnetic field, assumed to be force-free, from photospheric boundary data. Both methods rely on a well posed mathematical boundary value problem and are of the Grad-Rubin type, i.e., the couple $({\vec B},\alpha)$ is computed iteratively. They do differ from each other on the one hand by the way they address the zero-divergence of ${\vec B}$ issue, and on the other hand by the scheme they use for computing  α at each iteration. The comparison of the two methods is done by numerically computing two examples of nonlinear force-free fields associated to large scale strong electric current distributions, whose exact forms can be otherwise determined semi-analytically. In particular, the second solution has a large nonlinearity even in the weak field region – a feature which is not present in the actual magnetograms, but is interesting to consider as it does allow to push the methods to the limits of their range of validity. The best results obtained with those methods give a relative vector error smaller than 0.01. For the latter extreme case, our results show that higher resolution reconstructions with bounded convergence improve the approximated solution, which may be of some interest for the treatment of the data of future magnetographs.

[1]  T. Amari,et al.  3D Coronal magnetic field from vector magnetograms: non-constant-alpha force-free configuration of the active region NOAA 8151 , 2002 .

[2]  Patrick Ciarlet,et al.  On traces for functional spaces related to Maxwell's equations Part II: Hodge decompositions on the boundary of Lipschitz polyhedra and applications , 2001 .

[3]  G. Roumeliotis The “Stress-and-Relax” Method for Reconstructing the Coronal Magnetic Field from Vector Magnetograph Data , 1996 .

[4]  Parallel Construction of Nonlinear Force-Free Fields , 2004 .

[5]  G. Simon,et al.  Solar active region evolution : comparing models with observations : proceedings of Fourteenth International Summer Workshop National Solar Observatory/Sacramento Peak, Sunspot, New Mexico 88349, USA, 30 August - 3 September 1993 , 1994 .

[6]  B. Low Three-dimensional Structures of Magnetostatic Atmospheres. VII. Magnetic Flux Surfaces and Boundary Conditions , 2005 .

[7]  R. Kaiser,et al.  On the Existence of Force-Free Magnetic Fields¶with Small Nonconstant α in Exterior Domains , 2000 .

[8]  T. Sakurai Calculation of force-free magnetic field with non-constant α , 1981 .

[9]  S. Régnier,et al.  3D magnetic configuration of the Hα filament and X-ray sigmoid in NOAA AR 8151 , 2004 .

[10]  G. A. Gary,et al.  On the numerical computation of nonlinear force-free magnetic fields , 1985 .

[11]  B. Low,et al.  Modeling solar force-free magnetic fields , 1990 .

[12]  Zoran Mikic,et al.  RECONSTRUCTING THE SOLAR CORONAL MAGNETIC FIELD AS A FORCE-FREE MAGNETIC FIELD , 1997 .

[13]  Michael S. Wheatland,et al.  An Optimization Approach to Reconstructing Force-free Fields , 1997 .

[14]  R. Keppens,et al.  Extrapolation of a nonlinear force-free field containing a highly twisted magnetic loop , 2005 .

[15]  J. Aly,et al.  On the reconstruction of the nonlinear force-free coronal magnetic field from boundary data , 1989 .

[16]  H. Brezis Analyse fonctionnelle : théorie et applications , 1983 .

[17]  T. Amari,et al.  On the existence of non-linear force-free fields in three-dimensional domains , 2000 .

[18]  R. Chodura,et al.  A 3D Code for MHD Equilibrium and Stability , 1981 .

[19]  M. Bineau On the existence of force-free magnetic fields , 1972 .

[20]  T. Amari,et al.  A finite-element method for computing nonlinear force-free fields , 2001 .