LSFEM Implementation of MHD Numerical Solver
- 1 University J. E. Purkinje, ústí nad Labem, Czech Republic
- 2 Astronomical Institute of Czech Academy of Sciences, Ond?ejov, Czech Republic
Abstract
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.
- E. R. Priest, “Solar Magneto-Hydrodynamics,” Reidel, Dordrecht, 1984.
- K. Shibata and S. Tanuma, “Plasmoid-Induced-Reconnection and Fractal Reconnection,” Earth, Planets, and Space, Vol. 53, 2001, pp. 473-482. http://adsabs.harvard.edu/abs/2001EP%26S...53..473S
- M. Bárta, J. Büchner, M. Karlicky and J. Skála, “Spontaneous Current-Layer Fragmentation and Cascading Reconnection in Solar Flares. I. Model and Analysis,” Astrophysical Journal, Vol. 737, No. 1, 2011. doi:10.1088/0004-637X/737/1/24
- M. Bárta, J. Büchner, M. Karlicky and P. Kotr?, “Spontaneous Current-Layer Fragmentation and Cascading Reconnection in Solar Flares. II. Relation to Observations,” Astrophysical Journal, Vol. 730, No. 1, 2011, 6 pages. doi:10.1088/0004-637X/730/1/47
- O. Mátais, “Course 3: Large-Eddy Simulations of Turbulence,” In: M. Lesieur, A. Yaglom and F. David, Eds., New Trends in Turbulence, Vol. 730, No. 1, 2001, pp. 113-186.
- M. A. Leschziner, “Course 4: Statistical Turbulence Modelling for the Computation of Physically Complex Flows,” In: M. Lesieur, A. Yaglom and F. David, Eds., New Trends in Turbulence, Springer Verlag, Berlin, 2001, p. 187. doi:10.1007/3-540-45674-0_4
- M. J. Berger and J. Oliger, “Adaptive Mesh Refinement for Hyperbolic Partial Differential Equations,” Journal of Computational Physics, Vol. 53, No. 3, 1984, pp. 484-512. doi:10.1016/0021-9991(84)90073-1
- B. Fryxell, K. Olson, P. Ricker, F. X. Timmes, M. Zingale, D. Q. Lamb, P. MacNeice, R. Rosner, J. W. Truran, and H. Tufo, “FLASH: An Adaptive Mesh Hydrodynamics Code for Modeling Astrophysical Thermonuclear Flashes,” Astrophysical Journal Supplements, Vol. 131, No. 1, 2000, pp. 273-334.
- B. van der Holst and R. Keppens, “Hybrid Block-AMR in Cartesian and Curvilinear Coordinates: MHD Applications,” Journal of Computational Physics, Vol. 226, No. 1, 2007, pp. 925-946.
- B. Jiang, “The Least-Squares Finite Element Method,” Springer-Verlag, Berlin, 1998.
- P. P. Bochev and M. D. Gunzburger, “Least-Squares Finite Element Methods,” Springer Science+Business Media, New York, 2009.
- C. R. Sovinec, A. H. Glasser, T. A. Gianakon, D. C. Barnes, R. A. Nebel, S. E. Kruger, D. D. Schnack, S. J. Plimpton, A. Tarditi and M. S. Chu, “Nonlinear Magnetohydrodynamics Simulation Using High-Order Finite Elements,” Journal of Computational Physics, Vol. 195, No. 1, 2004, pp. 355-386. doi:10.1016/j.jcp.2003.10.004