Abstract
In this work, we propose a Crank-Nicolson-type scheme with variable steps for the time fractional Allen-Cahn equation. The proposed scheme is shown to be unconditionally stable (in a variational energy sense) and is maximum bound preserving. Interestingly, the discrete energy stability result obtained in this paper can recover the classical energy dissipation law when the fractional order α \rightarrow 1. That is, our scheme can asymptotically preserve the energy dissipation law in the α \rightarrow 1 limit. This seems to be the first work on a variable time-stepping scheme that can preserve both the energy stability and the maximum bound principle. Our Crank-Nicolson scheme is built upon a reformulated problem associated with the Riemann-Liouville derivative. As a byproduct, we build up a reversible transformation between the L1-type formula of the Riemann-Liouville derivative and a new L1-type formula of the Caputo derivative with the help of a class of discrete orthogonal convolution kernels. This is the first time such a discrete transformation is established between two discrete fractional derivatives. We finally present several numerical examples with an adaptive time-stepping strategy to show the effectiveness of the proposed scheme.
Original language | English |
---|---|
Pages (from-to) | A3503-A3526 |
Number of pages | 24 |
Journal | SIAM Journal on Scientific Computing |
Volume | 43 |
Issue number | 5 |
DOIs | |
Publication status | Published - 2021 |
Scopus Subject Areas
- Computational Mathematics
- Applied Mathematics
User-Defined Keywords
- adaptive time stepping
- asymptotic preserving
- energy stability
- maximum principle
- time-fractional Allen-Cahn equation