Abstract:In this paper, we present several unconditionally energy-stable invariant energy quadratization (IEQ) finite element methods (FEMs) with linear, first- and second-order accuracy for solving both the Cahn-Hilliard equation and the Allen-Cahn equation. For time discretization, we compare three distinct IEQ-FEM schemes that position the intermediate function introduced by the IEQ approach in different function spaces: finite element space, continuous function space, or a combination of these spaces. Rigorous proofs establishing the existence and uniqueness of the numerical solution, along with analyses of energy dissipation for both equations and mass conservation for the Cahn-Hilliard equation, are provided. The proposed schemes' accuracy, efficiency, and solution properties are demonstrated through numerical experiments.
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop unconditionally energy - stable finite element methods (FEMs) for solving the Cahn - Hilliard equation and the Allen - Cahn equation. Specifically, the paper focuses on the following points:
1. **Energy stability**: Ensure that the numerical scheme maintains the energy dissipation property at any time step, which is crucial for the accuracy of long - time simulations.
2. **Linearization treatment**: By introducing intermediate variables, the nonlinear terms are transformed into linear forms, thus avoiding the complexity brought by iterative techniques.
3. **High - order accuracy**: Propose numerical methods with linear, first - order, and second - order accuracy to improve computational efficiency and accuracy.
### Specific problem description
Both the Cahn - Hilliard equation and the Allen - Cahn equation are important models for describing phase - change processes and are widely used in fields such as materials science and fluid mechanics. These equations are essentially nonlinear partial differential equations, and it is difficult to obtain analytical solutions. Therefore, efficient numerical methods are required to solve them.
#### Cahn - Hilliard equation
The Cahn - Hilliard equation describes the phase separation process of a two - phase system, and its mathematical form is:
\[
u_t=\nabla\cdot(M(u)\nabla w)\quad\text{in }\Omega\times(0, T],
\]
\[
w = -\epsilon^2\Delta u+F'(u)\quad\text{in }\Omega\times(0, T],
\]
The boundary conditions are:
\[
\nabla u\cdot n = 0,\quad\nabla w\cdot n = 0\quad\text{on }\partial\Omega\times[0, T],
\]
The initial condition is:
\[
u = u_0\quad\text{on }\Omega\times\{t = 0\},
\]
where $\Omega\subseteq\mathbb{R}^d$ is a bounded region, $\epsilon$ is a small parameter, $M(u)\geq0$ represents the mobility function, $u_0(x)$ is the initial data, $n$ is the unit outer normal vector of the boundary $\partial\Omega$, and $F(u)$ is the nonlinear volume potential energy.
#### Allen - Cahn equation
The Allen - Cahn equation describes the interface evolution process, and its mathematical form is:
\[
u_t=\epsilon^2\Delta u - F'(u)\quad\text{in }\Omega\times(0, T],
\]
The boundary conditions are:
\[
\nabla u\cdot n = 0\quad\text{on }\partial\Omega\times[0, T],
\]
The initial condition is:
\[
u = u_0\quad\text{on }\Omega\times\{t = 0\}.
\]
### Solution
The paper proposes several finite element methods (FEMs) based on the invariant energy - quadraticization (IEQ) method. These methods can unconditionally maintain energy stability and have linear, first - order, and second - order accuracy. Specifically, the paper compares three different IEQ - FEM schemes, which handle the intermediate variables introduced by the IEQ method in different function spaces.
1. **Method 1 (BDF1 - IEQ - FEM1)**: Directly incorporate the intermediate variable $U$ into the finite element space $V_h$ and treat it the same as other variables. This method will lead to an augmented linear system and increase the computational cost.
2. **Method 2 (BDF1 - IEQ - FEM2)**: Place the intermediate variable $U$ in the continuous function space $C^0(\Omega)$ instead of the finite element space $V_h$. This method can calculate the intermediate variable independently and point - by - point, improving computational efficiency, but may lead to inconsistent decay of the approximate energy.