Numer. Math. Theor. Meth. Appl., 15 (2022), pp. 876-902.
Published online: 2022-10
Cited by
- BibTex
- RIS
- TXT
We develop an efficient and accurate spectral deferred correction (SDC) method for fractional differential equations (FDEs) by extending the algorithm in [14] for classical ordinary differential equations (ODEs). Specifically, we discretize the resulted Picard integral equation by the SDC method and accelerate the convergence of the SDC iteration by using the generalized minimal residual algorithm (GMRES). We first derive the correction matrix of the SDC method for FDEs and analyze the convergence region of the SDC method. We then present several numerical examples for stiff and non-stiff FDEs including fractional linear and nonlinear ODEs as well as fractional phase field models, demonstrating that the accelerated SDC method is much more efficient than the original SDC method, especially for stiff problems. Furthermore, we resolve the issue of low accuracy arising from the singularity of the solutions by using a geometric mesh, leading to highly accurate solutions compared to uniform mesh solutions at almost the same computational cost. Moreover, for long-time integration of FDEs, using the geometric mesh leads to great computational savings as the total number of degrees of freedom required is relatively small.
}, issn = {2079-7338}, doi = {https://doi.org/10.4208/nmtma.OA-2022-0012s}, url = {http://global-sci.org/intro/article_detail/nmtma/21084.html} }We develop an efficient and accurate spectral deferred correction (SDC) method for fractional differential equations (FDEs) by extending the algorithm in [14] for classical ordinary differential equations (ODEs). Specifically, we discretize the resulted Picard integral equation by the SDC method and accelerate the convergence of the SDC iteration by using the generalized minimal residual algorithm (GMRES). We first derive the correction matrix of the SDC method for FDEs and analyze the convergence region of the SDC method. We then present several numerical examples for stiff and non-stiff FDEs including fractional linear and nonlinear ODEs as well as fractional phase field models, demonstrating that the accelerated SDC method is much more efficient than the original SDC method, especially for stiff problems. Furthermore, we resolve the issue of low accuracy arising from the singularity of the solutions by using a geometric mesh, leading to highly accurate solutions compared to uniform mesh solutions at almost the same computational cost. Moreover, for long-time integration of FDEs, using the geometric mesh leads to great computational savings as the total number of degrees of freedom required is relatively small.