Papers by Kassem Mustapha
arXiv (Cornell University), Aug 7, 2022
Bookmarks Related papers MentionsView impact
Numerical Algorithms
We consider the time discretization of a linear parabolic problem by the discontinuous Galerkin (... more We consider the time discretization of a linear parabolic problem by the discontinuous Galerkin (DG) method using piecewise polynomials of degree at most r − 1 in t, for r ≥ 1 and with maximum step size k. It is well known that the spatial L2-norm of the DG error is of optimal order kr globally in time, and is, for r ≥ 2, superconvergent of order k2r− 1 at the nodes. We show that on the n th subinterval (tn− 1,tn), the dominant term in the DG error is proportional to the local right Radau polynomial of degree r. This error profile implies that the DG error is of order kr+ 1 at the right-hand Gauss–Radau quadrature points in each interval. We show that the norm of the jump in the DG solution at the left end point tn− 1 provides an accurate a posteriori estimate for the maximum error over the subinterval (tn− 1,tn). Furthermore, a simple post-processing step yields a continuous piecewise polynomial of degree r with the optimal global convergence rate of order kr+ 1. We illustrate thes...
Bookmarks Related papers MentionsView impact
IMA Journal of Numerical Analysis
A second-order accurate time-stepping scheme for solving a time-fractional Fokker–Planck equation... more A second-order accurate time-stepping scheme for solving a time-fractional Fokker–Planck equation of order $\alpha \in (0, 1)$, with a general driving force, is investigated. A stability bound for the semidiscrete solution is obtained for $\alpha \in (1/2,1)$ via a novel and concise approach. Our stability estimate is $\alpha $-robust in the sense that it remains valid in the limiting case where $\alpha $ approaches $1$ (when the model reduces to the classical Fokker–Planck equation), a limit that presents practical importance. Concerning the error analysis, we obtain an optimal second-order accurate estimate for $\alpha \in (1/2,1)$. A time-graded mesh is used to compensate for the singular behavior of the continuous solution near the origin. The time-stepping scheme scheme is associated with a standard spatial Galerkin finite element discretization to numerically support our theoretical contributions. We employ the resulting fully discrete computable numerical scheme to perform so...
Bookmarks Related papers MentionsView impact
Computational and Applied Mathematics, 2021
Bookmarks Related papers MentionsView impact
Communications on Applied Mathematics and Computation, 2020
Bookmarks Related papers MentionsView impact
Fractional Calculus and Applied Analysis, 2019
We establish the well-posedness of an initial-boundary value problem for a general class of linea... more We establish the well-posedness of an initial-boundary value problem for a general class of linear 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.
Bookmarks Related papers MentionsView impact
Computers & Mathematics with Applications, 2019
Bookmarks Related papers MentionsView impact
Computational Particle Mechanics, 2019
Bookmarks Related papers MentionsView impact
SIAM Journal on Numerical Analysis, 2016
Bookmarks Related papers MentionsView impact
SIAM Journal on Numerical Analysis, 2014
Bookmarks Related papers MentionsView impact
Numerical Algorithms, 2010
Bookmarks Related papers MentionsView impact
SIAM Journal on Numerical Analysis, 2020
Bookmarks Related papers MentionsView impact
Numerical Algorithms, 2021
Bookmarks Related papers MentionsView impact
Bookmarks Related papers MentionsView impact
K. MUSTAPHA School of Mathematics, University of New South Wales Sydney 2052, Australia kassim©ma... more K. MUSTAPHA School of Mathematics, University of New South Wales Sydney 2052, Australia kassim©maths, unsw. edu. au Abstract--In this work, we analyse the error of a discrete Petrov-Galerkin scheme for nonlin- ear ruth-order ordinary differential and integrodifferential equations on a finite interval subject to nonlinear side conditions. As a trial space we chose high-order Cm-splines. We prove optimal-order convergence and superconvergence in the knots for lower-order derivatives, where the range of deriva- tives for these enhanced convergences to hold is determined by the behaviour of the nonlocal part of the integrodifferential equation. Our results extend and simplify earlier results by Ganesh and Sloan [1]. The numerical experiments in this work, for several singularly perturbed ordinary differential equa- tions, demonstrate the power of our scheme that does not require any mesh restriction. (~) 2005 Elsevier Ltd. All rights reserved. Keywords--Discrete Petrov-Galerkin method, ...
Bookmarks Related papers MentionsView impact
Journal of Scientific Computing, 2021
Bookmarks Related papers MentionsView impact
Arabian Journal of Mathematics, 2020
In this paper, a weakly dissipative viscoelastic plate equation with an infinite memory is consid... more In this paper, a weakly dissipative viscoelastic plate equation with an infinite memory is considered. We show a general energy decay rate for a wide class of relaxation functions. To support our theoretical findings, some numerical illustrations are presented at the end. The numerical solution is computed using the popular finite element method in space, combined with time-stepping finite differences.
Bookmarks Related papers MentionsView impact
Mathematics of Computation, 2017
Bookmarks Related papers MentionsView impact
Journal of Scientific Computing, 2017
Bookmarks Related papers MentionsView impact
Computers & Mathematics with Applications, 2017
Bookmarks Related papers MentionsView impact
Uploads
Papers by Kassem Mustapha