Source author record

Joseph E. Pasciak

Explore this author’s publications, collaborators and research topics. This catalog record has not been claimed by its author.

ResearcherUnclaimed source record

math.NA math.AP

Catalog footprint

Research overview

3works

2topics

2close collaborators

Actions

Is this your work?

Open graph Browse papers

preprint2016arXiv

Numerical Approximation of Fractional Powers of Regularly Accretive Operators

We study the numerical approximation of fractional powers of accretive operators in this paper. Namely, if $A$ is the accretive operator associated with an accretive sesquilinear form $A(\cdot,\cdot)$ defined on a Hilbert space $\mathbb V$ contained in $L^2(Ω)$ , we approximate $A^{-β}$ for $β\in (0,1)$ . The fractional powers are defined in terms of the so-called Balakrishnan integral formula. Given a finite element approximation space $\mathbb V_h\subset \mathbb V$ , $A^{-β}$ is approximated by $A_h^{-β}π_h$ where $A_h$ is the operator associated with the form $A(\cdot,\cdot)$ restricted to $\mathbb V_h$ and $π_h$ is the $L^2(Ω)$ -projection onto $\mathbb V_h$ . We first provide error estimates for $(A^β-A_h^βπ_h)f$ in Sobolev norms with index in [0,1] for appropriate $f$ . These results depend on elliptic regularity properties of variational solutions involving the form $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 $A_h^βπ_h f$ . Finally, the results of numerical computations illustrating the proposed method are given.

preprint2016arXiv

The Approximation of Parabolic Equations Involving Fractional Powers of Elliptic Operators

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 $H^1_0(Ω)$ . The time dependent solution $u(x,t)$ is represented as a Dunford Taylor 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 $(z_qI-L)^{-1}v\in H^1_0(Ω)$ for a finite number of quadrature points $z_q$ lying along the contour. In turn, these quantities are approximated using complex valued continuous piecewise linear finite elements. Our main result provides $L^2(Ω)$ error estimates between the solution $u(\cdot,t)$ and its final approximation. Numerical results illustrating the behavior of the algorithms are provided.

preprint2013arXiv

Numerical Approximation of Fractional Powers of Elliptic Operators

We present and study a novel numerical algorithm to approximate the action of $T^β:=L^{-β}$ where $L$ is a symmetric and positive definite unbounded operator on a Hilbert space $H_0$ . The numerical method is based on a representation formula for $T^{-β}$ in terms of Bochner integrals involving $(I+t^2L)^{-1}$ for $t\in(0,\infty)$ . To develop an approximation to $T^β$ , we introduce a finite element approximation $L_h$ to $L$ and base our approximation to $T^β$ on $T_h^β:= L_h^{-β}$ . The direct evaluation of $T_h^β$ is extremely expensive as it involves expansion in the basis of eigenfunctions for $L_h$ . The above mentioned representation formula holds for $T_h^{-β}$ and we propose three quadrature approximations denoted generically by $Q_h^β$ . The two results of this paper bound the errors in the $H_0$ inner product of $T^β-T_h^βπ_h$ and $T_h^β-Q_h^β$ where $π_h$ is the $H_0$ orthogonal projection into the finite element space. We note that the evaluation of $Q_h^β$ involves application of $(I+(t_i)^2L_h)^{-1}$ with $t_i$ being either a quadrature point or its inverse. Efficient solution algorithms for these problems are available and the problems at different quadrature points can be straightforwardly solved in parallel. Numerical experiments illustrating the theoretical estimates are provided for both the quadrature error $T_h^β-Q_h^β$ and the finite element error $T^β-T_h^βπ_h$ .

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.