Spectral analysis and structure preserving preconditioners for fractional diffusion equations
Journal of Computational Physics(2016)
摘要
Fractional partial order diffusion equations are a generalization of classical partial differential equations, used to model anomalous diffusion phenomena. When using the implicit Euler formula and the shifted Grünwald formula, it has been shown that the related discretizations lead to a linear system whose coefficient matrix has a Toeplitz-like structure. In this paper we focus our attention on the case of variable diffusion coefficients. Under appropriate conditions, we show that the sequence of the coefficient matrices belongs to the Generalized Locally Toeplitz class and we compute the symbol describing its asymptotic eigenvalue/singular value distribution, as the matrix size diverges. We employ the spectral information for analyzing known methods of preconditioned Krylov and multigrid type, with both positive and negative results and with a look forward to the multidimensional setting. We also propose two new tridiagonal structure preserving preconditioners to solve the resulting linear system, with Krylov methods such as CGNR and GMRES. A number of numerical examples show that our proposal is more effective than recently used circulant preconditioners.
更多查看译文
关键词
Fractional diffusion equations,Toeplitz matrix,Locally Toeplitz sequence of matrices,Singular value/eigenvalue distribution,Preconditioning
AI 理解论文
溯源树
样例
生成溯源树,研究论文发展脉络