scholarly journals Rational Krylov methods for fractional diffusion problems on graphs

Author(s):  
Michele Benzi ◽  
Igor Simunec

AbstractIn this paper we propose a method to compute the solution to the fractional diffusion equation on directed networks, which can be expressed in terms of the graph Laplacian L as a product $$f(L^T) \varvec{b}$$ f ( L T ) b , where f is a non-analytic function involving fractional powers and $$\varvec{b}$$ b is a given vector. The graph Laplacian is a singular matrix, causing Krylov methods for $$f(L^T) \varvec{b}$$ f ( L T ) b to converge more slowly. In order to overcome this difficulty and achieve faster convergence, we use rational Krylov methods applied to a desingularized version of the graph Laplacian, obtained with either a rank-one shift or a projection on a subspace.

2021 ◽  
Vol 0 (0) ◽  
Author(s):  
Tobias Danczul ◽  
Clemens Hofreither

Abstract We establish an equivalence between two classes of methods for solving fractional diffusion problems, namely, Reduced Basis Methods (RBM) and Rational Krylov Methods (RKM). In particular, we demonstrate that several recently proposed RBMs for fractional diffusion can be interpreted as RKMs. This changed point of view allows us to give convergence proofs for some methods where none were previously available. We also propose a new RKM for fractional diffusion problems with poles chosen using the best rational approximation of the function 𝑧 −𝑠 with 𝑧 ranging over the spectral interval of the spatial discretization matrix. We prove convergence rates for this method and demonstrate numerically that it is competitive with or superior to many methods from the reduced basis, rational Krylov, and direct rational approximation classes. We provide numerical tests for some elliptic fractional diffusion model problems.


2020 ◽  
Vol 28 (3) ◽  
pp. 147-160
Author(s):  
Andrea Bonito ◽  
Diane Guignard ◽  
Ashley R. Zhang

AbstractWe consider the numerical approximation of the spectral fractional diffusion problem based on the so called Balakrishnan representation. The latter consists of an improper integral approximated via quadratures. At each quadrature point, a reaction–diffusion problem must be approximated and is the method bottle neck. In this work, we propose to reduce the computational cost using a reduced basis strategy allowing for a fast evaluation of the reaction–diffusion problems. The reduced basis does not depend on the fractional power s for 0 < smin ⩽ s ⩽ smax < 1. It is built offline once for all and used online irrespectively of the fractional power. We analyze the reduced basis strategy and show its exponential convergence. The analytical results are illustrated with insightful numerical experiments.


2021 ◽  
Vol 2021 (1) ◽  
Author(s):  
Nguyen Hoang Tuan ◽  
Nguyen Anh Triet ◽  
Nguyen Hoang Luc ◽  
Nguyen Duc Phuong

AbstractIn this work, we consider a fractional diffusion equation with nonlocal integral condition. We give a form of the mild solution under the expression of Fourier series which contains some Mittag-Leffler functions. We present two new results. Firstly, we show the well-posedness and regularity for our problem. Secondly, we show the ill-posedness of our problem in the sense of Hadamard. Using the Fourier truncation method, we construct a regularized solution and present the convergence rate between the regularized and exact solutions.


2020 ◽  
Vol 2020 (1) ◽  
Author(s):  
Abderrazak Nabti ◽  
Ahmed Alsaedi ◽  
Mokhtar Kirane ◽  
Bashir Ahmad

Abstract We prove the nonexistence of solutions of the fractional diffusion equation with time-space nonlocal source $$\begin{aligned} u_{t} + (-\Delta )^{\frac{\beta }{2}} u =\bigl(1+ \vert x \vert \bigr)^{ \gamma } \int _{0}^{t} (t-s)^{\alpha -1} \vert u \vert ^{p} \bigl\Vert \nu ^{ \frac{1}{q}}(x) u \bigr\Vert _{q}^{r} \,ds \end{aligned}$$ u t + ( − Δ ) β 2 u = ( 1 + | x | ) γ ∫ 0 t ( t − s ) α − 1 | u | p ∥ ν 1 q ( x ) u ∥ q r d s for $(x,t) \in \mathbb{R}^{N}\times (0,\infty )$ ( x , t ) ∈ R N × ( 0 , ∞ ) with initial data $u(x,0)=u_{0}(x) \in L^{1}_{\mathrm{loc}}(\mathbb{R}^{N})$ u ( x , 0 ) = u 0 ( x ) ∈ L loc 1 ( R N ) , where $p,q,r>1$ p , q , r > 1 , $q(p+r)>q+r$ q ( p + r ) > q + r , $0<\gamma \leq 2 $ 0 < γ ≤ 2 , $0<\alpha <1$ 0 < α < 1 , $0<\beta \leq 2$ 0 < β ≤ 2 , $(-\Delta )^{\frac{\beta }{2}}$ ( − Δ ) β 2 stands for the fractional Laplacian operator of order β, the weight function $\nu (x)$ ν ( x ) is positive and singular at the origin, and $\Vert \cdot \Vert _{q}$ ∥ ⋅ ∥ q is the norm of $L^{q}$ L q space.


Sign in / Sign up

Export Citation Format

Share Document