L²Roe: a low dissipation version of Roe's approximate Riemann solver for low Mach numbers

Abstract

Summary A modification of the Roe scheme called L2Roe for low dissipation low Mach Roe is presented. It reduces the dissipation of kinetic energy at the highest resolved wave numbers in a low Mach number test case of decaying isotropic turbulence. This is achieved by scaling the jumps in all discrete velocity components within the numerical flux function. An asymptotic analysis is used to show the correct pressure scaling at low Mach numbers and to identify the reduced numerical dissipation in that regime. Furthermore, the analysis allows a comparison with two other schemes that employ different scaling of discrete velocity jumps, namely, LMRoe and a method of Thornber et al. To this end, we present for the first time an asymptotic analysis of the last method. Numerical tests on cases ranging from low Mach number (M∞=0.001) to hypersonic (M∞=5) viscous flows are used to illustrate the differences between the methods and to show the correct behavior of L2Roe. No conflict is observed between the reduced numerical dissipation and the accuracy or stability of the scheme in any of the investigated test cases. Copyright © 2015 John Wiley & Sons, Ltd.

Spatial eigensolution analysis of discontinuous Galerkin schemes with practical insights for under-resolved computations and implicit LES

“…Also, some corrections of Riemann solvers aimed at low Mach number applications can effectively reduce them to nearly central fluxes around low-speed flow regions, cf. [47] . The dispersion curves of the physical modes shown in Fig.…”

An accurate and robust HLLC‐type Riemann solver for the compressible Euler system at various Mach numbers

“…Pressure fluctuations of order M 2 : Similar to the study in the work of Guillard and Viozat, 23 Equations (36), (37) and (38) or (41) imply the uniformity of the leading order pressure in space…”

“…where p * * is defined in Equation (31), p * is defined in Equation (10), f is defined in Equations (17) and (18). As argued in the work of Oßwald et al, 37 the all Mach correction term defined in Equation (31) helps to improve the behavior of the scheme at low Mach numbers, but it will introduce small disturbances in the vicinity of shocks. Thus, in the above formulation (34), we apply the pressure weight function in Equations (17) and (18) to turn off the all Mach correction at shocks.…”

“…[27][28][29] Unlike the Roe-Turkel scheme proposed by Guillard and Viozat 23 which modified both the eigenvalues and eigenvectors of the numerical dissipation, the all-speed Roe scheme only modified the nonlinear eigenvalues. [35][36][37] Although the low Mach number corrections for Godunov-type schemes are diverse, little attention has been devoted to the low Mach number problem of the HLLC scheme. Dellacherie et al 33,34 defined and analyzed the low Mach number problem through a linear analysis of a perturbed linear wave equation.…”

“…Furthermore, the theoretical framework of Dellacherie et al also justified other existing all Mach schemes, for example, the mentioned preconditioning methods of Guillard and Viozat, 23 Rieper's LMRoe scheme, 30 and other similar approaches. [35][36][37] Although the low Mach number corrections for Godunov-type schemes are diverse, little attention has been devoted to the low Mach number problem of the HLLC scheme. It was expected in the work of Rieper 30 that the low Mach number fix for Roe's method can be applied to other Godunov-type schemes including the HLLC approximate Riemann solver, but the extension was not straightforward due to the nonlinear formulations of its flux function.…”