High-Order Runge-Kutta Discontinuous Galerkin Methods with a New Type of Multi-Resolution WENO Limiters on Tetrahedral Meshes

In this paper, the second-order and third-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (WENO) limiters are proposed on tetrahedral meshes. The multi-resolution WENO limiter is an extension of a finite volume multi-resolution WENO scheme developed in [81], which serves as a limiter for RKDG methods on tetrahedral meshes. This new WENO limiter uses information of the DG solution essentially only within the troubled cell itself which is identified by a new modified version of the original KXRCF indicator [42], to build a sequence of hierarchical $L^2$ projection polynomials from zeroth degree to the second or third degree of the DG solution. The second-order and third-order RKDG methods with the associated multi-resolution WENO limiters are developed as examples for general high-order RKDG methods, which could maintain the original order of accuracy in smooth regions and keep essentially non-oscillatory property near strong discontinuities by gradually degrading from the optimal order to the first order. The linear weights inside the procedure of the new multi-resolution WENO limiters can be set as any positive numbers on the condition that they sum to one. A series of polynomials of different degrees within the troubled cell itself are applied in a WENO fashion to modify the DG solutions in the troubled cell on tetrahedral meshes. These new WENO limiters are very simple to construct, and can be easily implemented to arbitrary high-order accuracy on tetrahedral meshes. Such spatial reconstruction methodology improves the robustness in the simulation on the same compact spatial stencil of the original DG methods on tetrahedral meshes. Extensive one-dimensional (run as three-dimensional problems on tetrahedral meshes) and three-dimensional tests are performed to demonstrate the good performance of the RKDG methods with new multi-resolution WENO limiters.     

Adaptive Order WENO Reconstructions for the Semi-Lagrangian Finite Difference Scheme for Advection Problem

We present a new conservative semi-Lagrangian finite difference weighted essentially non-oscillatory scheme with adaptive order. This is an extension of the conservative semi-Lagrangian (SL) finite difference WENO scheme in [Qiu and Shu, JCP, 230 (4) (2011), pp. 863-889], in which linear weights in SL WENO framework were shown not to exist for variable coefficient problems. Hence, the order of accuracy is not optimal from reconstruction stencils. In this paper, we incorporate a recent WENO adaptive order (AO) technique [Balsara et al., JCP, 326 (2016), pp. 780-804] to the SL WENO framework. The new scheme can achieve an optimal high order of accuracy, while maintaining the properties of mass conservation and non-oscillatory capture of solutions from the original SL WENO. The positivity-preserving limiter is further applied to ensure the positivity of solutions. Finally, the scheme is applied to high dimensional problems by a fourth-order dimensional splitting. We demonstrate the effectiveness of the new scheme by extensive numerical tests on linear advection equations, the Vlasov-Poisson system, the guiding center Vlasov model as well as the incompressible Euler equations.     

Runge-Kutta Discontinuous Galerkin Method Using WENO-Type Limiters: Three-Dimensional Unstructured Meshes

This paper further considers weighted essentially non-oscillatory (WENO) and Hermite weighted essentially non-oscillatory (HWENO) finite volume methods as limiters for Runge-Kutta discontinuous Galerkin (RKDG) methods to solve problems involving nonlinear hyperbolic conservation laws. The application discussed here is the solution of 3-D problems on unstructured meshes. Our numerical tests again demonstrate this is a robust and high order limiting procedure, which simultaneously achieves high order accuracy and sharp non-oscillatory shock transitions.     

A Hermite WENO Method with Modified Ghost Fluid Method for Compressible Two-Medium Flow Problems

In this paper, we develop a novel approach by combining a new robust finite difference Hermite weighted essentially non-oscillatory (HWENO) method [51] with the modified ghost fluid method (MGFM) [25] to simulate the compressible two-medium flow problems. The main idea is that we first use the technique of the MGFM to transform a two-medium flow problem to two single-medium cases by defining the ghost fluids status based on the predicted interface status. Then the efficient and robust HWENO finite difference method is applied for solving the single-medium flow cases. By using immediate neighbor information to deal with both the solution and its derivatives, the fifth order finite difference HWENO scheme adopted in this paper is more compact and has higher resolution than the classical fifth order finite difference WENO scheme of Jiang and Shu [14]. Furthermore, by combining the HWENO scheme with the MGFM to simulate the two-medium flow problems, less ghost point information is needed than that in using the classical WENO scheme in order to obtain the same numerical accuracy. Various one-dimensional and two-dimensional two-medium flow problems are solved to illustrate the good performances of the proposed method.     

High Order Finite Difference WENO Methods with Unequal-Sized Sub-Stencils for the Degasperis-Procesi Type Equations

In this paper, we develop two finite difference weighted essentially non-oscillatory (WENO) schemes with unequal-sized sub-stencils for solving the Degasperis-Procesi (DP) and $\mu$-Degasperis-Procesi ($\mu$DP) equations, which contain nonlinear high order derivatives, and possibly peakon solutions or shock waves. By introducing auxiliary variable(s), we rewrite the DP equation as a hyperbolic-elliptic system, and the $\mu$DP equation as a first order system. Then we choose a linear finite difference scheme with suitable order of accuracy for the auxiliary variable(s), and two finite difference WENO schemes with unequal-sized sub-stencils for the primal variable. One WENO scheme uses one large stencil and several smaller stencils, and the other WENO scheme is based on the multi-resolution framework which uses a series of unequal-sized hierarchical central stencils. Comparing with the classical WENO scheme which uses several small stencils of the same size to make up a big stencil, both WENO schemes with unequal-sized sub-stencils are simple in the choice of the stencil and enjoy the freedom of arbitrary positive linear weights. Another advantage is that the final reconstructed polynomial on the target cell is a polynomial of the same degree as the polynomial over the big stencil, while the classical finite difference WENO reconstruction can only be obtained for specific points inside the target interval. Numerical tests are provided to demonstrate the high order accuracy and non-oscillatory properties of the proposed schemes.     

Moment-Based Multi-Resolution HWENO Scheme for Hyperbolic Conservation Laws

In this paper, a high-order moment-based multi-resolution Hermite weighted essentially non-oscillatory (HWENO) scheme is designed for hyperbolic conservation laws. The main idea of this scheme is derived from our previous work [J. Comput. Phys., 446 (2021) 110653], in which the integral averages of the function and its first order derivative are used to reconstruct both the function and its first order derivative values at the boundaries. However, in this paper, only the function values at the Gauss-Lobatto points in the one or two dimensional case need to be reconstructed by using the information of the zeroth and first order moments. In addition, an extra modification procedure is used to modify those first order moments in the troubled-cells, which leads to an improvement of stability and an enhancement of resolution near discontinuities. To obtain the same order of accuracy, the size of the stencil required by this moment-based multi-resolution HWENO scheme is still the same as the general HWENO scheme and is more compact than the general WENO scheme. Moreover, the linear weights are not unique and are independent of the node position, and the CFL number can still be 0.6 whether for the one or two dimensional case, which has to be 0.2 in the two dimensional case for other HWENO schemes. Extensive numerical examples are given to demonstrate the stability and resolution of such moment-based multi-resolution HWENO scheme.     

A Slope Constrained 4th Order Multi-Moment Finite Volume Method with WENO Limiter

This paper presents a new and better suited formulation to implement the limiting projection to high-order schemes that make use of high-order local reconstructions for hyperbolic conservation laws. The scheme, so-called MCV-WENO4 (multi-moment Constrained finite Volume with WENO limiter of 4th order) method, is an extension of the MCV method of Ii & Xiao (2009) by adding the 1st order derivative (gradient or slope) at the cell center as an additional constraint for the cell-wise local reconstruction. The gradient is computed from a limiting projection using the WENO (weighted essentially non-oscillatory) reconstruction that is built from the nodal values at 5 solution points within 3 neighboring cells. Different from other existing methods where only the cell-average value is used in the WENO reconstruction, the present method takes account of the solution structure within each mesh cell, and thus minimizes the stencil for reconstruction. The resulting scheme has 4th-order accuracy and is of significant advantage in algorithmic simplicity and computational efficiency. Numerical results of one and two dimensional benchmark tests for scalar and Euler conservation laws are shown to verify the accuracy and oscillation-less property of the scheme.     

Troubled-Cell Indicator Based on Mean Value Property for Hybrid WENO Schemes

In this paper, we introduce a new type of troubled-cell indicator to improve hybrid weighted essentially non-oscillatory (WENO) schemes for solving the hyperbolic conservation laws. The hybrid WENO schemes selectively adopt the high-order linear upwind scheme or the WENO scheme to avoid the local characteristic decompositions and calculations of the nonlinear weights in smooth regions. Therefore, they can reduce computational cost while maintaining non-oscillatory properties in non-smooth regions. Reliable troubled-cell indicators are essential for efficient hybrid WENO methods. Most of troubled-cell indicators require proper parameters to detect discontinuities precisely, but it is very difficult to determine the parameters automatically. We develop a new troubled-cell indicator derived from the mean value theorem that does not require any variable parameters. Additionally, we investigate the characteristics of indicator variable; one of the conserved properties or the entropy is considered as indicator variable. Detailed numerical tests for 1D and 2D Euler equations are conducted to demonstrate the performance of the proposed indicator. The results with the proposed troubled-cell indicator are in good agreement with pure WENO schemes. Also the new indicator has advantages in the computational cost compared with the other indicators.     

Finite Volume Hermite WENO Schemes for Solving the Hamilton-Jacobi Equation

In this paper, we present a new type of Hermite weighted essentially non-oscillatory (HWENO) schemes for solving the Hamilton-Jacobi equations on the finite volume framework. The cell averages of the function and its first one (in one dimension) or two (in two dimensions) derivative values are together evolved via time approaching and used in the reconstructions. And the major advantages of the new HWENO schemes are their compactness in the spacial field, purely on the finite volume framework and only one set of small stencils is used for different type of the polynomial reconstructions. Extensive numerical tests are performed to illustrate the capability of the methodologies.     

High Order Finite Difference Hermite WENO Fixed-Point Fast Sweeping Method for Static Hamilton-Jacobi Equations

In this paper, we combine the nonlinear HWENO reconstruction in [43] and the fixed-point iteration with Gauss-Seidel fast sweeping strategy, to solve the static Hamilton-Jacobi equations in a novel HWENO framework recently developed in [22]. The proposed HWENO frameworks enjoys several advantages. First, compared with the traditional HWENO framework, the proposed methods do not need to introduce additional auxiliary equations to update the derivatives of the unknown function $\phi$. They are now computed from the current value of $\phi$ and the previous spatial derivatives of $\phi$. This approach saves the computational storage and CPU time, which greatly improves the computational efficiency of the traditional HWENO scheme. In addition, compared with the traditional WENO method, reconstruction stencil of the HWENO methods becomes more compact, their boundary treatment is simpler, and the numerical errors are smaller on the same mesh. Second, the fixed-point fast sweeping method is used to update the numerical approximation. It is an explicit method and does not involve the inverse operation of nonlinear Hamiltonian, therefore any Hamilton-Jacobi equations with complex Hamiltonian can be solved easily. It also resolves some known issues, including that the iterative number is very sensitive to the parameter $ε$ used in the nonlinear weights, as observed in previous studies. Finally, to further reduce the computational cost, a hybrid strategy is also presented. Extensive numerical experiments are performed on two-dimensional problems, which demonstrate the good performance of the proposed fixed-point fast sweeping HWENO methods.     

A Robust WENO Type Finite Volume Solver for Steady Euler Equations on Unstructured Grids

A recent work of Li et al. [Numer. Math. Theor. Meth. Appl., 1(2008), pp. 92-112] proposed a finite volume solver to solve 2D steady Euler equations. Although the Venkatakrishnan limiter is used to prevent the non-physical oscillations nearby the shock region, the overshoot or undershoot phenomenon can still be observed. Moreover, the numerical accuracy is degraded by using Venkatakrishnan limiter. To fix the problems, in this paper the WENO type reconstruction is employed to gain both the accurate approximations in smooth region and non-oscillatory sharp profiles near the shock discontinuity. The numerical experiments will demonstrate the efficiency and robustness of the proposed numerical strategy.     

High-Order Local Discontinuous Galerkin Method with Multi-Resolution WENO Limiter for Navier-Stokes Equations on Triangular Meshes

In this paper, a new multi-resolution weighted essentially non-oscillatory (MR-WENO) limiter for high-order local discontinuous Galerkin (LDG) method is designed for solving Navier-Stokes equations on triangular meshes. This MR-WENO limiter is a new extension of the finite volume MR-WENO schemes. Such new limiter uses information of the LDG solution essentially only within the troubled cell itself, to build a sequence of hierarchical $L^2$ projection polynomials from zeroth degree to the highest degree of the LDG method. As an example, a third-order LDG method with associated same order MR-WENO limiter has been developed in this paper, which could maintain the original order of accuracy in smooth regions and could simultaneously suppress spurious oscillations near strong shocks or contact discontinuities. The linear weights of such new MR-WENO limiter can be any positive numbers on condition that their summation is one. This is the first time that a series of different degree polynomials within the troubled cell are applied in a WENO-type fashion to modify the freedom of degrees of the LDG solutions in the troubled cell. This MR-WENO limiter is very simple to construct, and can be easily implemented to arbitrary high-order accuracy and in higher dimensions on unstructured meshes. Such spatial reconstruction methodology improves the robustness in the numerical simulation on the same compact spatial stencil of the original LDG methods on triangular meshes. Some classicalviscous examples are given to show the good performance of this third-order LDG method with associated MR-WENO limiter.     

A High Order Sharp-Interface Method with Local Time Stepping for Compressible Multiphase Flows

In this paper, a new sharp-interface approach to simulate compressible multiphase flows is proposed. The new scheme consists of a high order WENO finite volume scheme for solving the Euler equations coupled with a high order path-conservative discontinuous Galerkin finite element scheme to evolve an indicator function that tracks the material interface. At the interface our method applies ghost cells to compute the numerical flux, as the ghost fluid method. However, unlike the original ghost fluid scheme of Fedkiw et al. [15], the state of the ghost fluid is derived from an approximate-state Riemann solver, similar to the approach proposed in [25], but based on a much simpler formulation. Our formulation leads only to one single scalar nonlinear algebraic equation that has to be solved at the interface, instead of the system used in [25]. Away from the interface, we use the new general Osher-type flux recently proposed by Dumbser and Toro [13], which is a simple but complete Riemann solver, applicable to general hyperbolic conservation laws. The time integration is performed using a fully-discrete one-step scheme, based on the approaches recently proposed in [5, 7]. This allows us to evolve the system also with time-accurate local time stepping. Due to the sub-cell resolution and the subsequent more restrictive time-step constraint of the DG scheme, a local evolution for the indicator function is applied, which is matched with the finite volume scheme for the solution of the Euler equations that runs with a larger time step. The use of a locally optimal time step avoids the introduction of excessive numerical diffusion in the finite volume scheme. Two different fluids have been used, namely an ideal gas and a weakly compressible fluid modeled by the Tait equation. Several tests have been computed to assess the accuracy and the performance of the new high order scheme. A verification of our algorithm has been carefully carried out using exact solutions as well as a comparison with other numerical reference solutions. The material interface is resolved sharply and accurately without spurious oscillations in the pressure field.     

A New Type of High-Order WENO Schemes for Hamilton-Jacobi Equations on Triangular Meshes

In this paper, a new type of third-order and fourth-order weighted essentially non-oscillatory (WENO) schemes is designed for simulating the Hamilton-Jacobi equations on triangular meshes. We design such schemes with the use of the nodal information defined on five unequal-sized spatial stencils, the application of monotone Hamiltonians as a building block, the artificial set of positive linear weights to make up high-order approximations in smooth regions simultaneously avoiding spurious oscillations nearby discontinuities of the derivatives of the solutions. The spatial reconstructions are convex combinations of the derivatives of a modified cubic/quartic polynomial defined on a big spatial stencil and four quadratic polynomials defined on small spatial stencils, and a third-order TVD Runge-Kutta method is used for the time discretization. The main advantages of these WENO schemes are their efficiency, simplicity, and can be easily implemented to higher dimensional unstructured meshes. Extensive numerical tests are performed to illustrate the good performance of such new WENO schemes.     

A New Approach of High Order Well-Balanced Finite Volume WENO Schemes and Discontinuous Galerkin Methods for a Class of Hyperbolic Systems with Source Terms

Hyperbolic balance laws have steady state solutions in which the flux gradients are nonzero but are exactly balanced by the source terms. In our earlier work [31–33], we designed high order well-balanced schemes to a class of hyperbolic systems with separable source terms. In this paper, we present a different approach to the same purpose: designing high order well-balanced finite volume weighted essentially non-oscillatory (WENO) schemes and RungeKutta discontinuous Galerkin (RKDG) finite element methods. We make the observation that the traditional RKDG methods are capable of maintaining certain steady states exactly, if a small modification on either the initial condition or the flux is provided. The computational cost to obtain such a well balanced RKDG method is basically the same as the traditional RKDG method. The same idea can be applied to the finite volume WENO schemes. We will first describe the algorithms and prove the well balanced property for the shallow water equations, and then show that the result can be generalized to a class of other balance laws. We perform extensive one and two dimensional simulations to verify the properties of these schemes such as the exact preservation of the balance laws for certain steady state solutions, the non-oscillatory property for general solutions with discontinuities, and the genuine high order accuracy in smooth regions.     

A Multi-Domain Hybrid DG and WENO Method for Hyperbolic Conservation Laws on Hybrid Meshes

In [SIAM J. Sci. Comput., 35(2)(2013), A1049–A1072], a class of multi-domain hybrid DG and WENO methods for conservation laws was introduced. Recent applications of this method showed that numerical instability may encounter if the DG flux with Lagrangian interpolation is applied as the interface flux during the moment of conservative coupling. In this continuation paper, we present a more robust approach in the construction of DG flux at the coupling interface by using WENO procedures of reconstruction. Based on this approach, such numerical instability is overcome very well. In addition, the procedure of coupling a DG method with a WENO-FD scheme on hybrid meshes is disclosed in detail. Typical testing cases are employed to demonstrate the accuracy of this approach and the stability under the flexibility of using either WENO-FD flux or DG flux at the moment of requiring conservative coupling.     

AUSM-Based High-Order Solution for Euler Equations

In this paper we demonstrate the accuracy and robustness of combining the advection upwind splitting method (AUSM), specifically AUSM⁺-UP [9], with high-order upwind-biased interpolation procedures, the weighted essentially non-oscillatory (WENO-JS) scheme [8] and its variations [2, 7], and the monotonicity preserving (MP) scheme [16], for solving the Euler equations. MP is found to be more effective than the three WENO variations studied. AUSM⁺-UP is also shown to be free of the so-called "carbuncle" phenomenon with the high-order interpolation. The characteristic variables are preferred for interpolation after comparing the results using primitive and conservative variables, even though they require additional matrix-vector operations. Results using the Roe flux with an entropy fix and the Lax-Friedrichs approximate Riemann solvers are also included for comparison. In addition, four reflective boundary condition implementations are compared for their effects on residual convergence and solution accuracy. Finally, a measure for quantifying the efficiency of obtaining high order solutions is proposed; the measure reveals that a maximum return is reached after which no improvement in accuracy is possible for a given grid size.     

Arbitrary-Lagrangian-Eulerian One-Step WENO Finite Volume Schemes on Unstructured Triangular Meshes

In this article we present a new class of high order accurate ArbitraryEulerian-Lagrangian (ALE) one-step WENO finite volume schemes for solving nonlinear hyperbolic systems of conservation laws on moving two dimensional unstructured triangular meshes. A WENO reconstruction algorithm is used to achieve high order accuracy in space and a high order one-step time discretization is achieved by using the local space-time Galerkin predictor proposed in [25]. For that purpose, a new element-local weak formulation of the governing PDE is adopted on moving space-time elements. The space-time basis and test functions are obtained considering Lagrange interpolation polynomials passing through a predefined set of nodes. Moreover, a polynomial mapping defined by the same local space-time basis functions as the weak solution of the PDE is used to map the moving physical space-time element onto a space-time reference element. To maintain algorithmic simplicity, the final ALE one-step finite volume scheme uses moving triangular meshes with straight edges. This is possible in the ALE framework, which allows a local mesh velocity that is different from the local fluid velocity. We present numerical convergence rates for the schemes presented in this paper up to sixth order of accuracy in space and time and show some classical numerical test problems for the two-dimensional Euler equations of compressible gas dynamics.     

Conservative Semi-Lagrangian Finite Difference WENO Formulations with Applications to the Vlasov Equation

In this paper, we propose a new conservative semi-Lagrangian (SL) finite difference (FD) WENO scheme for linear advection equations, which can serve as a base scheme for the Vlasov equation by Strang splitting [4]. The reconstruction procedure in the proposed SL FD scheme is the same as the one used in the SL finite volume (FV) WENO scheme [3]. However, instead of inputting cell averages and approximate the integral form of the equation in a FV scheme, we input point values and approximate the differential form of equation in a FD spirit, yet retaining very high order (fifth order in our experiment) spatial accuracy. The advantage of using point values, rather than cell averages, is to avoid the second order spatial error, due to the shearing in velocity (v) and electrical field (E) over a cell when performing the Strang splitting to the Vlasov equation. As a result, the proposed scheme has very high spatial accuracy, compared with second order spatial accuracy for Strang split SL FV scheme for solving the Vlasov-Poisson (VP) system. We perform numerical experiments on linear advection, rigid body rotation problem; and on the Landau damping and two-stream instabilities by solving the VP system. For comparison, we also apply (1) the conservative SL FD WENO scheme, proposed in [22] for incompressible advection problem, (2) the conservative SL FD WENO scheme proposed in [21] and (3) the non-conservative version of the SL FD WENO scheme in [3] to the same test problems. The performances of different schemes are compared by the error table, solution resolution of sharp interface, and by tracking the conservation of physical norms, energies and entropies, which should be physically preserved.     

High-Order Positivity-Preserving Well-Balanced Discontinuous Galerkin Methods for Euler Equations with Gravitation on Unstructured Meshes

In this paper, we propose a high-order accurate discontinuous Galerkin (DG) method for the compressible Euler equations under gravitational fields on unstructured meshes. The scheme preserves a general hydrostatic equilibrium state and provably guarantees the positivity of density and pressure at the same time. Comparing with the work on the well-balanced scheme for Euler equations with gravitation on rectangular meshes, the extension to triangular meshes is conceptually plausible but highly nontrivial. We first introduce a special way to recover the equilibrium state and then design a group of novel variables at the interface of two adjacent cells, which plays an important role in the well-balanced and positivity-preserving properties. One main challenge is that the well-balanced schemes may not have the weak positivity property. In order to achieve the well-balanced and positivity-preserving properties simultaneously while maintaining high-order accuracy, we carefully design DG spatial discretization with well-balanced numerical fluxes and suitable source term approximation. For the ideal gas, we prove that the resulting well-balanced scheme, coupled with strong stability preserving time discretizations, satisfies a weak positivity property. A simple existing limiter can be applied to enforce the positivity-preserving property, without losing high-order accuracy and conservation. Extensive one- and two-dimensional numerical examples demonstrate the desired properties of the proposed scheme, as well as its high resolution and robustness.