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.

No reviews yet.

No comments yet.