The suitability of Roe-type upwind schemes for the computation of steady quasi-1D flows of non-ideal fluids using explicit integration in the pseudo time is discussed. Based on the particular Roe linearization and entropy fix technique applied, several numerical difficulties can arise in modeling shock waves that are sonic either on the upstream or downstream side of the shock. These so-called sonic shocks typically occur away from stationary points of the cross-sectional area distribution, namely, where the geometrical source term does not vanish. The problem of selecting suitable formulations for the accurate simulation of non-classical quasi-1D steady flows is therefore addressed. Numerical experiments indicate a limited influence of the chosen Roe linearization technique, provided the so-called Property U, allowing the exact representation of steady shocks, is satisfied. Nevertheless, application of standard entropy fixes may either predict an incorrect steady-state transonic expansion neighbouring the sonic shock or even fail to attain a discrete steady state. In the latter case, lack of convergence is due to numerical unbalancing of the flux difference and source term integral over the transonic expansion which occurs in the close proximity of sonic shocks approaching their steady-state position. A simple modification to the synchronous splitting technique of [28] is proposed, which is able to produce the desired steady-state balance and allows substantial improvement in the resolution of sonic shocks.

On the numerical simulation of non-classical quasi-1D steady nozzle flows: Capturing sonic shocks

VIMERCATI, DAVIDE;GUARDONE, ALBERTO MATTEO ATTILIO
2018-01-01

Abstract

The suitability of Roe-type upwind schemes for the computation of steady quasi-1D flows of non-ideal fluids using explicit integration in the pseudo time is discussed. Based on the particular Roe linearization and entropy fix technique applied, several numerical difficulties can arise in modeling shock waves that are sonic either on the upstream or downstream side of the shock. These so-called sonic shocks typically occur away from stationary points of the cross-sectional area distribution, namely, where the geometrical source term does not vanish. The problem of selecting suitable formulations for the accurate simulation of non-classical quasi-1D steady flows is therefore addressed. Numerical experiments indicate a limited influence of the chosen Roe linearization technique, provided the so-called Property U, allowing the exact representation of steady shocks, is satisfied. Nevertheless, application of standard entropy fixes may either predict an incorrect steady-state transonic expansion neighbouring the sonic shock or even fail to attain a discrete steady state. In the latter case, lack of convergence is due to numerical unbalancing of the flux difference and source term integral over the transonic expansion which occurs in the close proximity of sonic shocks approaching their steady-state position. A simple modification to the synchronous splitting technique of [28] is proposed, which is able to produce the desired steady-state balance and allows substantial improvement in the resolution of sonic shocks.
2018
Non-classical gasdynamics; Roe scheme; Source terms; Steady-state flows; Computational Mathematics; Applied Mathematics
File in questo prodotto:
File Dimensione Formato  
VIMED01-18.pdf

Accesso riservato

Descrizione: Paper
: Publisher’s version
Dimensione 2.67 MB
Formato Adobe PDF
2.67 MB Adobe PDF   Visualizza/Apri
VIMED_OA_01-18.pdf

Open Access dal 05/08/2019

Descrizione: Paper open access
: Post-Print (DRAFT o Author’s Accepted Manuscript-AAM)
Dimensione 2.9 MB
Formato Adobe PDF
2.9 MB Adobe PDF Visualizza/Apri

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/11311/1031813
Citazioni
  • ???jsp.display-item.citation.pmc??? ND
  • Scopus 4
  • ???jsp.display-item.citation.isi??? 2
social impact