For a class of convection-diffusion equations with variable diffusivity, we construct third order accurate discontinuous Galerkin (DG) schemes on both one and two dimensional rectangular meshes. The DG method with an explicit time stepping can well be applied to nonlinear convection–diffusion equations. It is shown that under suitable time step restrictions, the scaling limiter proposed in Liu and Yu (2014) [23]when coupled with the present DG schemes preserves the solution bounds indicated by the initial data, i.e., the maximum principle, while maintaining uniform third order accuracy. These schemes can be extended to rectangular meshes in three dimension. The crucial for all model scenarios is that an effective test set can be identified to verify the desired bounds of numerical solutions. This is achieved mainly by taking advantage of the flexible form of the diffusive flux and the adaptable decomposition of weighted cell averages. Numerical results are presented to validate the numerical methods and demonstrate their effectiveness.

In this paper, an analysis of the superconvergence property of the semidiscrete discontinuous Galerkin (DG) method applied to one-dimensional time-dependent nonlinear scalar conservation laws is carried out. We prove that the error between the DG solution and a particular projection of the exact solution achieves $\left(k + \frac32\right)$-th order superconvergence when upwind fluxes are used. The results hold true for arbitrary nonuniform regular meshes and for piecewise polynomials of degree $k$ ($k \ge 1$), under the condition that $|f'(u)|$ possesses a uniform positive lower bound. Numerical experiments are provided to show that the superconvergence property actually holds true for nonlinear conservation laws with general flux functions, indicating that the restriction on $f(u)$ is artificial.

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.

No abstract uploaded!

We propose a structure-preserving doubling algorithm for a quadratic eigenvalue problem arising from the stability analysis of time-delay systems. We are particularly interested in the eigenvalues on the unit circle, which are difficult to estimate. The convergence and backward error of the algorithm are analyzed and three numerical examples are presented. Our experience shows that our algorithm is efficient in comparison to the few existing approaches for small to medium size problems.