A Second-Order, Weakly Energy-Stable Pseudo-spectral Scheme for the Cahn-Hilliard Equation and ...
发布时间:2022-06-16 | 发布人:知识服务与技术支持部标题: A Second-Order, Weakly Energy-Stable Pseudo-spectral Scheme for the Cahn-Hilliard Equation and Its Solution by the Homogeneous Linear Iteration Method
作者: Cheng, Kelong;Wang, Cheng);Wise, Steven M.;Yue, Xingye
来源出版物: JOURNAL OF SCIENTIFIC COMPUTING
出版年: DEC 2016
作者关键词: Cahn-Hilliard equation;Second order convex splitting;Energy stability;Fourier pseudo-spectral approximation;Linear iteration;Contraction mapping
研究方向: Mathematics
第一地址: 西南科技大学
入藏号:WOS:000387443600007
中国科学院文献情报中心期刊分区: 数学2区
摘要:
We present a second order energy stable numerical scheme for the two and three dimensional Cahn-Hilliard equation, with Fourier pseudo-spectral approximation in space. A convex splitting treatment assures the unique solvability and unconditional energy stability of the scheme. Meanwhile, the implicit treatment of the nonlinear term makes a direct nonlinear solver impractical, due to the global nature of the pseudo-spectral spatial discretization. We propose a homogeneous linear iteration algorithm to overcome this difficulty, in which an (where s the time step size) artificial diffusion term, a Douglas-Dupont-type regularization, is introduced. As a consequence, the numerical efficiency can be greatly improved, since the highly nonlinear system can be decomposed as an iteration of purely linear solvers, which can be implemented with the help of the FFT in a pseudo-spectral setting. Moreover, a careful nonlinear analysis shows a contraction mapping property of this linear iteration, in the discrete norm, with discrete Sobolev inequalities applied. Moreover, a bound of numerical solution in norm is also provided at a theoretical level. The efficiency of the linear iteration solver is demonstrated in our numerical experiments. Some numerical simulation results are presented, showing the energy decay rate for the Cahn-Hilliard flow with different values of .