A high‐resolution hybrid scheme for hyperbolic conservation laws

A new hybrid scheme is proposed, which combines the improved third‐order weighted essentially non‐oscillatory (WENO) scheme presented in this paper with a fourth‐order central scheme by a novel switch. Two major steps have been gone through for the construction of a high‐performance and stable hybrid scheme. Firstly, to enhance the WENO part of the hybrid scheme, a new reference smoothness indicator has been devised, which, combined with the nonlinear weighting procedure of WENO‐Z, can drive the third‐order WENO toward the optimal linear scheme faster. Secondly, to improve the hybridization with the central scheme, a hyperbolic tangent hybridization switch and its efficient polynomial counterpart are devised, with which we are able to fix the threshold value introduced by the hybridization. The new hybrid scheme is thus formulated, and a set of benchmark problems have been tested to verify the performance enhancement. Numerical results demonstrate that the new hybrid scheme achieves excellent performance in resolving complex flow features, even compared with the fifth‐order classical WENO scheme and WENO‐Z scheme. Copyright © 2015 John Wiley & Sons, Ltd.     

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.     

High‐order k‐exact WENO finite volume schemes for solving gas dynamic Euler equations on unstructured grids

This paper presents a family of High‐order finite volume schemes applicable on unstructured grids. The k‐exact reconstruction is performed on every control volume as the primary reconstruction. On a cell of interest, besides the primary reconstruction, additional candidate reconstruction polynomials are provided by means of very simple and efficient ‘secondary’ reconstructions. The weighted average procedure of the WENO scheme is then applied to the primary and secondary reconstructions to ensure the shock‐capturing capability of the scheme. This procedure combines the simplicity of the k‐exact reconstruction with the robustness of the WENO schemes and represents a systematic and unified way to construct High‐order accurate shock capturing schemes. To further improve the efficiency, an efficient problem‐independent shock detector is introduced. Several test cases are presented to demonstrate the accuracy and non‐oscillation property of the proposed schemes. The results show that the proposed schemes can predict the smooth solutions with uniformly High‐order accuracy and can capture the shock waves and contact discontinuities in high resolution. Copyright © 2011 John Wiley & Sons, Ltd.     

Robustness of the hybrid DRP‐WENO scheme for shock flow computations

This paper describes a new variant of hybrid scheme that is constructed by a wave‐capturing scheme and a nonoscillatory scheme for flow computations in the presence of shocks. The improved fifth‐order upwind weighted essentially nonoscillatory scheme is chosen to be conjugated with the seven‐point dispersion‐relation‐preserving scheme by means of an adaptive switch function of grid‐point type. The new hybrid scheme can achieve a better resolution than the hybrid scheme which is based on the classical weighted essentially scheme. Ami Harten's multiresolution analysis algorithm is applied to density field for detecting discontinuities and setting point values of the switch function adaptively. Moreover, the tenth‐order central filter is applied in smooth part of the flow field for damping dispersion errors. This scheme can promote overall computational efficiency and yield oscillation‐free results in shock flows. The resolution properties and robustness of the new hybrid scheme are tested in both 1D and 2D linear and nonlinear cases. It performs well for computing flow problems with rich structures of weak/strong shocks and large/small vortices, such as the shock‐boundary layer interaction problem in a shock tube, which illustrates that it is very robust and accurate for direct numerical simulation of gas‐dynamics flows. Copyright © 2011 John Wiley & Sons, Ltd.     

Application of central differencing and low‐dissipation weights in a weighted compact nonlinear scheme

This paper proposes WCNS‐CU‐Z, a weighted compact nonlinear scheme, that incorporates adapted central difference and low‐dissipative weights together with concepts of the adaptive central‐upwind sixth‐order weighted essentially non‐oscillatory scheme (WENO‐CU) and WENO‐Z schemes. The newly developed WCNS‐CU‐Z is a high‐resolution scheme, because interpolation of this scheme employs a central stencil constructed by upwind and downwind stencils. The smoothness indicator of the downwind stencil is calculated using the entire central stencil, and the downwind stencil is stopped around the discontinuity for stability. Moreover, interpolation of the sixth‐order WCNS‐CU‐Z exhibits sufficient accuracy in the smooth region through use of low‐dissipative weights. The sixth‐order WCNS‐CU‐Zs are implemented with a robust linear difference formulation (R‐WCNS‐CU6‐Z), and the resolution and robustness of this scheme were evaluated. These evaluations showed that R‐WCNS‐CU6‐Z is capable of achieving a higher resolution than the seventh‐order classical robust weighted compact nonlinear scheme and can provide a crisp result in terms of discontinuity. Among the schemes tested, R‐WCNS‐CU6‐Z has been shown to be robust, and variable interpolation type R‐WCNS‐CU6‐Z (R‐WCNS‐CU6‐Z‐V) provides a stable computation by modifying the first‐order interpolation when negative density or negative pressure arises after nonlinear interpolation. Copyright © 2016 John Wiley & Sons, Ltd.     

Developing a hybrid flux function suitable for hypersonic flow simulation with high‐order methods

In this paper, we develop a new hybrid Euler flux function based on Roe's flux difference scheme, which is free from shock instability and still preserves the accuracy and efficiency of Roe's flux scheme. For computational cost, only 5% extra CPU time is required compared with Roe's FDS. In hypersonic flow simulation with high‐order methods, the hybrid flux function would automatically switch to the Rusanov flux function near shock waves to improve the robustness, and in smooth regions, Roe's FDS would be recovered so that the advantages of high‐order methods can be maintained. Multidimensional dissipation is introduced to eliminate the adverse effects caused by flux function switching and further enhance the robustness of shock‐capturing, especially when the shock waves are not aligned with grids. A series of tests shows that this new hybrid flux function with a high‐order weighted compact nonlinear scheme is not only robust for shock‐capturing but also accurate for hypersonic heat transfer prediction. Copyright © 2015 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.     

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.     

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.     

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.     

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.     

A semi‐Lagrangian multi‐moment finite volume method with fourth‐order WENO projection

Numerical oscillation has been an open problem for high‐order numerical methods with increased local degrees of freedom (DOFs). Current strategies mainly follow the limiting projections derived originally for conventional finite volume methods and thus are not able to make full use of the sub‐cell information available in the local high‐order reconstructions. This paper presents a novel algorithm that introduces a nodal value‐based weighted essentially non‐oscillatory limiter for constrained interpolation profile/multi‐moment finite volume method (CIP/MM FVM) (Ii and Xiao, J. Comput. Phys., 222 (2007), 849–871) as an effort to pursue a better suited formulation to implement the limiting projection in schemes with local DOFs. The new scheme, CIP‐CSL‐WENO4 scheme, extends the CIP/MM FVM method by limiting the slope constraint in the interpolation function using the weighted essentially non‐oscillatory (WENO) reconstruction that makes use of the sub‐cell information available from the local DOFs and is built from the point values at the solution points within three neighboring cells, thus resulting a more compact WENO stencil. The proposed WENO limiter matches well the original CIP/MM FVM, which leads to a new scheme of high accuracy, algorithmic simplicity, and computational efficiency. We present the numerical results of benchmark tests for both scalar and Euler conservation laws to manifest the fourth‐order accuracy and oscillation‐suppressing property of the proposed scheme. Copyright © 2016 John Wiley & Sons, Ltd.     

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

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

A numerical study for WENO scheme-based on different lattice Boltzmann flux solver for compressible flows

In this paper, the finite difference weighted essentially non-oscillatory (WENO) scheme is incorporated into the recently developed four kinds of lattice Boltzmann flux solver (LBFS) to simulate compressible flows, including inviscid LBFS I, viscous LBFS II, hybrid LBFS III and hybrid LBFS IV. Hybrid LBFS can automatically realize the switch between inviscid LBFS I and viscous LBFS II through introducing a switch function. The resultant hybrid WENO–LBFS scheme absorbs the advantages of WENO scheme and hybrid LBFS. We investigate the performance of WENO scheme based on four kinds of LBFS systematically. Numerical results indicate that the devopled hybrid WENO–LBFS scheme has high accuracy, high resolution and no oscillations. It can not only accurately calculate smooth solutions, but also can effectively capture contact discontinuities and strong shock waves.     

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.     

Improvement of the WENO scheme smoothness estimator

This paper has analyzed the weights of the fifth‐order weighted essentially nonoscillatory scheme suggested by Jiang and Shu and a modified smoothness estimator is suggested to improve the accuracy. Using a function of the minimum and maximum of the original smoothness estimators, the new smoothness estimators have smaller variation among them than the original ones in smooth regions. Hence, as the new weights are closer to the optimal weights, the numerical accuracy is improved. Several numerical experiments are presented to demonstrate the accuracy and robustness of the new scheme. Copyright © 2009 John Wiley & Sons, Ltd.     

A new high-accuracy scheme for compressible turbulent flows

Shock-capturing and broad-bandwidth scale resolutions are two main challenges of compressible turbulent flow simulation. To meet the rigorous requests, a novel fifth-order hybrid scheme based on a uniform hybrid framework is designed. With the help of a continuous weight operator, the new scheme combines an upwind compact scheme for smooth regions and a compact-reconstruction weighted essentially non-oscillatory scheme for discontinuous regions. Numerical analyses and canonical numerical tests confirm that the new scheme has high accuracy, spectral-like resolution property and shock-capturing capability. Besides, the new scheme shows high computational efficiency compared to the related shock-capturing schemes and hybrid ones.     

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

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

A γ‐DGBGK scheme for compressible multi‐fluids

The paper presents a Discontinuous Galerkin γ‐BGK (γ‐DGBGK) method for compressible multicomponent flow simulations by coupling the discontinuous Galerkin method with a γ‐BGK scheme based on WENO limiters. In this γ‐DGBGK method, the construction of the flux in the DG method is based on the kinetic scheme which not only couples the convective and dissipative terms together, but also includes both discontinuous and continuous terms in the flux formulation at cell interfaces. WENO limiters are used to obtain uniform high‐order accuracy and sharp non‐oscillatory shock transition, and time accuracy obtained by integration for the flux function at the cell interface. Numerical examples in one and two space dimensions are presented to illustrate the robust and accuracy of the present scheme. Copyright © 2010 John Wiley & Sons, Ltd.