Abstract:Two new numerical schemes to approximate the Cahn-Hilliard equation with degenerate mobility (between stable values 0 and 1) are presented, by using two different non-centered approximation of the mobility. We prove that both schemes are energy stable and preserve the maximum principle approximately, i.e. the amount of the solution being outside of the interval [0,1] goes to zero in terms of a truncation parameter. Additionally, we present several numerical results in order to show the accuracy and the well behavior of the proposed schemes, comparing both schemes and the corresponding centered scheme.
What problem does this paper attempt to address?
The core problem that this paper attempts to solve is to develop new numerical methods for solving the Cahn - Hilliard equation with degenerate mobility. Specifically, the author proposes two new numerical schemes aiming to solve the following problems:
1. **Energy stability**: Ensure that the numerical scheme keeps the total free energy of the system from increasing with time during the simulation process, that is, it satisfies the energy - stable condition of the gradient - flow structure.
2. **Approximate preservation of the maximum principle**: Ensure that the value range of the numerical solution is as close as possible to the physically reasonable interval [0, 1]. It is mentioned in the paper that by introducing the truncation parameter ε, the part of the numerical solution outside the [0, 1] interval tends to zero as ε decreases.
3. **Handling degenerate mobility**: The mobility function \(M(\phi)\) in the Cahn - Hilliard equation degenerates to zero near the pure - phase regions (\(\phi = 0\) and \(\phi = 1\)), which brings difficulties to numerical solution. The author approximates the mobility function by non - centered difference, thereby avoiding directly dealing with singular terms and making the numerical scheme more stable and effective.
### Specific problem description
- **Equation background**: The Cahn - Hilliard equation is a fourth - order partial differential equation used to describe the interface dynamics in a two - phase system. The form of the equation is:
\[
\phi_t-\nabla\cdot\left(M(\phi)\nabla\frac{\delta E}{\delta\phi}\right) = 0,
\]
where \(M(\phi)\) is the mobility function and \(E(\phi)\) is the free - energy functional of the system.
- **Degenerate mobility**: When the mobility \(M(\phi)\) degenerates to zero in the pure - phase regions, traditional numerical methods may fail or be unstable. Therefore, new numerical schemes need to be designed to handle this degenerate situation.
- **Requirements for numerical schemes**: The new numerical schemes should not only ensure energy stability but also preserve the maximum principle of the solution as much as possible, that is, the value range of the solution should be limited within [0, 1] as much as possible.
### Main contributions of the paper
- **Proposing two new schemes**: Based on the non - centered difference approximation of the mobility function, two new numerical schemes (Gε - scheme and Jε - scheme) are proposed, and their energy stability and approximate preservation of the maximum principle are proved.
- **Theoretical analysis and numerical verification**: Through strict mathematical derivations, the energy stability and approximate preservation of the maximum principle of the new schemes are proved, and these properties are verified by numerical experiments.
In conclusion, this paper aims to develop efficient and stable numerical methods applicable to the Cahn - Hilliard equation with degenerate mobility, and solves the problems of numerical instability and unreasonable solution values that traditional methods may encounter when dealing with degenerate mobility.