Comparison of methods for modelling coronal magnetic fields
File(s)erin_fields.pdf (376.34 KB)
Published version
Author(s)
Goldstraw, EE
Hood, AW
Browning, PK
Cargill, PJ
Type
Journal Article
Abstract
Aims. Four different approximate approaches used to model the stressing of coronal magnetic fields due to an imposed photospheric motion are compared with each other and the results from a full time-dependent magnetohydrodynamic (MHD) code. The assumptions used for each of the approximate methods are tested by considering large photospheric footpoint displacements.
Methods. We consider a simple model problem, comparing the full non-linear MHD, determined with the Lare2D numerical code, with four approximate approaches. Two of these, magneto-frictional relaxation and a quasi-1D Grad-Shafranov approach, assume sequences of equilibria, whilst the other two methods, a second-order linearisation of the MHD equations and Reduced MHD, are time dependent.
Results. The relaxation method is very accurate compared to full MHD for force-free equilibria for all footpoint displacements, but has significant errors when the plasma β0 is of order unity. The 1D approach gives an extremely accurate description of the equilibria away from the photospheric boundary layers, and agrees well with Lare2D for all parameter values tested. The linearised MHD equations correctly predict the existence of photospheric boundary layers that are present in the full MHD results. As soon as the footpoint displacement becomes a significant fraction of the loop length, the RMHD method fails to model the sequences of equilibria correctly. The full numerical solution is interesting in its own right, and care must be taken for low β0 plasmas if the viscosity is too high.
Methods. We consider a simple model problem, comparing the full non-linear MHD, determined with the Lare2D numerical code, with four approximate approaches. Two of these, magneto-frictional relaxation and a quasi-1D Grad-Shafranov approach, assume sequences of equilibria, whilst the other two methods, a second-order linearisation of the MHD equations and Reduced MHD, are time dependent.
Results. The relaxation method is very accurate compared to full MHD for force-free equilibria for all footpoint displacements, but has significant errors when the plasma β0 is of order unity. The 1D approach gives an extremely accurate description of the equilibria away from the photospheric boundary layers, and agrees well with Lare2D for all parameter values tested. The linearised MHD equations correctly predict the existence of photospheric boundary layers that are present in the full MHD results. As soon as the footpoint displacement becomes a significant fraction of the loop length, the RMHD method fails to model the sequences of equilibria correctly. The full numerical solution is interesting in its own right, and care must be taken for low β0 plasmas if the viscosity is too high.
Date Issued
2018-02-26
Date Acceptance
2017-11-03
Citation
Astronomy and Astrophysics, 2018, 610
ISSN
0004-6361
Publisher
EDP Sciences
Journal / Book Title
Astronomy and Astrophysics
Volume
610
Copyright Statement
Reproduced with permission from Astronomy & Astrophysics © ESO 2018
Identifier
http://gateway.webofknowledge.com/gateway/Gateway.cgi?GWVersion=2&SrcApp=PARTNER_APP&SrcAuth=LinksAMR&KeyUT=WOS:000426130200002&DestLinkType=FullRecord&DestApp=ALL_WOS&UsrCustomerID=1ba7043ffcc86c417c072aa74d649202
Subjects
Science & Technology
Physical Sciences
Astronomy & Astrophysics
Sun: magnetic fields
Sun: corona
magnetohydrodynamics (MHD)
LARGE-SCALE CORONA
REDUCED MAGNETOHYDRODYNAMICS
SYNTHETIC MAGNETOGRAMS
SIMULATIONS
TURBULENCE
EQUATIONS
EVOLUTION
LOOPS
Publication Status
Published
Article Number
ARTN A48