A new smoothness indicator for third‐order WENO scheme

We present a new reference smoothness indicator for third‐order weighted essentially non‐oscillatory scheme to recover its design‐order convergence at critical points. This reference smoothness indicator, which involves both the candidate and global smoothness indicators in the weighted essentially non‐oscillatory framework, is devised according to a sufficient condition on the weights for third‐order convergence. The recovery of design‐order is verified by standard tests. Meanwhile, numerical results demonstrate that the present reference smoothness indicator produces sharper representation of the discontinuity owing to the combined effects of larger weight assignment to the discontinuous stencils and convergence rate recovery. Copyright © 2015 John Wiley & Sons, Ltd.     

Improved third‐order weighted essentially nonoscillatory scheme

A new third‐order WENO scheme is proposed to achieve the desired order of convergence at the critical points for scalar hyperbolic equations. A new reference smoothness indicator is introduced, which satisfies the sufficient condition on the weights for the third‐order convergence. Following the truncation error analysis, we have shown that the proposed scheme achieves the desired order accurate for smooth solutions with arbitrary number of vanishing derivatives if the parameter ε satisfies certain conditions. We have made a comparative study of the proposed scheme with the existing schemes such as WENO‐JS, WENO‐Z, and WENO‐N3 through different numerical examples. The result shows that the proposed scheme (WENO‐MN3) achieves better performance than these schemes.     

A low-dissipation third-order weighted essentially nonoscillatory scheme with a new reference smoothness indicator

The classical third-order weighted essentially nonoscillatory (WENO) scheme is notoriously dissipative as it loses the optimal order of accuracy at critical points and its two-point finite difference in the smoothness indicators is unable to differentiate the critical point from the discontinuity. In recent years, modifications to the smoothness indicators and weights of the classical third-order WENO scheme have been reported to reduce numerical dissipation. This article presents a new reference smoothness indicator for constructing a low-dissipation third-order WENO scheme. The new reference smoothness indicator is a nonlinear combination of the local and global stencil smoothness indicators. The resulting WENO-Rp3 scheme with the power parameter p=1.5 achieves third-order accuracy in smooth regions including critical points and has low dissipation, but numerical results show this scheme cannot keep the ENO property near discontinuities. The recommended WENO-R3 scheme (p=1) keeps the ENO property and performs better than several recently developed third-order WENO schemes.     

Improved third‐order weighted essentially nonoscillatory schemes with new smoothness indicators

In this article, we present two improved third‐order weighted essentially nonoscillatory (WENO) schemes for recovering their design‐order near first‐order critical points. The schemes are constructed in the framework of third‐order WENO‐Z scheme. Two new global smoothness indicators, τ_L3 and τ_L4, are devised by a nonlinear combination of local smoothness indicators (IS_k) and reference values (IS_G) based on Lagrangian interpolation polynomial. The performances of the proposed schemes are evaluated on several numerical tests governed by one‐dimensional linear advection equation or one‐ and two‐dimensional Euler equations. Numerical results indicate that the presented schemes provide less dissipation and higher resolution than the original WENO3‐JS and subsequent WENO3‐N scheme.     

Novel local smoothness indicators for improving the third‐order WENO scheme

The local smoothness indicators play an important role in the performance of a weighted essentially nonoscillatory (WENO) scheme. Due to having only 2 points available on each substencil, the local smoothness indicators calculated by conventional methods make the third‐order WENO scheme too dissipative. In this paper, we propose a different method to calculate the indicators by using all the 3 points on the global stencil of the third‐order WENO scheme. The numerical results demonstrate that the WENO scheme with the new indicators has less dissipation and better resolution than the conventional third‐order WENO scheme of Jiang and Shu for both smooth and discontinuous solutions.     

求解双曲守恒律的修正模板近似的五阶WENO格式

针对经典的五阶加权本质无振荡(WENO)格式在间断附近耗散过大以及临界点不能保精度的问题,本文提出了一种新的修正模板近似方法。改进了经典五阶WENO-JS格式中各候选子模板上数值通量的二阶多项式逼近,通过加入三次修正项使模板逼近达到四阶精度,并且通过引入可调函数φ使得新的格式具有ENO性质,理论分析新的格式具有保精度特性,通过一系列数值算例说明了新格式的高效性。     

求解双曲守恒律方程的WENO型熵相容格式

通过在单元交界面处进行高阶WENO重构,得到了一种求解双曲型守恒律方程的WENO型熵相容格式。用该格式对一维Burgers方程和Euler方程进行数值模拟,结果表明,该格式具有高精度、基本无振荡性等特点。

     

L1-type smoothness indicators based weighted essentially nonoscillatory scheme for Hamilton-Jacobi equations

This article presents an improved fifth-order finite difference weighted essentially nonoscillatory (WENO) scheme to solve Hamilton-Jacobi equations. A new type of nonlinear weights is introduced with the construction of local smoothness indicators on each local stencil that are measured with the help of generalized undivided differences in L¹-norm. A novel global smoothness measurement is also constructed with the help of local measurements from its linear combination. Numerical experiments are conducted in one- and two-dimensions to demonstrate the performance enhancement, resolution power, numerical accuracy for the proposed scheme, and compared it with the classical WENO scheme.     

An improved third-order finite difference weighted essentially nonoscillatory scheme for hyperbolic conservation laws

In this article, we present an improved third-order finite difference weighted essentially nonoscillatory (WENO) scheme to promote the order of convergence at critical points for the hyperbolic conservation laws. The improved WENO scheme is an extension of WENO-ZQ scheme. However, the global smoothness indicator has a little different from WENO-ZQ scheme. In this follow-up article, a convex combination of a second-degree polynomial with two linear polynomials in a traditional WENO fashion is used to compute the numerical flux at cell boundary. Although the same three-point information is adopted by the improved third-order WENO scheme, the truncation errors are smaller than some other third-order WENO schemes in L_∞ and L₂ norms. Especially, the convergence order is not declined at critical points, where the first and second derivatives vanish but not the third derivative. At last, the behavior of improved scheme is proved on a variety of one- and two-dimensional standard numerical examples. Numerical results demonstrate that the proposed scheme gives better performance in comparison with other third-order WENO schemes.     

非线性双曲型守恒律的高精度MmB差分格式

构造了一维非线性双曲型守恒律方程的一个高精度、高分辨率的广义G odunov型差分格式。其构造思想是:首先将计算区间划分为若干个互不相交的小区间,再根据精度要求等分小区间,通过各细小区间上的单元平均状态变量,重构各等分小区间交界面上的状态变量,并加以校正;其次,利用近似R iem ann解算子求解细小区间交界面上的数值通量,并结合高阶R unge-K u tta TVD方法进行时间离散,得到了高精度的全离散方法。证明了该格式的Mm B特性。然后,将格式推广到一、二维双曲型守恒方程组情形。最后给出了一、二维Eu ler方程组的几个典型的数值算例,验证了格式的高效性。     

Improvement of the weighted essentially nonoscillatory scheme based on the interaction of smoothness indicators

The weighted essentially nonoscillatory scheme is improved by introducing new smoothness indicators that evaluate the interactions among the classical smoothness indicators suggested by Jiang and Shu. The effect of the key parameters in the new smoothness indicators on the scheme is systematically investigated. The improved scheme has smaller dissipation with larger weight assignment to the discontinuous stencils and higher numerical accuracy with weights closer to the ideal weights. To verify the theory, benchmark problems governed by the linear transport equation, the 1‐dimensional nonlinear Burgers equation, and the Euler equations are conducted and analyzed, respectively. Better computational performances both on numerical resolution and accuracy are shown in the comparisons with other classical weighted essentially nonoscillatory schemes.     

An efficient low‐dissipative WENO filter scheme

This paper presents a simple and efficient procedure developed for tracing discontinuities in flow fields. Numerical experiments are carried out to test the new sensor coupled with the associated nonlinear WENO dissipation filter developed to suppress the numerical dissipation. The tests show that, for a problem containing shocks and vortices, the implementation of the new sensor and an optimized WENO scheme can obtain a stable and highly resolved solution. The numerical experiments demonstrated that the new filter scheme performs efficiently both in parallel and serial running for one‐dimensional inviscid flow problems. Direct numerical simulation of a Mach 5 turbulent boundary layer over a flat plate was carried out to demonstrate the applicability of the scheme to the DNS practices. Copyright © 2011 John Wiley & Sons, Ltd.     

一种提高极值点处收敛精度的三阶WENO-Z格式

为了提高三阶WENO-Z格式在极值点处的计算精度,通过理论推导给出三阶WENO格式满足收敛精度的充分条件。采用泰勒级数展开的方式,推导给出所构造格式非线性权重的计算公式,并综合权衡计算精度和计算稳定性确定所构造格式的参数。通过两个典型的精度测试,验证了改进格式在光滑流场极值点区域逼近三阶精度。进一步选用激波与熵波相互作用和Richtmyer-Meshkov不稳定性等经典算例,证实了本文提出的改进格式WENO-PZ3相较其他格式（WENO-JS3和WENO-Z3）不仅具有较高的精度,而且降低了格式的耗散,提高了对流场结构的分辨率。     

A new weighting method for improving the WENO‐Z scheme

The calculation of the weight of each substencil is very important for a weighted essentially nonoscillatory (WENO) scheme to obtain high‐order accuracy in smooth regions and keep the essentially nonoscillatory property near discontinuities. The weighting function introduced in the WENO‐Z scheme provides a straightforward method to analyze the accuracy order in smooth regions. In this paper, we construct a new sixth‐order global smoothness indicator (GSI‐6) and a function about GSI‐6 and the local smoothness indicators (IS_k) to calculate the weights. The analysis and numerical results show that, with the new weights, the scheme satisfies the sufficient condition for the fifth‐order convergence in smooth regions even at critical points. Meanwhile, it can also maintain low dissipation for discontinuous solutions due to relative large weights assigned to discontinuous substencils.     

Local maximum principle satisfying high‐order non‐oscillatory schemes

The main contribution of this work is to classify the solution region including data extrema for which high‐order non‐oscillatory approximation can be achieved. It is performed in the framework of local maximum principle (LMP) and non‐conservative formulation. The representative uniformly second‐order accurate schemes are converted in to their non‐conservative form using the ratio of consecutive gradients. Using the local maximum principle, these non‐conservative schemes are analyzed for their non‐linear LMP/total variation diminishing stability bounds which classify the solution region where high‐order accuracy can be achieved. Based on the bounds, second‐order accurate hybrid numerical schemes are constructed using a shock detector. The presented numerical results show that these hybrid schemes preserve high accuracy at non‐sonic extrema without exhibiting any induced local oscillations or clipping error. Copyright © 2015 John Wiley & Sons, Ltd.     

Improvement of multistep WENO scheme and its extension to higher orders of accuracy

This paper presents an efficient procedure for overcoming the deficiency of weighted essentially non‐oscillatory schemes near discontinuities. Through a thorough incorporation of smoothness indicators into the weights definition, up to ninth‐order accurate multistep methods are devised, providing weighted essentially non‐oscillatory schemes with enhanced order of convergence at transition points from smooth regions to a discontinuity, while maintaining stability and the essentially non‐oscillatory behavior. We also provide a detailed analysis of the resolution power and show that the solution enhancements of the new method at smooth regions come from their ability to render smoothness indicators closer to uniformity. The new scheme exhibits similar fidelity as other multistep schemes; however, with superior characteristics in terms of robustness and efficiency, as no logical statements or mapping function is needed. Extensions to higher orders of accuracy present no extra complexity. Numerical solutions of linear advection problems and nonlinear hyperbolic conservation laws are used to demonstrate the scheme's improved behavior for shock‐capturing problems. Copyright © 2016 John Wiley & Sons, Ltd.     

Numerical experiments with several variant WENO schemes for the Euler equations

Numerical experiments with several variants of the original weighted essentially non‐oscillatory (WENO) schemes (J. Comput. Phys. 1996; 126 :202–228) including anti‐diffusive flux corrections, the mapped WENO scheme, and modified smoothness indicator are tested for the Euler equations. The TVD Runge–Kutta explicit time‐integrating scheme is adopted for unsteady flow computations and lower–upper symmetric‐Gauss–Seidel (LU‐SGS) implicit method is employed for the computation of steady‐state solutions. A numerical flux of the variant WENO scheme in flux limiter form is presented, which consists of first‐order and high‐order fluxes and allows for a more flexible choice of low‐order schemes. Computations of unsteady oblique shock wave diffraction over a wedge and steady transonic flows over NACA 0012 and RAE 2822 airfoils are presented to test and compare the methods. Various aspects of the variant WENO methods including contact discontinuity sharpening and steady‐state convergence rate are examined. By using the WENO scheme with anti‐diffusive flux corrections, the present solutions indicate that good convergence rate can be achieved and high‐order accuracy is maintained and contact discontinuities are sharpened markedly as compared with the original WENO schemes on the same meshes. Copyright © 2008 John Wiley & Sons, Ltd.     

High‐order ADER‐WENO ALE schemes on unstructured triangular meshes—application of several node solvers to hydrodynamics and magnetohydrodynamics

In this paper, we present a class of high‐order accurate cell‐centered arbitrary Lagrangian–Eulerian (ALE) one‐step ADER weighted essentially non‐oscillatory (WENO) finite volume schemes for the solution of nonlinear hyperbolic conservation laws on two‐dimensional unstructured triangular meshes. High order of accuracy in space is achieved by a WENO reconstruction algorithm, while a local space–time Galerkin predictor allows the schemes to be high order accurate also in time by using an element‐local weak formulation of the governing PDE on moving meshes. The mesh motion can be computed by choosing among three different node solvers, which are for the first time compared with each other in this article: the node velocity may be obtained either (i) as an arithmetic average among the states surrounding the node, as suggested by Cheng and Shu, or (ii) as a solution of multiple one‐dimensional half‐Riemann problems around a vertex, as suggested by Maire, or (iii) by solving approximately a multidimensional Riemann problem around each vertex of the mesh using the genuinely multidimensional Harten–Lax–van Leer Riemann solver recently proposed by Balsara et al. Once the vertex velocity and thus the new node location have been determined by the node solver, the local mesh motion is then constructed by straight edges connecting the vertex positions at the old time level tⁿ with the new ones at the next time level t^{n + 1}. If necessary, a rezoning step can be introduced here to overcome mesh tangling or highly deformed elements. The final ALE finite volume scheme is based directly on a space–time conservation formulation of the governing PDE system, which therefore makes an additional remapping stage unnecessary, as the ALE fluxes already properly take into account the rezoned geometry. In this sense, our scheme falls into the category of direct ALE methods. Furthermore, the geometric conservation law is satisfied by the scheme by construction. We apply the high‐order algorithm presented in this paper to the Euler equations of compressible gas dynamics as well as to the ideal classical and relativistic magnetohydrodynamic equations. We show numerical convergence results up to fifth order of accuracy in space and time together with some classical numerical test problems for each hyperbolic system under consideration. Copyright © 2014 John Wiley & Sons, Ltd.     

High‐order ENO and WENO schemes for unstructured grids

This work describes the implementation and analysis of high‐order accurate schemes applied to high‐speed flows on unstructured grids. The class of essentially non‐oscillatory schemes (ENO), that includes weighted ENO schemes (WENO), is discussed in the paper with regard to the implementation of third‐ and fourth‐order accurate methods. The entire reconstruction process of ENO and WENO schemes is described with emphasis on the stencil selection algorithms. The stencils can be composed by control volumes with any number of edges, e.g. triangles, quadrilaterals and hybrid meshes. In the paper, ENO and WENO schemes are implemented for the solution of the dimensionless, 2‐D Euler equations in a cell centred finite volume context. High‐order flux integration is achieved using Gaussian quadratures. An approximate Riemann solver is used to evaluate the fluxes on the interfaces of the control volumes and a TVD Runge–Kutta scheme provides the time integration of the equations. Such a coupling of all these numerical tools, together with the high‐order interpolation of primitive variables provided by ENO and WENO schemes, leads to the desired order of accuracy expected in the solutions. An adaptive mesh refinement technique provides better resolution in regions with strong flowfield gradients. Results for high‐speed flow simulations are presented with the objective of assessing the implemented capability. Copyright © 2007 John Wiley & Sons, Ltd.     

一类高精度TVD差分格式及其应用

构造了一维非线性双曲型守恒律的一个新的高精度、高分辨率的守恒型TvD差分格式。其构造思想是：首先，将计算区间划分为若干个互不相交的小区间，再根据精度要求等分小区间，通过各细小区间上的单元平均状态变量，重构各细小区间交界面上的状态变量，并加以校正；其次，利用近似Riemann解计算细小区间交界面上的数值通量，并结合高阶Runge—Kutta TVD方法进行时间离散，得到了高精度的全离散方法。证明了该格式的TVD特性。该格式适合于使用分量形式计算而无须进行局部特征分解。通过计算几个典型的问题，验证了格式具有高精度、高分辨率且计算简单的优点。