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

He Yang, Fengyan Li

Research output: Contribution to journalArticlepeer-review

14 Scopus citations

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.

Original languageEnglish (US)
Pages (from-to)965-993
Number of pages29
JournalESAIM: Mathematical Modelling and Numerical Analysis
Volume50
Issue number4
DOIs
StatePublished - 2016

Keywords

  • Constrained transport
  • Discontinuous galerkin
  • Divergence-free
  • Error estimates
  • Ideal magnetohydrodynamic (MHD) equations
  • Magnetic induction equations
  • Stability

ASJC Scopus subject areas

  • Analysis
  • Numerical Analysis
  • Modeling and Simulation
  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Stability analysis and error estimates of an exactly divergence-free method for the magnetic induction equations'. Together they form a unique fingerprint.

Cite this