An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part I: Theory and numerical verification

This work presents an extension of discretely entropy stable discontinuous Galerkin (DG) methods to the resistive magnetohydrodynamics (MHD) equations. Although similar to the compressible Navier-Stokes equations at first sight, there are some important differences concerning the resistive MHD equations that need special focus. The continuous entropy analysis of the ideal MHD equations, which are the advective parts of the resistive MHD equations, shows that the divergence-free constraint on the magnetic field components must be incorporated as a non-conservative term in a form either proposed by Powell or Janhunen. Consequently, this non-conservative term needs to be discretized, such that the approximation is consistent with the entropy. As an extension of the ideal MHD system, we address in this work the continuous analysis of the resistive MHD equations and show that the entropy inequality holds. Thus, our first contribution is the proof that the resistive terms are symmetric and positive semi-definite when formulated in entropy space as gradients of the entropy variables. Moreover, this enables the construction of an entropy stable DG discretization for the resistive MHD equations. However, the resulting method suffers from large errors in the divergence-free constraint, since no particular treatment of divergence errors is included in the standard resistive MHD model. Hence, our second contribution is the extension of the resistive MHD equations with proper divergence cleaning based on a generalized Lagrange multiplier (GLM) strategy. We construct and analyze a DG method that is entropy stable for the resistive MHD equations and has a built-in GLM divergence cleaning mechanism. The theoretical derivations and proofs are then verified by several numerical examples...

[1]  M. Farrashkhalvat,et al.  Basic Structured Grid Generation: With an introduction to unstructured grid generation , 2003 .

[2]  P. Janhunen,et al.  A Positive Conservative Method for Magnetohydrodynamics Based on HLL and Roe Methods , 2000 .

[3]  A. Harten.  On the symmetric form of systems of conservation laws with entropy , 1983 .

[4]  David A. Kopriva,et al.  A Polynomial Spectral Calculus for Analysis of DG Spectral Element Methods , 2017, 1704.00709.

[5]  Leland Jameson,et al.  Numerical Convergence Study of Nearly Incompressible, Inviscid Taylor–Green Vortex Flow , 2005, J. Sci. Comput..

[6]  G. Taylor,et al.  Mechanism of the production of small eddies from large ones , 1937 .

[7]  Tatiana G. Elizarova,et al.  Numerical simulation of three-dimensional quasi-neutral gas flows based on smoothed magnetohydrodynamic equations , 2015, Computational Mathematics and Mathematical Physics.

[8]  S. Orszag,et al.  Small-scale structure of two-dimensional magnetohydrodynamic turbulence , 1979, Journal of Fluid Mechanics.

[9]  Tadmor.  Entropy functions for symmetric systems of conservation laws. Final Report , 1987 .

[10]  P. Dutt,et al.  Stable boundary conditions and difference schemas for Navier-Stokes equations , 1988 .

[11]  Gregor Gassner,et al.  The BR1 Scheme is Stable for the Compressible Navier–Stokes Equations , 2017, J. Sci. Comput..

[12]  Deep Ray,et al.  An entropy stable finite volume scheme for the two dimensional Navier-Stokes equations on triangular grids , 2017, Appl. Math. Comput..

[13]  Philip L. Roe,et al.  Affordable, entropy-consistent Euler flux functions II: Entropy production at shocks , 2009, J. Comput. Phys..

[14]  Gregor Gassner,et al.  An entropy stable nodal discontinuous Galerkin method for the two dimensional shallow water equations on unstructured curvilinear meshes with discontinuous bathymetry , 2015, J. Comput. Phys..

[15]  Eitan Tadmor,et al.  PERFECT DERIVATIVES, CONSERVATIVE DIFFERENCES AND ENTROPY STABLE COMPUTATION OF HYPERBOLIC CONSERVATION LAWS , 2016 .

[16]  Michael Dumbser,et al.  Space-time adaptive ADER-DG schemes for dissipative flows: Compressible Navier-Stokes and resistive MHD equations , 2016, Comput. Phys. Commun..

[17]  Gregor Gassner,et al.  Geometry effects in nodal discontinuous Galerkin methods on curved elements that are provably stable , 2016, Appl. Math. Comput..

[18]  Christopher A. Kennedy,et al.  Reduced aliasing formulations of the convective terms within the Navier-Stokes equations for a compressible fluid , 2008, J. Comput. Phys..

[19]  Christian Klingenberg,et al.  Entropy Stable Finite Volume Scheme for Ideal Compressible MHD on 2-D Cartesian Meshes , 2016, SIAM J. Numer. Anal..

[20]  Gregor Gassner,et al.  A Skew-Symmetric Discontinuous Galerkin Spectral Element Discretization and Its Relation to SBP-SAT Finite Difference Methods , 2013, SIAM J. Sci. Comput..

[21]  Matteo Parsani,et al.  Towards an Entropy Stable Spectral Element Framework for Computational Fluid Dynamics , 2016 .

[22]  Gregor Gassner,et al.  A uniquely defined entropy stable matrix dissipation operator for high Mach number ideal MHD and compressible Euler simulations , 2017, J. Comput. Phys..

[23]  Gregor Gassner,et al.  Affordable, entropy conserving and entropy stable flux functions for the ideal MHD equations , 2015, J. Comput. Phys..

[24]  Nail K. Yamaleev,et al.  Discretely conservative finite-difference formulations for nonlinear conservation laws in split form: Theory and boundary conditions , 2013, J. Comput. Phys..

[25]  P. Mininni,et al.  Paradigmatic flow for small-scale magnetohydrodynamics: properties of the ideal case and the collision of current sheets. , 2008, Physical review. E, Statistical, nonlinear, and soft matter physics.

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

[27]  Claus-Dieter Munz,et al.  Explicit Discontinuous Galerkin methods for unsteady problems , 2012 .

[28]  M. Carpenter,et al.  Fourth-order 2N-storage Runge-Kutta schemes , 1994 .

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

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

[31]  Timothy J. Barth,et al.  Numerical Methods for Gasdynamic Systems on Unstructured Meshes , 1997, Theory and Numerics for Conservation Laws.

[32]  Joachim Saur,et al.  Consistent boundary conditions at nonconducting surfaces of planetary bodies: Applications in a new Ganymede MHD model , 2014 .

[33]  George Em Karniadakis,et al.  A Discontinuous Galerkin Method for the Viscous MHD Equations , 1999 .

[34]  Sergio Pirozzoli,et al.  Numerical Methods for High-Speed Flows , 2011 .

[35]  Travis C. Fisher,et al.  High-order entropy stable finite difference schemes for nonlinear conservation laws: Finite domains , 2013, J. Comput. Phys..

[36]  Claus-Dieter Munz,et al.  Mesh Curving Techniques for High Order Discontinuous Galerkin Simulations , 2015 .

[37]  E. Tadmor.  Skew-selfadjoint form for systems of conservation laws , 1984 .

[38]  Chi-Wang Shu,et al.  A discontinuous Galerkin finite element method for directly solving the Hamilton-Jacobi equations , 2007, Journal of Computational Physics.

[39]  L. C. Woods.  An Introduction to the Kinetic Theory of Gases and magnetoplasmas. By L. C. WOODS. Oxford University Press, 1993. 300 pp. £27.50. , 1995, Journal of Fluid Mechanics.

[40]  Steven H. Frankel,et al.  Entropy Stable Spectral Collocation Schemes for the Navier-Stokes Equations: Discontinuous Interfaces , 2014, SIAM J. Sci. Comput..

[41]  J. Rossmanith.  High-Order Discontinuous Galerkin Finite Element Methods with Globally Divergence-Free Constrained Transport for Ideal MHD , 2013, 1310.4251.

[42]  S. Rebay,et al.  A High-Order Accurate Discontinuous Finite Element Method for the Numerical Solution of the Compressible Navier-Stokes Equations , 1997 .

[43]  Eitan Tadmor,et al.  Entropy stability theory for difference approximations of nonlinear conservation laws and related time-dependent problems , 2003, Acta Numerica.

[44]  Andrew R. Winters,et al.  A novel high-order, entropy stable, 3D AMR MHD solver with guaranteed positive pressure , 2016, J. Comput. Phys..

[45]  D. Rosenberg,et al.  Ideal evolution of MHD turbulence when imposing Taylor-Green symmetries , 2017 .

[46]  Christoph Altmann,et al.  Explicit discontinuous Galerkin methods for magnetohydrodynamics , 2012 .

[47]  Chi-Wang Shu,et al.  Entropy stable high order discontinuous Galerkin methods with suitable quadrature rules for hyperbolic conservation laws , 2017, J. Comput. Phys..

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

[49]  P. Roe,et al.  A Solution-Adaptive Upwind Scheme for Ideal Magnetohydrodynamics , 1999 .

[50]  H. C. Yee.  Skew-Symmetric Splitting and Stability of High Order Central Schemes , 2017 .

[51]  David A. Kopriva,et al.  Implementing Spectral Methods for Partial Differential Equations , 2009 .

[52]  M. Mock,et al.  Systems of conservation laws of mixed type , 1980 .

[53]  David A. Kopriva,et al.  Metric Identities and the Discontinuous Spectral Element Method on Curvilinear Meshes , 2006, J. Sci. Comput..

[54]  Claus-Dieter Munz,et al.  A RUNGE-KUTTA BASED DISCONTINUOUS GALERKIN METHOD WITH TIME ACCURATE LOCAL TIME STEPPING , 2011 .

[55]  Yong Liu,et al.  Entropy stable high order discontinuous Galerkin methods for ideal compressible MHD on structured meshes , 2018, J. Comput. Phys..

[56]  P. Mininni,et al.  Lack of universality in decaying magnetohydrodynamic turbulence. , 2009, Physical review. E, Statistical, nonlinear, and soft matter physics.

[57]  Antony Jameson,et al.  Formulation of Kinetic Energy Preserving Conservative Schemes for Gas Dynamics and Direct Numerical Simulation of One-Dimensional Viscous Compressible Flow in a Shock Tube Using Entropy and Kinetic Energy Preserving Schemes , 2008, J. Sci. Comput..

[58]  Gregor Gassner,et al.  Ideal GLM-MHD: About the entropy consistent nine-wave magnetic field divergence diminishing ideal magnetohydrodynamics equations , 2017, J. Comput. Phys..

[59]  Gregor Gassner,et al.  Split form nodal discontinuous Galerkin schemes with summation-by-parts property for the compressible Euler equations , 2016, J. Comput. Phys..

[60]  Dominik Derigs,et al.  A novel averaging technique for discrete entropy-stable dissipation operators for ideal MHD , 2016, J. Comput. Phys..

[61]  Eitan Tadmor,et al.  Well-balanced and energy stable schemes for the shallow water equations with discontinuous topography , 2011, J. Comput. Phys..

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

[63]  Michael Dumbser,et al.  High-Order Unstructured One-Step PNPM Schemes for the Viscous and Resistive MHD Equations , 2009 .

[64]  Florian Hindenlang,et al.  A Parallel 3D Discontinuous Galerkin Framework for Nonlinear Resistive MHD Simulations in Tokamak and Stellarator Geometries , 2016 .

[65]  Gregor Gassner,et al.  On the use of kinetic energy preserving DG-schemes for large eddy simulation , 2017, J. Comput. Phys..