kth.sePublications KTH
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
A bound preserving cut discontinuous Galerkin method for one dimensional hyperbolic conservation laws
Nanjing Univ Aeronaut & Astronaut, Sch Math, Nanjing 211106, Peoples R China..
Uppsala Univ, Dept Informat Technol, Div Sci Comp, SE-75105 Uppsala, Sweden..
KTH, School of Engineering Sciences (SCI), Mathematics (Dept.), Numerical Analysis, Optimization and Systems Theory.ORCID iD: 0000-0002-4911-467X
2024 (English)In: ESAIM: Mathematical Modelling and Numerical Analysis (ESAIM: M2AN), ISSN 2822-7840, E-ISSN 2804-7214, Vol. 58, no 5, p. 1651-1680Article in journal (Refereed) Published
Abstract [en]

In this paper we present a family of high order cut finite element methods with bound preserving properties for hyperbolic conservation laws in one space dimension. The methods are based on the discontinuous Galerkin framework and use a regular background mesh, where interior boundaries are allowed to cut through the mesh arbitrarily. Our methods include ghost penalty stabilization to handle small cut elements and a new reconstruction of the approximation on macro-elements, which are local patches consisting of cut and un-cut neighboring elements that are connected by stabilization. We show that the reconstructed solution retains conservation and order of convergence. Our lowest order scheme results in a piecewise constant solution that satisfies a maximum principle for scalar hyperbolic conservation laws. When the lowest order scheme is applied to the Euler equations, the scheme is positivity preserving in the sense that positivity of pressure and density are retained. For the high order schemes, suitable bound preserving limiters are applied to the reconstructed solution on macro-elements. In the scalar case, a maximum principle limiter is applied, which ensures that the limited approximation satisfies the maximum principle. Correspondingly, we use a positivity preserving limiter for the Euler equations, and show that our scheme is positivity preserving. In the presence of shocks additional limiting is needed to avoid oscillations, hence we apply a standard TVB limiter to the reconstructed solution. The time step restrictions are of the same order as for the corresponding discontinuous Galerkin methods on the background mesh. Numerical computations illustrate accuracy, bound preservation, and shock capturing capabilities of the proposed schemes.

Place, publisher, year, edition, pages
EDP Sciences , 2024. Vol. 58, no 5, p. 1651-1680
Keywords [en]
Hyperbolic conservation laws, discontinuous Galerkin method, cut elements, satisfying maximum principle, positivity preserving
National Category
Computational Mathematics
Identifiers
URN: urn:nbn:se:kth:diva-354606DOI: 10.1051/m2an/2024042ISI: 001317859200001Scopus ID: 2-s2.0-85205574724OAI: oai:DiVA.org:kth-354606DiVA, id: diva2:1904089
Note

QC 20241008

Available from: 2024-10-08 Created: 2024-10-08 Last updated: 2025-08-28Bibliographically approved

Open Access in DiVA

No full text in DiVA

Other links

Publisher's full textScopus

Authority records

Zahedi, Sara

Search in DiVA

By author/editor
Zahedi, Sara
By organisation
Numerical Analysis, Optimization and Systems Theory
In the same journal
ESAIM: Mathematical Modelling and Numerical Analysis (ESAIM: M2AN)
Computational Mathematics

Search outside of DiVA

GoogleGoogle Scholar

doi
urn-nbn

Altmetric score

doi
urn-nbn
Total: 62 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