Abstract:In this work, we present an efficient approach for the spatial and temporal discretization of the nonlocal Allen-Cahn equation, which incorporates various double-well potentials and an integrable kernel, with a particular focus on a non-smooth obstacle potential. While nonlocal models offer enhanced flexibility for complex phenomena, they often lead to increased computational costs and there is a need to design efficient spatial and temporal discretization schemes, especially in the non-smooth setting. To address this, we propose first- and second-order energy-stable time-stepping schemes combined with the Fourier collocation approach for spatial discretization. We provide energy stability estimates for the developed time-stepping schemes. A key aspect to our approach involves a representation of a solution via proximal operators. This together with the spatial and temporal discretizations enables direct evaluation of the solution that can bypass the solution of nonlinear, non-smooth, and nonlocal system. This method significantly improves computational efficiency, especially in the case of non-smooth obstacle potentials, and facilitates rapid solution evaluations in both two and three dimensions. We provide several numerical experiments to illustrate the effectiveness of our approach.
What problem does this paper attempt to address?
This paper attempts to solve the problem of efficient solution in the spatial and temporal discretization of the non - local Allen - Cahn equation. Specifically, the paper focuses on how to design efficient numerical methods to deal with non - local Allen - Cahn equations containing various double - well potential functions (such as regular potential, singular potential and obstacle potential). These equations are used to describe phase - transition phenomena in fields such as materials science, image processing and mean curvature flow.
### Main problems
1. **Computational efficiency**: Although non - local models provide greater modeling flexibility, they usually lead to higher computational costs. Therefore, it is necessary to design efficient spatio - temporal discretization schemes, especially when dealing with non - smooth potential functions.
2. **Energy stability**: The proposed time - stepping scheme needs to maintain energy stability to ensure the reliability and accuracy of the numerical solution.
3. **Solution of nonlinear and non - smooth systems**: Traditional numerical methods often require complex iterative solution processes when dealing with nonlinear and non - smooth systems, while the method proposed in this paper directly solves by using proximal operators, thus simplifying the calculation process.
### Solutions
1. **Time discretization**:
- First - order and second - order energy - stable time - stepping schemes are proposed.
- The semi - implicit Euler time - stepping method is used, combined with the fast Fourier transform (FFT) for explicit discretization of convolution.
- Represent the solution by approximate operators, thereby avoiding the complexity of solving nonlinear and non - smooth systems.
2. **Spatial discretization**:
- The Fourier collocation method is used for spatial discretization.
- Use the discrete Fourier transform (DFT) to efficiently calculate the convolution term.
3. **Application of approximate operators**:
- For the obstacle potential, the approximate operator is simplified to a point - by - point projection operator.
- For the regular potential, the approximate operator can be calculated by the Cardano formula.
- For the Flory - Huggins potential, since there is no known closed - form solution, Newton's method or the bisection method can be used for numerical solution.
### Numerical experiments
The paper verifies the effectiveness and accuracy of the proposed method through a series of numerical experiments, and shows the application effect of this method in 2D and 3D problems.
### Formula summary
- **Non - local diffusion operator**:
\[
Lu(x)=\int_{\Omega}(u(x)-u(y))\gamma(x - y)dy = c_{\gamma}u(x)-(\gamma * u)(x)
\]
where \( c_{\gamma}=\int_{\Omega}\gamma(x - y)dy\).
- **Non - local energy**:
\[
E(u)=\frac{1}{4}\int_{\Omega}\int_{\Omega}(u(x)-u(y))^{2}\gamma(x - y)dxdy+\int_{\Omega}F(u)dx
\]
- **Double - well potential function**:
\[
F(u)=f_{0}(u)+\psi(u)
\]
where \( f_{0}(u)=\frac{c_{F}}{2}(1 - u^{2})\), and \(\psi\) can be a regular potential, Flory - Huggins potential or obstacle potential.
- **Approximate operator**:
\[
\text{prox}_{\lambda\psi}(v)=\arg\min_s\left(\frac{\lambda\psi(s)}{2}+\frac{|s - v|^{2}}{2}\right)
\]
Through these methods and formulas, the paper provides an efficient and stable solution scheme for the non - local Allen - Cahn equation.