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.
DOI : 10.1051/m2an/2015061
Mots clés : Stability, error estimates, ideal magnetohydrodynamic (MHD) equations, constrained transport, divergence-free, discontinuous galerkin, magnetic induction equations
@article{M2AN_2016__50_4_965_0, author = {Yang, He and Li, Fengyan}, title = {Stability analysis and error estimates of an exactly divergence-free method for the magnetic induction equations}, journal = {ESAIM: Mathematical Modelling and Numerical Analysis }, pages = {965--993}, publisher = {EDP-Sciences}, volume = {50}, number = {4}, year = {2016}, doi = {10.1051/m2an/2015061}, zbl = {1348.78028}, mrnumber = {3521708}, language = {en}, url = {http://www.numdam.org/articles/10.1051/m2an/2015061/} }
TY - JOUR AU - Yang, He AU - Li, Fengyan TI - Stability analysis and error estimates of an exactly divergence-free method for the magnetic induction equations JO - ESAIM: Mathematical Modelling and Numerical Analysis PY - 2016 SP - 965 EP - 993 VL - 50 IS - 4 PB - EDP-Sciences UR - http://www.numdam.org/articles/10.1051/m2an/2015061/ DO - 10.1051/m2an/2015061 LA - en ID - M2AN_2016__50_4_965_0 ER -
%0 Journal Article %A Yang, He %A Li, Fengyan %T Stability analysis and error estimates of an exactly divergence-free method for the magnetic induction equations %J ESAIM: Mathematical Modelling and Numerical Analysis %D 2016 %P 965-993 %V 50 %N 4 %I EDP-Sciences %U http://www.numdam.org/articles/10.1051/m2an/2015061/ %R 10.1051/m2an/2015061 %G en %F M2AN_2016__50_4_965_0
Yang, He; Li, Fengyan. Stability analysis and error estimates of an exactly divergence-free method for the magnetic induction equations. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 50 (2016) no. 4, pp. 965-993. doi : 10.1051/m2an/2015061. http://www.numdam.org/articles/10.1051/m2an/2015061/
Divergence-free reconstruction of magnetic fields and WENO schemes for magne tohydrodynamics. J. Comput. Phys. 228 (2009) 5040–5056. | DOI | MR | Zbl
,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
and ,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
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
and ,The effect of nonzero on the numerical solution of the magnetohydrodynamic equations. J. Comput. Phys. 35 (1980) 426–430. | DOI | MR | Zbl
and ,Two Families of Mixed Finite Elements for Second Order Elliptic Problems. Numer. Math. 47 (1985) 217–235. | DOI | MR | Zbl
, and ,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
Simulation of magnetohydrodynamic flows: a constrained transport method. Astrophys. J. 322 (1988) 659–677. | DOI
and ,Stable upwind schemes for the magnetic induction equation. ESAIM: M2AN 43 (2009) 825–852. | DOI | Numdam | MR | Zbl
, , and ,An unsplit Godunov method for ideal MHD via constrained transport. J. Comput. Phys. 205 (2005) 509–539. | DOI | MR | Zbl
and ,Symmetric form of the equations of magnetohydrodynamics. Numer. Methods Mech. Continuum Medium 1 (1972) 26–34.
,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
,Arbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations. J. Comput. Phys. 231 (2012) 2655–2675. | DOI | MR | Zbl
and ,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
, and ,Central discontinuous Galerkin methods on overlapping cells with a nonoscillatory hierarchical reconstruction. SIAM J. Numer. Anal. 45 (2007) 2442–2467. | DOI | MR | Zbl
, , and ,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
, , and ,On stability of numerical schemes via frozen coefficients and the magnetic induction equations. BIT Num. Math. 50 (2010) 85–108. | DOI | MR | Zbl
and ,Vorticity-preserving Lax-Wendroff-type schemes for the system wave equation. SIAM J. Sci. Comput. 23 (2001) 170–192. | DOI | MR | Zbl
and ,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
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).
The constraint in shock-capturing magnetohydrodynamics codes. J. Comput. Phys. 161 (2000) 605–652. | DOI | MR | Zbl
,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
and ,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
and ,Cité par Sources :