A discontinuous Galerkin method for time fractional diffusion equations with variable coefficients

From MaRDI portal



Abstract: We propose a piecewise-linear, time-stepping discontinuous Galerkin method to solve numerically a time fractional diffusion equation involving Caputo derivative of order muin(0,1) with variable coefficients. For the spatial discretization, we apply the standard piecewise linear continuous Galerkin method. Well-posedness of the fully discrete scheme and error analysis will be shown. For a time interval~(0,T) and a spatial domain~Omega, our analysis suggest that the error in -norm is of order O(k2−fracmu2+h2) (that is, short by order fracmu2 from being optimal in time) where k denotes the maximum time step, and h is the maximum diameter of the elements of the (quasi-uniform) spatial mesh. However, our numerical experiments indicate optimal O(k2+h2) error bound in the stronger -norm. Variable time steps are used to compensate the singularity of the continuous solution near t=0.


The authors consider a piecewise-linear, time-stepping discontinuous Galerkin method to solve a time-fractional diffusion equation with variable coefficients involving a Caputo derivative of order \(\mu \in (0,1)\). A standard Galerkin method of total degree \(\leq 1\) on each spatial mesh elements is considered as a spatial discretization. Well-posedness of the fully discrete scheme and error analysis are proved. An error estimate is proved to be of order \(k^{2-\frac{\mu}{2}}+h^2\) in the \(L^2(L^2)\)-norm. Some numerical examples are presented and it is remarked that the order is \(k^ 2 +h^2 \) in the stronger \(L^\infty (L^2)\)-norm. Variable time steps are used to compensate the singularity of the continuous solution near \(t=0\).



Cites work


Cited in
(40)








This page was built for publication: A discontinuous Galerkin method for time fractional diffusion equations with variable coefficients

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