Convergence and superconvergence analyses of HDG methods for time fractional diffusion problems

DOI10.1007/S10444-015-9428-XMaRDI QIDQ273668zbMATH OpenOpenAlexFDO

Authors Kassem Mustapha, M. Nour, Bernardo Cockburn

Publication date 22 April 2016

Published in Advances in Computational Mathematics (Search for Journal in Brave)

Full work available at URL https://arxiv.org/abs/1412.2098

convergence analysis discontinuous Galerkin methods hybridization time fractional

Fractional derivatives and integrals (26A33) Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs (65M12)

Abstract: We study the hybridizable discontinuous Galerkin (HDG) method for the spatial discretization of time fractional diffusion models with Caputo derivative of order

0 < a l p h a < 1

. For each time

t i n [0, T]

, the HDG approximations are taken to be piecewise polynomials of degree

k g e 0

on the spatial domain~

O m e g a

, the approximations to the exact solution

u

in the

-norm and to

a b l a u

in the

-norm are proven to converge with the rate

h^{k + 1}

provided that

u

is sufficiently regular, where

h

is the maximum diameter of the elements of the mesh. Moreover, for

k g e 1

, we obtain a superconvergence result which allows us to compute, in an elementwise manner, a new approximation for

u

converging with a rate

h^{k + 2}

(ignoring the logarithmic factor), for quasi-uniform spatial meshes. Numerical experiments validating the theoretical results are displayed.

Review text

This paper studies a numerical method resulting from using exact integration in time and a hybridizable discontinuous Galerkin (HDG) method for the spatial discretization of the time fractional diffusion model problem: \(^cD^{1-\alpha}u-\Delta u = f\) in \(\Omega\times (0,T]\), subject to \(u=g(x)\) in \(\partial\Omega\times (0,T]\) and \(u(x,0)=u_0(x)\) in \(\Omega\). Here, \(\Omega\subset\mathbb R^d\) \((1\leq d\leq 3)\) is a convex polyhedral domain, \(f\), \(g\) and \(u_0\) are smooth enough to ensure the existence of a solution \(u\in W^{1,1}(0,T; H^2(\Omega))\). When the HDG approximations are piecewise polynomials of degree \(k\geq 0\), it is shown that the approximations to the exact solution \(u\) and to its gradient \(\nabla u\) in the \(L^\infty\) norm converge with the rate \(h^{k+1}\), where \(h\) is the maximum diameter of the elements of the mesh. This result is further improved by the authors in case \(k\geq 1\). Numerical experiments are also presented in order to support the theoretical findings.

Recommendations

Cites work

Cited in

(30)

This page was built for publication: Convergence and superconvergence analyses of HDG methods for time fractional diffusion problems

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q273668)