Stability and Robustness of Time-discretization Schemes for the Allen-Cahn Equation via Bifurcation and Perturbation Analysis

Wenrui Hao,Sun Lee,Xiaofeng Xu,Zhiliang Xu

2024-06-26

Abstract:The Allen-Cahn equation is a fundamental model for phase transitions, offering critical insights into the dynamics of interface evolution in various physical systems. This paper investigates the stability and robustness of frequently utilized time-discretization numerical schemes for solving the Allen-Cahn equation, with focuses on the Backward Euler, Crank-Nicolson (CN), convex splitting of modified CN, and Diagonally Implicit Runge-Kutta (DIRK) methods. Our stability analysis reveals that the Convex Splitting of the Modified CN scheme exhibits unconditional stability, allowing greater flexibility in time step selection, while the other schemes are conditionally stable. Additionally, our robustness analysis highlights that the Backward Euler method converges to correct physical solutions regardless of initial conditions. In contrast, the other methods studied in this work show sensitivity to initial conditions and may converge to incorrect physical solutions if the initial conditions are not carefully chosen. This study introduces a comprehensive approach to assessing stability and robustness in numerical methods for solving the Allen-Cahn equation, providing a new perspective for evaluating numerical techniques for general nonlinear differential equations.

Numerical Analysis

What problem does this paper attempt to address?

The main problem that this paper attempts to solve is to evaluate the stability and robustness of numerical methods for time discretization when solving the Allen - Cahn equation. Specifically, the paper focuses on the following commonly - used time - discretization methods:

1. **Backward Euler method**
2. **Crank - Nicolson scheme**
3. **Convex Splitting of Modified Crank - Nicolson**
4. **Diagonally Implicit Runge - Kutta (DIRK)**

### Research Background

The Allen - Cahn equation is a fundamental partial differential equation that describes the phase - transition process and is widely used in materials science, image processing, and other fields that require modeling of interface dynamics. Due to the non - linearity of this equation and the diffusive interface in the solution, developing robust and stable numerical methods has been a long - term challenge for accurate simulation.

### Research Objectives

The paper systematically studies the stability and robustness of the above - mentioned methods through bifurcation analysis and perturbation analysis. Specific objectives include:

1. **Stability**: Defined as the uniqueness of $\phi^{n + 1}$ given $\phi^n$, revealing the upper limit of the time step of the numerical scheme.
2. **Robustness**: Defined as the uniqueness of $\phi^n$ given $\phi^{n + 1}$, indicating the accuracy of the numerical scheme in converging to the physical solution.

### Main Findings

1. **Backward Euler method**:
- **Stability**: Unconditionally stable, allowing for a larger choice of time step.
- **Robustness**: Converges to the correct physical solution regardless of the initial conditions.

2. **Crank - Nicolson method**:
- **Stability**: Conditionally stable, with the time step $\Delta t\leq2\epsilon^2$.
- **Robustness**: Sensitive to the initial conditions and may converge to the wrong physical solution if the initial conditions are not properly chosen.

3. **Convex Splitting of Modified Crank - Nicolson**:
- **Stability**: Unconditionally stable, allowing for a larger choice of time step.
- **Robustness**: Similar to the Crank - Nicolson method, sensitive to the initial conditions.

4. **Diagonally Implicit Runge - Kutta method**:
- **Stability**: Conditionally stable, with the time step $\Delta t\leq\frac{\epsilon^2}{\max_i a_{ii}}$.
- **Robustness**: Sensitive to the initial conditions and may converge to the wrong physical solution.

### Conclusion

Through this research, the author provides a comprehensive framework for evaluating the stability and robustness of numerical methods for solving the Allen - Cahn equation. These findings are helpful for guiding the development of efficient and reliable computational methods for simulating phase - field models. In particular, the Backward Euler method performs well in dealing with stiff problems because of its unconditional stability and robustness.