On maximum-principle-satisfying high order schemes for scalar conservation laws

We construct uniformly high order accurate schemes satisfying a strict maximum principle for scalar conservation laws. A general framework (for arbitrary order of accuracy) is established to construct a limiter for finite volume schemes (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO) schemes) or discontinuous Galerkin (DG) method with first order Euler forward time discretization solving one-dimensional scalar conservation laws. Strong stability preserving (SSP) high order time discretizations will keep the maximum principle. It is straightforward to extend the method to two and higher dimensions on rectangular meshes. We also show that the same limiter can preserve the maximum principle for DG or finite volume schemes solving two-dimensional incompressible Euler equations in the vorticity stream-function formulation, or any passive convection equation with an incompressible velocity field. Numerical tests for both the WENO finite volume scheme and the DG method are reported.     

On positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations on rectangular meshes

We construct uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for Euler equations of compressible gas dynamics. The same framework also applies to high order accurate finite volume (e.g. essentially non-oscillatory (ENO) or weighted ENO (WENO)) schemes. Motivated by Perthame and Shu (1996) [20] and Zhang and Shu (2010) [26], a general framework, for arbitrary order of accuracy, is established to construct a positivity preserving limiter for the finite volume and DG methods with first order Euler forward time discretization solving one-dimensional compressible Euler equations. The limiter can be proven to maintain high order accuracy and is easy to implement. Strong stability preserving (SSP) high order time discretizations will keep the positivity property. Following the idea in Zhang and Shu (2010) [26], we extend this framework to higher dimensions on rectangular meshes in a straightforward way. Numerical tests for the third order DG method are reported to demonstrate the effectiveness of the methods.     

Local DG method using WENO type limiters for convection–diffusion problems

The local discontinuous Galerkin (LDG) method is a spatial discretization procedure for convection–diffusion equations, which employs useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers serving as numerical fluxes and limiters, which is termed as Runge–Kutta LDG (RKLDG) when TVD Runge–Kutta method is applied for time discretization. It has the advantage of flexibility in handling complicated geometry, h-p adaptivity, and efficiency of parallel implementation and has been used successfully in many applications. However, the limiters used to control spurious oscillations in the presence of strong shocks are less robust than the strategies of essentially non-oscillatory (ENO) and weighted ENO (WENO) finite volume and finite difference methods. In this paper, we investigated RKLDG methods with WENO and Hermite WENO (HWENO) limiters for solving convection–diffusion equations on unstructured meshes, with the goal of obtaining a robust and high order limiting procedure to simultaneously obtain uniform high order accuracy and sharp, non-oscillatory shock transition. Numerical results are provided to illustrate the behavior of these procedures.     

一种高精度的修正Hermite-ENO格式

设计一种基于三单元具有六阶精度的修正Hermite-ENO格式（CHENO）,求解一维双曲守恒律问题.CHENO格式利用有限体积法进行空间离散,在空间层上,使用ENO格式中的Newton差商法自适应选择模板.在重构半节点处的函数值及其一阶导数值时,利用Taylor展开给出修正Hermite插值使其提高到六阶精度,并设计了间断识别法与相应的处理方法以抑制间断处的虚假振荡;在时间层上采用三阶TVD Runge-Kutta法进行函数值及一阶导数值的推进.其主要优点是在达到高阶精度的同时具有紧致性.数值实验表明对一维双曲守恒律问题的求解达到了理论分析结果,是有效可行的.     

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.     

Conservative high order semi-Lagrangian finite difference WENO methods for advection in incompressible flow

In this paper, we propose a semi-Lagrangian finite difference formulation for approximating conservative form of advection equations with general variable coefficients. Compared with the traditional semi-Lagrangian finite difference schemes [5], [25], which approximate the advective form of the equation via direct characteristics tracing, the scheme proposed in this paper approximates the conservative form of the equation. This essential difference makes the proposed scheme naturally conservative for equations with general variable coefficients. The proposed conservative semi-Lagrangian finite difference framework is coupled with high order essentially non-oscillatory (ENO) or weighted ENO (WENO) reconstructions to achieve high order accuracy in smooth parts of the solution and to capture sharp interfaces without introducing spurious oscillations. The scheme is extended to high dimensional problems by Strang splitting. The performance of the proposed schemes is demonstrated by linear advection, rigid body rotation, swirling deformation, and two dimensional incompressible flow simulation in the vorticity stream-function formulation. As the information is propagating along characteristics, the proposed scheme does not have the CFL time step restriction of the Eulerian method, allowing for a more efficient numerical realization for many application problems.     

一类自适应运动网格的ALE方法

从积分形式的二维Lagrange流体力学方程组出发,使用ENO高阶插值多项式,推广了四边形结构网格下的一阶有限体积格式,构造一类结构网格下的高精度有限体积格式.结合有效的守恒重映方法,发展一类高精度的ALE方法,并结合自适应运动网格技术,进行ALE方法的数值模拟,得到预期的效果.     

ADER schemes for the shallow water equations in channel with irregular bottom elevation

This paper deals with the construction of high-order ADER numerical schemes for solving the one-dimensional shallow water equations with variable bed elevation. The non-linear version of the schemes is based on ENO reconstructions. The governing equations are expressed in terms of total water height, instead of total water depth, and discharge. The ENO polynomial interpolation procedure is also applied to represent the variable bottom elevation. ADER schemes of up to fifth order of accuracy in space and time for the advection and source terms are implemented and systematically assessed, with particular attention to their convergence rates. Non-oscillatory results are obtained for discontinuous solutions both for the steady and unsteady cases. The resulting schemes can be applied to solve realistic problems characterized by non-uniform bottom geometries.     

A High-Order WENO Finite Difference Scheme for the Equations of Ideal Magnetohydrodynamics

We present a high-order accurate weighted essentially non-oscillatory (WENO) finite difference scheme for solving the equations of ideal magnetohydrodynamics (MHD). This scheme is a direct extension of a WENO scheme, which has been successfully applied to hydrodynamic problems. The WENO scheme follows the same idea of an essentially non-oscillatory (ENO) scheme with an advantage of achieving higher-order accuracy with fewer computations. Both ENO and WENO can be easily applied to two and three spatial dimensions by evaluating the fluxes dimension-by-dimension. Details of the WENO scheme as well as the construction of a suitable eigen-system, which can properly decompose various families of MHD waves and handle the degenerate situations, are presented. Numerical results are shown to perform well for the one-dimensional Brio–Wu Riemann problems, the two-dimensional Kelvin–Helmholtz instability problems, and the two-dimensional Orszag–Tang MHD vortex system. They also demonstrate the importance of maintaining the divergence free condition for the magnetic field in achieving numerical stability. The tests also show the advantages of using the higher-order scheme. The new 5th-order WENO MHD code can attain an accuracy comparable with that of the second-order schemes with many fewer grid points.     

加权ENO格式的构造及数值模拟

将加权ENO格式推广到非结构三角形网格上,构造了一类加权ENO有限体积格式,提出的插值多项式的构造方式,可以减少计算时间.对于出现的病态方程组,给出了解决方法.此外还给出了插值点的选取方式及加权因子的构造方法.结合三阶TVD Runge Kutta时间离散,对二维欧拉方程组进行了数值试验.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

A high order ENO conservative Lagrangian type scheme for the compressible Euler equations

We develop a class of Lagrangian type schemes for solving the Euler equations of compressible gas dynamics both in the Cartesian and in the cylindrical coordinates. The schemes are based on high order essentially non-oscillatory (ENO) reconstruction. They are conservative for the density, momentum and total energy, can maintain formal high order accuracy both in space and time and can achieve at least uniformly second-order accuracy with moving and distorted Lagrangian meshes, are essentially non-oscillatory, and have no parameters to be tuned for individual test cases. One and two-dimensional numerical examples in the Cartesian and cylindrical coordinates are presented to demonstrate the performance of the schemes in terms of accuracy, resolution for discontinuities, and non-oscillatory properties.     

惯性约束聚变激光烧蚀驱动内界面不稳定性数值模拟

利用水平集方法和三阶本质无振荡(ENO)方法研究了流体力学界面不稳定性和激光烧蚀驱动内界面不稳定性数值模拟,并用数值算例与LARED-S程序的运算结果进行了比较.数值结果表明此方法计算惯性约束聚变激光烧蚀驱动内界面不稳定性的有效性和精确性.     

High order conservative Lagrangian schemes with Lax–Wendroff type time discretization for the compressible Euler equations

In this paper, we explore the Lax–Wendroff (LW) type time discretization as an alternative procedure to the high order Runge–Kutta time discretization adopted for the high order essentially non-oscillatory (ENO) Lagrangian schemes developed in 3 and 5. The LW time discretization is based on a Taylor expansion in time, coupled with a local Cauchy–Kowalewski procedure to utilize the partial differential equation (PDE) repeatedly to convert all time derivatives to spatial derivatives, and then to discretize these spatial derivatives based on high order ENO reconstruction. Extensive numerical examples are presented, for both the second-order spatial discretization using quadrilateral meshes [3] and third-order spatial discretization using curvilinear meshes [5]. Comparing with the Runge–Kutta time discretization procedure, an advantage of the LW time discretization is the apparent saving in computational cost and memory requirement, at least for the two-dimensional Euler equations that we have used in the numerical tests.     

Finite element form of FDV for widely varying flowfields

We present the Flowfield Dependent Variation (FDV) method for physical applications that have widely varying spatial and temporal scales. Our motivation is to develop a versatile numerical method that is accurate and stable in simulations with complex geometries and with wide variations in space and time scales. The use of a finite element formulation adds capabilities such as flexible grid geometries and exact enforcement of Neumann boundary conditions. While finite element schemes are used extensively by researchers solving computational fluid dynamics in many engineering fields, their use in space physics, astrophysical fluids and laboratory magnetohydrodynamic simulations with shocks has been predominantly overlooked. The FDV method is unique in that numerical diffusion is derived from physical parameters rather than traditional artificial viscosity methods. Numerical instabilities account for most of the difficulties when capturing shocks in these regimes. The first part of this paper concentrates on the presentation of our numerical method formulation for Newtonian and relativistic hydrodynamics. In the second part we present several standard simulation examples that test the method’s limitations and verify the FDV method. We show that our finite element formulation is stable and accurate for a range of both Mach numbers and Lorentz factors in one-dimensional test problems. We also present the converging/diverging nozzle which contains both incompressible and compressible flow in the flowfield over a range of subsonic and supersonic regions. We demonstrate the stability of our method and the accuracy by comparison with the results of other methods including the finite difference Total Variation Diminishing method. We explore the use of FDV for both non-relativistic and relativistic fluids (hydrodynamics) with strong shocks in order to establish the effectiveness in future applications of this method in astrophysical and laboratory plasma environments.     

自相似Euler方程的数值方法

Euler方程某些问题的解具有自相似特点，可以使用更准确的方法求解.提出了两种数值方法，分别称为自相似和准自相似方法，新方法可以使用现有守恒律方程的数值格式，无须设计特殊方法.对一维激波管问题、二维Riemann问题、激波反射以及激波折射问题进行了数值计算.对自相似Euler方程，一维计算结果显示数值解基本等同于精确解，二维结果也比现有文献计算的结果有更高的分辨率.对准自相似Euler方程，新方法可以求解不具有自相似性但接近自相似的问题，并在计算时间足够长时可以取得自相似Euler方程的效果.数值求解自相似Euler方程对自相似问题的研究，高分辨率、高精度格式的设计乃至Euler方程的精确解都有重要启示.      

一维Gross-Pitaevskii方程的高阶紧致分裂步多辛格式

首先把一维Gross-Pitaevskli方程改写成多辛Hamiltonian系统的形式,把形式通过分裂变成2个子哈密尔顿系统.然后,对这些子系统用辛或者多辛算法进行离散.通过对子系统数值算法的不同组合方式,得到不同精度的具有多辛算法特征数值格式.这些格式不仅具有多辛格式、分裂步方法和高阶紧致格式的特征,而且是质量守恒的.数值实验验证了新格式的数值行为.     

气相爆轰高阶中心差分-WENO组合格式自适应网格方法

研究一种高阶中心差分-WENO组合格式,并采用自适应网格方法进行二维和三维气相爆轰波的数值模拟.采用ZND爆轰模型的控制方程为包含化学反应源项的Euler方程组.组合格式在大梯度区采用WENO格式捕捉间断,在光滑区采用高阶中心差分格式提高计算效率.采用一种基于流场结构特征的自适应网格.计算结果,表明这种方法同时具有高精度、高分辨率和高效率的特点.     

不同层显式格式及在微尺度热传导中的应用

基于指数拟合和牛顿插值多项式构造了用于数值计算偏微分方程的不同层显式格式,应用不同层显式格式可获得不同精度的数值计算结果,并将该算法应用于微结构热传导方程和薄膜强瞬态热传导方程中．