Numerical study of the 2D Cahn-Hilliard model of phase separation with logarithmic potential
R. Abazari1, K. Yildirim2
1Department of Mathematics, University of Mohaghegh Ardabili, Ardabil, Iran
2Mus Alparslan University, Mus, Turkey
Keywords: small Phase separation, Cahn--Hilliard model, Splitting scheme, Crank--Nicolson method, Error estimation
Abstract
In this paper, we have approximated the solutions of the Cahn-Hilliard equation (CH) with the logarithmic potential function which arises in the modeling of phase separation of binary alloys. The CH equation is a high-order nonlinear equation,consequently, utilizing a common difference scheme on the CH equation causes long stencil schemes. To resolve the faults of long stencil schemes, we split the CH equation to a second-order system under the Neumann boundary condition and we applied a second-order scheme based on the 2D Crank-Nicolson method to discrete it. The uniqueness and error estimation of the approximated solution is proved. Also, preserving the conservation of mass and decreasing the total energy are investigated. Finally, to confirm the theoretical results, three examples with various initial conditions are presented.
|
