Joseph Pasciak

Bangti Jin, Raytcho Lazarov, Joseph Pasciak et al.

We investigate semi-discrete numerical schemes based on the standard Galerkin and lumped mass Galerkin finite element methods for an initial-boundary value problem for homogeneous fractional diffusion problems with non-smooth initial data. We assume that $Ω\subset \mathbb{R}^d$, $d=1,2,3$ is a convex polygonal (polyhedral) domain. We theoretically justify optimal order error estimates in $L_2$- and $H^1$-norms for initial data in $H^{-s}(Ω),~0\le s \le 1$. We confirm our theoretical findings with a number of numerical tests that include initial data $v$ being a Dirac $δ$-function supported on a $(d-1)$-dimensional manifold.

Stanislav Harizanov, Raytcho Lazarov, Pencho Marinov et al.

We discuss, study, and compare experimentally three methods for solving the system of algebraic equations $\mathbb{A}^α\bf{u}=\bf{f}$, $0< α<1$, where $\mathbb{A}$ is a symmetric and positive definite matrix obtained from finite difference or finite element approximations of second order elliptic problems in $\mathbb{R}^d$, $d=1,2,3$. The first method, introduced by Harizanov et.al, based on the best uniform rational approximation (BURA) $r_α(t)$ of $t^{1-α}$ for $0 \le t \le 1$, is used to get the rational approximation of $t^{-α}$ in the form $t^{-1}r_α(t)$. Here we develop another method, denoted by R-BURA, that is based on the best rational approximation $r_{1-α}(t)$ of $t^α$ on the interval $[0,1]$ and approximates $t^{-α}$ via $r^{-1}_{1-α}(t)$. The third method, introduced and studied by Bonito and Pasciak, is based on an exponentially convergent quadrature scheme for the Dundord-Taylor integral representation of the fractional powers of elliptic operators. All three methods reduce the solution of the system $\mathbb{A}^α\bf{u}=\bf{f}$ to solving a number of equations of the type $(\mathbb{A} +c\mathbb{I})\bf{u}= \bf{f}$, $c \ge 0$. Comprehensive numerical experiments on model problems with $\mathbb A$ obtained by approximation of elliptic equations in one and two spatial dimensions are used to compare the efficiency of these three algorithms depending on the fractional power $α$. The presented results prove the concept of the new R-BURA method, which performs well for $α$ close to $1$ in contrast to BURA, which performs well for $α$ close to $0$. As a result, we show theoretically and experimentally, that they have mutually complementary advantages.

Beiping Duan, Bangti Jin, Raytcho Lazarov et al.

We present and analyze a space-time Petrov-Galerkin finite element method for a time-fractional diffusion equation involving a Riemann-Liouville fractional derivative of order $α\in(0,1)$ in time and zero initial data. We derive a proper weak formulation involving different solution and test spaces and show the inf-sup condition for the bilinear form and thus its well-posedness. Further, we develop a novel finite element formulation, show the well-posedness of the discrete problem, and establish error bounds in both energy and $L^2$ norms for the finite element solution. In the proof of the discrete inf-sup condition, a certain nonstandard $L^2$ stability property of the $L^2$ projection operator plays a key role. We provide extensive numerical examples to verify the convergence of the method.