20.8NAMay 23
Spectral analysis and sine transform based preconditioning for a structure preserving stabilized scheme approximating the space-fractional Allen Cahn equation with logarithmic potentialDanyal Ahmad, Stefano Serra-Capizzano, Muhammad Sohaib et al.
We consider an initial boundary value problem of the space fractional Allen-Cahn equation with logarithmic Flory-Huggins potential. As an approximation technique, first-order weighted and shifted Grunwald difference formulae of the left and right fractional derivatives are used. The main focus of the present work is to study the spectral features of the underlying matrices and matrix sequences and to design proper preconditioners based on the spectral information. Then a computational and spectral analysis of the resulting preconditioned matrix sequences is performed. Numerical evidence and a short list of open problems complete the current study.
31.5NAMar 16
Structure-preserving preconditioning of discrete space-fractional diffusion equations with variable coefficient and θ-MethodMuhammad Faisal Khan, Asim Ilyas, Rolf Krause et al.
This paper studies the spectral properties of large matrices and the preconditioning of linear systems, arising from the finite difference discretization of a time-dependent space-fractional diffusion equation with a variable coefficient $a(x)$ defined on $Ω\subset \mathbb{R}^d$, $d=1,2$. The model involves a one-sided Riemann-Liouville fractional derivative multiplied by the function $a(x)$, discretized by the shifted Gr"unwald formula in space and the $θ$-method in time. The resulting all-at-once linear systems exhibit a $(d+1)$-level Toeplitz-like matrix structure, with $d=1,2$ denoting the space dimension, while the additional level is due to the time variable. A preconditioning strategy is developed based on the structural properties of the discretized operator. Using the generalized locally Toeplitz (GLT) theory, we analyze the spectral distribution of the unpreconditioned and preconditioned matrix sequences. The main novelty is that the analysis fully covers the case where the variable coefficient $a$ is nonconstant. Numerical results are provided to support the GLT based theoretical findings, and some possible extensions are briefly discussed.