Wave propagation problems arise in a wide range of applications. The energy conserving property is one of the guiding principles for numerical algorithms, in order to minimize the phase or shape errors after long time integration. In this paper, we develop and analyze a local discontinuous Galerkin (LDG) method for solving the wave equation. We prove optimal error estimates, superconvergence toward a particular projection of the exact solution, and the energy conserving property for the semi-discrete formulation. The analysis is extended to the fully discrete LDG scheme, with the centered second-order time discretization (the leap-frog scheme). Our numerical experiments demonstrate optimal rates of convergence and superconvergence. We also show that the shape of the solution, after long time integration, is well preserved due to the energy conserving property.

The convergence to steady state solutions of the Euler equations for weighted compact nonlinear schemes (WCNS) [Deng X. and Zhang H. (2000), J. Comput. Phys. 165, 22-44 and Zhang S., Jiang S and Shu C.-W. (2008), J. Comput. Phys. 227, 7294-7321] is studied through numerical tests. Like most other shock capturing schemes, WCNS also suffers from the problem that the residue can not settle down to machine zero for the computation of the steady state solution which contains shock waves but hangs at the truncation error level. In this paper, the techniques studied in [Zhang S. and Shu. C.-W. (2007), J. Sci. Comput. 31, 273-305 and Zhang S., Jiang S and Shu. C.-W. (2011), J. Sci. Comput. 47, 216-238], to improve the convergence to steady state solutions for WENO schemes, are generalized to the WCNS. Detailed numerical studies in one and two dimensional cases are performed. Numerical tests demonstrate the effectiveness of these techniques when applied to WCNS. The residue of various order WCNS can settle down to machine zero for typical cases while the small post-shock oscillations can be removed.

We further investigate the high order positivity-preserving discontinuous Galerkin (DG) methods for linear hyperbolic and radiative transfer equations developed in \cite{yuan}. The DG methods in \cite{yuan} can maintain positivity and high order accuracy, but they rely both on the scaling limiter in \cite{zhang1} and a rotational limiter, the latter may alter cell averages of the unmodulated DG scheme, thereby affecting conservation. Even though a Lax-Wendroff type theorem is proved in \cite{yuan}, guaranteeing convergence to weak solutions with correct shock speed when such rotational limiter is applied, it would still be desirable if a conservative DG method without changing the cell averages can be obtained which has both high order accuracy and positivity-preserving capability. In this paper, we develop and analyze such a DG method for both linear hyperbolic equations and radiative transfer equations. In the one-dimensional case, the method uses traditional DG space $P^k$ of piecewise polynomials of degree at most $k$. A key result is proved that the unmodulated DG solver in this case can maintain positivity of the cell average if the inflow boundary value and the source term are both positive, therefore the positivity-preserving framework in \cite{zhang1} can be used to obtain a high order conservative positivity-preserving DG scheme. Unfortunately, in two-dimensions this is no longer the case. We show that the unmodulated DG solver based either on $P^k$ or $Q^k$ spaces (piecewise $k$-th degree polynomials or piecewise tensor-product $k$-th degree polynomials) could generate negative cell averages. We augment the DG space with additional functions so that the positivity of cell averages from the unmodulated DG solver can be restored, thereby leading to high order conservative positivity-preserving DG scheme based on these augmented DG spaces following the framework in \cite{zhang1}. Computational results are provided to demonstrate the good performance of our DG schemes.

We present novel algorithms for compressible flows that\n\nare efficient for all Mach numbers. The approach is based on several\n\ningredients: semi-implicit schemes, the gauge decomposition\n\nof the velocity field and a second order formulation of the density\n\nequation (in the isentropic case) and of the energy equation (in\n\nthe full Navier-Stokes case). Additionally, we show that our approach\n\ncorresponds to a micro-macro decomposition of the model,\n\nwhere the macro field corresponds to the incompressible component\n\nsatisfying a perturbed low Mach number limit equation and\n\nthe micro field is the potential component of the velocity. Finally,\n\nwe also use the conservative variables in order to obtain a proper\n\nconservative formulation of the equations when the Mach number\n\nis order unity. We successively consider the isentropic case, the\n\nfull Navier-Stokes case, and the isentropic Navier-Stokes-Poisson\n\ncase. In this work, we only concentrate on the question of the\n\ntime discretization and show that the proposed method leads to\n\nAsymptotic Preserving schemes for compressible flows in the low\n\nMach number limit.

In this paper, we analyze the Lax-Wendroff discontinuous Galerkin (LWDG) method for solving linear conservation laws. The method was originally proposed by Guo et al., where they applied local discontinuous Galerkin (LDG) techniques to approximate high order spatial derivatives in the Lax-Wendroff time discretization. We show that, under the standard CFL condition $\tau\leq \lambda h$ (where $\tau$ and $h$ are the time step and the maximum element length respectively and $\lambda>0$ is a constant) and uniform or non-increasing time steps, the second order schemes with piecewise linear elements and the third order schemes with arbitrary piecewise polynomial elements are stable in the $L^2$ norm. The specific type of stability may differ with different choices of numerical fluxes. Our stability analysis includes multidimensional problems with divergence-free coefficients. Besides solving the equation itself, the LWDG method also gives approximations to its time derivative simultaneously. We obtain optimal error estimates for both the solution $u$ and its first order time derivative $u_t$ in one dimension, and numerical examples are given to validate our analysis.