Change search
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf
An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part I: Theory and numerical verification
Mathematisches Institut, Universität zu Köln, Köln, Germany.
Mathematisches Institut, Universität zu Köln, Köln, Germany.ORCID iD: 0000-0002-5902-1522
Mathematisches Institut, Universität zu Köln, Köln, Germany.
I. Physikalisches Institut, Universität zu Köln, Köln, Germany.
Show others and affiliations
2018 (English)In: Journal of Computational Physics, ISSN 0021-9991, E-ISSN 1090-2716Article in journal (Refereed) In press
Abstract [en]

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 non-conservative 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.

Place, publisher, year, edition, pages
Elsevier, 2018.
Keywords [en]
resistive magnetohydrodynamics, entropy stability, discontinuous Galerkin spectral element method, hyperbolic divergence cleaning, curvilinear hexahedral mesh, summation-by- parts
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:liu:diva-156851DOI: 10.1016/j.jcp.2018.06.027Scopus ID: 2-s2.0-85049958774OAI: oai:DiVA.org:liu-156851DiVA, id: diva2:1315781
Funder
EU, European Research Council, 714487Available from: 2019-05-14 Created: 2019-05-14 Last updated: 2019-05-23Bibliographically approved

Open Access in DiVA

An entropy stable nodal discontinuous Galerkin method for the resistive MHD equations. Part I: Theory and numerical verification(2911 kB)20 downloads
File information
File name FULLTEXT01.pdfFile size 2911 kBChecksum SHA-512
1aac42a1d357b9ab11bc6714d899b9f4b60e03a43e4e3f33b913696a082808f48b47f081c7a121c957fa5fa7305dec9c7271cc7df409a63f30d9ffd98cc37fba
Type fulltextMimetype application/pdf

Other links

Publisher's full textScopus

Search in DiVA

By author/editor
Winters, Andrew Ross
In the same journal
Journal of Computational Physics
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar
Total: 20 downloads
The number of downloads is the sum of all downloads of full texts. It may include eg previous versions that are now no longer available

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 38 hits
CiteExportLink to record
Permanent link

Direct link
Cite
Citation style
  • apa
  • ieee
  • modern-language-association-8th-edition
  • vancouver
  • Other style
More styles
Language
  • de-DE
  • en-GB
  • en-US
  • fi-FI
  • nn-NO
  • nn-NB
  • sv-SE
  • Other locale
More languages
Output format
  • html
  • text
  • asciidoc
  • rtf