In this paper, we investigate a simple limiter using weighted essentially non-oscillatory (WENO) methodology for the Runge-Kutta discontinuous Galerkin (RKDG) methods solving conservation laws, with the goal of obtaining a robust and high order limiting procedure to simultaneously achieve uniform high order accuracy and sharp, non-oscillatory shock transitions. The idea of this limiter is to reconstruct the entire polynomial, instead of reconstructing point values or moments in the classical WENO reconstructions. That is, the reconstruction polynomial on the target cell is a convex combination of polynomials on this cell and its neighboring cells and the nonlinear weights of the convex combination follow the classical WENO procedure. The main advantage of this limiter is its simplicity in implementation, especially for multi-dimensional meshes. Numerical results in one and two dimensions are provided to illustrate the behavior of this procedure.

For a class of convection-diffusion equations with variable diffusivity, we construct third order accurate discontinuous Galerkin (DG) schemes on both one and two dimensional rectangular meshes. The DG method with an explicit time stepping can well be applied to nonlinear convection–diffusion equations. It is shown that under suitable time step restrictions, the scaling limiter proposed in Liu and Yu (2014) [23]when coupled with the present DG schemes preserves the solution bounds indicated by the initial data, i.e., the maximum principle, while maintaining uniform third order accuracy. These schemes can be extended to rectangular meshes in three dimension. The crucial for all model scenarios is that an effective test set can be identified to verify the desired bounds of numerical solutions. This is achieved mainly by taking advantage of the flexible form of the diffusive flux and the adaptable decomposition of weighted cell averages. Numerical results are presented to validate the numerical methods and demonstrate their effectiveness.

This is an extension of our earlier work \cite{DS} in which a high order stable method was constructed for solving hyperbolic conservation laws on arbitrarily distributed point clouds. %The initial condition is given on such a point cloud, %and the algorithm solves for point values of the solution at %later time also on this point cloud. An algorithm of building a suitable polygonal mesh based on the random points was given and the traditional discontinuous Galerkin (DG) method was adopted on the constructed polygonal mesh. Numerical results in \cite{DS} show that the current scheme will generate spurious numerical oscillations when dealing with solutions containing strong shocks. In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for DG schemes on two-dimensional unstructured triangular meshes \cite{ZZSQ}, to our high order method on polygonal meshes. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the original method in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt is particularly simple to implement and will not harm the conservativeness and compactness of the original method. Moreover, we also extend the maximum-principle-satisfying limiter for the scalar case and the positivity-preserving limiter for the Euler system to our method. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good behavior of these limiters.

No abstract uploaded!

In this paper we present an error estimate for the explicit Runge-Kutta discontinuous Galerkin method to solve a linear hyperbolic equation in one dimension with discontinuous but piecewise smooth initial data. The discontinuous finite element space is made up of piecewise polynomials of arbitrary degree $k\geq1$, and time is advanced by the third order explicit total variation diminishing Runge-Kutta method under the standard CFL temporal-spatial condition. The $L^2(\mathbb{R}\backslash\mathcal{R}_T)$-norm error at the final time $T$ is optimal in both space and time, where $\mathcal{R}_T$ is the pollution region due to the initial discontinuity with the width $\mathcal{O}(\sqrt{T\beta}h^{1/2}\log(1/h))$. Here $h$ is the maximum cell length and $\beta$ is the flowing speed. These results are independent of the time step and hold also for the semi-discrete discontinuous Galerkin method.