English

Numerical Approximation of Fractional Powers of Regularly Accretive Operators

Numerical Analysis 2016-07-15 v3

Abstract

We study the numerical approximation of fractional powers of accretive operators in this paper. Namely, if AA is the accretive operator associated with an accretive sesquilinear form A(,)A(\cdot,\cdot) defined on a Hilbert space V\mathbb V contained in L2(Ω)L^2(\Omega), we approximate AβA^{-\beta} for β(0,1)\beta\in (0,1). The fractional powers are defined in terms of the so-called Balakrishnan integral formula. Given a finite element approximation space VhV\mathbb V_h\subset \mathbb V, AβA^{-\beta} is approximated by AhβπhA_h^{-\beta}\pi_h where AhA_h is the operator associated with the form A(,)A(\cdot,\cdot) restricted to Vh\mathbb V_h and πh\pi_h is the L2(Ω)L^2(\Omega)-projection onto Vh\mathbb V_h. We first provide error estimates for (AβAhβπh)f(A^\beta-A_h^{\beta}\pi_h)f in Sobolev norms with index in [0,1] for appropriate ff. These results depend on elliptic regularity properties of variational solutions involving the form A(,)A(\cdot,\cdot) and are valid for the case of less than full elliptic regularity. We also construct and analyze an exponentially convergent sinc quadrature approximation to the Balakrishnan integral defining AhβπhfA_h^{\beta}\pi_h f. Finally, the results of numerical computations illustrating the proposed method are given.

Keywords

Cite

@article{arxiv.1508.05869,
  title  = {Numerical Approximation of Fractional Powers of Regularly Accretive Operators},
  author = {Andrea Bonito and Joseph E. Pasciak},
  journal= {arXiv preprint arXiv:1508.05869},
  year   = {2016}
}

Comments

25 pages