应用于计算气动声学的优化有限紧致格式

对于含间断的计算气动声学问题,数值计算的格式不仅要求低耗散低色散的设计,对短波具有较高的分辨率,还要求能捕捉激波.中心紧致格式具有高精度,具有无耗散和低色散特征,但不能捕捉间断和激波;WENO格式处理间断较为成功,而耗散和色散误差相对较大.有限紧致格式可以将紧致格式与WENO格式相结合构造成混合格式,利用光滑因子之间的关系对激波区域进行自动判断,将传统的全域求解的紧致格式划分为有限的局部紧致求解,间断点上的激波捕捉铜梁自动作为局部紧致求解的边界通量,在在光滑区域具有紧致格式的高精度低耗散性质,在激波附近不产生非物理振荡.本文利用有限紧致格式思想,构造了新的适合于气动声学问题的优化有限紧致格式,将其应用于计算气动声学一维标准测试问题,对相关格式的模拟性能进行了评估,显示该格式在宽频声波传播和含有间断的声波传播模拟方面具有优势.     

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

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

二维溃坝绕流的WENO格式数值模拟

WENO(Weighted essentially non-oscillatory scheme)格式是一类新的高精度无振荡差分格式.本文将WENO格式和Runge-Kutta时间离散的思想应用于二维浅水方程组的求解,数值模拟矩形河道中大坝瞬间局部溃倒,下游有障碍物的洪水演进过程,并对模拟结果进行了分析,表明采用WENO格式所建立的高分辨率模型能够有效地模拟溃坝波的演进过程.     

二维双曲型守恒律的一种五阶松弛格式

松弛格式是Jin和Xin提出的无振荡有限差分方法,其主要思想是将守恒律转化为松弛方程组进行求解.本文用逐维五阶WENO重构和显隐式Runge-Kutta方法对松弛方程组的空间和时间进行离散,得到了一种求解二维双曲型守恒律五阶松弛格式.所得格式保持了松弛格式简单的优点,不用求解Riemann问题和计算通量函数的雅可比矩阵.通过二维Burgers方程和二维浅水方程的数值算例验证了格式的有效性.     

对流扩散方程的摄动有限体积(PFV)方法及讨论

在有限体积(FV)方法的重构近似中,引入数值摄动处理,即把界面数值通量摄动展开成网格间距的幂级数,并利用积分方程自身的性质求出幂级数的系数,同时获得高精度迎风和中心型摄动有限体积(PFV)格式.对标量输运方程给出积分近似为二阶、重构近似为二、三和四阶迎风和中心型PFV格式,这些PFV格式的结构形式及使用基点数与一阶迎风格式完全一致,迎风PFV格式满足对流有界准则;二阶和四阶中心PFV格式对网格Peclet数的任意值均为正型格式,比常用的二阶中心格式优越.用一维标量输运和方腔流动算例说明PFV格式的优良性能,并把PFV方法与性质相近的摄动有限差分(PFD)方法及相关的高精度方法作了对比分析.     

基于高阶和高效格式的悬停旋翼可压缩无粘绕流的计算

发展了一种基于高精度和高效格式计算悬停旋翼复杂绕流的隐式有限体积方法。控制方程为Euler方程,其中对流项通量的左右状态量采用五阶加权基本无振荡(WENO)格式重构,时间推进应用隐式LU-SGS算法,为进一步加速收敛,采用三层V循环多重网格松弛方法。考虑到多重网格方法的思想以及五阶WENO格式涉及6个网格单元,建议仅在最细网格上使用WENO格式。计算结果表明本文方法能有效捕捉激波,对尾迹也有较高分辨率,基于隐式LU-SGS算法的多重网格方法能有效提高计算效率。     

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

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

一种改进的中心-WENO混合格式

基于中心差分与WENO格式混合可以改善WENO格式耗散特性的思想,在理论推导的基础上,给出了一种用于激波捕捉计算的守恒型中心-WENO混合格式,该混合格式可视为三阶WENO格式和二阶中心差分格式的加权平均。在数值研究现有加权函数的基础上,给出了适用于该混合格式的加权函数,使其能够自适应地调整数值耗散以捕捉激波间断。数值结果表明:与三阶WENO格式相比,混合格式HY3_4能够降低数值耗散,更陡峭地捕捉间断,对复杂流场结构具有较高的分辨率;混合格式HY3_5对于包含高压比激波间断流场结构,能给出无振荡、低耗散的结果。     

加权基本无振荡格式研究进展

加权型基本无振荡WENO格式是近十年发展起来的一类高阶、高精度格式，它是在ENO格式的基础上采用加权思想构造的，对流场内的间断和细致结构具有较高的分辨率，适于求解包含激波、膨胀波以及接触间断等复杂结构的流场，目前已发展成为计算流体力学中一种重要的方法。本文针对加权型基本无振荡格式近年来的进展作一简要介绍。     

基于简单WENO-RKDG的溃坝问题数值模拟

对溃坝问题水流间断面的高精度、高分辨率数值模拟是水动力学的重要内容。简单加权本质无振荡(WENO)限制器由"问题单元"及其相邻单元的解重构"问题单元"的解,从而抑制数值解的非物理振荡,能够很好地模拟间断问题。本文详细介绍了简单WENO限制器的基本原理和过程。将简单WENO限制器-Runge-Kutta间断Galerkin方法应用于二维浅水控制方程的求解中,对二维矩形明渠中大坝瞬间全溃、局部溃塌所致的水流运动进行了数值模拟,并将数值计算结果与理论分析进行了比较。计算结果表明,方法能够清晰地捕捉到溃坝全过程中的间断,没有非物理的振荡现象发生,简单WENO限制器-RKDG方法能够很好地模拟溃坝波的演进过程。     

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

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

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.     

Linear stability of WENO schemes coupled with explicit Runge–Kutta schemes

This paper focuses on the results of the linear stability analysis of the finite‐difference weighted essentially non‐oscillatory (WENO) schemes with optimal weights. The standard WENO schemes between the third and 11th order, the order‐optimised WENO schemes of the sixth and eighth order and the bandwidth‐optimised WENO schemes of the third and fourth order are considered. Several explicit Runge–Kutta schemes including the recently published strong stability‐preserving explicit Runge–Kutta schemes are considered for time discretisation. The stability limits as well as dissipation and dispersion properties dependent on the Courant–Friedrichs–Lewy number are presented for a hyperbolic model equation. The different combinations of space and time discretisation schemes are compared in terms of their accuracy and efficiency. For a parabolic model equation, the viscous term is discretised with high‐order central differences. The stability limits for the parabolic problem are presented as well. Numerical results of linear test cases are shown. Copyright © 2011 John Wiley & Sons, Ltd.     

Preventing numerical oscillations in the flux‐split based finite difference method for compressible flows with discontinuities,II

Problems in the characteristic‐wise flux‐split based finite difference method when compressible flows with contact discontinuities or material interfaces are computed were presented and analyzed. The current analysis showed the following: (i) Even with the local characteristic decomposition technique, numerical errors could be caused by point‐wise flux vector splitting (FVS) methods, such as the Steger–Warming FVS or the van Leer FVS. Therefore, the Lax–Friedrichs type FVS method is required. (ii) If the isobars of a material are vertical lines, the combination of using the local characteristic decomposition and the global Lax–Friedrichs FVS can avoid velocity and pressure oscillations of contact discontinuities in this material for weighted essentially non‐oscillatory (WENO) schemes. (iii) For problems with material interfaces, the quasi‐conservative approach can be realized using characteristic‐wise flux‐split based finite difference WENO schemes if nonlinear WENO schemes in genuinely nonlinear characteristic fields can be guaranteed to be the same and the decomposition equation representing material interfaces is discretized properly. Copyright © 2015 John Wiley & Sons, Ltd.     

Well‐balanced finite difference weighted essentially non‐oscillatory schemes for the blood flow model

The blood flow model maintains the steady‐state solutions, in which the flux gradients are non‐zero but exactly balanced by the source term. In this paper, we design high order finite difference weighted essentially non‐oscillatory (WENO) schemes to this model with such well‐balanced property and at the same time keeping genuine high order accuracy. Rigorous theoretical analysis as well as extensive numerical results all indicate that the resulting schemes verify high order accuracy, maintain the well‐balanced property, and keep good resolution for smooth and discontinuous solutions. Copyright © 2016 John Wiley & Sons, Ltd.     

Third- and fourth-order well-balanced schemes for the shallow water equations based on the CWENO reconstruction

High-order finite volume schemes for conservation laws are very useful in applications, due to their ability to compute accurate solutions on quite coarse meshes and with very few restrictions on the kind of cells employed in the discretization. For balance laws, the ability to approximate up to machine precision relevant steady states allows the scheme to compute accurately, also on coarse meshes, small perturbations of such states, which are very relevant for many applications. In this paper, we propose third- and fourth-order accurate finite volume schemes for the shallow water equations. The schemes have the well-balanced property thanks to a path-conservative approach applied to an appropriate nonconservative reformulation of the equations. High-order accuracy is achieved by designing truly two-dimensional (2D) reconstruction procedures of the central WENO (CWENO ) type. The novel schemes are tested for accuracy and well-balancing and shown to maintain positivity of the water height on wet/dry transitions. Finally, they are applied to simulate the Tohoku 2011 tsunami event.     

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.     

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