Adaptive discontinuous Galerkin approximations to fourth order parabolic problems

An adaptive algorithm, based on residual type a posteriori indicators of errors measured in $L^{\infty}(L^2)$ and $L^2(L^2)$ norms, for a numerical scheme consisting of implicit Euler method in time and discontinuous Galerkin method in space for linear parabolic fourth order problems is presented. The a posteriori analysis is performed for convex domains in two and three space dimensions for local spatial polynomial degrees $r\ge 2$. The a posteriori estimates are then used within an adaptive algorithm, highlighting their relevance in practical computations, which results into substantial reduction of computational effort

A posteriori error control for discontinuous Galerkin methods for parabolic problems

We derive energy-norm a posteriori error bounds for an Euler time-stepping method combined with various spatial discontinuous Galerkin schemes for linear parabolic problems. For accessibility, we address first the spatially semidiscrete case, and then move to the fully discrete scheme by introducing the implicit Euler time-stepping. All results are presented in an abstract setting and then illustrated with particular applications. This enables the error bounds to hold for a variety of discontinuous Galerkin methods, provided that energy-norm a posteriori error bounds for the corresponding elliptic problem are available. To illustrate the method, we apply it to the interior penalty discontinuous Galerkin method, which requires the derivation of novel a posteriori error bounds. For the analysis of the time-dependent problems we use the elliptic reconstruction technique and we deal with the nonconforming part of the error by deriving appropriate computable a posteriori bounds for it.Comment: 6 figure

Krylov implicit integration factor discontinuous Galerkin methods on sparse grids for high dimensional reaction-diffusion equations

Computational costs of numerically solving multidimensional partial differential equations (PDEs) increase significantly when the spatial dimensions of the PDEs are high, due to large number of spatial grid points. For multidimensional reaction-diffusion equations, stiffness of the system provides additional challenges for achieving efficient numerical simulations. In this paper, we propose a class of Krylov implicit integration factor (IIF) discontinuous Galerkin (DG) methods on sparse grids to solve reaction-diffusion equations on high spatial dimensions. The key ingredient of spatial DG discretization is the multiwavelet bases on nested sparse grids, which can significantly reduce the numbers of degrees of freedom. To deal with the stiffness of the DG spatial operator in discretizing reaction-diffusion equations, we apply the efficient IIF time discretization methods, which are a class of exponential integrators. Krylov subspace approximations are used to evaluate the large size matrix exponentials resulting from IIF schemes for solving PDEs on high spatial dimensions. Stability and error analysis for the semi-discrete scheme are performed. Numerical examples of both scalar equations and systems in two and three spatial dimensions are provided to demonstrate the accuracy and efficiency of the methods. The stiffness of the reaction-diffusion equations is resolved well and large time step size computations are obtained

An a posteriori error analysis of a mixed finite element Galerkin approximation to second order linear parabolic problems

In this article, a posteriori error estimates are derived for a mixed finite element Galerkin approximation to second order linear parabolic initial and boundary value problems. Using mixed elliptic reconstruction method, a posteriori error estimates in $L^\infty(L^2)$ and $L^2(L^2)$-norms with optimal order of convergence for the solution as well as its flux are proved for the semidiscrete scheme. Finally, based on backward Euler method, a completely discrete scheme is analyzed and a posteriori bounds are derived, which improves earlier results on a posteriori estimates for mixed parabolic problems