Abstract:We consider a family of variable time-stepping Dahlquist-Liniger-Nevanlinna (DLN) schemes, which is unconditional non-linear stable and second order accurate, for the Allen-Cahn equation. The finite element methods are used for the spatial discretization. For the non-linear term, we combine the DLN scheme with two efficient temporal algorithms: partially implicit modified algorithm and scalar auxiliary variable algorithm. For both approaches, we prove the unconditional, long-term stability of the model energy under any arbitrary time step sequence. Moreover, we provide rigorous error analysis for the partially implicit modified algorithm with variable time-stepping. Efficient time adaptive algorithms based on these schemes are also proposed. Several one- and two-dimensional numerical tests are presented to verify the properties of the proposed time adaptive DLN methods.
What problem does this paper attempt to address?
### What problems does this paper attempt to solve?
This paper aims to solve several key problems in the numerical simulation of the Allen - Cahn equation, especially how to develop efficient and stable variable - time - step algorithms. Specifically:
1. **Unconditional Nonlinear Stability and Accuracy**:
- The paper proposes a class of variable - time - step schemes based on the Dahlquist - Liniger - Nevanlinna (DLN) method. These schemes have unconditional nonlinear stability and are second - order accurate on non - uniform time grids.
2. **Treatment of Nonlinear Terms**:
- For the nonlinear terms in the Allen - Cahn equation, the paper combines two efficient time algorithms: the partially implicit correction algorithm and the scalar auxiliary variable (SAV) algorithm. Both of these methods ensure the unconditional long - term stability of the model energy under any time - step sequence.
3. **Error Analysis**:
- The paper provides a strict error analysis of the partially implicit correction algorithm under variable - time - steps, proves its convergence at first - order time accuracy, and achieves second - order time accuracy on uniform time grids.
4. **Adaptive Time - Step Algorithm**:
- An adaptive - time algorithm based on the local truncation error (LTE) criterion is proposed to improve computational efficiency. These algorithms estimate LTE using certain explicit time - step methods and adjust the time - step size according to the estimation results.
5. **Numerical Verification**:
- The properties of the proposed adaptive DLN method are verified through one - dimensional and two - dimensional numerical tests, demonstrating the effectiveness and superiority of these methods in practical applications.
### Representation of Mathematical Formulas
- The form of the Allen - Cahn equation is:
\[
\begin{cases}
u_t-\epsilon^2\Delta u + f(u)=0, & x\in\Omega, \quad 0 < t\leq T, \\
\frac{\partial u}{\partial\vec{n}} = 0, & \text{on }\partial\Omega,
\end{cases}
\]
where \( f(u)=u^3 - u\), and the model energy \( E(u)\) is defined as:
\[
E(u):=\int_{\Omega}\left(\frac{\epsilon^2}{2}|\nabla u|^2+F(u)\right)dx,
\]
satisfying the energy dissipation law:
\[
\frac{d}{dt}E(u(t))=-\int_{\Omega}u_t^2dx\leq0.
\]
- The discrete form of the DLN method:
\[
\sum_{\ell = 0}^{2}\alpha_{\ell}y_{n - 1+\ell}=\hat{k}_n g\left(\sum_{\ell = 0}^{2}\beta(n)_{\ell}t_{n - 1+\ell},\sum_{\ell = 0}^{2}\beta(n)_{\ell}y_{n - 1+\ell}\right),
\]
where the coefficients \(\alpha_{\ell}\) and \(\beta(n)_{\ell}\) are as follows:
\[
\begin{bmatrix}
\alpha_2 \\ \alpha_1 \\ \alpha_0
\end{bmatrix}=\begin{bmatrix}
\frac{1}{2}(\theta + 1) \\ -\theta \\ \frac{1}{2}(\theta - 1)
\end{bmatrix},
\quad
\begin{bmatrix}
\beta(n)_2 \\ \beta(n)_1 \\ \beta(n)_0
\end{bmatrix}
\]