Abstract:In this paper, a positivity-preserving fifth-order finite volume compact-WENO scheme is proposed for solving compressible Euler equations. As we know conservative compact finite volume schemes have high resolution properties while WENO (Weighted Essentially Non-Oscillatory) schemes are essentially non-oscillatory near flow discontinuities. We extend the main idea of WENO schemes to some classical compact finite volume schemes [32], where lower order compact stencils are combined with WENO nonlinear weights to get a higher order finite volume compact-WENO scheme. The newly developed positivity-preserving limiter [46,44] is used to preserve positive density and internal energy for compressible Euler equations of fluid dynamics. The HLLC (Harten, Lax, and van Leer with Contact) approximate Riemann solver [39,2] is used to get the numerical flux at the cell interfaces. Numerical tests are presented to demonstrate the high-order accuracy, positivity-preserving, high-resolution and robustness of the proposed scheme.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is how to design a high - order finite - volume compact WENO (FVCW) scheme when solving the compressible Euler equations, so as to maintain positive density and internal energy, while having high resolution and robustness. Specifically, the paper focuses on how to combine the advantages of the compact scheme and the WENO scheme when dealing with discontinuous problems in fluid dynamics, and develop a method that can effectively avoid non - physical oscillations and maintain the stability of the numerical solution. The paper mentions that although the traditional linear compact scheme has high precision and good resolution, it will produce unacceptable non - physical oscillations when dealing with flows containing discontinuities. And although the WENO scheme can maintain the essentially non - oscillatory property near discontinuities, it usually has poor spectral resolution and excessive numerical dissipation. Therefore, this paper proposes a new high - order finite - volume compact WENO scheme, aiming to overcome these shortcomings and ensure the positive density and internal energy of the numerical solution, especially in extreme cases, such as vacuum or near - vacuum solutions.
### Main contributions
1. **High - order finite - volume compact WENO scheme**: Combining the high - resolution characteristics of the compact scheme and the essentially non - oscillatory characteristics of the WENO scheme, a fifth - order finite - volume compact WENO scheme is designed.
2. **Maintenance of positive density and internal energy**: A new positivity - preserving limiter is introduced to ensure that the density and internal energy are always positive during the numerical calculation process.
3. **HLLC approximate Riemann solver**: The HLLC approximate Riemann solver is used to calculate the numerical flux at the cell interface to reduce numerical dissipation and improve robustness.
4. **Verification by numerical experiments**: The high - order accuracy, high - resolution, positivity - preserving characteristics and robustness of the proposed FVCW scheme are verified through a series of numerical experiments.
### Mathematical model
The compressible Euler equations considered in the paper can be expressed as:
\[ \frac{\partial U}{\partial t}+\frac{\partial F(U)}{\partial x} = 0, \]
where \( U \) and \( F(U) \) are the conserved variables and the flux vector respectively, and are defined as:
\[ U=\begin{pmatrix} \rho \\ \rho u \\ E \end{pmatrix}, \quad F(U)=\begin{pmatrix} \rho u \\ \rho u^{2}+p \\ u(E + p) \end{pmatrix}, \]
where:
- \( \rho \) is the density,
- \( u \) is the particle velocity,
- \( p \) is the pressure,
- \( E \) is the total energy per unit volume,
- \( e \) is the specific internal energy, defined as \( e=\frac{p}{(\gamma - 1)\rho} \),
- \( \gamma \) is the specific heat ratio (for an ideal gas, \( \gamma = 1.4 \)).
The speed of sound \( a \) is defined as:
\[ a=\sqrt{\frac{\gamma p}{\rho}}. \]
### Numerical methods
1. **Finite - volume scheme**: By integrating the conservation equation on each cell, the finite - volume scheme is obtained.
2. **Compact WENO reconstruction**: At each cell interface, the compact WENO reconstruction is used to obtain the numerical flux with high - order accuracy.
3. **Positivity - preserving limiter**: The newly developed positivity - preserving rescaling limiter is used to ensure that the density and internal energy always remain positive.
4. **HLLC approximate Riemann solver**: It is used to calculate the numerical flux at the cell interface to reduce numerical dissipation and improve robustness.
### Conclusion
Through numerical experiments, the paper verifies that the proposed FVCW scheme has high - order accuracy, high - resolution, positivity - preserving characteristics and robustness when dealing with the compressible Euler equations. These characteristics make the scheme perform well in dealing with complex fluid - dynamic problems.