The main purpose of this paper is to analyze the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with carefully chosen implicit-explicit (IMEX) Runge-Kutta time discretization up to third order accuracy, for solving one-dimensional linear advection-diffusion equations. In the time discretization the advection term is treated explicitly and the diffusion term implicitly. There are three highlights of this work. The first is that we establish an important relationship between the gradient and interface jump of the numerical solution with the independent numerical solution of the gradient in the LDG methods. The second is that, by aid of the aforementioned relationship and the energy method, we show that the IMEX LDG schemes are unconditionally stable for the linear problems, in the sense that the time-step $\dt$ is only required to be upper-bounded by a constant which depends on the ratio of the diffusion and the square of the advection coefficients and is independent of the mesh-size $h$, even though the advection term is treated explicitly. The last is that under this time step condition, we obtain optimal error estimates in both space and time for the third order IMEX Runge-Kutta time-marching coupled with LDG spatial discretization. Numerical experiments are also given to verify the main results.

In this paper, we adapt a simple weighted essentially non-oscillatory (WENO) limiter, originally designed for discontinuous Galerkin (DG) schemes on two-dimensional unstructured triangular meshes, to the correction procedure via reconstruction (CPR) framework for solving nonlinear hyperbolic conservation laws on two-dimensional unstructured triangular meshes with straight or curved edges. This is an extension of our earlier work in which the WENO limiter was designed for the CPR framework on regular meshes. The objective of this simple WENO limiter is to simultaneously maintain uniform high order accuracy of the CPR framework in smooth regions and control spurious numerical oscillations near discontinuities. The WENO limiter we adopt in this paper uses information only from the target cell and its immediate neighbors. Hence, it is particularly simple to implement and will not harm the conservativeness and compactness of the CPR framework. Since the CPR framework with this WENO limiter does not in general satisfy the positivity preserving property, we also extend the positivity-preserving limiters to the CPR framework. Numerical results for both scalar equations and Euler systems of compressible gas dynamics are provided to illustrate the good behavior of this procedure.

In this paper, we develop a class of nonlinear compact schemes based on our previous linear central compact schemes with spectral-like resolution [X. Liu, S. Zhang, H. Zhang and C-W. Shu, A new class of central compact schemes with spectral-like resolution I: Linear schemes, Journal of Computational Physics 248 (2013) 235-256]. In our approach, we compute the flux derivatives on the cell-nodes by the physical fluxes on the cell nodes and numerical fluxes on the cell centers. To acquire the numerical fluxes on the cell centers, we perform a weighted hybrid interpolation of an upwind interpolation and a central interpolation. Through systematic analysis and numerical tests, we show that our nonlinear compact scheme has high order, high resolution and low dissipation, and has the same ability to capture strong discontinuities as regular weighted essentially non-oscillatory (WENO) schemes. It is a good choice for the simulation of multiscale problems with shock waves.

In this paper we present an {\em a priori} error estimate of the Runge-Kutta discontinuous Galerkin method for solving symmetrizable conservation laws, where the time is discretized with the third order explicit total variation diminishing Runge-Kutta method and the finite element space is made up of piecewise polynomials of degree $k\geq 2$. Quasi-optimal error estimate is obtained by energy techniques, for the so-called generalized E-fluxes under the standard temporal-spatial CFL condition $\tau\leq\gamma h$, where $h$ is the element length and $\tau$ is time step, and $\gamma$ is a positive constant independent of $h$ and $\tau$. Optimal estimates are also considered when the upwind numerical flux is used.

We consider an incompressible viscous flow without surface tension in a finite-depth domain of two dimensions, with free top boundary and fixed bottom boundary. This system is governed by the Navier-Stokes equations in this moving domain and the transport equation on the moving boundary. In this paper, we construct a stable numerical scheme to simulate the evolution of this system by discontinuous Galerkin method, and discuss the error analysis of the fluid under certain assumptions. Our formulation is mainly based on the geometric structure introduced earlier in the literature, and the natural energy estimate, which is rarely used in the numerical study of this system before.