Symmetric interior penalty Galerkin method for fractional-in-space phase-field equations

Stoll M., YÜCEL H.

AIMS MATHEMATICS, vol.3, no.1, pp.66-95, 2018 (Journal Indexed in ESCI) identifier identifier

  • Publication Type: Article / Article
  • Volume: 3 Issue: 1
  • Publication Date: 2018
  • Doi Number: 10.3934/math.2018.1.66
  • Title of Journal : AIMS MATHEMATICS
  • Page Numbers: pp.66-95


Fractional differential equations are becoming increasingly popular as a modelling tool to describe a wide range of non-classical phenomena with spatial heterogeneities throughout the applied sciences and engineering. However, the non-local nature of the fractional operators causes essential difficulties and challenges for numerical approximations. We here investigate the numerical solution of fractional-in-space phase-field models such as Allen-Cahn and Cahn-Hilliard equations via the contour integral method (CIM) for computing the fractional power of a matrix times a vector. Time discretization is performed by the first-and second-order implicit-explicit schemes with an adaptive time-step size approach, whereas spatial discretization is performed by a symmetric interior penalty Galerkin (SIPG) method. Several numerical examples are presented to illustrate the effect of the fractional power.