Search

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]

Summary: A robust, adaptive unstructured mesh refinement strategy for high‐order Runge‐Kutta discontinuous Galerkin method is proposed. The present work mainly focuses on accurate capturing of sharp gradient flow features like strong shocks in the simulations of two‐dimensional inviscid compressible flows. A posteriori finite volume subcell limiter is employed in the shock‐affected cells to control numerical spurious oscillations. An efficient cell‐by‐cell adaptive mesh refinement is implemented to increase the resolution of our simulations. This strategy enables to capture strong shocks without much numerical dissipation. A wide range of challenging test cases is considered to demonstrate the efficiency of the present adaptive numerical strategy for solving inviscid compressible flow problems having strong shocks. [ABSTRACT FROM AUTHOR]

SUMMARY A Cartesian grid method using immersed boundary technique to simulate the impact of body in fluid has become an important research topic in computational fluid dynamics because of its simplification, automation of grid generation, and accuracy of results. In the frame of Cartesian grid, one often uses finite volume method with second order accuracy or finite difference method. In this paper, an h-adaptive Runge-Kutta discontinuous Galerkin (RKDG) method on Cartesian grid with ghost cell immersed boundary method for arbitrarily complex geometries is developed. A ghost cell immersed boundary treatment with the modification of normal velocity is presented. The method is validated versus well documented test problems involving both steady and unsteady compressible flows through complex bodies over a wide range of Mach numbers. The numerical results show that the present boundary treatment to some extent reduces the error of entropy and demonstrate the efficiency, robustness, and versatility of the proposed approach. Copyright © 2013 John Wiley & Sons, Ltd. [ABSTRACT FROM AUTHOR]

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: 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]

• 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]

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 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]

• 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]

• A hybrid finite volume Hermite WENO method is presented for conservation laws. • The method is a combination of the DG framework with finite volume method. • The new method is higher efficiency than the existed HWENO method. In this paper, we propose a hybrid finite volume Hermite weighted essentially non-oscillatory (HWENO) scheme for solving one and two dimensional hyperbolic conservation laws, which would be the fifth order accuracy in the one dimensional case, while is the fourth order accuracy for two dimensional problems. The zeroth-order and the first-order moments are used in the spatial reconstruction, with total variation diminishing Runge-Kutta time discretization. Unlike the original HWENO schemes [28,29] using different stencils for spatial discretization, we borrow the thought of limiter for discontinuous Galerkin (DG) method to control the spurious oscillations, after this procedure, the scheme would avoid the oscillations by using HWENO reconstruction nearby discontinuities, and using linear approximation straightforwardly in the smooth regions is to increase the efficiency of the scheme. Moreover, the scheme still keeps the compactness as only immediate neighbor information is needed in the reconstruction. A collection of benchmark numerical tests for one and two dimensional cases are performed to demonstrate the numerical accuracy, high resolution and robustness of the proposed scheme. [ABSTRACT FROM AUTHOR]

It is an important and challenging issue for the numerical solution of radiative transfer equations to maintain both high order accuracy and positivity. For the two-dimensional radiative transfer equations, Ling et al. give a counterexample (Ling et al. (2018) [13]) showing that unmodulated discontinuous Galerkin (DG) solver based either on the P k or Q k polynomial spaces could generate negative cell averages even if the inflow boundary value and the source term are both positive (and, for time dependent problems, also a nonnegative initial condition). Therefore the positivity-preserving frameworks in Zhang and Shu (2010) [28] and Zhang et al. (2012) [29] which are based on the value of cell averages being positive cannot be directly used to obtain a high order conservative positivity-preserving DG scheme for the radiative transfer equations neither on rectangular meshes nor on triangular meshes. In Yuan et al. (2016) [26] , when the cell average of DG schemes is negative, a rotational positivity-preserving limiter is constructed which could keep high order accuracy and positivity in the one-dimensional radiative transfer equations with P k polynomials and could be straightforwardly extended to two-dimensional radiative transfer equations on rectangular meshes with Q k polynomials (tensor product polynomials). This paper presents an extension of the idea of the above mentioned one-dimensional rotational positivity-preserving limiter algorithm to two-dimensional high order positivity-preserving DG schemes for solving steady and unsteady radiative transfer equations on triangular meshes with P k polynomials. The extension of this method is conceptually plausible but highly nontrivial. We first focus on finding a special quadrature rule on a triangle which should satisfy some conditions. The most important one is that the quadrature points can be arranged on several line segments, on which we can use the one-dimensional rotational positivity-preserving limiter. Since the number of the quadrature points is larger than the number of basis functions of P k polynomial space, we determine a k -th polynomial by a L 2 -norm Least Square subject to its cell average being equal to the weighted average of the values on the quadrature points after using the rotational positivity-preserving limiter. Since the weights used here are the quadrature weights which are positive, then the cell average of the modified polynomial is nonnegative. And the final modified polynomial can be obtained by using the two-dimensional scaling positivity-preserving limiter on the triangular element. We theoretically prove that our rotational positivity-preserving limiter on triangular meshes could keep both high order accuracy and positivity. It is relatively simple to implement, and also does not affect convergence to weak solutions. The numerical results validate the high order accuracy and the positivity-preserving properties of our schemes. The advantage of the triangular meshes on handling complex domain is also presented in our numerical examples. • High order positivity-preserving limiter for discontinuous Galerkin methods to solve radiative transfer equations on triangular meshes. • We theoretically proved that the rotational positivity-preserving limiter on triangular meshes could keep both high order accuracy and positivity. • The numerical results validate the high order accuracy and the positivity-preserving property of our schemes. [ABSTRACT FROM AUTHOR]

Abstract: In this paper we generalize a new type of limiters based on the weighted essentially non-oscillatory (WENO) finite volume methodology for the Runge–Kutta discontinuous Galerkin (RKDG) methods solving nonlinear hyperbolic conservation laws, which were recently developed in [32] for structured meshes, to two-dimensional unstructured triangular meshes. The key idea of such limiters is to use the entire polynomials of the DG solutions from the troubled cell and its immediate neighboring cells, and then apply the classical WENO procedure to form a convex combination of these polynomials based on smoothness indicators and nonlinear weights, with suitable adjustments to guarantee conservation. The main advantage of this new limiter is its simplicity in implementation, especially for the unstructured meshes considered in this paper, as only information from immediate neighbors is needed and the usage of complicated geometric information of the meshes is largely avoided. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good performance of this procedure. [Copyright &y& Elsevier]