Search

The relativistic hydrodynamics (RHD) equations have three crucial intrinsic physical constraints on the primitive variables: positivity of pressure and density, and subluminal fluid velocity. However, numerical simulations can violate these constraints, leading to nonphysical results or even simulation failure. Designing genuinely physical-constraint-preserving (PCP) schemes is very difficult, as the primitive variables cannot be explicitly reformulated using conservative variables due to relativistic effects. In this paper, we propose three efficient Newton–Raphson (NR) methods for robustly recovering primitive variables from conservative variables. Importantly, we rigorously prove that these NR methods are always convergent and PCP, meaning they preserve the physical constraints throughout the NR iterations. The discovery of these robust NR methods and their PCP convergence analyses are highly nontrivial and technical. Our NR methods are versatile and can be seamlessly incorporated into any RHD schemes that require the recovery of primitive variables. As an application, we apply them to design PCP finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes for solving the RHD equations. Our PCP HWENO schemes incorporate high-order HWENO reconstruction, a PCP limiter, and strong-stability-preserving time discretization. We rigorously prove the PCP property of the fully discrete schemes using convex decomposition techniques. Moreover, we suggest the characteristic decomposition with rescaled eigenvectors and scale-invariant nonlinear weights to enhance the performance of the HWENO schemes in simulating large-scale RHD problems. Several demanding numerical tests are conducted to demonstrate the robustness, accuracy, and high resolution of the proposed PCP HWENO schemes and to validate the efficiency of our NR methods. • Propose three efficient Newton-Raphson (NR) methods for robustly recovering primitive variables. • Rigorously prove that these NR methods are always convergent. • Develop physical-constraint-preserving (PCP) Hermite WENO schemes for relativistic hydrodynamics. • Rigorously prove the PCP property of NR methods and fully discrete schemes. • Numerical examples demonstrate the robustness of PCP schemes and efficiency of NR methods. [ABSTRACT FROM AUTHOR]

In this paper, a new limiter using weighted essentially non-oscillatory (WENO) methodology is investigated for the Runge–Kutta discontinuous Galerkin (RKDG) methods for solving hyperbolic conservation laws. The idea is to use the high-order DG solution polynomial itself in the target cell and the linear polynomials which are reconstructed by the cell averages of solution in the target cell and its neighboring cells to reconstruct a new high-order polynomial in a manner of WENO methodology. Since only the linear polynomials need to be prepared for reconstruction, this limiter is very simple and compact with a stencil including only the target cell and its immediate neighboring cells. Numerical examples of various problems show that the new limiting procedure can simultaneously achieve uniform high-order accuracy and sharp, non-oscillatory shock transitions. [ABSTRACT FROM AUTHOR]

In this paper a third order finite volume weighted essentially non-oscillatory scheme is designed for solving hyperbolic conservation laws on tetrahedral meshes. Comparing with other finite volume WENO schemes designed on tetrahedral meshes, the crucial advantages of such new WENO scheme are its simplicity and compactness with the application of only six unequal size spatial stencils for reconstructing unequal degree polynomials in the WENO type spatial procedures, and easy choice of the positive linear weights without considering the topology of the meshes. The original innovation of such scheme is to use a quadratic polynomial defined on a big central spatial stencil for obtaining third order numerical approximation at any points inside the target tetrahedral cell in smooth region and switch to at least one of five linear polynomials defined on small biased/central spatial stencils for sustaining sharp shock transitions and keeping essentially non-oscillatory property simultaneously. By performing such new procedures in spatial reconstructions and adopting a third order TVD Runge–Kutta time discretization method for solving the ordinary differential equation (ODE), the new scheme's memory occupancy is decreased and the computing efficiency is increased. So it is suitable for large scale engineering requirements on tetrahedral meshes. Some numerical results are provided to illustrate the good performance of such scheme. [ABSTRACT FROM AUTHOR]

In this paper, a class of high-order central Hermite WENO (HWENO) schemes based on finite volume framework and staggered meshes is proposed for directly solving one- and two-dimensional Hamilton-Jacobi (HJ) equations. The methods involve the Lax-Wendroff type discretizations or the natural continuous extension of Runge-Kutta methods in time. This work can be regarded as an extension of central HWENO schemes for hyperbolic conservation laws (Tao et al. J. Comput. Phys. 318, 222-251, 2016) which combine the central scheme and the HWENO spatial reconstructions and therefore carry many features of both schemes. Generally, it is not straightforward to design a finite volume scheme to directly solve HJ equations and a key ingredient for directly solving such equations is the reconstruction of numerical Hamiltonians to guarantee the stability of methods. Benefited from the central strategy, our methods require no numerical Hamiltonians. Meanwhile, the zeroth-order and the first-order moments of the solution are involved in the spatial HWENO reconstructions which is more compact compared with WENO schemes. The reconstructions are implemented through a dimension-by-dimension strategy when the spatial dimension is higher than one. A collection of one- and two- dimensional numerical examples is performed to validate high resolution and robustness of the methods in approximating the solutions of HJ equations, which involve linear, nonlinear, smooth, non-smooth, convex or non-convex Hamiltonians. [ABSTRACT FROM AUTHOR]

In this continuing paper of (Zhu and Qiu, J Comput Phys 318 (2016), 110-121), a new fifth order finite difference weighted essentially non-oscillatory (WENO) scheme is designed to approximate the viscosity numerical solution of the Hamilton-Jacobi equations. This new WENO scheme uses the same numbers of spatial nodes as the classical fifth order WENO scheme which is proposed by Jiang and Peng (SIAM J Sci Comput 21 (2000), 2126-2143), and could get less absolute truncation errors and obtain the same order of accuracy in smooth region simultaneously avoiding spurious oscillations nearby discontinuities. Such new WENO scheme is a convex combination of a fourth degree accurate polynomial and two linear polynomials in aWENOtype fashion in the spatial reconstruction procedures. The linear weights of three polynomials are artificially set to be any random positive constants with a minor restriction and the new nonlinear weights are proposed for the sake of keeping the accuracy of the scheme in smooth region, avoiding spurious oscillations and keeping sharp discontinuous transitions in nonsmooth region simultaneously. The main advantages of such new WENO scheme comparing with the classical WENO scheme proposed by Jiang and Peng (SIAM J Sci Comput 21 (2000), 2126-2143) are its efficiency, robustness and easy implementation to higher dimensions. Extensive numerical tests are performed to illustrate the capability of the new fifth WENO scheme. [ABSTRACT FROM AUTHOR]

In this paper, we propose a high order conservative semi-Lagrangian (SL) finite difference Hermite weighted essentially non-oscillatory (HWENO) method for the Vlasov equation based on dimensional splitting. HWENO was first proposed for solving nonlinear hyperbolic problems by evolving both function values and its first derivative values (Qiu and Shu (2004) [23] ). The major advantage of HWENO, compared with the original WENO, lies in its compactness in reconstruction stencils. There are several new ingredients in this paper. Firstly we propose a mass-conservative SL HWENO scheme for a 1-D equation by working with a flux-difference form, following the work of Qiu and Christlieb (2010) [25] . Secondly, we propose a proper splitting for equations of partial derivatives in HWENO framework to ensure local mass conservation. The proposed fifth order SL HWENO scheme with dimensional splitting has been tested to work well in capturing filamentation structures without oscillations when the time step size is within the Eulerian CFL constraint. However, when the time stepping size becomes larger, numerical oscillations are observed for the ‘mass conservative’ dimensional splitting HWENO scheme, as there are extra source terms in equations of partial derivatives. In this case, we introduce WENO limiters to control oscillations. Classical numerical examples on linear passive transport problems, as well as the nonlinear Vlasov–Poisson system, have been tested to demonstrate the performance of the proposed scheme. [ABSTRACT FROM AUTHOR]

In this paper a new simple fifth order weighted essentially non-oscillatory (WENO) scheme is presented in the finite difference framework for solving the hyperbolic conservation laws. The new WENO scheme is a convex combination of a fourth degree polynomial with two linear polynomials in a traditional WENO fashion. This new fifth order WENO scheme uses the same five-point information as the classical fifth order WENO scheme [14,20] , could get less absolute truncation errors in L 1 and L ∞ norms, and obtain the same accuracy order in smooth region containing complicated numerical solution structures simultaneously escaping nonphysical oscillations adjacent strong shocks or contact discontinuities. The associated linear weights are artificially set to be any random positive numbers with the only requirement that their sum equals one. New nonlinear weights are proposed for the purpose of sustaining the optimal fifth order accuracy. The new WENO scheme has advantages over the classical WENO scheme [14,20] in its simplicity and easy extension to higher dimensions. Some benchmark numerical tests are performed to illustrate the capability of this new fifth order WENO scheme. [ABSTRACT FROM AUTHOR]

In this paper, we present a class of new Hermite weighted essentially non-oscillatory (HWENO) schemes based on finite volume framework to directly solve the Hamilton–Jacobi (HJ) equations. For HWENO reconstruction, both the cell average and the first moment of the solution are evolved, and for two dimensional case, HWENO reconstruction is based on a dimension-by-dimension strategy which is the first used in HWENO reconstruction. For spatial discretization, one of key points for directly solving HJ equation is the reconstruction of numerical fluxes. We follow the idea put forward by Cheng and Wang (2014) [3] to reconstruct the values of solution at Gauss–Lobatto quadrature points and numerical fluxes at the interfaces of cells, and for neither the convex nor concave Hamiltonian case, the monotone modification of numerical fluxes is added, which can guarantee the precision in the smooth region and converge to the entropy solution when derivative discontinuities come up. The third order TVD Runge–Kutta method is used for the time discretization. Extensive numerical experiments in one dimensional and two dimensional cases are performed to verify the efficiency of the methods. [ABSTRACT FROM AUTHOR]

In Zhu (2014) [1] , we presented a new type of the finite volume Hermite weighted essentially non-oscillatory (HWENO) schemes for solving the Hamilton–Jacobi equations on structured meshes, the key idea of the reconstruction in the HWENO schemes comes from the original WENO schemes, however the cell averages of the function and its two derivative values are together evolved via time approaching and used in the reconstructions, while only the function is evolved and used in the original WENO schemes which are nodal based approximations on finite difference framework. In this paper, we extend the method to solve the Hamilton–Jacobi equations on unstructured meshes. The major advantages of the new HWENO schemes on unstructured meshes are their compactness in the spatial field, purely on the finite volume framework and only one set of small stencils is needed for different types of the polynomial reconstruction. Extensive numerical tests are performed to illustrate the capability and high order of accuracy of these methodologies. [ABSTRACT FROM AUTHOR]

In { J. Comput. Phys. 229 (2010) 8105-8129}, we studied hybrid weighted essentially non-oscillatory (WENO) schemes with different indicators for hyperbolic conservation laws on uniform grids for Cartesian domains. In this paper, we extend the schemes to solve two-dimensional systems of hyperbolic conservation laws on curvilinear grids for non-Cartesian domains. Our goal is to obtain similar advantageous properties as those of the hybrid WENO schemes on uniform grids for Cartesian domains. Extensive numerical results strongly support that the hybrid WENO schemes with discontinuity indicators on curvilinear grids can also save considerably on computational cost in contrast to the pure WENO schemes. They also maintain the essentially non-oscillatory property for general solutions with discontinuities and keep the sharp shock transition. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, we extend a class of the Hermite weighted essentially non-oscillatory (HWENO) schemes for solving the Hamilton–Jacobi equations by Qiu and Shu (2005) [24] to two dimensional unstructured meshes. The idea of the reconstruction in the HWENO schemes comes from the original WENO schemes, however both the function and its first two derivative values are evolved via time advancing and used in the reconstructions, while only the function values are evolved and used in the original WENO schemes which are nodal based approximations. The third and fourth order HWENO schemes using the combinations of second order approximations with nonlinear weights and TVD Runge–Kutta time discretization method are used here. Comparing with the original WENO schemes for Hamilton–Jacobi equations, one major advantage of HWENO schemes presented here is its compactness in the reconstructions. Extensive numerical tests are performed to illustrate the capability and high order accuracy of the methodologies. [Copyright &y& Elsevier]

In Zhu and Qiu (J Comput Phys 228:6957-6976, ), we systematically investigated adaptive Runge-Kutta discontinuous Galerkin (RKDG) methods for hyperbolic conservation laws with different indicators, with an objective of obtaining efficient and reliable indicators to obtain better performance for adaptive computation to save computational cost. In this follow-up paper, we extend the method to solve two-dimensional problems. Although the main idea of the method for two-dimensional case is similar to that for one-dimensional case, the extension of the implementation of the method to two-dimensional case is nontrivial because of the complexity of the adaptive mesh with hanging nodes. We lay our emphasis on the implementation details including adaptive procedure, solution projection, solution reconstruction and troubled-cell indicator. Extensive numerical experiments are presented to show the effectiveness of the method. [ABSTRACT FROM AUTHOR]

Abstract: 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. [Copyright &y& Elsevier]

Abstract: In this paper, we continue on studying the Runge–Kutta discontinuous Galerkin (RKDG) methods to solve compressible multi-medium flow with conservative treatment of the moving material interface. Comparing with the paper by J. Qiu, T.G. Liu and B.C. Khoo [J. Comput. Phys. 222 (2007) 353–373], we adopt the HLLC flux instead of Lax-Friedrichs numerical flux, the finite volume weighted essentially nonoscillatory (WENO) and Hermite WENO (HWENO) reconstructions as limiter instead of TVB limiter for RKDG. The HLLC flux is based on the approximate Riemann solver with little numerical viscosity and can resolve the contact discontinuity and shear wave very well. For limiter procedure, first we use the KXRCF indicator to identify the troubled cell, then apply WENO or HWENO method to reconstruct the polynomial in the troubled cell, while maintaining the cell average. This limiter procedure is more accurate and less problem dependent than the TVB limiter. Numerical results in one dimension for multi-medium flows such as gas–gas and gas–water are provided to illustrate the capability of these procedures. [ABSTRACT FROM AUTHOR]

Abstract: 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. [ABSTRACT FROM AUTHOR]

Abstract: In this paper, we systematically investigate adaptive Runge–Kutta discontinuous Galerkin (RKDG) methods for hyperbolic conservation laws with different indicators which were based on the troubled cell indicators studied by Qiu and Shu [J. Qiu, C.-W. Shu, A comparison of troubled-cell indicators for Runge–Kutta discontinuous Galerkin mehtods using weighted essentially non-osillatory limiters, SIAM J. Sci. Comput. 27 (2005) 995–1013]. The emphasis is on comparison of the performance of adaptive RKDG method using different indicators, with an objective of obtaining efficient and reliable indicators to obtain better performance for adaptive computation to save computational cost. Both h-version and r-version adaptive methods are considered in the paper. The idea is to first identify “troubled cells” by different troubled-cell indicators, namely those cells where limiting might be needed and discontinuities might appear, then adopt an adaptive approach in these cells. A detailed numerical study in one-dimensional case is performed, addressing the issues of efficiency (less CPU cost and more accurate), non-oscillatory property, and resolution of discontinuities. [Copyright &y& Elsevier]

Abstract: In [J. Qiu, C.-W. Shu, Runge–Kutta discontinuous Galerkin method using WENO limiters, SIAM Journal on Scientific Computing 26 (2005) 907–929], Qiu and Shu investigated using weighted essentially non-oscillatory (WENO) finite volume methodology as limiters for the Runge–Kutta discontinuous Galerkin (RKDG) methods for solving nonlinear hyperbolic conservation law systems on structured meshes. In this continuation paper, we extend the method to solve two-dimensional problems on unstructured meshes, with the goal of obtaining a robust and high order limiting procedure to simultaneously obtain uniform high order accuracy and sharp, nonoscillatory shock transition for RKDG methods. Numerical results are provided to illustrate the behavior of this procedure. [Copyright &y& Elsevier]

Abstract: In this paper, a class of weighted essentially non-oscillatory (WENO) schemes with a Lax–Wendroff time discretization procedure, termed WENO-LW schemes, for solving Hamilton–Jacobi equations is presented. This is an alternative method for time discretization to the popular total variation diminishing (TVD) Runge–Kutta time discretizations. We explore the possibility in avoiding the nonlinear weights for part of the procedure, hence reducing the cost but still maintaining non-oscillatory properties for problems with strong discontinuous derivative. As a result, comparing with the original WENO with Runge–Kutta time discretizations schemes (WENO-RK) of Jiang and Peng [G. Jiang, D. Peng, Weighted ENO schemes for Hamilton–Jacobi equations, SIAM J. Sci. Comput. 21 (2000) 2126–2143] for Hamilton–Jacobi equations, the major advantages of WENO-LW schemes are more cost effective for certain problems and their compactness in the reconstruction. Extensive numerical experiments are performed to illustrate the capability of the method. [Copyright &y& Elsevier]

Abstract: The Runge–Kutta discontinuous Galerkin (RKDG) method for solving hyperbolic conservation laws is a high order finite element method, which utilizes the useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers, TVD Runge–Kutta time discretizations, and limiters. In this paper, we investigate using the RKDG finite element method for compressible two-medium flow simulation with conservative treatment of the moving material interfaces. Numerical results for both gas–gas and gas–water flows in one-dimension are provided to demonstrate the characteristic behavior of this approach. [Copyright &y& Elsevier]

Abstract: Runge–Kutta discontinuous Galerkin (RKDG) method is a high order finite element method for solving hyperbolic conservation laws employing useful features from high resolution finite volume schemes, such as the exact or approximate Riemann solvers serving as numerical fluxes, TVD Runge–Kutta time discretizations, and limiters. In most of the RKDG papers in the literature, the Lax–Friedrichs numerical flux is used due to its simplicity, although there are many other numerical fluxes which could also be used. In this paper, we systematically investigate the performance of the RKDG method based on different numerical fluxes, including the first-order monotone fluxes such as the Godunov flux, the Engquist–Osher flux, etc., and second-order TVD fluxes, with the objective of obtaining better performance by choosing suitable numerical fluxes. The detailed numerical study is mainly performed for the one dimensional system case, addressing the issues of CPU cost, accuracy, non-oscillatory property, and resolution of discontinuities. Numerical tests are also performed for two dimensional systems. [Copyright &y& Elsevier]

Abstract: In this paper we develop a Lax–Wendroff time discretization procedure for the discontinuous Galerkin method (LWDG) to solve hyperbolic conservation laws. This is an alternative method for time discretization to the popular total variation diminishing (TVD) Runge–Kutta time discretizations. The LWDG is a one step, explicit, high order finite element method. The limiter is performed once every time step. As a result, LWDG is more compact than Runge–Kutta discontinuous Galerkin (RKDG) and the Lax–Wendroff time discretization procedure is more cost effective than the Runge–Kutta time discretizations for certain problems including two-dimensional Euler systems of compressible gas dynamics when nonlinear limiters are applied. [Copyright &y& Elsevier]

Abstract: In this paper, a class of weighted essentially non-oscillatory (WENO) schemes based on Hermite polynomials, termed HWENO (Hermite WENO) schemes, for solving Hamilton–Jacobi equations is presented. The idea of the reconstruction in the HWENO schemes comes from the original WENO schemes, however both the function and its first derivative values are evolved in time and used in the reconstruction, while only the function values are evolved and used in the original WENO schemes. Comparing with the original WENO schemes of Jiang and Peng [Weighted ENO schemes for Hamilton–Jacobi equations, SIAM Journal on Scientific Computing 21 (2000) 2126] for Hamilton–Jacobi equations, one major advantage of HWENO schemes is its compactness in the reconstruction. Extensive numerical experiments are performed to illustrate the capability of the method. [Copyright &y& Elsevier]

In this paper, a class of fifth-order weighted essentially non-oscillatory (WENO) schemes based on Hermite polynomials, termed HWENO (Hermite WENO) schemes, for solving one-dimensional nonlinear hyperbolic conservation law systems is presented. The construction of HWENO schemes is based on a finite volume formulation, Hermite interpolation, and nonlinearly stable Runge–Kutta methods. The idea of the reconstruction in the HWENO schemes comes from the original WENO schemes, however both the function and its first derivative values are evolved in time and used in the reconstruction, while only the function values are evolved and used in the original WENO schemes. Comparing with the original WENO schemes of Liu et al. [J. Comput. Phys. 115 (1994) 200] and Jiang and Shu [J. Comput. Phys. 126 (1996) 202], one major advantage of HWENO schemes is its compactness in the reconstruction. For example, five points are needed in the stencil for a fifth-order WENO (WENO5) reconstruction, while only three points are needed for a fifth-order HWENO (HWENO5) reconstruction. For this reason, the HWENO finite volume methodology is more suitable to serve as limiters for the Runge–Kutta discontinuous Galerkin (RKDG) methods, than the original WENO finite volume methodology. Such applications in one space dimension is also developed in this paper. [Copyright &y& Elsevier]

In this paper, we review and construct fifth- and ninth-order central weighted essentially nonoscillatory (WENO) schemes based on a finite volume formulation, staggered mesh, and continuous extension of Runge–Kutta methods for solving nonlinear hyperbolic conservation law systems. Negative linear weights appear in such a formulation and they are treated using the technique recently introduced by Shi et al. (J. Comput. Phys.175, 108 (2002)). We then perform numerical computations and comparisons with the finite difference WENO schemes of Jiang and Shu (J. Comput. Phys.150, 97 (1999)) and Balsara and Shu (J. Comput. Phys.160, 405 (2000)). The emphasis is on the performance with or without a local characteristic decomposition. While this decomposition increases the computational cost, we demonstrate by our numerical experiments that it is still necessary to use it to control spurious oscillations when the order of accuracy is high, both for the central staggered grid and for the upwind nonstaggered grid WENO schemes. We use the shock entropy wave interaction problem to demonstrate the advantage of using higher order WENO schemes when both shocks and complex solution features coexist. [Copyright &y& Elsevier]

• A hybrid finite volume Hermite WENO method with artificial linear weights is presented for conservation laws. • The method is the fifth order accuracy in both one and two dimension. • The new method is higher efficiency than the existing HWENO method. In this paper, a fifth-order Hermite weighted essentially non-oscillatory (HWENO) scheme with artificial linear weights is proposed for one and two dimensional hyperbolic conservation laws, where the zeroth-order and the first-order moments are used in the spatial reconstruction. We construct the HWENO methodology using a nonlinear convex combination of a high degree polynomial with several low degree polynomials, and the associated linear weights can be any artificial positive numbers with only requirement that their summation equals one. The one advantage of the HWENO scheme is its simplicity and easy extension to multi-dimension in engineering applications for we can use any artificial linear weights which are independent on geometry of mesh. The another advantage is its higher order numerical accuracy using less candidate stencils for two dimensional problems. In addition, the HWENO scheme still keeps the compactness as only immediate neighbor information is needed in the reconstruction and has high efficiency for directly using linear approximation in the smooth regions. In order to avoid nonphysical oscillations nearby strong shocks or contact discontinuities, we adopt the thought of limiter for discontinuous Galerkin method to control the spurious oscillations. Some benchmark numerical tests are performed to demonstrate the capability of the proposed scheme. [ABSTRACT FROM AUTHOR]

• A new type of multi-resolution WENO limiter has been designed for Runge-Kutta discontinuous Galerkine (RKDG) schemes. • The limiter uses the information of the DG solution essentially only within the troubled cell itself, hence is very local. • Numerical examples are provided to demonstrate the good performance of this new limiter. In this paper, a new type of multi-resolution weighted essentially non-oscillatory (WENO) limiters for high-order Runge-Kutta discontinuous Galerkin (RKDG) methods is designed. This type of multi-resolution WENO limiters is an extension of the multi-resolution WENO finite volume and finite difference schemes developed in [43]. Such new limiters use information of the DG 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 RKDG method. The second-order, third-order, fourth-order, and fifth-order RKDG methods with these multi-resolution WENO limiters have been developed as examples, which could maintain the original order of accuracy in smooth regions and could simultaneously suppress spurious oscillations near strong discontinuities. The linear weights of such new multi-resolution WENO limiters can be any positive numbers on the condition that their sum equals one. This is the first time that 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. These new WENO limiters are very simple to construct, and can be easily implemented to arbitrary high-order accuracy and in higher dimensions. Such spatial reconstruction methodology improves the robustness in the numerical simulation on the same compact spatial stencil of the original DG methods. Benchmark examples are given to demonstrate the good performance of these RKDG methods with the associated multi-resolution WENO limiters. [ABSTRACT FROM AUTHOR]

In this paper, a compact and high order ADER (Arbitrary high order using DERivatives) scheme using the simple HWENO method (ADER-SHWENO) is proposed for hyperbolic conservation laws. The newly-developed method employs the Lax-Wendroff procedure to convert time derivatives to spatial derivatives, which provides the time evolution of the variables at the cell interfaces. This information is required for the simple HWENO reconstructions, which take advantages of the simple WENO and the classic HWENO. Compared with the original Runge-Kutta HWENO method (RK-HWENO), the new method has two advantages. Firstly, RK-HWENO method must solve the additional equations for reconstructions, which is avoided for the new method. Secondly, the SHWENO reconstruction is performed once with one stencil and is different from the classic HWENO methods, in which both the function and its derivative values are reconstructed with two different stencils, respectively. Thus the new method is more efficient than the RK-HWENO method. Moreover, the new method is more compact than the existing ADER-WENO method. Besides, the new method makes the best use of the information in the ADER method. Thus, the time evolution of the cell averages of the derivatives is simpler than that developed in the work (Li et al. (2021) [17]). Numerical tests indicate that the new method can achieve high order for smooth solutions both in space and time, keep non-oscillatory at discontinuities. • A compact fifth-order ADER scheme is proposed for hyperbolic conservation laws. • A compact reconstruction technique is designed by better utilizing the information in time evolution. • The computational stencil is more compact than the original ADERWENO scheme. • It is more efficient than the original multi-stage RK-HWENO scheme. [ABSTRACT FROM AUTHOR]

Since the classical WENO schemes [27] might suffer from slight post-shock oscillations (which are responsible for the numerical residual to hang at a truncation error level) and the new high-order multi-resolution WENO schemes [59] are successful to solve for steady-state problems, we apply these high-order finite volume multi-resolution WENO techniques to serve as limiters for high-order Runge-Kutta discontinuous Galerkin (RKDG) methods in simulating steady-state problems. Firstly, a new troubled cell indicator is designed to precisely detect the cells which would need further limiting procedures. Then the high-order multi-resolution WENO limiting procedures are adopted on a sequence of hierarchical L 2 projection polynomials of the DG solution within the troubled cell itself. By doing so, these RKDG methods with multi-resolution WENO limiters could gradually degrade from the optimal high-order accuracy to the first-order accuracy near strong discontinuities, suppress the slight post-shock oscillations, and push the numerical residual to settle down to machine zero in steady-state simulations. These new multi-resolution WENO limiters are very simple to construct and can be easily implemented to arbitrary high-order accuracy for solving steady-state problems in multi-dimensions. [ABSTRACT FROM AUTHOR]

In this paper, high-order Runge-Kutta discontinuous Galerkin (RKDG) methods with multi-resolution weighted essentially non-oscillatory (WENO) limiters are designed for solving hyperbolic conservation laws on triangular meshes. These multi-resolution WENO limiters are new extensions of the associated multi-resolution WENO finite volume schemes [49,50] which serve as limiters for RKDG methods from structured meshes [47] to triangular meshes. Such new WENO limiters use 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 [24] , to build a sequence of hierarchical L 2 projection polynomials from zeroth degree to the highest degree of the RKDG method. The second-order, third-order, and fourth-order RKDG methods with associated multi-resolution WENO limiters are developed as examples, which could maintain the original order of accuracy in smooth regions and keep essentially non-oscillatory property near strong shocks or contact discontinuities by gradually degrading from the highest order to the first order. The linear weights inside the procedure of the new multi-resolution WENO limiters can be any positive numbers on the condition that their sum equals one. This is the first time that 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 triangular meshes. These new WENO limiters are 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 simulation on the same compact spatial stencil of the original DG methods on triangular meshes. Extensive one-dimensional (run as two-dimensional problems on triangular meshes) and two-dimensional tests are performed to demonstrate the effectiveness of these RKDG methods with the new multi-resolution WENO limiters. [ABSTRACT FROM AUTHOR]

A moving mesh finite difference method based on the moving mesh partial differential equation is proposed for the numerical solution of the 2T model for multi-material, non-equilibrium radiation diffusion equations. The model involves nonlinear diffusion coefficients and its solutions stay positive for all time when they are positive initially. Nonlinear diffusion and preservation of solution positivity pose challenges in the numerical solution of the model. A coefficient-freezing predictor-corrector method is used for nonlinear diffusion while a cutoff strategy with a positive threshold is used to keep the solutions positive. Furthermore, a two-level moving mesh strategy and a sparse matrix solver are used to improve the efficiency of the computation. Numerical results for a selection of examples of multi-material non-equilibrium radiation diffusion show that the method is capable of capturing the profiles and local structures of Marshak waves with adequate mesh concentration. The obtained numerical solutions are in good agreement with those in the existing literature. Comparison studies are also made between uniform and adaptive moving meshes and between one-level and two-level moving meshes. [ABSTRACT FROM AUTHOR]

• A moving mesh discontinuous Galerkin method is presented for conservation laws. • The method is a combination of the DG method and the moving mesh PDE. • No interpolation is needed for physical variables from the old mesh to the new one. • The method achieves the theoretically predicted order of convergence for problems with smooth solutions and is able to capture shocks and concentrate mesh points in non-smooth regions. A moving mesh discontinuous Galerkin method is presented for the numerical solution of hyperbolic conservation laws. The method is a combination of the discontinuous Galerkin method and the mesh movement strategy which is based on the moving mesh partial differential equation approach and moves the mesh continuously in time and orderly in space. It discretizes hyperbolic conservation laws on moving meshes in the quasi-Lagrangian fashion with which the mesh movement is treated continuously and no interpolation is needed for physical variables from the old mesh to the new one. Two convection terms are induced by the mesh movement and their discretization is incorporated naturally in the DG formulation. Numerical results for a selection of one- and two-dimensional scalar and system conservation laws are presented. It is shown that the moving mesh DG method achieves the second and third order of convergence for P 1 and P 2 elements, respectively, for problems with smooth solutions and is able to capture shocks and concentrate mesh points in non-smooth regions. Its advantage over uniform meshes and its insensitiveness to mesh smoothness are also demonstrated. [ABSTRACT FROM AUTHOR]

This paper studies three high-order structure-preserving finite volume weighted essentially non-oscillatory (WENO) methods, which are not only well balanced (WB) for a general known hydrostatic equilibrium state but also preserve the positivity of density and pressure, for the compressible Euler equations under gravitational fields. These methods are built on a simple local scaling positivity-preserving (PP) limiter and a modified WENO-ZQ reconstruction exactly preserving the cell average value and scaling invariance. The WB properties of these three methods are achieved based on suitable numerical fluxes and approximation to the gravitational source terms. Based on some convex decomposition techniques as well as several critical properties of the admissible states and numerical flux, we carry out rigorous positivity-preserving analyses for these three WB schemes. We rigorously prove that the three WB methods, coupled with the PP limiter and a strong-stability-preserving time discretization, are always PP under suitable Courant-Friedrichs-Lewy conditions. Extensive numerical examples are provided to confirm WB and PP properties of three methods. [ABSTRACT FROM AUTHOR]

We propose a one-stage high-order gas-kinetic scheme (GKS) with interface-sharpening technique for the 4-equations model of compressible multi-component stiffened-gas flows. The scheme is constructed based on the framework of a one-stage efficient high-order GKS (EHGKS), which can achieve uniform high-order accuracy in both space and time. The main idea in the construction of the new scheme consists of three parts. Firstly, different from the original EHGKS, the non-oscillatory kinetic (NOK) flux is introduced to compute the leading term and the related flux evaluations, which enables the modified EHGKS suitable for stiffened gases, together with no pressure and velocity oscillations across the contact discontinuities. Then the modified scheme is utilized to discrete the conservative parts of the 4-equations model. Secondly, in order to avoid the oscillations across material interfaces, a high-order GKS is derived to solve the non-conservative parts of the 4-equations model following the idea of Abgrall [1]. Finally, to further reduce the over diffusion near a material interface, a simple interface-sharpening technique is developed, which combines the idea of interface construction and downwind scheme. The numerical results demonstrate that the proposed scheme can not only achieve the designed uniform accuracy in both space and time, but also sharply capture material interfaces. • A one-stage high-order GKS for the compressible multicomponent flows is developed. • The non-oscillatory kinetic flux is used to avoid the interface oscillation. • An interface-sharpening technique is developed to capture the interfaces sharply. • Numerical tests validate the good performance of the new scheme. [ABSTRACT FROM AUTHOR]

• A moving mesh DG method is proposed for the Ripa model. • The Ripa model has an extra relation η = C 1 h in the lake-at-rest steady state. • The new method preserves the lake-at-rest steady state and solution positivity. • The method is able to capture perturbations of the lake-at-rest steady state. A well-balanced moving mesh discontinuous Galerkin (DG) method is proposed for the numerical solution of the Ripa model – a generalization of the shallow water equations that accounts for effects of water temperature variations. Thermodynamic processes are important particularly in the upper layers of the ocean where the variations of sea surface temperature play a fundamental role in climate change. The well-balance property which requires numerical schemes to preserve the lake-at-rest steady-state is crucial to the simulation of perturbation waves over that steady state such as waves on a lake or tsunami waves in the deep ocean. To ensure the well-balance, positivity-preserving, and high-order properties, a DG-interpolation scheme (with or without scaling positivity-preserving limiter) and special treatments pertaining to the Ripa model are employed in the transfer of both the flow variables and bottom topography from the old mesh to the new one and in the TVB limiting process. Mesh adaptivity is realized using an MMPDE moving mesh approach and a metric tensor based on an equilibrium variable and water depth. A motivation is to adapt the mesh according to both the perturbations of the lake-at-rest steady-state and the water depth distribution (bottom topography structure). Numerical examples in one and two dimensions are presented to demonstrate the well-balance, high-order accuracy, and positivity-preserving properties of the method and its ability to capture small perturbations of the lake-at-rest steady-state. [ABSTRACT FROM AUTHOR]

In this paper, a new type of finite difference Hermite weighted essentially non-oscillatory (HWENO) schemes are constructed for solving Hamilton–Jacobi (HJ) equations. Point values of both the solution and its first derivatives are used in the HWENO reconstruction and evolved via time advancing. While the evolution of the solution is still through the classical numerical fluxes to ensure convergence to weak solutions, the evolution of the first derivatives of the solution is through a simple dimension-by-dimension non-conservative procedure to gain efficiency. The main advantages of this new scheme include its compactness in the spatial field and its simplicity in the reconstructions. Extensive numerical experiments in one and two dimensional cases are performed to verify the accuracy, high resolution and efficiency of this new scheme. [ABSTRACT FROM AUTHOR]

In this paper, a class of high-order central finite volume schemes is proposed for solving one- and two-dimensional hyperbolic conservation laws. Formulated on staggered meshes, the methods involve Hermite WENO (HWENO) spatial reconstructions, and Lax–Wendroff type discretizations or the natural continuous extension of Runge–Kutta methods in time. Differently from the central Hermite WENO methods we developed previously in Tao et al. (2015) [34] , the spatial reconstructions, a core ingredient of the methods, are based on the zeroth-order and the first-order moments of the solution, and are implemented through a dimension-by-dimension strategy when the spatial dimension is higher than one. This leads to much simpler implementation of the methods in higher dimension and better cost efficiency. Meanwhile, the proposed methods have the attractive features of the general central Hermite WENO methods such as being compact in reconstruction and requiring neither flux splitting nor numerical fluxes, while being accurate and essentially non-oscillatory. A collection of one- and two-dimensional numerical examples is presented to demonstrate high resolution and robustness of the methods in capturing smooth and non-smooth solutions. [ABSTRACT FROM AUTHOR]

A hybrid LDG-HWENO scheme is proposed for the numerical solution of KdV-type partial differential equations. It evolves the cell averages of the physical solution and its moments (a feature of Hermite WENO) while discretizes high order spatial derivatives using the local DG method. The new scheme has the advantages of both LDG and HWENO methods, including the ability to deal with high order spatial derivatives and the use of a small number of global unknown variables. The latter is independent of the order of the scheme and the spatial order of the underlying differential equations. One and two dimensional numerical examples are presented to show that the scheme can attain the same formal high order accuracy as the LDG method. [ABSTRACT FROM AUTHOR]

An efficient moving mesh finite difference method is developed for the numerical solution of equilibrium radiation diffusion equations in two dimensions. The method is based on the moving mesh partial differential equation approach and moves the mesh continuously in time using a system of meshing partial differential equations. The mesh adaptation is controlled through a Hessian-based monitor function and the so-called equidistribution and alignment principles. Several challenging issues in the numerical solution are addressed. Particularly, the radiation diffusion coefficient depends on the energy density highly nonlinearly. This nonlinearity is treated using a predictor–corrector and lagged diffusion strategy. Moreover, the nonnegativity of the energy density is maintained using a cutoff method which has been known in literature to retain the accuracy and convergence order of finite difference approximation for parabolic equations. Numerical examples with multi-material, multiple spot concentration situations are presented. Numerical results show that the method works well for radiation diffusion equations and can produce numerical solutions of good accuracy. It is also shown that a two-level mesh movement strategy can significantly improve the efficiency of the computation. [ABSTRACT FROM AUTHOR]

In this paper, we construct a high order weighted essentially non-oscillatory (WENO) finite difference discretization for compressible Navier-Stokes (NS) equations, which is rendered positivity-preserving of density and internal energy by a positivity-preserving flux splitting and a scaling positivity-preserving limiter. The novelty of this paper is WENO reconstruction performed on variables from a positivity-preserving convection diffusion flux splitting, which is different from conventional WENO schemes solving compressible NS equations. The core advantages of our proposed method are robustness and efficiency, which especially are suitable for solving tough demanding problems of both compressible Euler and NS equation including low density and low pressure flow regime. Moreover, in terms of computational cost, it is more efficient and easier to implement and extend to multi-dimensional problems than the positivity-preserving high order discontinuous Galerkin schemes and finite volume WENO scheme for solving compressible NS equations on rectangle domain. Benchmark tests demonstrate that the proposed positivity-preserving WENO schemes are high order accuracy, efficient and robust without excessive artificial viscosity for demanding problems involving with low density, low pressure, and fine structure. • A positivity-preserving high order accurate finite difference WENO scheme for compressible N-S equations. • The method is the high order accuracy and positivity-preserving. • WENO reconstruction is performed on variables from a positivity-preserving convection diffusion flux splitting. • The method is robust and suitable for solving the compressible flow with low density and low pressure. [ABSTRACT FROM AUTHOR]

Abstract: Ideal MHD equations arise in many applications such as astrophysical plasmas and space physics, and they consist of a system of nonlinear hyperbolic conservation laws. The exact density ρ and pressure p should be non-negative. Numerically, such positivity property is not always satisfied by approximated solutions. One can encounter this when simulating problems with low density, high Mach number, or much large magnetic energy compared with internal energy. When this occurs, numerical instability may develop and the simulation can break down. In this paper, we propose positivity-preserving discontinuous Galerkin and central discontinuous Galerkin methods for solving ideal MHD equations by following [X. Zhang, C.-W. Shu, Journal of Computational Physics 229 (2010) 8918–8934]. In one dimension, the positivity-preserving property is established for both methods under a reasonable assumption. The performance of the proposed methods, in terms of accuracy, stability and positivity-preserving property, is demonstrated through a set of one and two dimensional numerical experiments. The proposed methods formally can be of any order of accuracy. [Copyright &y& Elsevier]

Abstract: 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. [Copyright &y& Elsevier]

• A compact fifth-order gas-kinetic scheme is proposed for compressible Euler equations. • It is based on the framework of the one-stage efficient high-order gas-kinetic scheme. • A newly-developed simple and efficient compact reconstruction technique is used. • It is more efficient than the original one-stage EHGKS with WENO reconstruction. In this paper, a compact fifth-order gas-kinetic scheme (GKS) is proposed to solve the Euler equations of compressible flows. It is based on the framework of a one-stage efficient high-order GKS (EHGKS), which can achieve arbitrary high-order accuracy in both space and time, and the newly-developed compact reconstruction technique, combining both the ideas of the simple weighted essentially non-oscillatory scheme (WENO) and Hermite WENO (HWENO). The evolution model of EHGKS provides the time-dependent variables at cell interfaces, which can directly compute the cell-average slope. Thus the necessary information for a fifth-order compact HWENO reconstruction technique is obtained, and also making it possible to improve simple WENO as a compact reconstruction by determining the quartic polynomial based on the compact stencil. The one-stage time discretization technique of EHGKS makes it possible to utilize the reconstruction only once in each computing time step, which is very different from the present compact multi-stage GKS schemes. The newly-developed compact reconstruction provides the polynomial, not only the pointwise variables but also their spatial derivatives for the interface variables evaluations directly. Several numerical examples are given to validate the accuracy and effectiveness of the proposed compact one-stage scheme. And the numerical results also demonstrate that the new scheme is more efficient than the original EHGKS method with WENO reconstruction technique. [ABSTRACT FROM AUTHOR]

In this paper, a new type of high-order finite volume and finite difference multi-resolution Hermite weighted essentially non-oscillatory (HWENO) schemes are designed for solving hyperbolic conservation laws on structured meshes. Here we only use the information defined on a hierarchy of nested central spatial stencils but do not introduce any equivalent multi-resolution representation, the terminology of multi-resolution HWENO follows that of the multi-resolution WENO schemes (Zhu and Shu, 2018) [29]. The main idea of our spatial reconstruction is derived from the original HWENO schemes (Qiu and Shu, 2004) [19] , in which both the function and its first-order derivative values are evolved in time and used in the reconstruction. Our HWENO schemes use the same large stencils as the classical HWENO schemes which are narrower than the stencils of the classical WENO schemes for the same order of accuracy. Only the function values need to be reconstructed by our HWENO schemes, the first-order derivative values are obtained from the high-order linear polynomials directly. Furthermore, the linear weights of such HWENO schemes can be any positive numbers as long as their sum equals one, and there is no need to do any modification or positivity-preserving flux limiting in our numerical experiments. Extensive benchmark examples are performed to illustrate the robustness and good performance of such finite volume and finite difference HWENO schemes. • High-order multi-resolution Hermite weighted essentially non-oscillatory (HWENO) schemes are designed. • Both finite volume and finite difference schemes are presented. • The new methods are robust, there are no need any additional limiter. [ABSTRACT FROM AUTHOR]

• A high order finite difference scheme for compressible two-medium flows. • The scheme is conservative. • The new method can maintain high order accuracy when the solution is smooth. • WENO is applied to the primitive variables to maintain the equilibrium of velocity and pressure across the interface. In this paper, a high order finite difference conservative scheme is proposed to solve two-medium flows. Our scheme has four advantages: First, our scheme is conservative, which is important to ensure the numerical solution captures the main features properly. Second, our scheme directly applies the WENO interpolation method to the primitive variables so that it can maintain the equilibrium of velocity and pressure across the interface, which is very helpful to obtain a non-oscillatory solution. Third, the usage of nodal values enables us to manipulate algebraic functions easily. Fourth, the scheme can maintain high order accuracy when the solution is smooth. Extensive numerical experiments are performed to verify the high resolution and non-oscillatory performance of this new scheme. [ABSTRACT FROM AUTHOR]