We kindly invite you to the Miklós Farkas Seminar
24 October (Thursday) 10.15, BME, H306
Svetozar Margenov (Institute of Information and Communication Technologies, Bulgarian Academy of Sciences)
Numerical Solution of Fractional Diffusion Problems
This study is motivated by the recent achievements in fractional calculus and its numerous applications related to anomalous (super) diffusion. Let us consider a fractional power of a self-adjoint elliptic operator introduced through its spectral decomposition. It is self-adjoint but nonlocal. The nonlocal problems are computationally expensive. Several different techniques were recently proposed to localize the nonlocal elliptic operator, thus increasing the space dimension of the original computational domain.
An alternative approach [1, 2, 3] is discussed in this talk. Let A be a properly scaled symmetric and positive definite (SPD) sparse matrix, arising from finite element or finite difference discretization of the initial (standard, local) diffusion problem. A method for solving algebraic systems of linear equations involving Aα, 0<α<1, is presented. The solution methods are based on best uniform rational approximations (BURA) of the scalar function tα, 0≤t≤1. The method has exponential convergence rate with respect to the degree of rational approximation k. The error estimates of the last variant of BURA methods are robust with respect to the spectral condition number κ(A). A stabilized modification of the Remez algorithm is developed to compute the BURA of tα. Although the fractional power of A is a dense matrix, the algorithm has complexity of order O(N), where N is the number of unknowns. At this point we assume that some solver of optimal complexity (say multigrid or multilevel) is used for the involved systems with matrices A+djI, dj≥0, j=1,…,k. The comparative numerical tests confirm the advantages of the BURA method.
Acknowledgement: This research has been partially supported by the Bulgarian NSF Grant DN12/1.
(István Faragó, János Karátson, Róbert Horváth, Miklós Mincsovics, Gabriella Svantnerné Sebestyén)
Visit the homepage of the seminar: http://math.bme.hu/AlkAnalSzemi