Circumventing shock anomalies: a hybrid Riemann solver with high-order schemes

Ioannis Kokkinakis, Dimitris Drikakis

Research output: Contribution to conferenceAbstract

Abstract

Over the past 4 decades, there has been several publications on high-resolution methods and their applications to shock waves and turbulence. Despite the above efforts, the numerical behaviour of such schemes in the presence of very strong shockwaves, e.g., for Mach numbers greater than 10, is rather limited. In particular, the effect of grid misalignment to strong normal shockwaves is known to lead to unphysical phenomena that adversely affect the quality and even validity of the resulting flow (see [1, 2, 3] and references therein). The problem becomes particularly severe for Finite Volume (FV) methods when a contact resolving Riemann solver is used.

Riemann solvers that have a single formulation and claim to be carbuncle-free typically introduce a shear viscosity component in some form or another. This stabilizes the simulation, but it also results in inaccuracies in the prediction of steady shear waves, thus leading to a poor resolution of boundary layers. Hybrid Riemann solvers provide an alternative because the numerical dissipation can be better controlled within specific flow regions. In this study, we investigate a hybrid HLLC – HLL Riemann solver in conjunction with 2nd and higher-order schemes.

Several flow examples are presented, including numerical evidence showing that the hybrid Riemann solver does not generate artificial flow structures in the presence of very strong shocks. The effects of the spatial order of reconstruction on the accuracy and stability of the simulations is also discussed.

Conference

Conference6th European Conference on Computational Mechanics and 7th European Conference on Computational Fluid Dynamics 2018
Abbreviated titleECCM - ECFD 2018
CountryUnited Kingdom
CityGlasgow
Period11/06/1815/06/18

Fingerprint

Shear viscosity
Shear waves
Finite volume method
Flow structure
Shock waves
Mach number
Boundary layers
Turbulence

Keywords

  • high-order
  • finite volume
  • shockwave
  • Riemann solver
  • shock capturing

Cite this

Kokkinakis, I., & Drikakis, D. (2018). Circumventing shock anomalies: a hybrid Riemann solver with high-order schemes. 1. Abstract from 6th European Conference on Computational Mechanics and 7th European Conference on Computational Fluid Dynamics 2018, Glasgow, United Kingdom.
Kokkinakis, Ioannis ; Drikakis, Dimitris. / Circumventing shock anomalies : a hybrid Riemann solver with high-order schemes. Abstract from 6th European Conference on Computational Mechanics and 7th European Conference on Computational Fluid Dynamics 2018, Glasgow, United Kingdom.1 p.
@conference{260efab81d6c4f37bb8aaf93b362b4ca,
title = "Circumventing shock anomalies: a hybrid Riemann solver with high-order schemes",
abstract = "Over the past 4 decades, there has been several publications on high-resolution methods and their applications to shock waves and turbulence. Despite the above efforts, the numerical behaviour of such schemes in the presence of very strong shockwaves, e.g., for Mach numbers greater than 10, is rather limited. In particular, the effect of grid misalignment to strong normal shockwaves is known to lead to unphysical phenomena that adversely affect the quality and even validity of the resulting flow (see [1, 2, 3] and references therein). The problem becomes particularly severe for Finite Volume (FV) methods when a contact resolving Riemann solver is used. Riemann solvers that have a single formulation and claim to be carbuncle-free typically introduce a shear viscosity component in some form or another. This stabilizes the simulation, but it also results in inaccuracies in the prediction of steady shear waves, thus leading to a poor resolution of boundary layers. Hybrid Riemann solvers provide an alternative because the numerical dissipation can be better controlled within specific flow regions. In this study, we investigate a hybrid HLLC – HLL Riemann solver in conjunction with 2nd and higher-order schemes. Several flow examples are presented, including numerical evidence showing that the hybrid Riemann solver does not generate artificial flow structures in the presence of very strong shocks. The effects of the spatial order of reconstruction on the accuracy and stability of the simulations is also discussed.",
keywords = "high-order, finite volume, shockwave, Riemann solver, shock capturing",
author = "Ioannis Kokkinakis and Dimitris Drikakis",
year = "2018",
month = "6",
day = "13",
language = "English",
pages = "1",
note = "6th European Conference on Computational Mechanics and 7th European Conference on Computational Fluid Dynamics 2018, ECCM - ECFD 2018 ; Conference date: 11-06-2018 Through 15-06-2018",

}

Kokkinakis, I & Drikakis, D 2018, 'Circumventing shock anomalies: a hybrid Riemann solver with high-order schemes' 6th European Conference on Computational Mechanics and 7th European Conference on Computational Fluid Dynamics 2018, Glasgow, United Kingdom, 11/06/18 - 15/06/18, pp. 1.

Circumventing shock anomalies : a hybrid Riemann solver with high-order schemes. / Kokkinakis, Ioannis; Drikakis, Dimitris.

2018. 1 Abstract from 6th European Conference on Computational Mechanics and 7th European Conference on Computational Fluid Dynamics 2018, Glasgow, United Kingdom.

Research output: Contribution to conferenceAbstract

TY - CONF

T1 - Circumventing shock anomalies

T2 - a hybrid Riemann solver with high-order schemes

AU - Kokkinakis, Ioannis

AU - Drikakis, Dimitris

PY - 2018/6/13

Y1 - 2018/6/13

N2 - Over the past 4 decades, there has been several publications on high-resolution methods and their applications to shock waves and turbulence. Despite the above efforts, the numerical behaviour of such schemes in the presence of very strong shockwaves, e.g., for Mach numbers greater than 10, is rather limited. In particular, the effect of grid misalignment to strong normal shockwaves is known to lead to unphysical phenomena that adversely affect the quality and even validity of the resulting flow (see [1, 2, 3] and references therein). The problem becomes particularly severe for Finite Volume (FV) methods when a contact resolving Riemann solver is used. Riemann solvers that have a single formulation and claim to be carbuncle-free typically introduce a shear viscosity component in some form or another. This stabilizes the simulation, but it also results in inaccuracies in the prediction of steady shear waves, thus leading to a poor resolution of boundary layers. Hybrid Riemann solvers provide an alternative because the numerical dissipation can be better controlled within specific flow regions. In this study, we investigate a hybrid HLLC – HLL Riemann solver in conjunction with 2nd and higher-order schemes. Several flow examples are presented, including numerical evidence showing that the hybrid Riemann solver does not generate artificial flow structures in the presence of very strong shocks. The effects of the spatial order of reconstruction on the accuracy and stability of the simulations is also discussed.

AB - Over the past 4 decades, there has been several publications on high-resolution methods and their applications to shock waves and turbulence. Despite the above efforts, the numerical behaviour of such schemes in the presence of very strong shockwaves, e.g., for Mach numbers greater than 10, is rather limited. In particular, the effect of grid misalignment to strong normal shockwaves is known to lead to unphysical phenomena that adversely affect the quality and even validity of the resulting flow (see [1, 2, 3] and references therein). The problem becomes particularly severe for Finite Volume (FV) methods when a contact resolving Riemann solver is used. Riemann solvers that have a single formulation and claim to be carbuncle-free typically introduce a shear viscosity component in some form or another. This stabilizes the simulation, but it also results in inaccuracies in the prediction of steady shear waves, thus leading to a poor resolution of boundary layers. Hybrid Riemann solvers provide an alternative because the numerical dissipation can be better controlled within specific flow regions. In this study, we investigate a hybrid HLLC – HLL Riemann solver in conjunction with 2nd and higher-order schemes. Several flow examples are presented, including numerical evidence showing that the hybrid Riemann solver does not generate artificial flow structures in the presence of very strong shocks. The effects of the spatial order of reconstruction on the accuracy and stability of the simulations is also discussed.

KW - high-order

KW - finite volume

KW - shockwave

KW - Riemann solver

KW - shock capturing

M3 - Abstract

SP - 1

ER -

Kokkinakis I, Drikakis D. Circumventing shock anomalies: a hybrid Riemann solver with high-order schemes. 2018. Abstract from 6th European Conference on Computational Mechanics and 7th European Conference on Computational Fluid Dynamics 2018, Glasgow, United Kingdom.