Adaptive mesh refinement with spectral accuracy for magnetohydrodynamics in two space dimensions

Abstract

We examine the effect of accuracy of high-order spectral element methods, with or without adaptive mesh refinement (AMR), in the context of a classical configuration of magnetic reconnection in two space dimensions, the so-called Orszag-Tang vortex made up of a magnetic X-point centered on a stagnation point of the velocity. A recently developed spectral-element adaptive refinement incompressible magnetohydrodynamic (MHD) code is applied to simulate this problem. The MHD solver is explicit, and uses the Elsasser formulation on high-order elements. It automatically takes advantage of the adaptive grid mechanics that have been described elsewhere in the fluid context [Rosenberg, Fournier, Fischer, Pouquet, J. Comp. Phys. 215, 59-80 (2006)]; the code allows both statically refined and dynamically refined grids. Tests of the algorithm using analytic solutions are described, and comparisons of the Orszag-Tang solutions with pseudo-spectral computations are performed. We demonstrate for moderate Reynolds numbers that the algorithms using both static and refined grids reproduce the pseudo--spectral solutions quite well. We show that low-order truncation--even with a comparable number of global degrees of freedom--fails to correctly model some strong (sup--norm) quantities in this problem, even though it satisfies adequately the weak (integrated) balance diagnostics.

19 pages, 10 figures, 1 table. Submitted to New Journal of Physics

References in corpus (11)

• Generation of magnetic field by dynamo action in a turbulent flow of liquid sodium
• Magnetic Field Saturation in the Riga Dynamo Experiment
• The imprint of large-scale flows on turbulence
• Small scale structures in three-dimensional magnetohydrodynamic turbulence
• Adaptive Mesh Refinement for Singular Current Sheets in Incompressible Magnetohydrodynamic Flows
• Observation of a Turbulence-Induced Large Scale Magnetic Field
• Numerical solutions of the three-dimensional magnetohydrodynamic alpha-model
• A numerical study of the alpha model for two-dimensional magnetohydrodynamic turbulent flows
• Geophysical-astrophysical spectral-element adaptive refinement (GASpAR): Object-oriented h-adaptive code for geophysical fluid dynamics simulation
• Cancellation exponent and multifractal structure in two-dimensional magnetohydrodynamics: direct numerical simulations and Lagrangian averaged modeling
• Exact calculation of Fourier series in nonconforming spectral-element methods

Cited by in corpus (6)

• Kinetic dissipation and anisotropic heating in a turbulent collisionless plasma
• Ideal evolution of MHD turbulence when imposing Taylor-Green symmetries
• Fourier continuation method for incompressible fluids with boundaries
• A comparison of spectral element and finite difference methods using statically refined nonconforming grids for the MHD island coalescence instability problem
• Vector potential-based MHD solver for non-periodic flows using Fourier continuation expansions
• The role of field correlations on turbulent dissipation