N.H. Maruthi,S.V. Raghurama Rao
Abstract:An exact discontinuity capturing central solver developed recently, named MOVERS (Method of Optimal Viscosity for Enhanced Resolution of Shocks, J Computat Phys 2009;228:770-798), is analyzed and improved further to make it entropy stable. MOVERS, which is designed to capture steady shocks and contact discontinuities exactly by enforcing the Rankine-Hugoniot jump condition directly in the discretization process, is a low diffusive algorithm in a simple central discretization framework, free of complicated Riemann solvers and flux splittings. However, this algorithm needs an entropy fix to avoid nonsmoothness in the expansion regions. The entropy conservation equation is used as a guideline to introduce an optimal numerical diffusion in the smooth regions and a limiter based switchover is introduced for numerical diffusion based on jump conditions at the large gradients. The resulting new scheme is entropy stable, accurate and captures steady discontinuities exactly while avoiding an entropy fix.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop a numerical method that can accurately capture steady - state discontinuities (such as shock waves and contact discontinuities) while maintaining entropy stability in smooth regions. Specifically, the researchers hope to improve the existing MOVERS algorithm so that it can avoid non - smoothness problems without using entropy correction and reduce numerical dissipation in smooth regions.
### Main problems
1. **Accurately capture steady - state discontinuities**: Many existing numerical methods (such as Riemann solvers) can capture discontinuities, but usually need to introduce additional numerical dissipation to ensure stability, which may lead to less accurate results.
2. **Entropy stability**: When dealing with expansion regions, some algorithms may produce non - physical phenomena (such as expansion shock waves). To prevent this, entropy correction is usually introduced, but this will affect the algorithm's ability to accurately capture steady - state discontinuities.
3. **Reduce numerical dissipation in smooth regions**: Although the existing MOVERS - L algorithm can capture steady - state discontinuities without using entropy correction, it introduces too much numerical dissipation in smooth regions, affecting the calculation accuracy.
### Solutions
The author proposes a new hybrid algorithm - MOVERS - LE (L stands for limiter, E stands for entropy stability), aiming to combine the advantages of MOVERS and LLF methods. The specific methods are as follows:
- **In discontinuous regions**: Use numerical dissipation based on the Rankine - Hugoniot condition (from MOVERS) to ensure accurate capture of high - gradient regions.
- **In smooth regions**: Use numerical dissipation based on the entropy conservation equation (from the LLF method) to ensure entropy stability and reduce numerical dissipation.
In this way, MOVERS - LE can maintain accurate capture of steady - state discontinuities, provide less numerical dissipation in smooth regions, and ensure the entropy stability of the overall scheme.
### Mathematical formulas
- **Rankine - Hugoniot condition**:
\[
\Delta F = s\Delta U
\]
where \(\Delta F = F_R - F_L\) and \(\Delta U = U_R - U_L\), and \(s\) is the velocity of the discontinuity surface.
- **Numerical dissipation coefficient**:
\[
\alpha_{j + 1/2}=|s_{j + 1/2}|
\]
where \(s_{j + 1/2}\) is the velocity of the discontinuity surface.
- **Entropy conservation equation**:
\[
\frac{\partial(\rho S)}{\partial t}+\frac{\partial(\rho u S)}{\partial x} = 0
\]
where \(S\) is entropy, defined as:
\[
S = C_V\log\left(\frac{p}{\rho^\gamma}\right)+\text{constant}
\]
- **Numerical dissipation coefficient (entropy stability)**:
\[
a_{ES}=\frac{\partial(\rho u S)}{\partial(\rho S)}\approx u
\]
\[
a_{LLF - ES}=\max[|u_L|, |u_R|]
\]
Through these improvements, MOVERS - LE can perform well in various test cases, including steady - state shock waves, contact discontinuities, Sod - modified shock - tube problems, and strong discontinuities.