A family of functions for mass and energy flux splitting of the Euler equations

• Published in: Journal of computational physics Vol. 228; no. 23; pp. 8908 - 8918
• Main Authors: Raga, A.C., Cantó, J.
• Format: Journal Article
• Language: English
• Published: Kidlington Elsevier Inc 10.12.2009; Elsevier
• Subjects: Algorithms; Computational techniques; Contact; Diffusion; Discontinuity; Exact sciences and technology; Flux; Flux vector splitting; Flux vector splitting schemes; Fluxes; Mach number; Mathematical methods in physics; Physics; Shock capturing technique; Flux vector splitting schemes; Shock capturing technique; Discontinuity; Algorithms; Momentum transfer; Mass function; Calculation; Euler equations; Diffusion; Energy transfer; Calculation methods; Mach number
• ISSN: 0021-9991; 1090-2716
• DOI: 10.1016/j.jcp.2009.09.006

Flux vector splitting algorithms for the Euler equations are based on dividing the mass, momentum and energy fluxes into a “forward directed flux” F + and a “backward directed flux” F - (with F - = 0 for Mach numbers M > 1 and F + = 0 for M < - 1 ). van Leer (1979, 1982) [4,5] proposed using polynomials of the Mach number for computing F + and F - in the subsonic regime, and derived the lowest order polynomials that satisfy a set of chosen criteria. In this paper, we explore the possibility of increasing the order of these polynomials, with the purpose of reducing the diffusion across slow moving contact discontinuities of the flux vector splitting algorithm. We find that a moderate reduction of the diffusion, resulting in sharper shocks and contact discontinuities, can indeed be obtained with the higher order polynomials for the split fluxes.