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

Yang, He; Li, Fengyan

ESAIM: Mathematical Modelling and Numerical Analysis , Tome 50 (2016) no. 4, pp. 965-993.

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.

D.S. Balsara, Divergence-free reconstruction of magnetic fields and WENO schemes for magne tohydrodynamics. J. Comput. Phys. 228 (2009) 5040–5056. | DOI | MR | Zbl

D.S. Balsara and D.S. Spicer, A staggered mesh algorithm using high order Godunov fluxes to ensure solenoidal magnetic fields in magnetohydrodynamic simulations. J. Comput. Phys. 149 (1999) 270–292. | DOI | MR | Zbl

T. Barth, On the role of involutions in the discontinuous Galerkin discretization of Maxwell and magnetohydrodynamic systems, The IMA Volumes in Mathematics and its Applications, IMA volume on Compatible Spatial Discretizations, edited by Arnold, Bochev and Shashkov 142 (2005) 69–88. | MR | Zbl

N. Besse and D. Kröner, Convergence of the locally divergence free discontinuous Galerkin methods for induction equations for the 2D-MHD system. ESAIM: M2AN 39 (2005) 1117–1202. | DOI | Numdam | MR | Zbl

J.U. Brackbill and D.C. Barnes, The effect of nonzero $\nabla ·𝐁$ on the numerical solution of the magnetohydrodynamic equations. J. Comput. Phys. 35 (1980) 426–430. | DOI | MR | Zbl

F. Brezzi, J. Douglas, Jr. and L.D. Marini, Two Families of Mixed Finite Elements for Second Order Elliptic Problems. Numer. Math. 47 (1985) 217–235. | DOI | MR | Zbl

B. Cockburn, An introduction to the discontinuous Galerkin method for convection-dominated problems. Advanced Numerical Approximation of Nonlinear Hyperbolic Equations. Vol. 169 of Lect. Notes Math. (1998) 151–268. | MR | Zbl

C.R. Evans and J.F. Hawley, Simulation of magnetohydrodynamic flows: a constrained transport method. Astrophys. J. 322 (1988) 659–677. | DOI

F.G. Fuchs, K.H. Karlsen, S. Mishra and N.H. Risebro, Stable upwind schemes for the magnetic induction equation. ESAIM: M2AN 43 (2009) 825–852. | DOI | Numdam | MR | Zbl

T.A. Gardiner and J.M. Stone, An unsplit Godunov method for ideal MHD via constrained transport. J. Comput. Phys. 205 (2005) 509–539. | DOI | MR | Zbl

S.K. Godunov, Symmetric form of the equations of magnetohydrodynamics. Numer. Methods Mech. Continuum Medium 1 (1972) 26–34.

S. Li, High order central scheme on overlapping cells for magnetohydrodynamic flows with and without constrained transport method. J. Comput. Phys. 227 (2008) 7368–7393. | DOI | MR | Zbl

F. Li and L. Xu, Arbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations. J. Comput. Phys. 231 (2012) 2655–2675. | DOI | MR | Zbl

F. Li, L. Xu and S. Yakovlev, Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field. J. Comput. Phys. 230 (2011) 4828–4847. | DOI | MR | Zbl

Y. Liu, C.-W. Shu, E. Tadmor and M. Zhang, Central discontinuous Galerkin methods on overlapping cells with a nonoscillatory hierarchical reconstruction. SIAM J. Numer. Anal. 45 (2007) 2442–2467. | DOI | MR | Zbl

Y. Liu, C.-W. Shu, E. Tadmor and M. Zhang, $L^2$ stability analysis of the central discontinuous Galerkin method and a comparison between the central and regular discontinuous Galerkin methods. ESAIM: M2AN 42 (2008) 593–607. | DOI | Numdam | MR | Zbl

S. Mishra and M. Svärd, On stability of numerical schemes via frozen coefficients and the magnetic induction equations. BIT Num. Math. 50 (2010) 85–108. | DOI | MR | Zbl

K.W. Morton and R.L. Roe, Vorticity-preserving Lax-Wendroff-type schemes for the system wave equation. SIAM J. Sci. Comput. 23 (2001) 170–192. | DOI | MR | Zbl

P.A. Raviart and J.M. Thomas, A mixed finite element method for 2nd order elliptic problems, Mathematical Aspects of Finite Element Methods. Springer (1997) 292–315. | MR | Zbl

J.A. Rossmanith, An unstaggered, high-resolution constrained transport method for magnetohydrodynamic flows. SIAM J. Sci. Comput. 28 (2006) 1766–1797. | DOI | MR | Zbl

J.A. Rossmanith, High-order discontinuous Galerkin finite element methods with globally divergence-free constrained transport for ideal MHD. Preprint (2013).

G. Tóth, The $\nabla ·𝐁=0$ constraint in shock-capturing magnetohydrodynamics codes. J. Comput. Phys. 161 (2000) 605–652. | DOI | MR | Zbl

Q. Zhang and C-W. Shu, Error estimates to smooth solutions of Runge−Kutta discontinuous Galerkin methods for scalar conservation laws. SIAM J. Numer. Anal. 42 (2004) 641–666. | DOI | MR | Zbl

Q. Zhang and C-W. Shu, Stability analysis and a priori error estimates to the third order explicit Runge−Kutta discontinuous Galerkin method for scalar conservation laws, SIAM J. Numerical Anal. 48 (2010) 1038–1063. | DOI | MR | Zbl