Using Gauss-Jacobi quadrature rule to improve the accuracy of FEM for spatial fractional problems

Zongze Yang, Jungang Wang, Zhanbin Yuan, Yufeng Nie

Research output: Contribution to journalArticlepeer-review

2 Scopus citations

Abstract

Though the finite element method has been widely used in solving fractional differential equations, the effects of the Gaussian quadrature rule on the numerical results have rarely been considered. Since the fractional derivatives of the basis functions are not polynomials with integer power and always have weak singularities on some elements, the Gaussian quadrature rule (Gauss-Legendre quadrature rule) may not be suitable in assembling the fractional stiffness matrix. By splitting the integrand of the inner products into a weak singularity part and a smooth part and utilizing the Gauss-Jacobi quadrature rule for the weak singularity part, we present a modified algorithm to assemble the fractional stiffness matrix. The numerical results, conducted on 1D and 2D domains, show that our method can significantly improve the accuracy of the stiffness matrix as well as the accuracy of the numerical solution with much fewer Gaussian points.

Original languageEnglish
Pages (from-to)1389-1411
Number of pages23
JournalNumerical Algorithms
Volume89
Issue number3
DOIs
StatePublished - Mar 2022

Keywords

  • Finite element method
  • Fractional stiffness matrix
  • Gauss-Jacobi quadrature rule
  • Gaussian quadrature rule
  • Riemann-Liouville fractional derivatives

Fingerprint

Dive into the research topics of 'Using Gauss-Jacobi quadrature rule to improve the accuracy of FEM for spatial fractional problems'. Together they form a unique fingerprint.

Cite this