Source author record

Kassem Mustapha

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.NA math.AP Numerical Analysis math.CA

Catalog footprint

Research overview

14works

4topics

4close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2020arXiv

FEM for time-fractional diffusion equations, novel optimal error analyses

A semidiscrete Galerkin finite element method applied to time-fractional diffusion equations with time-space dependent diffusivity on bounded convex spatial domains will be studied. The main focus is on achieving optimal error results with respect to both the convergence order of the approximate solution and the regularity of the initial data. By using novel energy arguments, for each fixed time $t$ , optimal error bounds in the spatial $L^2$ - and $H^1$ -norms are derived for both cases: smooth and nonsmooth initial data.

preprint2020arXiv

Uniform stability for a spatially-discrete, subdiffusive Fokker-Planck equation

We prove stability estimates for the spatially discrete, Galerkin solution of a fractional Fokker-Planck equation, improving on previous results in several respects. Our main goal is to establish that the stability constants are bounded uniformly in the fractional diffusion exponent $α\in(0,1]$ . In addition, we account for the presence of an inhomogeneous term and show a stability estimate for the gradient of the Galerkin solution. As a by-product, the proofs of error bounds for a standard finite element approximation are simplified.

preprint2019arXiv

Regularity theory for time-fractional advection-diffusion-reaction equations

We investigate the behavior of the time derivatives of the solution to a linear time-fractional, advection-diffusion-reaction equation, allowing space- and time-dependent coefficients as well as initial data that may have low regularity. Our focus is on proving estimates that are needed for the error analysis of numerical methods. The nonlocal nature of the fractional derivative creates substantial difficulties compared with the case of a classical parabolic PDE. In our analysis, we rely on novel energy methods in combination with a fractional Gronwall inequality and certain properties of fractional integrals.

preprint2019arXiv

Well-posedness of time-fractional, advection-diffusion-reaction equations

We establish the well-posedness of an initial-boundary value problem for a general class of time-fractional, advection-diffusion-reaction equations, allowing space- and time-dependent coefficients as well as initial data that may have low regularity. Our analysis relies on novel energy methods in combination with a fractional Gronwall inequality and properties of fractional integrals.

preprint2016arXiv

An inverse source problem for a two-parameter anomalous diffusion with local time datum

We determine the space-dependent source term for a two-parameter fractional diffusion problem subject to nonlocal non-self-adjoint boundary conditions and two local time-distinct datum. A bi-orthogonal pair of bases is used to construct a series representation of the solution and the source term. The two local time conditions spare us from measuring the fractional integral initial conditions commonly associated with fractional derivatives. On the other hand, they lead to delicate $2\times 2$ linear systems for the Fourier coefficients of the source term and of the fractional integral of the solution at $t=0$ . The asymptotic behavior and estimates of the generalized Mittag-Leffler function are used to establish the solvability of these linear systems, and to obtain sufficient conditions for the existence of our construction. Analytical and numerical examples are provided.

preprint2015arXiv

Finite volume element method for two-dimensional fractional subdiffusion problems

In this paper, a semi-discrete spatial finite volume (FV) method is proposed and analyzed for approximating solutions of anomalous subdiffusion equations involving a temporal fractional derivative of order $α\in (0,1)$ in a two-dimensional convex polygonal domain. Optimal error estimates in $L^\infty(L^2)$ - norm is shown to hold. Superconvergence result is proved and as a consequence, it is established that quasi-optimal order of convergence in $L^{\infty}(L^{\infty})$ holds. We also consider a fully discrete scheme that employs FV method in space, and a piecewise linear discontinuous Galerkin method to discretize in temporal direction. It is, further, shown that convergence rate is of order $O(h^2+k^{1+α}),$ where $h$ denotes the space discretizing parameter and $k$ represents the temporal discretizing parameter. Numerical experiments indicate optimal convergence rates in both time and space, and also illustrate that the imposed regularity assumptions are pessimistic.

preprint2015arXiv

Numerical solution of the time-fractional Fokker-Planck equation with general forcing

We study two schemes for a time-fractional Fokker-Planck equation with space- and time-dependent forcing in one space dimension. The first scheme is continuous in time and is discretized in space using a piecewise-linear Galerkin finite element method. The second is continuous in space and employs a time-stepping procedure similar to the classical implicit Euler method. We show that the space discretization is second-order accurate in the spatial $L_2$ -norm, uniformly in time, whereas the corresponding error for the time-stepping scheme is $O(k^α)$ for a uniform time step $k$ , where $α\in(1/2,1)$ is the fractional diffusion parameter. In numerical experiments using a combined, fully-discrete method, we observe convergence behaviour consistent with these results.

preprint2014arXiv

A discontinuous Petrov-Galerkin method for time-fractional diffusion equations

We propose and analyze a time-stepping discontinuous Petrov-Galerkin method combined with the continuous conforming finite element method in space for the numerical solution of time-fractional subdiffusion problems. We prove the existence, uniqueness and stability of approximate solutions, and derive error estimates. To achieve high order convergence rates from the time discretizations, the time mesh is graded appropriately near~ $t=0$ to compensate the singular (temporal) behaviour of the exact solution near $t=0$ caused by the weakly singular kernel, but the spatial mesh is quasiuniform. In the $L_\infty((0,T);L_2(Ω))$ -norm ( $(0,T)$ is the time domain and $Ω$ is the spatial domain), for sufficiently graded time meshes, a global convergence of order $k^{m+α/2}+h^{r+1}$ is shown, where $0<α<1$ is the fractional exponent, $k$ is the maximum time step, $h$ is the maximum diameter of the spatial finite elements, and $m$ and $r$ are the degrees of approximate solutions in time and spatial variables, respectively. Numerical experiments indicate that our theoretical error bound is pessimistic. We observe that the error is of order ~ $k^{m+1}+h^{r+1}$ , that is, optimal in both variables.

preprint2014arXiv

A hybridizable discontinuous Galerkin method for fractional diffusion problems

We study the use of the hybridizable discontinuous Galerkin (HDG) method for numerically solving fractional diffusion equations of order $-α$ with $-1<α<0$ . For exact time-marching, we derive optimal algebraic error estimates {assuming} that the exact solution is sufficiently regular. Thus, if for each time $t \in [0,T]$ the approximations are taken to be piecewise polynomials of degree $k\ge0$ on the spatial domain~ $Ω$ , the approximations to $u$ in the $L_\infty\bigr(0,T;L_2(Ω)\bigr)$ -norm and to $\nabla u$ in the $L_\infty\bigr(0,T;{\bf L}_2(Ω)\bigr)$ -norm are proven to converge with the rate $h^{k+1}$ , where $h$ is the maximum diameter of the elements of the mesh. Moreover, for $k\ge1$ and quasi-uniform meshes, we obtain a superconvergence result which allows us to compute, in an elementwise manner, a new approximation for $u$ converging with a rate of $\sqrt{\log(T h^{-2/(α+1)})}\, \,h^{k+2}$ .

preprint2014arXiv

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

We study the hybridizable discontinuous Galerkin (HDG) method for the spatial discretization of time fractional diffusion models with Caputo derivative of order $0<α<1$ . For each time $t \in [0,T]$ , the HDG approximations are taken to be piecewise polynomials of degree $k\ge0$ on the spatial domain~ $Ω$ , the approximations to the exact solution $u$ in the $L_\infty\bigr(0,T;L_2(Ω)\bigr)$ -norm and to $\nabla u$ in the $L_\infty\bigr(0,T;{\bf L}_2(Ω)\bigr)$ -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\ge1$ , 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.

preprint2014arXiv

Time-stepping discontinuous Galerkin methods for fractional diffusion problems

Time-stepping $hp$ -versions discontinuous Galerkin (DG) methods for the numerical solution of fractional subdiffusion problems of order $-α$ with $-1<α<0$ will be proposed and analyzed. Generic $hp$ -version error estimates are derived after proving the stability of the approximate solution. For $h$ -version DG approximations on appropriate graded meshes near $t=0$ , we prove that the error is of order $O(k^{\max\{2,p\}+\fracα{2}})$ , where $k$ is the maximum time-step size and $p\ge 1$ is the uniform degree of the DG solution. For $hp$ -version DG approximations, by employing geometrically refined time-steps and linearly increasing approximation orders, exponential rates of convergence in the number of temporal degrees of freedom are shown. Finally, some numerical tests are given.

preprint2014arXiv

Time-stepping error bounds for fractional diffusion problems with non-smooth initial data

We apply the piecewise constant, discontinuous Galerkin method to discretize a fractional diffusion equation with respect to time. Using Laplace transform techniques, we show that the method is first order accurate at the \$n\$th time level \$t_n\$, but the error bound includes a factor \$t_n^{-1}\$ if we assume no smoothness of the initial data. We also show that for smoother initial data the growth in the error bound as \$t_n\$ decreases is milder, and in some cases absent altogether. Our error bounds generalize known results for the classical heat equation and are illustrated for a model problem.

preprint2013arXiv

A Superconvergent discontinuous Galerkin method for Volterra integro-differential equations, smooth and non-smooth kernels

We study the numerical solution for Volerra integro-differential equations with smooth and non-smooth kernels. We use a $h$ -version discontinuous Galerkin (DG) method and derive nodal error bounds that are explicit in the parameters of interest. In the case of non-smooth kernel, it is justified that the start-up singularities can be resolved at superconvergence rates by using non-uniformly graded meshes. Our theoretical results are numerically validated in a sample of test problems.

preprint2012arXiv

Superconvergence of a discontinuous Galerkin method for fractional diffusion and wave equations

We consider an initial-boundary value problem for $\partial_tu-\partial_t^{-α}\nabla^2u=f(t)$ , that is, for a fractional diffusion ( $-1<α<0$ ) or wave ( $0<α<1$ ) equation. A numerical solution is found by applying a piecewise-linear, discontinuous Galerkin method in time combined with a piecewise-linear, conforming finite element method in space. The time mesh is graded appropriately near $t=0$ , but the spatial mesh is quasiuniform. Previously, we proved that the error, measured in the spatial $L_2$ -norm, is of order $k^{2+α_-}+h^2\ell(k)$ , uniformly in $t$ , where $k$ is the maximum time step, $h$ is the maximum diameter of the spatial finite elements, $α_-=\min(α,0)\le0$ and $\ell(k)=\max(1,|\log k|)$ . Here, we generalize a known result for the classical heat equation (i.e., the case $α=0$ ) by showing that at each time level $t_n$ the solution is superconvergent with respect to $k$ : the error is of order $(k^{3+2α_-}+h^2)\ell(k)$ . Moreover, a simple postprocessing step employing Lagrange interpolation yields a superconvergent approximation for any $t$ . Numerical experiments indicate that our theoretical error bound is pessimistic if $α<0$ . Ignoring logarithmic factors, we observe that the error in the DG solution at $t=t_n$ , and after postprocessing at all $t$ , is of order $k^{3+α_-}+h^2$ .

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Kassem Mustapha

Research overview

Is this your work?

14 publications

FEM for time-fractional diffusion equations, novel optimal error analyses

Uniform stability for a spatially-discrete, subdiffusive Fokker-Planck equation

Regularity theory for time-fractional advection-diffusion-reaction equations

Well-posedness of time-fractional, advection-diffusion-reaction equations

An inverse source problem for a two-parameter anomalous diffusion with local time datum

Finite volume element method for two-dimensional fractional subdiffusion problems

Numerical solution of the time-fractional Fokker-Planck equation with general forcing

A discontinuous Petrov-Galerkin method for time-fractional diffusion equations

A hybridizable discontinuous Galerkin method for fractional diffusion problems

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

Time-stepping discontinuous Galerkin methods for fractional diffusion problems

Time-stepping error bounds for fractional diffusion problems with non-smooth initial data

A Superconvergent discontinuous Galerkin method for Volterra integro-differential equations, smooth and non-smooth kernels

Superconvergence of a discontinuous Galerkin method for fractional diffusion and wave equations

See the researcher in context

Building this map preview