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.

In this paper, we discuss high order finite difference weighted essentially non-oscillatory (WENO) schemes, coupled with total variation diminishing (TVD) Runge-Kutta (RK) temporal integration, for solving the semilinear hyperbolic system of a correlated random walk model describing movement of animals and cells in biology. Since the solutions to this system are non-negative, we discuss a positivity-preserving limiter without compromising accuracy. Analysis is performed to justify the maintanance of third order spatial / temporal accuracy when the limiters are applied to a third order finite difference scheme and third order TVD-RK time discretization for solving this model. Numerical results are also provided to demonstrate these methods up to fifth order accuracy.

In this paper, we study the superconvergence of the error for the local discontinuous Galerkin (LDG) finite element method for one-dimensional linear parabolic equations when the alternating flux is used. We prove that if we apply piecewise $k$-th degree polynomials, the error between the LDG solution and the exact solution is ($k+2$)-th order superconvergent at the Radau points with suitable initial discretization. Moreover, we also prove the LDG solution is ($k+2$)-th order superconvergent for the error to a particular projection of the exact solution. Even though we only consider periodic boundary condition, this boundary condition is not essential, since we do not use Fourier analysis. Our analysis is valid for arbitrary regular meshes and for $\mathcal{P}^k$ polynomials with arbitrary $k\geq1$. We perform numerical experiments to demonstrate that the superconvergence rates proved in this paper are sharp.

Fourier spectral methods achieve exponential accuracy both on the approximation level and for solving partial differential equations (PDEs), if the solution is analytic. If the solution is discontinuous but piecewise analytic up to the discontinuities, Fourier spectral methods produce poor pointwise accuracy, but still maintains exponential accuracy after post-processing. In earlier work, an extended technique is provided to recover exponential accuracy for functions which have end-point singularities, from the knowledge of point values on standard collocation points. In this paper, we develop a technique to recover exponential accuracy from the first $N$ Fourier coefficients of functions which are analytic in the open interval but have unbounded derivative singularities at end points. With this post-processing method, we are able to obtain exponential accuracy of spectral methods applied to linear transport equations involving such functions.