TY - JOUR
T1 - Using Gauss-Jacobi quadrature rule to improve the accuracy of FEM for spatial fractional problems
AU - Yang, Zongze
AU - Wang, Jungang
AU - Yuan, Zhanbin
AU - Nie, Yufeng
N1 - Publisher Copyright:
© 2021, The Author(s), under exclusive licence to Springer Science+Business Media, LLC, part of Springer Nature.
PY - 2022/3
Y1 - 2022/3
N2 - 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.
AB - 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.
KW - Finite element method
KW - Fractional stiffness matrix
KW - Gauss-Jacobi quadrature rule
KW - Gaussian quadrature rule
KW - Riemann-Liouville fractional derivatives
UR - http://www.scopus.com/inward/record.url?scp=85110551912&partnerID=8YFLogxK
U2 - 10.1007/s11075-021-01158-x
DO - 10.1007/s11075-021-01158-x
M3 - 文章
AN - SCOPUS:85110551912
SN - 1017-1398
VL - 89
SP - 1389
EP - 1411
JO - Numerical Algorithms
JF - Numerical Algorithms
IS - 3
ER -