Arbitrary order exactly divergence-free central discontinuous Galerkin methods for ideal MHD equations

Ideal magnetohydrodynamic (MHD) equations consist of a set of nonlinear hyperbolic conservation laws, with a divergence-free constraint on the magnetic field. Neglecting this constraint in the design of computational methods may lead to numerical instability or nonphysical features in solutions. In our recent work [F. Li, L. Xu, S. Yakovlev, Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field, Journal of Computational Physics 230 (2011) 4828-4847], second and third order exactly divergence-free central discontinuous Galerkin methods were proposed for ideal MHD equations. In this paper, we further develop such methods with higher order accuracy. The novelty here is that the well-established H(div)-conforming finite element spaces are used in the constrained transport type framework, and the magnetic induction equations are extensively explored in order to extract sufficient information to uniquely reconstruct an exactly divergence-free magnetic field. The overall algorithm is local, and it can be of arbitrary order of accuracy. Numerical examples are presented to demonstrate the performance of the proposed methods especially when they are fourth order accurate.

[1]  Chi-Wang Shu.  TVB uniformly high-order schemes for conservation laws , 1987 .

[2]  C. Munz,et al.  Hyperbolic divergence cleaning for the MHD equations , 2002 .

[3]  Chi-Wang Shu,et al.  Strong Stability-Preserving High-Order Time Discretization Methods , 2001, SIAM Rev..

[4]  L. D. Marini,et al.  Two families of mixed finite elements for second order elliptic problems , 1985 .

[5]  Chi-Wang Shu,et al.  L 2 STABILITY ANALYSIS OF THE CENTRAL DISCONTINUOUS GALERKIN METHOD AND A COMPARISON BETWEEN THE CENTRAL AND REGULAR , 2008 .

[6]  J. Stone,et al.  An unsplit Godunov method for ideal MHD via constrained transport , 2005, astro-ph/0501557.

[7]  J. Brackbill,et al.  The Effect of Nonzero ∇ · B on the numerical solution of the magnetohydrodynamic equations☆ , 1980 .

[8]  Chi-Wang Shu,et al.  The Runge-Kutta Discontinuous Galerkin Method for Conservation Laws V , 1998 .

[9]  J. T..  ON THE ROLE OF INVOLUTIONS IN THE DISCONTINUOUS GALERKIN DISCRETIZATION OF MAXWELL AND MAGNETOHYDRODYNAMIC SYSTEMS , 2005 .

[10]  Mark Ainsworth,et al.  Dispersive and dissipative behaviour of high order discontinuous Galerkin finite element methods , 2004 .

[11]  Chi-Wang Shu,et al.  Runge-Kutta Discontinuous Galerkin Method Using WENO Limiters , 2005, SIAM J. Sci. Comput..

[12]  Sergey Yakovlev,et al.  Central discontinuous Galerkin methods for ideal MHD equations with the exactly divergence-free magnetic field , 2011, J. Comput. Phys..

[13]  Dinshaw Balsara,et al.  Second-Order-accurate Schemes for Magnetohydrodynamics with Divergence-free Reconstruction , 2003, astro-ph/0308249.

[14]  J. Hawley,et al.  Simulation of magnetohydrodynamic flows: A Constrained transport method , 1988 .

[15]  P. Raviart,et al.  A mixed finite element method for 2-nd order elliptic problems , 1977 .

[16]  Eitan Tadmor,et al.  Non-Oscillatory Hierarchical Reconstruction for Central and Finite Volume Schemes , 2006 .

[17]  Chi-Wang Shu,et al.  L2 Stability Analysis of the Central Discontinuous Galerkin Method and a Comparison between the Central and Regular Discontinuous Galerkin Methods , 2008 .

[18]  Chi-Wang Shu,et al.  Locally Divergence-Free Discontinuous Galerkin Methods for MHD Equations , 2005, J. Sci. Comput..

[19]  Dinshaw S. Balsara.  Divergence-free reconstruction of magnetic fields and WENO schemes for magnetohydrodynamics , 2009, J. Comput. Phys..

[20]  Chi-Wang Shu,et al.  Central Discontinuous Galerkin Methods on Overlapping Cells with a Nonoscillatory Hierarchical Reconstruction , 2007, SIAM J. Numer. Anal..

[21]  Shengtai Li,et al.  A fourth-order divergence-free method for MHD flows , 2010, J. Comput. Phys..

[22]  Kenneth G. Powell,et al.  AN APPROXIMATE RIEMANN SOLVER FOR MAGNETOHYDRODYNAMICS (That Works in More than One Dimension) , 1994 .

[23]  T. Barth.  On the Role of Involutions in the Discontinuous Galerkin Discretization of Maxwell and Magnetohydrodynamic Systems , 2013 .

[24]  Dinshaw Balsara,et al.  Divergence-free adaptive mesh refinement for Magnetohydrodynamics , 2001 .

[25]  Shengtai Li,et al.  High order central scheme on overlapping cells for magneto-hydrodynamic flows with and without constrained transport method , 2008, J. Comput. Phys..

[26]  D. Balsara,et al.  A Staggered Mesh Algorithm Using High Order Godunov Fluxes to Ensure Solenoidal Magnetic Fields in Magnetohydrodynamic Simulations , 1999 .

[27]  Chi-Wang Shu,et al.  Runge–Kutta Discontinuous Galerkin Methods for Convection-Dominated Problems , 2001, J. Sci. Comput..

[28]  Sergey Yakovlev,et al.  Locally divergence-free central discontinuous Galerkin methods for ideal MHD equations , 2013, J. Comput. Sci..

[29]  Paul R. Woodward,et al.  Extension of the Piecewise Parabolic Method to Multidimensional Ideal Magnetohydrodynamics , 1994 .

[30]  Chi-Wang Shu,et al.  Locally divergence-free discontinuous Galerkin methods for the Maxwell equations , 2004, Journal of Computational Physics.

[31]  G. Tóth.  The ∇·B=0 Constraint in Shock-Capturing Magnetohydrodynamics Codes , 2000 .

[32]  M. Fortin,et al.  E cient rectangular mixed fi-nite elements in two and three space variables , 1987 .