Asymptotically compatible energy of variable-step fractional BDF2 formula for time-fractional Cahn-Hilliard model

Abstract

A new discrete energy dissipation law of the variable-step fractional BDF2 (second-order backward differentiation formula) scheme is established for time-fractional Cahn-Hilliard model with the Caputo's fractional derivative of order $\alpha \in (0,1)$, under a weak step-ratio constraint $0.4753\le {\tau }_k/{\tau }_{k-1}<r^{\ast }(\alpha )$, where ${\tau }_k$ is the $k$-th time-step size and $r^{\ast }(\alpha )\ge 4.660$ for $\alpha \in (0,1)$.We propose a novel discrete gradient structure by a local-nonlocal splitting technique, that is, the fractional BDF2 formula is split into a local part analogue to the two-step backward differentiation formula of the first derivative and a nonlocal part analogue to the L1-type formula of the Caputo's derivative. More interestingly, in the sense of the limit $\alpha \to 1^{-}$, the discrete energy and the corresponding energy dissipation law are asymptotically compatible with the associated discrete energy and the energy dissipation law of the variable-step BDF2 method for the classical Cahn-Hilliard equation, respectively. Numerical examples with an adaptive stepping procedure are provided to demonstrate the accuracy and the effectiveness of our proposed method.

21 pages, 22 figures

References in corpus (7)

• Mesh-robustness of an energy stable BDF2 scheme with variable steps for the Cahn-Hilliard model
• Stability of variable-step BDF2 and BDF3 methods
• The variable-step L1 scheme preserving a compatible energy law for time-fractional Allen-Cahn equation
• Energy Stable L2 Schemes for Time-Fractional Phase-Field Equations
• Time-fractional Cahn-Hilliard equation: Well-posedness, degeneracy, and numerical solutions
• A decreasing upper bound of energy for time-fractional phase-field equations
• Energy stability of variable-step L1-type schemes for time-fractional Cahn-Hilliard model