On new strategies to control the accuracy of WENO algorithms close to discontinuities
This paper is devoted to the construction and analysis of new nonlinear optimal weights for weighted ENO (WENO) interpolation capable of raising the order of accuracy close to discontinuities. The new nonlinear optimal weights are constructed using a strategy inspired by the original WENO algorithm,...
| Autores: | , , |
|---|---|
| Tipo de recurso: | artículo |
| Estado: | Versión enviada para evaluación y publicación |
| Fecha de publicación: | 2019 |
| País: | España |
| Institución: | Universidad Politécnica de Cartagena(UPCT) |
| Repositorio: | Repositorio Digital UPCT |
| OAI Identifier: | oai:repositorio.upct.es:10317/10656 |
| Acceso en línea: | http://hdl.handle.net/10317/10656 https://epubs.siam.org/doi/10.1137/18M1214937 |
| Access Level: | acceso abierto |
| Palabra clave: | WENO schemes New optimal weights Improved adaption to discontinuities Signal processing Matemática Aplicada 12 Matemáticas |
| Sumario: | This paper is devoted to the construction and analysis of new nonlinear optimal weights for weighted ENO (WENO) interpolation capable of raising the order of accuracy close to discontinuities. The new nonlinear optimal weights are constructed using a strategy inspired by the original WENO algorithm, and they work very well for corner or jump singularities, leading to optimal theoretical accuracy. This is the first part of a series of two papers. In this first part we analyze the performance of the new algorithms proposed for univariate function approximation in the point values (interpolation problem). In the second part, we will extend the analysis to univariate function approximation in the cell averages (reconstruction problem). Our aim is twofold: to raise the order of accuracy of the WENO type interpolation schemes both near discontinuities and in the interval which contains the singularity. The first problem can be solved using the new nonlinear optimal weights, but the second one requires a new strategy that locates the position of the singularity inside the cell in order to attain adaption. This new strategy is inspired by the ENO-SR schemes proposed by Harten [J. Comput. Phys., 83 (1989), pp. 148--184]. Thus, we will introduce two different algorithms in the point values. The first one can deal with corner singularities and jump discontinuities for intervals not containing the singularity. The second algorithm can also deal with intervals containing corner singularities, as they can be detected from the point values, but jump discontinuities cannot, as the information of their position is lost during the discretization process. As mentioned before, the second part of this work will be devoted to the cell averages and, in this context, it will be possible to work with jump discontinuities as well. |
|---|