Search

In this paper, continuous research is undertaken to explore the underlying mechanism of numerical shock instabilities of Godunov-type schemes for strong shocks. By conducting dissipation analysis of Godunov-type schemes and a sequence of numerical experiments, we are able to clarify that the instability may be attributed to insufficient entropy production inside the numerical shock structure. As a result, a general entropy-control technique for improving the robustness of various Godunov-type schemes at strong shocks is developed. It plays a part in guaranteeing that enough entropy is produced inside the numerical shock structure. Furthermore, such a modified approach does not introduce any additional numerical dissipation on linear degenerate waves to suppress the shock instability. Numerical results that are obtained for various test cases indicate that the proposed methods have a good performance in terms of accuracy and robustness. [ABSTRACT FROM AUTHOR]

A hybrid numerical flux scheme is proposed by adapting the carbuncle-free modified Harten-Lax-van Leer contact (HLLCM) scheme to smoothly revert to the Harten-Lax-van Leer contact (HLLC) scheme in regions of shear. This hybrid scheme, referred to as the HLLCT scheme, employs a novel, velocity-based shear sensor. In contrast to the non-local pressure-based shock sensors often used in carbuncle cures, the proposed shear sensor can be computed in a localized manner meaning that the HLLCT scheme can be easily introduced into existing codes without having to implement additional data structures. Through numerical experiments, it is shown that the HLLCT scheme is able to resolve shear layers accurately without succumbing to the shock instability. [ABSTRACT FROM AUTHOR]

In this study, we investigate the so called carbuncle phenomenon by means of numerical experiments and heuristic considerations. We identify two main sources for the carbuncle: instability of the 1d shock position and low numerical viscosity on shear waves. We also describe how higher order stabilizes the 1d shock position and, thus, reduces the carbuncle. [ABSTRACT FROM AUTHOR]

It is well known that low diffusion Riemann solvers with minimal smearing on contact and shear waves are vulnerable to shock instability problems, including the carbuncle phenomenon. In the present study, we concentrate on exploring where the instability grows out and how the dissipation inherent in Riemann solvers affects the unstable behaviors. With the help of numerical experiments and a linearized analysis method, it has been found that the shock instability is strongly related to the unstable modes of intermediate states inside the shock structure. The consistency of mass flux across the normal shock is needed for a Riemann solver to capture strong shocks stably. The famous carbuncle phenomenon is interpreted as the consequence of the inconsistency of mass flux across the normal shock for a low diffusion Riemann solver. Based on the results of numerical experiments and the linearized analysis, a robust Godunov-type scheme with a simple cure for the shock instability is suggested. With only the dissipation corresponding to shear waves introduced in the vicinity of strong shocks, the instability problem is circumvented. Numerical results of several carefully chosen strong shock wave problems are investigated to demonstrate the robustness of the proposed scheme. [ABSTRACT FROM AUTHOR]

In this paper, we propose new Euler flux functions for use in a finite-volume Euler/Navier–Stokes code, which are very simple, carbuncle-free, yet have an excellent boundary-layer-resolving capability, by combining two different Riemann solvers into one based on a rotated Riemann solver approach. We show that very economical Euler flux functions can be devised by combining the Roe solver (a full-wave solver) and the Rusanov/HLL solver (a fewer-wave solver), based on a rotated Riemann solver approach: a fewer-wave solver automatically applied in the direction normal to shocks to suppress carbuncles and a full-wave solver applied, again automatically, across shear layers to avoid an excessive amount of dissipation. The resulting flux functions can be implemented in a very simple and economical manner, in the form of the Roe solver with modified wave speeds, so that converting an existing Roe flux code into the new fluxes is an extremely simple task. They require only 7–14% extra CPU time and no problem-dependent tuning parameters. These new rotated fluxes are not only robust for shock-capturing, but also accurate for resolving shear layers. This is demonstrated by an extensive series of numerical experiments with standard finite-volume Euler and Navier–Stokes codes, including various shock instability problems and also an unstructured grid case. [Copyright &y& Elsevier]