Search

The aim of this paper is to develop and analyze high-order time stepping schemes for approximately solving semilinear subdiffusion equations. We apply the convolution quadrature generated by k-step backward differentiation formula (BDFk) to discretize the time-fractional derivative with order α ∈ (0, 1) and modify the starting steps in order to achieve optimal convergence rate. This method has already been well-studied for the linear fractional evolution equations in Jin, Li, and Zhou [SIAM J. Sci. Comput., 39 (2017), pp. A3129-A3152], while the numerical analysis for the nonlinear problem is still missing in the literature. By splitting the nonlinear potential term into an irregular linear part and a smoother nonlinear part and using the generating function technique, we prove that the convergence order of the corrected BDFk scheme is O(τmin(k,1+2α-ε)) without imposing further assumption on the regularity of the solution. Numerical examples are provided to support our theoretical results. [ABSTRACT FROM AUTHOR]

An implicit-explicit multistep method based on the backward difference formulae (BDF) is proposed for time discretization of parabolic equations with a non-self-adjoint operator. Implicit and explicit schemes are used for the self-adjoint and anti-self-adjoint parts of the operator, respectively. For a k-step method, some correction terms are added to the starting k - 1 steps to maintain kth-order convergence without imposing further compatibility conditions at the initial time. Long-time kth-order convergence for the numerical method is proved under the assumptions that the operator is coercive and that the non-self-adjoint part is low order. Such an operator often appears in practical computation (such as the Stokes--Darcy system) but may violate the standard sectorial angle condition used in the literature for analysis of BDF. In particular, the proposed method and analysis in this paper extend the long-time energy error analysis of the Stokes--Darcy system in Chen et al. [SIAM J. Numer. Anal., 51 (2013), pp. 2563--2584; Numer. Math., 134 (2016), pp. 857--879] to general symmetrized and decoupled BDF methods up to order 6 by using the generating function technique. [ABSTRACT FROM AUTHOR]