In this paper, we introduce a general approach for proving optimal $L^2$ error estimates for the semi-discrete local discontinuous Galerkin (LDG) methods solving linear high order wave equations. The optimal order of error estimates hold not only for the solution itself but also for the auxiliary variables in the LDG method approximating the various order derivatives of the solution. Several examples including the one-dimensional third order wave equation, one-dimensional fifth order wave equation, and multi-dimensional Schr\"{o}dinger equation are explored to demonstrate this approach. The main idea is to derive energy stability for the various auxiliary variables in the LDG discretization, via using the scheme and its time derivatives with different test functions. Special projections are utilized to eliminate the jump terms at the cell boundaries in the error estimate in order to achieve the optimal order of accuracy.

No abstract uploaded!

In this paper, we develop third-order conservative sign-preserving and steady-state-preserving time integrations and seek their applications in multispecies and multireaction chemical reactive flows. In this problem, the density and pressure are nonnegative, and the mass fraction for the ith species, denoted as zi, 1<=i<=M, should be between 0 and 1, where M is the total number of species. There are four main difficulties in constructing high-order bound-preserving techniques for multispecies and multireaction detonations. First of all, most of the bound-preserving techniques available are based on Euler forward time integration. Therefore, for problems with stiff source, the time step will be significantly limited. Secondly, the mass fraction does not satisfy a maximum-principle and hence it is not easy to preserve the upper bound 1. Thirdly, in most of the previous works for gaseous denotation, the algorithm relies on second-order Strang splitting methods where the ux and stiff source terms can be solved separately, and the extension to high-order time discretization seems to be complicated. Finally, most of the previous ODE solvers for stiff problems cannot preserve the total mass and the positivity of the numerical approximations at the same time. In this paper, we will construct third-order conservative sign-preserving Rugne-Kutta and multistep methods to overcome all these difficulties. The time integrations do not depend on the Strang splitting, i.e. we do not split the flux and the stiff source terms. Moreover, the time discretization can handle the stiff source with large time step and preserves the steady-state. Numerical experiments will be given to demonstrate the good performance of the bound-preserving technique and the stability of the scheme for problems with stiff source terms.

In this paper, we apply discontinuous Galerkin (DG) methods to solve two model equations: Krause's consensus models and pressureless Euler equations. These two models are used to describe the collisions of particles, and the distributions can be identified as density functions. If the particles are placed at a single point, then the density function turns out to be a $\delta$-function and is difficult to be well approximated numerically. In this paper, we use DG method to approximate such a singularity and demonstrate the good performance of the scheme. Since the density functions are always positive, we apply a positivity-preserving limiter to them. Moreover, for pressureless Euler equations, the velocity satisfies the maximum principle. We also construct special limiters to fulfill this requirement.

No abstract uploaded!