A Mixed FEM for a Time-Fractional Fokker–Planck Model

Samir Karaa, Kassem Mustapha*, Naveed Ahmed

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review


We propose and analyze a mixed finite element method for the spatial approximation of a time-fractional Fokker–Planck equation in a convex polyhedral domain, where the given driving force is a function of space. Taking into account the limited smoothing properties of the model, and considering an appropriate splitting of the errors, we employed a sequence of clever energy arguments to show optimal convergence rates with respect to both approximation properties and regularity results. In particular, error bounds for both primary and secondary variables are derived in L2-norm for cases with smooth and nonsmooth initial data. We further investigate a fully implicit time-stepping scheme based on a convolution quadrature in time generated by the backward Euler method. Our main result provides pointwise-in-time optimal L2-error estimates for the primary variable. Numerical examples are then presented to illustrate the theoretical contributions.

Original languageEnglish
Article number59
JournalJournal of Scientific Computing
Issue number3
Publication statusPublished - Jan 2024


  • 35S10
  • 65M12
  • 65M15
  • 65M60
  • 65M70
  • Convolution quadrature
  • Error analysis
  • Fractional Fokker–Planck equation
  • Mixed finite elements
  • Smooth and nonsmooth initial data

ASJC Scopus subject areas

  • Software
  • Theoretical Computer Science
  • Numerical Analysis
  • General Engineering
  • Computational Mathematics
  • Computational Theory and Mathematics
  • Applied Mathematics

Cite this