We present an implicit-explicit well-balanced finite volume scheme for the Euler equations with a gravitational source term which is able to deal also with low Mach flows. To visualize the different scales we use the non-dimensionalized equations on which we apply a pressure splitting and a Suliciu relaxation. On the resulting model, we apply a splitting of the flux into a linear implicit and an non-linear explicit part that leads to a scale independent time-step. The explicit step consists of a Godunov type method based on an approximative Riemann solver where the source term is included in the flux formulation. We develop the method for a first order scheme and give an extension to second order. Both schemes are designed to be well-balanced, preserve the positivity of density and internal energy and have a scale independent diffusion. We give the low Mach limit equations for well-prepared data and show that the scheme is asymptotic preserving. These properties are numerically validated by various test cases.

An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity / Thomann, A.; Puppo, G.; Klingenberg, C.. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 420:(2020), p. 109723. [10.1016/j.jcp.2020.109723]

An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity

Puppo G.;
2020

Abstract

We present an implicit-explicit well-balanced finite volume scheme for the Euler equations with a gravitational source term which is able to deal also with low Mach flows. To visualize the different scales we use the non-dimensionalized equations on which we apply a pressure splitting and a Suliciu relaxation. On the resulting model, we apply a splitting of the flux into a linear implicit and an non-linear explicit part that leads to a scale independent time-step. The explicit step consists of a Godunov type method based on an approximative Riemann solver where the source term is included in the flux formulation. We develop the method for a first order scheme and give an extension to second order. Both schemes are designed to be well-balanced, preserve the positivity of density and internal energy and have a scale independent diffusion. We give the low Mach limit equations for well-prepared data and show that the scheme is asymptotic preserving. These properties are numerically validated by various test cases.
2020
Asymptotic preserving; Euler equations with gravity; IMEX scheme; positivity preserving; Suliciu relaxation; well-balanced
01 Pubblicazione su rivista::01a Articolo in rivista
An all speed second order well-balanced IMEX relaxation scheme for the Euler equations with gravity / Thomann, A.; Puppo, G.; Klingenberg, C.. - In: JOURNAL OF COMPUTATIONAL PHYSICS. - ISSN 0021-9991. - 420:(2020), p. 109723. [10.1016/j.jcp.2020.109723]
File allegati a questo prodotto
File Dimensione Formato  
Thomann_An-all-speed-second-order_2020.pdf

solo gestori archivio

Note: 10.1016/j.jcp.2020.109723
Tipologia: Versione editoriale (versione pubblicata con il layout dell'editore)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 1.09 MB
Formato Adobe PDF
1.09 MB Adobe PDF   Contatta l'autore
Thomann_preprint_An-all-speed-second-order_2020.pdf

accesso aperto

Tipologia: Documento in Pre-print (manoscritto inviato all'editore, precedente alla peer review)
Licenza: Tutti i diritti riservati (All rights reserved)
Dimensione 1.34 MB
Formato Adobe PDF
1.34 MB Adobe PDF

I documenti in IRIS sono protetti da copyright e tutti i diritti sono riservati, salvo diversa indicazione.

Utilizza questo identificativo per citare o creare un link a questo documento: https://hdl.handle.net/11573/1467287
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 24
  • ???jsp.display-item.citation.isi??? 20
social impact