基于非结构四边形网格的WENO格式

基于非结构四边形网格发展求解双曲守恒律的三阶加权基本无振荡（WENO）格式.针对任意非结构四边形网格选取重构模板，并给出基于线性多项式的三阶线性重构.但对于一般的非结构四边形网格，会出现非常大的线性权和负权，使得非线性重构的WENO格式对光滑问题也不稳定.本文给出一个处理非常大的线性权的优化重构方法，对优化后得到的负线性权采用分裂方法进行处理.对于非线性权，提出一种考虑局部网格和物理量间断的新光滑度量因子.采用优化重构方法和新的非线性权，当前的三阶WENO格式在质量很差的网格上也具有很好的稳定性.理论的三阶精度在数值精度测试算例中得到验证，同时一范数和无穷范数的误差绝对值不依赖于网格质量；具有强间断的数值结果证明了当前格式的有效性.     

基于WENO重构的熵稳定格式求解浅水方程

针对浅水方程,提出一种数值求解格式:空间方向采用满足熵稳定条件的数值通量,并在单元交界面处进行高阶WENO重构,时间上的推进采用强稳定的Runge-Kutta方法.模拟一维和二维经典问题,结果表明,该格式具有分辨率高、基本无振荡性等特点.     

多车种LWR交通流模型的半离散中心迎风格式

对多车种LWR交通流模型,给出一种半离散中心迎风格式,该格式以五阶WENO-Z重构和半离散中心迎风数值通量为基础.WENO-Z重构方法的引入提高了格式的精度,并保证格式具有基本无振荡的性质.时间的离散采用保持强稳定性的Runge-Kutta方法.通过数值算例验证了格式的有效性.     

求解粘性Hamilton-Jacobi方程的高阶方法

提出了求解具有粘性项的Hamilton-Jacobi方程的二阶、四阶方法.该方法以加权基本无振荡(WENO)格式为基础,通过修正数值通量函数和构造右端粘性项的基于非线性限制器的二阶近似、基于Taylor展开的四阶近似,成功地求解了一维、二维的粘性Hamilton-Jacobi方程.给出的算例说明了本方法具有高分辨率、鲁棒性和无振荡特性.     

求解双曲守恒律方程的高分辨率熵稳定格式

熵稳定格式从物理概念出发,保证总熵关于时间耗散,在计算过程中无需进行熵修正,有效避免如膨胀激波,负压力等非物理现象,显示出独特的优点.通过插入限制器和在单元交界面处进行高阶重构,得到一类高分辨率的熵稳定格式.算例结果表明,格式具有可靠性,高精度和基本无振荡性等特点.     

求解双曲守恒律方程的高分辨率熵相容格式

为提高熵相容格式的精度,利用限制器机制构造高分辨率格式,将构造的通量限制器插入熵相容格式,得到一类高分辨率熵相容格式.构造Euler方程高分辨率熵相容格式时,对熵相容格式中的几个参数做简单调整,提高了接触间断处的分辨率.将所得格式的数值结果与熵相容格式的数值结果比较表明,构造的高分辨率熵相容格式具有稳健和基本无振荡等特性.     

一种新的可计算可压缩流动的预处理方法

针对使用可压缩流动数值方法求解不可压缩流动存在的刚性问题,基于虚拟压缩法思想,构造了一种以Mach数、速度、密度、温度等变量为元素的预处理矩阵,改变了控制方程组的特征根并使其量级更接近.通过理论推导与分析,证明新方法相比Weiss, Pletcher, Dailey和Choi的方法而言,不仅能降低方程组的刚性,提高了数值求解效率,而且拥有更好的稳定性,此外还能实现低速流动和高速流动之间的光滑过渡.采用有限差分格式进行离散,对流项的Roe格式作为基本加权无振荡(WENO)格式的求解器,黏性项则使用中心型紧致差分格式来计算,与预处理矩阵相结合展开数值实验,结果表明新预处理方法可以实现对无黏和有黏不可压缩流动问题的高精度模拟,且拥有比Weiss和Pletcher等提出的方法更好的收敛性和稳定性.     

双曲型方程无振荡选取(NOS)有限差分格式

从几何观点解释了双曲型方程差分格式的TVD条件,导出了常用二阶差分格式的无振荡条件,发展了一种具有时空三阶精度的无振荡选取NOS差分格式.从单个双曲型方程的一些典型算例,显示了该格式高精度、无振荡和逻辑简单的特点,并能有效避免通常使用维数分裂法向二维推广时带来的空间耗散不对称性.     

计算流体力学中的高精度数值方法回顾

在过去的二、三十年中,计算流体力学(CFD)领域的高精度数值方法的设计和应用研究非常活跃.高精度数值方法主要针对具有复杂解结构流场的模拟而设计.回顾CFD中主要用于可压缩流模拟的几类高精度格式的发展与应用.可压缩流的一个重要特征是流场中存在激波、界面以及其它间断,同时还常常在解的光滑区域包含复杂结构.这对设计既不振荡又保持高阶精度的格式带来特别的挑战.重点讨论本质无振荡(ENO)、加权本质无振荡(WENO)有限差分与有限体积格式、间断Galerkin有限元(DG)方法,描述它们各自的特点、长处与不足,简要回顾这些方法的发展和应用,重点介绍它们近五年来的最新进展.     

三维不可压N-S方程的多重网格求解

应用全近似存储(Full Approximation Storage,FAS)多重网格法和人工压缩性方法求解了三维不可压Navi-er-Stokes方程.在解粗网格差分方程时,对Neumann边界条件采用增量形式进行更新,离散方程用对角化形式的近似隐式因子分解格式求解,其中空间无粘项分别用MUSCL格式和对称TVD格式进行离散.对90°弯曲的方截面管道流动和4:1椭球体层流绕流的数值模拟表明,多重网格的计算时间比单重网格节省一半以上,且无限制函数的MUSCL格式比TVD格式对流动结构有更好的分辨能力.     

Hybrid weighted essentially non-oscillatory schemes with different indicators

A key idea in finite difference weighted essentially non-oscillatory (WENO) schemes is a combination of lower order fluxes to obtain a higher order approximation. The choice of the weight to each candidate stencil, which is a nonlinear function of the grid values, is crucial to the success of WENO schemes. For the system case, WENO schemes are based on local characteristic decompositions and flux splitting to avoid spurious oscillation. But the cost of computation of nonlinear weights and local characteristic decompositions is very high. In this paper, we investigate hybrid schemes of WENO schemes with high order up-wind linear schemes using different discontinuity indicators and explore the possibility in avoiding the local characteristic decompositions and the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong shocks. The idea is to identify discontinuity by an discontinuity indicator, then reconstruct numerical flux by WENO approximation in discontinuous regions and up-wind linear approximation in smooth regions. These indicators are mainly based on the troubled-cell indicators for discontinuous Galerkin (DG) method which are listed in the paper by Qiu and Shu (J. Qiu, C.-W. Shu, A comparison of troubled-cell indicators for Runge–Kutta discontinuous Galerkin methods using weighted essentially non-oscillatory limiters, SIAM Journal of Scientific Computing 27 (2005) 995–1013). The emphasis of the paper is on comparison of the performance of hybrid scheme using different indicators, with an objective of obtaining efficient and reliable indicators to obtain better performance of hybrid scheme to save computational cost. Detail numerical studies in one- and two-dimensional cases are performed, addressing the issues of efficiency (less CPU time and more accurate numerical solution), non-oscillatory property.     

An Eulerian–Lagrangian WENO finite volume scheme for advection problems

We develop a locally conservative Eulerian–Lagrangian finite volume scheme with the weighted essentially non-oscillatory property (EL–WENO) in one-space dimension. This method has the advantages of both WENO and Eulerian–Lagrangian schemes. It is formally high-order accurate in space (we present the fifth order version) and essentially non-oscillatory. Moreover, it is free of a CFL time step stability restriction and has small time truncation error. The scheme requires a new integral-based WENO reconstruction to handle trace-back integration. A Strang splitting algorithm is presented for higher-dimensional problems, using both the new integral-based and pointwise-based WENO reconstructions. We show formally that it maintains the fifth order accuracy. It is also locally mass conservative. Numerical results are provided to illustrate the performance of the scheme and verify its formal accuracy.     

WENO schemes on arbitrary mixed-element unstructured meshes in three space dimensions

The paper extends weighted essentially non-oscillatory (WENO) methods to three dimensional mixed-element unstructured meshes, comprising tetrahedral, hexahedral, prismatic and pyramidal elements. Numerical results illustrate the convergence rates and non-oscillatory properties of the schemes for various smooth and discontinuous solutions test cases and the compressible Euler equations on various types of grids. Schemes of up to fifth order of spatial accuracy are considered.     

High-order conservative finite difference GLM–MHD schemes for cell-centered MHD

We present and compare third- as well as fifth-order accurate finite difference schemes for the numerical solution of the compressible ideal MHD equations in multiple spatial dimensions. The selected methods lean on four different reconstruction techniques based on recently improved versions of the weighted essentially non-oscillatory (WENO) schemes, monotonicity preserving (MP) schemes as well as slope-limited polynomial reconstruction. The proposed numerical methods are highly accurate in smooth regions of the flow, avoid loss of accuracy in proximity of smooth extrema and provide sharp non-oscillatory transitions at discontinuities.     

Positivity-preserving high order finite difference WENO schemes for compressible Euler equations

In Zhang and Shu (2010) [20], Zhang and Shu (2011) [21] and Zhang et al. (in press) [23], we constructed uniformly high order accurate discontinuous Galerkin (DG) and finite volume schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics. In this paper, we present an extension of this framework to construct positivity-preserving high order essentially non-oscillatory (ENO) and weighted essentially non-oscillatory (WENO) finite difference schemes for compressible Euler equations. General equations of state and source terms are also discussed. Numerical tests of the fifth order finite difference WENO scheme are reported to demonstrate the good behavior of such schemes.     

An improved weighted essentially non-oscillatory scheme for hyperbolic conservation laws

In this article we develop an improved version of the classical fifth-order weighted essentially non-oscillatory finite difference scheme of [G.S. Jiang, C.W. Shu, Efficient implementation of weighted ENO schemes, J. Comput. Phys. 126 (1996) 202–228] (WENO-JS) for hyperbolic conservation laws. Through the novel use of a linear combination of the low order smoothness indicators already present in the framework of WENO-JS, a new smoothness indicator of higher order is devised and new non-oscillatory weights are built, providing a new WENO scheme (WENO-Z) with less dissipation and higher resolution than the classical WENO. This new scheme generates solutions that are sharp as the ones of the mapped WENO scheme (WENO-M) of Henrick et al. [A.K. Henrick, T.D. Aslam, J.M. Powers, Mapped weighted essentially non-oscillatory schemes: achieving optimal order near critical points, J. Comput. Phys. 207 (2005) 542–567], however with a 25% reduction in CPU costs, since no mapping is necessary. We also provide a detailed analysis of the convergence of the WENO-Z scheme at critical points of smooth solutions and show that the solution enhancements of WENO-Z and WENO-M at problems with shocks comes from their ability to assign substantially larger weights to discontinuous stencils than the WENO-JS scheme, not from their superior order of convergence at critical points. Numerical solutions of the linear advection of discontinuous functions and nonlinear hyperbolic conservation laws as the one dimensional Euler equations with Riemann initial value problems, the Mach 3 shock–density wave interaction and the blastwave problems are compared with the ones generated by the WENO-JS and WENO-M schemes. The good performance of the WENO-Z scheme is also demonstrated in the simulation of two dimensional problems as the shock–vortex interaction and a Mach 4.46 Richtmyer–Meshkov Instability (RMI) modeled via the two dimensional Euler equations.     

高阶FD-WENO格式在数值求解Rayleigh-Taylor不稳定性问题中的应用

用高阶加权本质上无振荡有限差分格式(FD-WENO),求解重力作用下高密度比二维流体界面Rayleigh-Taylor不稳定性问题及激光烧蚀Rayleigh-Taylor不稳定性问题,均获得较为理想的数值结果.     

High-order kinetic flux vector splitting schemes in general coordinates for ideal quantum gas dynamics

A class of high-order kinetic flux vector splitting schemes are presented for solving ideal quantum gas dynamics based on quantum statistical mechanics. The collisionless quantum Boltzmann equation approach is adopted and both Bose–Einstein and Fermi–Dirac gases are considered. The formulas for the split flux vectors are derived based on the general three-dimensional distribution function in velocity space and formulas for lower dimensions can be directly deduced. General curvilinear coordinates are introduced to treat practical problems with general geometry. High-order accurate schemes using weighted essentially non-oscillatory methods are implemented. The resulting high resolution kinetic flux splitting schemes are tested for 1D shock tube flows and shock wave diffraction by a 2D wedge and by a circular cylinder in ideal quantum gases. Excellent results have been obtained for all examples computed.     

基于气动声学问题的高阶非线性紧致格式

文章基于线性中心紧致差分格式, 通过非线性加权插值的方法来求解网格中心处的函数值.这类格式保持了原有中心紧致差分格式的高阶精度和低耗散特性, 同时其分辨率也非常高, 由于其非线性插值的机制, 使得这类格式能够捕捉强激波, 所以这类新的高阶非线性紧致格式是一种较好的模拟湍流和气动声学等多尺度问题的方法.      

Monotonicity Preserving Weighted Essentially Non-oscillatory Schemes with Increasingly High Order of Accuracy

In this paper we design a class of numerical schemes that are higher-order extensions of the weighted essentially non-oscillatory (WENO) schemes of G.-S. Jiang and C.-W. Shu (1996) and X.-D. Liu, S. Osher, and T. Chan (1994). Used by themselves, the schemes may not always be monotonicity preserving but coupled with the monotonicity preserving bounds of A. Suresh and H. T. Huynh (1997) they perform very well. The resulting monotonicity preserving weighted essentially non-oscillatory (MPWENO) schemes have high phase accuracy and high order of accuracy. The higher-order members of this family are almost spectrally accurate for smooth problems. Nevertheless, they, have robust shock capturing ability. The schemes are stable under normal CFL numbers. They are also efficient and do not have a computational complexity that is substantially greater than that of the lower-order members of this same family of schemes. The higher accuracy that these schemes offer coupled with their relatively low computational complexity makes them viable competitors to lower-order schemes, such as the older total variation diminishing schemes, for problems containing both discontinuities and rich smooth region structure. We describe the MPWENO schemes here as well as show their ability to reach their designed accuracies for smooth flow. We also examine the role of steepening algorithms such as the artificial compression method in the design of very high order schemes. Several test problems in one and two dimensions are presented. For multidimensional problems where the flow is not aligned with any of the grid directions it is shown that the present schemes have a substantial advantage over lower-order schemes. It is argued that the methods designed here have great utility for direct numerical simulations and large eddy simulations of compressible turbulence. The methodology developed here is applicable to other hyperbolic systems, which is demonstrated by showing that the MPWENO schemes also work very well on magnetohydrodynamical test problems.