Author(s)

Jan Skála, Miroslav Bárta

Astronomical Institute of Czech Academy of Sciences, Ond?ejov, Czech Republic.
University J. E. Purkinje, ústí nad Labem, Czech Republic.

Many problems in physics are inherently of multi-scale nature. The issues of MHD turbulence or magnetic reconnection, namely in the hot and sparse, almost collision-less astrophysical plasmas, can stand as clear examples. The Finite Element Method (FEM) with adaptive gridding appears to be the appropriate numerical implementation for handling the broad range of scales contained in such high Lundquist-number MHD problems. In spite the FEM is now routinely used in engineering practice in solid-state and fluid dynamics, its usage for MHD simulations has recently only begun and only few implementations exist so far. In this paper we present our MHD solver based on the Least-Square FEM (LSFEM) formulation. We describe the transformation of the MHD equations into form required for finding the LSFEM functional and some practical issues in implementation of the method. The algorithm was tested on selected problems of ideal (non-resistive) and resistive MHD. The tests show the usability of LSFEM for solving MHD equations.

Magnetohydrodynamics (MHD); Least-Squares Finite Element Method; Adaptive Mesh Refinement; Magnetic Reconnection; Solar Eruptions; MHD Turbulence

Skála, J. and Bárta, M. (2012) LSFEM Implementation of MHD Numerical Solver. Applied Mathematics, 3, 1842-1850. doi: 10.4236/am.2012.331250.