zbMATH — the first resource for mathematics

On the unconditionally gradient stable scheme for the Cahn-Hilliard equation and its implementation with Fourier method. (English) Zbl 1305.65250
Summary: We implement a nonlinear unconditionally gradient stable scheme by Eyre, within the Fourier method framework for the long-time numerical integration of the Allen-Cahn and Cahn-Hilliard equations, which are gradient flows of the Allen-Cahn energy. We propose a new iterative procedure to solve the nonlinear scheme. When the iterative scheme is applied to the Allen-Cahn equation, we show that the nonlinear iteration is a contractive mapping in the \(L^2\) norm for large time steps. For the Cahn-Hilliard equation, we establish that the proposed iterative scheme converges with a time step constraint. Further, we numerically demonstrate that the iterative scheme converges for large time steps. The scheme allows for spectral accuracy in space and fast simulation of the dynamics in high dimensions while preserving the discrete form of the energy law. For the general potential well, we present the gradient stable scheme without introducing extra stabilizing terms. Therefore, the numerical error introduced by the operator splitting is reduced to its minimum.

65T50 Numerical methods for discrete and fast Fourier transforms
82C26 Dynamic and nonequilibrium phase transitions (general) in statistical mechanics
65M12 Stability and convergence of numerical methods for initial value and initial-boundary value problems involving PDEs
80A22 Stefan problems, phase changes, etc.
Full Text: DOI