Positivity-preserving high order discontinuous Galerkin schemes for compressible Euler equations with source terms

In [16], [17], we constructed uniformly high order accurate discontinuous Galerkin (DG) schemes which preserve positivity of density and pressure for the Euler equations of compressible gas dynamics with the ideal gas equation of state. The technique also applies to high order accurate finite volume schemes. For the Euler equations with various source terms (e.g., gravity and chemical reactions), it is more difficult to design high order schemes which do not produce negative density or pressure. In this paper, we first show that our framework to construct positivity-preserving high order schemes in [16], [17] can also be applied to Euler equations with a general equation of state. Then we discuss an extension to Euler equations with source terms. Numerical tests of the third order Runge–Kutta DG (RKDG) method for Euler equations with different types of source terms are reported.     

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.     

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.     

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.     

Positivity preserving semi-Lagrangian discontinuous Galerkin formulation: Theoretical analysis and application to the Vlasov–Poisson system

Semi-Lagrangian (SL) methods have been very popular in the Vlasov simulation community , , , , , ,  and . In this paper, we propose a new Strang split SL discontinuous Galerkin (DG) method for solving the Vlasov equation. Specifically, we apply the Strang splitting for the Vlasov equation [6], as a way to decouple the nonlinear Vlasov system into a sequence of 1-D advection equations, each of which has an advection velocity that only depends on coordinates that are transverse to the direction of propagation. To evolve the decoupled linear equations, we propose to couple the SL framework with the semi-discrete DG formulation. The proposed SL DG method is free of time step restriction compared with the Runge–Kutta DG method, which is known to suffer from numerical time step limitation with relatively small CFL numbers according to linear stability analysis. We apply the recently developed positivity preserving (PP) limiter [37], which is a low-cost black box procedure, to our scheme to ensure the positivity of the unknown probability density function without affecting the high order accuracy of the base SL DG scheme. We analyze the stability and accuracy properties of the SL DG scheme by establishing its connection with the direct and weak formulations of the characteristics/Lagrangian Galerkin method [23]. The quality of the proposed method is demonstrated via basic test problems, such as linear advection and rigid body rotation, and via classical plasma problems, such as Landau damping and the two stream instability.     

On a new defect of shock-capturing methods

This paper proposes an explanation and a cure (or avoidance) to the new defect found of Eulerian shock-capturing methods in “A note on the conservative schemes for the Euler equations” by Tang and Liu [H. Tang, Tiegang Liu, A note on the conservative schemes for the Euler equations, J. Comput. Phys. 218 (2006) 451–459]. The latter gives a numerical investigation using several popular high resolution conservative schemes applied to Riemann problems of inviscid, compressible, perfect gas flows in Eulerian and Lagrangian coordinates with an initial high density ratio as well as a high pressure ratio. The results show that these methods work very inefficiently when applied to such problems and may give inaccurate numerical results, especially in shock location (or speed), even with a very fine grid.We have found that in problems of this type a strong rarefaction wave (SRW) is present adjacent to a contact line. Godunov averaging over the wave then produces large errors which, when the wave is strong, also persist for a long time. The cumulative error is thus very large which violates the strength of the contact line adjacent to it which, in turn, affects the speed and hence the location of the shock on the other side of the contact. We confirm this numerically using a method based on the unified coordinates with the shock-adaptive Godunov scheme plus contact strength preserving. The method, when applied to the Examples 2.1 and 2.2 of Tang and Liu [H. Tang, Tiegang Liu, A note on the conservative schemes for the Euler equations, J. Comput. Phys. 218 (2006) 451–459], produces high quality results even for comparatively coarse grids.     

Positivity-preserving high-resolution schemes for systems of conservation laws

A new class of flux-limited schemes for systems of conservation laws is presented that is both high-resolution and positivity-preserving. The schemes are obtained by extending the Steger–Warming method to second-order accuracy through the use of component-wise TVD flux limiters while ensuring that the coefficients of the discretization equation are positive. A coefficient is considered positive if it has all-positive eigenvalues and has the same eigenvectors as those of the convective flux Jacobian evaluated at the corresponding node. For certain systems of conservation laws, such as the Euler equations for instance, this condition is sufficient to guarantee positivity-preservation. The method proposed is advantaged over previous positivity-preserving flux-limited schemes by being capable to capture with high resolution all wave types (including contact discontinuities, shocks, and expansion fans). Several test cases are considered in which the Euler equations in generalized curvilinear coordinates are solved in 1D, 2D, and 3D. The test cases confirm that the proposed schemes are positivity-preserving while not being significantly more dissipative than the conventional TVD methods. The schemes are written in general matrix form and can be used to solve other systems of conservation laws, as long as they are homogeneous of degree one.     

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

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

Linear high-resolution schemes for hyperbolic conservation laws: TVB numerical evidence

The Osher–Chakrabarthy family of linear flux-modification schemes is considered. Improved lower bounds on the compression factors are provided, which suggest the viability of using the unlimited version. The LLF flux formula is combined with these schemes in order to obtain efficient finite-difference algorithms. The resulting schemes are applied to a battery of numerical tests, going from advection and Burgers equations to Euler and MHD equations, including the double Mach reflection and the Orszag–Tang 2D vortex problem. Total-variation-bounded (TVB) behavior is evident in all cases, even with time-independent upper bounds. The proposed schemes, however, do not deal properly with compound shocks, arising from non-convex fluxes, as shown by Buckley–Leverett test simulations.     

A class of hybrid DG/FV methods for conservation laws I: Basic formulation and one-dimensional systems

By comparing the discontinuous Galerkin (DG) and the finite volume (FV) methods, a concept of ‘static reconstruction’ and ‘dynamic reconstruction’ is introduced for high-order numerical methods. Based on the new concept, a class of hybrid DG/FV schemes is presented for one-dimensional conservation law using a ‘hybrid reconstruction’ approach. In the hybrid DG/FV schemes, the lower-order derivatives of a piecewise polynomial solution are computed locally in a cell by the DG method based on Taylor basis functions (called as ‘dynamic reconstruction’), while the higher-order derivatives are re-constructed by the ‘static reconstruction’ of the FV method, using the known lower-order derivatives in the cell itself and its adjacent neighboring cells. The hybrid DG/FV methods can greatly reduce CPU time and memory required by the traditional DG methods with the same order of accuracy on the same mesh, and they can be extended directly to unstructured and hybrid grids in two and three dimensions similar to the DG and/or FV methods. The hybrid DG/FV methods are applied to one-dimensional conservation law, including linear and non-linear scalar equation and Euler equations. In order to capture the strong shock waves without spurious oscillations, a simple shock detection approach is developed to mark ‘trouble cells’, and a moment limiter is adopted for higher-order schemes. The numerical results demonstrate the accuracy, and the super-convergence property is shown for the third-order hybrid DG/FV schemes. In addition, by analyzing the eigenvalues of the semi-discretized system in one dimension, we discuss the spectral properties of the hybrid DG/FV schemes to explain the super-convergence phenomenon.     

A comment on the computation of non-conservative products

We are interested in the solution of non-conservative hyperbolic systems, and consider in particular the so-called path-conservative schemes (see e.g. [2], [3]) which rely on the theoretical work in [1]. The example of the standard Euler equations for a perfect gas is used to illuminate some computational issues and shortcomings of this approach.     

Generalized conservative approximations of split convective derivative operators

We propose a strategy to design locally conservative finite-difference approximations of convective derivatives for shock-free compressible flows with arbitrary order of accuracy, that generalizes the approach of Ducros et al. (2000) [1], and that can be applied as a building block of low-dissipative, hybrid shock-capturing methods. The approximations stem from application of standard central difference formulas to split forms of the convective terms in the compressible Euler equations, which guarantee strong numerical stability and (near) energy preservation in the inviscid limit. A convenient implementation of the high-order fluxes is suggested, which guarantees improved computational efficiency over existing methods. Numerical tests performed for isotropic turbulence at zero viscosity show stability of schemes with order of accuracy up to ten, and effectiveness of convective splitting of Kennedy and Gruber (2008) [2] in providing extra stability in the presence of strong density variations. Numerical simulations of compressible turbulent boundary layer flow indicate suitability of the method for non-uniform grids, and overall support superior computational efficiency of high-order schemes.     

A hybrid,center-difference,limiter method for simulations of compressible multicomponent flows with Mie-Grüneisen equation of state

We develop an efficient spatially high-order, Cartesian-mesh, hybrid, center-difference, limiter methodology for numerical simulations of compressible multicomponent flows with isotropic Mie-Grüneisen equation of state. Effective switching between center-difference and limiter schemes is achieved by a set of robust tolerance and Lax-entropy based criterion [18]. Oscillations that could result from a mixed stencil scheme are minimized by requiring that the limiter method approaches the center-difference method in smooth regions. To achieve this the limiter is based on a norm of the deviation of WENO reconstruction weights from ideal. Results from a spatially 4th order version of the methodology are presented in one and two dimensions utilizing the California Institute of Technology’s VTF (Virtual Test Facility) AMROC [7] software.     

基于静动态重构的高阶DG/FV混合格式在二维非结构网格中的推广

通过比较间断Galerkin有限元方法(DGM)和有限体积方法(FVM),提出"静态重构"和"动态重构"的概念,进一步建立基于静动态"混合重构"算法的三阶DG/FV混合格式.在DG/FV混合格式中,单元平均值和一阶导数由DGM方法"动态重构",二阶导数利用FVM方法"静态重构";在此基础上,构造高阶多项式插值函数,得到...     

一种基于高阶节点间断Galerkin方法的新型限制器设计

使用高阶间断Galerkin(discontinuous Galerkin, DG)方法求解双曲守恒律方程组时, 非物理效应常常导致计算过程的中断, 这在很大程度上制约着该方法在计算流体力学中的应用.文章结合局部单元上原始流动变量的Taylor展开, 设计了一种新型的限制器, 通过对各阶空间导数的重构, 有效地消除了非物理振荡的不利影响.对二维Euler方程的计算结果表明, 该限制器不仅能够捕捉高质量的激波, 而且能够保证残值的有效收敛.      

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.     

A fully second order implicit/explicit time integration technique for hydrodynamics plus nonlinear heat conduction problems

We present a fully second order implicit/explicit time integration technique for solving hydrodynamics coupled with nonlinear heat conduction problems. The idea is to hybridize an implicit and an explicit discretization in such a way to achieve second order time convergent calculations. In this scope, the hydrodynamics equations are discretized explicitly making use of the capability of well-understood explicit schemes. On the other hand, the nonlinear heat conduction is solved implicitly. Such methods are often referred to as IMEX methods [2], [1], [3]. The Jacobian-Free Newton Krylov (JFNK) method (e.g. [10], [9]) is applied to the problem in such a way as to render a nonlinearly iterated IMEX method. We solve three test problems in order to validate the numerical order of the scheme. For each test, we established second order time convergence. We support these numerical results with a modified equation analysis (MEA) [21], [20]. The set of equations studied here constitute a base model for radiation hydrodynamics.     

求解Euter方程组的一种简单的守恒误差

本文提出了一种解一维不定常Euler方程组的较简单的一阶,二阶守恒差分格式,一阶格式简称为FUDE(First-Order UPwind Difference Scheme for Euler Equations),二阶格式简记为SUDE。本文的格式避免了大的数值振荡及膨胀激波,同时计算量较MacCormack^[1]格式仅略有增加,文中若干数值试验表明本文格式的分辨率是令人满意的。     

高精度TVD有限差分格式的自由梯度条件

1自由流条件守恒性气动方程组的高精度TVD有限差分格式一般通过使用通量限制器加入非线性的人工粘性来构成。在曲线坐标系中，通量限制器比较的两个因子与各自所在单元的界面雅可比转换行列式相关。为了消除畸变网格单元界面雅可比转换行列式的变化对通量限制器比较结果的影响，本文提出了高精度TVD格式的自由梯度性质，并分别就Osher和HartenTVD格式讨论了同时满足自由流和自由梯度的度量张量计算方法，使高精度TVD格式在畸变网格上得以实现。任意曲线坐标系中守恒型二维欧拉方程表示为其中，E、F为笛卡尔坐标系中的通量；J为坐标变…