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

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

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

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

三阶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格式有较高的精度和收敛速度,有较低的数值耗散,能有效捕捉间断、激波和涡等复杂流动。     

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

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

求解双曲型守恒律的半离散三阶中心迎风格式

给出了求解一维双曲型守恒律的一种半离散三阶中心迎风格式，并利用逐维进行计算的方法将格式推广到二维守恒律。构造格式时利用了波传播的单侧局部速度，三阶重构方法的引入保证了格式的精度。时间方向的离散采用三阶TVD Runge—Kutta方法。本文格式保持了中心差分格式简单的优点，即不需用Riemann解算器，避免了进行特征分解过程。用该格式对一维和二维守恒律进行了大量的数值试验，结果表明本文格式是高精度、高分辨率的。     

带有沙漏控制的相对自由度壳元非线性动力分析

针对八节点相对自由度壳单元，给出了单元内坐标和位移的插值公式，利用HuWashinzu变分原理，基于拟应变法，在大变形情况下推导了拟应变的表达式，构造了带有沙漏控制的动力问题的有限元求解格式。通过算例表明该文提出的基于相对自由度壳元的沙漏控制算法能够很好地解决非线性动力问题，可改善计算精度和计算效率。     

AN HLLC SCHEME FOR THE SEVEN-EQUATION MULTIPHASE MODEL AND ITS APPLICATION TO COMPRESSIBLE MULTICOMPONENT

针对Saurel和Abgrall提出的两速度两压力的七方程可压缩多相流模型,改进了其数值解法并应用于模拟可压缩多介质流动问题．在Saurel等的算子分裂法基础上,根据Abgrall的多相流系统应满足速度和压力的均匀性不随时间改变的思想,推导了与HLLC格式一致的非守恒项离散格式以及体积分数发展方程的迎风格式．进一步,通过改变分裂步顺序,构造了稳健的结合算子分裂的三阶TVD龙格一库塔方法．最后通过几个一维和二维高密度比高压力比气液两相流算例,显示了该方法在计算精度和稳健性上的改进效果．     

七方程可压缩多相流模型的HLLC格式及应用

针对Saurel和Abgrall提出的两速度两压力的七方程可压缩多相流模型,改进了其数值解法并应用于模拟可压缩多介质流动问题.在Saurel等的算子分裂法基础上,根据Abgrall的多相流系统应满足速度和压力的均匀性不随时间改变的思想,推导了与HLLC格式一致的非守恒项离散格式以及体积分数发展方程的迎风格式.进一步,通过改变分裂步顺序,构造了稳健的结合算子分裂的三阶TVD龙格-库塔方法.最后通过几个一维和二维高密度比高压力比气液两相流算例,显示了该方法在计算精度和稳健性上的改进效果.      

通量限制器对算法性能的影响

在非正交曹线坐标系下,本文给出了求解非线性双曲型Ｅｕｌｅｒ方程的ＬＵ－ＡＵＳＭＬＷ算法。为了改进该算法的性能,将高分辨率ＡＵＳＭＰＷ格式的空间精度由一阶精度扩展到三阶精度。分析了选择通量限制器对算法稳定性、收敛性和精度的影响,并构造了一种新的能量限制器。本文数值结果表明,通量限制器是决定ＬＵ－ＡＵＳＭＰＷ算法性能的关键因素,并且该算法采用本文构造的通量限制器,求解非线性双曲型Ｅｕｌｅｒ方程,具有较     

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.     

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

改进的五阶WENO-Z+格式

WENO-ZWENO-Z$+\!$格式的性能提升依赖于新增项的作用,该项的作用是在WENO-Z格式的基础上进一步增大欠光滑子模板上的权重. 系数$\lambda$被设置用来调控该项的作用, 以避免负耗散. 本文指出了WENO-Z$+\!$格式的缺陷,其所采用$\lambda $的取值方式既不能充分发挥格式的潜力, 也未完全消除负耗散;提出$\lambda $的值应随当地流场数据变化,方能充分发挥新增项在降低数值耗散、提高分辨率上的潜力. 基于此,本文重新设计了$\lambda $的计算公式,该公式能自适应地调控新增项的作用: 只在欠光滑子模板上的权重容易过度增大的地方削弱该项的作用,以避免负耗散; 在其他地方则充分发挥新增项的作用,最大限度增大欠光滑子模板上的权重, 提高格式的分辨率.将使用该系数公式的新格式命名为WENO-Z++, 并对其数值性能进行了系统的研究.理论分析表明, 新格式在间断处具有基本无振荡(essentially non-oscillatory,ENO)特性和更低的数值耗散. 对近似色散关系(approximate dispersion relation,ADR)的研究表明,新格式有效地避免了因过度增大欠光滑子模板上的权重而带来的负耗散,其频谱特性也得到了显著提升.本文还推导了使新格式在极值点处也能保持最优阶的精度的参数设置.一系列求解Euler方程的数值试验表明,新格式的激波捕捉能力和对复杂流场结构的分辨率都显著优于原WENO-Z$+\!$格式.}     

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.     

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.     

用隐式方法和WENO格式计算悬停旋翼跨声速无粘流场

寻找一种能够准确计算以涡为主要特征的复杂流场和克服尾迹耗散问题的数值方法,一直是旋翼空气动力学研究的热点和难点。本文发展了一种基于高阶迎风格式计算悬停旋翼无粘流场的隐式数值方法。无粘通量采用Roe通量差分分裂格式,为提高精度,使用五阶WENO格式进行左右状态插值,并与MUSCL插值进行比较。为提高收敛到定常解的效率,时间推进采用LU-SGS隐式方法。用该方法对一跨声速悬停旋翼无粘流场进行了数值计算,数值结果表明WENO-Roe的激波分辨率高于MUSCL-Roe,体现出了格式精度的提高对计算结果的改善,LU-SGS隐式方法的计算效率比5步Runge-Kutta显式方法的高。