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.     

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

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

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 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.     

Third‐order WENO scheme with a new smoothness indicator

In this article, we have devised a new reference smoothness indicator for third‐order weighted essentially non‐oscillatory (WENO) scheme to achieve desired order of convergence at critical points. In the context of the weighted essentially non‐oscillatory scheme, reference smoothness indicator is constructed in such a way that it satisfies the sufficient condition on the weights for the third‐order convergence. The goal is to construct a reference smoothness indicator such that the resulted scheme have to achieve the required order of accuracy even if the first two derivatives vanish but not the third derivative. The construction of such reference smoothness indicator is not possible through a linear combination of local smoothness indicators only. We have proposed a reference smoothness indicator to be of the fourth order of accuracy on three‐point stencil that contains the linear combination of the first derivative information of the local and global stencils. The performance enhancement of the WENO scheme through this reference smoothness indicator is verified through the standard numerical experiments. Numerical results indicate that the new scheme provides better results in comparison with the earlier third‐order WENO schemes like WENO‐JS and WENO‐Z. Copyright © 2017 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.     

三阶WENO格式精度分析与改进

高分辨率激波捕捉格式对含激波流场的数值模拟具有重要意义。通过理论推导分析了WENO-JS3格式和WENO-Z3格式的精度,发现两种格式在光滑流场区域（包含极值点处）具有相同的理论精度且均低于三阶设计精度,WENO-Z3格式由于增大了非光滑模板的非线性权重使其计算精度有所提高。在理论推导的基础上,提出了WENO-Z+3格式及其改进格式（WENO-Z+3P1和WENO-Z+3P32）,且改进格式在光滑流场区域能满足所设计的三阶精度要求。选用一维平面黎曼问题及双马赫反射等经典算例,验证了本文提出的WENO-Z+3格式及其改进格式相较其他格式具有耗散低和对流场结构分辨率高的特性。     

无粘可压缩流动的改进高精度方法

为更准确捕捉复杂流场的流动细节,通过对WENO格式的光滑因子进行改进,发展了一种新的五阶WENO格式。对三阶ENO格式进行加权可以得到五阶WENO格式,但是不同的加权处理,WENO格式在极值处保持加权基本无振荡的效果不同,本文构造了二阶精度的局部光滑因子,及不含一阶二阶导数的高阶全局光滑因子,从而实现WENO格式在极值处有五阶精度。基于改进五阶WENO格式,对一维对流方程、一维和二维可压缩无粘问题进行算例验证,并与传统WENO-JS格式和WENO-Z格式进行比较。计算结果表明,改进五阶WENO格式有较高的精度和收敛速度,有较低的数值耗散,能有效捕捉间断、激波和涡等复杂流动。     

Convergence issues in using high‐resolution schemes and lower–upper symmetric Gauss–Seidel method for steady shock‐induced combustion problems

This paper reports numerical convergence study for simulations of steady shock‐induced combustion problems with high‐resolution shock‐capturing schemes. Five typical schemes are used: the Roe flux‐based monotone upstream‐centered scheme for conservation laws (MUSCL) and weighted essentially non‐oscillatory (WENO) schemes, the Lax–Friedrichs splitting‐based non‐oscillatory no‐free parameter dissipative (NND) and WENO schemes, and the Harten–Yee upwind total variation diminishing (TVD) scheme. These schemes are implemented with the finite volume discretization on structured quadrilateral meshes in dimension‐by‐dimension way and the lower–upper symmetric Gauss–Seidel (LU–SGS) relaxation method for solving the axisymmetric multispecies reactive Navier–Stokes equations. Comparison of iterative convergence between different schemes has been made using supersonic combustion flows around a spherical projectile with Mach numbers M = 3.55 and 6.46 and a ram accelerator with M = 6.7. These test cases were regarded as steady combustion problems in literature. Calculations on gradually refined meshes show that the second‐order NND, MUSCL, and TVD schemes can converge well to steady states from coarse through fine meshes for M = 3.55 case in which shock and combustion fronts are separate, whereas the (nominally) fifth‐order WENO schemes can only converge to some residual level. More interestingly, the numerical results show that all the schemes do not converge to steady‐state solutions for M = 6.46 in the spherical projectile and M = 6.7 in the ram accelerator cases on fine meshes although they all converge on coarser meshes or on fine meshes without chemical reactions. The result is based on the particular preconditioner of LU–SGS scheme. Possible reasons for the nonconvergence in reactive flow simulation are discussed.Copyright © 2012 John Wiley & Sons, Ltd.     

A high performance fifth-order multistep WENO scheme

Many efforts have been made to improve the accuracy of the conventional weighted essentially nonoscillatory (WENO) scheme at transition points (connecting a smooth region and a discontinuity point). This paper analyzes these works and further develops a more effective multistep WENO scheme. Theoretical analysis and numerical results show that the new scheme not only improves the accuracy by one order higher than the traditional fifth-order WENO schemes at transition point but also maintains the fifth-order accuracy in smooth regions even at critical point where the first derivative vanishes.     

加权型紧致格式与加权本质无波动格式的比较

线性紧致格式和加权本质无波动格式是两种典型的高阶精度数值格式,它们各有优缺点.线性紧致格式在具有高阶精度的同时,格式的分辨率也比较高,耗散低,是计算多尺度流场结构的较好格式,但是不能计算具有强激波的流场.加权本质无波动格式是一种高阶精度捕捉激波格式,鲁棒性好,但耗散比较高,分辨率也不理想.近年来,在莱勒的线性紧致格式基础上,采用加权本质无波动格式捕捉激波思想,发展了一系列加权型紧致格式.本文较全面地比较了加权型紧致格式和加权本质无波动格式,包括构造方法、鲁棒性、分辨率、耗散特性、收敛特性以及并行计算效率.结果表明,现有的加权型紧致格式基本保持了加权本质无波动格式的性质,对于气动力等宏观量的计算,比加权本质无波动格式没有明显的优势.      

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.     

A parameter-free ε-adaptive algorithm for improving weighted compact nonlinear schemes

In this paper, we propose a parameter-free algorithm to calculate ε, a parameter of small quantity initially introduced into the nonlinear weights of weighted essentially nonoscillatory (WENO) scheme to avoid denominator becoming zero. The new algorithm, based on local smoothness indicators of fifth-order weighted compact nonlinear scheme (WCNS), is designed in a manner to adaptively increase ε in smooth areas to reduce numerical dissipation and obtain high-order accuracy, and decrease ε in discontinuous areas to increase numerical dissipation and suppress spurious numerical oscillations. We discuss the relation between critical points and discontinuities and illustrate that, when large gradient areas caused by high-order critical points are not well resolved with sufficiently small grid spacing, numerical oscillations arise. The new algorithm treats high-order critical points as discontinuities to suppress numerical oscillations. Canonical numerical tests are carried out, and computational results indicate that the new adaptive algorithm can help improve resolution of small scale flow structures, suppress numerical oscillations near discontinuities, and lessen susceptibility to flux functions and interpolation variables for fifth-order WCNS. The new adaptive algorithm can be conveniently generalized to WENO/WCNS with different orders.     

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.     

A novel weighting switch function for uniformly high‐order hybrid shock‐capturing schemes

Hybrid schemes are very efficient for complex compressible flow simulation. However, for most existing hybrid schemes in literature, empirical problem‐dependent parameters are always needed to detect shock waves and hence greatly decrease the robustness and accuracy of the hybrid scheme. In this paper, based on the nonlinear weights of the weighted essentially non‐oscillatory (WENO) scheme, a novel weighting switch function is proposed. This function approaches 1 with high‐order accuracy in smooth regions and 0 near discontinuities. Then, with the new weighting switch function, a seventh‐order hybrid compact‐reconstruction WENO scheme (HCCS) is developed. The new hybrid scheme uses the same stencil as the fifth‐order WENO scheme, and it has seventh‐order accuracy in smooth regions even at critical points. Numerical tests are presented to demonstrate the accuracy and robustness of both the switch function and HCCS. Comparisons also reveal that HCCS has lower dissipation and less computational cost than the seventh‐order WENO scheme. Copyright © 2016 John Wiley & Sons, Ltd.     

On the performance of WENO/TENO schemes to resolve turbulence in DNS/LES of high‐speed compressible flows

High‐speed compressible turbulent flows typically contain discontinuities and have been widely modeled using Weighted Essentially Non‐Oscillatory (WENO) schemes due to their high‐order accuracy and sharp shock capturing capability. However, such schemes may damp the small scales of turbulence and result in inaccurate solutions in the context of turbulence‐resolving simulations. In this connection, the recently developed Targeted Essentially Non‐Oscillatory (TENO) schemes, including adaptive variants, may offer significant improvements. The present study aims to quantify the potential of these new schemes for a fully turbulent supersonic flow. Specifically, DNS of a compressible turbulent channel flow with M = 1.5 and Re_τ = 222 is conducted using OpenSBLI, a high‐order finite difference computational fluid dynamics framework. This flow configuration is chosen to decouple the effect of flow discontinuities and turbulence and focus on the capability of the aforementioned high‐order schemes to resolve turbulent structures. The effect of the spatial resolution in different directions and coarse grid implicit LES are also evaluated against the WALE LES model. The TENO schemes are found to exhibit significant performance improvements over the WENO schemes in terms of the accuracy of the statistics and the resolution of the three‐dimensional vortical structures. The sixth‐order adaptive TENO scheme is found to produce comparable results to those obtained with nondissipative fourth‐ and sixth‐order central schemes and reference data obtained with spectral methods. Although the most computationally expensive scheme, it is shown that this adaptive scheme can produce satisfactory results if used as an implicit LES model.     

Local discontinuous Galerkin method for a nonlocal viscous conservation laws

The purpose of this article is to investigate high‐order numerical approximations of scalar conservation laws with nonlocal viscous term. The viscous term is given in the form of convolution in space variable. With the help of the characteristic of viscous term, we design a semidiscrete local discontinuous Galerkin (LDG) method to solve the nonlocal model. We prove stability and convergence of semidiscrete LDG method in L² norm. The theoretical analysis reveals that the present numerical scheme is stable with optimal convergence order for the linear case, and it is stable with sub‐optimal convergence order for nonlinear case. To demonstrate the validity and accuracy of our scheme, we test the Burgers equation with two typical nonlocal fractional viscous terms. The numerical results show the convergence order accuracy in space for both linear and nonlinear cases. Some numerical simulations are provided to show the robustness and effectiveness of the present numerical scheme.     

非稳态Hamilton-Jacobi方程的7阶加权紧致非线性格式

Hamilton-Jacobi (HJ) 方程是一类重要的非线性偏微分方程, 在物理学、流体力学、图像处理、微分几何、金融数学、最优化控制理论等方面有着广泛的应用. 由于HJ方程的弱解存在但不唯一, 且解的导数可能出现间断, 导致其数值求解具有一定的难度. 本文提出了非稳态HJ方程的7阶精度加权紧致非线性格式 (WCNS). 该格式结合了Hamilton函数的Lax-Friedrichs型通量分裂方法和一阶空间导数左、右极限值的高阶精度混合节点和半节点型中心差分格式. 基于7点全局模板和4个4点子模板推导了半节点函数值的高阶线性逼近和4个低阶线性逼近, 以及全局模板和子模板的光滑度量指标. 为避免间断附近数值解产生非物理振荡以及提高格式稳定性, 采用WENO型非线性插值方法计算半节点函数值. 时间离散采用3阶TVD型Runge-Kutta方法. 通过理论分析验证了WCNS格式对于光滑解具有最佳的7阶精度. 为方便比较, 经典的7阶WENO格式也被推广用于求解HJ方程. 数值结果表明, 本文提出的WCNS格式能够很好地模拟HJ方程的精确解, 且在光滑区域能够达到7阶精度; 与经典的同阶WENO格式相比, WCNS格式在精度、收敛性和分辨率方面更优, 计算效率略高.      

基于混合型WENO格式数值模拟四类Riemann问题

本文基于三阶WENO格式和三阶WENO-Z格式,利用有限体积法研究了同一格式在不同方向使用不同的求解器,以及在不同方向采用相同精度的不同格式数值模拟4类Riemann问题,分析各向异性对数值计算结果的影响．数值结果表明,无论是使用不同的格式,还是使用不同的求解器进行数值模拟,都会导致数值结果不同程度地失去图像对称性．