Stability analysis and error estimates of an exactly divergence-free method for the magnetic induction equations^∗

¹ Department of Radiology, Stanford University, Stanford, CA 94305-5105, USA.
hyang2@stanford.edu. Present address: Department of Mathematics, Ohio State University, Columbus, OH 43210-1174, USA. yang.1671@osu.edu
² Department of Mathematical Sciences, Rensselaer Polytechnic Institute, Troy, NY 12180-3590, USA.
lif@rpi.edu

Received: 23 March 2015
Revised: 7 July 2015

Abstract

In this paper, we consider an exactly divergence-free scheme to solve the magnetic induction equations. This problem is motivated by the numerical simulations of ideal magnetohydrodynamic (MHD) equations, a nonlinear hyperbolic system with a divergence-free condition on the magnetic field. Computational methods without satisfying such condition may lead to numerical instability. One class of methods, constrained transport schemes, is widely used as divergence-free treatments. So far there is not much analysis available for such schemes. In this work, we take an exactly divergence-free scheme proposed by [Li and Xu, J. Comput. Phys. 231 (2012) 2655–2675] as a candidate of the constrained transport schemes, and adapt it to solve the magnetic induction equations. For the resulting scheme applied to the equations with a constant velocity field, we carry out von Neumann analysis for numerical stability on uniform meshes. We also establish the stability and error estimates based on energy methods. In particular, we identify the stability mechanism due to the spatial and temporal discretizations, and the role of the exactly divergence-free property of the numerical solution for stability. The analysis based on energy methods can be extended to non-uniform meshes, and they can also be applied to the magnetic induction equations with a variable velocity field, which is more relevant to the MHD simulations.