An hp-Local Discontinuous Galerkin method for Parabolic\ud Integro-Differential Equations

In this article, a priori error analysis is discussed for an hp-local discontinuous Galerkin (LDG) approximation to a parabolic integro-differential equation. It is shown that the L2 -norm of the gradient and the L2 -norm of the potential are optimal in the discretizing parameter h and suboptimal in the degree of polynomial p. Due to the presence of the integral term, an introduction of an expanded mixed type Ritz-Volterra projection helps to achieve optimal estimates. Further, it is observed that a negative norm estimate of the gradient plays a crucial role in our convergence analysis. As in the elliptic case, similar results on order of convergence are established for the semidiscrete method after suitably modifying the numerical fluxes. The optimality of these theoretical results is tested in a series of numerical experiments on two dimensional domains

Analysis of optimal error estimates and superconvergence of the discontinuous Galerkin method for convection-diffusion problems in one space dimension

In this paper, we study the convergence and superconvergence properties of the discontinuous Galerkin (DG) method for a linear convection-diffusion problem in one-dimensional setting. We prove that the DG solution and its derivative exhibit optimal O(hp+1) and O(hp) convergence rates in the L 2 -norm, respectively, when p-degree piecewise polynomials with p ≥ 1 are used. We further prove that the p-degree DG solution and its derivative are O(h2p) superconvergent at the downwind and upwind points, respectively. Numerical experiments demonstrate that the theoretical rates are optimal and that the DG method does not produce any oscillation. We observed optimal rates of convergence and superconvergence even in the presence of boundary layers when Shishkin meshes are used

hp-Version discontinuous Galerkin finite element methods for semilinear parabolic problems

We consider the hp-version interior penalty discontinuous Galerkin finite element method (hp-DGFEM) for semilinear parabolic equations with mixed Dirichlet and Neumann boundary conditions. Our main concern is the error analysis of the hp--DGFEM on shape--regular spatial meshes. We derive error bounds under various hypotheses on the regularity of the solution, for both the symmetric and non--symmetric versions of DGFEM

A posteriori discontinuous Galerkin error estimates for transient convection–diffusion equations

AbstractA posteriori error estimates are derived for unsteady convection–diffusion equations discretized with the non-symmetric interior penalty and the local discontinuous Galerkin methods. First, an error representation formula in a user specified output functional is derived using duality techniques. Then, an Lt2(Lx2) a posteriori estimate consisting of elementwise residual-based error indicators is obtained by eliminating the dual solution. Numerical experiments are performed to assess the convergence rates of the various error indicators on a model problem

Discontinuous Galerkin Methods for Mass Transfer through Semi-Permeable Membranes

A discontinuous Galerkin (dG) method for the numerical solution of initial/boundary value multi-compartment partial differential equation (PDE) models, interconnected with interface conditions, is presented and analysed. The study of interface problems is motivated by models of mass transfer of solutes through semi-permeable membranes. More specifically, a model problem consisting of a system of semilinear parabolic advection-diffusion-reaction partial differential equations in each compartment, equipped with respective initial and boundary conditions, is considered. Nonlinear interface conditions modelling selective permeability, congestion and partial reflection are applied to the compartment interfaces. An interior penalty dG method is presented for this problem and it is analysed in the space-discrete setting. The a priori analysis shows that the method yields optimal a priori bounds, provided the exact solution is sufficiently smooth. Numerical experiments indicate agreement with the theoretical bounds and highlight the stability of the numerical method in the advection-dominated regime