Enhanced HLLEM and HLL-CPS schemes for all Mach number flows based using anti-diffusion coefficients

A. Gogoi,J.C. Mandal

2024-11-14

Abstract:This paper compares the HLLEM and HLL-CPS schemes for Euler equations and proposes improvements for all Mach number flows. Enhancements to the HLLEM scheme involve adding anti-diffusion terms in the face normal direction and modifying anti-diffusion coefficients for linearly degenerate waves near shocks. The HLL-CPS scheme is improved by adjusting anti-diffusion coefficients for the face normal direction and linearly degenerate waves. Matrix stability, linear perturbation, and asymptotic analyses demonstrate the robustness of the proposed schemes and their ability to capture low Mach flow features. Numerical tests confirm that the schemes are free from shock instabilities at high speeds and accurately resolve low Mach number flow features.

Numerical Analysis,Computational Physics

What problem does this paper attempt to address?

This paper attempts to solve the following two main problems:

1. **Numerical shock instability in high - Mach - number flows**:
- Existing approximate Riemann solvers (such as HLLEM and HLL - CPS) are prone to numerical shock instability when dealing with high - speed flows. These instabilities include the carbuncle phenomenon, odd - even decoupling, kinked Mach stem, etc.
- For example, although the HLLEM scheme can accurately capture contact discontinuities and shear waves, it is prone to produce numerical shock instability at high Mach numbers; while the HLL - CPS scheme can avoid the carbuncle phenomenon in some cases, it still shows slight instability in supersonic corner tests.

2. **Insufficient accuracy in low - Mach - number flows**:
- Existing approximate Riemann solvers have difficulty in accurately capturing flow - field characteristics when dealing with low - Mach - number flows, especially on structured grids. These problems include the checkerboard - type decoupling of the pressure and velocity fields.
- For example, the HLLEM and HLL - CPS schemes perform poorly in low - Mach - number flows and cannot accurately resolve boundary layers and shear layers.

To solve the above problems, this paper proposes an improved HLLEM and HLL - CPS scheme, which enhances their stability and accuracy in the full Mach - number range by introducing an anti - diffusion coefficient. Specifically:

- For the HLLEM scheme, the improvements include adding an anti - diffusion term in the face - normal direction and modifying the anti - diffusion coefficient near the linearly degenerate wave.
- For the HLL - CPS scheme, the improvements include adjusting the anti - diffusion coefficients in the face - normal direction and for the linearly degenerate wave.

Through matrix stability analysis, linear perturbation analysis and asymptotic analysis, the stability of the proposed improved scheme at high Mach numbers and the accuracy at low Mach numbers are verified. Numerical test results also show that the improved scheme has no shock instability in high - speed flows and can accurately resolve low - Mach - number flow characteristics.

### Key formulas

1. **Conservative form of the Euler equations**:
\[
\frac{\partial \mathbf{U}}{\partial t}+\frac{\partial \mathbf{F}(\mathbf{U})}{\partial x}+\frac{\partial \mathbf{G}(\mathbf{U})}{\partial y} = 0
\]
where:
\[
\mathbf{U}=\begin{bmatrix}
\rho\\
\rho u\\
\rho v\\
\rho E
\end{bmatrix}, \quad
\mathbf{F}(\mathbf{U})=\begin{bmatrix}
\rho u\\
\rho u^{2}+p\\
\rho uv\\
(\rho E + p)u
\end{bmatrix}, \quad
\mathbf{G}(\mathbf{U})=\begin{bmatrix}
\rho v\\
\rho uv\\
\rho v^{2}+p\\
(\rho E + p)v
\end{bmatrix}
\]

2. **Equation of state**:
\[
p = (\gamma - 1)\left(\rho E-\frac{1}{2}\rho(u^{2}+v^{2})\right)
\]
where \(\gamma\) is the specific heat ratio, which is assumed to be 1.4 in this paper.

3. **HLLEM flux formula**:
\[
F(U_R,