14 research outputs found
Numerical solving unsteady space-fractional problems with the square root of an elliptic operator
An unsteady problem is considered for a space-fractional equation in a
bounded domain. A first-order evolutionary equation involves the square root of
an elliptic operator of second order. Finite element approximation in space is
employed. To construct approximation in time, regularized two-level schemes are
used. The numerical implementation is based on solving the equation with the
square root of the elliptic operator using an auxiliary Cauchy problem for a
pseudo-parabolic equation. The scheme of the second-order accuracy in time is
based on a regularization of the three-level explicit Adams scheme. More
general problems for the equation with convective terms are considered, too.
The results of numerical experiments are presented for a model two-dimensional
problem.Comment: 21 pages, 7 figures. arXiv admin note: substantial text overlap with
arXiv:1412.570
Positive approximations of the inverse of fractional powers of SPD M-matrices
This study is motivated by the recent development in the fractional calculus
and its applications. During last few years, several different techniques are
proposed to localize the nonlocal fractional diffusion operator. They are based
on transformation of the original problem to a local elliptic or
pseudoparabolic problem, or to an integral representation of the solution, thus
increasing the dimension of the computational domain. More recently, an
alternative approach aimed at reducing the computational complexity was
developed. The linear algebraic system , is considered, where is a properly normalized (scalded) symmetric
and positive definite matrix obtained from finite element or finite difference
approximation of second order elliptic problems in ,
. The method is based on best uniform rational approximations (BURA)
of the function for and natural .
The maximum principles are among the major qualitative properties of linear
elliptic operators/PDEs. In many studies and applications, it is important that
such properties are preserved by the selected numerical solution method. In
this paper we present and analyze the properties of positive approximations of
obtained by the BURA technique. Sufficient conditions for
positiveness are proven, complemented by sharp error estimates. The theoretical
results are supported by representative numerical tests
A Gauss-Laguerre approach for the resolvent of fractional powers
This paper introduces a very fast method for the computation of the resolvent of fractional powers of
operators. The analysis is kept in the continuous setting of (potentially unbounded) self-adjoint positive operators in Hilbert spaces. The method is based on the Gauss-Laguerre rule, exploiting a particular integral representation of the resolvent. We provide sharp error estimates that can be used to a priori select the number of nodes to achieve a prescribed tolerance