Neumann fractional diffusion problems: BURA solution methods and algorithms

Stanislav Harizanov, Nikola Kosturski, Svetozar Margenov and Yavor Vutov

Mathematics and Computers in Simulation (MATCOM), 2021, vol. 189, issue C, 85-98

Abstract: Let us consider the non-local problem −Lαu=f, α∈(0,1), L is a second order self-adjoint elliptic operator in Ω⊂Rd with Neumann boundary conditions on ∂Ω. The problem is discretized by finite difference or finite element method, thus obtaining the linear system Aαu=f, A is sparse symmetric and positive semidefinite matrix. The proposed method is based on best uniform rational approximations (BURA) of degree k, rα,k, of the scalar function tα, t∈[0,1]. Then, the approximate solution ur of the fractional power linear system is defined as u≈ur=λ2−αrα,k(λ2A†)f, where λ2 is the first positive eigenvalue of A, and A† stands for the Moore–Penrose pseudo inverse of A. The BURA method reduces the non-local problem to solution of k linear systems with matrices A+diI, di>0, i=1,…,k. An exponential convergence rate with respect to k is proven. The error estimates are robust with respect to the condition number of A in the subspace orthogonal to the constant vectors. The algorithm has almost optimal computational complexity assuming that optimal iterative solvers are applied to the auxiliary sparse linear systems. The first group of numerical tests illustrates in detail the obtained theoretical results. Finite difference discretization of a model 2D problem is used for this purpose. The second part of numerical tests demonstrates the applicability of BURA methods for 3D problems in domains with general geometry. Linear finite elements on unstructured tetrahedral meshes with local mesh refinement are used in the presented large-scale experiments, confirming also the almost optimal computational complexity.

Keywords: BURA; Fractional diffusion; Neumann boundary value problems; unstructured meshes; computational complexity (search for similar items in EconPapers)
Date: 2021
References: View references in EconPapers View complete reference list from CitEc
Citations: View citations in EconPapers (1)

Downloads: (external link)
http://www.sciencedirect.com/science/article/pii/S0378475420302469
Full text for ScienceDirect subscribers only

Related works:
This item may be available elsewhere in EconPapers: Search for items with the same title.

Export reference: BibTeX RIS (EndNote, ProCite, RefMan) HTML/Text

Persistent link: https://EconPapers.repec.org/RePEc:eee:matcom:v:189:y:2021:i:c:p:85-98

DOI: 10.1016/j.matcom.2020.07.018

Access Statistics for this article

Mathematics and Computers in Simulation (MATCOM) is currently edited by Robert Beauwens

More articles in Mathematics and Computers in Simulation (MATCOM) from Elsevier
Bibliographic data for series maintained by Catherine Liu ().