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

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.Comment: 19 pages, 10 figures, 1 table. Submitted to New Journal of Physic

Approximation of Parametric Derivatives by the Empirical Interpolation Method

We introduce a general a priori convergence result for the approximation of parametric derivatives of parametrized functions. We consider the best approximations to parametric derivatives in a sequence of approximation spaces generated by a general approximation scheme, and we show that these approximations are convergent provided that the best approximation to the function itself is convergent. We also provide estimates for the convergence rates. We present numerical results with spaces generated by a particular approximation scheme—the Empirical Interpolation Method—to confirm the validity of the general theory