Applying the Lagrange multiplier technique to reconstruct a force-free magnetic field

被引:3
作者
Fatholahzadeh, S. [1 ]
Jafarpour, M. H. [1 ]
Nasiri, S. [1 ]
机构
[1] Shahid Beheshti Univ, Dept Phys, Tehran 1983969411, Iran
关键词
Sun: corona; Sun: magnetic fields; Sun: activity; Sun: atmosphere; Sun: flares; methods: data analysis; ACTIVE-REGION; PLASMA BETA; EXTRAPOLATION; EQUILIBRIA;
D O I
10.1051/0004-6361/202243471
中图分类号
P1 [天文学];
学科分类号
0704 ;
摘要
Context. The magnetic field in the solar corona can be reconstructed by extrapolating the data obtained by measuring the magnetic field in the photosphere. It is widely thought that the dynamics of the solar corona is governed by the magnetic field. The magnetic field in the corona can therefore be reconstructed using the physics of force-free fields. Several models have been developed so far to reconstruct the magnetic field, but shortcomings in this regard remain. Therefore, alternative models in this respect can still be proposed. Here we apply the Lagrange multiplier technique to render an optimization model.Aims. The main aim of the present paper is to propose a method of constrained optimization using the vector Lagrange multiplier and compare the results with those of preexisting models.Methods. Our main focus is on the conceptual modification of the optimization models. Since these models are computationally efficient and easy to implement, any possible progressive step would be welcome. The Lagrange multiplier technique is a powerful mathematical tool that has been successfully applied to many areas in physics. It may serve this purpose. In the absence of nonmagnetic forces such as pressure, gravity, and dissipative forces, the coronal medium is dominated by magnetic force. Thus, the function that is considered to be minimized may include a divergence term subject to the constraint force-free term, which yields a solenoidal and force-free magnetic field by an iterative process.Results. The numerical analysis of the proposed model was conducted through an artificial magnetogram and an observational vector magnetogram obtained by SDO/HMI images. The results obtained confirm that (i) the Lagrangian to be minimized in the present model converges slightly faster toward zero, at least for initial iteration steps, (ii) the energy variation during the optimization is compatible with the variations in previous studies, (iii) the numerical results seem to be compatible with a semi-analytical solution as the test case, and (iv) the model is applied to a real magnetogram, and relevant quantities such as magnetic energy content, the current-weighted angle between the current density vector and magnetic field, and the fractional flux errors are computed.Conclusions. The methods and techniques that convert a constrained problem into an unconstrained one are the penalty method, the barrier method, the augmented Lagrangian method, and the Lagrange multiplier technique, for instance. We have employed the Lagrange multiplier technique, by which any constrained condition could be added to the Lagrangian by an appropriate Lagrange multiplier. In our case, the constraint is the force-freeness of the magnetic field and is therefore a special case. The method has the following advantages: (i) The convergence rate is slightly higher for the initial iteration steps, which may help us for time-series events, while several magnetograms must be considered and limited steps of iteration may be needed. (ii) The Lagrangian is introduced based on the Lagrange multiplier technique, which facilitates first fixing a physical compromise such as a divergence-free condition, and subsequently adding any given constraint term. (iii) The quantities obtained by the constrained optimization with the vector Lagrange multiplier model, that is, the relative magnetic energy, the ratio of the magnetic energy to its initial value, the angle between the electric current and the magnetic field, fractional flux errors, the normalized vector error, the vector correlation, and the Cauchy-Schwarz indicator, are comparable with those of the comparison models considered in this article.
引用
收藏
页数:11
相关论文
共 50 条
  • [31] A fluid dynamics approach for the computation of nonlinear force-free magnetic field
    Li, JQ
    Wang, JX
    Wei, FS
    CHINESE JOURNAL OF ASTRONOMY AND ASTROPHYSICS, 2003, 3 (03): : 247 - 256
  • [32] A SELF-CONSISTENT NONLINEAR FORCE-FREE SOLUTION FOR A SOLAR ACTIVE REGION MAGNETIC FIELD
    Wheatland, M. S.
    Reginer, S.
    ASTROPHYSICAL JOURNAL LETTERS, 2009, 700 (02) : L88 - L91
  • [33] Non-constant-α force-free field of active region NOAA 8210
    Régnier, S
    Amari, T
    Canfield, RC
    SOLMAG 2002: PROCEEDINGS OF THE MAGNETIC COUPLING OF THE SOLAR ATMOSPHERE EUROCONFERENCE AND IAU COLLOQUIUM 188, 2002, 505 : 65 - 68
  • [34] MODELING MAGNETIC FIELD STRUCTURE OF A SOLAR ACTIVE REGION CORONA USING NONLINEAR FORCE-FREE FIELDS IN SPHERICAL GEOMETRY
    Guo, Y.
    Ding, M. D.
    Liu, Y.
    Sun, X. D.
    DeRosa, M. L.
    Wiegelmann, T.
    ASTROPHYSICAL JOURNAL, 2012, 760 (01)
  • [35] Nonlinear Force-free Modeling of Flare-related Magnetic Field Changes at the Photosphere and Chromosphere
    Kleint, Lucia
    Wheatland, Michael S.
    Mastrano, Alpha
    McCauley, Patrick I.
    ASTROPHYSICAL JOURNAL, 2018, 865 (02)
  • [36] ON THE NUMERICAL COMPUTATION OF NONLINEAR FORCE-FREE MAGNETIC-FIELDS
    WU, ST
    SUN, MT
    CHANG, HM
    HAGYARD, MJ
    GARY, GA
    ASTROPHYSICAL JOURNAL, 1990, 362 (02) : 698 - 708
  • [37] Reconstructing spherical nonlinear force-free field in the solar corona
    Song, M. T.
    Fang, C.
    Zhang, H. Q.
    Tang, Y. H.
    Wu, S. T.
    Zhang, Y. A.
    ASTROPHYSICAL JOURNAL, 2007, 666 (01) : 491 - 500
  • [38] Extrapolation of a nonlinear force-free field containing a highly twisted magnetic loop
    Valori, G
    Kliem, B
    Keppens, R
    ASTRONOMY & ASTROPHYSICS, 2005, 433 (01): : 335 - 347
  • [39] A Nonlinear Force-Free Magnetic Field Approximation Suitable for Fast Forward-Fitting to Coronal Loops. III. The Free Energy
    Aschwanden, Markus J.
    SOLAR PHYSICS, 2013, 287 (1-2) : 369 - 389
  • [40] HOW MUCH ENERGY CAN BE STORED IN A 3-DIMENSIONAL FORCE-FREE MAGNETIC-FIELD
    ALY, JJ
    ASTROPHYSICAL JOURNAL, 1991, 375 (02) : L61 - L64