A comparison of spectral element and finite difference methods using statically refined nonconforming grids for the MHD island coalescence instability problem

Abstract

A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)] is applied to simulate the problem of MHD island coalescence instability (MICI) in two dimensions. MICI is a fundamental MHD process that can produce sharp current layers and subsequent reconnection and heating in a high-Lundquist number plasma such as the solar corona [Ng and Bhattacharjee, Phys. Plasmas, 5, 4028 (1998)]. Due to the formation of thin current layers, it is highly desirable to use adaptively or statically refined grids to resolve them, and to maintain accuracy at the same time. The output of the spectral-element static adaptive refinement simulations are compared with simulations using a finite difference method on the same refinement grids, and both methods are compared to pseudo-spectral simulations with uniform grids as baselines. It is shown that with the statically refined grids roughly scaling linearly with effective resolution, spectral element runs can maintain accuracy significantly higher than that of the finite difference runs, in some cases achieving close to full spectral accuracy.

19 pages, 17 figures, submitted to Astrophys. J. Suppl

References in corpus (6)

• On Solving the Coronal Heating Problem
• A Chandra X-ray Investigation of the Violent Interstellar Medium: From Dwarf Starbursts to Ultraluminous Infrared Galaxies
• A Comparison of Two Shallow Water Models with Non-Conforming Adaptive Grids: classical tests
• Geophysical-astrophysical spectral-element adaptive refinement (GASpAR): Object-oriented h-adaptive code for geophysical fluid dynamics simulation
• Adaptive mesh refinement with spectral accuracy for magnetohydrodynamics in two space dimensions
• A Constrained Tectonics Model for Coronal Heating

Cited by in corpus (8)

• Reconnection and particle acceleration in interacting flux ropes -- II. 3D effects on test particles in magnetically dominated plasmas
• High-Lundquist Number Scaling in Three-Dimensional Simulations of Parker's Model of Coronal Heating
• Interacting tilt and kink instabilities in repelling current channels
• A paradigmatic flow for small-scale magnetohydrodynamics: properties of the ideal case and the collision of current sheets
• FINMHD: an adaptive finite element code for magnetic reconnection and plasmoid chains formation in Magnetohydrodynamics
• Formation of plasmoid chains and fast magnetic reconnection during nonlinear evolution of the tilt instability
• Petschek-type reconnection in the high-Lundquist-number regime during nonlinear evolution on the tilt instability
• The role of field correlations on turbulent dissipation