ENO schemes with subcell resolution

doi:10.1016/0021-9991(89)90226-X

Journal of Computational Physics

Volume 83, Issue 1, July 1989, Pages 148-184

Dedicated to Eugene Isaacson on his 70th birthday

https://doi.org/10.1016/0021-9991(89)90226-X Get rights and content

Abstract

In this paper, we introduce the notion of subcell resolution, which is based on the observation that unlike point values, cell-averages of a discontinuous piecewise-smooth function contain information about the exact location of the discontinuity within the cell. Using this observation we design an essentially non-oscillatory (ENO) reconstruction technique which is exact for cell averages of discontinuous piecewise-polynomial functions of the appropriate degree. Later on we incorporate this new reconstruction technique into Godunov-type schemes in order to produce a modification of the ENO schemes which prevents the smearing of contact discontinuities.

References (11)

ColellaP. et al.
J. Comput. Phys.
(1984)
HartenA.
J. Comput. Phys.
(1983)
HartenA. et al.
J. Appl. Numer. Math.
(1986)
HartenA. et al.
RoeP.L.
J. Comput. Phys.
(1981)
GodunovS.K.
Math. Sb.
(1959)

There are more references available in the full text version of this article.

Cited by (360)

An improved fifth-order WENO scheme with symmetry-preserving smoothness indicators for hyperbolic conservation laws
2023, Journal of Computational Physics
Ha et al. [1] have constructed a novel set of smoothness indicators by using the $L_{1}$ -norm measure and introducing a user-tunable parameter to propose a new WENO scheme, dubbed WENO-NS. Kim et al. [2] have improved WENO-NS with another user-tunable parameter to balance the contribution of substencils, leading to the WENO-P scheme. Recently, Rathan et al. [3] have modified WENO-P by devising a higher-order global smoothness indicator. The associated MWENO-P scheme aimed to fix the issue of WENO-P and WENO-NS that they cannot achieve the optimal convergence orders at the critical points where the first and second derivatives are zero. The above-said schemes outperform the existing many WENO schemes in some senses. However, their smoothness indicators violate the symmetry-preserving property preserved by those of many other well-established WENO schemes and rely on one or two user-tunable parameters. In the present study, a new set of smoothness indicators satisfying the symmetry-preserving property without using any user-tunable parameters is constructed. The resultant scheme is called the improved MWENO-P scheme and abbreviated as IMWENO-P. We provide detailed theoretical analysis and numerical examples to verify that the new scheme can attain optimal convergence orders even in the presence of second-order critical points. Euler equations are simulated to demonstrate the enhancements of the new scheme, such as the lower dissipation and better resolution.
Locally order-preserving mapping for WENO methods
2023, Journal of Computational and Applied Mathematics
In our previous studies (Li and Zhong, 2021a; Li and Zhong, 2021b), the commonly reported issue that most of the existing mapped WENO schemes suffer from either losing high resolutions or generating spurious oscillations in long-run simulations of hyperbolic problems has been successfully addressed, by devising the improved mapped WENO schemes, namely MOP-WENO-X, where “X” stands for the version of the existing mapped WENO scheme. However, all the MOP-WENO-X schemes bring about the serious deficiency that their resolutions in the region with smoothly varying high-frequency waves are dramatically decreased compared to their associated WENO-X schemes. The purpose of this paper is to overcome this drawback. We first present the locally order-preserving (LOP) mapping. Then, by using a new proposed posteriori adaptive technique, we apply this LOP property to obtain the new mappings from those of the WENO-X schemes. The essential idea of the posteriori adaptive technique is to identify the global stencil in which the existing mappings fail to preserve the LOP property, and then replace the mapped weights with the weights of the classic WENO-JS scheme to recover the LOP property. We build the resultant mapped WENO schemes and denote them as LOP-WENO-X. The numerical experiments demonstrate that the resolutions in the region with high-frequency but smooth waves of the LOP-WENO-X schemes are similar to or even better than those of their associated WENO-X schemes and naturally much higher than the MOP-WENO-X schemes. Furthermore, the LOP-WENO-X schemes gain all the great advantages of the MOP-WENO-X schemes, such as attaining high resolutions and in the meantime preventing spurious oscillations near discontinuities when solving the one-dimensional linear advection problems with long output times, and significantly reducing the post-shock oscillations in the simulations of the two-dimensional problems with shock waves.
On the supremum of the steepness parameter in self-adjusting discontinuity-preserving schemes
2022, Computers and Fluids
Self-adjusting steepness (SAS)-based schemes preserve various structures in the compressible flows. These schemes provide a range of desired behaviors depending on the steepness-adjustable limiters with the steepness measured by a steepness parameter. These properties include either second-order accuracy with exact steepness infima that are theoretically given or having anti-diffusive/compression properties with a larger steepness parameter. Nevertheless, the supremum of the steepness parameter has not been determined theoretically yet. In this study, we demonstrate that any anti-diffusive limiter should be limited by Ultra-bee limiter according to Sweby's total variation diminishing (TVD) condition. Two typical steepness-adjustable limiters are analyzed in detail including the tangent of hyperbola for interface capturing (THINC) limiter and the steepness-adjustable harmonic (SAH) limiter. Applying this constraint, we derive for the first time two inequalities which the steepness parameters much satisfy. Furthermore, we obtain the analytical expression of the Courant-Friedrichs-Lewy (CFL) number-dependent supremum of the steepness parameter. Using this solution, we then propose supremum-determined SAS schemes. These schemes are further extended to solve the compressible Euler equations. The results of typical numerical tests confirm our theoretical conclusions and show that the final schemes are capable of sharply capturing contact discontinuities and minimizing numerical oscillations.
Self-adjusting steepness-based schemes that preserve discontinuous structures in compressible flows
2022, Journal of Computational Physics
Although conventional total-variation-diminishing schemes can capture discontinuities with essentially non-oscillatory properties and are nominal of second-order accuracy in smooth regions, they cannot generally maintain sharp resolution for various discontinuities (especially contact discontinuity for long-time simulations) owing to excessive numerical dissipation. Designing schemes to overcome this disadvantage is a difficult problem. In this paper, we propose to solve this problem by introducing a novel framework of self-adjusting steepness-based schemes, which are mainly constructed in four steps: (1) design a slope limiter containing a steepness parameter β that provides a mechanism to enable the scheme to accurately solve both smooth and discontinuous problems with proper values of β; (2) determine the infimum of β such that the scheme is order-optimized; (3) determine the supremum of β such that the scheme has a nonlinear stable anti-diffusion/compression effect; (4) calculate the steepness parameter β using an adaptive algorithm in terms of the infimum and supremum to ensure that the final scheme obtains essentially non-oscillatory and sharp resolutions for various discontinuities while maintaining the nominal second-order accuracy for smooth regions. Moreover, such a framework has been implemented on the tangent-of-hyperbola-for-interface-capturing scheme (Xiao et al. (2005) [27]) and the classical harmonic limiter (van Leer (1974) [37]), resulting in two specific self-adjusting steepness-based schemes. Finally, these two schemes were tested using a series of numerical examples in one and two dimensions. It was found that the schemes cannot only obtain second-order accuracy in smooth regions but also preserve discontinuous flow structures, especially contact discontinuities, even after long computation times.
Hybrid high-order steepness-adjustable harmonic scheme based on combined discontinuity sensors
2022, Computers and Fluids
Extending discontinuity-capturing schemes to high resolution using a hybrid approach has become popular in recent times. The hybrid schemes reported in literature usually combine high-order compact/explicit-schemes in smooth regions and the weighted essentially non-oscillatory (WENO) scheme in discontinuities by a discontinuity sensor. However, there are two major concerns: (1) the significant dissipation of the WENO scheme will smear the contact discontinuities, and (2) the accuracy of the discontinuity sensor will significantly influence the performance of the hybrid scheme. However, these two aspects are typically coupled. Therefore, in this study, we attempt to replace the discontinuity-capturing schemes (WENO) with discontinuity-preserving schemes (steepness-adjustable harmonic (SAH) scheme). Then, we extend the SAH scheme to high-order accuracy by hybridizing it with a linear fourth-degree polynomial (P4) scheme. The new P4-SAH scheme was employed to evaluate the nine different discontinuity sensors. The numerical results revealed that a single discontinuity sensor can hardly adapt to different numerical cases. Therefore, we propose a new class of discontinuity sensors that combines MP-WGVC discontinuity sensors. Several classical one- and two-dimensional Euler equation benchmark tests were simulated. The numerical results demonstrate that the new hybrid P4-SAH scheme with a combined MP-WGVC sensor not only sharply captures shock waves and contact discontinuities, but also capture small vortex structures.
High order recovery of geometric interfaces from cell-average data
2024, arXiv

View all citing articles on Scopus

: Research was supported under NSF Grant DMS85-03294, DARPA Grant in the ACMP Program, ONR Grant N00014-86-K-0691, NASA Ames Interchange NCA2-185, and NASA Langley Grant NAG1-270. Also research was partially supported under the National Aeronautics and Space Administration under NASA Contract NAS1-18107 while the author was in residence at the Institute for Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, Hampton, VA 23665-5225.

View full text

ENO schemes with subcell resolution

Abstract

J. Comput. Phys.

J. Comput. Phys.

J. Appl. Numer. Math.

J. Comput. Phys.

Math. Sb.