Abstract
In this paper, a fast finite volume method is proposed for the initial and boundary value problems of spatial fractional diffusion equations on nonuniform meshes. The discretizations of the Riemann-Liouville fractional derivatives lead to unstructured dense coefficient matrices, differing from the Toeplitz-like structure under the uniform mesh. The fast algorithm is proposed by using the sum-of-exponentials (SOE) technique to the spatial kernel xα−1,α∈(0,1). Then, the matrix-vector multiplications of the resulting coefficient matrices could be implemented in O(mlog2m) operations, where m denotes the size of matrices. Iterative solvers are preferably applied to obtain the numerical solution. The proposed fast scheme is proved to be unconditionally stable for sufficiently accurate SOE approximation. Meanwhile, a banded preconditioner is exploited to accelerate the Krylov subspace method. Numerical experiments are provided to demonstrate the efficiency of the proposed fast algorithm.
Talk to us
Join us for a 30 min session where you can share your feedback and ask us any queries you have
Disclaimer: All third-party content on this website/platform is and will remain the property of their respective owners and is provided on "as is" basis without any warranties, express or implied. Use of third-party content does not indicate any affiliation, sponsorship with or endorsement by them. Any references to third-party content is to identify the corresponding services and shall be considered fair use under The CopyrightLaw.