Evaluation of Euler Fluxes for Hypersonic Heating Computations

Peer Reviewedhttp://deepblue.lib.umich.edu/bitstream/2027.42/83575/1/AIAA-41605-439.pd

A hybrid numerical flux for supersonic flows with application to rocket nozzles

The numerical simulation of shock waves in supersonic flows is challenging because of several instabilities which can affect the solution. Among them, the carbuncle phenomenon can introduce nonphysical perturbations in captured shock waves. In the present work, a hybrid numerical flux is proposed for the evaluation of the convective fluxes that avoids carbuncle and keeps high-accuracy on shocks and boundary layers. In particular, the proposed flux is a combination between an upwind approximate Riemann problem solver and the Local Lax-Friedrichs scheme. A simple strategy to mix the two fluxes is proposed and tested in the framework of a discontinuous Galerkin discretisation. The approach is investigated on the subsonic flow in a channel, on the supersonic flow around a cylinder, on the supersonic flow on a flat plate and on the flow in a overexpanded rocket nozzle

Low Mach number effect in simulation of high Mach number flow

In this note, we relate the two well-known difficulties of Godunov schemes: the carbuncle phenomena in simulating high Mach number flow, and the inaccurate pressure profile in simulating low Mach number flow. We introduced two simple low-Mach-number modifications for the classical Roe flux to decrease the difference between the acoustic and advection contributions of the numerical dissipation. While the first modification increases the local numerical dissipation, the second decreases it. The numerical tests on the double-Mach reflection problem show that both modifications eliminate the kinked Mach stem suffered by the original flux. These results suggest that, other than insufficient numerical dissipation near the shock front, the carbuncle phenomena is strongly relevant to the non-comparable acoustic and advection contributions of the numerical dissipation produced by Godunov schemes due to the low Mach number effect.Comment: 9 pages, 1 figur

A carbuncle cure for the Harten-Lax-van Leer contact (HLLC) scheme using a novel velocity-based sensor

AbstractA 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.</jats:p

Advanced numerical methods for dissipative and non-dissipative relativistic hydrodynamics

High-energy physical phenomena such as astrophysical events and heavy-ion collisions contain a hydrodynamic aspect in which a branch of fluid dynamics called relativistic hydrodynamics (RHD) is required for its mathematical description. The resulting equations must be, more often than not, solved numerically for scientists to ascertain useful information regarding the fluid system in question. This thesis describes and presents a twodimensional computational fluid dynamics (CFD) solver for dissipative and non-dissipative relativistic hydrodynamics, i.e. in the presence and absence of physically resolved viscosity and heat conduction. The solver is based on a finite volume, Godunov-type, HighResolution Shock-Capturing (HRSC) framework, containing a plethora of numerical implementations such as high-order Weighted-Essentially Non-Oscillatory (WENO) spatial reconstruction, approximate Riemann solvers and a third-order Total Variation Diminishing (TVD) Runge–Kutta method. The base numerical solver for the solution of non-dissipative RHD is extensively tested using a series of one-dimensional test cases, namely, a smooth flow problem and shock-tube configurations as well as the two-dimensional vortex sheet and Riemann problem test cases. For the case of non-dissipative relativistic hydrodynamics the relativistic CFD solver is found to perform well in terms of the orders of accuracy achieved and its ability to resolve shock wave patterns. Numerical pathologies have been identified when the relativistic HLLC Riemann solver is used in multi-dimensions for problems exhibiting strong shock waves. This is attributed to the so-called Carbuncle problem which is shown to occur because of pressure differencing within the process of restoring the missing contact discontinuity of its predecessor, the HLL Riemann solver. To avoid this numerical pathology and improve the robustness of numerical solutions that make use of the HLLC Riemann solver, the development of a rotated-hybrid Riemann solver arising from the hybridisation of the HLL and HLLC (or Rusanov and HLLC) approximate Riemann solvers is presented. A standalone application of the HLLC Riemann solver can produce spurious numerical artefacts when it is employed in conjunction with Godunov-type high-order methods in the presence of discontinuities. It has been found that a rotated-hybrid Riemann solver with the proposed HLL/HLLC (Rusanov/HLLC) scheme could overcome the difficulty of the spurious numerical artefacts and presents a robust solution for the Carbuncle problem. The proposed rotated-hybrid Riemann solver provides sufficient numerical dissipation to capture the behaviour of strong shock waves for relativistic hydrodynamics. Therefore, focus is placed on two benchmark test cases (odd-even decoupling and double-Mach reflection problems) and the investigation of two astrophysical phenomena, the relativistic Richtmyer– Meshkov instability and the propagation of a relativistic jet. In all presented test cases, the Carbuncle problem is shown to be eliminated by employing the proposed rotated-hybrid Riemann solver. This strategy is problem-independent, straightforward to implement and provides a consistent robust numerical solution when combined with Godunov-type highorder schemes for relativistic hydrodynamics...[cont.]Aerospac