The approximation of parabolic equations involving fractional powers of elliptic operators Academic Article uri icon

abstract

  • 2016 Elsevier B.V. We study the numerical approximation of a time dependent equation involving fractional powers of an elliptic operator L defined to be the unbounded operator associated with a Hermitian, coercive and bounded sesquilinear form on H01(). The time dependent solution u(x,t) is represented as a DunfordTaylor integral along a contour in the complex plane. The contour integrals are approximated using sinc quadratures. In the case of homogeneous right-hand-sides and initial value v, the approximation results in a linear combination of functions (zqIL)1vH01() for a finite number of quadrature points zq lying along the contour. In turn, these quantities are approximated using complex valued continuous piecewise linear finite elements. Our main result provides L2() error estimates between the solution u(,t) and its final approximation. Numerical results illustrating the behavior of the algorithms are provided.

published proceedings

  • JOURNAL OF COMPUTATIONAL AND APPLIED MATHEMATICS

author list (cited authors)

  • Bonito, A., Lei, W., & Pasciak, J. E.

citation count

  • 19

complete list of authors

  • Bonito, Andrea||Lei, Wenyu||Pasciak, Joseph E

publication date

  • January 2017