M.C. Martí,P. Mulet,D.F. Yáñez,D. Zorío
Abstract:A set of arbitrarily high-order WENO schemes for reconstructions on nonuniform grids is presented. These non-linear interpolation methods use simple smoothness indicators with a linear cost with respect to the order, making them easy to implement and computationally efficient. The theoretical analysis to verify the accuracy and the essentially non-oscillatory properties are presented together with some numerical experiments involving algebraic problems in order to validate them. Also, these general schemes are applied for the solution of conservation laws and hyperbolic systems in the context of finite volume methods.
What problem does this paper attempt to address?
The main problem that this paper attempts to solve is to develop an efficient WENO (Weighted Essentially Non - Oscillatory) scheme for dealing with the numerical reconstruction problem on non - uniform grids. Specifically, the authors aim to:
1. **Improve computational efficiency**: By introducing a simplified and less computationally costly smoothness indicator, the newly proposed WENO scheme can operate efficiently on non - uniform grids and maintain high - order accuracy.
2. **Ensure accuracy**: Verify that the proposed WENO scheme can achieve the highest - order accuracy when the data does not cross singular points, and exhibit the essentially non - oscillatory property when encountering discontinuities.
3. **Wide applicability**: Apply these improved WENO schemes in the finite - volume method to solve the numerical solutions of hyperbolic conservation laws and hyperbolic systems, especially their performance in shock problems.
### Key technical points
- **Smoothness indicator**: A new simplified smoothness indicator is proposed, whose computational complexity grows linearly with the order of accuracy, rather than the traditional quadratic growth. This significantly reduces the computational cost.
\[
I_{R,i,h} = \sum_{j = i}^{r + i - 1}\left(\frac{f_{j + 1,h}-f_{j,h}}{c_{j + 1}-c_j}\right)^2
\]
- **Weight design**: Based on the simplified smoothness indicator, nonlinear weights are designed to ensure that appropriate reconstruction results can be obtained in both smooth and discontinuous regions.
\[
\omega_{R,i,h}=\frac{\alpha_{R,i,h}}{\sum_{j = 0}^{r'}\alpha_{R,j,h}}, \quad \text{where} \quad \alpha_{R,i,h}=\frac{1}{r' + 1}\left(1+\left(\frac{d_{R,h}}{I_{R,i,h}+\varepsilon}\right)^s\right)^{-1}
\]
- **Global average weight**: A global average weight is introduced to further improve the accuracy of the scheme near discontinuities.
\[
\omega_{R,h}=\frac{1}{1 + d_{R,h}J_{R,h,s}}, \quad \text{where} \quad J_{R,h,s}=\sum_{i = 0}^{r'}\frac{1}{I_{R,i,h}^s+\varepsilon}
\]
### Application and verification
The authors verify the effectiveness of this scheme through a series of numerical experiments, including the numerical solutions of algebraic problems and hyperbolic conservation laws. The experimental results show that the new scheme not only achieves the expected high - order accuracy in smooth regions, but also exhibits good non - oscillatory properties in discontinuous regions.
In conclusion, the main contribution of this paper is to propose an efficient WENO scheme suitable for non - uniform grids, which solves the problems of high computational complexity and accuracy loss of traditional WENO schemes on non - uniform grids.