With the help of the asymptotic expansion for the classic Li formula and based on the L1-type compact difference scheme,we propose a temporal Richardson extrapolation method for the fractional sub-diffusion equation.T...With the help of the asymptotic expansion for the classic Li formula and based on the L1-type compact difference scheme,we propose a temporal Richardson extrapolation method for the fractional sub-diffusion equation.Three extrapolation formulas are presented,whose temporal convergence orders in L_(∞)-norm are proved to be 2,3-α,and 4-2α,respectively,where 0<α<1.Similarly,by the method of order reduction,an extrapola-tion method is constructed for the fractional wave equation including two extrapolation formulas,which achieve temporal 4-γ and 6-2γ order in L_(∞)-norm,respectively,where1<γ<2.Combining the derived extrapolation methods with the fast algorithm for Caputo fractional derivative based on the sum-of-exponential approximation,the fast extrapolation methods are obtained which reduce the computational complexity significantly while keep-ing the accuracy.Several numerical experiments confirm the theoretical results.展开更多
基金supported by the National Natural Science Foundation of China(grant number 11671081).
文摘With the help of the asymptotic expansion for the classic Li formula and based on the L1-type compact difference scheme,we propose a temporal Richardson extrapolation method for the fractional sub-diffusion equation.Three extrapolation formulas are presented,whose temporal convergence orders in L_(∞)-norm are proved to be 2,3-α,and 4-2α,respectively,where 0<α<1.Similarly,by the method of order reduction,an extrapola-tion method is constructed for the fractional wave equation including two extrapolation formulas,which achieve temporal 4-γ and 6-2γ order in L_(∞)-norm,respectively,where1<γ<2.Combining the derived extrapolation methods with the fast algorithm for Caputo fractional derivative based on the sum-of-exponential approximation,the fast extrapolation methods are obtained which reduce the computational complexity significantly while keep-ing the accuracy.Several numerical experiments confirm the theoretical results.