In this paper, we develop and analyze a new multiscale discontinuous Galerkin (DG) method for one-dimensional stationary Schr\"{o}dinger equations with open boundary conditions \bo{which have highly oscillating solutions}. Our method uses a smaller finite element space than the WKB local DG (WKB-LDG) method while achieving the same order of accuracy with no resonance errors. We prove that the DG approximation converges optimally with respect to the mesh size $h$ in $L^2$ norm without the typical constraint \bo{that $h$ has to be smaller than the wave length}. Numerical experiments were carried out to verify the second order optimal convergence rate of the method and to demonstrate its ability to capture oscillating solutions on coarse meshes in the applications to Schr\"{o}dinger equations.

We study semi-Lagrangian discontinuous Galerkin (SLDG) and Runge-Kutta discontinuous Galerkin (RKDG) schemes for some front propagation problems in the presence of an obstacle term, modeled by a nonlinear Hamilton-Jacobi equation of the form $\min(u_t + c u_x, u - g(x))=0$, in one space dimension. New convergence results and error bounds are obtained for Lipschitz regular data. These ``low regularity" assumptions are the natural ones for the solutions of the studied equations. Numerical tests are given to illustrate the behavior of our schemes.

In this paper we consider the drift-diffusion (DD) model of one dimensional semiconductor devices, which is a system involving not only first derivative convection terms but also second derivative diffusion terms and a coupled Poisson potential equation. Optimal error estimates are obtained for both the semi-discrete and fully discrete local discontinuous Galerkin (LDG) schemes with smooth solutions. In the fully discrete scheme, we couple the implicit-explicit (IMEX) time discretization with the LDG spatial discretization, in order to allow larger time steps and to save computational cost. The main technical difficulty in the analysis is to treat the inter-element jump terms which arise from the discontinuous nature of the numerical method and the nonlinearity and coupling of the models. A simulation is also performed to validate the analysis.

In this paper, we propose a new type of weighted essentially non-oscillatory (WENO) limiter, which belongs to the class of Hermite WENO (HWENO) limiters, for the Runge-Kutta discontinuous Galerkin (RKDG) methods solving hyperbolic conservation laws. This new HWENO limiter is a modification of the simple WENO limiter proposed recently by Zhong and Shu. Both limiters use information of the DG solutions only from the target cell and its immediate neighboring cells, thus maintaining the original compactness of the DG scheme. The goal of both limiters is to obtain high order accuracy and non-oscillatory properties simultaneously. The main novelty of the new HWENO limiter in this paper is to reconstruct the polynomial on the target cell in a least square fashion while the simple WENO limiter is to use the entire polynomial of the original DG solutions in the neighboring cells with an addition of a constant for conservation. The modification in this paper improves the robustness in the computation of problems with strong shocks or contact discontinuities, without changing the compact stencil of the DG scheme. Numerical results for both one and two dimensional equations including Euler equations of compressible gas dynamics are provided to illustrate the viability of this modified limiter.

The main purpose of this paper is to analyze the stability and error estimates of the local discontinuous Galerkin (LDG) methods coupled with implicit-explicit (IMEX) time discretization schemes, for solving multi-dimensional convection-diffusion equations with nonlinear convection. By establishing the important relationship between the gradient and the interface jump of the numerical solution with the independent numerical solution of the gradient in the LDG method, on both rectangular and triangular elements, we can obtain the same stability results as in the one-dimensional case, i.e, the IMEX LDG schemes are unconditionally stable for the multi-dimensional convection-diffusion problems, in the sense that the time-step $\dt$ is only required to be upper-bounded by a positive constant independent of the spatial mesh size $h$. Furthermore, by the aid of the so-called elliptic projection and the adjoint argument, we can also obtain optimal error estimates in both space and time, for the corresponding fully discrete IMEX LDG schemes, under the same condition, i.e, if piecewise polynomial of degree $k$ is adopted on either rectangular or triangular meshes, we can show the convergence accuracy is of order $\Ocal(h^{k+1}+\dt^s)$ for the $s$-th order IMEX LDG scheme ($s=1,2,3$) under consideration. Numerical experiments are also given to verify our main results.