Universität zu Köln

Bohm, Marvin, Winters, Andrew R., Gassner, Gregor J., Derigs, Dominik, Hindenlang, Florian and Saur, Joachim (2020). An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part I: Theory and numerical verification. J. Comput. Phys., 422. SAN DIEGO: ACADEMIC PRESS INC ELSEVIER SCIENCE. ISSN 1090-2716

Full text not available from this repository.

Abstract

The first paper of this series presents a discretely entropy stable discontinuous Galerkin (DG) method for the resistive magnetohydrodynamics (MHD) equations on three-dimensional curvilinear unstructured hexahedral meshes. Compared to other fluid dynamics systems such as the shallow water equations or the compressible Navier-Stokes equations, the resistive MHD equations need special considerations because of the divergence-free constraint on the magnetic field. For instance, it is well known that for the symmetrization of the ideal MHD system as well as the continuous entropy analysis a non-conservative term proportional to the divergence of the magnetic field, typically referred to as the Powell term, must be included. As a consequence, the mimicry of the continuous entropy analysis in the discrete sense demands a suitable DG approximation of the nonconservative terms in addition to the ideal MHD terms. This paper focuses on the resistive MHD equations: Our first contribution is a proof that the resistive terms are symmetric and positive-definite when formulated in entropy space as gradients of the entropy variables, which enables us to show that the entropy inequality holds for the resistive MHD equations. This continuous analysis is the key for our DG discretization and guides the path for the construction of an approximation that discretely mimics the entropy inequality, typically termed entropy stability. Our second contribution is a detailed derivation and analysis of the discretization on three-dimensional curvilinear meshes. The discrete analysis relies on the summation-by-parts property, which is satisfied by the DG spectral element method (DGSEM) with Legendre-Gauss-Lobatto (LGL) nodes. Although the divergence-free constraint is included in the non-conservative terms, the resulting method has no particular treatment of the magnetic field divergence errors, which might pollute the solution quality. Our final contribution is the extension of the standard resistive MHD equations and our DG approximation with a divergence cleaning mechanism that is based on a generalized Lagrange multiplier (GLM). As a conclusion to the first part of this series, we provide detailed numerical validations of our DGSEM method that underline our theoretical derivations. In addition, we show a numerical example where the entropy stable DGSEM demonstrates increased robustness compared to the standard DGSEM. (C) 2018 Elsevier Inc. All rights reserved.

Item Type:	Journal Article
Creators:	Creators Email ORCID ORCID Put Code Bohm, Marvin UNSPECIFIED UNSPECIFIED UNSPECIFIED Winters, Andrew R. UNSPECIFIED UNSPECIFIED UNSPECIFIED Gassner, Gregor J. UNSPECIFIED UNSPECIFIED UNSPECIFIED Derigs, Dominik UNSPECIFIED UNSPECIFIED UNSPECIFIED Hindenlang, Florian UNSPECIFIED UNSPECIFIED UNSPECIFIED Saur, Joachim UNSPECIFIED UNSPECIFIED UNSPECIFIED
URN:	urn:nbn:de:hbz:38-310521
DOI:	10.1016/j.jcp.2018.06.027
Journal or Publication Title:	J. Comput. Phys.
Volume:	422
Date:	2020
Publisher:	ACADEMIC PRESS INC ELSEVIER SCIENCE
Place of Publication:	SAN DIEGO
ISSN:	1090-2716
Language:	English
Faculty:	Unspecified
Divisions:	Unspecified
Subjects:	no entry
Uncontrolled Keywords:	Keywords Language NONLINEAR CONSERVATION-LAWS; COMPRESSIBLE NAVIER-STOKES; SHALLOW-WATER EQUATIONS; FINITE-ELEMENT-METHOD; HIGH-ORDER; BOUNDARY-CONDITIONS; DIFFERENCE-SCHEMES; VOLUME SCHEME; IDEAL MHD; SYSTEMS Multiple languages Computer Science, Interdisciplinary Applications; Physics, Mathematical Multiple languages
URI:	http://kups.ub.uni-koeln.de/id/eprint/31052