Paper detail

Optimal long-time decay rate of solutions of complete monotonicity-preserving schemes for nonlinear time-fractional evolutionary equations

The solution of the nonlinear initial-value problem $\mathcal{D}_{t}^αy(t)=-λy(t)^γ$ for $t>0$ with $y(0)>0$, where $\mathcal{D}_{t}^α$ is a Caputo derivative of order $α\in (0,1)$ and $λ, γ$ are positive parameters, is known to exhibit $O(t^{α/γ})$ decay as $t\to\infty$. No corresponding result for any discretisation of this problem has previously been proved. In the present paper it is shown that for the class of complete monotonicity-preserving ($\mathcal{CM}$-preserving) schemes (which includes the L1 and Grünwald-Letnikov schemes) on uniform meshes $\{t_n:=nh\}_{n=0}^\infty$, the discrete solution also has $O(t_{n}^{-α/γ})$ decay as $t_{n}\to\infty$. This result is then extended to $\mathcal{CM}$-preserving discretisations of certain time-fractional nonlinear subdiffusion problems such as the time-fractional porous media and $p$-Laplace equations. For the L1 scheme, the $O(t_{n}^{-α/γ})$ decay result is shown to remain valid on a very general class of nonuniform meshes. Our analysis uses a discrete comparison principle with discrete subsolutions and supersolutions that are carefully constructed to give tight bounds on the discrete solution. Numerical experiments are provided to confirm our theoretical analysis.