Second-order LOD multigrid method for multidimensional Riesz fractional diffusion equation

From MaRDI portal




Abstract: We propose a locally one dimensional (LOD) finite difference method for multidimensional Riesz fractional diffusion equation with variable coefficients on a finite domain. The numerical method is second-order convergent in both space and time directions, and its unconditional stability is strictly proved. Comparing with the popular first-order finite difference method for fractional operator, the form of obtained matrix algebraic equation is changed from (IA)uk+1=uk+bk+1 to (IwidetildeA)uk+1=(I+widetildeB)uk+ildebk+1/2; the three matrices A, widetildeA and widetildeB are all Toeplitz-like, i.e., they have completely same structure and the computational count for matrix vector multiplication is mathcalO(NlogN); and the computational costs for solving the two matrix algebraic equations are almost the same. The LOD-multigrid method is used to solve the resulting matrix algebraic equation, and the computational count is mathcalO(NlogN) and the required storage is mathcalO(N), where N is the number of grid points. Finally, the extensive numerical experiments are performed to show the powerfulness of the second-order scheme and the LOD-multigrid method.



Cites work


Cited in
(23)


Describes a project that uses

Uses Software





This page was built for publication: Second-order LOD multigrid method for multidimensional Riesz fractional diffusion equation

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q466801)