A third-order compact nonlinear scheme for compressible flow simulations

Reynolds-averaged Navier-Stokes simulations based on second-order numerical methods are widely used by commercial codes and work as dominating tools for most industrial applications. They, however, suffer from limitations in accurate and reliable predictions of skin-friction drag and aerodynamic heating, as well as in simulations of complex flows such as large-scale separation and transition. A remedy for this is the development of high-order schemes, by which numerically induced dissipation and dispersion errors of low-order schemes can be effectively reduced. Weighted compact nonlinear schemes (WCNSs) are a family of high-resolution nonlinear shock-capturing methods. A stencil-selection procedure is introduced in the proposed work with an aim to improve the nonlinear weight of the third-order WCNS. By using the approximate dispersion relation analysis, it is demonstrated that the new scheme has reduced dissipation and dispersion errors, compared with WCNSs using two typical nonlinear weights. Improvements are also achieved by the new scheme in numerical tests such as the double Mach reflection problem and the Rayleigh-Taylor instability simulation, which are characterized by strong shock discontinuities and rich small scales, respectively. The new scheme is therefore highly favored in the simulation of flow problems involving strong discontinuities and multiscales phenomena.     

非稳态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格式在精度、收敛性和分辨率方面更优, 计算效率略高.      

Implicit large eddy simulation of a supersonic flat-plate boundary layer flow by weighted compact nonlinear scheme

The performance of implicit large eddy simulation (ILES) of a supersonic flat-plate turbulent boundary layer flow by weighted compact nonlinear scheme (WCNS) has been investigated. In view of features of WCNS and ILES, it was expected that ILES by WCNS could be an efficient approach to perform LES of supersonic turbulent flows. The flowfield calculated by WCNS was of lower turbulent intensity compared with an explicit LES data obtained by a numerical scheme of the same order of accuracy on a computational grid of similar resolution. It was concluded that the numerical dissipation inherent in WCNS is so large that applying WCNS to ILES of this flowfield is inefficient compared with explicit LES.     

尺度无关高阶WCNS格式

高速流场的数值模拟中, 既要保证对小尺度结构的高保真分辨, 又要实现对激波稳定、无振荡地捕捉.当前工程中广泛应用的高精度数值格式虽然都能一定程度地满足上述两种要求, 但仍与理想目标存在较大差距.例如, 模拟雷诺应力模型等小尺度问题时, 高精度格式在间断解附近易产生数值振荡.基于高精度格式所存在的上述问题, 本文引入去尺度函数, 探索了一种更加简单稳定的非线性权重构造方法, 并将其应用于7阶精度加权紧致非线性格式WCNS, 提出了一种尺度无关的7阶WCNS格式.该格式的性能与灵敏度参数和尺度因子的选择无关, 并且在小尺度下仍可以有效捕捉流场激波.同时, 该格式在间断处具有基本无振荡性质, 且在任意尺度函数下保持尺度无关, 并且在极值点处也能保持最优精度.本文还推导了7阶D权函数的形式.最后, 在一维线性对流方程中验证了新格式在流场光滑区能够达到设计精度, 并通过一系列数值实验证明了尺度无关的7阶WCNS格式在激波捕捉能力上具有良好表现, 为WCNS格式改进和解决可压缩湍流等非线性问题提供了一种新途径.      

High-order maximum-principle-preserving and positivity-preserving weighted compact nonlinear schemes for hyperbolic conservation laws

In this paper,the maximum-principle-preserving(MPP)and positivitypreserving(PP)flux limiting technique will be generalized to a class of high-order weighted compact nonlinear schemes(WCNSs)for scalar conservation laws and the compressible Euler systems in both one and two dimensions.The main idea of the present method is to rewrite the scheme in a conservative form,and then define the local limiting parameters via case-by-case discussion.Smooth test problems are presented to demonstrate that the proposed MPP/PP WCNSs incorporating a third-order Runge-Kutta method can attain the desired order of accuracy.Other test problems with strong shocks and high pressure and density ratios are also conducted to testify the performance of the schemes.     

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.     

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.     

基于HWCNS格式的紧致插值方法研究

在保证良好间断捕捉能力的前提下,能够达到更高的分辨率,一直是有限差分方法努力的方向。基于HWCNS格式构造思想,发展了一种高精度非线性紧致插值方法,构造了紧致七阶HWCNS格式,分析了其频谱特性,利用多个典型算例对所构造的格式性能进行了考察。结果表明,在模拟包含间断和多尺度流动结构的长时间演化问题中,本文发展的方法在计算结果精度和综合计算效率方面优于显式五阶HWCNS和七阶WENO格式,与频谱分析结论一致。     

高精度加权紧致非线性格式的研究进展

综述了高精度加权紧致非线性格式WCNS在理论分析以及复杂流动应用方面的研究进展. 首先回顾了国内外高精度格式研究的概况, 然后介绍了WCNS的研究与发展历程. 在对WCNS进行了Fourier分析和渐近稳定性分析后, 给出了WCNS求解多维复杂流动的算例.      

基于非结构/混合网格模拟黏性流的高阶精度DDG/FV混合方法

间断Galerkin有限元方法(discontinuous Galerkin method, DGM) 因具有计算精度高、模板紧致、易于并行等优点, 近年来已成为非结构/混合网格上广泛研究的高阶精度数值方法. 但其计算量和内存需求量巨大, 特别是对于网格规模达到百万甚至数千万的大型三维实际复杂外形问题, 其计算量和存储量对计算资源的消耗是难以承受的. 基于“混合重构”的DG/FV 格式可以有效降低DGM 的计算量和存储量. 本文将DDG 黏性项离散方法推广应用于DG/FV 混合算法, 得到新的DDG/FV混合格式, 以进一步提高DG/FV混合算法对于黏性流动模拟的计算效率. 通过Couette流动、层流平板边界层、定常圆柱绕流, 非定常圆柱绕流和NACA0012 翼型绕流等二维黏性流算例, 优化了DDG 通量公式中的参数选择, 验证了DDG/FV 混合格式对定常和非定常黏性流模拟的精度和计算效率, 并与广泛使用的BR2-DG 格式的计算结果和效率进行对比研究. 一系列数值实验结果表明, 本文构造的DDG/FV混合格式在二维非结构/混合网格的Navier-Stokes 方程求解中, 在达到相同的数值精度阶的前提下, 相比BR2-DG格式, 对于隐式时间离散的定常问题计算效率提高了2 倍以上, 对于显式时间离散的非定常问题计算效率提高1.6 倍, 并且在一些算例中, 混合格式具有更优良的计算稳定性. DDG/FV 混合格式提升了计算效率和稳定性, 具有良好的应用前景.      

Efficient high-order immersed interface methods for heat equations with interfaces

An efficient high-order immersed interface method （IIM） is proposed to solve two-dimensional （2D） heat problems with fixed interfaces on Cartesian grids, which has the fourth-order accuracy in the maximum norm in both time and space directions. The space variable is discretized by a high-order compact （HOC） difference scheme with correction terms added at the irregular points. The time derivative is integrated by a Crank-Nicolson and alternative direction implicit （ADI） scheme. In this case, the time accuracy is just second-order. The Richardson extrapolation method is used to improve the time accuracy to fourth-order. The numerical results confirm the convergence order and the efficiency of the method.     

Compact Third-Order Multidimensional Upwind Scheme for Navier–Stokes Simulations

A new compact third-order scheme for the solution of the unsteady Navier--Stokes equations on unstructured grids is proposed. The scheme is a cell-based algorithm, belonging to the class of Multidimensional Upwind schemes, which uses a finite-element reconstruction procedure over the cell to achieve third order (spatial) accuracy. Derivation of the scheme is given. The asymptotic accuracy, for steady/unsteady inviscid or viscous flow situations, is proved using numerical experiments. These results are compared with the performances of a second-order multidimensional upwind scheme. The new compact high-order discretization proves to have excellent parallel scalability, which makes it well suited for large-scale computations on parallel supercomputers. Our studies show clearly the advantages of the new compact third-order scheme compared with the classical second-order Multidimensional Upwind scheme. Received 29 October 2001 and accepted 21 March 2002     

Large-eddy simulation of wing tip vortex in the near field

Large-eddy simulation with filtered-structure-function subgrid model and implicit large-eddy simulation (ILES without explicit subgrid model) using high-order accuracy and high resolution compact scheme have been performed on the tip vortex shedding from a rectangular half-wing with a NACA 0012 airfoil section and a rounded wing tip. The formation of the tip vortex and its initial development in the boundary layer and the near field wake are investigated and analysed in detail. The physics, why the tip vortex, which is originally turbulent in the boundary layer, is re-laminarised and becomes stable and laminar rapidly after shedding in the near field, is revealed by this simulation. The computation also shows the widely used second-order subgrid model is not consistent to six-order compact scheme and would degenerate the six-order LES results to second-order. Therefore, high-order schemes, grid refinement and six-order subgrid models are critical to LES approaches.     

薄板弯曲分析的高阶高效无网格法

与传统有限元法相比,无网格法具有节点形函数高度光滑、易于形成高阶近似等优势,更适合于以薄板弯曲问题为代表的高阶偏微分方程的数值求解。然而,高阶无网格法的形函数是非多项式的有理函数,导致弱形式的区域积分难以得到精确计算,通常采用的高阶高斯积分方法需使用大量积分点,计算效率低且精度不高。本文针对薄板弯曲问题的高阶（三阶）无网格法分析,首次发展了与该高阶近似相一致的曲率光顺方案,并基于背景三角形积分单元建立了相应的数值积分格式,大幅度减少了所需的积分点数目。所发展方法的关键在于计算刚度阵所需的形函数的二阶导数由形函数及其一阶导数通过散度定理确定,而非对形函数直接求导获得。数值结果表明,基于标准的高斯积分方案的高阶无网格法精度不高,不能精确再现纯弯曲和线性弯曲模式,且得到的弯矩场分布存在严重的虚假数值振荡。而本文所建议的基于曲率光顺方案的高阶无网格法能够方便高效地求解薄板弯曲问题,尤其是它能精确反映纯弯曲和线性弯曲模式。与标准的高斯积分方法和目前主流的常曲率光顺方法相比,本文方法在计算效率、精度、弯矩分布等方面均展现出显著优势,因而具有较好的应用价值。     

New optimized flux difference schemes for improving high-order weighted compact nonlinear scheme with applications

To improve the spectral characteristics of the high-order weighted compact nonlinear scheme(WCNS), optimized flux difference schemes are proposed. The disadvantages in previous optimization routines, i.e., reducing formal orders, or extending stencil widths, are avoided in the new optimized schemes by utilizing fluxes from both cell-edges and cell-nodes. Optimizations are implemented with Fourier analysis for linear schemes and the approximate dispersion relation(ADR) for nonlinear schemes. Clas...     

高精度格式WCNS-E-5的Fourier分析与应用

对高精度加权紧致非线性格式WCNS-E-5进行了理论分析和应用研究。首先采用Fourier方法分析了WCNS-E-5与WCNS-5及其它高阶格式(迎风偏置格式EUW-5与标准格式)的差分误差特性,研究了它们在分辨效率方面的特性和相位误差在多维方向上的各向异性特性。分析结果显示WCNS-E-5与WCNS-5在色散与耗散特性方面优于EUW-5,分辨效率也普遍高于EUW-5和格式,而且它们的相位速度在多维方向上传播所表现的各向异性特性与其它高阶格式一致。WCNS-E-5的高精度特性与WCNS-5的一致,但在计算过程中少了三对角矩阵求解因而它的效率更高,于是采用WCNS-E-5数值模拟了二维/三维复杂流场,算例包括单涡斜向运动、二维Rie-mann问题以及存在分离的大攻角钝锥流动问题。计算结果体现了WCNS-E-5具有很低的数值耗散误差,它对激波、膨胀波和接触面等间断具有很好的捕捉能力,得到的图像清晰光滑,准确再现了真实流动现象。     

Development and assessment of a variable-order non-oscillatory scheme for convection term discretization

A new scheme for convection term discretization is developed, called VONOS (variable-order non-oscillatory scheme). The development of the scheme is based on the behaviour of well-known non-oscillatory schemes in the pure convection of a step profile test case. The new scheme is a combination of the QUICK and BSOU (bounded second-order upwind) schemes. These two schemes do not have the same formal order of accuracy and for that reason the formal order of accuracy of the new scheme is variable. The scheme is conservative, bounded and accurate. The performance of the new scheme was assessed in three test cases. The results showed that it is more accurate than currently used higher-order schemes, so it can be used in a general purpose algorithm in order to save computational time for the same level of accuracy. © 1998 John Wiley & Sons, Ltd.     

A simultaneous space-time wavelet method for nonlinear initial boundary value problems

A high-precision and space-time fully decoupled numerical method is developed for a class of nonlinear initial boundary value problems. It is established based on a proposed Coiflet-based approximation scheme with an adjustable high order for the functions over a bounded interval, which allows the expansion coefficients to be explicitly expressed by the function values at a series of single points. When the solution method is used, the nonlinear initial boundary value problems are first spatially discretized into a series of nonlinear initial value problems by combining the proposed wavelet approximation and the conventional Galerkin method, and a novel high-order step-by-step time integrating approach is then developed for the resulting nonlinear initial value problems with the same function approximation scheme based on the wavelet theory. The solution method is shown to have the N th-order accuracy, as long as the Coiflet with [0, 3 N-1]compact support is adopted, where N can be any positive even number. Typical examples in mechanics are considered to justify the accuracy and efficiency of the method.     

Hermite梯度重构核近似配点法及其在功能梯度材料板的静力分析中的应用

由于直接配点法在求解边值问题时边界上的求解精度较低,本文提出了Hermite梯度重构核近似配点法（HGCM）来改进边界求解精度。重构核近似是无网格法中一种常用的近似函数,但是其在求解高阶导数时格式复杂且非常耗时。HGCM采用梯度重构核近似构建形函数的任意高阶导数,提高了计算效率;通过Hermite配点法构建离散方程,提高了边界求解精度。这种方法在求解对应变系数四阶偏微分方程的功能梯度材料板的静力问题时精度高,计算效率高,并可进一步推广应用于高阶偏微分方程描述的边值问题。     

THE INFLUENCE OF GLOBAL-DIRECTION STENCIL ON GRADIENT AND HIGH-ORDER DERIVATIVES RECONSTRUCTION OF UNSTRUCTURED FINITE VOLUME METHODS 1)

模板选择方式对非结构有限体积方法的计算准确性会产生显著影响. 在之前的工作中, 基于局部方向模板存在的问题, 我们探索了一种更加简单有效的全局方向模板选择方法, 并将其应用于二阶精度非结构有限体积求解器. 基于该方法找到的模板单元均沿着壁面法向与流向, 可有效捕捉流场变化, 反映流动的各向异性, 并且模板选择过程脱离了对网格拓扑的依赖, 避免了局部方向模板选择方法中复杂的阵面推进与方向判断过程, 克服了在大压缩比三角形网格上模板单元偏离壁面法向的现象, 同时在二阶精度求解器上得到了较高的计算精度与计算准确性. 为了进一步验证全局方向模板在高阶精度非结构有限体积方法中应用的可行性, 本文初步测试了该模板对变量梯度及高阶导数重构的影响. 经检验, 在不同类型的网格上, 采用全局方向模板得到的变量梯度与高阶导数误差明显低于局部方向模板, 同时也低于共点模板的计算误差. 此外, 在高斯积分点处由全局方向模板得到的变量点值与导数误差同样在三种模板中最低. 因此该模板选择方法在非结构有限体积梯度与高阶导数重构方面具有较好的数值表现, 具备在高阶精度非结构有限体积求解器中应用并推广的可行性.